Characterization of mitochondrial genome of Chinese wild mulberry silkworm, Bomyx mandarina (Lepidoptera: Bombycidae)

Mitochondrial Genome Architecture and Evolutionary Origin of the Yao Silkworm, a Living Fossil of the Domestic Silkworm Bombyx mori (Lepidoptera: Bombycidae)

The Yao silkworm is a unique silkworm resource producing yellow flat plate silk that has only been reared by the Baiku Yao ethnic group in Nandan County, Guangxi Province, China for a thousand years. Here, we report the mitochondrial genomes (mitogenomes) of five Yao silkworm strains and 10 local Guangxi strains of the domestic silkworm (Bombyx mori) L. (Lepidoptera: Bombycidae), and use the resulting mitogenomes and the available Bombyx mitogenomes to characterize their genome architecture and trace the evolutionary origin of the Yao silkworm. The five Yao silkworm mitogenomes exhibited genome architectures identical to typical set of 37 mitochondrial genes (13 protein-coding genes, 22 transfer RNAs, and two ribosomal RNAs) and a high level of genome sequence similarity with the domestic silkworm. Mitogenome-based phylogenetic reconstruction provided solid evidence that the Yao silkworm shares a common ancestor with the domestic silkworm. Sliding window analysis uncovered a distinct variation pattern in the mitogenome between the Yao silkworm and the other domestic silkworm strains. The phylogenetic analyses revealed a basal placement of the Yao silkworm among all available domestic silkworm strains, indicating that the Yao silkworm is an ancient population of the domestic silkworm. Our data indicated that the Yao silkworm (B. mori) is a lineage of the domestic silkworm, which for the first time provides insights into the origin of the Yao silkworm.

Phylogeographic Relationships among Bombyx mandarina (Lepidoptera: Bombycidae) Populations and Their Relationships to B. mori Inferred from Mitochondrial Genomes

We report 37 mitochondrial genome (mitogenome) sequences of Bombyx mori strains (Lepidoptera: Bombycidae) and four of B. mandarina individuals, each preserved and collected, respectively, in South Korea. These mitogenome sequences combined with 45 public data showed a substantial genetic reduction in B. mori strains compared to the presumed ancestor B. mandarina, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, concordant to the Qinling-Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses, as has previously been detected. However, one individual that was in close proximity to the south Qinling-Huaihe line was exceptional, belonging to the northern group. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. Although a substantial genetic divergence is present between B. mandarina in southern China and Japan, a highly supported sister relationship between the two regional populations may suggest the potential origin of Japanese B. mandarina from southern China instead of the Korean peninsula.

Characterization of the complete mitochondrial genome of Orthaga olivacea Warre (Lepidoptera Pyralidae) and comparison with other Lepidopteran insects

Orthaga olivacea Warre (Lepidoptera: Pyralidae) is an important agricultural pest of camphor trees (Cinnamomum camphora). To further supplement the known genome-level features of related species, the complete mitochondrial genome of Orthaga olivacea is amplified, sequenced, annotated, analyzed, and compared with 58 other species of Lepidopteran. The complete sequence is 15,174 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a putative control region. Base composition is biased toward adenine and thymine (79.02% A+T) and A+T skew are slightly negative. Twelve of the 13 PCGs use typical ATN start codons. The exception is cytochrome oxidase 1 (cox1) that utilizes a CGA initiation codon. Nine PCGs have standard termination codon (TAA); others have incomplete stop codons, a single T or TA nucleotide. All the tRNA genes have the typical clover-leaf secondary structure, except for trnS(AGN), in which dihydrouridine (DHU) arm fails to form a stable stem-loop structure. The A+T-rich region (293 bp) contains a typical Lepidopter motifs 'ATAGA' followed by a 17 bp poly-T stretch, and a microsatellite-like (AT)13 repeat. Codon usage analysis revealed that Asn, Ile, Leu2, Lys, Tyr and Phe were the most frequently used amino acids, while Cys was the least utilized. Phylogenetic analysis suggested that among sequenced lepidopteran mitochondrial genomes, Orthaga olivacea Warre was most closely related to Hypsopygia regina, and confirmed that Orthaga olivacea Warre belongs to the Pyralidae family.

Comparative Transcriptomics Provides Insights into Reticulate and Adaptive Evolution of a Butterfly Radiation

Butterfly eyes are complex organs that are composed of a diversity of proteins and they play a central role in visual signaling and ultimately, speciation, and adaptation. Here, we utilized the whole eye transcriptome to obtain a more holistic view of the evolution of the butterfly eye while accounting for speciation events that co-occur with ancient hybridization. We sequenced and assembled transcriptomes from adult female eyes of eight species representing all major clades of the Heliconius genus and an additional outgroup species, Dryas iulia. We identified 4,042 orthologous genes shared across all transcriptome data sets and constructed a transcriptome-wide phylogeny, which revealed topological discordance with the mitochondrial phylogenetic tree in the Heliconius pupal mating clade. We then estimated introgression among lineages using additional genome data and found evidence for ancient hybridization leading to the common ancestor of Heliconius hortense and Heliconius clysonymus. We estimated the Ka/Ks ratio for each orthologous cluster and performed further tests to demonstrate genes showing evidence of adaptive protein evolution. Furthermore, we characterized patterns of expression for a subset of these positively selected orthologs using qRT-PCR. Taken together, we identified candidate eye genes that show signatures of adaptive molecular evolution and provide evidence of their expression divergence between species, tissues, and sexes. Our results demonstrate: 1) greater evolutionary changes in younger Heliconius lineages, that is, more positively selected genes in the cydno-melpomene-hecale group as opposed to the sara-hortense-erato group, and 2) suggest an ancient hybridization leading to speciation among Heliconius pupal-mating species.

The Complete Mitochondrial Genome and Expression Profile of Mitochondrial Protein-Coding Genes in the Bisexual and Parthenogenetic Haemaphysalis longicornis

