Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Complete mitochondrial genomes of Dactylogyrus crucifer and Dactylogyrus zandti reveal distinct patterns of codon usage within Dactylogyrus

Monogeneans of the genus Dactylogyrus Diesing, 1850, the largest genus in the family Dactylogyridae, mostly parasitize the gills of cyprinoid hosts; however, only 3 Dactylogyrus' mitochondrial genomes (mitogenomes) are studied so far. The aim of this research is to extend our understanding of the mitogenomes of Dactylogyrus. We sequenced the mitogenomes of D. crucifer and D. zandti isolated from Rutilus rutilus and Abramis brama orientalis in northwest China, and then we compared these mitogenomes with other monogeneans. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of D. crucifer and D. zandti and characterized the mitogenomes to understand the gene structure, gene identity, the secondary structures of the 22 tRNA genes, and relative synonymous codon usage. We used the analytic Bayesian Information and Maximum Likelihood methods to determine their associated phylogenetic trees. The mitogenomes of D. crucifer and D. zandti were 14,403 and 18,584 bp, respectively. Organization and positioning of these genes were in accordance with Dactylogyrus lamellatus and Dactylogyrus tuba. The nucleotide composition of Dactylogyridae was different from other families of Monogenea, and the A+T count of genus Dactylogyrus (54 - 58.4 %) was lower than other genus species of the family Dactylogyridea (63.9 - 78.4 %) in protein-coding genes. Dactylogyrus members displayed a codon usage bias. The relative synonymous codon used by Dactylogyrus was not conserved and was lower than other monogeneans. The codon use patterns of closely-related species isolated from closely-related hosts were identical. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrus isolated from host subfamily Leuciscinae formed a sister-group. Our results contributed significantly to an increased database of mitogenomes, more than 50 %, for Dactylogyrus that may help future studies of mitochondrial genes and codon uses for the analysis of monogenean phylogenetics.

Characterization of the Complete Mitochondrial Genome of the Central Highland Grey-Shanked Douc Langur (Pygathrix cinerea), a Critically Endangered Species Endemic to Vietnam (Mammalia: Primates)

The grey-shanked douc langur (Pygathrix cinerea) is a recently described, critically endangered primate, endemic to Vietnam. In this study, we describe the Central Highland species' complete mitochondrial genome (mitogenome-mtDNA). It is a circular molecule with a length of 16,541 base pairs (bp). The genome consists of 37 genes, consistent with those found in most other vertebrates, including 13 protein coding genes, 22 transfer RNAs, and two ribosomal RNAs. A comparison with the mitogenomes of more than 50 primates showed that the mitogenome of Vietnamese Central Highland Pygathrix cinerea has a conservative gene order. We identified 43 nucleotide differences when comparing this genome with a previously published mitogenome of Pygathrix cinerea. It is evident that there are distinct differences between the Pygathrix cinerea we are currently studying and other Pygathrix cinerea specimens. These differences are unlikely to be solely the result of sequencing errors, as the mitogenomes were generated using high-quality methods. The genetic divergence observed between the two Pygathrix cinerea mitogenomes implies the potential existence of at least two distinct lineages or forms of this primate species within its native range in Vietnam.

Evolutionary Insights from the Mitochondrial Genome of Oikopleura dioica: Sequencing Challenges, RNA Editing, Gene Transfers to the Nucleus, and tRNA Loss

Sequencing the mitochondrial genome of the tunicate Oikopleura dioica is a challenging task due to the presence of long poly-A/T homopolymer stretches, which impair sequencing and assembly. Here, we report on the sequencing and annotation of the majority of the mitochondrial genome of O. dioica by means of combining several DNA and amplicon reads obtained by Illumina and MinIon Oxford Nanopore Technologies with public RNA sequences. We document extensive RNA editing, since all homopolymer stretches present in the mitochondrial DNA correspond to 6U-regions in the mitochondrial RNA. Out of the 13 canonical protein-coding genes, we were able to detect eight, plus an unassigned open reading frame that lacked sequence similarity to canonical mitochondrial protein-coding genes. We show that the nad3 gene has been transferred to the nucleus and acquired a mitochondria-targeting signal. In addition to two very short rRNAs, we could only identify a single tRNA (tRNA-Met), suggesting multiple losses of tRNA genes, supported by a corresponding loss of mitochondrial aminoacyl-tRNA synthetases in the nuclear genome. Based on the eight canonical protein-coding genes identified, we reconstructed maximum likelihood and Bayesian phylogenetic trees and inferred an extreme evolutionary rate of this mitochondrial genome. The phylogenetic position of appendicularians among tunicates, however, could not be accurately determined.

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Comparative mitogenomics of marine angelfishes (F: Pomacanthidae)

The targeted capture of ultraconserved elements (UCEs) has substantially increased the amount of genetic data available for phylogenomic reconstructions. These capture datasets frequently contain mitochondrial DNA as a by-product, often in the form of complete mitogenomes. These can be efficiently harvested to expand existing datasets without additional costs. Here, we present new mitochondrial genomes for six marine angelfish species (F: Pomacanthidae), assembled and annotated from off-target UCE reads. We provide the first comparative analysis of all mitochondrial genomes available for the Pomacanthidae. Results showed that the average length of pomacanthid mitogenomes is 16.8 kbp. Total GC and AT content varied between 44.5% and 46.3%, and 53.7% and 55.5%, respectively. The architecture of angelfish mitogenomes was comparable to that seen in other fish species with 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and the control region. All 13 PCGs evolved under purifying selection, highlighting a high level of selection pressure and gene expression to preserve genetic integrity. The ND6 and ATP8 genes had the highest ratio of non-synonymous (dN) to synonymous (dS) substitutions, indicating a relaxation of purifying selection constraints. Finally, these newly assembled mitogenomes will allow further investigations of the population genetics, systematics and evolutionary biology of one of the most prominent reef fish family in the aquarium trade.

Comparative Analysis of Complete Mitochondrial Genomes of Five Chromodorididae Species (Nudibranchia:Doridina)

