Comparison of base composition and codon usage in insect mitochondrial genomes

The systematic codon usage bias has an important effect on genetic adaption in native species

Random mutations increase genetic variety and natural selection enhances adaption over generations. Codon usage biases (CUB) provide clues about the genome adaptation mechanisms of native species and extremophile species. Significant numbers of gene (CDS) of nine classes of endangered, native species, including extremophiles and mesophiles were utilised to compute CUB. Codon usage patterns differ among the lineages of endangered and extremophiles with native species. Polymorphic usage of nucleotides with codon burial suggests parallelism of native species within relatively confined taxonomic groups. Utilizing the deviation pattern of CUB of endangered and native species, I present a calculation parameter to estimate the extinction risk of endangered species. Species diversity and extinction risk are both positively associated with the propensity of random mutation in CDS (Coding DNA sequence). Codon bias tenet profoundly selected and it governs to adaptive evolution of native species.

Organellar Genomes of Sargassum hemiphyllum var. chinense Provide Insight into the Characteristics of Phaeophyceae

Sargassum hemiphyllum var. chinense, a prevalent seaweed along the Chinese coast, has economic and ecological significance. However, systematic positions within Sargassum and among the three orders of Phaeophyceae, Fucales, Ectocarpales, and Laminariales are in debate. Here, we reported the organellar genomes of S. hemiphyllum var. chinense (34,686-bp mitogenome with 65 genes and 124,323 bp plastome with 173 genes) and the investigation of comparative genomics and systematics of 37 mitogenomes and 22 plastomes of Fucales (including S. hemiphyllum var. chinense), Ectocarpales, and Laminariales in Phaeophyceae. Whole genome collinearity analysis showed gene number, type, and arrangement were consistent in organellar genomes of Sargassum with 360 SNP loci identified as S. hemiphyllum var. chinense and two genes (rps7 and cox2) identified as intrageneric classifications of Sargassum. Comparative genomics of the three orders of Phaeophyceae exhibited the same content and different types (petL was only found in plastomes of the order Fucales and Ectocarpales) and arrangements (most plastomes were rearranged, but trnA and trnD in the mitogenome represented different orders) in genes. We quantified the frequency of RNA-editing (canonical C-to-U) in both organellar genomes; the proportion of edited sites corresponded to 0.02% of the plastome and 0.23% of the mitogenome (in reference to the total genome) of S. hemiphyllum var. chinense. The repetition ratio of Fucales was relatively low, with scattered and tandem repeats (nine tandem repeats of 14-24 bp) dominating, while most protein-coding genes underwent negative selection (Ka/Ks < 1). Collectively, these findings provide valuable insights to guide future species identification and evolutionary status of three important Phaeophyceae order species.

JW, Vieira DBR, Brandão RCF, Medeiros DBDA, Nunes Neto JP. Sequencing and Description of the Mitochondrial Genome of Orthopodomyia fascipes (Diptera: Culicidae)

The genus Orthopodomyia Theobald, 1904 (Diptera: Culicidae) comprises 36 wild mosquito species, with distribution largely restricted to tropical and temperate areas, most of which are not recognized as vectors of epidemiological importance due to the lack of information related to their bionomy and involvement in the cycle transmission of infectious agents. Furthermore, their evolutionary relationships are not completely understood, reflecting the scarcity of genetic information about the genus. Therefore, in this study, we report the first complete description of the mitochondrial genome of a Neotropical species representing the genus, Orthopodomyia fascipes Coquillet, 1906, collected in the Brazilian Amazon region. Using High Throughput Sequencing, we obtained a mitochondrial sequence of 15,598 bp, with an average coverage of 418.5×, comprising 37 functional subunits and a final portion rich in A + T, corresponding to the control region. The phylogenetic analysis, using Maximum Likelihood and Bayesian Inference based on the 13 protein-coding genes, corroborated the monophyly of Culicidae and its two subfamilies, supporting the proximity between the tribes Orthopodomyiini and Mansoniini, partially disagreeing with previous studies based on the use of molecular and morphological markers. The information generated in this study contributes to a better understanding of the taxonomy and evolutionary history of the genus and other groups of Culicidae.

