‘Mutiny on the Bounty’: the genetic history of Norfolk Island reveals extreme gender-biased admixture
© Benton et al. 2015
Received: 3 June 2015
Accepted: 28 August 2015
Published: 3 September 2015
The Pacific Oceania region was one of the last regions of the world to be settled via human migration. Here we outline a settlement of this region that has given rise to a uniquely admixed population. The current Norfolk Island population has arisen from a small number of founders with mixed Caucasian and Polynesian ancestry, descendants of a famous historical event. The ‘Mutiny on the Bounty’ has been told in history books, songs and the big screen, but recently this story can be portrayed through comprehensive molecular genetics. Written history details betrayal and murder leading to the founding of Pitcairn Island by European mutineers and the Polynesian women who left Tahiti with them. Investigation of detailed genealogical records supports historical accounts.
Using genetics, we show distinct maternal Polynesian mitochondrial lineages in the present day population, as well as a European centric Y-chromosome phylogeny. These results comprehensively characterise the unique gender-biased admixture of this genetic isolate and further support the historical records relating to Norfolk Island.
Our results significantly refine previous population genetic studies investigating Polynesian versus Caucasian diversity in the Norfolk Island population and add information that is beneficial to future disease and gene mapping studies.
KeywordsNorfolk Island Mitochondria Y-chromosome Genetic isolate Population history
The Bounty set for Tahiti in 1787, with the goal of securing and distributing breadfruit for the slave trade. After an unexpectedly long stopover in Tahiti, the majority of the Bounty crew overthrew Captain William Bligh. Led by Lieutenant Fletcher Christian, 18 mutineers returned to Tahiti. Expecting to be discovered by the British, Christian along with 8 other mutineers, 6 Tahitian men and 18 Tahitian women (plus a baby girl) fled to the Pitcairn Islands, where the Bounty was burnt to avoid detection and to eliminate the only viable escape route.
Conflict became rife, with numerous murders as well as suicide, leaving all Tahitian men and 7 mutineers dead. When Ned Young succumbed to respiratory failure in 1800, John Adams became the sole surviving male on Pitcairn. Slowly the island’s population grew, later bolstered by three further male European settlers—John Buffett and John Evans in 1823 and George Nobbs in 1828 . These three men were the only outsiders to settle permanently on the island and subsequently married into the community.
Continual population growth and the concomitant pressure on resources led the British Government to allow the Pitcairn Islanders to settle Norfolk Island. On the 8th June 1856, a total of 193 inhabitants of Pitcairn Island were relocated. Today Norfolk Island has a population of about 1600 permanent residents, with the 2006 census reporting that approximately half of the permanent population are descendants of the Pitcairn founders .
Previously our group have illustrated the unique genomic properties of the Norfolk Island pedigree [5, 6], as well as its use in disease gene mapping studies [7–10]; establishing that this population isolate is able to identify both common and rare variants associated with complex phenotypes that may be difficult to study in larger but more heterogeneous populations. An important part of this process is identifying and understanding genetic admixture, it is critical to account for population substructure when conducting disease gene mapping studies. In this short report, we detail significant improvements to our previous estimates for both maternal and paternal ancestry in the unique population of Norfolk Island. These findings are interesting from a population genetics point of view and in turn highlight a key part of Oceanic history. Our results further comprehensively characterise the unique gender-biased admixture of this genetic isolate which is an important consideration when studying complex disease traits.
Cohort collection and ethics
Accurate and detailed historical accounts have been used by genealogists to create and maintain a well-documented database of the entire Norfolk Island population, spanning all the way back to the original founders. This pedigree has been drawn up and is maintained in a genealogy program known as Brother’s Keeper. The pedigree includes ~5700 individuals coalescing over 11 generations or 200 years back to the original 9 European sailors and 12 Tahitian women. The Norfolk Island Health Study which has already been well established in previous research sampled individuals from the lower four generations of the pedigree and included 386 (64 %) individuals possessing lineages back to the founders and 216 individuals (36 %) who were considered to be new founders and did not show direct ancestral links. An updated core pedigree was constructed using this information and genetic information as it became available through genetic studies. Currently, the core pedigree structure contains those individuals that are most closely related and coalesce directly back to the original founders. The Norfolk Island Health Study (NIHS) has already been well established in previous research [7–10]. In this study, we used a representative group of Norfolk Island individuals selected from the pedigree, meaning that they relate back to the original founders, and we have phenotype and genotype information for them. The total number of core pedigree members selected was 330 (this was adjusted to exclude individuals under the age of 18 years), which consisted of 152 males and 178 females. All individuals gave written informed consent. Ethical approval was granted prior to the commencement of the study by the Griffith University Human Research Ethics Committee (ethical approval no: 1300000485) and the project was carried out in accordance with the relevant guidelines, which complied with the Helsinki Declaration for human research.
Orthographic map creation
A custom script was written by David Eccles to convert shapefile information into SVG images (https://github.com/gringer/bioinfscripts/blob/master/perlshaper.pl). This script was used with both 1:50 M resolution world political boundaries from NED (http://www.naturalearthdata.com/) and 1:1 M boundaries from GADM (http://www.gadm.org/) to create the local map and orthographic world projection in Fig. 1. Maps were combined together and overlaid with additional text using Inkscape (https://inkscape.org/en/).
Detection of the mitochondrial 9 bp repeat/deletion
To explore the presence of Polynesian mitochondrial genomes in the Norfolk Island population 23 individuals with direct maternal relationship to the Polynesian founders were selected. Primer sequences were designed to target the mitochondrial COII/tRNAlys region known to contain the 9 bp deletion in haplogroup B carriers; Forward primer 5′ AGGGCCCGTATTTACCCTATAG 3′, Reverse primer 5′ ATTTAGTTGGGGCATTTCACTG 3′. The target amplicon had an expected size of 133-bp for samples without the 9-bp deletion and 124-bp for samples containing the deletion. Briefly, 20 ng of genomic DNA was amplified with PCR Premix Optimisation Buffer E supplied by Epicentre Technologies, 0.3 μM of each mtDEL primer and AmpliTaq polymerase in a 20 μL reaction volume. Amplification was performed on a PE Applied Biosytems GeneAmp System 9700 with a 95 °C initial denaturation for 4 min followed by 30 cycles of 95 °C for 30 s, 60 °C for 1 min and 72 °C for 1 min, with a final extension at 72 °C for 7 min. Genotyping of the 9-bp deletion involved the use of fluorescently labelled primers. PCR products were incorporated into a formamide and ROX size standard mix prior to loading on an Applied Biosystems 310 Genetic Analyser. Fragment lengths were accurately determined after detection by automated gel capillary electrophoresis. Full details of the above method and RFLP data are available online: https://www120.secure.griffith.edu.au/rch/items/0e72bbad-feb0-c2bf-54d7-fef9adf7f5e3/1/.
Preliminary mitochondrial genome resequencing
Initially, the mitochondrial DNA of the selected 23 samples was resequenced with the VariantSEQr resequencing system (Applied Biosystems, Foster City, CA). The primer set (RS000056015_01 mitoALL) included a total of 46 amplicons. Each primer pair is tailed with universal M13 sequence and is designed for universal PCR and sequencing conditions. PCR was performed with a final volume of 10 μL containing 1.0 μL of genomic DNA, 2.0 μL mitoSEQr primer mix, 5.0 μL of 2X AmpliTaq Gold PCR Master Mix and 1.6 μL of 50 % UltraPure glycerol on a AB9700 with the following cycling conditions: Heat activation at 96 °C for 5 min, followed by 40 cycles of 94 °C for 30 s, 60 °C for 45 s and 72 °C for 45 s. A final extension of 72 °C for 10 min was included before cooling to hold at 4 °C. Prior to sequencing, PCR products were cleaned-up with 2 μL of ExoSAP-IT and 37 °C incubation for 30 min followed by heat inactivation at 80 °C for 15 min. A final volume of 10 μL contained 4.0 μL of BigDye Terminator Ready Reaction Mix v3.1, 1.0 μL of M13 primer, 3.0 μL of deionised water and 2.0 μL of PCR product. Cycling conditions were initiated with heat activation of 96 °C for 1 min, followed by 25 cycles of 96 °C for 10 s, 50 °C for 5 min and 60 °C for 4 min, a final hold step of 4 °C was included. To clean the sequence prior to electrophoresis, 2.5 μL of 125 mM EDTA was added to each tube and mixed followed by 30 μL 100 % ethanol. This solution was allowed to incubate at room temperature for 15 min. Following incubation, samples were centrifuged at 2000×g for 45 min with removal of supernatant prior to addition of 30 μL of 70 % ethanol and a further centrifugation of 15 min at 2000×g. Again, supernatant was removed, and samples were allowed to air dry before being resuspended in 10 μL of Hi-Di Formamide. Electrophoresis was performed on an AB3730xl (POP7) DNA Analyser running buffer with EDTA using a 36-cm array.
Sample preparation and complete mitochondrial genome sequencing
To facilitate more accurate and high throughput sequencing of complete mitochondrial genomes, we developed an in house method using a combination of long-range PCR and the Life Technologies Ion Torrent Next Generation Sequencing platform. This was performed on a representative cross section of the Norfolk Island pedigree (n = 322).
Mitochondrial DNA was enriched and purified before undergoing library preparation. Samples were amplified by long-range PCR, utilising two primer pairs which produce overlapping fragments covering the entire human mitochondrial genome. Primer sequences used were as follows: fragment 1 from 569 to 9819 b.p fwrd AAC CAA ACC CCA AAG ACA CC and rvrs GCC AAT AAT GAC GTG AAG TCC fragment 2 from 9611 to 626 b.p fwrd TCC CAC TCC TAA ACA CAT CC and rvrs TTT ATG GGG TGA TGT GAG CC. Primer sequences are given in the 5′ to 3′ orientation. Final reagent concentrations used were 1x Buffer, 0.5 mM dNTPs, 200 nM each primer, 3 % DMSO and 0.05 units/μL Roche expand long-range enzyme with 100 ng genomic DNA input in a 50 μL reaction volume. The optimum thermocycling conditions consisted of a two-stage PCR with an initial activation step at 92 °C for 2 min followed by 9 cycles of 92 °C for 10 s, 55 °C for 15 s, 68 °C for 9 min and 19 cycles of 92 °C for 15 s, 60 °C for 20 s, 68 °C for 9 min and a final extension of 7 min at 68 °C followed by an 8 °C hold. PCR products were run on 1 % agarose gels to accommodate for large fragment sizes at 70 V for 90 min and visualised under UV light. Ethidium bromide was used for staining at a 4 % concentration. Negative controls were included in each PCR run to check for contamination and ensure adequate quality control.
Mitochondrial PCR fragments were cleaned using QIAquick post PCR cleanup columns to remove any excess dNTPs or other reagents and were then quantified using Agilent DNA 12000 chips on a bioanalyser. Overlapping long-range PCR fragments were pooled together in equimolar amounts for each sample. After pooling and mixing 100 ng was aliquoted for library preparation. Samples underwent library preparation using NEBNext® Fast DNA Library Prep Set for Ion Torrent following the manufacturer’s protocol. Physical shearing using a Biorupter was found to produce a much more even distribution of fragments than an enzymatic approach. All samples were physically sheared post fragment pooling for 15 cycles (30 s ON, 30 s OFF) on high at 4 °C. Barcoding of samples was performed using Bioo Scientific NEXTflex DNA Barcodes for Ion PGM. The E-Gel SizeSelect system was used to select a target size of 350-bp libraries for use with 200-bp sequencing chemistry. Agencourt AMPure XP beads from Beckman and Coulter were used for all cleaning steps. After the final amplification and purification step all libraries were quantified on a Bioanalyser using DNA1000 chips. Samples for each multiplex were pooled in equimolar amounts and diluted to 26 pM, with 48 samples being plexed into a single sequencing reaction on 316 chips. Template preparation and sequencing was undertaken according to manufacturer’s instructions for the Ion Torrent platform.
A custom pipeline was developed to accommodate QC, alignment, assembly and analysis, with scripts available for download at a GitHub repository (https://github.com/gringer/phylotreescripts). Briefly, raw sequence reads were aligned to the Reconstructed Sapiens Reference Sequence (RSRS)  using Bowtie2 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml). SAMtools (http://www.htslib.org/) was used to produce Binary Alignment (BAM) files, which were processed into VCF files, then converted to consensus FASTA sequences. These FASTA sequences were compared with variant reference information from Phylotree in order to find mitochondrial haplotypes that most closely represented the consensus sequences for that individual.
All phylogenetic reconstructions were generated in R (http://www.r-project.org/), using the APE package (http://cran.r-project.org/web/packages/ape/). Reconstruction and labelling of the Norfolk Island pedigree diagram was performed using the Kinship2 R package (http://cran.r-project.org/web/packages/kinship2/).
Y chromosome analysis
DNA extraction and genome-wide genotyping procedures for the NIHS are previously documented . All SNPs located on the Y chromosome were QC filtered and extracted from genotype files using PLINK v1.07 (http://pngu.mgh.harvard.edu/~purcell/plink/); this is available as Additional file 1, Build36 map information is available as Additional file 2. The resultant SNP data were imported to R with haplotype and phylogenetic analysis performed using APE. SNP labelling was performed in accordance with the minimal skeleton Y tree .
Results and discussion
Mitochondrial evidence for Polynesian female founders
Y-chromosome evidence for European male founders
The results presented here show excellent concordance with written historical accounts and greatly improve upon prior findings, refining estimates of sex-dependent admixture bias and specific haplotypic resolution within the Norfolk Island population. Sequencing of the mitochondrial genomes from a large majority of the pedigree members allowed us to ascertain that approximately 40 % of the current day population have Polynesian derived mtDNA. This is a large increase over our previous estimate of 24 % described in McEvoy et al. . Similarly, we have been able to use information obtained from dense SNP arrays to infer the haplotypic status of 240 males based on Y-chromsome genotype data, which is a substantial improvement upon previous work that was limited to a small number of STRs. Overall, these new results add much more clarity to the genetic ancestral picture of the Norfolk Island population.
Knowledge of the population substructure caused by the divergence of genetically distinct ancestral populations is important when conducting disease gene mapping studies. For example, at the very least admixture requires adjustment when performing statistical tests of association within and between populations. The data presented here, unlike previous results, will allow accurate assessment and correction of sex-specific lineages when conducting disease gene mapping studies. In particular, when performing disease gene mapping studies accurate assignment of the mitochondrial haplogroups that distinguish between and within Polynesian and European ancestry subgroups should allow real pathological point variants in the mitochondrial genome to be discerned from those increasing in frequency due to population genetic factors such as genetic drift.
The geographic isolation of Norfolk Island, alongside a small number of European and Polynesian founders—including the additional 3 European men settling Pitcairn later—have established a unique genetic isolate. Consequently, extreme gender-biased admixture contributes to this unique genomic structure and provides exciting opportunities for future population genetics and disease gene mapping studies.
Her Majesty’s Ship
Norfolk Island Health Study
polymerase chain reaction
restriction fragment length polymorphism
Single Nucleotide Polymorphism
short tandem repeat
This research was supported by funding from a National Health and Medical Research Council of Australia project grant and utilised infrastructure purchased with Australian Government EIL SuperScience Funds as part of the Therapeutic Innovation Australia—Queensland Node project. MB was supported by a Corbett Postgraduate Research Scholarship, and RL was partially supported by a Corbett Research Fellowship. The authors would like to acknowledge Dr Peter Underhill for his comment and suggestions aiding in the refinement of the Y-chromosome phylogeny, and David Eccles for providing the perl scripts to generate the world map. Finally, we extend our appreciation to the Norfolk Islanders who volunteered for this study.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Edgecombe J. Norfolk Island - South Pacific: island of history and many delights. Thornleigh: Australia: J.M. Edgecombe; 1999.Google Scholar
- Hoare M. Norfolk Island: a revised and enlarged history 1774–1998. Rockhampton: Central Queensland University Press; 1999.Google Scholar
- Refshauge WF, Walsh RJ. Pitcairn Island - fertility and population-growth, 1790–1856. Ann Hum Biol. 1981;8:303–12.View ArticlePubMedGoogle Scholar
- Mathews SP. NORFOLK ISLAND: Census of population and housing, 8 August 2006. Statistical report on characteristics of population and dwellings. 2006.Google Scholar
- Macgregor S, Bellis C, Lea RA, Cox H, Dyer T, Blangero J, et al. Legacy of mutiny on the Bounty: founder effect and admixture on Norfolk Island. Eur J Hum Genet. 2010;18:67–72.PubMed CentralView ArticlePubMedGoogle Scholar
- McEvoy BP, Zhao ZZ, Macgregor S, Bellis C, Lea RA, Cox H, et al. European and Polynesian admixture in the Norfolk Island population. Heredity (Edinb). 2010;229–234.
- Bellis C, Hughes RM, Begley KN, Quinlan S, Lea RA, Heath SC, et al. Phenotypical characterisation of the isolated Norfolk Island population focusing on epidemiological indicators of cardiovascular disease. Hum Hered. 2006;60:211–9.View ArticleGoogle Scholar
- Bellis C, Cox HC, Ovcaric M, Begley KN, Lea RA, Quinlan S, et al. Linkage disequilibrium analysis in the genetically isolated Norfolk Island population. Heredity (Edinb). 2008;100:366–73.View ArticleGoogle Scholar
- Cox HC, Bellis C, Lea RA, Quinlan S, Hughes R, Dyer T, et al. Principal component and linkage analysis of cardiovascular risk traits in the Norfolk Isolate. Hum Hered. 2009;68:55–64.PubMed CentralView ArticlePubMedGoogle Scholar
- Benton MC, Lea RA, Macartney-Coxson D, Carless MA, Göring HH, Bellis C, et al. Mapping eQTLs in the Norfolk Island genetic isolate identifies candidate genes for CVD Risk Traits. Am J Hum Genet. 2013;93:1087–99.PubMed CentralView ArticlePubMedGoogle Scholar
- Behar DM, Van Oven M, Rosset S, Metspalu M, Loogväli EL, Silva NM, et al. A “copernican” reassessment of the human mitochondrial DNA tree from its root. Am J Hum Genet. 2012;90:675–84.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Oven M, Van Geystelen A, Kayser M, Decorte R, Larmuseau MHD. Seeing the wood for the trees: a minimal reference phylogeny for the human Y chromosome. Hum Mutat. 2014;35:187–91.View ArticlePubMedGoogle Scholar
- Bellis C. Use of the Isolated Norfolk Island Population for Cardiovascular Disease Risk Trait Genetic Analysis. PhD thesis. Griffith University: School of Medical Sciences; 2009.Google Scholar
- Redd AJ, Takezaki N, Sherry ST, McGarvey ST, Sofro AS, Stoneking M. Evolutionary history of the COII/tRNALys intergenic 9 base pair deletion in human mitochondrial DNAs from the Pacific. Mol Biol Evol. 1995;12:604–15.PubMedGoogle Scholar
- Kayser M, Brauer S, Cordaux R, Casto A, Lao O, Zhivotovsky LA, et al. Melanesian and Asian origins of Polynesians: mtDNA and Y chromosome gradients across the Pacific. Mol Biol Evol. 2006;23:2234–44.View ArticlePubMedGoogle Scholar
- Van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009;30:E386–94.View ArticlePubMedGoogle Scholar
- Chiaroni J, Underhill PA, Cavalli-Sforza LL. Y chromosome diversity, human expansion, drift, and cultural evolution. Proc Natl Acad Sci U S A. 2009;106:20174–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Underhill PA, Passarino G, Lin AA, Shen P, Mirazón Lahr M, Foley RA, et al. The phylogeography of Y chromosome binary haplotypes and the origins of modern human populations. Ann Hum Genet. 2001;65(Pt 1):43–62.View ArticlePubMedGoogle Scholar
- Sanchez JJ, Hallenberg C, Børsting C, Hernandez A, Morling N. High frequencies of Y chromosome lineages characterized by E3b1, DYS19-11, DYS392-12 in Somali males. Eur J Hum Genet. 2005;13:856–66.View ArticlePubMedGoogle Scholar