The tick Haemaphysalis longicornis is widely distributed in eastern Asia, New Zealand and Australia, and is well-known as a vector of multiple zoonotic pathogens. This species exhibits two reproductive strategies, bisexual and obligate parthenogenetic reproduction. Hence, in the current study, the complete mitochondrial genomes of the bisexual and parthenogenetic populations were assembled and analyzed, and the expression of the mitochondrial protein-coding genes was evaluated and compared between the two reproductive populations. The results indicated that the length of the mitochondrial genomes of the two reproductive populations is 14,694 and 14,693 bp in the bisexual and parthenogenetic populations, respectively. The AT content in the mitochondrial genome of the bisexual and obligate parthenogenetic population reached 77.22 and 77.34%, respectively. The phylogenetic tree was constructed combining 13 protein-coding genes, which showed that the genetic distance between the bisexual and parthenogenetic populations was less than that between the subspecies. The expression of the mitochondrial protein-coding genes was quantitatively analyzed at different feeding status for the bisexual and parthenogenetic populations, and the results showed significant differences in the expression patterns of these genes, suggesting that they might trigger specific energy utilization mechanisms due to their different reproductive strategies and environmental pressures.

The mitochondrial genome of Muga silkworm (Antheraea assamensis) and its comparative analysis with other lepidopteran insects

Muga (Antheraea assamensis) is an economically important silkmoth endemic to the states of Assam and Meghalaya in India and is the producer of the strongest known commercial silk. However, there is a scarcity of genomic and proteomic data for understanding the organism at a molecular level. Our present study is on decoding the complete mitochondrial genome (mitogenome) of A. assamensis using next generation sequencing technology and comparing it with other available lepidopteran mitogenomes. Mitogenome of A. assamensis is an AT rich circular molecule of 15,272 bp (A+T content ~80.2%). It contains 37 genes comprising of 13 protein coding genes (PCGs), 22 tRNA and 2 rRNA genes along with a 328 bp long control region. Its typical tRNA^Met-tRNA^Ile-tRNA^Gln arrangement differed from ancestral insects (tRNA^Ile-tRNA^Gln-tRNA^Met). Two PCGs cox1 and cox2 were found to have CGA and GTG as start codons, respectively as reported in some lepidopterans. Interestingly, nad4l gene showed higher transversion mutations at intra-species than inter-species level. All PCGs evolved under strong purifying selection with highest evolutionary rates observed for atp8 gene while lowest for cox1 gene. We observed the typical clover-leaf shaped secondary structures of tRNAs with a few exceptions in case of tRNA^Ser1 and tRNA^Tyr where stable DHU and TΨC loop were absent. A significant number of mismatches (35) were found to spread over 19 tRNA structures. The control region of mitogenome contained a six bp (CTTAGA/G) deletion atypical of other Antheraea species and lacked tandem repeats. Phylogenetic position of A. assamensis was consistent with the traditional taxonomic classification of Saturniidae. The complete annotated mitogenome is available in GenBank (Accession No. KU379695). To the best of our knowledge, this is the first report on complete mitogenome of A. assamensis.

Mitochondrial Genomes of Two Bombycoidea Insects and Implications for Their Phylogeny

The mitochondrial genome (mt genome) provides important information for understanding molecular evolution and phylogenetics. As such, the two complete mt genomes of Ampelophaga rubiginosa and Rondotia menciana were sequenced and annotated. The two circular genomes of A. rubiginosa and R. menciana are 15,282 and 15,636 bp long, respectively, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The nucleotide composition of the A. rubiginosa mt genome is A + T rich (81.5%) but is lower than that of R. menciana (82.2%). The AT skew is slightly positive and the GC skew is negative in these two mt genomes. Except for cox1, which started with CGA, all other 12PCGs started with ATN codons. The A + T-rich regions of A. rubiginosa and R. menciana were 399 bp and 604 bp long and consist of several features common to Bombycoidea insects. The order and orientation of A. rubiginosa and R. menciana mitogenomes with the order trnM-trnI-trnQ-nad2 is different from the ancestral insects in which trnM is located between trnQ and nad2 (trnI-trnQ-trnM-nad2). Phylogenetic analyses indicate that A. rubiginosa belongs in the Sphingidae family, and R. menciana belongs in the Bombycidae family.

The first complete mitochondrial genome for the subfamily Limacodidae and implications for the higher phylogeny of Lepidoptera

The mitochondrial genome (mitogenome) provides important information for understanding molecular evolution and phylogeny. To determine the systematic status of the family Limacodidae within Lepidoptera, we infer a phylogenetic hypothesis based on the complete mitogenome of Monema flavescens (Lepidoptera: Limacodidae). The mitogenome of M. flavescens is 15,396 base pairs (bp), and includes 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region (CR). The AT skew of this mitogenome is slightly negative and the nucleotide composition is also biased towards A + T nucleotides (80.5%). All PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is initiated by CGA. All tRNAs display the typical clover-leaf structure characteristic of mitochondrial tRNAs, with the exception of trnS1 (AGN). The mitogenome CR is 401 bp and consists of several features common to Lepidoptera. Phylogenetic analysis using Bayesian Inference (BI) and Maximum Likelihood (ML) based on nucleotide and amino acid sequences of 13 mitochondrial PCGs indicates that M. flavescens belongs to Zygaenoidea. We obtain a well-supported phylogenetic tree consisting of Yponomeutoidea + (Tortricoidea + Zygaenoidea + (Papilionoidea + (Pyraloidea + (Noctuoidea + (Geometroidea + Bombycoidea))))).

Mitochondrial genome of the sweet potato hornworm, Agrius convolvuli (Lepidoptera: Sphingidae), and comparison with other Lepidoptera species

In the present study, we sequenced the complete mitochondrial genome (mitogenome) of Agrius convolvuli (Lepidoptera: Sphingidae) and compared it with previously sequenced mitogenomes of lepidopteran species. The mitogenome was a circular molecule, 15 349 base pairs (bp) long, containing 37 genes. The order and orientation of genes in the A. convolvuli mitogenome were similar to those in sequenced mitogenomes of other lepidopterans. All 13 protein-coding genes (PCGs) were initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which seemed to be initiated by the codon CGA, as observed in other lepidopterans. Three of the 13 PCGs had the incomplete termination codon T, while the remainder terminated with TAA. Additionally, the codon distributions of the 13 PCGs revealed that Asn, Ile, Leu2, Lys, Phe, and Tyr were the most frequently used codon families. All transfer RNAs were folded into the expected cloverleaf structure except for tRNA^Ser(AGN), which lacked a stable dihydrouridine arm. The length of the adenine (A) + thymine (T)-rich region was 331 bp. This region included the motif ATAGA followed by a 19-bp poly-T stretch and a microsatellite-like (TA)₈ element next to the motif ATTTA. Phylogenetic analyses (maximum likelihood and Bayesian methods) showed that A. convolvuli belongs to the family Sphingidae.