Mitochondrial genome comparison and phylogenetic position of Fannia pusio among the Calyptratae flies

Fannia pusio, the chicken dung fly species, remains unexplored despite its forensic, sanitary, and veterinary importance in the Nearctic and Neotropical regions. In this study, we obtained the complete mitochondrial genome of Fannia pusio for the first time using next-generation sequencing. We compared it with previously published mitogenomes of the genus from the Palearctic region, and its phylogenetic position was studied based on the concatenated protein-coding genes (PCGs) dataset of Calyptratae flies. The circular mitochondrial genome of F. pusio is 16,176 bp in length, with a high A + T content (78.3%), whose gene synteny, codon usage analysis, and amino acid frequency are similar to previously reported Fannia mitogenomes. All PCGs underwent purifying selection except the nad2 gene. Interspecific K2P distances of PCGs of Fannia yielded an average of 12.4% (8.1%-21.1%). The Fannia genus is monophyletic and closely related to Muscidae based on molecular data. Further taxonomic sampling is required to deep into the phylogenetic relationships of the originally proposed species-groups and subgroups within the genus. These results provide a valuable dataset for studying the mitochondrial genome evolution and a resource for the taxonomy and systematics of Fannia.

Comparative genome and phylogenetic analysis revealed the complex mitochondrial genome and phylogenetic position of Conopomorpha sinensis Bradley

Conopomorpha sinensis Bradley is a destructive pest that causes severe economic damage to litchi and longan. Previous C. sinensis research has focused on population life tables, oviposition selectivity, pest population prediction, and control technology. However, there are few studies on its mitogenome and phylogenetic evolution. In this study, we sequenced the whole mitogenome of C. sinensis by the third-generation sequencing, and analyzed the characteristics of its mitogenome by comparative genome. The complete mitogenome of C. sinensis is a typical circular and double-stranded structure. The ENC-plot analyses revealed that natural selection could affect the information of codon bias of the protein-coding genes in the mitogenome of C. sinensis in the evolutionary process. Compared with 12 other Tineoidea species, the trnA-trnF gene cluster of tRNA in the C. sinensis mitogenome appears to have a new arrangement pattern. This new arrangement has not been found in other Tineoidea or other Lepidoptera, which needs further exploration. Meanwhile, a long AT repeated sequence was inserted between trnR and trnA, trnE and trnF, ND1 and trnS in the mitogenome of C. sinensis, and the reason for this sequence remains to be further studied. Furthermore, the results of phylogenetic analysis showed that the litchi fruit borer belonged to Gracillariidae, and Gracillariidae was monophyletic. The results will contribute to an improved understanding of the complex mitogenome and phylogeny of C. sinensis. It also will provide a molecular basis for further research on the genetic diversity and population differentiation of C. sinensis.

Complete mitochondrial genome assembly and comparison of Camellia sinensis var. Assamica cv. Duntsa