Mitochondrial genome is an important molecular marker for exploring the phylogenetic relationships of species and revealing molecular evolution. In the present study, 5 mitogenomes of Chromodorididae (Chromodoris lochi, Chromodoris colemani, Chromodoris elisabethina, Chromodoris annae and Hypselodoris whitei) were systemically investigated. The lengths of the mitogenomes sequences were 14248 bp, 14257 bp, 14252 bp, 14254 bp and 14856 bp, respectively. Most protein-coding genes (PCGs) were initiated with the common ATG codon and terminated with the TAA and TAG. We calculated Ka/Ks values for all 13 PCGs of Chromodorididae species, all ratios were less than 1, indicating selection by purification. Phylogenetic relationships were constructed by Bayesian inference (BI) and maximum likelihood (ML) methods based on all complete genomes of 50 species, primarily from the family Chromodorididae (Doridina) and 2 outgroups. This phylogenetic tree provided further additional references for the classification of the suborder Doridina. Gene rearrangement suggested a more conserved pattern of gene sequences in the superfamily Chromodoridoidea. These results and newly sequenced will contribute to a better understanding of Chromodorididae and provide reference for further phylogenetic studies.

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Effect of locomotor preference on the evolution of mitochondrial genes in Bovidae

Locomotor preferences and habitat types may drive animal evolution. In this study, we speculated that locomotor preference and habitat type may have diverse influences on Bovidae mitochondrial genes. We used selection pressure and statistical analysis to explore the evolution of mitochondrial DNA (mtDNA) protein-coding genes (PCGs) from diverse locomotor preferences and habitat types. Our study demonstrates that locomotor preference (energy demand) drives the evolution of Bovidae in mtDNA PCGs. The habitat types had no significant effect on the rate of evolution in Bovidae mitochondrial genes. Our study provides deep insight into the adaptation of Bovidae.

Complete Mitochondrial Genome of Four Peristediidae Fish Species: Genome Characterization and Phylogenetic Analysis

The systematic revision of the family Peristediidae remains an unresolved issue due to their diverse and unique morphology. Despite the popularity of using mitochondrial genome research to comprehensively understand phylogenetic relationships in fish, genetic data for peristediid fish need to be included. Therefore, this study aims to investigate the mitochondrial genomic characteristics and intra-family phylogenetic relationships of Peristediidae by utilizing mitochondrial genome analysis. Therefore, this study aims to investigate the phylogenetic relationship of Peristediidae by utilizing mitochondrial genome analysis. The mitochondrial genome of four species of Peristediidae (Peristedion liorhynchus, Satyrichthys welchi, Satyrichthys rieffeli, and Scalicus amiscus) collected in the East China Sea was studied. The mitochondrial gene sequence lengths of four fish species were 16,533 bp, 16,526 bp, 16,527 bp, and 16,526 bp, respectively. They had the same mitochondrial structure and were all composed of 37 genes and one control region. Most PCGs used ATG as the start codon, and a few used GTG as the start codon. An incomplete stop codon (TA/T) occurred. The AT-skew and GC-skew values of 13 PCGs from four species were negative, and the GC-skew amplitude was greater than that of AT-skew. All cases of D-arm were found in tRNA-Ser (GCT). The Ka/Ks ratio analysis indicated that 13 PCGs were suffering purifying selection. Based on 12 PCGs (excluding ND6) sequences, a phylogenetic tree was constructed using Bayesian inference (BI) and maximum likelihood (ML) methods, providing a further supplement to the scientific classification of Peristediidae fish. According to the results of divergence time, the four species of fish had apparent divergence in the Early Cenozoic, which indicates that the geological events at that time caused the climax of species divergence and evolution.

Characterization of the complete mitochondrial genome of Desmaulus extinctorium (Littorinimorpha, Calyptraeoidea, Calyptraeidae) and molecular phylogeny of Littorinimorpha

For the purpose of determining the placement of Calyptraeidae within the Littorinimorpha, we hereby furnish a thorough analysis of the mitochondrial genome (mitogenome) sequence of Desmaulus extinctorium. This mitogenome spans 16,605 base pairs and encompasses the entire set of 37 genes, including 13 PCGs, 22 tRNAs and two rRNAs, with an evident AT bias. Notably, tRNASer1 and tRNASer2 lack dihydrouracil (DHU) arms, resulting in an inability to form a secondary structure. Similarly, tRNAAla lacks a TΨC arm, rendering it incapable of forming a secondary structure. In contrast, the remaining tRNAs demonstrate a characteristic secondary structure reminiscent of a cloverleaf. A comparison with ancestral gastropods reveals distinct differences in three gene clusters (or genes), encompassing 15 tRNAs and eight PCGs. Notably, inversions and translocations represent the major types of rearrangements observed in D. extinctorium. Phylogenetic analysis demonstrates robust support for a monophyletic grouping of all Littorinimorpha species, with D. extinctorium representing a distinct Calyptraeoidea clade. In summary, this investigation provides the first complete mitochondrial dataset for a species of the Calyptraeidae, thus providing novel insights into the phylogenetic relationships within the Littorinimorpha.

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

MtDNA deletions and aging

Aging is the major risk factor in most of the leading causes of mortality worldwide, yet its fundamental causes mostly remain unclear. One of the clear hallmarks of aging is mitochondrial dysfunction. Mitochondria are best known for their roles in cellular energy generation, but they are also critical biosynthetic and signaling organelles. They also undergo multiple changes with organismal age, including increased genetic errors in their independent, circular genome. A key group of studies looking at mice with increased mtDNA mutations showed that premature aging phenotypes correlated with increased deletions but not point mutations. This generated an interest in mitochondrial deletions as a potential fundamental cause of aging. However, subsequent studies in different models have yielded diverse results. This review summarizes the research on mitochondrial deletions in various organisms to understand their possible roles in causing aging while identifying the key complications in quantifying deletions across all models.

Comparative mitogenomics of Brachiopods reveals conservatism in articulate species and unusualness in inarticulate species

BACKGROUND

Brachiopods are a phylum of marine invertebrates with over 10,000 fossil species. Today, there are fewer than 500 extant species assigned to the class Articulata or Inarticulata and for which knowledge of evolutionary genetics and genomics is still poor. Until now, complete mitogenome sequences of two inarticulate species and four articulate species were available.

METHODS AND RESULTS

The complete mitogenome of the inarticulate brachiopod species Lingula reevii (20,778 bp) was obtained by using next generation sequencing. It contains 12 protein-coding genes (the annotation of atp8 is unsure), two ribosomal RNA genes, 26 transfer RNA genes, and one supernumerary ORF that is also conserved in the inarticulate species Lingula anatina. It is hypothesized that this ORF could represent a Lingula-specific mtORFan gene (without obvious homology to other genes). Comparative mitogenomics indicate the mitochondrial gene order of L. reevii is unique among brachiopods, and that compared to articulate species, inarticulate species exhibit massive mitogenome rearrangements, deviant ATP8 protein sequences and supernumerary ORFs, possibly representing species- or lineage-specific mtORFan genes.

