Pedigree Derived Mutation Rate Across the Entire Mitochondrial Genome of the Norfolk Island Population

doi:10.21203/rs.3.rs-823709/v1

Download PDF

Research Article

Pedigree Derived Mutation Rate Across the Entire Mitochondrial Genome of the Norfolk Island Population

https://doi.org/10.21203/rs.3.rs-823709/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Estimates of mutation rates for various regions of the mitochondrial genome (mtGenome) vary widely, depending on whether they are inferred using a phylogenetic approach or obtained directly from pedigrees. Traditionally, only the control region, or small portions of the coding region have been targeted for analysis due to the cost and effort required to produce whole mtGenome Sanger profiles. Here, we report one of the first pedigree derived mutation rates for the entire human mtGenome. The entire mtGenome from 225 individuals originating from Norfolk Island was analysed to estimate the population mutation rate and compared against published mutation rates. These individuals were from 45 maternal lineages spanning 345 generational events. Mutation rates for various portions of the mtGenome were calculated. Nine mutations (including two transitions and seven cases of heteroplasmy) were observed, resulting in a rate of 0.063 mutations/site/million years (95% confidence interval: 0.033 – 0.118). These mutation rates are approximately 17 times higher than estimates derived from phylogenetic analysis with heteroplasmy detected in 13 samples (n=225, 5.8% individuals). Providing one of the first pedigree derived estimates for the entire mtGenome, this study provides a better understanding of mtGenome evolution and has relevance to many research fields, including medicine, anthropology and forensics.

Environmental Policy

Epigenetics & Genomics

Mutation Rate

Whole mtGenome Sequencing

Mitochondrial DNA

Norfolk Island

Pedigree