Comparative Mitochondrial Genome Analysis of Eligma narcissus and other Lepidopteran Insects Reveals Conserved Mitochondrial Genome Organization and Phylogenetic Relationships

In this study, we sequenced the complete mitochondrial genome of Eligma narcissus and compared it with 18 other lepidopteran species. The mitochondrial genome (mitogenome) was a circular molecule of 15,376 bp containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and an adenine (A) + thymine (T) − rich region. The positive AT skew (0.007) indicated the occurrence of more As than Ts. The arrangement of 13 PCGs was similar to that of other sequenced lepidopterans. All PCGs were initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which was initiated by the CGA sequence, as observed in other lepidopterans. The results of the codon usage analysis indicated that Asn, Ile, Leu, Tyr and Phe were the five most frequent amino acids. All tRNA genes were shown to be folded into the expected typical cloverleaf structure observed for mitochondrial tRNA genes. Phylogenetic relationships were analyzed based on the nucleotide sequences of 13 PCGs from other insect mitogenomes, which confirmed that E. narcissus is a member of the Noctuidae superfamily.

Mitochondrial genome characterization of Tecia solanivora (Lepidoptera: Gelechiidae) and its phylogenetic relationship with other lepidopteran insects

The complete mitogenome of the potato tuber moth Tecia solanivora (Lepidoptera: Gelechiidae) was sequenced, annotated, characterized and compared with 140 species of the order Lepidoptera. The circular genome is 15,251 bp, containing 37 genes (13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A+T-rich region). The gene arrangement was identical to other lepidopteran mitogenomes but different from the ancestral arrangement found in most insects for the tRNA-Met gene (A+T-region, tRNA-I, tRNA-Q, tRNA-M). The mitogenome of T. solanivora is highly A+T-biased (78.2%) and exhibits negative AT- and GC-skews. All PCGs are initiated by canonical ATN start codons, except for Cytochrome Oxidase subunit 1 (COI), which is initiated by CGA. Most PCGs have a complete typical stop codon (TAA). Only NAD1 has a TAG stop codon and the COII and NAD5 genes have an incomplete stop codon consisting of just a T. The A+T-rich region is 332 bp long and contains common features found in lepidopteran mitogenomes, including the 'ATAGA' motif, a 17 bp poly (T) stretch and a (AT)8 element preceded by the 'ATTTA' motif. Other tandem repeats like (TAA)4 and (TAT)7 were found, as well as (T)6 and (A)10 mononucleotide repeat elements. Finally, this mitogenome has 20 intergenic spacer regions. The phylogenetic relationship of T. solanivora with 28 other lepidopteran families (12 superfamilies) showed that taxonomic classification by morphological features coincides with the inferred phylogeny. Thus, the Gelechiidae family represents a monophyletic group, suggesting that T. solanivora and Pectinophora gossypiella have a recent common ancestor.

Characterization of the Complete Mitochondrial Genome of Cerura menciana and Comparison with Other Lepidopteran Insects

The complete mitochondrial genome (mitogenome) of Cerura menciana (Lepidoptera: Notodontidae) was sequenced and analyzed in this study. The mitogenome is a circular molecule of 15,369 bp, containing 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and a A+T-rich region. The positive AT skew (0.031) indicated that more As than Ts were present. All PCGs were initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which was initiated by CAG. Two of the 13 PCGs contained the incomplete termination codon T or TA, while the others were terminated with the stop codon TAA. The A+T-rich region was 372 bp in length and consisted of an ‘ATAGA’ motif followed by an 18 bp poly-T stretch, a microsatellite-like (AT)₈ and a poly-A element upstream of the trnM gene. Results examining codon usage indicated that Asn, Ile, Leu2, Lys, Tyr and Phe were the six most frequently occurring amino acids, while Cys was the rarest. Phylogenetic relationships, analyzed based on the nucleotide sequences of the 13 PCGs from other insect mitogenomes, confirmed that C. menciana belongs to the Notodontidae family.

The mitochondrial genome of the peach fruit fly, Bactrocera zonata (Saunders) (Diptera: Tephritidae): Complete DNA sequence, genome organization, and phylogenetic analysis with other tephritids using next generation DNA sequencing

Mitochondrial genome can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. The complete 15,935 bp mitochondrial genome of Bactrocera zonata (Diptera: Tephritidae), is assembled from Illumina MiSeq read data. The mitogenome information for B. zonata was compared to the homologous sequences of other tephritids. Annotation indicated that the structure and orientation of 13 protein coding genes (PCGs), 22 tRNA and 2 rRNA sequences were typical of, and similar to, the ten closely related tephritid species. The nucleotide composition shows heavily biased toward As and Ts accounting 73.34% and exhibits a slightly positive AT skew, which is similar to other known tephritid species. All PCGs are initiated by ATN codons, except for cox1 with TCG and atp8 with GTG. Nine PCGs use a common stop codon of TAA or TAG, whereas the remaining four use an incomplete termination codon T or TA likely to be completed by adenylation. All tRNAs have the typical clover-leaf structure, with an exception for trnS((AGN)). Four short intergenic spacers showed high degree of conservation among B. zonata and other ten tephritids. A poly(T) stretch at the 5' end followed by [TA(A)]n-like stretch and a tandem repeats of 39 bp has been observed in CR. The analysis of gene evolutionary rate revealed that the cox1 and atp6 exhibits lowest and highest gene substitution rates, respectively than other genes. The phylogenetic relationships based on Maximum Likelihood method using all protein-coding genes and two ribosomal RNA genes confirmed that B. zonata is closely related to Bactrocera correcta, Bactrocera carambolae, Bactrocera papayae, and Bactrocera philippinensis and Bactrocera dorsalis belonging to B. dorsalis species complex forms a monophyletic clade, which is in accordance with the traditional morphological classification and recent molecular works.