CONCLUSION

The results of this study enrich genetics knowledge of extant brachiopods, which may eventually help to test hypotheses about their decline.

Comparative mitochondrial genome brings insights to slight variation in gene proportion and large intergenic spacer and phylogenetic relationship of mudskipper species

Fish mitochondrial genome have been largely studied worldwide for evolutionary and other genetic purposes and the structure and gene organization are commonly conservative. However, several studies have demonstrated that this scenario may present variations in some taxa, showing differentiation on the gene rearrangement. In this study, the complete mitogenome of terrestrial fish Boleophthalmus dussumieri was generated and compared with other species of the Exudercidae fishes. The newly complete mitogenome generated is circular and 16,685 bp of length, and it contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA genes (tRNAs), and one control region (CR), with high conservative structure, like other Mudskippers. Most of the PCG showed similar codon usage bias. The gene length was found to be different specially for the CR, 12S rRNA gene and ND5 gene in some taxon. All the Boleophthalmus species showed a gene duplication in the CR, except for B. dussumieri, and they presented a long intergenic spacer specially on the tRNA-Pro/ OH Tandem duplication/random loss (TDRL) and dimer-mitogenome and nonrandom loss (DMNL) are suitable to explain the mitogenome rearrangement observed in this study. The phylogenetic analysis well supported the monophyly of all mudskipper species and the analysis positioned the Periophthalmus clade as the most basal of the terrestrial fishes. This finding provides basis and brings insights for gene variation, gene rearrangements and replications showing evidence for variety of mitochondrial structure diversity within mudskippers.

Variability and evolution of gene order rearrangement in mitochondrial genomes of arthropods (except Hexapoda)

In the species-rich Phylum Arthropoda, the mitochondrial genome is relatively well conserved both in terms of number and order of genes. However, specific clades have a 'typical' gene order that differs from the putative arthropod ancestral arrangement. The aim of this work was to compare the rate of mitochondrial gene rearrangements at inter- and intra-taxonomic levels in the Arthropoda and to postulate the most parsimonious ancestral orders representing the four major arthropod lineages. For this purpose, we performed a comparative genomic analysis of arthropod mitochondrial genomes available in the NCBI database. Using a combination of bioinformatics methods that examined mitochondrial gene rearrangements in 464 species of arthropods from three subphyla (Chelicerata, Myriapoda, and Crustacea [except Hexapoda, previously analyzed]), we observed differences in the rate of rearrangement within major lineages. A higher rate of mitochondrial genome rearrangement was observed in Crustacea and Chelicerata compared to Myriapoda. Likewise, early branching clades exhibit less variability in mitochondrial genome order than late branching clades, within each subphylum. We identified 'hot regions' in the mitochondrial genome of each studied subphylum, and postulated the most likely ancestral gene order in each subphylum and taxonomic order. Our work provides new evidence on the evolutionary dynamics of mitochondrial genome gene order in arthropods and new mitochondrial genome architectures in different taxonomic divisions within each major lineage of arthropods.

Complete mitochondrial genome of Stethoconus japonicus (Hemiptera: Miridae): Insights into the evolutionary traits within the family Miridae

The mitochondrial (mt) genome sequence of insects possesses numerous evolutionary traits. To better understand the evolution of mt genomes within the family Miridae, the complete mt genome of the predatory Japanese plant bug Stethoconus japonicus Schumacher was sequenced before undertaking a comparative analysis of all reported plant bug mt genomes. The mt genome of S. japonicus is a closed-circular and double-stranded DNA molecule of 16,274 bp (GenBank: MK341530), which consists of 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and a putative control region (CR). Consistent with other plant bugs, the mt genome of S. japonicus is strongly AT-biased (73.49 %) with A-skewed (0.202) and C-skewed (-0.248). All 13 PCGs initiate translation using ATN codons and TAA served as complete stop codons for eight PCGs, which as incomplete stop codon "T-" for cox1, nad1, nad5-6 and "TA-" for cox2. Regarding other features, all 22 tRNAs could be folded into typical cloverleaf secondary structures. The control region is 1,717 bp and contains a long tandem repeat sequence of a 165 bp unit repeated six times. Similar sequence with variable number of tandemly repeated units from intra-genus CRs is a distinct characteristic of plant bug mt genomes. Phylogenetic relationships of 15 bugs were eventually analyzed based on Maximum likelihood (ML) and Bayesian inference (BI) methods using 17 mt genome sequences. In the phylogenetic trees, species from a same genus or subfamily are clustered into a branch with high supporting values.and the result suggest that Deraeocorinae is more closely related to Mirinae than Bryocorinae. Finally, this study revealed that mutation of tRNA anticodon is a useful phylogenetic marker that could be of significance for studies of evolutionary patterns.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

The complete mitochondrial genome of Echinolaelaps fukienensis provide insights into phylogeny and rearrangement in the superfamily Dermanyssoidea

BACKGROUND

Echinolaelaps fukienensis is the dominant mite species parasitic on the body surface of the genus Niviventer. The mitochondrial genome (mitogenome) has its own independent genetic material and genetic system, and is now widely used in population genetics, genealogical biogeography, phylogeny and molecular evolution studies. Species diversity of the superfamily Dermanyssoidea is very rich, but its mitogenomes AT content is high, and it is difficult to amplify the complete mitogenome by routine PCR. To date, we have only obtained the mitogenomes of 6 species, scarcity on sequence data has greatly impeded the studies in the superfamily Dermanyssoidea.

METHODS

Echinolaelaps fukienensis were collected in 2019 from the body surface of Niviventer confucianus (Rodentia, Muridae) in Yunnan Province. The E. fukienensis mitogenome was determined and analyzed for the first time using the Illumina Novoseq 6000 platform. Phylogenetic analyses of the superfamily Dermanyssoidea were conducted based on the entire mitogenome sequences.

RESULTS

