Genetic Diversity in Mitochondria

Genetic Diversity in Mitochondria

We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

As I was studying for my GRE exam, I was reading on Mitochondria. I then remembered that they undergo replication similar to bacteria. So I googled for info and found on Wikipedia that there is a " near-absence of genetic recombination in mitochondrial DNA"

So then, if there is a lack of genetic diversity, how can these mitochondria still be functional? If we assume they are affected by random mutations as DNA can be. Is it the case that they replicate infrequently enough to not suffer ill-effects of low diversity?

Your claim that "there is a lack of genetic diversity" in mitochondria is not correct. There is no connection between recombination and mutation rates per se. In fact, mitochondria have much higher mutation rates than nucleus and different types of mitochondrial DNA can even co-exist in one organism, a phenomenon called heteroplasmy (which is found at a frequency of ~90% in human).

But: there is a connection between mutation fixation rate and effective population size. Because mitochondria are normally inherited through maternal line, effective population size for their genome is smaller than that of the nuclear genome, and thus any mutations in the mitochondrial DNA reach fixation faster, but again this is counteracted by high mutation rates.

There is no problem with functionality: as always if something gets broken (and it is easy to break something in the mitochondrial genome, since genes are normally tightly packed in it) respective genotype is eliminated or reduced in frequency by natural selection. In the case of mitochondria this selection acts both upon the level of individual organisms (intra-individual somatic pool of mitochondria: positive and negative selection) and the level of populations (mostly purifying selection in humans). It is nevertheless true, that given the virtual lack of recombination, selection becomes the only force of purging deleterious mutations.

Population biology and genetic diversity of two adjacent shrimp ( Parapenaeopsis coromandelica ) populations exploited under different fishing pressures in the coastal waters of Sri Lanka

Parapenaeopsis coromandelica shrimp populations along the western coast of Sri Lanka have supported coastal trawling for the last hundreds of years. Two non-overlapping adjacent fishing grounds (Hendala and Negombo) sustain different fishing intensities. In order to obtain information on the population structure and genetic diversity of P. coromandelica in these two regions, differences in length–weight relationships, growth, spawning seasons, sex-ratios, gonadosomatic index, length at 50% maturity (L 50 ) and sequence variation of the mitochondrial cytochrome oxidase subunit I were examined. Significant differences in population biology and genetic diversity were revealed from the two fishing grounds. Samples of Hendala showed lower L 50 and genetic diversity which are considered as potential effects and symptoms of extensive selective harvesting. Further, the shrimps' behaviour seems to be triggering the separation through low mixing of individuals at the two fishing grounds resulting in significant divergence based on haplotype frequencies. Management of P. coromandelica should consider the revealed biological and genetic evidences on existence of two sub-populations/stocks together with a routine monitoring of genetic effects due to harvesting.

Wallace, D. C. A mitochondrial bioenergetic etiology of disease. J. Clin. Invest. 123, 1405–1412 (2013).

Pham, T. D. et al. Cristae remodeling causes acidification detected by integrated graphene sensor during mitochondrial outer membrane permeabilization. Sci. Rep. 6, 35907 (2016).

Walker, J. E. The ATP synthase: the understood, the uncertain and the unknown. Biochem. Soc. Trans. 41, 1–16 (2013).

Wallace, D. C. & Fan, W. Energetics, epigenetics, mitochondrial genetics. Mitochondrion 10, 12–31 (2010).

Wallace, D. C., Fan, W. & Procaccio, V. Mitochondrial energetics and therapeutics. Annu. Rev. Pathol. 5, 297–348 (2010).

Picard, M. et al. Progressive increase in mtDNA 3243A>G heteroplasmy causes abrupt transcriptional reprogramming. Proc. Natl. Acad. Sci. USA 111, E4033–E4042 (2014).

Couvillion, M. T., Soto, I. C., Shipkovenska, G. & Churchman, L. S. Synchronized mitochondrial and cytosolic translation programs. Nature 533, 499–503 (2016).

Giles, R. E., Blanc, H., Cann, H. M. & Wallace, D. C. Maternal inheritance of human mitochondrial DNA. Proc. Natl. Acad. Sci. USA 77, 6715–6719 (1980).

Pyle, A. et al. Extreme-depth re-sequencing of mitochondrial DNA finds no evidence of paternal transmission in humans. PLoS Genet. 11, e1005040 (2015).

Brown, W. M., Prager, E. M., Wang, A. & Wilson, A. C. Mitochondrial DNA sequences of primates: tempo and mode of evolution. J. Mol. Evol. 18, 225–239 (1982).

Neckelmann, N., Li, K., Wade, R. P., Shuster, R. & Wallace, D. C. cDNA sequence of a human skeletal muscle ADP/ATP translocator: lack of a leader peptide, divergence from a fibroblast translocator cDNA, and coevolution with mitochondrial DNA genes. Proc. Natl. Acad. Sci. USA 84, 7580–7584 (1987).

Wallace, D. C. et al. Sequence analysis of cDNAs for the human and bovine ATP synthase beta subunit: mitochondrial DNA genes sustain seventeen times more mutations. Curr. Genet. 12, 81–90 (1987).

Wallace, D. C. Mitochondrial DNA variation in human radiation and disease. Cell 163, 33–38 (2015).

Sorrentino, V., Menzies, K. J. & Auwerx, J. Repairing mitochondrial dysfunction in disease. Annu. Rev. Pharmacol. Toxicol. 58, 353–389 (2018).

Gorman, G. S. et al. Prevalence of nuclear and mitochondrial DNA mutations related to adult mitochondrial disease. Ann. Neurol. 77, 753–759 (2015).

Lightowlers, R. N., Taylor, R. W. & Turnbull, D. M. Mutations causing mitochondrial disease: what is new and what challenges remain? Science 349, 1494–1499 (2015).

Elliott, H. R., Samuels, D. C., Eden, J. A., Relton, C. L. & Chinnery, P. F. Pathogenic mitochondrial DNA mutations are common in the general population. Am. J. Hum. Genet. 83, 254–260 (2008).

Chinnery, P. F., Elliott, H. R., Hudson, G., Samuels, D. C. & Relton, C. L. Epigenetics, epidemiology and mitochondrial DNA diseases. Int. J. Epidemiol. 41, 177–187 (2012).

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