The complete mitochondrial genome of Antheraea pernyi strain Yu6 (Lepidoptera: Saturniidae)

In this study, the complete mitochondrial genome of Antheraea pernyi strain Yu6 (Lepidoptera: Saturniidae) has been reported for the first time. It is a circular molecule of 15,569 bp in length, containing 37 typical coding genes and 1 non-coding AT-rich region. The overall composition of the mitogenome is A (39.27%), G (7.71%), C (12.01%), and T (41.01%). Its gene order and content are identical to the common type found in most insect mitogenomes. All protein coding genes (PCGs) start with a typical ATN initiation codon, except for the cox1 gene, which begins with TTAG codon. Nine genes used standard complete termination codon TAA, whereas the cox1, cox2, nad3, and nad5 genes end with single T. All tRNAs display typical secondary cloverleaf structures as those of other insects. Additionally, the 552 bp long AT-rich region contained several structures common to the other lepidopterons, such as some structures of repeated motifs and microsatellite-like elements. Phylogenetic analysis showed that the Antheraea pernyi Yu6 was close to Saturniidae.

The complete mitochondrial genome of Antheraea pernyi strain 731 (Lepidoptera: Saturniidae)

In this study, the complete mitochondrial genome of Antheraea pernyi strain 731 (Lepidoptera: Saturniidae) is determined for the first time. It is a circular molecule of 15,570 bp in length, with 37 typical coding genes and one non-coding A T-rich region. Its gene order and content are identical to the common type found in most insect mitogenomes. All protein-coding genes (PCGs) start with a typical ATN initiation codon, except for the cox1 gene, which begins with TTAG codon. Nine genes used standard complete termination codon TAA, whereas the cox1, cox2, nad3, and nad5 genes end with single T. All tRNAs display typical secondary cloverleaf structures as those of other insects. Additionally, the non-coding AT-rich region is 553 bp long, located between rrnS and trnM genes. It contains some structures of repeated motifs and microsatellite-like elements characteristic of the other lepidopterons. Phylogenetic analysis showed that the Antheraea pernyi 731 was close to Saturniidae.

Morphological characters are compatible with mitogenomic data in resolving the phylogeny of nymphalid butterflies (lepidoptera: papilionoidea: nymphalidae)

Nymphalidae is the largest family of butterflies with their phylogenetic relationships not adequately approached to date. The mitochondrial genomes (mitogenomes) of 11 new nymphalid species were reported and a comparative mitogenomic analysis was conducted together with other 22 available nymphalid mitogenomes. A phylogenetic analysis of the 33 species from all 13 currently recognized nymphalid subfamilies was done based on the mitogenomic data set with three Lycaenidae species as the outgroups. The mitogenome comparison showed that the eleven new mitogenomes were similar with those of other butterflies in gene content and order. The reconstructed phylogenetic trees reveal that the nymphalids are made up of five major clades (the nymphaline, heliconiine, satyrine, danaine and libytheine clades), with sister relationship between subfamilies Cyrestinae and Biblidinae, and most likely between subfamilies Morphinae and Satyrinae. This whole mitogenome-based phylogeny is generally congruent with those of former studies based on nuclear-gene and mitogenomic analyses, but differs considerably from the result of morphological cladistic analysis, such as the basal position of Libytheinae in morpho-phylogeny is not confirmed in molecular studies. However, we found that the mitogenomic phylogeny established herein is compatible with selected morphological characters (including developmental and adult morpho-characters).

