Gene characteristics of the complete mitochondrial genomes of Paratoxodera polyacantha and Toxodera hauseri (Mantodea: Toxoderidae)

Comparative Analysis of the Complete Mitochondrial Genomes of Three Sisoridae (Osteichthyes, Siluriformes) and the Phylogenetic Relationships of Sisoridae

The family Sisoridae is one of the largest and most diverse Asiatic catfish families, with most species occurring in the water systems of the Qinhai-Tibetan Plateau and East Himalayas. At present, the phylogenetic relationship of the Sisoridae is relatively chaotic. In this study, the mitochondrial genomes (mitogenomes) of three species Creteuchiloglanis kamengensis, Glaridoglanis andersonii, and Exostoma sp. were systematically investigated, the phylogenetic relationships of the family were reconstructed and to determine the phylogenetic position of Exostoma sp. within Sisoridae. The lengths of the mitogenomes' sequences of C. kamengensis, G. andersonii, and Exostoma sp. were 16,589 bp, 16,531 bp, and 16,529 bp, respectively. They all contained one identical control region (D-loop), two ribosomal RNAs (rRNAs), 13 protein-coding genes (PCGs) and 22 transfer RNA (tRNA) genes. We applied two approaches, Bayesian Inference (BI) and Maximum Likelihood (ML), to construct phylogenetic trees. Our findings revealed that the topological structure of both ML and BI trees exhibited significant congruence. Specifically, the phylogenetic tree strongly supports the monophyly of Sisorinae and Glyptosternoids and provides new molecular biological data to support the reconstruction of phylogenetic relationships with Sisoridae. This study is of great scientific value for phylogenetic and genetic variation studies of the Sisoridae.

Codon usage bias analysis in the mitochondrial genomes of five Rhingia Scopoli (Diptera, Syrphidae, Eristalinae) species

This study presents the sequencing and annotation of mitochondrial genomes from five Rhingia species of the family Syrphidae, focusing on codon bias. Each species possessed 22 tRNAs genes, 13 protein-coding genes, 2 rRNAs genes, and a control region, without any observed gene rearrangements. Nucleotide composition analysis revealed a higher AT content compared with GC content, indicating AT enrichment. Neutrality plot, Parity rule 2 bias, and effective number of codons plot analyses collectively indicated that natural selection primarily influences the codon usage bias in the five Rhingia species. Relative synonymous codon usage analysis identified the optimal codons for Rhingia binotata, R. fromosana, R. campestris, R. louguanensis, and R. xanthopoda as 10, 14, 10, 11, and 12, respectively, all ending with A/U and exhibiting AT preference. Phylogenetic analysis, based on maximum likelihood and Bayesian inference methods applied to three datasets, confirmed the monophyly of Rhingia. In conclusion, this research establishes a foundation for understanding the phylogenetic evolution and codon usage patterns in Rhingia, offering valuable for future studies.

Molecular insights into codon usage analysis of mitochondrial fission and fusion gene: relevance to neurodegenerative diseases

Mitochondrial dysfunction is the leading cause of neurodegenerative disorders like Alzheimer's disease and Parkinson's disease. Mitochondria is a highly dynamic organelle continuously undergoing the process of fission and fusion for even distribution of components and maintaining proper shape, number, and bioenergetic functionality. A set of genes governs the process of fission and fusion. OPA1, Mfn1, and Mfn2 govern fusion, while Drp1, Fis1, MIEF1, and MIEF2 genes control fission. Determination of specific molecular patterns of transcripts of these genes revealed the impact of compositional constraints on selecting optimal codons. AGA and CCA codons were over-represented, and CCC, GTC, TTC, GGG, ACG were under-represented in the fusion gene set. In contrast, CTG was over-represented, and GCG, CCG, and TCG were under-represented in the fission gene set. Hydropathicity analysis revealed non-polar protein products of both fission and fusion gene set transcripts. AGA codon repeats are an integral part of translational regulation machinery and present a distinct pattern of over-representation and under-representation in different transcripts within the gene sets, suggestive of selective translational force precisely controlling the occurrence of the codon. Out of six synonymous codons, five synonymous codons encoding for leucine were used differently in both gene sets. Hence, forces regulating the occurrence of AGA and five synonymous leucine-encoding codons suggest translational selection. A correlation of mutational bias with gene expression and codon bias and GRAVY and AROMA signifies the selection pressure in both gene sets, while the correlation of compositional bias with gene expression, codon bias, protein properties, and minimum free energy signifies the presence of compositional constraints. More than 25% of codons of both gene sets showed a significant difference in codon usage. The overall analysis shed light on molecular features of gene sets involved in fission and fusion.

Comparative analysis of codon usage patterns and phylogenetic implications of five mitochondrial genomes of the genus Japanagallia Ishihara, 1955 (Hemiptera, Cicadellidae, Megophthalminae)

Japanagallia is a genus of Cicadomorpha in the family of leafhoppers that are plant piercing-sucking insects, and it is difficult to distinguish by morphological characteristics. So far, only one complete mitochondrial genome data has been reported for the genus Japanagallia. Therefore, in order to better understand this group, we assembled and annotated the complete mitochondrial genomes of five Japanagallia species, and analyzed their codon usage patterns. Nucleotide composition analysis showed that AT content was higher than GC content, and the protein-coding sequences preferred to end with A/T at the third codon position. Relative synonymous codon usage analysis revealed most over-represented codon ends with A or T. Parity plot analysis revealed the codon usage bias of mitochondrial genes was influenced by both natural selection and mutation pressure. In the neutrality plot, the slopes of regression lines were < 0.5, suggesting that natural selection was playing a major role while mutation pressure was of minor importance. The effective number of codons showed that the codon usage bias between genes and genomes was low. Correspondence analysis revealed that the codon usage pattern differed among 13 protein-coding genes. Phylogenetic analyses based on three datasets using two methods (maximum likelihood and Bayesian inference), restored the Megophthalminae monophyly with high support values (bootstrap support values (BS) = 100, Bayesian posterior probability (PP) = 1). In the obtained topology, the seven Japanagallia species were clustered into a monophyletic group and formed a sister group with Durgade. In conclusion, our study can provide a reference for the future research on organism evolution, identification and phylogeny relationships of Japanagallia species.