Wallace, D. C. et al. Mitochondrial DNA mutation associated with Leber’s hereditary optic neuropathy. Science 242, 1427–1430 (1988).

Holt, I. J., Harding, A. E., Petty, R. K. & Morgan-Hughes, J. A. A new mitochondrial disease associated with mitochondrial DNA heteroplasmy. Am. J. Hum. Genet. 46, 428–433 (1990).

Wallace, D. C. et al. Familial mitochondrial encephalomyopathy (MERRF): genetic, pathophysiological, and biochemical characterization of a mitochondrial DNA disease. Cell 55, 601–610 (1988).

Shoffner, J. M. et al. Myoclonic epilepsy and ragged-red fiber disease (MERRF) is associated with a mitochondrial DNA tRNA (Lys) mutation. Cell 61, 931–937 (1990).

Goto, Y., Nonaka, I. & Horai, S. A mutation in the tRNA (Leu)(UUR) gene associated with the MELAS subgroup of mitochondrial encephalomyopathies. Nature 348, 651–653 (1990).

Huoponen, K., Vilkki, J., Aula, P., Nikoskelainen, E. K. & Savontaus, M. L. A new mtDNA mutation associated with Leber hereditary optic neuroretinopathy. Am. J. Hum. Genet. 48, 1147–1153 (1991).

Johns, D. R., Neufeld, M. J. & Park, R. D. An ND-6 mitochondrial DNA mutation associated with Leber hereditary optic neuropathy. Biochem. Biophys. Res. Commun. 187, 1551–1557 (1992).

Fischel-Ghodsian, N., Prezant, T. R., Bu, X. & Oztas, S. Mitochondrial ribosomal RNA gene mutation in a patient with sporadic aminoglycoside ototoxicity. Am. J. Otolaryngol. 14, 399–403 (1993).

Hutchin, T. et al. A molecular basis for human hypersensitivity to aminoglycoside antibiotics. Nucleic Acids Res. 21, 4174–4179 (1993).

Pfeffer, G. & Chinnery, P. F. Diagnosis and treatment of mitochondrial myopathies. Ann. Med. 45, 4–16 (2013).

Andreu, A. L. et al. A nonsense mutation (G15059A) in the cytochrome b gene in a patient with exercise intolerance and myoglobinuria. Ann. Neurol. 45, 127–130 (1999).

Andreu, A. L. et al. Exercise intolerance due to mutations in the cytochrome b gene of mitochondrial DNA. N. Engl. J. Med. 341, 1037–1044 (1999).

Andreu, A. L., Checcarelli, N., Iwata, S., Shanske, S. & DiMauro, S. A missense mutation in the mitochondrial cytochrome b gene in a revisited case with histiocytoid cardiomyopathy. Pediatr. Res. 48, 311–314 (2000).

Mancuso, M. et al. Mitochondrial myopathy and complex III deficiency in a patient with a new stop-codon mutation (G339X) in the cytochrome b gene. J. Neurol. Sci. 209, 61–63 (2003).

DiMauro, S. Mitochondrial encephalomyopathies—at fifty years on: the Robert Wartenberg lecture. Neurology 81, 281–291 (2013).

Emmanuele, V. et al. A novel mutation in the mitochondrial DNA cytochrome b gene (MTCYB) in a patient with mitochondrial encephalomyopathy, lactic acidosis, and strokelike episodes syndrome. J. Child Neurol. 28, 236–242 (2013).

Wallace, D.C., Lott, M.T. & Procaccio, V. Mitochondrial medicine: the mitochondrial biology and genetics of metabolic and degenerative diseases, cancer, and aging. in Emery and Rimoin’s Principles and Practice of Medical Genetics Vol. 1 (eds. Rimoin, D. L. et al.) Ch. 11, 1–153 (Churchill Livingstone Elsevier, Philadelphia, 2013).

Wallace, D.C., Lott, M.T. & Procaccio, V. Mitochondrial biology—mitochondrial medicine. in Principles and Practice of Medical Genetics and Genomics (eds. Pyeritz, R. E. et al.) (Elsevier, Philadelphia, 2018).

MITOMAP. A Human Mitochondrial Genome Database. (2018).

Jenuth, J. P., Peterson, A. C., Fu, K. & Shoubridge, E. A. Random genetic drift in the female germline explains the rapid segregation of mammalian mitochondrial DNA. Nat. Genet. 14, 146–151 (1996).

Wai, T., Teoli, D. & Shoubridge, E. A. The mitochondrial DNA genetic bottleneck results from replication of a subpopulation of genomes. Nat. Genet. 40, 1484–1488 (2008).

Cree, L. M. et al. A reduction of mitochondrial DNA molecules during embryogenesis explains the rapid segregation of genotypes. Nat. Genet. 40, 249–254 (2008).

Li, M. et al. Transmission of human mtDNA heteroplasmy in the genome of the Netherlands families: support for a variable-size bottleneck. Genome Res. 26, 417–426 (2016).

Wilson, I. J. et al. Mitochondrial DNA sequence characteristics modulate the size of the genetic bottleneck. Hum. Mol. Genet. 25, 1031–1041 (2016).

Carling, P. J., Cree, L. M. & Chinnery, P. F. The implications of mitochondrial DNA copy number regulation during embryogenesis. Mitochondrion 11, 686–692 (2011).

Jokinen, R. & Battersby, B. J. Insight into mammalian mitochondrial DNA segregation. Ann. Med. 45, 149–155 (2013).

Floros, V. I. et al. Segregation of mitochondrial DNA heteroplasmy through a developmental genetic bottleneck in human embryos. Nat. Cell Biol. 20, 144–151 (2018).

Fan, W. et al. A mouse model of mitochondrial disease reveals germline selection against severe mtDNA mutations. Science 319, 958–962 (2008).

Stewart, J. B. et al. Strong purifying selection in transmission of mammalian mitochondrial DNA. PLoS Biol. 6, e10 (2008).

Sharpley, M. S. et al. Heteroplasmy of mouse mtDNA is genetically unstable and results in altered behavior and cognition. Cell 151, 333–343 (2012).

Shoffner, J. M. et al. Mitochondrial DNA variants observed in Alzheimer disease and Parkinson disease patients. Genomics 17, 171–184 (1993).