Camellia sinensis var. Assamica cv. Duntsa (C.duntsa), a valuable Theaceae from Hunan Province, has been looked at as a precious tea resource by local farmers because of its economic and ecological value. Genomics study on C.duntsa is essential for the domestication and enhancement of tea tree varieties. In the present study, we used a hybrid approach based on Illumina and PacBio data to sequence and assemble the mitochondrial genome of C.duntsa. The mitochondrial genome of C.duntsa was estimated to be 1,081,996 base pairs (bp) and eighty-one genes consisting of one pseudogene, three ribosomal RNA (rRNA) genes, thirty transfer RNA (tRNA) genes, and forty-seven protein-coding genes (PCGs). Tetramer repetitions made up 43.90% of simple sequence repeats (SSRs). The codon usage bias of the Theaceae mitochondrial gene atp9 was altered by mutation, but the codon usage of other genes was shaped by natural selection. Besides, there are eighteen gene-containing homologous regions between the chloroplast and mitochondrial genomes of C. duntsa.Some genomes including atp8, cox1, cox3, nad7, nad9, rpl16, rpl2, rps19, rps4, and sdh4 are absent in the mitochondrial genome of several Theaceae plant. However, C. duntsa maintains these genes integrity and functionality. Another gene, rps16, is either lacking from the mitochondrial genome of C. duntsa or is present as a pseudogene. C. duntsa and C. sinensis (OM809792) are very similar, as shown by a collinear match across four species of Theaceae; the most conservative genes are nad5, atp9, cox2, rps3, trnA-TGC, trnI-GAT, rrn18, trnV-GAC, and ccmFN. Similarly, the genome's phylogenetic trees revealed that C. duntsa was the sister species to C. sinensis. The results confirmed that the C. duntsa and C. sinensis (OM809792) mitochondrial genome underwent gene rearrangement.In general, our results shows that genomic information from organelles can help us understand plant phylogeny and can also be used to make molecular markers and study how genetic traits change over time. Our research will contribute to the population genetics and evolution of tea plant.

Characterization, comparative phylogenetic, and gene transfer analyses of organelle genomes of Rhododendron × pulchrum

Rhododendron × pulchrum, an important horticultural species, is widely distributed in Europe, Asia, and North America. To analyze the phylogenetic and organelle genome information of R. × pulchrum and its related species, the organelle genome of R. × pulchrum was sequenced and assembled. The complete mitochondrial genome showed lineage DNA molecules, which were 816,410 bp long and contained 64 genes, namely 24 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 37 protein-coding genes. The chloroplast genome of R. × pulchrum was reassembled and re-annotated; the results were different from those of previous studies. There were 42 and 46 simple sequence repeats (SSR) identified from the mitochondrial and chloroplast genomes of R. × pulchrum, respectively. Five genes (nad1, nad2, nad4, nad7, and rps3) were potentially useful molecular markers. The R. × pulchrum mitochondrial genome collinear alignment among five species of the Ericaceae showed that the mitochondrial genomes of these related species have a high degree of homology with R. × pulchrum in this gene region, and the most conservative genes were trnC-GCA, trnD-GUC, trnM-CAU, trnN-GUU, trnY-GUA, atp4, nad4, nad2, nad5, ccmC, and rrn26. The phylogenetic trees of mitochondrial genome showed that R. simsii was a sister to R. × pulchrum. The results verified that there was gene rearrangement between R. × pulchrum and R. simsii mitochondrial genomes. The codon usage bias of 10 Ericaceae mitochondrial genes and 7 Rhododendron chloroplast genes were influenced by mutation, while other genes codon usages had undergone selection. The study identified 13 homologous fragments containing gene sequences between the chloroplast and mitochondrial genomes of R. × pulchrum. Overall, our results illustrate the organelle genome information could explain the phylogenetics of plants and could be used to develop molecular markers and genetic evolution. Our study will facilitate the study of population genetics and evolution in Rhododendron and other genera in Ericaceae.

Description of the mitogenome and phylogeny of Aedes spp. (Diptera: Culicidae) from the Amazon region

The genus Aedes (Diptera: Culicidae) includes species of great epidemiological relevance, particularly involved in transmission cycles of leading arboviruses in the Brazilian Amazon region, such as the Zika virus (ZIKV), Dengue virus (DENV), Yellow fever virus (YFV), and Chikungunya virus (CHIKV). We report here the first putatively complete sequencing of the mitochondrial genomes of Brazilian populations of the species Aedes albopictus, Aedes scapularis and Aedes serratus. The sequences obtained showed an average length of 14,947 bp, comprising 37 functional subunits, typical in animal mitochondria (13 PCGs, 22 tRNA, and 2 rRNA). The phylogeny reconstructed by Maximum likelihood method, based on the concatenated sequences of all 13 PCGs produced at least two non-directly related groupings, composed of representatives of the subgenus Ochlerotatus and Stegomyia of the genus Aedes. The data and information produced here may be useful for future taxonomic and evolutionary studies of the genus Aedes, as well as the Culicidae family.