The E. fukienensis mitogenome was 14,402 bp, which is known the smallest genome of the superfamily Dermanyssoidea, encoding a total of 37 genes, including 13 PCGs, 22 tRNAs, 2 rRNAs and 1 control region. Most protein-coding genes use ATN as the start codon and TAN as the stop codon. AT and GC skew of atp8 genes in E. fukienensis were both 0. The average length of 22 tRNA genes of E. fukienensis was 64 bp, and secondary structures of tRNAs showed base mismatches and missing D-arms in many places. Compared with gene arrangement pattern of the hypothetical ancestor of arthropods, the E. fukienensis mitogenome shows a novel arrangement pattern. Phylogenetic tree supported the monophyly of the superfamily Dermanyssoidea. Echinolaelaps fukienensis being the least genetic distant (0.2762) and most closely related to Varroa destructor.

CONCLUSIONS

This study analyzed comprehensive the structure and evolution of the E. fukienensis mitogenome for the first time, enriches molecular data of the genus Echinolaelaps, which will contribute to further understand phylogeny and rearrangement patterns of the superfamily Dermanyssoidea.

The evolutionary characterization of Gamasida based on mitochondrial genes codon usage pattern

Mites belonging to the suborder Gamasida are species-rich and habitat-diverse, with a worldwide distribution. To adapt to the environment and obtain better living conditions, all species of the suborder Gamasida have been undergoing constant evolution. The complete mitochondrial genome (mitogenome) is an invaluable molecular marker for studying the origin of species, genetic differentiation between closely related species, and between intraspecific groups. In some species of the suborder Gamasida, mitochondrial tRNA genes are truncated and carried unstable genetic information. This study presents a comparative analysis of codon usage pattern and preference of 13 protein-coding genes of 24 species in 17 genera and 10 families of the suborder Gamasida. Results showed that have an obvious AT preference (0.664-0.829) for codon usage in the suborder Gamasida. Most of the optimal and high-frequency codons also end in A/T. The degree of natural selection varies between the same protein-coding genes of different gamasid mites or among different protein-coding genes within the same gamasid mites. Base and codon usage pattern and preference are very similar between the same species and genus, namely the closer species, the more similar their bases and codons usage patterns and preference are. T bases and C bases were the preference bases for codon usage of 24 species in the suborder Gamasida. Evolution of the suborder Gamasida was dominated by natural selection (64.1%). This study provides the first comprehensive analysis of codon usage in the suborder Gamasida, which will greatly improve our understanding of codon usage patterns and preference, genetics, and evolution of the suborder Gamasida. It will help to evaluate the degree of molecular adaptation in the suborder Gamasida and to further explore evolutionary features of the suborder Gamasida.

The complete mitochondrial genome of Morishitium polonicum (Trematoda, Cyclocoelidae) and its phylogenetic implications

Trematodes can adversely impact the health and survival of wild animals. The trematode family Cyclocoelidae, which includes large digenean bird parasites, lacks molecular analysis, and reclassifications have not been supported. This study produced the first fully assembled and annotated mitochondrial genome sequence for the trematode Morishitium polonicum. The whole length of the M. polonicum (GenBank accession number: OP930879) mitogenome is 14083 bp, containing 22 transfer ribonucleic acids (tRNAs), 2 ribosomal RNAs (rRNAs, rrnL and rrnS), and a noncoding control section (D-loop) 13777 to 13854 bp in length. The 12 PCG areas have 3269 codons and a total length of 10053 bp, which makes up 71.38% of the mitochondrial genome's overall sequence. Most (10/12) of the PCGs that code for proteins begin with ATG, while the nad4L and nad1 genes have a GTG start codon. Phylogenetic analysis using the concatenated nucleotide sequences of 12 PCGs, and the ML tree analysis results showed that M. polonicum is more closely related to with Echinostomatidae and Fasciolidae, which indicates that the family Cyclocoelidae is more closely associated with Echinochasmidae. This study provides mtDNA information, and analysis of mitogenomic structure and evolution. Moreover, we aimed to understand the phylogenetic relationships of this fluke.

Phylomitogenomics elucidates the evolution of symbiosis in Thoracotremata (Decapoda: Cryptochiridae, Pinnotheridae, Varunidae)

Background

Thoracotremata belong to the large group of "true" crabs (infraorder Brachyura), and they exhibit a wide range of physiological and morphological adaptations to living in terrestrial, freshwater and marine habitats. Moreover, the clade comprises various symbiotic taxa (Aphanodactylidae, Cryptochiridae, Pinnotheridae, some Varunidae) that are specialised in living with invertebrate hosts, but the evolutionary history of these symbiotic crabs is still partially unresolved.

Methods

Here we assembled and characterised the complete mitochondrial genomes (hereafter mitogenomes) of three gall crab species (Cryptochiridae): Kroppcarcinus siderastreicola, Opecarcinus hypostegus and Troglocarcinus corallicola. A phylogenetic tree of the Thoracotremata was reconstructed using 13 protein-coding genes and two ribosomal RNA genes retrieved from three new gall crab mitogenomes and a further 72 available thoracotreme mitogenomes. Furthermore, we applied a comparative analysis to characterise mitochondrial gene order arrangement, and performed a selection analysis to test for selective pressure of the protein-coding genes in symbiotic Cryptochiridae, Pinnotheridae, and Varunidae (Asthenognathus inaequipes and Tritodynamia horvathi).

Results

The results of the phylogenetic reconstruction confirm the monophyly of Cryptochiridae, which clustered separately from the Pinnotheridae. The latter clustered at the base of the tree with robust branch values. The symbiotic varunids A. inaequipes and T. horvathi clustered together in a clade with free-living Varunidae species, highlighting that symbiosis in the Thoracotremata evolved independently on multiple occasions. Different gene orders were detected in symbionts and free-living species when compared with the ancestral brachyuran gene order. Lastly, the selective pressure analysis detected two positively selected sites in the nad6 gene of Cryptochiridae, but the evidence for positive selection in Pinnotheridae and A. inaequipes and T. horvathi was weak. Adaptive evolution of mitochondrial protein-coding genes is perhaps related to the presumably higher energetic demands of a symbiotic lifestyle.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

Origin matters: Using a local reference genome improves measures in population genomics

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

Insights into the Deep Phylogeny and Novel Gene Rearrangement of Mytiloidea from Complete Mitochondrial Genome