Novel Gene Rearrangements in Mitochondrial Genomes of four families of Praying Mantises (Insecta, Mantodea) and Phylogenetic Relationships of Mantodea

The mitochondrial genome (mitogenome) plays an important role in phylogenetic studies of many species. The mitogenomes of many praying mantis groups have been well-studied, but mitogenomes of special mimic praying mantises, especially Acanthopoidea and Galinthiadoidea species, are still sorely lacking in the NCBI database. The present study analyzes five mitogenomes from four species of Acanthopoidea (Angela sp., Callibia diana, Coptopteryx sp., Raptrix fusca) and one of Galinthiadoidea (Galinthias amoena) that were sequenced by the primer-walking method. Among Angela sp. and Coptopteryx sp., three gene rearrangements were detected in ND3-A-R-N-S-E-F and COX1-L2-COX2 gene regions, two of which were novel. In addition, individual tandem repeats were found in control regions of four mitogenomes (Angela sp., C. diana, Coptopteryx sp., G. amoena). For those, plausible explanations were derived from the tandem duplication-random loss (TDRL) model and the slipped-strand mispairing model. One potential motif was found in Acanthopidae that was seen as a synapomorphy. Several conserved block sequences (CBSs) were detected within Acanthopoidea that paved the way for the design of specific primers. Via BI and ML analysis, based on four datasets (PCG12, PCG12R, PCG123, PCG123R), the merged phylogenetic tree within Mantodea was reconstructed. This showed that the monophyly of Acanthopoidea was supported and that the PCG12R dataset was the most suitable for reconstructing the phylogenetic tree within Mantodea.

Exploring the Mitogenomes of Mantodea: New Insights from Structural Diversity and Higher-Level Phylogenomic Analyses

The recently reorganized classification of Mantodea has made significant progress in resolving past homoplasy problems, although some relationships among higher taxa remain uncertain. In the present study, we utilized newly sequenced mitogenomes and nuclear gene sequences of 23 mantid species, along with published data of 53 mantises, to perform familial-sampling structural comparisons of mantodean mitogenomes and phylogenomic studies. Our rstructural analysis revealed generally conserved mitogenome organizations, with a few cases of tRNA gene rearrangements, including the detection of trnL₂ duplication for the first time. In our phylogenetic analysis, we found a high degree of compositional heterogeneity and lineage-specific evolutionary rates among mantodean mitogenomes, which frequently corresponded to several unexpected groupings in the topologies under site-homogeneous models. In contrast, the topologies obtained using the site-heterogeneous mixture model fit the currently accepted phylogeny of Mantodea better. Topology tests and four-cluster likelihood mapping analyses further determined the preferred topologies. Our phylogenetic results confirm the monophyly of superfamilial groups Schizomantodea, Amerimantodea, Heteromantodea, Promantidea, and Mantidea and recover the early-branching relationships as (Mantoidoidea + (Amerimantodea + (Metallyticoidea + Cernomantodea))). Additionally, the results suggest that the long-unresolved phylogenetic position of Majangidae should be placed within Mantidea, close to Mantoidea, rather than within Epaphroditoidea. Our findings contribute to understanding the compositional and structural diversity in mantodean mitogenomes, underscore the importance of evolutionary model selection in phylogenomic studies, and provide new insights into the high-level phylogeny of Mantodea.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

Evolution of Gene Arrangements in the Mitogenomes of Ensifera and Characterization of the Complete Mitogenome of Schizodactylus jimo

Gene arrangement (relative location of genes) is another evolutionary marker of the mitogenome that can provide extensive information on the evolutionary mechanism. To explore the evolution of gene arrangements in the mitogenome of diversified Ensifera, we sequenced the mitogenome of the unique dune cricket species found in China and used it for phylogenetic analysis, in combination with 84 known Ensiferan mitogenomes. The mitogenome of Schizodactylus jimo is a 16,428-bp circular molecule that contains 37 genes. We identified eight types of gene arrangement in the 85 ensiferan mitogenomes. The gene location changes (i.e., gene translocation and duplication) were in three gene blocks: I-Q-M-ND2, rrnl-rns-V, and ND3-A-R-N-S-E-F. From the phylogenetic tree, we found that Schizodactylus jimo and most other species share a typical and ancient gene arrangement type (Type I), while Grylloidea has two types (Types II and III), and the other five types are rare and scattered in the phylogenetic tree. We deduced that the tandem replication–random loss model is the evolutionary mechanism of gene arrangements in Ensifera. Selection pressure analysis revealed that purifying selection dominated the evolution of the ensiferan mitochondrial genome. This study suggests that most gene rearrangements in the ensiferan mitogenome are rare accidental events.

Mitochondrial genomes of 10 Mantidae species and their phylogenetic implications

This article aims to present a phylogenetic evaluation of Mantidae based on a mitochondrial genome (mitogenome) data set. The mitogenome of 10 Mantidae species were sequenced using next-generation sequencing. The length of nine the complete mitogenomes ranged from 15,371 bp in Tenodera aridifolia to 16,063 bp in Hierodula longa. Mantidae mitogenomes have 37 genes and control region with two exceptions: five trnR copies in Statilia maculata, and H. zhangi was incomplete missing trnI, trnQ, trnM and a portion of the control region. There was a large noncoding region (LNC) between trnM and nad2 in H. chinensis, H. longa, H. maculata and Titanodula sp. Most of protein-coding genes (PCGs) used the typical start ATN codon and TAA/TAG stop codons. All tRNAs fold into the typical clover-leaf secondary structure except trnS1 which lacks a dihydrouracil (DHU) arm. Nucleotide diversity and Ka/Ks analysis of 13 PCGs showed that atp8 had the highest variability and fastest evolutionary rate. Phylogenetic relationships among 42 Mantidae species were reconstructed using the 13 PCGs and two rRNA genes using Bayesian Inference (BI) and Maximum Likelihood (ML) methods. Of the seven mantid subfamilies included in this analysis, only four had multiple exemplars, and of those only Mantinae and Vatinae formed monophyletic groups in BI and ML trees. Consistent with previous studies, the monophyly of the Hierudulinae and Tenoderinae were not been supported. The present results imply that it is necessary to combine nuclear molecular markers and external characteristic to understand the phylogenetic relationships within Mantidae.