Compositional features and pattern of codon usage for mitochondrial CO genes among reptiles

The phenomenon of non-random occurrence of synonymous nucleotide triplets (codons) in the coding sequences of genes is the codon usage bias (CUB). In this study, we used bioinformatic tool kit to analyze the compositional pattern and CUB of mitogenes namely COI, COII and COIII across different orders of reptiles. Estimation of overall base composition in the protein-coding sequences of COI, COII and COIII genes of the reptilian orders revealed an uneven usage of nucleotides. The overall count of A nucleotide was found to be the highest while the overall count of G nucleotide was the least. The CO genes across the three reptilian orders were prominently AT biased. Comparison of the GC proportion at each codon position displayed that GC1 percentage ranked the highest in all the three CO genes of the reptilian orders. SCUO values indicated weaker CUB, while considerable variation of SCUO values existed in the three CO genes across the studied reptiles. Relative synonymous codon usage (RSCU) values indicated that mostly the A ending codons were preferred. Based on the parameters namely neutrality plot, mutational responsive index and translational selection, we could conclude that natural selection was the major evolutionary force in COI, COII and COIII genes in the studied reptilian orders. However, correspondence analysis, parity plot and correlation studies indicated the existence of mutation pressure as well on the CO genes.

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))))).

Mitochondrial genome sequencing and phylogeny of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis (Diptera: Culicidae)

The genus Haemagogus (Diptera: Culicidae) comprises species of great epidemiological relevance, involved in transmission cycles of the Yellow fever virus and other arboviruses in South America. So far, only Haemagogus janthinomys has complete mitochondrial sequences available. Given the unavailability of information related to aspects of the evolutionary biology and molecular taxonomy of this genus, we report here, the first sequencing of the mitogenomes of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis. The mitogenomes showed an average length of 15,038 bp, average AT content of 79.3%, positive AT-skews, negative GC-skews, and comprised 37 functional subunits (13 PCGs, 22 tRNA, and 02 rRNA). The PCGs showed ATN as start codon, TAA as stop codon, and signs of purifying selection. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogenetic analyzes of Bayesian inference and Maximum Likelihood, based on concatenated sequences from all 13 PCGs, produced identical topologies and strongly supported the monophyletic relationship between the Haemagogus and Conopostegus subgenera, and corroborated with the known taxonomic classification of the evaluated taxa, based on external morphological aspects. The information produced on the mitogenomes of the Haemagogus species evaluated here may be useful in carrying out future taxonomic and evolutionary studies of the genus.

Towards the implementation of a DNA barcode library for the identification of Peruvian species of Anastrepha (Diptera: Tephritidae)

The genus Anastrepha is a diverse lineage of fruit-damaging tephritid flies widespread across the Neotropical Region. Accurate taxonomic identification of these flies is therefore of paramount importance in agricultural contexts. DNA barcoding libraries are molecular-based tools based on a short sequence of the mitochondrial COI gene enabling rapid taxonomic identification of biological species. In this study, we evaluate the utility of this method for species identification of Peruvian species of Anastrepha and assemble a preliminary barcode profile for the group. We obtained 73 individual sequences representing the 15 most common species, 13 of which were either assigned to previously recognized or newly established BINs. Intraspecific genetic divergence between sampled species averaged 1.01% (range 0-3.3%), whereas maximum interspecific values averaged 8.67 (range 8.26-17.12%). DNA barcoding was found to be an effective method to discriminate between many Peruvian species of Anastrepha that were tested, except for most species of the fraterculus species group, which were all assigned to the same BIN as they shared similar and, in some cases, identical barcodes. We complemented this newly produced dataset with 86 published sequences to build a DNA barcoding library of 159 sequences representing 56 Peruvian species of Anastrepha (approx. 58% of species reported from that country). We conclude that DNA barcoding is an effective method to distinguish among Peruvian species of Anastrepha outside the fraterculus group, and that complementary methods (e.g., morphometrics, additional genetic markers) would be desirable to assist sensu stricto species identification for phytosanitary surveillance and management practices of this important group of pestiferous flies.