The complete mitochondrial genome of Bombyx mori strain Yu39 (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Bombyx mori strain Yu39 (Lepidoptera: Bombycidae) is a circular molecule of 15,652 bp in length, containing 37 typical mitochondrial genes: 13 protein-coding genes (PCGs), 22 transfer RNAs, 2 ribosomal RNAs and a non-coding AT-rich region. Its gene order and arrangement are identical to the common type found in most insect mitogenomes. All PCGs start with a typical ATN codon, except for the cox1 gene, which begins with uncertained codon. All PCGs terminate in the common stop codon TAA, except for the cox1 and cox2, which use single T as their stop codons. The non-coding AT-rich region is 494-bp long, located between rrnS and trnM genes. It contains some structures of repeated motifs and microsatellite-like elements characteristic of the other lepidopterons.

The complete mitochondrial genome of Bombyx mori strain Yu5 (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Bombyx mori strain Yu5 (Lepidoptera: Bombycidae) is a circular molecule of 15,644  bp in length, containing 37 typical coding genes and one non-coding AT-rich region. The overall composition of the mitogenome was A (43.07%), G (7.34%), C (11.33%), and T (38.26%). Its gene order and content were identical to the common type found in most insect mitogenomes. All protein coding genes (PCGs) started with a typical ATN initiation codon, except for the cox1 gene, which began with CGA codon. Eleven genes used standard complete termination codon TAA, whereas the cox1 and cox2 genes ended with single T. All tRNA genes displayed typical secondary cloverleaf structures as those of other insects. Additionally, the 494  bp long AT-rich region contained several structures common to the other lepidopterons, such as some structures of repeated motifs and microsatellite-like elements and a poly-A element upstream of transfer RNA M (trnM) gene.

The complete mitochondrial genome of Bombyx mori strain Chunyun (Lepidoptera: Bombycidae)

In this study, the complete mitochondrial genome sequence of Bombyx mori strain Chunyun (Lepidoptera: Bombycidae) has been reported for the first time. The total length of the mitogenome was 15,659 bp with 37 typical animal mitochondrial genes and 1 non-coding A + T-rich region. Its gene content and order were identical to those of other lepidopteran mitochondrial genomes. All protein-coding genes (PCGs) were initiated by ATN codons except for the COI gene, which began with CGA codon. Eleven PCGs stopped with termination codon TAA, whereas the COI and COII genes ended with single T. All tRNAs have typical structures of insect mitochondrial tRNAs. The 494 bp AT-rich region contains several features common to other lepidopterans, such as some structures of repeated motifs and microsatellite-like elements and a poly-A element upstream of transfer RNA M (trnM) gene.

The complete mitochondrial genome of Bombyx mori strain Jin6 (Lepidoptera: Bombycidae)

In this study, the complete mitochondrial genome of Bombyx mori strain Jin6 (Lepidoptera: Bombycidae) has been reported for the first time. It was a circular molecule of 15,648 bp in length, containing 37 typical coding genes and one non-coding AT-rich region. The overall composition of the mitogenome was A (43.05%), G (7.30%), C (11.35%), and T (38.29%). Its gene order and content were identical to the common type found in most insect mitogenomes. All protein coding genes (PCGs) started with a typical ATN initiation codon, except for the cox1 gene, which began with CGA codon. Eleven genes used standard complete termination codon TAA, whereas the cox1 and cox2 genes ended with single T. All tRNA genes displayed typical secondary cloverleaf structures as those of other insects. Additionally, the 494 bp long AT-rich region contained several structures common to the other lepidopterons, such as some structures of repeated motifs and microsatellite-like elements.

The complete mitochondrial genome of Bombyx mori strain Huayu (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Bombyx mori strain Huayu (Lepidoptera: Bombycidae) is determined in this study. The genome was 15,666 bp long, with 37 typical animal mitochondrial genes and 1 non-coding A + T-rich region. Its gene content and order were identical to those of other lepidopteran mitochondrial genomes. All protein-coding genes (PCGs) were initiated by ATN codons except for the COI gene, which began with uncertained codon. Eleven PCGs stopped with termination codon TAA, whereas the COI and COII genes ended with single T. All tRNAs have typical structures of insect mitochondrial tRNAs. The 494 bp AT-rich region contains several features common to other lepidopterans, such as the motif ATAGA followed by an 18 bp poly-T stretch and an 11 bp poly-A element upstream of transfer RNA M (trnM) gene.

The complete mitochondrial genome of the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

The complete mitochondrial genome (mitogenome) of Plutella xylostella (Lepidoptera: Plutellidae) was determined (GenBank accession No. KM023645). The length of this mitogenome is 16,014 bp with 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and an A + T-rich region. It presents the typical gene organization and order for completely sequenced lepidopteran mitogenomes. The nucleotide composition of the genome is highly A + T biased, accounting for 81.48%, with a slightly positive AT skewness (0.005). All PCGs are initiated by typical ATN codons, except for the gene cox1, which uses CGA as its start codon. Some PCGs harbor TA (nad5) or incomplete termination codon T (cox1, cox2, nad2 and nad4), while others use TAA as their termination codons. The A + T-rich region is located between rrnS and trnM with a length of 888 bp.

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A+T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA(Ser)(AGN). The A. atlas A+T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.

Genome-wide characterization of adaptation and speciation in tiger swallowtail butterflies using de novo transcriptome assemblies

Hybrid speciation appears to be rare in animals, yet characterization of possible examples offers to shed light on the genomic consequences of this unique phenomenon, as well as more general processes such as the role of adaptation in speciation. Here, we first generate transcriptome assemblies for a putative hybrid butterfly species, Papilio appalachiensis, its parental species, P. glaucus and P. canadensis, and an outgroup, P. polytes. Then, we use these data to infer genome-wide patterns of introgression and genomic mosaicism using both phylogenetic and population genetic approaches. Our results reveal that there is little genetic divergence among all three of the focal species, but the subset of gene trees that strongly support a specific tree topology suggest widespread sharing of genetic variation between P. appalachiensis and both parental species, likely as a result of hybrid speciation. We also find evidence for substantial shared genetic variation between P. glaucus and P. canadensis, which may be due to gene flow or ancestral variation. Consistent with previous work, we show that P. applachiensis is more similar to P. canadensis at Z-linked genes and more similar to P. glaucus at mitochondrial genes. We also identify a variety of targets of adaptive evolution, which appear to be enriched for traits that are likely to be important in the evolution of this butterfly system, such as pigmentation, hormone sensitivity, developmental processes, and cuticle formation. Overall, our results provide a genome-wide portrait of divergence and introgression associated with adaptation and speciation in an iconic butterfly radiation.

Biogeography revealed by mariner-like transposable element sequences via a Bayesian coalescent approach

Genetic diversity of natural populations is useful in biogeographical studies. Here, we apply a Bayesian method based on the coalescent model to dating biogeographical events by using published DNA sequences of wild silkworms, Bombyx mandarina, and the domesticated model organisms B. mori, both of which categorized into the order of Lepidoptera, sampled from China, Korea, and Japan. The sequences consist of the BmTNML locus and the flanking intergenic regions. The BmTNML locus is composed of cecropia-type mariner-like element (MLE) with inverted terminal repeats, and three different transposable elements (TE), including L1BM, BMC1 retrotransposons, and BmamaT1, are inserted into the MLE. Based on the genealogy defined by TE insertions/deletions (indels), we estimated times to the most recent common ancestor and these indels events using the flanking, MLE, and indels sequences, respectively. These estimates by using MLE sequences strongly correlated with those by using flanking sequences, implying that cecropia-type MLEs can be used as a molecular clock. MLEs are thought to have transmitted horizontally among different species. By using a pair of published cecropia-type MLE sequences from lepidopteran insect, an emperor moth, and a coral in Ryukyu Islands, we demonstrated dating of horizontal transmission between species which are distantly related but inhabiting geographically close region.

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A+T-rich region. The A+T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.

The complete mitochondrial genome of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae) was determined. The circular genome is 15,345 bp long, and presents a typical gene organization and order for sequenced mitogenomes of Bombycidea species. The nucleotide composition of the genome is highly A+T biased, accounting for 79.86%. The AT skew of the genome is slightly negative, indicating the occurrence of more Ts than As, as found in other Saturniidae species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI and COII, which are tentatively designated by CGA and GTG, respectively, as observed in other insects. Four of 13 PCGs, including COI, COII, ATP6, and ND3, harbor the incomplete termination codons, T or TA. With an exception for tRNASer(AGN), all other tRNAs can form a typical clover-leaf structure of mitochondrial tRNA. The 359 bp A+T-rich region of S. c. cynthia contains non-repetitive sequences, but harbors several features common to the Bombycidea insects, including the motif ATAGA followed by a poly-T stretch of 19 bp, a microsatellite-like (AT)7 element preceded by the ATTTA motif, and a poly-A element upstream tRNAMet. The phylogenetic analyses support the morphology-based current hypothesis that Bombycidae and Saturniidae are monophyletic. Our result confirms that Saturniini and Attacini form a reciprocal monophyletic group within Saturniidae.

Little gene flow between domestic silkmoth Bombyx mori and its wild relative Bombyx mandarina in Japan, and possible artificial selection on the CAD gene of B. mori

We analyzed PCR-amplified carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase (CAD) gene fragments from 146 Bombyx mori native strains and found extremely low levels of DNA polymorphism. Two haplotypes were identified, one of which was predominant. CAD haplotype analysis of 42 samples of Japanese B. mandarina revealed four haplotypes. No common haplotype was shared between the two species and at least five base substitutions were detected. This result was suggestive of low levels of gene flow between the two species. The nucleotide diversity (π) scores of the two samples differed markedly: lower π values were estimated for B. mori native strains than Japanese B. mandarina. We further analyzed 12 Chinese B. mandarina derived from seven areas of China, including Taiwan. The results clearly indicated that the π score was ~80-fold greater in Chinese B. mandarina than in B. mori. The extremely low level of DNA polymorphism in B. mori compared to its wild relatives suggested that the CAD gene itself or its tightly linked regions are possible targets for silkworm domestication.

Complete sequence of the mitochondrial genome of the Japanese buff-tip moth, Phalera flavescens (Lepidoptera: Notodontidae)

We sequenced the complete mitochondrial genome of Phalera flavescens. The mitogenome is 15,659 bp in length, including 13 protein-coding genes (atp6, atp8, cox1-3, nad1-6, nad4L, cob), two ribosomal RNAs (rrnS and rrnL), 22 transfer RNAs and an AT-rich region, a putative control region (D-loop). Gene order and orientation were found to be identical to those of other completely sequenced lepidopteran mitogenomes. All 13 protein-coding genes start with the common codon ATN, except for the cox1 gene, which uses CGA as the initial codon. Nine of the 13 protein-coding genes stop with codon TAA, while the cox1, cox2, nad5, and nad4 genes stop with the single nucleotide T. All tRNA genes can be folded into canonical cloverleaf secondary structure, except for trnS1, which loses the ''DHU'' arm. Six overlapping sequences totaling 20 bp (1-8 bp for each sequence) and 16 intergenic spacer sequences, totaling 276 bp (1-58 bp for each sequence) are scattered throughout the genome; the largest intergenic spacer is located between the trnQ and nad2 genes. A microsatellite-like structure (AT)(6)ACC(AT)(6) and 16-bp poly-T elements preceded by the ATTTA motif are present in the D-loop region. Additionally, unexpectedly, an extra 190-bp insertion, with unknown function, was found in the small subunit rRNA gene (rrnS); this gene is the longest known (1020 bp) among all of the Lepidoptera.

The complete mitogenome of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) and comparison with other lepidopteran insects

The complete mitochondrial genome (mitogenome) of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) was determined to be 15,653bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A+T-rich region. It has the typical gene organization and order of mitogenomes from lepidopteran insects. 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 subunit 1 (cox1) gene which was initiated by CGA. The cox1 and cox2 genes had incomplete stop codons consisting of just a T. All the tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome was 495bp in length and consisted of several features common to the lepidopteras. Phylogenetic analysis showed that the B. mori Dazao was close to Bombycidae.

