Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Yu-Ping Li^I; Wu Song^I; Sheng-Lin Shi^I; Yan-Qun Liu^{I, II}; Min-Hui Pan^II; Fang-Yin Dai^II; Cheng Lu^II; Zhong-Huai Xiang^II

^ICollege of Bioscience and Biotechnology, Shenyang Agricultural University, Liaoning, Shenyang, China

^IIThe Key Sericultural Laboratory of Agricultural Ministry, Southwest University, Chongqing, China

^{Send correspondence to}

ABSTRACT

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.

Key words: nucleotide substitution pattern, mitochondrial genome, Bombyx mori, Chinese Bombyx mandarina, Japanese Bombyx mandarina.

The mitochondrial (mt) DNA of insects is a self-replicating, circular DNA molecule about 14-20 kb in size, encoding a conserved set of 37 genes (13 protein genes, 22 tRNA genes, and two rRNA genes) and an A+T-rich segment known as control region (Wolstenholme, 1992; Boore, 1999). Several comparative mt genomics studies for investigation of the nucleotide variation and evolutionary patterns have been carried out, including Drosophila melanogaster subgroup members (Ballard, 2000), Ostrinia nubilalis and O. furnicalis (Coates et al., 2005), as well as Bombyx mori and its close relative Japanese B. mandarina (Yukuhiro et al., 2002). Generally, the A+T-rich region shows a higher level of sequence variability than other regions of the genome (Fauron and Wolstenholme, 1980; Shao et al., 2001)

The mulberry silkworm Bombyx mori is the only truly domesticated insect. It is thought to be domesticated from the wild mulberry silkworm, B. mandarina, about 5,000-10,000 years ago (Goldsmith et al., 2005). B. mori and B. mandarina are morphologically and physiologically similar. However, two types of wild mulberry silkworms with uniform morphology and different chromosome numbers per haploid genome are currently found in mulberry fields. The Japanese B. mandarina (JaBm) living in Japan and Korea has 27 chromosomes in its haploid genome, whereas the Chinese B. mandarina (ChBm), distributed over China, has 28 chromosomes per haploid genome, like its domesticated counterpart B. mori that also carries 28 chromosomes (Astaurov et al., 1959). JaBm was thought to be a probable wild ancestor of B. mori, and two sets of biased patterns in the mt genome nucleotide substitution between the two have been confirmed, i.e., a difference in preference of transitional changes between major and minor DNA strands and a variation in the degree of nucleotide substitution by genomic regions (Yukuhiro et al., 2002).

Phylogenetic analysis (Arunkumar et al., 2006; Pan et al., 2008) and historical records (Goldsmith et al., 2005) make it clear that the domesticated silkworm was directly domesticated from the Chinese wild silkworm rather than from the Japanese wild silkworm. We determined the complete mt genome of ChBm (GenBank accession number: AY301620) that contains a typical gene complement, order, and arrangement identical to that of B. mori and JaBm (Pan et al., 2008). In the present study, we examined the nucleotide variation pattern of the complete mt genome between B. mori and ChBm. The purpose of this study was to determine whether the nucleotide variation pattern between B. mori and its wild ancestor, ChBm, is consistent with that between B. mori and its close relative, JaBm.

The mt genome sequences of B. mori C108 (AB070264), Backokjam (NC_002355), Aojuku (AB083339), Xiafang (AY048187), and JaBm (AB070263) were downloaded from GenBank, aligned by Clustal X (Thompson et al., 1997), and the alignment was analyzed by MEGA version 4.1 (Tamura et al., 2007), including the variable nucleotide sites and the substitution sites.

The four mt genome sequences of B. mori C108 (15,656 bp), Backokjam (15,643 bp), Aojuku (15,635 bp), and Xiafang (15,664 bp) were aligned to produce a 15,685 nt consensus alignment, of which only 89 (0.57%) nucleotide sites are variable, including insertions-deletions (indels) (data not shown). Within the four B. mori strains, only 33 substitutions were detected, including 24 transition (ts) and 9 transversion (tv) mutations. The variable nucleotide sites ranged from 63 to 74 nt when each two of the four B. mori mt genome sequences were compared, showing no significant differences in variable nucleotide rate (χ² = 1.25, d.f. = 5, p > 0.05).