Mitochondrial Phylogenomics Suggests Complex Evolutionary Pattern of Pronotal Foliaceous Mimicry in Hierodulinae (Mantodea: Mantidae), with Description of a New Species of Rhombodera Burmeister, 1838 from China

Hierodulinae is a species-rich mantid subfamily, with some species bearing a notable leaf-like pronotum. However, the evolutionary pattern and taxonomic significance of the leaf-like pronotum are largely unknown. Here, we present a phylogenomic analysis of the Hierodulinae genera Rhombodera Burmeister, 1838, and Hierodula Burmeister, 1838 based on mitochondrial genomes. We also describe a new species, namely Rhombodera hyalina sp. nov. from Guangxi, China. Our phylogenetic result, together with the evidence from male genitalia, suggests the division of the Oriental Hierodula and Rhombodera complex into three clades. We find a complex pattern on the evolution of the leaf-like pronotum, which is present in at least five lineages, respectively, of the above three clades.

Novel Gene Rearrangement in the Mitochondrial Genome of Three Garra and Insights Into the Phylogenetic Relationships of Labeoninae

Complete mitochondrial genomes (mitogenomes) can provide valuable information for phylogenetic relationships, gene rearrangement, and molecular evolution. Here, we report the mitochondrial whole genomes of three Garra species and explore the mechanisms of rearrangements that occur in their mitochondrial genomes. The lengths of the mitogenomes’ sequences of Garra dengba, Garra tibetana, and Garra yajiangensis were 16,876, 16,861, and 16,835, respectively. They contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genomes of three Garra species were rearranged compared to other fish mitochondrial genomes. The tRNA-Thr, tRNA-Pro and CR (T-P-CR) genes undergo replication followed by random loss of the tRNA-Thr and tRNA-Pro genes to form tRNA-Thr, CR1, tRNA-Pro and CR2 (T-CR-P-CR). Tandem duplication and random loss best explain this mitochondrial gene rearrangement. These results provide a foundation for future characterization of the mitochondrial gene arrangement of Labeoninae and further phylogenetic studies.

Phylogenetic Implication of Large Intergenic Spacers: Insights from a Mitogenomic Comparison of Prosopocoilus Stag Beetles (Coleoptera: Lucanidae)

Simple Summary

Insect mitochondrial genomes (mitogenomes) show high diversity in some lineages. In the mitogenome of some Coleoptera species, a large intergenic spacer (IGS) has been identified. However, very little is known about mitogenomes of lucanid beetles. In this work, to provide further insight into the phylogenic relationships among species in lucanid beetles (genus Prosopocoilus), two Prosopocoilus species (Prosopocoilus castaneus and Prosopocoilus laterotarsus) were newly sequenced and comparatively analyzed. Significantly, the two newly sequenced Prosopocoilus species contained a large IGS located between trnI and trnQ. Our phylogenomic analyses showed that P. castaneus and P. laterotarsus were clustered in a clade with typical Prosopocoilus species (Prosopocoilus confucius, Prosopocoilus blanchardi, and Prosopocoilusastacoides). These results provide valuable data for the future study of the phylogenetic relationships in this genus.

Abstract

To explore the characteristics of mitogenomes and discuss the phylogenetic relationships within the genus Prosopocoilus, the mitogenomes of two species (P. castaneus and P. laterotarsus) were newly sequenced and comparatively analyzed. The arrangement of the mitogenome in these two lucanid beetles was the same as that in the inferred ancestral insect, and the nucleotide composition was highly biased towards A + T as in other lucanids. The evolutionary rates of 13 protein-coding genes (PCGs) suggested that their evolution was based on purifying selection. Notably, we found evidence of the presence of a large IGS between trnI and trnQ genes, whose length varied from 375 bp (in P. castaneus) to 158 bp (in P. laterotarsus). Within the large IGS region, a short sequence (TAAAA) was found to be unique among these two species, providing insights into phylogenomic reconstruction. Phylogenetic analyses were performed using the maximum likelihood (IQ-TREE) and Bayesian (PhyloBayes) methods based on 13 protein-coding genes (PCGs) in nucleotides and amino acids (AA) from published mitogenomes (n = 29). The genus Prosopocoilus was found to constitute a distinct clade with high nodal support. Overall, our findings suggested that analysis of the characteristics of the large IGS (presence or absence, size, and location) in mitogenomes of the genus Prosopocoilus may be informative for the phylogenetic and taxonomic analyses and for evaluation of the genus Prosopocoilus, despite the dense sampling materials needed.

Novel Mitochondrial Gene Rearrangement and Intergenic Regions Exist in the Mitochondrial Genomes from Four Newly Established Families of Praying Mantises (Insecta: Mantodea)

Simple Summary

Mantodea is regarded as an excellent material to study the gene rearrangements and large non-coding regions (LNCRs) in mitochondrial genomes. Meanwhile, as a result of the convergent evolution and parallelism, the gene rearrangements and LNCRs are specific to some taxonomic groups within Mantodea, which play an important role in phylogenetic relationship research. Nine mitochondrial genomes (mitogenomes) from four newly established families of praying mantises are obtained and annotated. Eight types of gene rearrangements, including four novel types of gene rearrangements in Mantodea, are detected, which can be explained by the tandem replication-random loss (TDRL) model. Moreover, one conserved motif between trnI-trnQ is detected in Toxoderidae. This study shed light on the formation mechanisms of these gene rearrangements and LNCRs in four newly established families of praying mantises.

Abstract