Hutchin, T. & Cortopassi, G. A mitochondrial DNA clone is associated with increased risk for Alzheimer disease. Proc. Natl. Acad. Sci. USA 92, 6892–6895 (1995).

Santoro, A. et al. Evidence for sub-haplogroup h5 of mitochondrial DNA as a risk factor for late onset Alzheimer’s disease. PLoS One 5, e12037 (2010).

Chagnon, P. et al. Phylogenetic analysis of the mitochondrial genome indicates significant differences between patients with Alzheimer disease and controls in a French-Canadian founder population. Am. J. Med. Genet. 85, 20–30 (1999).

Carrieri, G. et al. Mitochondrial DNA haplogroups and APOE4 allele are non-independent variables in sporadic Alzheimer’s disease. Hum. Genet. 108, 194–198 (2001).

van der Walt, J. M. et al. Analysis of European mitochondrial haplogroups with Alzheimer disease risk. Neurosci. Lett. 365, 28–32 (2004).

Lakatos, A. et al. Association between mitochondrial DNA variations and Alzheimer’s disease in the ADNI cohort. Neurobiol. Aging 31, 1355–1363 (2010).

van der Walt, J. M. et al. Mitochondrial polymorphisms significantly reduce the risk of Parkinson disease. Am. J. Hum. Genet. 72, 804–811 (2003).

Ghezzi, D. et al. Mitochondrial DNA haplogroup K is associated with a lower risk of Parkinson’s disease in Italians. Eur. J. Hum. Genet. 13, 748–752 (2005).

Khusnutdinova, E. et al. A mitochondrial etiology of neurodegenerative diseases: evidence from Parkinson’s disease. Ann. NY Acad. Sci. 1147, 1–20 (2008).

Fachal, L. et al. No evidence of association between common European mitochondrial DNA variants in Alzheimer, Parkinson, and migraine in the Spanish population. Am. J. Med. Genet. B. Neuropsychiatr. Genet. 168B, 54–65 (2015).

Hudson, G. et al. Two-stage association study and meta-analysis of mitochondrial DNA variants in Parkinson disease. Neurology 80, 2042–2048 (2013).

Chalkia, D. et al. Association between mitochondrial DNA haplogroup variation and autism spectrum disorders. JAMA Psychiatry 74, 1161–1168 (2017).

Wang, Y., Picard, M. & Gu, Z. Genetic evidence for elevated pathogenicity of mitochondrial DNA heteroplasmy in Autism Spectrum Disorder. PLoS Genet. 12, e1006391 (2016).

Fuku, N. et al. Mitochondrial haplogroup N9a confers resistance against type 2 diabetes in Asians. Am. J. Hum. Genet. 80, 407–415 (2007).

Nishigaki, Y. et al. Mitochondrial haplogroup N9b is protective against myocardial infarction in Japanese males. Hum. Genet. 120, 827–836 (2007).

Chalkia, D. et al. Mitochondrial DNA associations with East Asian metabolic syndrome. Biochim. Biophys. Acta 1859, 878–892 (2018).

Latorre-Pellicer, A. et al. Mitochondrial and nuclear DNA matching shapes metabolism and healthy ageing. Nature 535, 561–565 (2016).

Yu, X. et al. Dissecting the effects of mtDNA variations on complex traits using mouse conplastic strains. Genome Res. 19, 159–165 (2009).

Wallace, D. C. Genetics: mitochondrial DNA in evolution and disease. Nature 535, 498–500 (2016).

Wallace, D. C. & Chalkia, D. Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease. Cold Spring Harb. Perspect. Biol. 5, a021220 (2013).

Craven, L. et al. Pronuclear transfer in human embryos to prevent transmission of mitochondrial DNA disease. Nature 465, 82–85 (2010).

Hyslop, L. A. et al. Towards clinical application of pronuclear transfer to prevent mitochondrial DNA disease. Nature 534, 383–386 (2016).

Tachibana, M. et al. Mitochondrial gene replacement in primate offspring and embryonic stem cells. Nature 461, 367–372 (2009).

Tachibana, M. et al. Towards germline gene therapy of inherited mitochondrial diseases. Nature 493, 627–631 (2013).

Kang, E. et al. Mitochondrial replacement in human oocytes carrying pathogenic mitochondrial DNA mutations. Nature 540, 270–275 (2016).

Paull, D. et al. Nuclear genome transfer in human oocytes eliminates mitochondrial DNA variants. Nature 493, 632–637 (2013).

Wolf, D. P., Hayama, T. & Mitalipov, S. Mitochondrial genome inheritance and replacement in the human germline. EMBO J. 36, 2177–2181 (2017).

Ma, H. et al. Functional human oocytes generated by transfer of polar body genomes. Cell Stem Cell 20, 112–119 (2017).

Slone, J., Zhang, J. & Huang, T. Experience from the first live-birth derived from oocyte nuclear transfer as a treatment strategy for mitochondrial diseases. J. Mol. Genet. Med. 11, 1000258 (2017).

Yamada, M. et al. Genetic drift can compromise mitochondrial replacement by nuclear transfer in human oocytes. Cell Stem Cell 18, 749–754 (2016).

Barritt, J. A., Brenner, C. A., Malter, H. E. & Cohen, J. Mitochondria in human offspring derived from ooplasmic transplantation. Hum. Reprod. 16, 513–516 (2001).

Zhang, J. et al. Live birth derived from oocyte spindle transfer to prevent mitochondrial disease. Reprod. Biomed. Online 34, 361–368 (2017).

Kujoth, G. C. et al. Mitochondrial DNA mutations, oxidative stress, and apoptosis in mammalian aging. Science 309, 481–484 (2005).

Trifunovic, A. et al. Premature ageing in mice expressing defective mitochondrial DNA polymerase. Nature 429, 417–423 (2004).

Vermulst, M. et al. DNA deletions and clonal mutations drive premature aging in mitochondrial mutator mice. Nat. Genet. 40, 392–394 (2008).

Kang, E. et al. Age-related accumulation of somatic mitochondrial DNA mutations in adult-derived human iPSCs. Cell Stem Cell 18, 625–636 (2016).

Schriner, S. E. et al. Extension of murine life span by overexpression of catalase targeted to mitochondria. Science 308, 1909–1911 (2005).

Lee, H. Y. et al. Targeted expression of catalase to mitochondria prevents age-associated reductions in mitochondrial function and insulin resistance. Cell Metab. 12, 668–674 (2010).