However, the mt genome sequences of ChBm (15,682 bp) and JaBm (15,928 bp) were aligned to produce a 15,970 nt consensus alignment, of which 639 (4.00%) nucleotide sites are variable including indels, and 325 substitutions were detected. Compared to these data observed within the four B. mori strains, the sequence divergence within B. mandarina is significantly greater (4.00% vs 0.57%; χ² = 415.26, d.f. = 1, p < 0.001). Moreover, we identified 484 variable nucleotide sites including indels (3.08% sequence divergence) from a 15,722 nt consensus alignment between the mt genomes of B. mori C108 and ChBm, whereas 855 variable nucleotide sites including indels (5.35% sequence divergence) were identified from a 15,970 nt consensus alignment between B. mori C108 and JaBm. Compared to the sequence divergence observed between B. mori C108 and JaBm, that between B. mori and ChBm was significantly smaller (3.08% vs 5.35%; χ² = 101.36, d.f. = 1, p < 0.001). These results led us to further compare the nucleotide substitution pattern between B. mori C108 and ChBm with that between B. mori C108 and JaBm.

A total of 381 substitutions were detected between the mt genomes of B. mori C108 and ChBm: 230 transition and 151 transversion mutations (ts:tv = κ ≈ 230/151 = 1.52). This value is significantly deviatory from neutral expectation (1:2, χ² = 125.30, d.f. = 1, p < 0.001), indicating that evolutionary pressures are acting upon the two mt genomes. However, between B. mori C108 and JaBm, 514 substitutions were detected: 414 transition and 100 transversion mutations (κ ≈ 4.14), also significantly deviated from neutral expectation (1:2, χ² = 517.47, d.f. = 1, p < 0.001). Excess transition mutation was also reported in Drosophila (κ = 761/180 ≈ 4.23) and Ostrinia (κ = 134/48 ≈ 2.88) mt genomes, and attributed to non-neutral evolutionary forces or population effects (Ballard, 2000; Coates et al., 2005). Compared to the data observed between B. mori and JaBm, the transition mutation ratio between B. mori and ChBm was significantly smaller (1.52 vs 4.14) (χ² = 44.14, d.f. = 1, p < 0.001), suggesting that their nucleotide substitution pattern is different.

The nucleotide differences between B. mori and ChBm mt genes except for tRNA genes are shown in Table 1. The numbers and rates of nucleotide difference between the mt genes of B. mori and ChBm were found to be decreased compared to those seen between B. mori and JaBm, except for the srRNA gene and A+T-rich region, in which the numbers and rates of nucleotide difference were increased. The full srRNA sequences of B. mori C108 and ChBm were aligned to produce a 785 nt consensus alignment, of which 12 substitutions were detected, showing a higher nucleotide sequence divergence than the lrRNA gene. Comparing B. mori C108 and ChBm, the NADH dehydrogenase subunit encoding nad6 gene was the most conserved among genes or regions, followed by the ATP synthase F0 subunit encoding atp8 gene. The cytochrome b (cob) gene was the most variable, followed by the cytochrome c oxidase subunit encoding cox3 gene. The nucleotide sequence divergence in the A+T-rich region was moderate among these genes or regions. However, comparing the mt genes of B. mori and JaBm, the srRNA gene was the most conserved among genes or regions, and nucleotide sequence divergences in the lrRNA gene and the A+T-rich region were also low, compared to those of protein-coding genes.

In this study, we found no significant differences in transition rates for the protein-coding genes encoded by the major strand (22 A-G and 101 C-T transitions in 6,909 bp) between the mt genomes of B. mori C108 and ChBm, compared to those encoded by the minor strand (48 A-G and 13 C-T transitions in 4296 bp) (χ² = 2.13, d.f. = 1, p > 0.05). However, on the major strand there were significantly more changes for T-C transitions in protein-coding genes than observed on the minor strand (χ² = 35.35, d.f. = 1, p < 0.001), while on the minor strand there were significantly more changes for A-G transitions than observed on the major strand (χ² = 27.23, d.f. = 1, p < 0.001). These data are consistent with those previously reported (Ballard 2000; Yukuhiro et al., 2002).