Long non-coding regions (NCRs) and gene rearrangements are commonly seen in mitochondrial genomes of Mantodea and are primarily focused on three regions: CR-I-Q-M-ND2, COX2-K-D-ATP8, and ND3-A-R-N-S-E-F-ND5. In this study, eight complete and one nearly complete mitochondrial genomes of praying mantises were acquired for the purpose of discussing mitochondrial gene rearrangements and phylogenetic relationships within Mantodea, primarily in the newly established families Haaniidae and Gonypetidae. Except for Heterochaeta sp. JZ-2017, novel mitochondrial gene arrangements were detected in Cheddikulama straminea, Sinomiopteryx graham, Pseudovates chlorophaea, Spilomantis occipitalis. Of note is the fact that one type of novel arrangement was detected for the first time in the Cyt b-S2-ND1 region. This could be reliably explained by the tandem replication-random loss (TDRL) model. The long NCR between trnT and trnP was generally found in Iridopteryginae and was similar to the ND4L or ND6 gene. Combined with gene rearrangements and intergenic regions, the monophyly of Haaniidae was supported, whereas the paraphyly of Gonypetidae was recovered. Furthermore, several synapomorphies unique to some clades were detected that conserved block sequences between trnI and trnQ and gaps between trnT and trnP in Toxoderidae and Iridopteryginae, respectively.

The complete mitochondrial genome of Leptomantella tonkinae (Hebard, 1920) (Mantodea: Leptomantellidae) and its phylogeny

The complete mitochondrial (mt) genome of Leptomantella tonkinae (Hebard, 1920) was 15,527 bp in length and contained 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and one control region. The gene arrangement of mt genome of L. tonkinae was identical to the primitive mantis. The overall AT content of the mt genome was 74%. In ML and BI phylogenetic analyses, the monophyly of Leptomantellidae was robustly supported and the clade of Leptomantellidae is a sister clade to the group of (Gonypetidae+(Leptomantellidae+(Amorphoscelidae+Nanomantidae))).

Cryptic Species Exist in Vietnamella sinensis Hsu, 1936 (Insecta: Ephemeroptera) from Studies of Complete Mitochondrial Genomes

Ephemeroptera (Insecta: Pterygota) are widely distributed all over the world with more than 3500 species. During the last decade, the phylogenetic relationships within Ephemeroptera have been a hot topic of research, especially regarding the phylogenetic relationships among Vietnamellidae. In this study, three mitochondrial genomes from three populations of Vienamella sinensis collected from Tonglu (V. sinensis TL), Chun’an (V. sinensis CN), and Qingyuan (V. sinensis QY) in Zhejiang Province, China were compared to discuss the potential existence of cryptic species. We also established their phylogenetic relationship by combining the mt genomes of 69 Ephemeroptera downloaded from NCBI. The mt genomes of V. sinensis TL, V. sinensis CN, and V. sinensis QY showed the same gene arrangement with lengths of 15,674 bp, 15,674 bp, and 15,610 bp, respectively. Comprehensive analyses of these three mt genomes revealed significant differences in mt genome organization, genetic distance, and divergence time. Our results showed that the specimens collected from Chun’an and Tonglu in Zhejiang Province, China belonged to V. sinensis, and the specimens collected from Qingyuan, Zhejiang Province, China were a cryptic species of V. sinensis. In maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees, the monophyly of the family Vietnamellidae was supported and Vietnamellidae has a close relationship with Ephemerellidae.

Gaining insights into the compositional constraints and molecular phylogeny of five silkworms mitochondrial genome

This study was performed to identify codon usage bias (CUB), genetic similarity and phylogenetic analysis of complete mitochondrial genomes along with separate sequences of 13 protein coding genes per each genome from five types of silkworm including Bombyx mori, Bombyx mandarina, Samia cynthia ricini, Antheraea pernyi and Antheraea assama. Nucleotide composition analysis suggested that AT content was higher than GC content and t-test analysis revealed significance difference (p < 0.01) between AT and GC content. Relative synonymous CUB analysis revealed most over-represented codon ends with A or T. Parity plot analysis revealed both natural selection and mutation pressure influenced CUB of mitochondrial genes while neutrality plot analysis suggested that role of natural selection was higher than mutation pressure. The effective number of codons (ENC) revealed the CUB was low among genes and genomes. In phylogenetic analysis of complete mitochondrial genomes, the B. mori fell in a same cluster with Bombyx mandarina and showed the most similarity (96.7%). In terms of protein coding genes, COX1, COX2 and COX3 showed the most obvious differences. In conclusion, comparative analysis of mitochondrial genomes could be used to identify differences in gene organization, accurate phylogenetic analysis and clustering of different types of silkworms.

First Report on Mitochondrial Gene Rearrangement in Non-Biting Midges, Revealing a Synapomorphy in Stenochironomus Kieffer (Diptera: Chironomidae)

Simple Summary

Gene rearrangement is an additional type of data to support relationships of taxa, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. The concept to use mitochondrial gene rearrangements as phylogenetic markers has been proposed since the mid-1980s, the synapomorphic gene rearrangements have been identified from many lineages. However, mitochondrial gene rearrangement has never been observed in the non-biting midges (Diptera: Chironomidae). Here, seven new mitogenomes of the genus Stenochironomus were sequenced and analyzed. Coupled with published data, phylogenetic analyses were performed within Chironominae. The present study showed that mitogenomes of Stenochironomus are showing a higher A+T bias than other chironomid species. A synapomorphic gene rearrangement that the gene order rearranges from trnI-trnQ-trnM to trnI-trnM-trnQ was identified within Stenochironomus, which is the first instance of mitochondrial gene rearrangement discovered in the Chironomidae. The monophyly of the genus Stenochironomus was strongly supported by mitogenomes. Our study provides new insights into the mitochondrial gene order of Chironomidae, and provides a valuable resource for understanding synapomorphic gene rearrangements.

Abstract