Woo, D. K. et al. Mitochondrial genome instability and ROS enhance intestinal tumorigenesis in APC(Min/+) mice. Am. J. Pathol. 180, 24–31 (2012).

Arnheim, N. & Cortopassi, G. Deleterious mitochondrial DNA mutations accumulate in aging human tissues. Mutat. Res. 275, 157–167 (1992).

Soong, N. W., Hinton, D. R., Cortopassi, G. & Arnheim, N. Mosaicism for a specific somatic mitochondrial DNA mutation in adult human brain. Nat. Genet. 2, 318–323 (1992).

Corral-Debrinski, M. et al. Mitochondrial DNA deletions in human brain: regional variability and increase with advanced age. Nat. Genet. 2, 324–329 (1992).

Corral-Debrinski, M. et al. Marked changes in mitochondrial DNA deletion levels in Alzheimer brains. Genomics 23, 471–476 (1994).

Coskun, P. E., Beal, M. F. & Wallace, D. C. Alzheimer’s brains harbor somatic mtDNA control-region mutations that suppress mitochondrial transcription and replication. Proc. Natl. Acad. Sci. USA 101, 10726–10731 (2004).

Coskun, P. E. et al. Systemic mitochondrial dysfunction and the etiology of Alzheimer’s disease and down syndrome dementia. J. Alzheimers Dis. 20(Suppl 2), S293–S310 (2010).

Coskun, P. et al. A mitochondrial etiology of Alzheimer and Parkinson disease. Biochim. Biophys. Acta 1820, 553–564 (2012).

Horton, T. M. et al. Marked increase in mitochondrial DNA deletion levels in the cerebral cortex of Huntington’s disease patients. Neurology 45, 1879–1883 (1995).

Corral-Debrinski, M. et al. Hypoxemia is associated with mitochondrial DNA damage and gene induction. Implications for cardiac disease. J. Am. Med. Assoc. 266, 1812–1816 (1991).

Ji, F. et al. Mitochondrial DNA variant associated with Leber hereditary optic neuropathy and high-altitude Tibetans. Proc. Natl. Acad. Sci. USA 109, 7391–7396 (2012).

Zheng, S. et al. Role of mtDNA haplogroups in COPD susceptibility in a southwestern Han Chinese population. Free Radic. Biol. Med. 53, 473–481 (2012).

Sherer, T. B. et al. Mechanism of toxicity of pesticides acting at complex I: relevance to environmental etiologies of Parkinson’s disease. J. Neurochem. 100, 1469–1479 (2007).

Cannon, J. R. & Greenamyre, J. T. Neurotoxic in vivo models of Parkinson’s disease recent advances. Prog. Brain Res. 184, 17–33 (2010).

Bulstrode, H. et al. Mitochondrial DNA and traumatic brain injury. Ann. Neurol. 75, 186–195 (2014).

Kilbaugh, T. J. et al. Peripheral blood mitochondrial DNA as a biomarker of cerebral mitochondrial dysfunction following traumatic brain injury in a porcine model. PLoS One 10, e0130927 (2015).

Baudouin, S. V. et al. Mitochondrial DNA and survival after sepsis: a prospective study. Lancet 366, 2118–2121 (2005).

Gomez, R. et al. Association of mitochondrial allele 4216C with increased risk for complicated sepsis and death after traumatic injury. J. Trauma 66, 850–857 (2009). discussion 857–858.

Hendrickson, S. L. et al. Mitochondrial DNA haplogroups influence. AIDS progression. AIDS 22, 2429–2439 (2008).

Zhang, Q. et al. Circulating mitochondrial DAMPs cause inflammatory responses to injury. Nature 464, 104–107 (2010).

Zhang, Q., Itagaki, K. & Hauser, C. J. Mitochondrial DNA is released by shock and activates neutrophils via p38 map kinase. Shock 34, 55–59 (2010).

Krysko, D. V. et al. Emerging role of damage-associated molecular patterns derived from mitochondria in inflammation. Trends Immunol. 32, 157–164 (2011).

Dan Dunn, J., Alvarez, L. A., Zhang, X. & Soldati, T. Reactive oxygen species and mitochondria: A nexus of cellular homeostasis. Redox Biol. 6, 472–485 (2015).

Nadeau-Vallée, M. et al. A critical role of interleukin-1 in preterm labor. Cytokine Growth Factor Rev. 28, 37–51 (2016).

Angelin, A. et al. Foxp3 reprograms T cell metabolism to function in low-glucose, high-lactate environments. Cell Metab. 25, 1282–1293.e7 (2017).

Ravi, S., Mitchell, T., Kramer, P., Chacko, B. & Darley-Usmar, V. M. Mitochondria in monocytes and macrophages-implications for translational and basic research. Int. J. Biochem. Cell Biol. 53, 202–207 (2014).


Whole genome resequencing of Japanese wild silkworm

The pupa of a single individual of Japanese wild silkworm was used to extract genomic DNA, based on a standard protocol for genomic DNA extraction. We followed the manufacturer's instructions (Illumina) to prepare DNA library. Briefly, we made it by following the workflow on Illumina Genome Analyzer II system: cluster generation, template hybridization, isothermal amplification, linearization, blocking, and denaturization and hybridization of the sequencing primers. Then the Illumina base-calling pipeline (SolexaPipeline-0.3) was applied to get sequences from the fluorescent images. Short read data have been deposited in the NCBI Short Read Archive under the accession number SRA009886.

Public data availability