The complete mitochondrial genome of Leucoptera malifoliella Costa (Lepidoptera: Lyonetiidae)

The mitochondrial genome (mitogenome) of Leucoptera malifoliella (=L. scitella) (Lepidoptera: Lyonetiidae) was sequenced. The size was 15,646 bp with gene content and order the same as those of other lepidopterans. The nucleotide composition of L. malifoliella mitogenome is highly A+T biased (82.57%), ranked just below Coreana raphaelis (82.66%) (Lepidoptera: Lycaenidae). All protein-coding genes (PCGs) start with the typical ATN codon except for the cox1 gene, which uses CGA as the initiation codon. Nine PCGs have the common stop codon TAA, four PCGs have the common stop codon T as incomplete stop codons, and nad4l and nad6 have TAG as the stop codon. Cloverleaf secondary structures were inferred for 22 tRNA genes, but trnS1(AGN) was found to lack the DHU stem. The secondary structure of rrnL and rrnS is generally similar to other lepidopterans but with some minor differences. The A+T-rich region includes the motif ATAGA, but the poly (T) stretch is replaced by a stem-loop structure, which may have a similar function to the poly (T) stretch. Finally, there are three long repeat (154 bp) sequences followed by one short repeat (56 bp) with four (TA)(n) intervals, and a 10-bp poly-A is present upstream of trnM. Phylogenetic analysis shows that the position of Yponomeutoidea, as represented by L. malifoliella, is the same as traditional classifications. Yponomeutoidea is the sister to the other lepidopteran superfamilies covered in the present study.

Nucleotide sequence variation in mitochondrial COI gene among 147 silkworm (Bombyx mori) strains from Japanese, Chinese, European and moltinism classes

We characterized the nucleotide sequences of PCR-amplified mitochondrial COI fragments of 147 silkworm (Bombyx mori) strains that have been maintained in the National Institute of Agrobiological Sciences. Coding sequences (714 bp) of the 147 COI fragments were classified into eight haplotypes based on nucleotide differences at eight segregating sites. No length variation was identified in this region. The 5'-noncoding region showed different features, wherein changes in the number of Ts in the T-stretch, together with two base substitutions, were observed. As a result, the 147 COI noncoding sequences were classified into six haplotypes. Combining the coding and noncoding regions, we identified 14 haplotypes. One of the 14 haplotypes, Hap1A was exclusively abundant in the Japanese native strain class, while this haplotype was less frequent in the other three native strain classes. This finding suggests that the Japanese strain class underwent significant genetic differentiation from the Chinese, European, and moltinism classes, when the each class is regarded as a population. Comparison of the nucleotide sequences to those of B. mandarina (which inhabits Japan) revealed changes that are significantly larger than those within either B. mori or B. mandarina. Furthermore, we detected no common haplotypes between them, which suggests the concept of suppressed gene flow between the two species.