Characteristic and variability of five complete aphid mitochondrial genomes: Aphis fabae mordvilkoi, Aphis craccivora, Myzus persicae, Therioaphis tenera and Appendiseta robiniae (Hemiptera; Sternorrhyncha; Aphididae)

The complete mitochondrial genomes of aphids Aphis fabae mordvilkoi, A. craccivora, Myzus persicae from Aphidinae as well as Therioaphis tenera and Appendiseta robiniae from Calaphidinae were sequenced and compared with the genomes of other aphid species. A. fabae mordvilkoi, Th. tenera and A. robiniae mitogenomes were sequenced and analyzed for the first time. The annotation of A. craccivora and M. persicae were corrected compared to what was previously published. According to our data there is no translocation of tRNA-Tyr gene in A. craccivora mitogenome and this aphid species has an ancestral type of mitochondrial gene order. A + T content in all 5 mitogenomes was higher than 80%. A + T content in the Th. tenera CR was 59.5% which is untypically low. CRs of all 5 studied mitogenomes had 2 conserved motifs at their ends and extended G + C rich region. A. craccivora, M. persicae and Th. tenera had large tandem repeats inside the CRs. Detailed molecular analysis of all 5 aphid mitochondrial genomes showed the importance of a deep understanding of the molecular organization of all the functional regions of the mitochondrial DNA, which helps to avoid mistakes during genome annotation.

Mitochondrial genomes and polymorphic regions of Gonimbrasia belina and Gynanisa maja (Lepidoptera: Saturniidae), two important edible caterpillars of Southern Africa

Mopane worms are the vernacular designation for the edible caterpillars of the African emperor moths Gonimbrasia belina and Gynanisa maja. Both species, particularly G. belina, are widely harvested in Southern Africa, and their populations are declining. Despite their commercial, nutritional, and cultural importance, their genetic data are currently unavailable. We sequenced two complete mitogenomes from each species using Ion Torrent technology, and identified informative markers in the complete mitogenomes of the two species for use in future studies. Comparing the conspecific mitogenomes allowed the identification of regions with high nucleotide diversity in ATP6, ND1, ND4, ND5, ND6, and CYTB genes. The final panels of markers will allow for the survey of 3117 bp in G. belina, and 3990 bp in Gy. maja. Phylogenetic reconstruction within the family Saturniidae recovered the tribe Bunaeini as monophyletic and basal to Saturniidae, and the tribe Attacini as a monophyletic clade nested within the tribe Saturniini. The G. belina and Gy. maja mitogenomes are the first representatives of African Saturniidae, a taxonomic group with relevance as a food resource on the continent. This study represents the first step towards assessing the genetic diversity, population structure, and phylogeography of African edible caterpillars.

Compositional properties and codon usage pattern of mitochondrial ATP gene in different classes of Arthropoda

Codon usage bias (CUB) is defined as the usage of synonymous codons unequally for an amino acid in a gene transcript. It is influenced by both mutation pressure and natural selection and is a species-specific property. In our current study, we used bioinformatic methods to investigate the coding sequences of mitochondrial adenosine triphosphate gene (MT-ATP) in different classes of arthropoda to know the codon usage pattern of the gene as no work was described earlier. The analysis of compositional properties suggested that the gene is AT rich. The effective number of codons revealed the CUB of both ATP6 and ATP8 gene was moderate. Heat map showed that the codons ending with AT were negatively associated with GC3 while the codons ending with GC were positively associated with GC3 in all the classes of arthropoda. Correspondence study revealed that the pattern of codon usage of ATP6 and ATP8 genes differed across classes. Neutrality plot suggested the codon usage bias of these two genes in phylum arthropoda was influenced by both mutation pressure and natural selection.

Comparative analysis of the complete mitochondrial genomes of five Achilidae species (Hemiptera: Fulgoroidea) and other Fulgoroidea reveals conserved mitochondrial genome organization