Through NCBI Short Read Archive with accession number SRA009208, we downloaded all of the 40 silkworm resequencing data (see Table ​ Table1), 1 ), including 29 domesticated silkworms and 11 Chinese wild ones, detailed information of which can be found elsewhere [7], the silkworm nuclear reference genome sequence was derived from the SilkDB [11] ( or From GenBank we obtained three mt genome sequences (accession numbers <"type":"entrez-nucleotide","attrs":<"text":"AB070264","term_id":"18640062","term_text":"AB070264">> AB070264, <"type":"entrez-nucleotide","attrs":<"text":"AY301620","term_id":"60391177","term_text":"AY301620">> AY301620, and <"type":"entrez-nucleotide","attrs":<"text":"NC_003395","term_id":"18644896","term_text":"NC_003395">> NC_003395) and their annotations, which were serving as references. We used these reference genomes in all of our consensus assembly below.

Table 1

The detailed information of samples and sequencing summary

Sample IDName of strainSystem or locationEffective Depth (X)Genome coverage (%)DifferenceRate (%)
D7Lao IILaos12.8897.421.32
D8India M3India51.299.220.44
D10Italy16Italy, Europe24.1299.460.35
D11Soviet Union NO.1Former SU, Europe10.9894.750.41
D15Mutation M2Mutation22.0399.621.32
D16A06EGuangdong province, China75.0499.711.52
D17DamaoSichuan province, China28.6699.110.64
D18Ankang NO.4Shanxi province, China79.9399.910.62
D19ZT500Gansu province, China165.4299.990.61
D20ZhuguiZhejiang province, China77.6199.770.43
D21BilianJiangsu province, China62.8399.491.4
D22ZT900Sichuan province, China80.8199.611.05
D23ZT100Hunan province, China89.799.830.41
D24SihongJiangsu province, China105.0399.750.31
D25XiaoshiwanZhejiang province, China20.4698.310.59
D26C108Chongqing, China27.5792.660.4
D27Sichuang M3Sichuan province, China51.399.670.58
D28ZT000Guizhou province, China166.8599.960.53
D29HandanHebei province, China51.9499.900.32
W1B. mandarina ZiyangSichuan province, China27.696.072.49
W2B. mandarina NanchongSichuan province, China106.1598.182.46
W3B. mandarina HongyaSichuan province, China25.2794.181.94
W4B. mandarina PengshanSichuan province, China76.0395.451.95
W5B. mandarina AnkangShanxi province, China30.5492.302.01
W6B. mandarina YichangHubei province, China108.9896.632.77
W7B. mandarina YanchengJiangsu province, China72.0697.102.03
W8B. mandarina LuzhouSichuan province, China40.5492.452.33
W9B. mandarina HunanHunan province, China114.5796.912.03
W10B. mandarina SuzhouJiangsu province, China49.1195.402.05
W11B. mandarina RongchangChongqing, China79.3195.892.04
W12B. mandarina JapanHokkaido, Japan13.5281.470.81

Reads alignment and consensus assembly

Although we only included the mt genome sequence in the analysis, we prefer to map raw short reads onto the whole genome [including nuclear DNA sequence and mitochondrial DNA (mtDNA) sequence], in order to make alignment more reliable. Due to the fact that three sorts of resequencing data of silkworm have been in hand - which are the domesticated group, the Chinese wild population, and the Japanese wild population - there must be three different references to be used by program SOAP v1.09 [13]. So we built three different reference genomes which are the silkworm nuclear reference genome sequence with C108 mtDNA sequence, with Chinese wild mtDNA sequence and with Japanese wild mtDNA sequence. The mapping results showed that Chinese wild silkworms have a higher mismatch rate than the domesticated ones (see Table ​ Table1), 1 ), which suggests a high sequence diversity in Chinese wild group, even in the intraspecific comparison between Chinese wild variety Ankang and the Chinese wild reference (which is also from Ankang in Shanxi Province of China). Followed by SOAPsnp [14], which based on Bayesian theory, we calculate the posterior probability of each possible genotype at every genome position, from the alignment results for each sample. Then the consensus was structured by the highest probability. 41 consensus sequences have been put into the Genbank under accession number <"type":"entrez-nucleotide","attrs":<"text":"GU966591","term_id":"291575480","term_text":"GU966591">> GU966591- <"type":"entrez-nucleotide","attrs":<"text":"GU966631","term_id":"291576040","term_text":"GU966631">> GU966631.

SNP detection and experimental validation

We used SOAPsnp [14] to call SNPs for each variety. After setting

six steps were used to filter out the unreliable variants: 1) we set Q20 quality cutoff 2) two unique reads were allowed at least 3) the SNPs had to be at least 5 bp away from each other 4) the approximate copy number of flanking sequences had to be no more than 2 5) P value of the rank sum test had to be no less than 0.05 6) the number of unique mapped reads had to be greater than or equal to half of the number of total mapped reads.

To evaluate our SNP calling strategy, we randomly selected 50 sites for PCR-Sanger dideoxy sequencing validation using the AB 3730XL. After manually checking all the intensity trace files, we found that all the sites were confirmed by the PCR-sequencing results.

Linkage Disequilibrium (LD) measure calculation

To measure LD level in the silkworm mt population, we used the normalized measure of allelic association estimate D' [15], which can not be easily influenced by rare alleles examined [16]. We set the parameters in the software Haploview [17] as follows:

-maxdistance 200-dprime-minGeno 0.6-minMAF 0.1-hwcutoff 0.001. Then spot chart was plotted with R scripts which drew averaged D' against pairwise marker distance.

Phylogeny reconstruction

Five silkworm mt sequences (with accession number <"type":"entrez-nucleotide","attrs":<"text":"AB070264","term_id":"18640062","term_text":"AB070264">> AB070264, <"type":"entrez-nucleotide","attrs":<"text":"AY301620","term_id":"60391177","term_text":"AY301620">> AY301620, <"type":"entrez-nucleotide","attrs":<"text":"NC_003395","term_id":"18644896","term_text":"NC_003395">> NC_003395, <"type":"entrez-nucleotide","attrs":<"text":"NC_004622","term_id":"162279939","term_text":"NC_004622">> NC_004622, and <"type":"entrez-nucleotide","attrs":<"text":"NC_012727","term_id":"238563960","term_text":"NC_012727">> NC_012727) were firstly aligned using MUSCLE v3.7 [18] with default settings, adjusting coordinates for the three resequencing groups. From this, we got 46 mt genomes, after integrating these five sequences and 41 consensus that we have, to perform phylogenetic reconstruction following MEGA v4 [19] by using Neighbour-Joining method and Mrbayes v3.1.2 [20] under the Bayesian theory. In Mrbayes, we chose the GTR+gamma+I model, and set the chain length to 50,000,000 (1 sample/1000 generations) and burned in the first 10,000 samples. Almost identical results were obtained in two independent runs. The quality of being dependable of the NJ trees was bootstrapped with 1000 replicate estimates. To evaluate the confidence level in the tree selection, we applied statistical tests with CONSEL [21]. The site-wise likelihood file is derived from PhyML calculatation [22].