The extant marine mussels which belong to the Mytiloidea are widespread species inhabiting mostly coastal waters, with some distributed in the deep sea. To clarify the classification systems and phylogenetic relationships range from genus to family level within Mytiloidea, new sequence was used in a phylogenetic analysis including all the available Mytiloidea mitochondrial genomes. In this study, the complete mitochondrial genome of Vignadula atrata is 15,624 bp in length and contains 12 protein-coding genes (PCGs, atp8 is absent), two ribosomal RNA genes, and 22 transfer RNA genes. Phylogenetic analysis based on 12 PCGs showed that it has a close relationship to Bathymodiolus. The analysis of gene rearrangements in the Pteriomorphia showed that the arrangements are highly variable across species, novel gene rearrangements were found within Mytiloidea. The V. atrata mitogenome was provided in detail, with notes on the sequence and a key to the species of Vignadula. This study provides a perspective on the taxonomic histories of the marine mussels and refines the unclear relationship between the origin and evolution of species in Mytiloidea within Bivalvia.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

A computational study on mitogenome-encoded proteins of Pavo cristatus and Pavo muticus identifies key genetic variations with functional implications

BACKGROUND

The Pavo cristatus population, native to the Indian subcontinent, is thriving well in India. However, the Pavo muticus population, native to the tropical forests of Southeast Asia, has reduced drastically and has been categorised as an endangered group. To understand the probable genetic factors associated with the decline of P. muticus, we compared the mitogenome-encoded proteins (13 proteins) between these two species.

RESULTS

Our data revealed that the most frequent variant between these two species was mtND1, which had an alteration in 9.57% residues, followed by mtND5 and mtATP6. We extended our study on the rest of the proteins and observed that cytochrome c oxidase subunits 1, 2, and 3 do not have any change. The 3-dimensional structure of all 13 proteins was modeled using the Phyre2 programme. Our data show that most of the proteins are alpha helical, and the variations observed in P. muticus reside on the surface of the respective proteins. The effect of variation on protein function was also predicted, and our results show that amino acid substitution in mtND1 at 14 sites could be deleterious. Similarly, destabilising changes were observed in mtND1, 2, 3, 4, 5, and 6 and mtATP6-8 due to amino acid substitution in P. muticus. Furthermore, protein disorder scores were considerably altered in mtND1, 2, and 5 of P. muticus.

CONCLUSIONS

The results presented here strongly suggest that variations in mitogenome-encoded proteins of P. cristatus and P. muticus may alter their structure and functions. Subsequently, these variations could alter energy production and may correlate with the decline in the population of P. muticus.

Evolutionary Rates, Divergence Rates, and Performance of Individual Mitochondrial Genes Based on Phylogenetic Analysis of Copepoda

Copepoda is a large and diverse group of crustaceans, which is widely distributed worldwide. It encompasses roughly 9 orders, whose phylogeny remains unresolved. We sequenced the complete mitochondrial genome (mitogenome) of Sinergasilus major (Markevich, 1940) and used it to explore the phylogeny and mitogenomic evolution of Copepoda. The mitogenome of S. major (14,588 bp) encodes the standard 37 genes as well as a putative control region, and molecular features are highly conserved compared to other Copepoda mitogenomes. Comparative analyses indicated that the nad2 gene has relatively high nucleotide diversity and evolutionary rate, as well as the largest amount of phylogenetic information. These results indicate that nad2 may be a better marker to investigate phylogenetic relationships among closely related species in Copepoda than the commonly used cox1 gene. The sister-group relationship of Siphonostomatoida and Cyclopoida was recovered with strong support in our study. The only topological ambiguity was found within Cyclopoida, which might be caused by the rapid evolution and sparse taxon sampling of this lineage. More taxa and genes should be used to reconstruct the Copepoda phylogeny in the future.

Detecting gene breakpoints in noisy genome sequences using position-annotated colored de-Bruijn graphs

BACKGROUND

Identifying the locations of gene breakpoints between species of different taxonomic groups can provide useful insights into the underlying evolutionary processes. Given the exact locations of their genes, the breakpoints can be computed without much effort. However, often, existing gene annotations are erroneous, or only nucleotide sequences are available. Especially in mitochondrial genomes, high variations in gene orders are usually accompanied by a high degree of sequence inconsistencies. This makes accurately locating breakpoints in mitogenomic nucleotide sequences a challenging task.

RESULTS

This contribution presents a novel method for detecting gene breakpoints in the nucleotide sequences of complete mitochondrial genomes, taking into account possible high substitution rates. The method is implemented in the software package DeBBI. DeBBI allows to analyze transposition- and inversion-based breakpoints independently and uses a parallel program design, allowing to make use of modern multi-processor systems. Extensive tests on synthetic data sets, covering a broad range of sequence dissimilarities and different numbers of introduced breakpoints, demonstrate DeBBI 's ability to produce accurate results. Case studies using species of various taxonomic groups further show DeBBI 's applicability to real-life data. While (some) multiple sequence alignment tools can also be used for the task at hand, we demonstrate that especially gene breaks between short, poorly conserved tRNA genes can be detected more frequently with the proposed approach.

CONCLUSION

The proposed method constructs a position-annotated de-Bruijn graph of the input sequences. Using a heuristic algorithm, this graph is searched for particular structures, called bulges, which may be associated with the breakpoint locations. Despite the large size of these structures, the algorithm only requires a small number of graph traversal steps.

Comparative Mitogenomics of Jumping Spiders with First Complete Mitochondrial Genomes of Euophryini (Araneae: Salticidae)

Salticidae is the most species-rich family of spiders with diverse morphology, ecology and behavior. However, the characteristics of the mitogenomes within this group are poorly understood with relatively few well-characterized complete mitochondrial genomes. In this study, we provide completely annotated mitogenomes for Corythalia opima and Parabathippus shelfordi, which represent the first complete mitogenomes of the tribe Euophryini of Salticidae. The features and characteristics of the mitochondrial genomes are elucidated for Salticidae by thoroughly comparing the known well-characterized mitogenomes. The gene rearrangement between trnL2 and trnN was found in two jumping spider species, Corythalia opima and Heliophanus lineiventris Simon, 1868. Additionally, the rearrangement of nad1 to between trnE and trnF found in Asemonea sichuanensis Song & Chai, 1992 is the first protein-coding gene rearrangement in Salticidae, which may have an important phylogenetic implication for the family. Tandem repeats of various copy numbers and lengths were discovered in three jumping spider species. The codon usage analyses showed that the evolution of codon usage bias in salticid mitogenomes was affected by both selection and mutational pressure, but selection may have played a more important role. The phylogenetic analyses provided insight into the taxonomy of Colopsus longipalpis (Żabka, 1985). The data presented in this study will improve our understanding of the evolution of mitochondrial genomes within Salticidae.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Aythya marila