In the present study, the complete mitochondrial genomes (mitogenomes) of five Achilidae (Hemiptera: Fulgoroidea), Betatropis formosana, two new species (Magadhaideus luodiana sp. nov and Peltatavertexalis horizontalis sp. nov), Plectoderini sp. and Paracatonidia sp., were sequenced for the first time through next-generation sequencing. The five mitogenomes ranged from 15,214 to 16,216 bp in length, with the typical gene content and arrangement usually observed in Hexapods. The motif "ATGATAA" between atp8 and atp6 was found in all the analyzed species. An overlap "AAGCTTA" between trnW and trnC was observed in the mitogenomes of most Fulgoroidea. The structural and compositional analyses of 26 Fulgoroidea mitogenomes, including the gene rearrangement of five tRNAs (trnW, trnC and trnY; trnT and trnP), the A + T content and AT-skew of the whole mitogenomes, and the nuclear acid and amino acid compositions of the protein-coding genes (PCGs), revealed family-level differences between Delphacidae and other families (Achilidae, Flatidae, Fulgoridae, Issidae and Ricaniidae). Phylogenetic analyses of 13 protein-coding genes from 26 Fulgoroidea species by maximum likelihood and Bayesian Inference were consistent and well supported the basal position of Delphacidae, a close affinity among the families Flatidae, Issidae and Ricaniidae, and a close relationship between Achilidae and Fulgoridae.

Comparative analysis of the mitochondrial genome of the fungus Colletotrichum lindemuthianum, the causal agent of anthracnose in common beans

Fungi of the genus Colletotrichum are economically important and are used as models in plant-pathogen interaction studies. In this study, the complete mitochondrial genomes of two Colletotrichum lindemuthianum isolates were sequenced and compared with the mitochondrial genomes of seven species of Colletotrichum. The mitochondrial genome of C. lindemuthianum is a typical circular molecule 37,446 bp (isolate 89 A₂ 2-3) and 37,440 bp (isolate 83.501) in length. The difference of six nucleotides between the two genomes is the result of a deletion in the ribosomal protein S3 (rps3) gene in the 83.501 isolate. In addition, substitution of adenine for guanine within the rps3 gene in the mitochondrial genome of the 83.501 isolate was observed. Compared to the previously sequenced C. lindemuthianum mitochondrial genome, an exon no annotated in the cytochrome c oxidase I (cox1) gene and a non-conserved open reading frame (ncORF) were observed. The size of the mitochondrial genomes of the seven species of Colletotrichum was highly variable, being attributed mainly to the ncORF, ranging from one to 10 and also from introns ranging from one to 11 and which encode a total of up to nine homing endonucleases. This paper reports for the first time by means of transcriptome that then ncORFs are transcribed in Colletotrichum spp. Phylogeny data revealed that core mitochondrial genes could be used as an alternative in phylogenetic relationship studies in Colletotrichum spp. This work contributes to the genetic and biological knowledge of Colletotrichum spp., which is of great economic and scientific importance.

More evolution underground: Accelerated mitochondrial substitution rate in Australian burrowing freshwater crayfishes (Decapoda: Parastacidae)

To further understand the evolutionary history and mitogenomic features of Australia's highly distinctive freshwater crayfish fauna, we utilized a recently described rapid mitogenome sequencing pipeline to generate 24 new crayfish mitogenomes including a diversity of burrowing crayfish species and the first for Astacopsis gouldi, the world's largest freshwater invertebrate. Whole mitogenome-based phylogeny estimates using both Bayesian and Maximum Likelihood methods substantially strengthen existing hypotheses for systematic relationships among Australian freshwater crayfish with evidence of pervasive diversifying selection and accelerated mitochondrial substitution rate among the members of the clade representing strongly burrowing crayfish that may reflect selection pressures for increased energy requirement for adaptation to terrestrial environment and a burrowing lifestyle. Further, gene rearrangements are prevalent in the burrowing crayfish mitogenomes involving both tRNA and protein coding genes. In addition, duplicated control regions were observed in two closely related Engaeus species, together with evidence for concerted evolution. This study significantly adds to the understanding of Australian freshwater crayfish evolutionary relationships and suggests a link between mitogenome evolution and adaptation to terrestrial environments and a burrowing lifestyle in freshwater crayfish.