Tests of effective population size

First, using the SNP information of each group, we inferred the mt sequence mismatch distributions for domesticated group and Chinese wild group, where a ragged distribution implies the stable population size, and where a bell-shaped pattern is often related to population growth [23,24]. Further, Tajima's D and Fu and Li's D tests for population size were also performed with DnaSP v5 [25].

Neutrality tests

It became clear to us that, during the domestication process of the silkworm, 354 genes in the nuclear genome bear a strong human selection footprint [7]. The mt genome may therefore also experience selective pressure. In order to examine mtDNA protein evolution, we compare the rate of nonsynonymous to synonymous mutation within and between species in each gene. If these ratios differ significantly, they provide evidence of selection pressure [26]. We then applied Williams' correction to calculate the G statistic [27]. Moreover, we also preformed the two likelihood ratio tests (LRTs), based on widely used branch-site models from PAML 4.2 [28], to detect any positive selection.

Materials and Methods

DNA Samples

A total of 74 samples were collected, including blood samples (16), hair samples (16), and dried leather samples (42). Due to degradation, DNA extractions were successful for only 21 of the 42 dried leather samples ( table 1 ). Therefore, the total number of DNA samples was reduced to 53. Both of the two subspecies were included, with five of them being Ailurus fulgens fulgens and the others being Ailurus fulgens styani ( table 1 ). The blood and hair samples were obtained from the Chongqing Zoo and Chengdu Zoos of China, and their wild origins were known. Blood samples were anticoagulated with heparin and stored at −70°C before DNA extraction. The hair samples were collected by plucking and stored at −70°C. The dried leather samples were obtained from collections of the Kunming Institute of Zoology, Chinese Academy of Sciences, and stored at −70°C after sampling. The 53 red pandas were originally from 8 different geographic locations in the Sichuan and Yunnan provinces of China ( fig. 1 ). Although efforts were made to avoid sampling related individuals, the relationships among animals in the sample were generally unknown.

DNA Extraction, Polymerase Chain Reaction, and Sequencing

DNA extractions from blood samples follow the standard phenol-chloroform method. The fresh hair and dried leather samples were first treated with proteinase K at 56°C for 2 h and then incubated with 10% Chelex 100 (Bio-Rad) at 98°C for 30 min. After centrifugation at a high speed (10,000 rpm) for 10 min, the supernatants were collected and directly used as DNA templates for PCR ( Walsh 1990 ). The PCR was conducted by predenaturing at 94°C for 2 min, cycling at 94°C for 1 min, 56°C for 1 min, and 72°C for 1 min for 35–40 cycles, and a final extension at 72°C for 5 min. The primer sequences are CAC CAT CAA CAC CCA AAG CTG (forward) and TTC ATG GGC CCG GAG CGA G (reverse), which amplify a 276-bp fragment located upstream of the mtDNA D-loop region. The PCR products were purified through low-melting-point agarose gel electrophoresis. Sequencing was conducted on an ABI377 automatic sequencer with both forward and reverse primers.

Phylogenetic Analysis and Statistical Tests of Neutrality

For phylogenetic analysis, parsimony (PAUP, version 3.1.1 Swofford 1993 ) and median-joining network analyses ( Bandelt, Forster, and Röhl 1999 ) were used. The homologous sequence of the raccoon (Procyon lotor), the closest living relative of the red panda, was included as an outgroup. The pairwise mismatch distribution was generated using Arlequin, version 2.000 (Schneider, Roessli, and Excoffier 2000). The essential population parameter θ was estimated using Watterson's (1975) estimate, Tajima's (1983) estimate, and Fu's (1994) UPBLUE estimate. Watterson's estimate is based on the number of segregating sites among the sequences. Tajima's estimate is based on the calculation of the mean number of pairwise differences of the sequences, while Fu's UPBLUE estimate is done by incorporating the genealogical information of the sequences. A statistical test of neutrality was carried out using Fu's (1997)FS test. Strictly speaking, all three of these estimators of θ are based on the infinite-sites model ( Watterson 1975 Tajima 1983 Fu 1997 ). Since the sequences generated in this study are from the D-loop region that has mutation hot spots, the infinite-sites model is violated to some extent. To minimize the effect of violation of the model on the estimation of θ, as well as statistical tests of neutrality, we inferred all the required information for parameter estimation and neutrality testing from the parsimony analysis. This was done by first reconstructing a parsimony tree from the sequences and then inferring the required information from the tree. For example, to infer the total number of mutations in the sample, we counted the total number of steps in the parsimony tree. For each pair of sequences, the distance needed for UPBLUE could easily be computed from the parsimony tree as well.

Fu's FS test of neutrality was used to infer the population history of the red panda. The FS value tends to be negative when there is an excess of recent mutations, and therefore a large negative value of FS will be taken as evidence against the neutrality of mutations, an indication of deviation caused by population growth and/or selection.

Jodi Nunnari

Jodi Nunnari received her BA from the College of Wooster and her PhD from Vanderbilt University. She began her studies of mitochondria as a post-doctoral fellow at the University of California, San Francisco in the lab of Peter Walter. Nunnari is currently a Professor of Molecular and Cellular Biology at the University of California, Davis… Continue Reading

More Talks in Cell Biology


Restriction fragment length polymorphisms in the D-loop,tRNA Met+Glu+Il and ND3 gene fragments

MtDNA D-loop, tRNA Met+Glu+Il and ND3 fragments from T739, TA2,615, BALB/C, C3H, C57BL/6J and DBA/2, were cleaved by 16 restriction endonucleases, including HaeIII, HinfI, EcoRV, HindIII, HpaI, BamHI, ApaI, NdeII, XhoI, XbaI, AluI, RsaI, StuI, DraI, AvaI and HaeII, respectively. No differences in restriction maps were observed, nor were there any variations in the 46 restriction sites in mtDNA D-loop,tRNA Met+Glu+Ile and ND3 fragments in these strains(Fig. 1, Table 2).

Restriction patterns of mouse mt DNA D-Loop, tRN Ile+GlN+Met and ND3 fragments