Aythya marila is a large diving duck belonging to the family Anatidae. However, the phylogenetic relationship among these Aythya species remains unclear due to the presence of extensive interspecific hybridization events within the Aythya genus. Here, we sequenced and annotated the complete mitochondrial genome of A. marila, which contained 22 tRNAs, 13 protein-coding genes (PCGs), 2 ribosomal RNAs, and 1 D-loop, with a length of 16,617 bp. The sizes of the PCGs ranged from 297 to 1824 bp and were all, except for ND6, located on the heavy chain (H). ATG and TAA were the most common start and termination codons of the 13 PCGs, respectively. The fastest- and slowest-evolving genes were ATP8 and COI, respectively. Codon usage analysis indicated that CUA, AUC, GCC, UUC, CUC, and ACC were the six most frequent codons. The nucleotide diversity values indicated a high level of genetic diversity in A. marila. F_ST analysis suggested a widespread gene exchange between A. baeri and A. nyroca. Moreover, phylogenetic reconstructions using the mitochondrial genomes of all available Anatidae species showed that, in addition to A. marila, four major clades among the Anatidae (Dendrocygninae, Oxyurinae, Anserinae, and Anatinae) were closely related to A. fuligula. Overall, this study provides valuable information on the evolution of A. marila and new insights into the phylogeny of Anatidae.

Genome skimming elucidates the evolutionary history of Octopoda

Phylogenies for Octopoda have, until now, been based on morphological characters or a few genes. Here we provide the complete mitogenomes and the nuclear 18S and 28S ribosomal genes of twenty Octopoda specimens, comprising 18 species of Cirrata and Incirrata, representing 13 genera and all five putative families of Cirrata (Cirroctopodidae, Cirroteuthidae, Grimpoteuthidae, Opisthoteuthidae and Stauroteuthidae) and six families of Incirrata (Amphitretidae, Argonautidae, Bathypolypodidae, Eledonidae, Enteroctopodidae, and Megaleledonidae) which were assembled using genome skimming. Phylogenetic trees were built using Maximum Likelihood and Bayesian Inference with several alignment matrices. All mitochondrial genomes had the 'typical' genome composition and gene order previously reported for octopodiforms, except Bathypolypus ergasticus, which appears to lack ND5, two tRNA genes that flank ND5 and two other tRNA genes. Argonautoidea was revealed as sister to Octopodidae by the mitochondrial protein-coding gene dataset, however, it was recovered as sister to all other incirrate octopods with strong support in an analysis using nuclear rRNA genes. Within Cirrata, our study supports two existing classifications suggesting neither is likely in conflict with the true evolutionary history of the suborder. Genome skimming is useful in the analysis of phylogenetic relationships within Octopoda; inclusion of both mitochondrial and nuclear data may be key.

The First Complete Mitochondrial Genome of Genus Isocapnia (Plecoptera: Capniidae) and Phylogenetic Assignment of Superfamily Nemouroidea

Capniidae are a family of stoneflies, also known as snow flies, who emerge in winter. The phylogeny of Capniidae is widely accepted to be based on morphological analysis. Until now, only five Capniidae mitochondrial genomes have been sequenced so far. In addition, sampling is required to determine an accurate phylogenetic association because the generic classification of this family is still controversial and needs to be investigated further. In this study, the first mitogenome of genus Isocapnia was sequenced with a length of 16,200 bp and contained 37 genes, including a control region, two rRNAs, 22 tRNAs, and 13 PCGs, respectively. Twelve PCGs originated with the common start codon ATN (ATG, ATA, or ATT), while nad5 used GTG. Eleven PCGs had TAN (TAA or TAG) as their last codon; however, cox1 and nad5 had T as their final codon due to a shortened termination codon. All tRNA genes demonstrated the cloverleaf structure, which is distinctive for metazoans excluding the tRNASer1 (AGN) that missed the dihydrouridine arm. A Phylogenetic analysis of the superfamily Nemouroidea was constructed using thirteen PCGs from 32 formerly sequenced Plecoptera species. The Bayesian inference and maximum likelihood phylogeny tree structures derived similar results across the thirteen PCGs. Our findings strongly supported Leuctridae + ((Capniidae + Taeniopterygidae) + (Nemouridae + Notonemouridae)). Ultimately, the best well-supported generic phylogenetic relationship within Capniidae is as follows; (Isocapnia + (Capnia + Zwicknia) + (Apteroperla + Mesocapnia)). These findings will enable us to better understand the evolutionary relationships within the superfamily Nemouroidea and the generic classification and mitogenome structure of the family Capniidae.

The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

The complete mitochondrial genome of Hemigrapsus sinensis (Brachyura, Grapsoidea, Varunidae) and its phylogenetic position within Grapsoidea

BACKGROUND

In this study, the complete mitogenome of Hemigrapsus sinensis was the first identified and analyzed.

OBJECTIVE

The complete mitochondrial genome of Hemigrapsus sinensis (Brachyura, Grapsoidea, Varunidae) and its phylogenetic position within Grapsoidea.

METHODS

The sample of Hemigrapsus sinensis was collected and DNA was extracted. After sequencing, NOVOPlasty was used for sequence assembly. Annotate sequences with MITOS WebServer, tRNAscan-SE2.0, and NCBI database. MEGA was used for sequence analysis and Phylosuite was used for phylogenetic tree construction. DnaSP was used to calculate Ka/Ks.

RESULTS

This mitochondrial genome shows that it was 15,900 bp and encoded 13 PCGs, 22 tRNA genes, two rRNA genes, and one control region. The genome composition tends to A + T (74.34%) and presents a negative GC-skew (- 0.22) and AT-skew (- 0.03). The PCGs initiation codon was the typical ATN and termination codon was the typical TAN, incomplete T or missing. The ML and BI trees showed that H. sinensis was most closely related to Hemigrapsus and clustered together with the Varunidae. And our phylogenetic trees provide proof that Ocypodoidea and Grapsoidea may be of common origin. Meanwhile, in the phylogenetic tree, parallel mixing of Chiromantes and Orisarma raised doubts over the traditional classification system. Besides, Incomplete Lineage sorting (ILS) was observed in Varunidae. In the subsequent analysis of evolution rate, we found that all of the PCGs (NAD4 was not calculated) had undergone negative selections, indicating the conservation of mitochondrial genes of H. sinensis during the evolution.

CONCLUSION

Therefore, researching the complete mitogenome of H. sinensis would be contributing to molecular taxonomy, phylogenetic relationship, and breeding optimization within the Grapsoidea superfamily.

