- Research article
- Open Access
Genetic diversity, molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes
© Li et al; licensee BioMed Central Ltd. 2010
- Received: 16 November 2009
- Accepted: 24 March 2010
- Published: 24 March 2010
Mitochondria are a valuable resource for studying the evolutionary process and deducing phylogeny. A few mitochondria genomes have been sequenced, but a comprehensive picture of the domestication event for silkworm mitochondria remains to be established. In this study, we integrate the extant data, and perform a whole genome resequencing of Japanese wild silkworm to obtain breakthrough results in silkworm mitochondrial (mt) population, and finally use these to deduce a more comprehensive phylogeny of the Bombycidae.
We identified 347 single nucleotide polymorphisms (SNPs) in the mt genome, but found no past recombination event to have occurred in the silkworm progenitor. A phylogeny inferred from these whole genome SNPs resulted in a well-classified tree, confirming that the domesticated silkworm, Bombyx mori, most recently diverged from the Chinese wild silkworm, rather than from the Japanese wild silkworm. We showed that the population sizes of the domesticated and Chinese wild silkworms both experience neither expansion nor contraction. We also discovered that one mt gene, named cytochrome b, shows a strong signal of positive selection in the domesticated clade. This gene is related to energy metabolism, and may have played an important role during silkworm domestication.
We present a comparative analysis on 41 mt genomes of B. mori and B. mandarina from China and Japan. With these, we obtain a much clearer picture of the evolution history of the silkworm. The data and analyses presented here aid our understanding of the silkworm in general, and provide a crucial insight into silkworm phylogeny.
- Effective Population Size
- Nuclear Genome
- Cytb Gene
- Domestication Event
- Haploid Chromosome Number
Mitochondria have a genetic system which is different and distinct from the nuclear genome. They have a rapid tempo of nucleotide substitution and present a special form of maternal inheritance . The domesticated silkworm, B. mori, is the model organism of Lepidoptera, which has been producing industrial silk for humans over a period of more than 5,000 years . It embraces more than 1,000 inbred strains worldwide , and has recently been serving as a bioreactor . There are two major populations of B. mandarina: notably the Chinese, with a haploid chromosome number of 28 (equal to those of B. mori), and the Japanese, with a haploid chromosome number of 27. They have different geographic distributions . In 2006 Arunkumar et al.  analyzed six complete mt genomes to construct the phylogeny for the Bombycidae family, and although the concurrence of their methods found the Chinese B. mandarina to be the progenitor species of B. mori, some of their methods (ML phylogeny) also tentatively proposed that Chinese and Japanese B. mandarina share a clade, belonging to a single, monophyletic group on the phylogenic tree, distinct from B. mori. Based on the same data, Pan et al.  concurred with the consensus of Arunkumar et al., that the Japanese B. mandarina is a more distant taxon from the domesticated silkworm than the Chinese B. mandarina.
A clearer picture for the Bombycidae phylogeny in confirmation of this became desirable. Here, we take advantage of 41 mt genomes (covering both wild silkworm groups and the domesticated silkworms) to again infer the ancestor of B. mori, and arrive at compelling evidence confirming the conclusions of that the apparent progenitor of B. mori is the Chinese B. mandarina, in confirmation of both Arunkumar et al and Pan et al.
The domesticated silkworm has evolved to be completely dependent on humans for survival, and many domestication-related genes have been found to be selected for, when comparing the domesticated silkworm with the Chinese wild silkworm, accounting for 2.9% of the genome-wide genes . This remarkable footprint by selection effect implied a potential selection imprint left in the silkworm mt genome. Moreover, previous studies have demonstrated that variants in mt genome were often subjected to selection [8–10]. With availability of genome-wide resequencing data from 40 silkworms, here we have the opportunity to make numerous comparisons and to investigate the phylogeny of silkworm mitochondria, in order to study whether the silkworm mt population has experienced any expansion or contraction, and to detect the selection effect shaping on the mt genome.
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
The detailed information of samples and sequencing summary
Name of strain
System or location
Effective Depth (X)
Genome coverage (%)
Soviet Union NO.1
Former SU, Europe
Guangdong province, China
Sichuan province, China
Shanxi province, China
Gansu province, China
Zhejiang province, China
Jiangsu province, China
Sichuan province, China
Hunan province, China
Jiangsu province, China
Zhejiang province, China
Sichuan province, China
Guizhou province, China
Hebei province, China
B. mandarina Ziyang
Sichuan province, China
B. mandarina Nanchong
Sichuan province, China
B. mandarina Hongya
Sichuan province, China
B. mandarina Pengshan
Sichuan province, China
B. mandarina Ankang
Shanxi province, China
B. mandarina Yichang
Hubei province, China
B. mandarina Yancheng
Jiangsu province, China
B. mandarina Luzhou
Sichuan province, China
B. mandarina Hunan
Hunan province, China
B. mandarina Suzhou
Jiangsu province, China
B. mandarina Rongchang
B. mandarina Japan
AB070264 represents the strain C108, which served as the domesticated population reference. AY301620 was used as the Chinese wild mitochondria reference, which is from Ankang in Shanxi Province of China . NC_003395 is a Japanese wild individual from Tsukuba, Ibaraki, Japan , and we map our resequencing data (Japanese wild silkworm) onto it. Mt genome sequence of Antheraea pernyi (NC_004622) and Eriogyna pyretorum (NC_012727) were also downloaded to be applied as outgroup when reconstructing phylogeny.
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 . 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 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 , 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 GU966591-GU966631.
SNP detection and experimental validation
We used SOAPsnp  to call SNPs for each variety. After setting
-u -r 0.01 -m,
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' , which can not be easily influenced by rare alleles examined . We set the parameters in the software Haploview  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.
Five silkworm mt sequences (with accession number AB070264, AY301620, NC_003395, NC_004622, and NC_012727) were firstly aligned using MUSCLE v3.7  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  by using Neighbour-Joining method and Mrbayes v3.1.2  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 . The site-wise likelihood file is derived from PhyML calculatation .
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 .
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 . 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 . We then applied Williams' correction to calculate the G statistic . Moreover, we also preformed the two likelihood ratio tests (LRTs), based on widely used branch-site models from PAML 4.2 , to detect any positive selection.
Determining mapping strategy and mapping summary
Before applying the strategy of mapping raw reads onto different references, first round mapping was initiated, which only took a single mt genome C108 as a reference to be mapped. After SNP calling, there was an obvious separation among these three groups in phylogeny analysis (data not shown). This impressed us that we should use different reference datasets based on silkworm classification as a reference for each group, so that we can get more realistic genetic variation. Before we clearly acknowledge the separation among these three populations, it is essential to perform these alignments. Here, we just take the mt sequence into consideration; however, we must introduce the nuclear DNA sequence, because there are some similar regions between the mt sequence and nuclear genome. A sequencing read would have more than one hit mapped to the whole genome (the mt genome plus the nuclear genome), while it would be with unique hits to the mt sequence when we excluded the nuclear DNA sequence. The latter, exclusive mt genome situation would interrupt meaningful SNP detection .
After aligning short reads onto the reference genome by using the SOAP alignment program , we obtained a 2319.56 fold effective depth for all the 41 varieties and each ranged from 8.99 to 165.42 fold. The minimal depth 8.99 fold covering 93.71% of the mt genome (see Table 1) suggests that there was sufficient content to detect the variation. Though this fluctuation in depth was apparent, we found there was no significant difference [Mann Whitney U (MWU), P = 0.36] when one takes the depth between domesticated silkworms and wild ones into comparison.
Mitochondria diversity and the higher polymorphism level in wild silkworm
Statistical summary for three different groups.
Effective Length (bp)
1.14 × 10-3
Chinese wild silkworms
6.20 × 10-3
Japanese wild silkworms
2.02 × 10-3
No recombination event in silkworm mt genome
Furthermore, there was also no correlation between genome distance and |D'|, and distribution of site pairs with |D'| < 1 were randomly distributed with physical distance (Pearson's correlation coefficient, ρ = 5.31 × 10-3; P = 0.80). Although this is the first investigation for silkworm mitochondria recombination, there seems to be only limited occurrence, making it unnecessary to consider its effect on population structure or evolution.
To evaluate the level of confidence in tree selection, we test whether alternative placements of Japanese wild silkworm among the Bombycidae could be accepted on basis of the approximately unbiased (AU) test and the weighted Shimodaira-Hasegawa (WSH) test . With our dataset, they reject all the alternative scenarios (P < 1 × 10-6). As mentioned in the Method section, the reference of the domesticated and Chinese wild group is from C108 and Ankang, respectively, but another individual from these two lines did not cluster together with its references (Figure 2A and Figure 2B). This result suggests the intraspecific mtDNA sequence diversity would be higher than the interspecific diversity.
Stable population size in domesticated and Chinese wild population
Positive selection detection
Number of replacement and synonymous substitutions within and between B. mori and B. mandarina
Fixed between species
Polymorphic within species
Gadj* = 9.75
P = 1.8 × 10-3
We also looked for positive selection using branch-site likelihood ratio tests for this gene, showing significant difference (P = 0.014) between the foreground and background. These two kinds of independent tests both show statistical power in this scenario, so we can believe the cytb gene is important to the domesticated silkworm, which will be elucidated in the discussion.
The silkworm is the best-characterized model for biochemical, genetic, and genomic studies of the order Lepidoptera, and for insect domestication. Here we report a full-scale genome-wide mitochondrial map, which includes all the major silkworm classes in Bombycidae, and unveils - at a single base-pair resolution - the distinction between the two wild silkworm subgroups. This map will provide a valuable resource for further study on silkworm mitochondria evolution, and will accelerate the research of functional identification among silkworm genes.
We estimated θπ values for SNPs in mitochondria and found they were 1.14 × 10-3 and 6.20 × 10-3 in B. mori and B. mandarina, respectively, and we also know that θπ values in the nuclear genome were 1.36 × 10-2 and 1.53 × 10-2 in B. mori and B. mandarina, respectively . The relatively larger polymorphism level increase in nuclear genome is mainly due to the background transposable element content  and rapid decay of LD . The pattern in humans, whose mutation rate is 2.8 × 10-3  and 0.8 × 10-3  in our mt and nuclear genomes, respectively, provides an interesting contrast to that of the silkworm. But compared to Drosophila simulans, we found a similar phenomenon of a nuclear genome mutation rate (1.8 × 10-2)  being higher than the mitochondrial rate (1.1 × 10-3) . This indicated that the insect nuclear genome could tolerate more alteration, and the differential of effective population size between nuclear genome and mt genome would be more distinct in the insect. Although the diversity of mitochondria was smaller than that of the nuclear, the difference between the domesticated silkworms and the wild silkworms was apparent, which suggested they had exclusive genetic background, and local inbreeding had produced more effect on the domesticated population.
In recent years, scientists have only used limited samples or partial genome datasets to infer the phylogenic relationship amongst silkworm clades. The main features of these analyses have been on the basis of incomplete genetic information. Although it is possible for them to reach the conclusion that the true ancestor of the domesticated silkworm is the Chinese wild group, representatives of several basal lineages remain missing. Here, taking advantage of 46 complete mt genomes, we can confidently verify the consensus phylogeny of Arunkumar et al and Pan et al that B. mori was domesticated from the Chinese B. mandarina. Both the AU test and WSH test also provided further support on this argument. We also endeavored to study the major genetic characteristics left by silkworm domestication process, such as the different level of polymorphisms within each group, and the selection force shaping on the mt genome. In the mitochondrion of eukaryotes and in aerobic prokaryotes, the cytb protein is a component of respiratory chain complex III, and is related to electron transfer. Mutation in the site 318 of human cytb gene causes severe deficiency in the respiratory chain enzyme involving in patient exercise intolerance . In silkworms, previous studies also showed that nine midgut-enriched genes related to energy metabolism under selection play a critical role in food digestion and nutrient absorption . The identification of cytb gene indicates that mt genes probably contribute equally to energy metabolism process, and is important to the domesticated silkworm. This is the first time we sort out a gene under positive selection in silkworm mtDNA, and it may provide a new insight into the silkworm mtDNA evolution process.
The constant population size inference for B. mori and B. mandarina in China indicated that the Chinese wild silkworm led diverse and distinct lifestyles in the wild mulberry, and the domesticated silkworm retained and conserved their abundant interspecific polymorphism. These results would pave the way for breeding strategies later.
B. mori is not only well adapted to human handling, but is wholly dependent on humans for survival, in addition to being well differentiated trait-wise from its wild cousin. This makes silkworm domestication a distinctive event in agricultural history, deserving the same attention as domestication of livestock and crop plants; especially as it took place in a different geographical region (Asia vs. the Fertile Crescent) , and in a distinctly different culture from the earliest known, and better studied domestication events.
Based on whole genome comparative analysis, with 41 silkworm mitochondrion and 5 other available mt genome sequences, we identified 347 SNPs in the silkworm mt genome. Using this data, we infer that the silkworm mt genome did not undergo recombination, and find it very strongly implied that the Chinese B. mandarina is the most recent ancestor of domesticated B. mori, in confirmation of the conclusions of Arunkumar et al and Pan et al. It is found that, after silkworm domestication and group divergence, the effective population sizes in both groups remained constant, and that identifiable genes underwent powerful selection pressures in B. mori.
These comprehensive overviews will provide new insight into the evolution of the silkworm. It is believed that the silkworm domestication event is of comparative significance equal to the concurrent domestication events of crop plants and animal livestock, during the same time period of human history.
We thank Dr. Kazuei Mita at National Institute of Agrobiological Sciences in Japan for his kindly providing the Japanese wild silkworm genomic DNA. We thank two anonymous referees providing valuable comments. This work is supported by grants from the National Basic Research Program of China (No. 2005CB121000), the Program for Changjiang Scholars and Innovative Research Team in University (IRT0750), and the Shenzhen Bureau of Science Technology & Information, China (JC200903190764A).
- Avise JC, Arnold J, Ball RM, Bermingham E, Lamb T, Neigel JE, Reeb CA, Saunders NC: Intraspecific Phylogeography - the Mitochondrial-DNA Bridge between Population-Genetics and Systematics. Annu Rev Ecol Syst. 1987, 18: 489-522.View ArticleGoogle Scholar
- Xiang ZH, Huang JT, Xia JG, Lu C: Biology of sericulture. 2005, Beijing: China Forestry Publishing HouseGoogle Scholar
- Goldsmith MR, Shimada T, Abe H: The genetics and genomics of the silkworm, Bombyx mori. Annu Rev Entomol. 2005, 50: 71-100. 10.1146/annurev.ento.50.071803.130456.View ArticlePubMedGoogle Scholar
- Tomita M, Munetsuna H, Sato T, Adachi T, Hino R, Hayashi M, Shimizu K, Nakamura N, Tamura T, Yoshizato K: Transgenic silkworms produce recombinant human type III procollagen in cocoons. Nat Biotechnol. 2003, 21 (1): 52-56. 10.1038/nbt771.View ArticlePubMedGoogle Scholar
- Arunkumar KP, Metta M, Nagaraju J: Molecular phylogeny of silkmoths reveals the origin of domesticated silkmoth, Bombyx mori from Chinese Bombyx mandarina and paternal inheritance of Antheraea proylei mitochondrial DNA. Mol Phylogenet Evol. 2006, 40 (2): 419-427. 10.1016/j.ympev.2006.02.023.View ArticlePubMedGoogle Scholar
- Pan M, Yu Q, Xia Y, Dai F, Liu Y, Lu C, Zhang Z, Xiang Z: Characterization of mitochondrial genome of Chinese wild mulberry silkworm, Bomyx mandarina (Lepidoptera: Bombycidae). Sci China C Life Sci. 2008, 51 (8): 693-701. 10.1007/s11427-008-0097-6.View ArticlePubMedGoogle Scholar
- Xia Q, Guo Y, Zhang Z, Li D, Xuan Z, Li Z, Dai F, Li Y, Cheng D, Li R, et al: Complete resequencing of 40 genomes reveals domestication events and genes in silkworm (Bombyx). Science. 2009, 326 (5951): 433-436. 10.1126/science.1176620.PubMed CentralView ArticlePubMedGoogle Scholar
- Ballard JW, Kreitman M: Unraveling selection in the mitochondrial genome of Drosophila. Genetics. 1994, 138 (3): 757-772.PubMed CentralPubMedGoogle Scholar
- Ballard JW, Kreitman M: Is mitochondrial DNA a strictly neutral marker?. Trends Ecol Evol. 1995, 10: 485-488. 10.1016/S0169-5347(00)89195-8.View ArticlePubMedGoogle Scholar
- Rand DM: The units of selection on mitochondrial DNA. Annu Rev Ecol Syst. 2001, 32: 415-448. 10.1146/annurev.ecolsys.32.081501.114109.View ArticleGoogle Scholar
- Wang J, Xia Q, He X, Dai M, Ruan J, Chen J, Yu G, Yuan H, Hu Y, Li R, et al: SilkDB: a knowledgebase for silkworm biology and genomics. Nucleic Acids Res. 2005, D399-402. 33 DatabaseGoogle Scholar
- Yukuhiro K, Sezutsu H, Itoh M, Shimizu K, Banno Y: Significant levels of sequence divergence and gene rearrangements have occurred between the mitochondrial genomes of the wild mulberry silkmoth, Bombyx mandarina, and its close relative, the domesticated silkmoth, Bombyx mori. Mol Biol Evol. 2002, 19 (8): 1385-1389.View ArticlePubMedGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24 (5): 713-714. 10.1093/bioinformatics/btn025.View ArticlePubMedGoogle Scholar
- Li R, Li Y, Fang X, Yang H, Wang J, Kristiansen K: SNP detection for massively parallel whole-genome resequencing. Genome Res. 2009, 19 (6): 1124-1132. 10.1101/gr.088013.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Lewontin RC: The Interaction of Selection and Linkage. I. General Considerations; Heterotic Models. Genetics. 1964, 49 (1): 49-67.PubMed CentralPubMedGoogle Scholar
- Wang WY, Barratt BJ, Clayton DG, Todd JA: Genome-wide association studies: theoretical and practical concerns. Nat Rev Genet. 2005, 6 (2): 109-118. 10.1038/nrg1522.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17 (12): 1246-1247. 10.1093/bioinformatics/17.12.1246.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9 (3): 552-569.PubMedGoogle Scholar
- Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129 (2): 555-562.PubMed CentralPubMedGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.View ArticlePubMedGoogle Scholar
- McDonald JH, Kreitman M: Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991, 351 (6328): 652-654. 10.1038/351652a0.View ArticlePubMedGoogle Scholar
- Sokal RR, Rohlf FJ: Biometry: the principles and practice of statistics in biological research. 1994, New York: Freeman, 3Google Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Awadalla P, Eyre-Walker A, Smith JM: Questioning Evidence for Recombination in Human Mitochondrial DNA. Science. 2000, 288 (5473): 1931a-10.1126/science.288.5473.1931a.View ArticlePubMedGoogle Scholar
- Guo SW: Linkage disequilibrium measures for fine-scale mapping: a comparison. Hum Hered. 1997, 47 (6): 301-314. 10.1159/000154430.View ArticlePubMedGoogle Scholar
- Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Syst Biol. 2002, 51 (3): 492-508. 10.1080/10635150290069913.View ArticlePubMedGoogle Scholar
- Watson E, Bauer K, Aman R, Weiss G, von Haeseler A, Paabo S: mtDNA sequence diversity in Africa. Am J Hum Genet. 1996, 59 (2): 437-444.PubMed CentralPubMedGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMed CentralPubMedGoogle Scholar
- Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 133 (3): 693-709.PubMed CentralPubMedGoogle Scholar
- Xia Q, Zhou Z, Lu C, Cheng D, Dai F, Li B, Zhao P, Zha X, Cheng T, Chai C, et al: A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004, 306 (5703): 1937-1940. 10.1126/science.1102210.View ArticlePubMedGoogle Scholar
- Ingman M, Kaessmann H, Paabo S, Gyllensten U: Mitochondrial genome variation and the origin of modern humans. Nature. 2000, 408 (6813): 708-713. 10.1038/35047064.View ArticlePubMedGoogle Scholar
- Wang J, Wang W, Li R, Li Y, Tian G, Goodman L, Fan W, Zhang J, Li J, Guo Y, et al: The diploid genome sequence of an Asian individual. Nature. 2008, 456 (7218): 60-65. 10.1038/nature07484.PubMed CentralView ArticlePubMedGoogle Scholar
- Begun DJ, Holloway AK, Stevens K, Hillier LW, Poh YP, Hahn MW, Nista PM, Jones CD, Kern AD, Dewey CN, et al: Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans. PLoS Biol. 2007, 5 (11): e310-10.1371/journal.pbio.0050310.PubMed CentralView ArticlePubMedGoogle Scholar
- Andreu AL, Hanna MG, Reichmann H, Bruno C, Penn AS, Tanji K, Pallotti F, Iwata S, Bonilla E, Lach B, et al: Exercise intolerance due to mutations in the cytochrome b gene of mitochondrial DNA. N Engl J Med. 1999, 341 (14): 1037-1044. 10.1056/NEJM199909303411404.View ArticlePubMedGoogle Scholar
- Driscoll CA, Macdonald DW, O'Brien SJ: From wild animals to domestic pets, an evolutionary view of domestication. Proc Natl Acad Sci USA. 2009, 106 (Suppl 1): 9971-9978. 10.1073/pnas.0901586106.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.