. D-loop (1100 bp) . . tRNA Ile+GlN+Met (1126 bp) . . ND3 (534 bp) . .
Enzyme . Sites . Fragment length . Sites . Fragment length . Sites . Fragment length .
HaeIII 4 456,445,122. 4 470,290,136. 1 392,142
HinfI 1 700,400 0 1126 4 253,84.
ApaI 1 655,445 1 658,469 0 534
AluI 4 545,269,102. 3 665,318,135. 1 310,224
BamHI 0 1100 2 710,252,164 0 534
DraI 0 1100 1 634,493 1 451,83
StuI 0 1100 2 605,297,225 1 391,143
NdeII 2 845,118,37 0 1126 0 534
RsaI 6 720,178,125. 4 392,259,228. 2 353,103,78
XbaI 1 680,420 0 1126 0 524
. D-loop (1100 bp) . . tRNA Ile+GlN+Met (1126 bp) . . ND3 (534 bp) . .
Enzyme . Sites . Fragment length . Sites . Fragment length . Sites . Fragment length .
HaeIII 4 456,445,122. 4 470,290,136. 1 392,142
HinfI 1 700,400 0 1126 4 253,84.
ApaI 1 655,445 1 658,469 0 534
AluI 4 545,269,102. 3 665,318,135. 1 310,224
BamHI 0 1100 2 710,252,164 0 534
DraI 0 1100 1 634,493 1 451,83
StuI 0 1100 2 605,297,225 1 391,143
NdeII 2 845,118,37 0 1126 0 534
RsaI 6 720,178,125. 4 392,259,228. 2 353,103,78
XbaI 1 680,420 0 1126 0 524

Note: There are no other restriction enzyme sites.

Single-stranded conformation polymorphisms in D-loop 5′ and 3′ end fragments

To further analyze the genetic variations in these strains, D-loop 5′and 3′ end fragments, both of high variability, were subjected to PCR-SSCP analysis. No differences in the SSCP electrophoresis bands were observed (Fig. 2).

Genetic Diversity in Mitochondria - Biology

All articles published by MDPI are made immediately available worldwide under an open access license. No special permission is required to reuse all or part of the article published by MDPI, including figures and tables. For articles published under an open access Creative Common CC BY license, any part of the article may be reused without permission provided that the original article is clearly cited.

Feature Papers represent the most advanced research with significant potential for high impact in the field. Feature Papers are submitted upon individual invitation or recommendation by the scientific editors and undergo peer review prior to publication.

The Feature Paper can be either an original research article, a substantial novel research study that often involves several techniques or approaches, or a comprehensive review paper with concise and precise updates on the latest progress in the field that systematically reviews the most exciting advances in scientific literature. This type of paper provides an outlook on future directions of research or possible applications.

Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world. Editors select a small number of articles recently published in the journal that they believe will be particularly interesting to authors, or important in this field. The aim is to provide a snapshot of some of the most exciting work published in the various research areas of the journal.

Genetic differences between people across the world are no greater than differences between pigeons

In terms of genetics, the difference between any two humans is about the same as the difference between two pigeons.

These results from a massive “DNA barcoding” study will come as a surprise to anyone who thinks of themselves as truly unique, or who places great store by differences between groups.

The typical genetic distinctions within a species, including humans, is one in 1,000 of the "letters" that make up a DNA sequence.

"Culture, life experience and other things can make people very different but in terms of basic biology, we're like the birds,” said Dr Mark Stoeckle from the University of Basel, who co-led the study.


"By determining the genetic variety within species of the animal kingdom, made possible only recently by the burgeoning number of DNA sequences, we've documented the absence of human exceptionalism."

Dr Stoeckle and his colleagues made use of a database containing over five million genetic barcodes from more than 100,000 animal species, assembled by scientists worldwide over the past 15 years.

They looked specifically at mitochondrial sequences – pieces of DNA found in the mitochondria, which are tiny structures providing cells with power.

"If a Martian landed on Earth and met a flock of pigeons and a crowd of humans, one would not seem more diverse than the other according to the basic measure of mitochondrial DNA," said Jesse Ausubel, director of the program for the human environment at The Rockefeller University.

In a paper published in the journal Human Evolution, the researchers concluded that in fact there is very little genetic diversity within most animal species, and very clear genetic distinctions between species.

Allio R, Donega S, Galtier N, Nabholz B (2017) Large variation in the ratio of mitochondrial to nuclear mutation rate across animals: implications for genetic diversity and the use of mitochondrial DNA as a molecular marker. Mol Biol Evol 34(11):2762–2772.

Ballard JWO, Whitlock MC (2004) The incomplete natural history of mitochondria. Mol Ecol 13(4):729–744.

Bazin E, Glémin S, Galtier N (2006) Population size does not influence mitochondrial genetic diversity in animals. Science 312(5773):570–572.

Blanchet S, Prunier JG, De Kort H (2017) Time to go bigger: emerging patterns in macrogenetics. Trends Genet 33(9):579–580.

Bohonak AJ, Vandergast AG (2011) The value of DNA sequence data for studying landscape genetics. Mol Ecol 20(12):2477–2479.

CBD and UNEP (2010) Strategic plan for biodiversity 2011–2020 and the Aichi targets

Currie DJ (1991) Energy and large-scale patterns of animal- and plant-species richness. Am Nat 137(1):27–49

DeBalsi KL, Hoff KE, Copeland WC (2017) Role of the mitochondrial DNA replication machinery in mitochondrial DNA mutagenesis, aging and age-related diseases. Ageing Res Rev 33:89–104.

DiBattista JD (2008) Patterns of genetic variation in anthropogenically impacted populations. Conserv Genet 9(1):141–156.

Edwards S, Bensch S (2009) Looking forwards or looking backwards in avian phylogeography? A comment on Zink and Barrowclough 2008. Mol Ecol 18:2930–2933

Falconer DS, Mackay TFC (1996) Introduction to quantitative genetics, 4th edn. Longman, Harlow

Fisher RA (1930) The genetical theory of natural selection. Oxford University Press, Oxford

Frankham R (1996) Relationship of genetic variation to population size in wildlife. Conserv Biol 10(6):1500–1508.

Galtier N, Jobson RW, Nabholz B, Glémin S, Blier PU (2009a) Mitochondrial whims: metabolic rate, longevity and the rate of molecular evolution. Biol Lett 5(3):413–416.

Galtier N, Nabholz B, Glémin S, Hurst GDD (2009b) Mitochondrial DNA as a marker of molecular diversity: a reappraisal. Mol Ecol 18(22):4541–4550.