(1) Background: Gene rearrangement of mitochondrial genome, especially those with phylogenetic signals, has long fascinated evolutionary biologists. The synapomorphic gene rearrangements have been identified across multiple orders and at many different taxonomic levels, supporting the monophyletic or systematic relationships of related lineages. However, mitochondrial gene rearrangement has never been observed in the non-biting midges (Diptera: Chironomidae); (2) methods: in this study, the complete mitogenomes of seven Stenochironomus species were sequenced and analyzed for the first time; (3) results: each mitogenome of Stenochironomus contains 37 typical genes and a control region. The whole mitogenomes of Stenochironomus species exhibit a higher A+T bias than other published chironomid species. The gene order rearranges from trnI-trnQ-trnM to trnI-trnM-trnQ in all the seven mitogenomes of Stenochironomus, which might be act as a synapomorphy of the genus, supporting the monophyletic of Stenochironomus species. In addition, another derived gene cluster: trnA-trnG-ND3-trnR exists in Stenochironomus tobaduodecimus. The derived gene orders described above are the first case of mitochondrial gene rearrangement in Chironomidae. Coupled with published data, phylogenetic relationships were reconstructed within Chironominae, and strongly supported the monophyly of Stenochironomus; (4) conclusions: our study provides new insights into the mitochondrial gene order of Chironomidae, and provides a valuable resource for understanding the synapomorphic gene rearrangements.

Strategies and Patterns of Codon Bias in Molluscum Contagiosum Virus

Trends associated with codon usage in molluscum contagiosum virus (MCV) and factors governing the evolution of codon usage have not been investigated so far. In this study, attempts were made to decipher the codon usage trends and discover the major evolutionary forces that influence the patterns of codon usage in MCV with special reference to sub-types 1 and 2, MCV-1 and MCV-2, respectively. Three hypotheses were tested: (1) codon usage patterns of MCV-1 and MCV-2 are identical; (2) SCUB (synonymous codon usage bias) patterns of MCV-1 and MCV-2 slightly deviate from that of human host to avoid affecting the fitness of host; and (3) translational selection predominantly shapes the SCUB of MCV-1 and MCV-2. Various codon usage indices viz. relative codon usage value, effective number of codons and codon adaptation index were calculated to infer the nature of codon usage. Correspondence analysis and correlation analysis were performed to assess the relative contribution of silent base contents and significance of codon usage indices in defining bias in codon usage. Among the tested hypotheses, only the second and third hypotheses were accepted.

Comprehensive Draft Genome Analyses of Three Rockfishes (Scorpaeniformes, Sebastiscus) via Genome Survey Sequencing

Sebastiscus species, marine rockfishes, are of essential economic value. However, the genomic data of this genus is lacking and incomplete. Here, whole genome sequencing of all species of Sebastiscus was conducted to provide fundamental genomic information. The genome sizes were estimated to be 802.49 Mb (S. albofasciatus), 786.79 Mb (S. tertius), and 776.00 Mb (S. marmoratus) by using k-mer analyses. The draft genome sequences were initially assembled, and genome-wide microsatellite motifs were identified. The heterozygosity, repeat ratios, and numbers of microsatellite motifs all suggested possibly that S. tertius is more closely related to S. albofasciatus than S. marmoratus at the genetic level. Moreover, the complete mitochondrial genome sequences were assembled from the whole genome data and the phylogenetic analyses genetically supported the validation of Sebastiscus species. This study provides an important genome resource for further studies of Sebastiscus species.

Study on the Characteristic Codon Usage Pattern in Porcine Epidemic Diarrhea Virus Genomes and Its Host Adaptation Phenotype

Porcine epidemic diarrhea virus (PEDV), which classified in the genus Alphacoronavirus, family Coronaviridae, is one of the most important pathogens that cause heavy economic losses in pig industry. Although intensive mutation and recombination analysis of PEDV strains were provided, systematic genome analysis were needed to elucidate the evolution mechanism and codon usage adaptation profiles of the pathogen. Here, a comprehensive investigation was carried out to reveal the systematic evolutionary processes of synonymous codon usage and host-adapted evolution phenotype of PEDV genome. We found a low codon usage bias (CUB) in PEDV genome and that nucleotide compositions, natural selection, mutation pressure and geographical diversity shapes the codon usage patterns of PEDV, with natural selection dominated the overall codon usage bias in PEDV than the others. By using the relative codon deoptimization index (RCDI) and similarity index (SiD) analysis, we observed that genotype II PEDV strains showed the highest level of adaptation phenotype to Sus scrofa than another divergent clade. To the best of our knowledge, this is the first comprehensive report elaborating the codon usage and host adaptation of PEDV. The findings offer an insight into our understanding of factors involved in PEDV evolution, adaptation and fitness toward their hosts.

Complete mitochondrial genomes of four species of praying mantises (Dictyoptera, Mantidae) with ribosomal second structure, evolutionary and phylogenetic analyses

Praying mantises are distributed all over the world. Though some Mantodea mitogenomes have been reported, an evolutionary genomic and phylogenetic analysis study lacks the latest taxonomic system. In the present study, four new mitogenomes were sequenced and annotated. Deroplatys truncate, D. lobate, Amorphoscelis chinensis and Macromantis sp. belong to Deroplatyidae, Amorphoscelidae and Photinaidae family, respectively. Our results indicated that the ATP8 gene may be lost in D. truncate and D. lobata mt genome, and four tRNA genes have not been found in D. truncate, D. lobata and Macromantis sp. A dN/dS pair analysis was conducted and it was found that all genes have evolved under purifying selection. Furthermore, we tested the phylogenetic relationships between the eight families of the Mantodea, including 35 species of praying Mantis. Based on the complete mitochondrial genome data, it was also suggested as sister to Deroplatyidae + Mantidae, Metallyticus sp., the only representative of Metallyticidae, is sister to the remaining mantises. Our results support the taxonomic system of Schwarz and Roy and are consistent with previous studies.

Insight into the Phylogenetic Relationships among Three Subfamilies within Heptageniidae (Insecta: Ephemeroptera) along with Low-Temperature Selection Pressure Analyses Using Mitogenomes