Digging deeper: new gene order rearrangements and distinct patterns of codons usage in mitochondrial genomes among shrimps from the Axiidea, Gebiidea and Caridea (Crustacea: Decapoda)

BACKGROUND

Whole mitochondrial DNA is being increasingly utilized for comparative genomic and phylogenetic studies at deep and shallow evolutionary levels for a range of taxonomic groups. Although mitogenome sequences are deposited at an increasing rate into public databases, their taxonomic representation is unequal across major taxonomic groups. In the case of decapod crustaceans, several infraorders, including Axiidea (ghost shrimps, sponge shrimps, and mud lobsters) and Caridea (true shrimps) are still under-represented, limiting comprehensive phylogenetic studies that utilize mitogenomic information.

METHODS

Sequence reads from partial genome scans were generated using the Illumina MiSeq platform and mitogenome sequences were assembled from these low coverage reads. In addition to examining phylogenetic relationships within the three infraorders, Axiidea, Gebiidea, and Caridea, we also investigated the diversity and frequency of codon usage bias and mitogenome gene order rearrangements.

RESULTS

We present new mitogenome sequences for five shrimp species from Australia that includes two ghost shrimps, Callianassa ceramica and Trypaea australiensis, along with three caridean shrimps, Macrobrachium bullatum, Alpheus lobidens, and Caridina cf. nilotica. Strong differences in codon usage were discovered among the three infraorders and significant gene order rearrangements were observed. While the gene order rearrangements are congruent with the inferred phylogenetic relationships and consistent with taxonomic classification, they are unevenly distributed within and among the three infraorders.

DISCUSSION

Our findings suggest potential for mitogenome rearrangements to be useful phylogenetic markers for decapod crustaceans and at the same time raise important questions concerning the drivers of mitogenome evolution in different decapod crustacean lineages.

The rate of adaptive evolution in animal mitochondria

We have investigated whether there is adaptive evolution in mitochondrial DNA, using an extensive data set containing over 500 animal species from a wide range of taxonomic groups. We apply a variety of McDonald-Kreitman style methods to the data. We find that the evolution of mitochondrial DNA is dominated by slightly deleterious mutations, a finding which is supported by a number of previous studies. However, when we control for the presence of deleterious mutations using a new method, we find that mitochondria undergo a significant amount of adaptive evolution, with an estimated 26% (95% confidence intervals: 5.7-45%) of nonsynonymous substitutions fixed by adaptive evolution. We further find some weak evidence that the rate of adaptive evolution is correlated to synonymous diversity. We interpret this as evidence that at least some adaptive evolution is limited by the supply of mutations.

Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera:Culicidae), members of the coronator group

BACKGROUND

The Coronator Group currently encompasses six morphologically similar species (Culex camposi Dyar, Culex coronator Dyar and Knab, Culex covagarciai Forattini, Culex usquatus Dyar, Culex usquatissimus Dyar, and Culex ousqua Dyar). Culex coronator has been incriminated as a potential vector of West Nile Virus (WNV), Saint Louis Encephalitis Virus (SLEV), and Venezuelan Equine Encephalitis Virus (VEEV). The complete mitochondrial genome of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi was sequenced, annotated, and analyzed to provide genetic information about these species.

RESULTS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi varied from 15,573 base pairs in Cx. usquatus to 15,576 in Cx. coronator. They contained 37 genes (13 protein-encoding genes, 2 rRNA genes, and 22 tRNA genes) and the AT-rich control region. Comparative analyses of the 37 genes demonstrated the mitochondrial genomes to be composed of variable and conserved genes. Despite the small size, the ATP8, ATP6 plus NADH5 protein-encoding genes were polymorphic, whereas tRNAs and rRNAs were conserved. The control region contained some poly-T stretch. The Bayesian phylogenetic tree corroborated that both the Coronator Group and the Culex pipens complex are monophyletic taxa.