A Blind-Identification Test on Criconema annuliferum (de Man, 1921) Micoletzky, 1925 Species Complex Corroborate the Hyper-Cryptic Species Diversity Using Integrative Taxonomy

Ring nematodes are obligate ectoparasites on crops and natural herbaceous and woody plants, and some species are of economic importance and cause damage to roots of several crops. Recent integrative taxonomical analyses recognized the existence of two cryptic species within the Criconema annuliferum morphotype in Spain. In this study, we corroborated that morphometric, morphological and a multi-locus analysis (including the ribosomal markers D2-D3 expansion segments of 28S rRNA, ITS rRNA, 18S RNA, and the mitochondrial DNA cytochrome oxidase I gene) identified a new lineage clearly separated from C. annuliferum, C. paraannuliferum and C. plesioannuliferum. The new lineage was described herein as Criconema pseudoannuliferum sp. nov., confirming that C. annuliferum species complex species complex comprises a hyper-cryptic species complex. This research analysed soil samples from the rhizosphere of maritime pine (Pinus pinaster Ait.) forests in Bermeja-Crestellina Mountain, located at the western part of Málaga province, southern Spain. The integrative taxonomical analyses revealed the occurrence of a new cryptic species identified using females, males and juveniles with detailed morphology, morphometry and molecular markers, described herein as Criconema pseudoannuliferum sp. nov. All molecular markers (D2-D3, ITS, 18S and COI) were obtained from the same individual that was also used for morphological and morphometric analyses. This research demonstrated the hidden diversity within the C. annuliferum species complex species complex can reach to four lineages under ribosomal and mitochondrial gene markers for one morphospecies group, which includes four species, viz. C. annuliferum, C. paraannuliferum, C. plesioannuliferum, and C. pseudoannuliferum sp. nov. Criconema pseudoannuliferum sp. nov. was detected in moderate soil density in two maritime pine forests (5 and 25 nematodes/500 cm³ of soil) suggesting that does not cause damage to maritime pine.

mtR_find: A Parallel Processing Tool to Identify and Annotate RNAs Derived from the Mitochondrial Genome

RNAs originating from mitochondrial genomes are abundant in transcriptomic datasets produced by high-throughput sequencing technologies, primarily in short-read outputs. Specific features of mitochondrial small RNAs (mt-sRNAs), such as non-templated additions, presence of length variants, sequence variants, and other modifications, necessitate the need for the development of an appropriate tool for their effective identification and annotation. We have developed mtR_find, a tool to detect and annotate mitochondrial RNAs, including mt-sRNAs and mitochondria-derived long non-coding RNAs (mt-lncRNA). mtR_find uses a novel method to compute the count of RNA sequences from adapter-trimmed reads. When analyzing the published datasets with mtR_find, we identified mt-sRNAs significantly associated with the health conditions, such as hepatocellular carcinoma and obesity, and we discovered novel mt-sRNAs. Furthermore, we identified mt-lncRNAs in early development in mice. These examples show the immediate impact of miR_find in extracting a novel biological information from the existing sequencing datasets. For benchmarking, the tool has been tested on a simulated dataset and the results were concordant. For accurate annotation of mitochondria-derived RNA, particularly mt-sRNA, we developed an appropriate nomenclature. mtR_find encompasses the mt-ncRNA transcriptomes in unpreceded resolution and simplicity, allowing re-analysis of the existing transcriptomic databases and the use of mt-ncRNAs as diagnostic or prognostic markers in the field of medicine.

Nanopore adaptive sampling for targeted mitochondrial genome sequencing and bloodmeal identification in hematophagous insects

BACKGROUND

Blood-feeding insects are important vectors for an array of zoonotic pathogens. While previous efforts toward generating molecular resources have largely focused on major vectors of global medical and veterinary importance, molecular data across a large number of hematophagous insect taxa remain limited. Advancements in long-read sequencing technologies and associated bioinformatic pipelines provide new opportunities for targeted sequencing of insect mitochondrial (mt) genomes. For engorged hematophagous insects, such technologies can be leveraged for both insect mitogenome genome assembly and identification of vertebrate blood-meal sources.

METHODS

We used nanopore adaptive sampling (NAS) to sequence genomic DNA from four species of field-collected, blood-engorged mosquitoes (Aedes and Culex spp.) and one deer fly (Chrysops sp.). NAS was used for bioinformatical enrichment of mtDNA reads of hematophagous insects and potential vertebrate blood-meal hosts using publically available mt genomes as references. We also performed an experimental control to compare results of traditional non-NAS nanopore sequencing to the mt genome enrichment by the NAS method.

RESULTS

Complete mitogenomes were assembled and annotated for all five species sequenced with NAS: Aedes trivittatus, Aedes vexans, Culex restuans, Culex territans and the deer fly, Chrysops niger. In comparison to data generated during our non-NAS control experiment, NAS yielded a substantially higher proportion of reference-mapped mtDNA reads, greatly streamlining downstream mitogenome assembly and annotation. The NAS-assembled mitogenomes ranged in length from 15,582 to 16,045 bp, contained between 78.1% and 79.0% A + T content and shared the anticipated arrangement of 13 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs. Maximum likelihood phylogenies were generated to further characterize each insect species. Additionally, vertebrate blood-meal analysis was successful in three samples sequenced, with mtDNA-based phylogenetic analyses revealing that blood-meal sources for Chrysops niger, Culex restuans and Aedes trivittatus were human, house sparrow (Passer domesticus) and eastern cottontail rabbit (Sylvilagus floridanus), respectively.

CONCLUSIONS

Our findings show that NAS has dual utility to simultaneously molecularly identify hematophagous insects and their blood-meal hosts. Moreover, our data indicate NAS can facilitate a wide array of mitogenomic systematic studies through novel 'phylogenetic capture' methods. We conclude that the NAS approach has great potential for broadly improving genomic resources used to identify blood-feeding insects, answer phylogenetic questions and elucidate complex pathways for the transmission of vector-borne pathogens.

Next-Generation Sequencing of the Complete Huaibei Grey Donkey Mitogenome and Mitogenomic Phylogeny of the Equidae Family