We determined 15 complete and two nearly complete mitogenomes of Heptageniidae belonging to three subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae) and six genera (Afronurus, Epeorus, Leucrocuta, Maccaffertium, Stenacron, and Stenonema). Species of Rhithrogeninae and Ecdyonurinae had the same gene rearrangement of CR-I-M-Q-M-ND2, whereas a novel gene rearrangement of CR-I-M-Q-NCR-ND2 was found in Heptageniinae. Non-coding regions (NCRs) of 25-47 bp located between trnA and trnR were observed in all mayflies of Heptageniidae, which may be a synapomorphy for Heptageniidae. Both the BI and ML phylogenetic analyses supported the monophyly of Heptageniidae and its subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae). The phylogenetic results combined with gene rearrangements and NCR locations confirmed the relationship of the subfamilies as (Heptageniinae + (Rhithrogeninae + Ecdyonurinae)). To assess the effects of low-temperature stress on Heptageniidae species from Ottawa, Canada, we found 27 positive selection sites in eight protein-coding genes (PCGs) using the branch-site model. The selection pressure analyses suggested that mitochondrial PCGs underwent positive selection to meet the energy requirements under low-temperature stress.

Novel tRNA gene rearrangements in the mitochondrial genomes of praying mantises (Mantodea: Mantidae): Translocation, duplication and pseudogenization

Gene rearrangements have been found in several mitochondrial genomes of Mantodea, located in the gene blocks CR-I-Q-M-ND2, COX1-K-D-ATP8 and ND3-A-R-N-S-E-F-ND5. We have sequenced one mitogenome of Amelidae (Yersinia mexicana) and six mitogenomes of Mantidae to discuss the mitochondrial gene rearrangement and the phylogenetic relationship within Mantidae. These mitogenomes showed rearrangements of tRNA genes except for Asiadodis yunnanensis and Hierodula zhangi. These novel gene rearrangements of Mantidae were primarily concentrated in the region of CR-I-Q-M-ND2, including gene translocation, duplication and pseudogenization. For the occurrences of these rearrangements, the tandem duplication-random loss (TDRL) model and slipped-strand mispairing model were suitable to explain. Large non-coding regions (LNCRs) located in the region of CR-I-Q-M-ND2 were detected in most Mantidae species, whereas some LNCRs had high similarity to the control region (CR). Both BI and ML phylogenetic analyses supported the monophyly of Mantidae and the paraphyly of Mantinae. The phylogenetic results with the gene order and the location of NCRs acted as forceful evidence that specific gene rearrangements and special LNCRs may be synapomorphies for several groups of mantises.

Insight into Codon Utilization Pattern of Tumor Suppressor Gene EPB41L3 from Different Mammalian Species Indicates Dominant Role of Selection Force

Simple Summary

The present study envisaged the codon usage pattern analysis of tumor suppressor gene EPB41L3 for the human, brown rat, domesticated cattle, and Sumatran orangutan. Most amino acids are coded by more than one synonymous codon, but they are used in a biased manner. The codon usage bias results from multiple factors like compositional properties, dinucleotide abundance, neutrality, parity, tRNA pool, etc. Understanding codon bias is central to fields as diverse as molecular evolution, gene expressivity, protein translation, and protein folding. This kind of studies is important to see the effects of various evolutionary forces on codon usage. The present study indicated that the selection force is dominant over other forces shaping codon usage in the envisaged organisms.

Abstract

Uneven codon usage within genes as well as among genomes is a usual phenomenon across organisms. It plays a significant role in the translational efficiency and evolution of a particular gene. EPB41L3 is a tumor suppressor protein-coding gene, and in the present study, the pattern of codon usage was envisaged. The full-length sequences of the EPB41L3 gene for the human, brown rat, domesticated cattle, and Sumatran orangutan available at the NCBI were retrieved and utilized to analyze CUB patterns across the selected mammalian species. Compositional properties, dinucleotide abundance, and parity analysis showed the dominance of A and G whilst RSCU analysis indicated the dominance of G/C-ending codons. The neutrality plot plotted between GC12 and GC3 to determine the variation between the mutation pressure and natural selection indicated the dominance of selection pressure (R = 0.926; p < 0.00001) over the three codon positions across the gene. The result is in concordance with the codon adaptation index analysis and the ENc-GC3 plot analysis, as well as the translational selection index (P2). Overall selection pressure is the dominant pressure acting during the evolution of the EPB41L3 gene.

Mitogenome Analysis of Four Lamiinae Species (Coleoptera: Cerambycidae) and Gene Expression Responses by Monochamus alternatus When Infected with the Parasitic Nematode, Bursaphelenchus mucronatus

We determined the mitochondrial gene sequence of Monochamus alternatus and three other mitogenomes of Lamiinae (Insect: Coleoptera: Cerambycidae) belonging to three genera (Aulaconotus, Apriona and Paraglenea) to enrich the mitochondrial genome database of Lamiinae and further explore the phylogenetic relationships within the subfamily. Phylogenetic trees of the Lamiinae were built using the Bayesian inference (BI) and maximum likelihood (ML) methods and the monophyly of Monochamus, Anoplophora, and Batocera genera was supported. Anoplophora chinensis, An. glabripennis and Aristobia reticulator were closely related, suggesting they may also be potential vectors for the transmission of the pine wood pathogenic nematode (Bursaphelenchus xylophilus) in addition to M. alternatus, a well-known vector of pine wilt disease. There is a special symbiotic relationship between M. alternatus and Bursaphelenchus xylophilus. As the native sympatric sibling species of B. xylophilus, B. mucronatus also has a specific relationship that is often overlooked. The analysis of mitochondrial gene expression aimed to explore the effect of B. mucronatus on the energy metabolism of the respiratory chain of M. alternatus adults. Using RT-qPCR, we determined and analyzed the expression of eight mitochondrial protein-coding genes (COI, COII, COIII, ND1, ND4, ND5, ATP6, and Cty b) between M. alternatus infected by B. mucronatus and M. alternatus without the nematode. Expression of all the eight mitochondrial genes were up-regulated, particularly the ND4 and ND5 gene, which were up-regulated by 4-5-fold (p < 0.01). Since longicorn beetles have immune responses to nematodes, we believe that their relationship should not be viewed as symbiotic, but classed as parasitic.