Phylogeny and evolutionary history of the silkworm

The silkworm, Bombyx mori, played an important role in the old Silk Road that connected ancient Asia and Europe. However, to date, there have been few studies of the origins and domestication of this species using molecular methods. In this study, DNA sequences of mitochondrial and nuclear loci were used to infer the phylogeny and evolutionary history of the domesticated silkworm and its relatives. All of the phylogenetic analyses indicated a close relationship between the domesticated silkworm and the Chinese wild silkworm. Domestication was estimated to have occurred about 4100 years ago (ya), and the radiation of the different geographic strains of B. mori about 2000 ya. The Chinese wild silkworm and the Japanese wild silkworm split about 23600 ya. These estimates are in good agreement with the fossil evidence and historical records. In addition, we show that the domesticated silkworm experienced a population expansion around 1000 ya. The divergence times and the population dynamics of silkworms presented in this study will be useful for studies of lepidopteran phylogenetics, in the genetic analysis of domestic animals, and for understanding the spread of human civilizations.

The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

Background

Lepidoptera encompasses more than 160,000 described species that have been classified into 45–48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to the non-ditrysian lineage Hepialoidea and conducted a comparative mitogenomic analysis across Lepidoptera.

Results

The mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,816 bp long with an A + T content of 81.28 % and 82.34 %, respectively. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and the A + T-rich region. Different tandem repeats in the A + T-rich region mainly account for the size difference between the two mitogenomes. All the protein-coding genes start with typical mitochondrial initiation codons, except for cox1 (CGA) and nad1 (TTG) in both mitogenomes. The anticodon of trnS(AGN) in T. renzhiensis and T. yunnanensis is UCU instead of the mostly used GCU in other sequenced Lepidoptera mitogenomes. The 1,584-bp sequence from rrnS to nad2 was also determined for an unspecified ghost moth (Thitarodes sp.), which has no repetitive sequence in the A + T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A + T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif ‘ATAGA’ and poly-T stretch in the A + T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes.

Conclusion

The mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event(s) likely occurred after Hepialoidea diverged from other lepidopteran lineages. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.

The complete mitochondrial genome of the wild silkworm moth, Actias selene

The complete mitochondrial genome (mitogenome) of Actias selene (Lepidoptera: Saturniidae) was determined to be 15,236 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The arrangement of 13 PCGs was similar to that of other sequenced lepidopterans. The AT skew of the mitogenome of A. selene was slightly negative, indicating a higher number of T compared to A nucleotides. The nucleotide composition of the mitogenome of A. selene was also biased toward A+T nucleotides (78.91%). All PCGs were initiated by ATN codons, except for the gene encoding cytochrome c oxidase subunit 1 (cox1), which may be initiated by the TTAG, as observed in other lepidopterans. Three genes, including cox1, cox2, and nad5, had incomplete stop codons consisting of just a T. With an exception for trnS1(AGN), all the other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome of A. selene was 339 bp in length, and contains several features common to the Lepidopteras, including non-repetitive sequences, a conserved structure combining the motif ATAGA and an 18-bp poly-T stretch and a poly-A element upstream of trnM gene. Phylogenetic analysis showed that A. selene was close to Saturniidae.

The Complete Mitochondrial Genome of Ctenoptilum vasava (Lepidoptera: Hesperiidae: Pyrginae) and Its Phylogenetic Implication

We here report the first complete mitochondrial (mt) genome of a skipper, Ctenoptilum vasava Moore, 1865 (Lepidoptera: Hesperiidae: Pyrginae). The mt genome of the skipper is a circular molecule of 15,468 bp, containing 2 ribosomal RNA genes, 24 putative transfer RNA (tRNA), genes including an extra copy of trnS (AGN) and a tRNA-like insertion trnL (UUR), 13 protein-coding genes and an AT-rich region. All protein-coding genes (PCGs) are initiated by ATN codons and terminated by the typical stop codon TAA or TAG, except for COII which ends with a single T. The intergenic spacer sequence between trnS (AGN) and ND1 genes also contains the ATACTAA motif. The AT-rich region of 429 bp is comprised of nonrepetitive sequences, including the motif ATAGA followed by an 19 bp poly-T stretch, a microsatellite-like (AT)(3) (TA)(9) element next to the ATTTA motif, an 11 bp poly-A adjacent to tRNAs. Phylogenetic analyses (ML and BI methods) showed that Papilionoidea is not a natural group, and Hesperioidea is placed within the Papilionoidea as a sister to ((Pieridae + Lycaenidae) + Nymphalidae) while Papilionoidae is paraphyletic to Hesperioidea. This result is remarkably different from the traditional view where Papilionoidea and Hesperioidea are considered as two distinct superfamilies.

The complete mitochondrial genome of the rice moth, Corcyra cephalonica

The complete mitochondrial genome (mitogenome) of the rice moth, Corcyra cephalonica Stainton (Lepidoptera: Pyralidae) was determined as a circular molecular of 15,273 bp in size. The mitogenome composition (37 genes) and gene order are the same as the other lepidopterans. Nucleotide composition of the C. cephalonica mitogenome is highly A+T biased (80.43%) like other insects. Twelve protein-coding genes start with a typical ATN codon, with the exception of coxl gene, which uses CGA as the initial codon. Nine protein-coding genes have the common stop codon TAA, and the nad2, cox1, cox2, and nad4 have single T as the incomplete stop codon. 22 tRNA genes demonstrated cloverleaf secondary structure. The mitogenome has several large intergenic spacer regions, the spacer1 between trnQ gene and nad2 gene, which is common in Lepidoptera. The spacer 3 between trnE and trnF includes microsatellite-like repeat regions (AT)18 and (TTAT)(3). The spacer 4 (16 bp) between trnS2 gene and nad1 gene has a motif ATACTAT; another species, Sesamia inferens encodes ATCATAT at the same position, while other lepidopteran insects encode a similar ATACTAA motif. The spacer 6 is A+T rich region, include motif ATAGA and a 20-bp poly(T) stretch and two microsatellite (AT)(9), (AT)(8) elements.