The rates of nucleotide substitutions between the mt genomes of B. mori C108 and JaBm vary by genomic region: the five genes surrounding the A+T-rich region and two rRNA genes (Class 1: nad2, cox1, cox2, atp6, and nad1) are more conserved than the remaining genes (Class 2: cox3, nad4, nad5, and cob), whereas no significant difference in substitution rates within Classes 1 and Classes 2 has been detected (Yukuhiro et al., 2002). In this study, we found that, between the mt genomes of B. mori and ChBm, there is a significant difference in the nucleotide substitution rate for the 13 protein-coding genes (χ² = 39.68, d.f. = 12, p < 0.005), which is consistent with that found between B. mori and JaBm. Further analysis showed that between the mt genomes of B. mori and ChBm the average substitution rate of Class 1 (0.0195) was significantly smaller than that of Class 2 (0.0308) (χ² = 12.70, d.f. = 1, p < 0.005); moreover, a significant difference in the substitution rate within Classes 2 was detected (χ² = 15.30, d.f. = 3, p < 0.005), but there was no significant difference in substitution rates within Classes 1 (χ² = 0.77, d.f. = 4, p > 0.9).

The A+T-rich region is the only major non-coding region in the mt genomes of animals. Based on the assumption that the A+T-rich region is the fastest evolving region of the mtDNA (Zhang and Hewitt, 1997; Giuffra et al., 2000), it has become popular as a molecular marker for population genetic and phylogeographic studies of animals, including insects. However, according to Zhang and Hewitt (1997), in terms of nucleotide substitution, the A+T-rich region might not be the fastest evolving region of the mtDNA of insects. Between the mt genes of B. mori and JaBm, the nucleotide sequence divergence in the A+T-rich region is low, higher only than that of the srRNA gene, whereas between B. mori and ChBm, the A+T-rich region presents a moderate sequence divergence, lower only than those of the three protein-coding genes cob, cox3, and nad4L. These observations of nucleotide sequence divergence in the A+T-rich region showed that this is not the fastest evolving gene or region in Bombyx mt genomes, providing direct evidence for the above mentioned assumption of Zhang and Hewitt (1997).

The primary finding of this study was that the nucleotide substitution pattern between the mt genomes of B. mori and ChBm is different from that between B. mori and JaBm. Although the preference of transition by DNA strands between B. mori and ChBm was consistent with that between B. mori and JaBm, the regional variation in nucleotide substitution rate showed a different feature. Comparative mt genomics revealed a lower level of transition mutation (ts:tv = 1.52) between B. mori and ChBm, but a significantly higher level of transition mutation (ts:tv = 4.14) between B. mori and JaBm. It has been suggested that the A+T richness in the mt genome would cause an apparently lower transition bias ratio in closely related species (Tamura, 1992; Yu et al., 1999; Liao and Lu, 2000; Arunkumar et al., 2006). However, the present investigation of the A+T content of the mt genome from the three samples showed that they were almost identical (from 81.36% in B. mori to 81.68% in ChBm), indicating that the heavy bias towards A+T is not suitable to explain this case. Therefore, the distinct nucleotide substitution pattern suggests that the mt genomes of ChBm and JaBm are the result of different evolutionary forces. The mechanism responsible for that difference remains unclear. A possible explanation is that the environmental selection effect could lead to different mutational biases. In line with previous reports (Messier and Stewart, 1997; Wise et al., 1998; Creevey and McInerney, 2002; Jansa et al., 2003), this study presents evidence for mtDNA selection.

Acknowledgments

This work was supported in part by grants from the National Natural Science Foundation of China (No. 30800803), the Scientific Research Project for Commonwealth Industry of the Agricultural Ministry (No. nyhyzx07-020-17), the Scientific Research Project for High School of the Department of Education of the Liaoning Province (No. 2008643), and the Young Scholar Foundation of the Shenyang Agricultural University (No. 20070112), respectively.