Hoban S, Bruford M, D’Urban Jackson J, Lopes-Fernandes M, Heuertz M, Hohenlohe PA, Paz-Vinas I, Sjögren-Gulve P, Segelbacher G, Vernesi C, Aitken S, Bertola LD, Bloomer P, Breed M, Rodríguez-Correa H, Funk WC, Grueber CE, Hunter ME, Jaffe R, Liggins L, Mergeay J, Moharrek F, O’Brien D, Ogden R, Palma-Silva C, Pierson J, Ramakrishnan U, Simo-Droissart M, Tani N, Waits L, Laikre L (2020) Genetic diversity targets and indicators in the CBD post-2020 Global Biodiversity Framework must be improved. Biol Conserv 248:108654.

Hoffmann S, Spitkovsky D, Radicella JP, Epe B, Wiesner RJ (2004) Reactive oxygen species derived from the mitochondrial respiratory chain are not responsible for the basal levels of oxidative base modifications observed in nuclear DNA of mammalian cells. Free Radic Biol Med 36(6):765–773.

James J, Eyre-Walker A (2020) Mitochondrial DNA sequence diversity in mammals: a correlation between the effective and census population sizes. Genome Biol Evol.

Johnson MTJ, Munshi-South J (2017) Evolution of life in urban environments. Science 358(6363):eaam8327.

Lanfear R, Thomas JA, Welch JJ, Brey T, Bromham L (2007) Metabolic rate does not calibrate the molecular clock. Proc Natl Acad Sci USA 104(39):15388–15393.

Lawrence ER, Fraser DJ (2020) Latitudinal biodiversity gradients at three levels: linking species richness, population richness and genetic diversity. Glob Ecol Biogeogr.

Leigh DM, Hendry AP, Vázquez-Domínguez E, Friesen VL (2019) Estimated six per cent loss of genetic variation in wild populations since the industrial revolution. Evol Appl 12(8):1505–1512.

Li H, Xiang-Yu J, Dai G, Gu Z, Ming C, Yang Z, Ryder OA, Li WH, Fu YX, Zhang YP (2016) Large numbers of vertebrates began rapid population decline in the late 19th century. Proc Natl Acad Sci USA 113(49):14079–14084.

Loveless M, Hamrick J (1984) Ecological determinants of genetic structure in plant populations. Annu Rev Ecol Syst 51:65–95

Manel S, Guerin PE, Mouillot D, Blanchet S, Velez L, Albouy C, Pellissier L (2020) Global determinants of freshwater and marine fish genetic diversity. Nat Commun 11(1):1–9.

Millette KL, Fugère V, Debyser C, Greiner A, Chain FJJ, Gonzalez A (2020) No consistent effects of humans on animal genetic diversity worldwide. Ecol Lett 23(1):55–67.

Miraldo A, Li S, Borregaard MK, Florez-Rodriguez A, Gopalakrishnan S, Rizvanovic M, Wang Z, Rahbek C, Marske KA, Nogues-Bravo D (2016) An Anthropocene map of genetic diversity. Science 353(6307):1532–1535.

Mittelbach GG, Schemske DW, Cornell HV, Allen AP, Brown JM, Bush MB, Harrison SP, Hurlbert AH, Knowlton N, Lessios HA, McCain CM, McCune AR, McDade LA, McPeek MA, Near TJ, Price TD, Ricklefs RE, Roy K, Sax DF, Schluter D, Sobel JM, Turelli M (2007) Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. Ecol Lett 10(4):315–331.

Mittell EA, Nakagawa S, Hadfield JD (2015) Are molecular markers useful predictors of adaptive potential? Ecol Lett 18(8):772–778.

Nabholz B, Glémin S, Galtier N (2008a) Strong variations of mitochondrial mutation rate across mammals—the longevity hypothesis. Mol Biol Evol 25(1):120–130.

Nabholz B, Mauffrey JF, Bazin E, Galtier N, Glemin S (2008b) Determination of mitochondrial genetic diversity in mammals. Genetics 178(1):351–361.

Nevo E, Beiles A, Ben-Shlomo R (1984) The evolutionary significance of genetic diversity: ecological, demographic, and life history correlates. In: Mani GS (ed) Lecture Notes in Biomathematics vol 53, evolutionary dynamics of genetic diversity. Springer, Berlin, pp 13–213.

Pontarp M, Bunnefeld L, Cabral JS, Etienne RS, Fritz SA, Gillespie R, Graham CH, Hagen O, Hartig F, Huang S, Jansson R, Maliet O, Münkemüller T, Pellissier L, Rangel TF, Storch D, Wiegand T, Hurlbert AH (2019) The latitudinal diversity gradient: novel understanding through mechanistic eco-evolutionary models. Trends Ecol Evol 34(3):211–223.

Rubinoff D, Holland BS (2005) Between two extremes: mitochondrial DNA is neither the panacea nor the nemesis of phylogenetic and taxonomic inference. Syst Biol 54(6):952–961.

Schmidt C, Domaratzki M, Kinnunen RP, Bowman J, Garroway CJ (2020a) Continent-wide effects of urbanization on bird and mammal genetic diversity. Proc R Soc B Biol Sci 287(1920):20192497.

Schmidt C, Dray S, Garroway CJ (2020b) Genetic and species-level biodiversity patterns are linked by demography and ecological opportunity. bioRxiv.

Soulé ME (1976) In: Ayala FJ (ed) Allozyme variation, its determinants in space and time. Sinauer Associates, Sunderland, Massachusetts. pp. 60–77

Stearns SC (1992) The evolution of life histories. Oxford University Press, New York

Theodoridis S, Fordham DA, Brown SC, Li S, Rahbek C, Nogues-Bravo D (2020) Evolutionary history and past climate change shape the distribution of genetic diversity in terrestrial mammals. Nat Commun 11(1):2557.

Zink RM, Barrowclough GF (2008) Mitochondrial DNA under siege in avian phylogeography. Mol Ecol 17(9):2107–2121.


  1. Gasar

    I think you are not right. Write in PM, we will talk.

  2. Nabi Ulmalhamah

    Bravo, one sentence ... another idea

  3. Bondig

    you can look at that infinitely.

  4. Son

    It seems that I have already seen in another blog about this topic.

  5. Kirisar

    In my opinion you are not right. I am assured. Let's discuss.

Write a message