The complete mitochondrial genome of the oriental fruit moth Grapholita molesta (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae) currently is one of the economically most destructive pest species of stone and pome fruits worldwide. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,776 bp long, with an A + T content of 81.24%, containing 37 typical animal mitochondrial genes and an A + T-rich region. All gene are arranged as hypothesized ancestral gene order of insects except for trnM, which was shuffled from 3′ downstream of trnQ to 5′ upstream of trnI. cox1 gene uses unusual CGA start codon, as that in all other sequenced lepidopteran mitochondrial genome. The secondary structures for the two rRNA genes were predicted. All helices typically present in insect mitochondrial rRNA genes are generated. A microsatellite sequence was inserted into the region of H2347 in rrnL in G. molesta and two other sequenced tortricid mitochondrial genomes, indicating that the insertion event in this helix might occurred anciently in family Tortricidae. All of the 22 typical animal tRNA genes have a typical cloverleaf structure except for trnS2, in which the D-stem pairings in the DHU arm are absent. An intergenic sequence is present between trnQ and nad2 as well as in other sequenced lepidopteran mitochondrial genomes, which was presumed to be a remnant of trnM gene and its boundary sequences after the duplication of trnM to the upstream of trnI in Lepidoptera. The A + T-rich region is 836 bp, containing six repeat sequences of “TTATTATTATTATTAAATA(G)TTT.”

Mitochondrial genome of the cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) and comparison with other Lepidopterans

We present the complete sequence of the mitochondrial genome (mitogenome) of the cotton bollworm Helicoverpa armigera. The 15,347-bp mitogenome of H. armigera was arranged in the same order described for all other sequenced lepidopterans, which differs from the most common type found in insects, due to the movement of trnM to a position 5'-upstream of trnI. The gene overlap in the H. armigera mitogenome is totally 23 bp in six locations. The H. armigera mitogenome has a total of 175 bp of intergenic spacer sequences spread over 14 regions ranging in size from 1 to 45 bp. The nucleotide composition of the whole mitogenome of H. armigera is highly A+T biased, accounting for 80.97%, with a slightly positive AT skewness and negative GC skewness, indicating the occurrence of more A than T, C more than G. The protein-encoding genes have typical mitochondrial start codons, except for cox1, which contains the unusual CGA. The cox1, cox2, and nad4 genes have incomplete stop codons (T). The lrRNA and srRNA genes are 1395 and 794-bp long, respectively. All tRNAs have a typical cloverleaf structure of mitochondrial tRNAs, except for trnS1(AGN), the dihydrouridine arm of which could not form a stable stem-loop structure. The H. armigera A+T-rich region contains a conserved structure combining the motif ATAGA and a 19-bp poly-T stretch, but absence of the 9-bp poly-A element upstream of trnM.

Mitochondrial genome sequence and expression profiling for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

We report the assembly of the 14,054 bp near complete sequencing of the mitochondrial genome of the legume pod borer (LPB), Maruca vitrata (Lepidoptera: Crambidae), which we subsequently used to estimate divergence and relationships within the lepidopteran lineage. The arrangement and orientation of the 13 protein-coding, 2 rRNA, and 19 tRNA genes sequenced was typical of insect mitochondrial DNA sequences described to date. The sequence contained a high A+T content of 80.1% and a bias for the use of codons with A or T nucleotides in the 3rd position. Transcript mapping with midgut and salivary gland ESTs for mitochondrial genome annotation showed that translation from protein-coding genes initiates and terminates at standard mitochondrial codons, except for the coxI gene, which may start from an arginine CGA codon. The genomic copy of coxII terminates at a T nucleotide, and a proposed polyadenylation mechanism for completion of the TAA stop codon was confirmed by comparisons to EST data. EST contig data further showed that mature M. vitrata mitochondrial transcripts are monocistronic, except for bicistronic transcripts for overlapping genes nd4/nd4L and nd6/cytb, and a tricistronic transcript for atp8/atp6/coxIII. This processing of polycistronic mitochondrial transcripts adheres to the tRNA punctuated cleavage mechanism, whereby mature transcripts are cleaved only at intervening tRNA gene sequences. In contrast, the tricistronic atp8/atp6/coxIII in Drosophila is present as separate atp8/atp6 and coxIII transcripts despite the lack of an intervening tRNA. Our results indicate that mitochondrial processing mechanisms vary between arthropod species, and that it is crucial to use transcriptional information to obtain full annotation of mitochondrial genomes.

The complete mitochondrial genome of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae)

The complete mitochondrial genome (mitogenome) of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae) was determined. The genome is a circular molecule 15 481 bp long. It presents a typical gene organization and order for completely sequenced lepidopteran mitogenomes, but differs from the insect ancestral type for the placement of tRNA(Met). The nucleotide composition of the genome is also highly A + T biased, accounting for 80.38%, with a slightly positive AT skewness (0.010), indicating the occurrence of more As than Ts, as found in the Noctuoidea species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI, which is tentatively designated by the CGA codon as observed in other lepidopterans. Four of 13 PCGs harbor the incomplete termination codon, T or TA. All tRNAs have a typical clover-leaf structure of mitochondrial tRNAs, except for tRNA(Ser)(AGN), the DHU arm of which could not form a stable stem-loop structure. The intergenic spacer sequence between tRNA(Ser)(AGN) and ND1 also contains the ATACTAA motif, which is conserved across the Lepidoptera order. The H. cunea A+T-rich region of 357 bp is comprised of non-repetitive sequences, but harbors several features common to the Lepidoptera insects, including the motif ATAGA followed by an 18 bp poly-T stretch, a microsatellite-like (AT)(8) element preceded by the ATTTA motif, an 11 bp poly-A present immediately upstream tRNA(Met). The phylogenetic analyses support the view that the H. cunea is closerly related to the Lymantria dispar than Ochrogaster lunifer, and support the hypothesis that Noctuoidea (H. cunea, L. dispar, and O. lunifer) and Geometroidea (Phthonandria atrilineata) are monophyletic. However, in the phylogenetic trees based on mitogenome sequences among the lepidopteran superfamilies, Papillonoidea (Artogeia melete, Acraea issoria, and Coreana raphaelis) joined basally within the monophyly of Lepidoptera, which is different to the traditional classification.

Genetic diversity, molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes

Background

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.

Results

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.

Conclusions

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.