First complete mitogenomes of Diamesinae, Orthocladiinae, Prodiamesinae, Tanypodinae (Diptera: Chironomidae) and their implication in phylogenetics

BACKGROUND

The mitochondrial genome (mitogenome) has been extensively used for phylogenetic and evolutionary analysis in Diptera, but the study of mitogenome is still scarce in the family Chironomidae.

METHODS

Here, the first complete mitochondrial genomes of four Chironomid species representing Diamesinae, Orthocladiinae, Prodiamesinae and Tanypodinae are presented. Coupled with published mitogenomes of two, a comparative mitochondrial genomic analysis between six subfamilies of Chironomidae was carried out.

RESULTS

Mitogenomes of Chironomidae are conserved in structure, each contains 37 typical genes and a control region, and all genes arrange the same gene order as the ancestral insect mitogenome. Nucleotide composition is highly biased, the control region displayed the highest A + T content. All protein coding genes are under purifying selection, and the ATP8 evolves at the fastest rate. In addition, the phylogenetic analysis covering six subfamilies within Chironomidae was conducted. The monophyly of Chironomidae is strongly supported. However, the topology of six subfamilies based on mitogenomes in this study is inconsistent with previous morphological and molecular studies. This may be due to the high mutation rate of the mitochondrial genetic markers within Chironomidae. Our results indicate that mitogenomes showed poor signals in phylogenetic reconstructions at the subfamily level of Chironomidae.

Comparative analysis of the complete mitochondrial genomes of three rockfishes (Scorpaeniformes, Sebastiscus) and insights into the phylogenetic relationships of Sebastidae

Mitochondrial genome is a powerful molecule marker to provide information for phylogenetic relationships and revealing molecular evolution in ichthyological studies. Sebastiscus species, a marine rockfish, are of essential economic value. However, the taxonomic status and phylogenetic relationships of Sebastidae have been controversial so far. Here, the mitochondrial genomes (mitogenomes) of three species, S. tertius, S. albofasciatus, and S. marmoratus, were systemically investigated. The lengths of the mitogenomes' sequences of S. tertius, S. albofasciatus, and S. marmoratus were 16910, 17056, and 17580 bp, respectively. It contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNA (tRNA) genes, and one identical control region (D-loop) among the three species. The genetic distance and Ka/Ks ratio analyses indicated 13 PCGs were suffering purifying selection and the selection pressures were different from certain deep-sea fishes, which were most likely due to the difference in their living environment. The phylogenetic tree was constructed by Bayesian Inference (BI) and Maximum Likelihood (ML). Most interestingly, the results indicated that Sebastidae and Scorpaenidae were grouped into a separate branch, so the taxonomic status of Sebastidae should be classified into subfamily Sebastinae. Our results may lead to a taxonomic revision of Scorpaenoidei.

Six complete mitochondrial genomes of mayflies from three genera of Ephemerellidae (Insecta: Ephemeroptera) with inversion and translocation of trnI rearrangement and their phylogenetic relationships

As a small order of Pterygota (Insecta), Ephemeroptera has almost 3,500 species around the world. Ephemerellidae is a widely distributed common group of Ephemeroptera. However, the relationship among Ephemerellidae, Vietnamellidae and Teloganellidae is still in dispute. In this study, we sequenced six complete mitogenomes of three genera from Ephemerellidae (Insecta: Ephemeroptera): Ephemerella sp. Yunnan-2018, Serratella zapekinae, Serratella sp. Yunnan-2018, Serratella sp. Liaoning-2019, Torleya grandipennis and T. tumiforceps. These mitogenomes were employed to reveal controversial phylogenetic relationships among the Ephemeroptera, with emphasis on the phylogenetic relationships among Ephemerellidae. The lengths of the six mayfly mitogenomes ranged from 15,134 bp to 15,703 bp. Four mitogenomes of Ephemerella sp. Yunnan-2018, Serratella zapekinae, Serratella sp. Yunnan-2018 and Serratella sp. Liaoning-2019 had 22 tRNAs including an inversion and translocation of trnI. By contrast, the mitogenomes of T. tumiforceps and T. grandipennis had 24 tRNAs due to an extra two copies of inversion and translocation of trnI. Within the family Ephemerellidae, disparate gene rearrangement occurred in the mitogenomes of different genera: one copy of inversion and translocation trnI in the genera Ephemerella and Serratella, and three repeat copies of inversion and translocation of trnI in the genus Torleya. A large non-coding region (≥200 bp) between trnS1 (AGN) and trnE was detected in T. grandipennis and T. tumiforceps. Among the phylogenetic relationship of the Ephemeroptera, the monophyly of almost all families except Siphlonuridae was supported by BI and ML analyses. The phylogenetic results indicated that Ephemerellidae was the sister clade to Vietnamellidae whereas Teloganellidae was not a sister clade of Ephemerellidae and Vietnamellidae.

Analysis of codon usage of severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) and its adaptability in dog

•

Codon analysis reveal natural selection and other undefined factors dominates the overall codon usage bias in SARS-CoV-2 rather than mutational pressure.

•

The host adaptation potential of SARS-CoV-2 is more in human as compared to dog.

Severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) is recognized as one of the life-threatening viruses causing the most destructive pandemic in this century. The genesis of this virus is still unknown. To elucidate its molecular evolution and regulation of gene expression, the knowledge of codon usage is a pre-requisite. In this study, an attempt was made to document the genome-wide codon usage profile and the various factors influencing the codon usage patterns of SARS-CoV-2 in human and dog. The SARS-CoV-2 genome showed relative abundance of A and U nucleotides and relative synonymous codon usage analysis revealed that the preferred synonymous codons mostly end with A/U. The analysis of ENc-GC3s, Neutrality and Parity rule 2 plots indicated that natural selection and other undefined factors dominate the overall codon usage bias in SARS-CoV-2 whereas the impact of mutation pressure is comparatively minor. The codon adaptation index and relative codon deoptimization index of SARS-CoV-2 deciphered that human is more favoured host for adaptation compared to dog. These results enhance our understanding of the factors involved in evolution of the novel human SARS-CoV-2 and its adaptability in dog.

