Complete mitochondrial genome of theGrapholita dimorphaKomai (Lepidoptera: Tortricidae)

Complete mitochondrial genome analysis and molecular phylogenetic implications of Kennelia xylinana (Lepidoptera: Tortricidae)

Kennelia is a small genus in Tortricidae that is distributed in the Oriental and Palaearctic regions, and its taxonomic position within the subfamily Olethreutinae is controversial. For a comprehensive understanding of the genus, we sequenced the mitogenome of Kennelia xylinana, the type species of Kennelia, and Ancylis unculana, a species of Enarmoniini; analyzed the mitogenome characteristics of K. xylinana; and explored its phylogenetic position. Similar to other members of Lepidoptera, the mitogenome of K. xylinana is 15,762-bp long and consists of 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, and a noncoding control region. In particular, we found a structure (TATAATTAATAA)₁₁ in the middle of the AT-rich region. Based on the Bayesian inference and maximum likelihood analyses of the 13 PCGs of 40 tortricid species, representing 8 tribes of 2 subfamilies, K. xylinana was clustered with two members of Enarmoniini, A. unculana and Loboschiza koenigiana, and formed highly supported monophyly. The results indicate that Kennelia should be placed in the tribe Enarmoniini.

Reads in a haystack: extracting complete mitogenome sequences hidden in baculovirus datasets

Thaumatotibia leucotreta (Lepidoptera, Tortricidae) is one of many economically important insect pests for which no complete mitogenome sequence is available. The complete mitochondrial sequences for this species and other key pests could assist in the development of novel molecular techniques, such as enabling the identification of population-specific markers which could assist in improved monitoring of populations. The objective of this study was to determine whether NGS datasets generated for entomopathogenic viruses contain reads originating from host mitochondrial DNA. A total of 28 NGS datasets generated for the baculovirus Cryptophlebia leucotreta granulovirus (CrleGV) were analysed in this study. Three datasets contained sufficient reads providing adequate coverage for the assembly of complete mitogenomes. All 13 protein-coding genes, 22 tRNAs and both rRNAs present in the mitogenomes of other species within the Grapholitini tribe, were identified. Phylogenetic analysis of the mitogenomes at both an intrafamilial and interspecies level grouped the sequences within the Olethreutinae and T. leucotreta clades, respectively. Analysis of single nucleotide variations (SNVs) between each T. leucotreta sequence indicated up to 75 differences across the mitogenome. The methodology used in this study could be expanded to other baculovirus NGS datasets enabling the generation of novel lepidopteran mitogenome sequences.

The complete mitochondrial genome of Cochylimorpha cultana (Lederer, 1855) (Lepidoptera: Tortricidae)

The complete mitochondrial genome of Cochylimorpha cultana (Lederer) (Lepidoptera: Tortricidae) was 15,348 bps in size, and comprised 37 genes, which were 13 PCGs, 22 tRNA genes and two rRNA genes. Most PCGs used the conventional ATN start codon, except for cox1 initiating with CGA. Four genes (cox1, cox2, nad4 and nad5) used single T residue as stop condon. 21 out of 22 tRNAs are folded into the cloverleaf secondary structure, except for trnS1. The phylogenetic analysis based on maximum-likelihood (ML) method revealed that the evolutionary status of C. cultana in Tortricinae at the molecular level, which agrees well with the classical taxonomy.

The complete mitochondrial genome of Trifida elongate and comparative analysis of 43 leafhoppers

In this study, the mitochondrial genome of Trifida elongate was sequenced, and comparative analysis of T. elongate and other 43 leafhoppers was performed based on the mitochondrial genome. The mitochondrial genome sequence length of T. elongate was 14,924 bp. It comprised 13 protein-coding genes, 2 ribosomal RNA (rRNA) genes, 22 transport RNA (tRNA) genes, and 1 non-coding control region. The control region is located between the rrnS and trnI genes, is characterized by two tandem repeats and three simple sequence repeats. Phylogenetic analysis showed that T. elongate is closely related with Bolanusoides shaanxiensis and Limassolla lingchuanensis (bootstrap value = 92% and posterior probabilities = 1). Analysis of synonymous and non-synonymous nucleotide substitutions showed that Ka/Ks value of the 13 protein-coding genes of 8 subfamily leafhoppers were less than 1 ranging from 0.0315 to 0.9928. atp8 had the highest Ka/Ks value whereas cox1 had the lowest Ka/Ks value. This study provides information on the structure and sequence characteristics of the mitochondrial genome of T. elongata. Typhlocybinae is clustered with (Cicadellinae+(Idiocerinae+(Mileewinae+(Nirvaninae+(Evacanthinae+Ledrinae))))).

Multiple refugia from penultimate glaciations in East Asia demonstrated by phylogeography and ecological modelling of an insect pest

BACKGROUND

Refugial populations in Quaternary glaciations are critical to understanding the evolutionary history and climatic interactions of many extant species. Compared with the well-studied areas of Europe and Northern America, refugia of species in eastern Asia remain largely unknown. Here, we investigated the phylogeographic history of a globally important insect pest, the oriental fruit moth Grapholita molesta, in its native range of China.

RESULTS

