Comparative mitogenomic analysis provides evolutionary insights into Formica (Hymenoptera: Formicidae).

Formica is a large genus in the family Formicidae with high diversity in its distribution, morphology, and physiology. To better understand evolutionary characteristics of Formica, the complete mitochondrial genomes (mitogenomes) of two Formica species were determined and a comparative mitogenomic analysis for this genus was performed. The two newly sequenced Formica mitogenomes each included 37 typical mitochondrial genes and a large non-coding region (putative control region), as observed in other Formica mitogenomes. Base composition, gene order, codon usage, and tRNA secondary structure were well conserved among Formica species, whereas diversity in sequence size and structural characteristics was observed in control regions. We also observed several conserved motifs in the intergenic spacer regions. These conserved genomic features may be related to mitochondrial function and their highly conserved physiological constraints, while the diversity of the control regions may be associated with adaptive evolution among heterogenous habitats. A negative AT-skew value on the majority chain was presented in each of Formica mitogenomes, indicating a reversal of strand asymmetry in base composition. Strong codon usage bias was observed in Formica mitogenomes, which was predominantly determined by nucleotide composition. All 13 mitochondrial protein-coding genes of Formica species exhibited molecular signatures of purifying selection, as indicated by the ratio of non-synonymous substitutions to synonymous substitutions being less than 1 for each protein-coding gene. Phylogenetic analyses based on mitogenomic data obtained fairly consistent phylogenetic relationships, except for two Formica species that had unstable phylogenetic positions, indicating mitogenomic data are useful for constructing phylogenies of ants. Beyond characterizing two additional Formica mitogenomes, this study also provided some key evolutionary insights into Formica.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation.

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Differential gene expression patterns between the head and thorax of Gynaephora aureata are associated with high-altitude adaptation.

Grassland caterpillars (Lepidoptera: Erebidae: Gynaephora) are important pests in alpine meadows of the Qinghai-Tibetan Plateau (QTP). These pests have morphological, behavioral, and genetic adaptations for survival in high-altitude environments. However, mechanisms underlying high-altitude adaptation in QTP Gynaephora species remain largely unknown. Here, we performed a comparative analysis of the head and thorax transcriptomes of G. aureata to explore the genetic basis of high-altitude adaptation. We detected 8,736 significantly differentially expressed genes (sDEGs) between the head and thorax, including genes related to carbohydrate metabolism, lipid metabolism, epidermal proteins, and detoxification. These sDEGs were significantly enriched in 312 Gene Ontology terms and 16 KEGG pathways. We identified 73 pigment-associated genes, including 8 rhodopsin-associated genes, 19 ommochrome-associated genes, 1 pteridine-associated gene, 37 melanin-associated genes, and 12 heme-associated genes. These pigment-associated genes were related to the formation of the red head and black thorax of G. aureata. A key gene, yellow-h, in the melanin pathway was significantly upregulated in the thorax, suggesting that it is related to the formation of the black body and contributed to the adaptation of G. aureata to low temperatures and high ultraviolet radiation in the QTP. Another key gene, cardinal, in the ommochrome pathway was significantly upregulated in the head and may be related to red warning color formation. We also identified 107 olfactory-related genes in G. aureata, including genes encoding 29 odorant-binding proteins, 16 chemosensory proteins, 22 odorant receptor proteins, 14 ionotropic receptors, 12 gustatory receptors, 12 odorant degrading enzymes, and 2 sensory neuron membrane proteins. Diversification of olfactory-related genes may be associated with the feeding habits of G. aureata, including larvae dispersal and searching for plant resources available in the QTP. These results provide new insights into high-altitude adaptation of Gynaephora in the QTP and may contribute to the development of new control strategies for these pests.

Comparative Mitogenomic Analyses of Darkling Beetles (Coleoptera: Tenebrionidae) Provide Evolutionary Insights into tRNA-like Sequences.

Tenebrionidae is widely recognized owing to its species diversity and economic importance. Here, we determined the mitochondrial genomes (mitogenomes) of three Tenebrionidae species (Melanesthes exilidentata, Anatolica potanini, and Myladina unguiculina) and performed a comparative mitogenomic analysis to characterize the evolutionary characteristics of the family. The tenebrionid mitogenomes were highly conserved with respect to genome size, gene arrangement, base composition, and codon usage. All protein-coding genes evolved under purifying selection. The largest non-coding region (i.e., control region) showed several unusual features, including several conserved repetitive fragments (e.g., A+T-rich regions, G+C-rich regions, Poly-T tracts, TATA repeat units, and longer repetitive fragments) and tRNA-like structures. These tRNA-like structures can bind to the appropriate anticodon to form a cloverleaf structure, although base-pairing is not complete. We summarized the quantity, types, and conservation of tRNA-like sequences and performed functional and evolutionary analyses of tRNA-like sequences with various anticodons. Phylogenetic analyses based on three mitogenomic datasets and two tree inference methods largely supported the monophyly of each of the three subfamilies (Stenochiinae, Pimeliinae, and Lagriinae), whereas both Tenebrioninae and Diaperinae were consistently recovered as polyphyletic. We obtained a tenebrionid mitogenomic phylogeny: (Lagriinae, (Pimeliinae, ((Tenebrioninae + Diaperinae), Stenochiinae))). Our results provide insights into the evolution and function of tRNA-like sequences in tenebrionid mitogenomes and contribute to our general understanding of the evolution of Tenebrionidae.

Comparative mitogenomic analyses provide evolutionary insights into the retrolateral tibial apophysis clade (Araneae: Entelegynae).

The retrolateral tibial apophysis (RTA) clade is the largest spider lineage within Araneae. To better understand the diversity and evolution, we newly determined mitogenomes of ten RTA species from six families and performed a comparative mitogenomics analysis by combining them with 40 sequenced RTA mitogenomes available on GenBank. The ten mitogenomes encoded 37 typical mitochondrial genes and included a large non-coding region (putative control region). Nucleotide composition and codon usage were well conserved within the RTA clade, whereas diversity in sequence length and structural features was observed in control region. A reversal of strand asymmetry in nucleotide composition, i.e., negative AT-skews and positive GC-skews, was observed in each RTA species, likely resulting from mitochondrial gene rearrangements. All protein-coding genes were evolving under purifying selection, except for atp8 whose Ka/Ks was larger than 1, possibly due to positive selection or selection relaxation. Both mutation pressure and natural selection might contribute to codon usage bias of 13 protein-coding genes in the RTA lineage. Phylogenetic analyses based on mitogenomic data recovered a family-level phylogeny within the RTA; {[(Oval calamistrum clade, Dionycha), Marronoid clade], Sparassidae}. This study characterized RTA mitogenomes and provided some new insights into the phylogeny and evolution of the RTA clade.