The Huaibei grey donkey (HGD) is an endangered species and a vital native breed in Anhui Province, China. However, its complete mitogenome, phylogeny, and maternal origin remain unclear. The objectives of this study were to detect the genetic diversity of the HGD and investigate its phylogenetic relationship with other breeds to inform conservation management. The complete mitogenome of the HGD was sequenced through next-generation sequencing, and the most variable region in the mitochondrial DNA displacement-loop (D-loop) was amplified via a polymerase chain reaction (PCR). Next, we used the median-joining network (MJN) to calculate the genetic relationships among populations and the neighbor-jointing method to build a phylogenetic tree and speculate as to its origin. The results showed that the mitogenome contains 22 tRNAs, 2 rRNAs, 13 PCGs, and 1 D-loop region. Analyzing the D-loop region of the HGDs, we identified 23 polymorphic sites and 11 haplotypes. The haplotype and nucleotide diversity were 0.87000 (Hd) and 0.02115 (Pi), respectively. The MJN analysis indicated that the HGD potentially has two maternal lineages, and phylogenetic analysis indicated that the Somali lineage could be the most probable domestication center for this breed. Therefore, our mitogenome analysis highlights the high genetic diversity of the HGD, which may have originated from the Somali wild ass, as opposed to the Asian wild ass. This study will provide a useful resource for HGD conservation and breeding.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei

The enigmatic scaphopods, or tusk shells, are a small and rare group of molluscs whose phylogenomic position among the Conchifera is undetermined, and the taxonomy within this class also needs revision. Such work is hindered by there only being a very few mitochondrial genomes in this group that are currently available. Here, we present the assembly and annotation of the complete mitochondrial genome from Dentaliida Pictodentalium vernedei, whose mitochondrial genome is 14,519 bp in size, containing 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The nucleotide composition was skewed toward A-T, with a 71.91% proportion of AT content. Due to the mitogenome-based phylogenetic analysis, we defined P. vernedei as a sister to Graptacme eborea in Dentaliida. Although a few re-arrangements occurred, the mitochondrial gene order showed deep conservation within Dentaliida. Yet, such a gene order in Dentaliida largely diverges from Gadilida and other molluscan classes, suggesting that scaphopods have the highest degree of mitogenome arrangement compared to other molluscs.

Characterization of four mitochondrial genomes of Crambidae (Lepidoptera, Pyraloidea) and phylogenetic implications

Loxostege turbidalis, Loxostege aeruginalis, Pyrausta despicata, and Crambus perlellus belong to Crambidae, Pyraloidea. Their mitochondrial genomes (mitogenomes) were successfully sequenced. The mitogenomes of L. turbidalis, L. aeruginalis, P. despicata, and C. perlellus are 15 240 bp, 15 339 bp, 15 389 bp, and 15 440 bp. The four mitogenomes all have a typical insect mitochondrial gene order, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and one A + T rich region (control region). The PCGs are initiated by the typical ATN codons, except CGA for the cox1 gene. Most PCGs terminate with common codon TAA or TAG, the incomplete codon T is found as the stop codon for cox2, nad4, and nad5. Most tRNA genes exhibit typical cloverleaf structure, except trnS1 (AGN) lacking the dihydrouridine (DHU) arm. The secondary structure of rRNA of four mitogenomes were predicted. Poly-T structure and micro-satellite regions are conserved in control regions. The phylogenetic analyses based on 13 PCGs showed the relationships of subfamilies in Pyraloidea. Pyralidae, and Crambidae are monophyletic, respectively. Pyralidae comprises four subfamilies, which form the following topology with high support values: (Galleriinae + ((Pyralinae + Epipaschiinae)+ Phycitinae)). Crambidae includes seven subfamilies and is divided into two lineages. Pyraustinae and Spilomelinae are sister groups of each other, and form the "PS clade." Other five subfamilies (Crambinae, Acentropinae, Scopariinae, Schoenobiinae, and Glaphyriinae) form the "non-PS clade" in the Bayesian inference tree. However, Schoenobiinae is not grouped with the other four subfamilies and located at the base of Crambidae in two maximum likelihood trees.

Characterization of mitochondrial genome of Indian Ocean blue-spotted maskray, Neotrygon indica and its phylogenetic relationship within Dasyatidae Family

The present study characterized complete mitochondrial genome of Blue-spotted maskray, Neotrygon indica and studied the evolutionary relationship of the species within the Dasyatidae family. The total length of the mitogenome was 17,974 bp including 37 genes and a non-coding control region. The average frequency of nucleotides in protein-coding genes was A: 29.1 %, T: 30.2 %, G: 13.0 % and C: 27.7 % with AT content of 59.3 %. The values of AT and GC skewness were -0.018 and -0.338, respectively. Comparative analyses showed a large number of average synonymous substitutions per synonymous site (Ks) in gene NADH4 (5.07) followed by NADH5 (4.72). High values of average number of non-synonymous substitutions per non-synonymous site (Ka) were observed in genes ATPase8 (0.54) and NADH2 (0.44). Genes NADH4L and NADH2 showed high interspecific genetic distance values of 0.224 ± 0.001 and 0.213 ± 0.002, respectively. Heat map analysis showed variation in codon usage among different species of the Dasyatidae family. The phylogenetic tree showed a sister relationship between the Dasyatinae and the Neotrygoninae subfamilies. Neotrygon indica formed as a sister species to the clade consisting of N. varidens and N. orientalis. Based on the present results, Neotrygon indica could have diverged from the common ancestor of the two latter in the Plio-Pleistocene. The present study showed distinct characteristics of N. indica from its congeners through comparative mitogenomics.

Mitogenomic Features and Evolution of the Nile River Dominant Tilapiine Species (Perciformes: Cichlidae)

To better understand the diversity and evolution of cichlids, we sequenced, assembled, and annotated the complete mitochondrial genomes of three Nile tilapiine species (Coptodon zillii, Oreochromis niloticus, and Sarotherodon galilaeus) dominating the Nile River waters. Our results showed that the general mitogenomic features were conserved among the Nile tilapiine species. The genome length ranged from 16,436 to 16,631 bp and a total of 37 genes were identified (two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and 1 control region). The ND6 was the only CDS that presented a negative AT skew and a positive GC skew. The most extended repeat sequences were in the D-loop followed by the pseudogenes (trnS^GCU). The ND5 showed relatively high substitution rates whereas ATP8 had the lowest substitution rate. The codon usage bias displayed a greater quantity of NNA and NNC at the third position and anti-bias against NNG. The phylogenetic relationship based on the complete mitogenomes and CDS was able to differentiate the three species as previously reported. This study provides new insight into the evolutionary connections between various subfamilies within cichlids while providing new molecular data that can be applied to discriminate between Nile tilapiine species and their populations.