Genetic structure analyses unveiled three distinct groups and a set of populations with admixture. Approximate Bayesian Computation (ABC) analyses support range expansion of this moth from southwest groups of Yunnan and Sichuan to northern and eastern China. A set of admixed populations was found around these two ancestral groups. This pattern of genetic structure points to two refugia located in the Yunnan region and Sichuan Basin. The split of the two refugia was dated to 329.2 thousand years ago in the penultimate glacial period. One of the lineages was exclusively found around the Sichuan Basin, indicating the formation of endemic populations in this refugium. Ecological niche model analysis suggested a shrinking distribution from the LIG period to the MID period in the Sichuan lineage but a wide and stable distribution in the other lineage.

CONCLUSIONS

Our results for the first time suggest that Yunnan and Sichuan jointly served as two large-scale refugia in eastern Asia in Quaternary glaciations, helping to maintain genetic diversity overall.

Convergent herbivory on conifers by Choristoneura moths after boreal forest formation

Mitogenomes are useful markers for phylogenetic studies across a range of taxonomic levels. Here, we focus on mitogenome variation across the tortricid moth genus Choristoneura and particularly the spruce budworm (Choristoneura fumiferana) species complex, a notorious pest group of North American conifer forests. Phylogenetic relationships of Tortricidae, representing two subfamilies, four tribes and nine genera, were analyzed using 21 mitogenomes. These included six newly-sequenced mitogenomes for species in the spruce budworm complex plus three additional Choristoneura species and 12 previously published mitogenomes from other tortricids and one from the Cossidae. We evaluated the phylogenetic informativeness of the mitogenomes and reconstructed a time-calibrated tree with fossil and secondary calibrations. We found that tortricid mitogenomes had conserved protein and ribosomal regions, and analysis of all protein-coding plus ribosomal genes together provided an efficient marker at any taxonomic rank. The time-calibrated phylogeny showed evolutionary convergence of conifer feeding within Choristoneura, with two independent lineages, the Nearctic spruce budworm complex and the Palearctic species Choristoneura murinana, both shifting onto conifers about 11 million years ago from angiosperms. These two host-plant shifts both occurred after the formation of boreal forest in the late Miocene. Haplotype diversification within the spruce budworm complex occurred in the last 4 million years, and is probably linked to the initial cooling cycles of the Northern Hemisphere in the Pliocene.

The mitochondrial genome of the lepidopteran host cadaver (Thitarodes sp.) of Ophiocordyceps sinensis and related phylogenetic analysis

To understand the phylogeny of the host insect (Thitarodes sp.) of the fungus Ophiocordyceps sinensis, we sequenced, annotated and characterized the complete mitochondrial (mt) genome of the host cadaver of a natural O. sinensis. Further, we compared the Thitarodes sp. mt genome with those of the other 7 sequenced Hepialidae and examined the phylogenetic relationships using a constructed Maximum Likelihood (ML) phylogenetic tree and mt genomic features (genetic distances and intergenic spacers). The mt genome is a circular molecule of 16,280bp in length with a high A+T content (81.20%) and contains 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and an AT-rich region. The gene arrangement is identical to the ancestral arrangement but differs from those of other lepidopteran mt genomes because of the arrangement of tRNA genes. The tRNA region, which is located between the AT-rich region and nad2, is trnI/trnQ/trnM (IQM) in Thitarodes sp., rather than the trnM/trnI/trnQ (MIQ) of the Lepidoptera-specific rearrangement. All PCGs begin with the canonical start codons ATN or NTG, except for cox1, which starts with CGA. Most PCGs terminate with the typical stop codon TAA, although some have an incomplete stop codon (T). The 1473bp AT-rich region is located between the rrnS (12S rRNA) and trnI, which is the longest sequenced in a Thitarodes mt genome to date, containing nine 112bp copies and one partial copy of a 55bp sequence. The results derived from the phylogenetic tree, the genetic distances and the intergenic spacers of the mt genome show that the host insect of O. sinensis belongs to the Thitarodes, while Endoclita signifer and Napialus hunanensis form a relatively distinct lineage from Thitarodes. The sequence and full annotation of this moth mt genome will provide more molecular information about the Exoporia within the Lepidoptera, and the clarification of its phylogeny will improve the management of this insect resource and the conservation and sustainable use of this endangered medicinal species in China.

The Complete Mitochondrial Genome of Brachmia macroscopa (Lepidoptera: Gelechiidae) and Its Related Phylogenetic Analysis

The sweet potato leaf folder, Brachmia macroscopa, is an important pest in China. The complete mitogenome, which consists of 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region, was sequenced and found to be 15,394 bp in length (GeneBank no. KT354968). The gene order and orientation of the B. macroscopa mitogenome were similar to those of other sequenced lepidopteran species. All of the PCGs started with ATN as the canonical start codon except for cox1, which started with CGA. In regard to stop codons, most PCGs stopped at TAA except for cox2, which stopped at TA, and nad4, which stopped at a single T. Thirteen PCGs of the available species (33 taxa) were used to demonstrate phylogenetic relationships. The ditrysian cluster was supported as a monophyletic clade at high levels by using maximum likelihood and Bayesian methods. The apoditrysian group, covering the Gelechioidea, formed a monophyletic clade with a bootstrap value of 88% and a posterior probability of 1.00. The superfamily Gelechioidea was supported as a monophyletic lineage by a posterior probability of 1.00.

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.