The mitochondrial genome (mtGenome), and the control region in particular, is the most analysed DNA sequence in human evolutionary studies. As a result, the rate of change in mitochondrial DNA (mtDNA) is relatively well understood. However, since the 1990s, evidence has emerged that mtDNA mutation rates vary widely and are dependent upon the estimation method used and the population studied (e.g. when comparing Swedes and Icelanders, respectively, in Cavelier et al. ¹ and Siguroardottir et al. ².

The rate at which mtDNA mutates is predominantly estimated using phylogenetic or pedigree methods. Phylogenetic analysis of mtDNA relies upon haplotype trees and phylogenetically derived divergence rates, while human pedigree-based mutation rates are estimated by comparison of parent/offspring pairs or using deep-rooted familial lineages at particular loci and counting the number of novel mutations per pair, and this value is then divided by the number of meioses ³. Discrepancy exists between the estimates produced using these two methods, with no consensus on the cause of this discrepancy. Several plausible causes have been suggested, namely: differences in the mutation rate at different locations within the mtGenome; sample size and selection; the effect of natural selection and genetic drift; the occurrence of somatic mutations; the unintended sequencing of nuclear mitochondrial pseudogenes (NUMTs); and the leakage of paternal mtDNA and recombination ⁴. Thus, the choice of which estimate to use in population studies is determined by the specific purpose of the investigation. Phylogenetically based estimations may be more suited to investigations of deep history as this takes into account mutations that have reached a considerable frequency in the population. In contrast, pedigree-based estimations may be more suited to studies of recent history as newly acquired mutations that have not had time to become fixed within the population are considered ^5,6.

The human mtGenome is split into two sections: a large coding region that is responsible cellular energy production and for the gene production for transfer RNA (tRNA), ribosomal RNA (rRNA) translation; and a smaller control region ⁷. The control region also contains the polymorphism-rich hypervariable regions (HVI and HVII), which are traditionally used for comparisons in missing persons cases, criminal cases and for historical military cases ⁸. While the mitochondrial control region has been extensively studied, minimal literature explores the mutation rate for the entire mtGenome. Typically only the control region (or portions of) are examined due to the cost and effort required to produce whole mtGenome Sanger profiles and the perceived stability of the remainder of the genome. The advent of Next Generation Sequencing (NGS) technologies has provided impetus and the potential to expand the current body of knowledge surrounding the prediction of mutation rates. Since it was first introduced in 2005 ⁹, NGS has transformed the application and use of genetic data generation across a number of fields. The advantages of NGS over Sanger Sequencing have been exploited for analysis, particularly in forensics where the increased discrimination power that results from sequencing the entire mtGenome has provided valuable information for human identification ^10–15.

In an effort to estimate the mutation rate of the entire human mtGenome using a non-phylogenetic approach, we utilised 254 members of the most recent four generations of the Norfolk Island (NI) core pedigree ^16–20 and compared entire mtGenome sequences of maternal relatives from 45 maternal lineages. NI is a small, remote island in the Pacific with a unique admixture of paternal European ancestry in combination with Polynesian maternal origins ²¹. The population has a well-documented history and genealogy. Accurate and detailed historical accounts have been used by genealogists to create and maintain a well-documented, large multigenerational NI pedigree which consists of 5,742 individuals spanning 11 generations, and 200 years to the original founders ²². Due to the many inbreeding loops found in the population’s early generations, and the size and complexity, the NI pedigree was reconstructed to include all core individuals relating back to the original founders (n = 1388) ²⁰.

This study defines one of the first pedigree derived mutation rates encompassing the entire mtGenome. The expansion of the current paradigms and evidence of mutation rates allows researchers to more fully understand the processes that have shaped the evolution of the mtGenome. Our findings have application and implications for various research fields including clinical genetics, human evolution and forensic identity testing.

Ethics

Ethical clearance for the Norfolk Island mitochondrial DNA analysis portion of this study was provided originally by the Griffith University Human Research Ethics Committee (Approval MSC/04/09/HREC). Ethical clearance was transferred to and is now provided by the Queensland University of Technology Human Research Ethics Committee (Approval Number: 1400000749). No other ethical clearance was required. All Norfolk participants provided informed consent for research involvement

Sample Selection

This research used individuals from the NI Health Study and the associated NI core pedigree for research investigations at Queensland University of Technology. The NI Health Study has been well described in previous research ^16–20. The NI core pedigree contains individuals that are most closely related to the original founders. For illustration, the NI core pedigree is shown in Fig. 1. Using the NI core pedigree, maternal pedigrees were generated by establishing a list of founding mothers and tracing their maternal lines. Resultant pedigrees were generated using the Pedigree ²³, Kinship2 ²⁴ and Tidyverse ²⁵ packages. All pedigree analysis was conducted in RStudio ²⁶. The code utilised for this study is available in the GitHub repository http://sirselim.github.io/presentations/mt_tracing.html. In total, 45 resultant pedigrees (families) were chosen (Fig. 2). From these families, 225 individuals (including 125 females and 100 males) were chosen for sequencing, corresponding to 345 mtDNA transmissions. All individuals chosen were from the last four generations of the NI core pedigree and relate back to the original founders.

Whole mtGenome Sequencing

The entire mtGenome was amplified using long-range PCR with two overlapping primer sets (fragment 1: mt.569–9819, F: 5’-AACCAAACCCCAAAGACACC-3’ and R: 5’-GCCAATAATGACGTGAAGTCC-3’; and fragment 2: mt.9611–626, F: 5’-TCCCACTCCTAAACACATCC-3’ and R: 5’-TTTATGGGGTGATGTGAGCC-3’). Library preparation and sequencing was performed using an Ion Torrent™ high throughout sequencing protocol established in-house ²⁷.

Data Analysis

Data analysis was performed using the bioinformatics pipeline outlined in Harvey et al. ²⁷. Sequences were obtained for the entire mtGenome and were aligned in relation to the revised Cambridge Reference Sequence (rCRS) ²⁸ using an online tool Mitomaster to annotate variants and call haplotype information for each sample ²⁹.

Length variants that are known hotspots for insertion and deletions (indels) were ignored, including positions 309, 455, 463, 573, 960, 5899, 8276, 8285 and 16193 relative to the rCRS ³⁰. Variants in positions 16180 to 16183 were also ignored as reliable sequence data could not be obtained (see results and discussion).

In cases of lengthy heteroplasmy, a single dominant variant was identified. In cases of point heteroplasmy, the position was deemed heteroplasmic if the following requirements were achieved:

1. For NGS data, the minor allele frequency was greater than 20% (of total coverage). This threshold was chosen to reliably differentiate between signal and noise and reduce false positives.

2. When sequenced using Sanger Sequencing, the minor allele was readily visible (a distinct peak of normal morphology was evident and white space beneath the peak could be observed without changing the chromatogram view to examine the signal closer to the baseline) ¹².

3. The minor allele was evident in two different high-quality Sanger sequences (for example, when using both forward and reverse primers).

All mutations were confirmed using Sanger Sequencing. Amplification was performed in a 25µL total reaction volume using 1.6µL magnesium chloride, 5µL 5X GoTaq® Colourless Flexi Reaction Buffer, 0.5µL of 10mM dNTPs, 9.4µL deionised water, 2µL GoTaq® Hot Start Polymerase, 2µL of each 10µM amplification primer and 50ng DNA extract. Thermal cycling conditions were 95°C for 30 seconds; 30 cycles of 95°C for 30 seconds, 57°C for 40 seconds, 72°C for 30 seconds; and a 68°C hold for 3 minutes. Purification of amplified products prior to sequencing was performed using ExoSAP-IT™ PCR Product Clean-up Reagent, using 5µL PCR product, 4µL deionised water, and 1µL ExoSAP-IT™. Thermal cycling conditions were as per the manufacturer’s instructions. Sequencing was performed via the BigDye™ Terminator v3.1 Cycle Sequencing Kit. Reactions consisted of 9.64µL deionized water; 4µL BigDye™ Terminator v3.1 Sequencing Buffer; 0.5µL BigDye™ Terminator v3.1 Ready Reaction Mix; 3.2µL sequencing primer at 10µM; and 2.66µL PCR product for a total reaction volume of 20µL. Thermal cycling conditions used were: 96°C hold for 1 minute; 30 cycles of 96°C for 10 seconds, 50°C for 5 seconds, and 60°C for 4 minutes; followed by 1 cycle of 4°C for 5 minutes and 10°C for 5 minutes. Sequence product purification was performed via ethanol precipitation. Sequence detection was performed by capillary electrophoresis on a 3500 Genetic Analyser using a 50cm array, the FastSeq instrument protocol with default instrument settings. Amplification and sequencing primers are outlined in Supplementary Table S1.

The mutation rate was derived from the number of detected mutations per number of ‘meioses’ or transmission events, which is the number of cumulative generations tracing back to the maternal ancestor. The mutation rate is expressed as mutations per base pair (site) per million years (mutations/site/Myr), where the generational time was assumed to be 25 years. Different values of the mutation rate were estimated according to different assumptions: (1) all mutations that were detected (including heteroplasmy) were considered for the mutation rate calculation, and (2) all mutations that were detected (excluding heteroplasmy) were considered for the mutation rate calculation.

Confidence intervals (CIs: 95%) were calculated using Epitools, an online tool provided by AusVet Animal Health Services ³¹. The program outputs intervals using five alternative calculation methods as described in Brown et al. ³². The Wilson and Clopper Pearson methods were reported.

Whole mtGenome sequencing to identify mtDNA variants in the Norfolk Island pedigree

The whole mtGenome was sequenced for 225 individuals representing 45 independent mtDNA lineages from the core pedigree of the NI population isolate on an Ion Torrent™ NGS platform (Fig. 1). Sequence quality (Phred) scores remained consistent at > 25 for all samples at the median read length (140bp), with data reaching a median depth of approximately 370X for all 225 samples. Sequencing reads were aligned in relation to the rCRS and variants annotated using Mitomaster, with called variants subsequently validated by Sanger Sequencing. In the rCRS, position T16189 is flanked by cytosine residues between positions 16183 and 16194. In the samples examined in this study, this sequence stretch was consistently and reproducibly represented in NGS and Sanger data analysis when no length or point variation was present. However, in a number of samples, a combination of the T16189C transition and a deletion at position 16189 resulted in low coverage values and unreliable Sanger data for this region. Consequently, it was difficult to determine the exact number of cytosine residues present. Similar results were obtained for the HVII C-track region between positions 302 to 310 and 310 to 316. In this instance, difficulties in sequence confirmation were observed when there was an insertion of one or more cytosine (resulting in a stretch of ≥8 cytosine residues). Mutations in these regions were omitted during the comparison of maternally related individuals, and hence were not included in the mutation rate calculations performed in this study.

Calculation of mutation rates across the mtGenome

When the individual pedigrees were examined, the number of meioses and mutations were observed (Table 1). Homoplasmic mutations were found at two different sites in one individual each. Heteroplasmy was confirmed at 7 positions across 13 individuals (n = 225, 5.78% individuals). Analysis of the individuals belonging to 45 families from NI revealed 13.3% of the families contained at least one individual with mtDNA heteroplasmy as a result of mutation across the entire mtGenome. The frequency of point heteroplasmy in the NI samples did not show significant differences in relation to gender (Chi-squared test: X² = 0.472, df = 1, P = 0.492), nor was it found to be associated with a specific mtDNA haplogroup (Fisher exact test: P = 0.080).

Table 1

Mutations observed in 45 maternal lineages
Pedigree ID	Number of meioses in the pedigree	Number of individuals analysed	Number of mutations observed	Number of individuals with mutation	Mutations
13	14	10	2	1, 3	A8470R, A16280R
23	10	7	0	-	-
167	3	3	0	-	-
208	7	7	1	3	T9012Y
219	4	3	0	-	-
242	13	5	0	-	-
490	8	4	1	1	T146Y
691	130	59	2	2, 1	A2833R, A8817G
694	63	29	0	-	-
833	2	2	0	-	-
1088	3	4	0	-	-
1597	15	7	1	1	A16247G
1878	10	7	0	-	-
2071	5	4	0	-	-
2212	2	2	0	-	-
2328	2	2	0	-	-
2462	4	2	0	-	-
2802	3	3	0	-	-
2913	2	2	0	-	-
3310	4	3	0	-	-
3363	2	2	0	-	-
3963	3	3	0	-	-
5581	2	2	0	-	-
110390	1	2	0	-	-
124570	1	2	1	1	C16344Y
125750	2	3	0	-	-
169040	1	2	0	-	-
169810	1	2	0	-	-
218290	3	2	0	-	-
301005	1	2	0	-	-
310881	1	2	0	-	-
311021	2	3	0	-	-
313181	1	2	0	-	-
313211	1	2	0	-	-
313851	1	2	0	-	-
315560	1	2	0	-	-
317611	1	2	0	-	-
317671	2	3	0	-	-
317741	1	2	0	-	-
319011	1	2	0	-	-
319271	2	3	0	-	-
320321	2	3	0	-	-
321870	2	3	0	-	-
329180	4	5	1	2	C16320Y
400099	2	2	0	-	-
Total	345	225	9
Heteroplasmy are designated with an International Union of Biochemistry code ³³.

Further analysis of the 345 mtDNA transmissions identified 9 mutations (2 transitions and 7 cases of heteroplasmy) across the entire mtGenome, suggesting that 1.57x10^-6mutations occur (at a detectable level) in the mtGenome in each generation. Assuming a generation time of 25 years, the mutation rate for the entire mtGenome would be 0.0630 mutations/site/Myr (95% CI: 0.033 – 0.118) (Table 2). In order to make an estimate of separate mutation rates for the HVI and HVII regions, HVI was delimited by positions 16024 – 16383 (360 bp) and HVII delimited by positions 57 – 371 (315 bp) in accordance with Sigurgardottir et al. ².

Table 2 Mutation Rates for the NI extended pedigree

		95% CI Clopper Pearson		95% CI Wilson
	Point Estimate	Lower	Upper	Lower	Upper
mtGenome (1-16596)	0.063	0.029	0.118	0.033	0.118
HVI (16024–16383)	1.292	0.357	3.276	0.501	3.276
HVII (57–371)	0.369	0.013	2.038	0.064	2.064
HVI/HVII (16024–16383 and 57–371)	0.861	0.279	1.991	0.368	1.991
Control Region (1-576 and 16024-16596)	0.505	0.164	1.167	0.216	1.167
Coding Region (576 – 16024)	0.030	0.008	0.076	0.012	0.076

Pedigree derived mutation rate calculated for various regions of the mtGenome. Rates are expressed in mutations per site per million years (25 year generation time). mtGenome – mitochondrial genome. HVI – hypervariable region 1. HVII – hypervariable region II. All the mutations (including heteroplasmy) that were detected were considered

With the introduction of NGS, rapid sequencing of the entire mtGenome has become widely accessible for applications such as medical diagnostics, disease aetiology, and in forensic science with aged and degraded remains ^34–36. As NGS technologies are further developed, the benefits of high throughput, high quality, and low-cost methodologies continues to rise. With the advent of NGS it has also become feasible to sequence the entire mtGenome in large pedigrees and population cohorts. In this study we utilised NGS mtGenome data representing four generations of the NI pedigree to calculate mutation rates across the entire mtGenome.

Homopolymeric C-tracks in HVI (Position 16183 to 16194) and HVII (Position 302 to 316)

In this study, NGS analysis of homopolymer stretches (such as poly C regions) were a challenge for the Ion Torrent™ platform resulting in sequencing slippage in regions containing long homopolymer stretches. This is due to the degree of change in voltage loss resolution (above 6 to 8 bp) inherent to the platform. As a result, these regions have poor mapping quality and can lead to false-positive and false-negative calls and typically require assessment using alternative methods. The identification and calling INDELS also presented a bioinformatic challenge in the samples due to sequencing error bias and required a higher level of manual assessment. Issues in sequencing data in homopolymeric regions are not unique to NGS. Sanger Sequencing is also susceptible to such errors ³⁷, with all sequencing platforms prone to sequencing errors corresponding to the chemistry and platforms used ³⁸.

In 134 samples (n = 225, 59.6%), a combination of the T16189C transition and a deletion at position 16189 resulted in low coverage values and unreliable Sanger data for this region. Consequently, it became difficult to determine the number of cytosine residues present. In accordance with Interpretation Guidelines for mtDNA Sequencing outlined by the Laboratory Division of the FBI ³⁹, no attempt was made to count the number of residues for interpretation purposes and all comparisons assumed the same number were present. Similar results were obtained for the HVII C-track between position 302 to 310 and 310 to 316. Difficulties were also observed when an insertion of one or more Cs (resulting in a stretch of ≥8 C residues) were present. In an attempt to reduce the number of changes identified as false mutations, these regions were omitted during each comparison of maternally related individuals.

Heteroplasmy

Several previous studies have demonstrated Sanger Sequencing to be valid for quantification of heteroplasmy greater than 10% and that NGS could detect and quantify heteroplasmy as low as 1% ⁴⁰. In this study, to ensure robustness of the data and subsequent analysis, a threshold of 20% or greater was used. Previous studies have also performed independent DNA extraction, PCR amplification and sequencing to authenticate heteroplasmy results and exclude any contamination ^4,41,42. Here, heteroplasmy was confirmed via Sanger Sequencing.

Point heteroplasmy was confirmed at 7 positions across 13 samples (n = 225, 5.8% individuals), a rate found to be lower than other studies that have previously sequenced the entire mtGenome. Ramos et al. ⁴¹ observed point heteroplasmy in 12.8% of samples (13 of 101 individuals), while Santos et al. ⁴ observed point heteroplasmy in 12.1% of samples (28 of 232 individuals). Neither of these studies used NGS and both methodologies relied upon the same method to quantify levels of heteroplasmy. Specifically, these studies measured the height of the two peaks directly from the sequencing electropherogram and calculated the proportion between the heights of each peak with respect to the sum of the height of the two peaks. For comparison purposes, the obtained numerical proportions were then averaged. Another method previously used included cloning PCR products encompassing the D-loop region and subsequent sequencing of the 26–66 clones to determine the number containing each mtDNA variant ^4,41. The increased percentage of heteroplasmy observed in those studies when compared to ours could be due to the lower minor variant threshold used. The threshold used in Ramos et al. ⁴¹ was 10%, and although no threshold was stated in Santos et al. ⁴, that study included cases of heteroplasmy where the minor variant was as low as 2.5%. In addition to those outlined in Table 1, heteroplasmic mixtures below the 20% threshold were observed in three additional individuals, with MAFs of 12, 15 and 15% for A16280R, T9012Y and C16320, respectively. As these were below the 20% threshold, they were not considered in any further analysis. Further heteroplasmic variants may be evident at other positions across the mtGenome if our threshold was reduced to one comparable with Ramos et al. ⁴¹ or Santos et al. ⁴. In our study, 6.7% (n = 45) of the NI pedigrees examined, at least one individual within the lineage presented mtDNA heteroplasmy produced by mutations across the coding region. This value is consistent with the percentage reported by Santos et al. ⁴³ (6.5% of families). Interestingly, no individuals from the NI pedigrees examined in this study showed heteroplasmy at the positions examined by Santos et al. ⁴³.

Length heteroplasmy is typically observed in every individual where a transition at position 16,189 results in a homopolymer of nine or more cytosine residues, with no length heteroplasmy observed when seven or fewer cytosine residues are present ^12,44−46. In the NI samples, consistent and reliable Sanger Sequencing was not achieved for individuals that exhibited the 16,189 transition and therefore, length heteroplasmy could not be confirmed. Hence, to ensure consistent reporting, the decision was made to disregard variants in positions 16180 to 16183 for all samples. The poly-C tracks of both the HVI and HVII regions are known to have high insertion/deletion rates, resulting in length heteroplasmy ⁴⁷. For example, Santos et al. ⁴ reported length heteroplasmy produced by the insertion of cytosine residues in the poly-C tract of HVRI and HVRII respectively in 22.92% and 54.16% of the families studied. It is accepted that the general mechanism for generating length heteroplasmy in these regions is replication slippage ⁴⁸. These regions were also excluded during the analysis of the NI individuals and no other length heteroplasmy in any region was observed in the NI samples.

Mutation Rate

In the analysis of 345 genetic transmissions, 9 mutations (2 transitions and 7 cases of heteroplasmy) were detected across the entire mtGenome, suggesting the mutation rate for the entire mtGenome to be 0.063 mutations/site/Myr (95% CI: 0.033–0.118). Table 3 compares the derived mutation rate for each region calculated using the NI samples with previous studies using closed populations. For consistency, mutation rates for each study were converted to mutations/site/Myr with a defined generation time of 25 years. Note that although rates are grouped based on region (such as control region, coding region), the sequence range varies from study to study. Furthermore, inclusion requirements varied for each study impacting the rate calculated (for example, only substitutions for which there was evidence of a germinal origin were considered in Santos et al. ⁴, or all substitutions detected were considered in Madrigal et al. ⁴⁹). As expected, the mutation rate for the coding region was found to be smaller than that reported for the control region (0.030 mutations/site/Myr for coding region verses 0.505 mutations/site/Myr for control region), and as with other studies, the estimation reported here is much higher than those obtained by phylogenetic methods (Table 3).

In the NI samples, with the exception of the coding region, the mutation rate calculated for various regions of the mtGenome were found to be relatively high when compared to previous published studies. This may be attributed to differences in the population examined, study design or analysis. In this study, the rates calculated for the NI population included all mutations (heteroplasmy, germline, somatic). In contrast, for the HVI region, Santos et al. ⁴ reported a mutation rate of 0.074 mutations/site/Myr (17x lower than the rate for NI: 1.292 mutations/site/Myr), however their rate only included substitutions (including heteroplasmy) with a germinal origin present in women fixed at the individual level. While inclusion criteria alone does not fully explain the variations identified via pedigree-derived mutation rates, it likely contributes to the observed difference.

Interestingly, the mutation rate for the coding region was one of the smallest rates reported, despite our study being the only one to examine the entire region (bases 576 to 16024). This variation could be due to population differences. NI is a small, isolated population, and therefore the evolution of extremely high mutation rates is unlikely to occur unless organisms are under special circumstances ⁵⁰, with beneficial mutations rarely compensating for deleterious mutations.

In addition, heteroplasmic variants were not included in all studies assessing the mutation rate of mtDNA (for example, ¹). It is reasonable to suggest that heteroplasmy may resolve in favour of a ‘new’ base, and therefore, exclusion of heteroplasmic variants may underestimate the true mutation rate. However, the inclusion of heteroplasmic mutations introduces the limitation that only those that have reached a level of 20% of the mtDNA population are deemed detected. This raises the question of whether the discrepancy between phylogenetic and pedigree derived rates could be answered by simply excluding heteroplasmic variants, on the presumption that they will not become homoplasmic and/or will not be transmitted at the level of the population. Heteroplasmic mutations that reached 20%, the minimum threshold required to be included in the pedigree derived mutation rate reported here, should have a higher probability of becoming homoplasmic than those reaching only 1–2%.

Table 3

Summary of the derived mtDNA mutation rate for various published studies
Region	Study	Type	Sequence Range	mtDNA transmissions	Sub^a	Hetero^b	Total^c	Reported Mutation Rate	Standardised Mutation Rate*
HVI	⁴	Pedigree	16024–16383	321				1.84 x 10^− 6 mut/site/generation^d	0.074^d
	⁵¹	Phylogenetic	16024–16401					10.3x10^-8 mut/site/year^m	0.103
	¹	Pedigree	Amplification range: L15997 and H202	292	0			0.45 mut/site/Myr^f,g	0.360^f
	²	Pedigree	16024–16383	705	3	3	6	-	0.948^e
	Current	Pedigree	16024–16383	345	1	3	4	-	1.292^e
	⁴	Pedigree	16024–16383	321	0	6	6	-	2.083^e
	⁴⁹	Pedigree	16024–16383	299			7	71.1 x 10^− 6 mut/site/generation^h	2.844^h
HVII	⁵¹	Phylogenetic	29–408					7.39x10^-8 mut/site/year^m	0.074
	²	Pedigree	57–371	705	1	0	1	-	0.181^f
	¹	Pedigree	Amplification range: L16483 and H580	291	0	0	0	0.42 mut/site/Myr^f,g	0.336^f
	Current	Pedigree	57–371	345	0	1	1	-	0.369^e
	⁴	Pedigree	57–371	321				10.83 x 10^− 6 mut/site/generation^d	0.433^d
	⁴	Pedigree	57–371	321	0	5	5	-	1.984^e
HVI/HVII	¹	Pedigree	Amplification range: L16483 and H580	291	0	0	0	0.21mut/site/Myr^f,g	0.168^f
	⁴	Pedigree	57–371 and 16024–16383	321				6.04 x 10^− 6 mut/site/generation^d	0.241^d
	²	Pedigree	57–371 and 16024–16383	705				0.32 mut/site/Myr^f,g	0.256^f
	²	Pedigree	57–371 and 16024–16383	705	3	3	6	-	0.506^e
	Current	Pedigree	57–371 and 16024–16383	345	1	4	5	-	0.861^e
	⁵²	Pedigree	Undefined	327			10	2.5/site/Myr^e,g	2.000^e
	⁴	Pedigree	57–371 and 16024–16383	321	0	11	11	-	2.037^e
Control Region	⁵³	Pedigree	Undefined					1.5 x 10^− 6 mut/site/generation	0.060
	⁵¹	Phylogenetic	Undefined					7.00x10^-8 mut/site/year^m	0.070
	⁴	Pedigree	1–400 and 16024–16569	321				4.19 x 10^− 6 mut/site/generation^d	0.168^d
	¹	Pedigree	Amplification Range: L15997 and H202, with L16483 and H580	292	0	0	0	0.21 mut/site/Myr^f,g	0.168^f
	⁵⁴	Pedigree	1–576 and 16024–16596	185			1	0.24 mut/site/Myr^g,i	0.188ⁱ
	Current	Pedigree	1-576 and 16024–16569	345	1	4	5	-	0.505^e
	⁴	Pedigree	1–400 and 16024–16569	321				1.28 x 10^− 5 mut/site/generation^j	0.512^j
	⁴	Pedigree	1–400 and 16024–16596	321	0	6	6	1.92 x 10^− 5 mut/site/generationⁱ	0.768ⁱ
	⁴	Pedigree	1–400 and 16024–16569	321	0	11	11	3.52 x 10^− 5 mut/site/generation^e	1.194^e
Coding Region	⁵⁵	Phylogenetic	577–16023					1.25 ± 0.68x10^-8 sub/site/year^o	0.0125 ± 0.0068
	⁵⁶	Phylogenetic	Undefined					1.70x10^− 8 sub/site/yearⁿ	0.017
	⁴³	Pedigree	3230–4331	311				5.89 x 10^− 7 mut/site/generation^d	0.024^d
	Current	Pedigree	576–16024	345	1	3	4	-	0.030^e
	⁴³	Pedigree	3230–4331	311				1.03 x 10^− 6 mut/site/generation^l	0.041^l
	⁴³	Pedigree	3230–4331	311				5.84 x 10^− 6 mut/site/generationⁱ	0.234ⁱ
	⁴³	Pedigree	3230–4331	311				5.84 x 10^− 6 mut/site/generation^j	0.234^j
	⁴³	Pedigree	3230–4331	311	0	3	3	8.75 x 10^− 6 mut/site/generation^e	0.350^e
	⁴³	Pedigree	3230–4331	311				8.75 x 10^− 6 mut/site/generation^k	0.350^k
	¹	Pedigree	5550–6550	256	0	0	0	0.54 mut/site/Myr^f,g	0.432
Complete mtGenome	⁵³	Pedigree	1–16569					2.7 x 10^− 7 mut/site/generation^o	0.011
	⁵⁵	Phylogenetic						1.92x10^-8 sub/site/year	0.0192
	Current	Pedigree	1–16569	345	2	7	9	-	0.063^e
mut/site/Myr – mutations per nucleotide per million years. Mutation rates have been adjusted for consistent comparison.

Sub – substitutions

Hetero – heteroplasmic mutations

(a) Number of homoplasmic mutations observed

(b) Number of heteroplasmic mutations observed

(d) Only the substitutions (including heteroplasmy) with a germinal origin present in women that would become fixed at the individual level were considered

(e) All the substitutions (including heteroplasmy) that were detected were considered

(f) Only homoplasmic mutations were considered

(g) One generation is 20 years

(h) Unclear if heteroplasmy included in rate

(i) Only the substitutions (including for which there was evidence of a germinal origin were considered

(j) Only the substitutions (including heteroplasmy) present in women for whom there was evidence of a germinal origin were considered

(k) Only the substitutions (including heteroplasmy) present in women were considered

(l) Only the substitutions (including heteroplasmy) present in women that would become fixed at the individual level considering neutrality were considered

(m) Date of divergence human-chimpanzee used to calibrate evolutionary rate: 4.9 million years

(n) Date of divergence human-chimpanzee used to calibrate evolutionary rate: 5 million years

(o) Date of divergence human-chimpanzee used to calibrate evolutionary rate: 6.5 million year

The mtDNA control region is commonly used to assess human evolution and population movements. Many of these studies rely on phylogenetic analysis of mtDNA control region haplotype trees and phylogenetically derived rates of divergence. Estimates of the mutation rate with a non-phylogenetic approach (namely pedigree analysis) are reported to be approximately ten-fold higher than phylogenetically derived rates. Using this approach, 345 genetic transmissions, 9 mutations (2 transitions and 7 cases of heteroplasmy) were identified across the entire mtGenome, resulting in a mutation rate (obtained by employing the same definition of mutation used by other authors) of 0.0630 mutations/site/Myr (95% CI: 0.033–0.118). The mutation rate produced from this study is one of the first to use extended pedigree analysis of the entire mtGenome in combination with NGS. If mutations arise randomly in the mtGenome, these results suggest that newly arising mutations in the mitochondrial coding region are eliminated before they reach a detectable frequency. Defining this mechanism may provide further insight into the interpretation of mtGenome data across numerous fields of study.

Data availability

The data that support the findings of this study have restricted access due to ethical requirements and agreements with the Norfolk Island community, and so are not publicly available. Aggregate data may, however, be made available from the authors upon reasonable request. The Norfolk genetics steering committee will assess data access requests via our GRC computational genetics group (interested researchers should contact [email protected]).

Code availability

The code utilised for this study is available in the GitHub repository https://github.com/GRC-CompGen/mitochondrial_seq_pipeline

and http://sirselim.github.io/presentations/mt_tracing.html

Acknowledgements

The Norfolk Island mitochondrial analysis was supported by NHMRC APP1083450. This research was also supported by infrastructure purchased with Australian Government EIF Super Science Funds as part of the Therapeutic Innovation Australia - Queensland Node project (LRG). In addition, Jasmine Connell was the recipient of an Institute of Health and Biomedical Innovation (IHBI) Queensland University of Technology (QUT) postgraduate student scholarship and Miles Benton was supported by a Corbett Postgraduate Research Grant. The authors wish to acknowledge the Norfolk Island individuals who have generously donated their DNA samples and time to participate in this research study and for the ongoing community support which aids our research.

Author contribution statement

JRC - Conceptualization, formal analysis, investigation, methodology, validation, writing – original draft, review and editing; MCB – methodology, writing – review & editing; RAL – methodology, analysis review, writing – review & editing; HGS – sample preparation, quality control, writing – review & editing; JC - conceptualization, supervision, writing – review & editing; LMH - conceptualization, project administration, supervision, writing – review & editing; KMW - conceptualization, project administration, resources, supervision, writing – review & editing; LRG - conceptualization, project administration, resources, supervision, writing – review & editing

Conflicts of interest/Competing interests

The authors declare that they have no competing interests.

Statement on adherence to relevant ethical guidelines and informed consent obtained while conducting the experiment.

This study was conducted in accordance with the Australian National Statement on Ethical conduct in Human Research, with approval of the Queensland University of Technology Human Research Ethics Committee (approval number 1600000464). All participants were over 18 and provided informed consent.

1. Cavelier, L., Jazin, E., Jalonen, P. & Gyllensten, U. MtDNA substitution rate and segregation of heteroplasmy in coding and noncoding regions. Hum. Genet.107, 45–50 (2000).

2. Sigurğardóttir, S., Helgason, A., Gulcher, J. R., Stefansson, K. & Donnelly, P. The mutation rate in the human mtDNA control region. Am. J. Hum. Genet.66, 1599–1609 (2000).

3. Henn, B. M., Gignoux, C. R., Feldman, M. & Mountain, J. L. Characterizing the time dependency of human mitochondrial DNA mutation rate estimates. Mol. Biol. Evol.26, 217–230 (2009).

4. Santos, C. et al. Understanding differences between phylogenetic and pedigree-derived mtDNA mutation rate: A model using families from the Azores Islands (Portugal). Mol. Biol. Evol.22, 1490–1505 (2005).

5. Macaulay, V. et al. mtDNA mutation rates - No need to panic. Am. J. Hum. Genet.61, 983–990 (1997).

6. Pääbo, S. Mutational hot spots in the mitochondrial microcosm. Am. J. Hum. Genet.59, 493–496 (1996).

7. Kareem, M. A., Abdulzahra, A. I., Hameed, I. H. & Jebor, M. A. A new polymorphic positions discovered in mitochondrial DNA hypervariable region HVIII from central and north-central of Iraq. Mitochondrial DNA Part A27, 3250–3254 (2016).

8. Melton, T., Holland, C. & Holland, M. Forensic mitochondrial DNA analysis: Current practice and future potential. Forensic Sci. Rev.24, 101–122 (2012).

9. Margulies, M. et al. Genome sequencing in microfabricated high-density picolitre reactors. Nature437, 376 (2005).

10. Bodner, M. et al. Helena, the hidden beauty: Resolving the most common West Eurasian mtDNA control region haplotype by massively parallel sequencing an Italian population sample. Forensic Sci. Int. Genet.15, 21–26 (2015).

11. Coble, M. et al. Single nucleotide polymorphisms over the entire mtDNA genome that increase the power of forensic testing in Caucasians. Int. J. Legal Med.118, 137–146 (2004).

12. Just, R. et al. Full mtGenome reference data: Development and characterization of 588 forensic-quality haplotypes representing three U.S. populations. Forensic Sci. Int. Genet.14, 141–155 (2015).

13. King, J. L. et al. High-quality and high-throughput massively parallel sequencing of the human mitochondrial genome using the Illumina MiSeq. Forensic Sci. Int. Genet.12, 128–135 (2014).

14. Parsons, T. & Coble, M. Increasing the forensic discrimination of mitochondrial DNA testing through analysis of the entire mitochondrial DNA genome. Croat. Med. J.42, 304–309 (2001).

15. Stoljarova, M., King, J. L., Takahashi, M., Aaspõllu, A. & Budowle, B. Whole mitochondrial genome genetic diversity in an Estonian population sample. Int. J. Legal Med.130, 67–71 (2016).

16. Bellis, C. et al. Phenotypical Characterisation of the Isolated Norfolk Island Population Focusing on Epidemiological Indicators of Cardiovascular Disease. Hum. Hered.60, 211–219 (2005).

17. Bellis, C. et al. Linkage disequilibrium analysis in the genetically isolated Norfolk Island population. Heredity100, 366–373 (2008).

18. Cox, H. C. et al. Principal Component and Linkage Analysis of Cardiovascular Risk Traits in the Norfolk Isolate. Hum. Hered.68, 55–64 (2009).

19. Benton, M. C. et al. Mapping eQTLs in the Norfolk Island Genetic Isolate Identifies Candidate Genes for CVD Risk Traits. Am. J. Hum. Genet.93, 1087–1099 (2013).

20. Benton, M. et al. ‘Mutiny on the Bounty’: the genetic history of Norfolk Island reveals extreme gender-biased admixture. Investig. Genet.6, (2015).

21. Mcevoy, B. P. et al. European and Polynesian admixture in the Norfolk Island population. Hered. Cardiff105, 229–34 (2010).

22. Macgregor, S. et al. Legacy of mutiny on the Bounty: founder effect and admixture on Norfolk Island. Eur. J. Hum. Genet. EJHG18, 67–72 (2010).

23. Coster, A. pedigree. (2013).

24. Therneau, T. M., Daniel, S., Sinnwell, J. & Atkinson, E. kinship2. (2015).

25. Wickham, H. & RStudio. tidyverse. (2017).

26. RStudio. (2017).

27. Harvey, N. R. et al. Ion torrent high throughput mitochondrial genome sequencing (HTMGS). PLOS ONE14, e0224847 (2019).

28. Andrews, R. M. et al. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat. Genet.23, 147 (1999).

29. MITOMAP. MITOMAP A human mitochondrial genome database. http://www.mitomap.org/MITOMAP (2017).

30. Parson, W. EMPOP mtDNA Database Directions for Use. (2019).

31. AusVet Animal Health Services. Epi Tools - Calculate confidence limits for a sample proportion. http://epitools.ausvet.com.au/content.php?page=CIProportion (2016).

32. Brown, L. D., Cai, T. T. & DasGupta, A. Interval estimation for a binomial proportion. Stat. Sci.16, 101–133 (2001).

33. Nomenclature Committee of the International Union of Biochemistry. Nomenclature for incompletely specified bases in nucleic acid sequences. Recommendations 1984. Biochem. J.229, 281–286 (1985).

34. Schirmer, M. et al. Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform. Nucleic Acids Res.43, e37–e37 (2015).

35. Templeton, J. E. L. et al. DNA capture and next-generation sequencing can recover whole mitochondrial genomes from highly degraded samples for human identification. Investig. Genet.4, 26 (2013).

36. Yang, Y., Xie, B. & Yan, J. Application of next generation sequencing technology in forensic science. Genomics Proteomics Bioinformatics12, 190–197 (2014).

37. Yohe, S. & Thyagarajan, B. Review of Clinical Next-Generation Sequencing. Arch. Pathol. Lab. Med.141, 1544–1557 (2017).

38. Ratan, A. et al. Comparison of Sequencing Platforms for Single Nucleotide Variant Calls in a Human Sample. PLOS ONE8, e55089 (2013).

39. Budowle, B., DiZinno, J. A. & Wilson, M. R. Interpretation guidelines for mitochondrial DNA sequencing. in 10th International Symposium on Human Identification (10th International Symposium on Human Identification, 1999).

40. Kloss-Brandstätter, A. et al. Validation of next-generation sequencing of entire mitochondrial genomes and the diversity of mitochondrial DNA mutations in oral squamous cell carcinoma. PLoS ONE10, (2015).

41. Ramos, A. et al. Frequency and Pattern of Heteroplasmy in the Complete Human Mitochondrial Genome. PLOS ONE8, e74636 (2013).

42. Santos, C. et al. Frequency and pattern of heteroplasmy in the control region of human mitochondrial DNA. J. Mol. Evol.67, 191–200 (2008).

43. Santos, C. et al. Mutation patterns of mtDNA: Empirical inferences for the coding region. BMC Evol. Biol.8, 167 (2008).

44. Melton, T. Mitochondrial DNA Heteroplasmy. Forensic Sci. Rev.16, 1–20 (2004).

45. Lee, H. Y., Chung, U., Yoo, J.-E., Park, M. J. & Shin, K.-J. Quantitative and qualitative profiling of mitochondrial DNA length heteroplasmy. Electrophoresis25, 28–34 (2004).

46. Irwin, J. et al. Investigation of heteroplasmy in the human mitochondrial DNA control region: A synthesis of observations from more than 5000 global population samples. J. Mol. Evol.68, 516–527 (2009).

47. Hauswirth, W. W. & Clayton, D. A. Length heterogeneity of a conserved displacement-loop sequence in human mitochondrial DNA. Nucleic Acids Res.13, 8093–8104 (1985).

48. Hauswirth, W., Van De Walle, M., Olivo, P. & Laipis, P. Heterogeneous mitochondrial DNA D-loop sequences in bovine tissue. Cell37, 1001–1007 (1984).

49. Madrigal, L. et al. High mitochondrial mutation rates estimated from deep-rooting Costa Rican pedigrees. Am. J. Phys. Anthropol.148, 327–333 (2012).

50. Sniegowski, P. D., Gerrish, P. J., Johnson, T. & Shaver, A. The evolution of mutation rates: separating causes from consequences. BioEssays22, 1057–1066 (2000).

51. Horai, S., Hayasaka, K., Kondo, R., Tsugane, K. & Takahata, N. Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs. Proc. Natl. Acad. Sci.92, 532–536 (1995).

52. Parsons, T. et al. A high observed substitution rate in the human mitochondrial DNA control region. Nat. Genet.15, 363–368 (1997).

53. Rebolledo-Jaramillo, B. et al. Maternal age effect and severe germ-line bottleneck in the inheritance of human mitochondrial DNA. Proc. Natl. Acad. Sci.111, 15474–15479 (2014).

54. Howell, N. et al. The pedigree rate of sequence divergence in the human mitochondrial genome: There is a difference between phylogenetic and pedigree rates. Am. J. Hum. Genet.72, 659–670 (2003).

55. Fu, Q. et al. A revised timescale for human evolution based on ancient mitochondrial genomes. Curr. Biol. CB23, 553–559 (2013).

56. Ingman, M., Kaessmann, H., Pääbo, S. & Gyllensten, U. Mitochondrial genome variation and the origin of modern humans. Nature408, 708–713 (2000).

No competing interests reported.

SupplementaryTableS.docx
Supplementary Table S1

Download PDF

Editorial decision: Major revision
01 Dec, 2021
Reviews received at journal
29 Nov, 2021
Reviews received at journal
26 Sep, 2021
Reviewers agreed at journal
08 Sep, 2021
Reviewers invited by journal
02 Sep, 2021
Editor assigned by journal
25 Aug, 2021
Editor invited by journal
25 Aug, 2021
Submission checks completed at journal
25 Aug, 2021
First submitted to journal
18 Aug, 2021

You are reading this latest preprint version

Pedigree Derived Mutation Rate Across the Entire Mitochondrial Genome of the Norfolk Island Population

Status:

Version 1

Abstract

Figures

Introduction

Materials And Methods

Ethics

Sample Selection

Whole mtGenome Sequencing

Data Analysis

Results

Whole mtGenome sequencing to identify mtDNA variants in the Norfolk Island pedigree

Calculation of mutation rates across the mtGenome

Discussion

Homopolymeric C-tracks in HVI (Position 16183 to 16194) and HVII (Position 302 to 316)

Heteroplasmy

Mutation Rate

Declarations

Data availability

Code availability

Author contribution statement

Conflicts of interest/Competing interests

References

Additional Declarations

Supplementary Files

Status:

Version 1