The complete mitochondrial genome of Xanthomantis bimaculata (Mantodea: Iridopterygidae) and its phylogeny

The mitochondrial genome sequence of Xanthomantis bimaculata (Mantodea: Iridopterygidae) from Yunnan, China is a circular molecule with the typical insect mitochondrial gene arrangement, which is 15,941 bp in length and contains 22 tRNAs, two rRNAs, 13 protein-coding genes, and one control region. The overall AT content of the mitogenome is 73.12% (A = 37.58%, T = 35.54%, C = 16.54%, G = 10.34%). In BI and ML phylogenetic analyses, X. bimaculata was a sister clade to Sceptuchus simplex. The monophyly of the families Iridopterygidae, Thespidae and Liturgusidae were supported.

Selection constraints determine preference for A/U-ending codons in Taxus contorta

Unequal utilization of synonymous codons is a well-known phenomenon among living organisms. This phenomenon plays a major role in the enhancement of the accuracy and efficiency of translation. Gymnosperms are rarely paid attention in this aspect. Understanding the degree of and determining the forces influencing codon usage bias (CUB) in Taxus contorta, an endangered Himalayan gymnosperm, will prove useful in interpreting the evolutionary characteristics of this species. Using RNAseq data, 93 790 assembled transcripts were clustered into 32 701 unigenes. Around 13 061 full-length sequences were utilized for the analysis of CUB. Compositional properties showed that GC-content ranged from 28.76% to 65.22%, with an average value of 44.28%, suggesting an AT-rich genome. The mean effective number of codons (ENC) value revealed that CUB is not strong in T. contorta. The preferred codons tended to be A/U ending, whereas the avoided codons tended to be G/C ending. A P2 index of 0.54 and a Mutation Responsive Index (MRI) value of -0.02 in addition to the results revealed by the neutrality, ENC, and parity plots showed that natural selection is a predominating factor governing CUB. Mutational pressure, gene length, hydropathiciy, aromaticity, and nucleotide composition influence CUB weakly.

The complete mitochondrial genome of Annamanum lunulatum (Coleoptera: Lamiinae) and its phylogeny

The complete mitochondrial genome of the Annamanum lunulatum is 15,610 bp in length, which contains 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and the A + T-rich region. The arrangement of genes is identical to all know longhorn beetles mitochondrial genomes. The overall AT content of the mitochondrial genome is 75.3%, whereas the AT content of A + T-rich region is 84.3%. In ML and BI phylogenetic analyses, A. lunulatum is a sister clade to Blepephaeus succinctor, and the monophyly of Lamiinae is supported.

The complete mitochondrial genome of Mantis religiosa (Mantodea: Mantidae) from Canada and its phylogeny

The complete mitochondrial genome of Mantis religiosa (Mantodea: Mantidae) from Canada was successfully sequenced. The mitochondrial genome was a circular molecule of 15,521 bp in length, containing 13 protein-coding genes, two rRNA genes, 23 tRNA genes (including an extra tRNAArg gene), and the control region. The AT content of the whole genome was 76.9% and the length of the control region was 636 bp with 81.9% AT content. The structure of the M. religiosa mitochondrial genome from Canada was almost identical to M. religiosa from China and their genetic distance was just 0.017. In Bayesian inference (BI) and maximum likelihood (ML) analyses, we found that M. religiosa was a sister clade to Statilia sp. and the monophyly of the genera Hierodula and Rhombodera was not supported.

The complete mitochondrial genomes of five longicorn beetles (Coleoptera: Cerambycidae) and phylogenetic relationships within Cerambycidae

Cerambycidae is one of the most diversified groups within Coleoptera and includes nearly 35,000 known species. The relationships at the subfamily level within Cerambycidae have not been convincingly demonstrated and the gene rearrangement of mitochondrial genomes in Cerambycidae remains unclear due to the low numbers of sequenced mitogenomes. In the present study, we determined five complete mitogenomes of Cerambycidae and investigated the phylogenetic relationship among the subfamilies of Cerambycidae based on mitogenomes. The mitogenomic arrangement of all five species was identical to the ancestral Cerambycidae type without gene rearrangement. Remarkably, however, two large intergenic spacers were detected in the mitogenome of Pterolophia sp. ZJY-2019. The origins of these intergenic spacers could be explained by the slipped-strand mispairing and duplication/random loss models. A conserved motif was found between trnS2 and nad1 gene, which was proposed to be a binding site of a transcription termination peptide. Also, tandem repeat units were identified in the A + T-rich region of all five mitogenomes. The monophyly of Lamiinae and Prioninae was strongly supported by both MrBayes and RAxML analyses based on nucleotide datasets, whereas the Cerambycinae and Lepturinae were recovered as non-monophyletic.

Analysis of Nipah Virus Codon Usage and Adaptation to Hosts

A recent outbreak of Nipah virus (NiV) in India has caused 17 deaths in people living in districts of Kerala state. Its zoonotic nature, as well as high rate of human-to-human transmission, has led researchers worldwide to work toward understanding the different aspects of the NiV. We performed a codon usage analysis, based on publicly available nucleotide sequences of NiV and its host adaptation, along with other members of the Henipavirus genus in ten hosts. NiV genome encodes nine open reading frames; and overall, no significant bias in codon usage was observed. Aromaticity of proteins had no impact on codon usage. An analysis of preferred codons used by NiV and the tRNA pool in human cells indicated that NiV prefers codons from a suboptimal anticodon tRNA pool. We observed that codon usage by NiV is mainly constrained by compositional and selection pressures, not by mutational forces. Parameters that define NiV and host relatedness in terms of codon usage were analyzed, with a codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index calculations; which indicated that, of all hosts analyzed, NiV was best adapted to African green monkeys. A comparative analysis based on the relative codon deoptimization index (RCDI) for host adaptation of NiV, Hendra virus (HeV), Cedar virus (CedV), and Hendra like Mojiang virus (MojV) revealed that except for dogs and ferrets, all evaluated hosts were more susceptible to HeV than NiV.