Received: February 16, 2009; Accepted: July 1, 2009.

Associate Editor: Louis Bernard Klaczko

• Arunkumar KP, Metta M and Nagaraju J (2006) 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 40:419-427.
• Astaurov BL, Golisheva MD and Roginskaya IS (1959) Chromosome complex of Ussuri geographical race of Bombyx mandarina M. with special reference to the problem of the origin of the domesticated silkworm, Bombyx mori L. Cytology 1:327-332.
• Ballard JW (2000) Comparative genomics of mitochondrial DNA in members of the Drosophila melanogaster subgroup. J Mol Evol 51:48-63.
• Boore JL (1999) Animal mitochondrial genomes. Nucleic Acids Res 27:1767-1780.
• Coates BS, Sumerford DV, Hellmich RL and Lewis LC (2005) Partial mitochondrial genome sequences of Ostrinia nubilalis and Ostrinia furnicalis Int J Biol Sci 1:13-18.
• Creevey C and McInerney JO (2002) An algorithm for detecting directional and non-directional positive selection, neutrality and negative selection in protein coding DNA sequences. Gene 300:43-51.
• Fauron CMR and Wolstenholme DR (1980) Extensive diversity among Drosophila species with respect to nucleotide sequences within the adenine + thymine rich region of mitochondrial DNA molecules. Nucleic Acids Res 11:2439-2453.
• Giuffra E, Kijas JMH, Amarger V, Carlborg O, Jeon JT and Andersson L (2000) The origin of the domestic pig: Independent domestication and subsequent introgression. Genetics 154:1785-179l.
• Goldsmith MR, Shimada T and Abe H (2005) The genetics and genomics of the silkworm, Bombyx mori Annu Rev Entomol 50:71-100.
• Jansa SA, Lundrigan BL and Tucker PK (2003) Tests for positive selection on immune and reproductive genes in closely related species of the murine genus Mus J Biol Evol 39:123-128.
• Liao SY and Lu C (2000) Progress on animal mitochondrial genome. Progr Biochem Biophys 27:508-512 (in Chinese).
• Messier W and Stewart CB (1997) Episodic adaptive evolution of primate lysozymes. Nature 385:151-154.
• Pan MH, Yu QY, Xia YL, Dai FY, Liu YQ, Lu C, Zhang Z and Xiang ZH (2008) Characterization of the mitochondrial genome of the Chinese wild mulberry silkworm, Bomyx mandarina (Lepidoptera, Bombycidae). Sci China Ser C-Life Sci 51:693-701.
• Shao R, Campbell NJH and Barker SCB (2001) Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera). Mol Biol Evol 18:858-865.
• Tamura K (1992) The rate and pattern of nucleotide substitution in Drosophila mitochondrial DNA. Mol Biol Evol 9:814-825.
• Tamura K, Dudley J, Nei M and Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software v. 4.0. Mol Biol Evol 24:1596-1599.
• Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F and Higgins DG (1997) The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25:4876-4882.
• Wise CA, Sraml M and Easteal S (1998) Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans, but not in chimpanzees. Genetics 148:409-421.
• Wolstenholme DR (1992) Animal mitochondrial DNA: Structure and evolution. Int Rev Cytol 141:173-216.
• Yu H, Wang W, Fang S, Zhang YP, Lin FJ and Geng ZC (1999) Phylogeny and evolution for the Drosophila nasuta subgroup based on mitochondrial ND4 and ND4L gene sequences. Mol Phylogenet Evol 13:556-565.
• Yukuhiro K, Sezutsu H, Itoh M, Shimizu K and Banno Y (2002) 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 19:1385-1389.
• Zhang DX and Hewitt GM (1997) Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol 25:99-120.

Publication Dates

• Publication in this collection
  18 Dec 2009
• Date of issue
  2010

History

• Received
  16 Feb 2009
• Accepted
  01 July 2009