CONCLUSIONS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx. usquatissimus and Cx. camposi share the same gene composition and arrangement features that match to those reported for most Culicidae species. They are composed of the same 37 genes and the AT-rich control region, which contains poly-T stretches that may be involved in the functional role of the mitochondrial genome. Taken together, results of the dN/dS ratios, the sliding window analyses and the Bayesian phylogenetic analyses suggest that ATP6, ATP8 and NADH5 are promising genes to be employed in phylogenetic studies involving species of the Coronator Group, and probably other species groups of the subgenus Culex. Bayesian topology corroborated the morphological hypothesis of the Coronator Group as monophyletic lineage within the subgenus Culex.

Analysis of codon usage bias of mitochondrial genome in Bombyx mori and its relation to evolution

BACKGROUND

Synonymous codon usage bias (SCUB) is an inevitable phenomenon in organismic taxa, generally referring to differences in the occurrence frequency of codons across different species or within the genome of the same species. SCUB happens in various degrees under pressure from nature selection, mutation bias and other factors in different ways. It also attaches great significance to gene expression and species evolution, however, a systematic investigation towards the codon usage in Bombyx mori (B. mori) has not been reported yet. Moreover, it is still indistinct about the reasons contributing to the bias or the relationship between the bias and the evolution of B. mori.

RESULTS

The comparison of the codon usage pattern between the genomic DNA (gDNA) and the mitochondrial DNA (mtDNA) from B. mori suggests that mtDNA has a higher level of codon bias. Furthermore, the correspondence analysis suggests that natural selection, such as gene length, gene function and translational selection, dominates the codon preference of mtDNA, while the composition constraints for mutation bias only plays a minor role. Additionally, the clustering results of the silkworm superfamily suggest a lack of explicitness in the relationship between the codon usage of mitogenome and species evolution.

CONCLUSIONS

Among the complicated influence factors leading to codon bias, natural selection is found to play a major role in shaping the high bias in the mtDNA of B. mori from our current data. Although the cluster analysis reveals that codon bias correlates little with the species evolution, furthermore, a detailed analysis of codon usage of mitogenome provides better insight into the evolutionary relationships in Lepidoptera. However, more new methods and data are needed to investigate the relationship between the mtDNA bias and evolution.

Influence of certain forces on evolution of synonymous codon usage bias in certain species of three basal orders of aquatic insects

Forces that influence the evolution of synonymous codon usage bias are analyzed in six species of three basal orders of aquatic insects. The rationale behind choosing six species of aquatic insects (three from Ephemeroptera, one from Plecoptera, and two from Odonata) for the present analysis is based on phylogenetic position at the basal clades of the Order Insecta facilitating the understanding of the evolution of codon bias and of factors shaping codon usage patterns in primitive clades of insect lineages and their subtle differences in some of their ecological and environmental requirements in terms of habitat-microhabitat requirements, altitudinal preferences, temperature tolerance ranges, and consequent responses to climate change impacts. The present analysis focuses on open reading frames of the 13 protein-coding genes in the mitochondrial genome of six carefully chosen insect species to get a comprehensive picture of the evolutionary intricacies of codon bias. In all the six species, A and T contents are observed to be significantly higher than G and C, and are used roughly equally. Since transcription hypothesis on codon usage demands A richness and T poorness, it is quite likely that mutation pressure may be the key factor associated with synonymous codon usage (SCU) variations in these species because the mutation hypothesis predicts AT richness and GC poorness in the mitochondrial DNA. Thus, AT-biased mutation pressure seems to be an important factor in framing the SCU variation in all the selected species of aquatic insects, which in turn explains the predominance of A and T ending codons in these species. This study does not find any association between microhabitats and codon usage variations in the mitochondria of selected aquatic insects. However, this study has identified major forces, such as compositional constraints and mutation pressure, which shape patterns of codon usage in mitochondrial genes in the primitive clades of insect lineages.