The Complete Mitogenome of Bombyx Mori Strain Dazao (Lepidoptera: Bombycidae) and Comparison with Other Lepidopteran Insects
Total Page:16
File Type:pdf, Size:1020Kb
View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Genomics 101 (2013) 64–73 Contents lists available at SciVerse ScienceDirect Genomics journal homepage: www.elsevier.com/locate/ygeno The complete mitogenome of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) and comparison with other lepidopteran insects Qiu-Ning Liu 1, Bao-Jian Zhu 1, Li-Shang Dai, Chao-Liang Liu ⁎ College of Life Science, Anhui Agricultural University, 130 Changjiang West Road 230036, PR China article info abstract Article history: The complete mitochondrial genome (mitogenome) of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) Received 31 August 2012 was determined to be 15,653 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes Accepted 6 October 2012 and a A+T-rich region. It has the typical gene organization and order of mitogenomes from lepidopteran in- Available online 13 October 2012 sects. The AT skew of this mitogenome was slightly positive and the nucleotide composition was also biased toward A+T nucleotides (81.31%). All PCGs were initiated by ATN codons, except for cytochrome c oxidase Keywords: subunit 1 (cox1) gene which was initiated by CGA. The cox1 and cox2 genes had incomplete stop codons Mitochondrial genome Lepidoptera consisting of just a T. All the tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. Bombyx mori The A+T-rich region of the mitogenome was 495 bp in length and consisted of several features common Phylogeny to the lepidopteras. Phylogenetic analysis showed that the B. mori Dazao was close to Bombycidae. © 2012 Elsevier Inc. All rights reserved. 1. Introduction phylogenetic analyses to the selected species from Lepidoptera, Orthoptera and Diptera based on the mitogenome sequences were Insect mitogenome DNA (mtDNA) is a circular molecule with performed using the neighbor-joining (NJ) method. 14–19 kb in length and contains 13 protein coding genes (PCGs), sub- units 6 and 8 of the ATPase (atp6 and atp8), cytochrome c oxidase subunits 1–3(cox1-cox3), cytochrome B (cob), NADH dehydrogenase Table 1 The complete mitochondrial genome of Lepidoptera. subunits 1–6 and 4 L (nad1-6 and nad4L), two ribosomal RNA genes encoding the small and large subunit rRNAs (rrnL and rrnS), 22 trans- Species Length (bp) Accession number References fer RNA (tRNA) genes and a A+T-rich region which contains initia- Bombyx mori Dazao 15, 653 This study tion sites for transcription and replication [1–3]. Mitochondrial DNA Bombyx mori Xiafang 15, 664 AY048187 [4] has been widely used as an informative molecular marker for diverse Bombyx mori C108 15, 656 AB070264 [5] Japanese Bombyx mandarina 15, 928 NC_003395 [5] evolutionary studies among species [3]. Chinese Bombyx mandarina 15, 682 AY301620 [6] Although Lepidoptera is the 2nd most numerous order of insects, only Antheraea pernyi 15, 575 AY242996 [7] a few complete or near-complete mitogenomes are currently available in Antheraea yamamai 15, 338 EU726630 [8] GenBank (Table 1). The silk-producing insects with economic value in Eriogyna pyretorum 15, 327 FJ685653 [9] Lepidoptera belong to two families of moth, Bombycidae and Saturniidae. Samia cynthia ricini 15, 366 JN215366 [10] Actias selene 15, 236 JX186589 [11] The complete mitogenomes of three species of Bombycidae and six spe- Caligula boisduvalii 15, 360 NC_010613 [12] cies of Saturniidae were sequenced: Bombyx mori Xiafang [4], Bombyx Manduca sexta 15, 516 EU286785 [13] mori C108 [5], Japanese B. mandarina [5] and Chinese B. mandarina [6] be- Hyphantria cunea 15, 481 GU592049 [14] longing to family Bombycidae; Antheraea pernyi [7], Antheraea yamamai Helicoverpa armigera 15, 347 GU188273 [15] Ochrogaster lunifer 15, 593 AM946601 [16] [8], Eriogyna pyretorum [9], Samia cynthia ricini [10], Actias selene [11] Ostrinia nubilalis 14, 535 NC_003367 [17] and Caligula boisduvalii [12] belonging to the family Saturniidae. Chilo suppressalis 15, 456 HQ860290 [18] Bombyx mori strain Dazao, a member of the Bombycidae family, Diatraea saccharalis 15,490 FJ240227 [19] is an important model organism and used as a bioreactor for the Phthonandria atrilineata 15, 499 EU569764 [20] production of recombinant proteins [28]. In the present study, the Adoxophyes honmai 15, 680 DQ073916 [21] Grapholita molesta 15, 776 HQ116416 [22] complete mitogenome of B. mori strain Dazao was sequenced. The Spilonota lechriaspis 15, 368 HM204705 [23] Artogeia melete 15, 140 EU597124 [24] Coreana raphaelis 15, 314 NC_007976 [25] ⁎ Corresponding author. Fax: +86 551 5786201. Acraea issoria 15, 245 GQ376195 [26] E-mail address: [email protected] (C.-L. Liu). Fabriciana nerippe 15, 140 JF504707 [27] 1 These authors contributed equally to this work. 0888-7543/$ – see front matter © 2012 Elsevier Inc. All rights reserved. http://dx.doi.org/10.1016/j.ygeno.2012.10.002 Q.-N. Liu et al. / Genomics 101 (2013) 64–73 65 Table 2 2.2. PCR amplification and sequencing Primers used for amplification of the mitogenome of B. mori Dazao. Primer pair Primer sequence (5′→3′) For amplification of the entire mitogenome of B. mori strain Dazao, nine primer sets were designed according to the conserved nucleo- F1 GCTTTTGGGCTCATACCTCA R1 GATGAAATACCTGCAAGATGAAG tide sequences of known mitochondrial sequences from Lepidoptera F2 TGGAGCAGGAACAGGATGAAC insects [4–12] and synthesized in Sunbiotech Co., Ltd. (Beijing, R2 GAGACCADTACTTGCTTTCAG China) (Table 2). The fragments were amplified using Aidlab Taq Kit F3 ATTTGTGGAGCTAATCATAG (Aidlab Co., Beijing, China) according to the manufacturer's instruc- R3 GGTCAGGGACTATAATCTAC F4 TCGACCTGGAACTTTAGC tions. The PCR was performed under the following conditions: R4 GCAGCTATAGCCGCTCCTACT 3 min at 94 °C, followed by 35 cycles of 30 s at 94 °C, 1–3 min at F5 TAAAGCAGAAACAGGAGTAG 55–60 °C, 10 min at 72 °C. The PCR products were separated by aga- R5 ATTGCGATATTATTTCTTTTG rose gel electrophoresis (1% w/v) and purified using a DNA gel extrac- F6 ACATTCTTAGGTGGATTA tion kit (Aidlab Co., Beijing, China). The purified PCR products were R6 GTTAAAGTGGCATTATCT F7 GGAGCTTCTACATGAGCTTTTGG ligated into the T-vector (TaKaRa Co., Dalian, China) and sequenced R7 GTTTGCGACCTCGATGTTG at least three times at Sunbiotech Co., Ltd., (Beijing, China). F8 GGTCCCTTACGAATTTGAATATATCCT R8 AAACTAGGATTAGATACCCTATTAT 2.3. Sequence assembly and gene annotation F9 CTCTACTTTGTTACGACTTATT R9 TCTAGGCCAATTCAACAACC Sequence annotation was performed using the blast tools in NCBI web site (http://blast.ncbi.nlm.nih.gov/Blast)andDNAStarpackage 2. Materials and methods (DNAStar Inc. Madison, USA). Identification of tRNA genes was verified using the tRNAscan-SE program. The potential stem–loop secondary 2.1. DNA extraction structures within these tRNA gene sequences were calculated using the tRNAscan-SE Search Server (http://lowelab.ucsc.edu/tRNAscan-SE/) The silkworm strain Dazao was maintained in our laboratory. All [29]. The secondary structures of tRNA genes that could not be predicted larvae were fed with fresh mulberry leaves under normal conditions. using the tRNAscan-SE were analyzed by comparison with the Total DNA was isolated from the third day of fifth instar larvae using nucleotide sequences of other insect tRNA sequences [4–12].The the Aidlab Genomic DNA Extraction Kit (Aidlab Co., Beijing, China) PCGs were identified by sequence similarity with B. mori C108 [5]. according to the manufacturer's instructions. Extracted DNA was The nucleotide sequences of PCGs were translated with the inverte- used for PCR amplification of the complete mitogenome. brate mitogenome genetic code. Alignments of PCGs from various Fig. 1. Map of the mitogenome of B. mori strain Dazao. tRNA genes are labeled according to the IUPAC-IUB. Single letter amino acids above the bar indicate coding sequence on major strand whereas below the bar indicates coding on minor strand. Anti-clockwise PCGs or rRNA genes are located on L strand and others are located on H strand. 66 Q.-N. Liu et al. / Genomics 101 (2013) 64–73 lepidopteran mitogenomes were performed using Clustal X [30]. A+T-rich region (Fig. 1). The order of genes and the orientation of Composition skewness was calculated according to the formulas: the mitogenome of B. mori Dazao are identical to the sequenced lep- AT skew=[A− T]/[A+T]; GC skew=[G− C]/[G+ C] [31].Tandem idopteran mitogenomes. The placement of the trnM tRNA gene in repeats in the A+ T-rich region were predicted using the Tandem the B. mori Dazao mitogenome differs from the ancestral gene order Repeats Finder program (http://tandem.bu.edu/trf/trf.html) [32]. observed in insects. The ancestral order is: A+T-rich region, trnI, trnQ, trnM, nad2 [39]. The placement of trnM may therefore represent 2.4. Phylogenetic analysis a molecular feature exclusive to lepidopteran mtDNAs [3,37]. The nu- cleotide composition of the mitogenome of B. mori Dazao is: A=6732 To reconstruct the phylogenetic relationship among lepidopteran (43.01%), T=5995 (38.30%), G=1154 (7.37%) and C=1772 insects, the complete mitogenomes of 26 lepidopteran species were (11.32%). As observed in other lepidopterans, the nucleotide compo- obtained from the GenBank database (Table 1). These mitogenomes sition of the B. mori Dazao mitogenome is biased toward A+T were divided into 11 lepidopteran superfamilies within the lepidopter- (81.31%) (Table 4). an suborder. The mitogenomes of Locusta migratoria (NC_001712), Drosophila yakuba (NC_001322) and Anopheles gambiae (NC_002084) were used as outgroups [33–35]. The amino acid sequences of each of 3.2. Overlapping and intergenic spacer regions the 13 mitochondrial PCGs were aligned using default settings and concatenated. The concatenated set of amino acid sequences from the The B. mori Dazao mitogenome has a total of 31 bp overlap be- 13 PCGs was used in phylogenetic analysis, which was performed using tween genes in eight locations. Two longest eight nucleotides over- neighbor-joining (NJ) method with the MEGA version 5.0 program [36]. laps located between trnW and trnC as well as trnF and nad5 were found (Table 3), however, the eight nucleotides overlap between 3.