Genetic aspects of mitochondrial genome evolution

Daphnia sp. (Branchiopoda: Cladocera) Mitochondrial Genome Gene Rearrangement and Phylogenetic Position Within Branchiopoda

In high-altitude (4500 m) freshwater lakes, Daphnia is the apex species and the dominant zooplankton. It frequently dwells in the same lake as the Gammarid. Branchiopoda, a class of Arthropoda, Crustacea, is a relatively primitive group in the subphylum Crustacea, which originated in the Cambrian period of the Paleozoic. The complete mitogenome sequence of Daphnia sp. (Branchiopoda: Cladocera) was sequenced and annotated in this study and deposited in GenBank. The sequence structure of this species was studied by comparing the original sequences with BLAST. In addition, we have also researched the mechanisms of their mitochondrial gene rearrangement by establishing a model. We have used the Bayesian inference [BI] and maximum likelihood [ML] methods to proceed with phylogenetic analysis inference, which generates identical phylogenetic topology that reveals the phylogenetic state of Daphnia. The complete mitogenome of Daphnia sp. shows that it was 15,254 bp in length and included two control regions (CRs) and 37 genes (13 protein-coding genes, 22 tRNAs and two ribosomal RNAs [16S and 12S]). In addition to tRNA-Ser (GCT), other tRNAs have a typical cloverleaf secondary structure. Meanwhile, the mitogenome of Daphnia sp. was clearly rearranged when compared to the mitogenome of typical Daphnia. In a word, we report a newly sequenced mitogenome of Daphnia sp. with a unique rearrangement phenomenon. These results will be helpful for further phylogenetic research and provide a foundation for future studies on the characteristics of the mitochondrial gene arrangement process in Daphnia.

The absence of canonical respiratory complex I subunits in male-type mitogenomes of three Donax species

Bivalves are an extraordinary class of animals in which species with a doubly uniparental inheritance (DUI) of mitochondrial DNA have been described. DUI is characterized as a mitochondrial homoplasmy of females and heteroplasmy of male individuals where F-type mitogenomes are passed to the progeny with mother egg cells and divergent M-type mitogenomes are inherited with fathers sperm cells. However, in most cases only male individuals retain divergent mitogenome inherited with spermatozoa. Additionally, in many of bivalves, unique mitochondrial features, like additional genes, gene duplication, gene extensions, mitochondrial introns, and recombination, were observed. In this study, we sequenced and assembled male-type mitogenomes of three Donax species. Comparative analysis of mitochondrial sequences revealed a lack of all seven NADH dehydrogenase subunits as well as the presence of three long additional open reading frames lacking identifiable homology to any of the existing genes.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes

BACKGROUND

Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes.

METHODS

Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis.

RESULTS

The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%.

CONCLUSION

The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Characterization of the complete mitochondrial genome of Spirobolus grahami (Diplopoda: Spirobolidae) with phylogenetic analysis

Diplopoda is one of the most diverse and important groups of soil arthropods, but little research has been done on their phylogenetic relationship and evolution. Here, we sequenced and annotated the complete mitochondrial genomes of Spirobolus grahami. The total mitogenome of S. grahami was typical circular, double-stranded molecules, with 14,875 bp in length, including 13 protein-coding genes, 22 tRNAs, two rRNAs, and one control region. Base composition analysis suggested that the mitochondrial sequences were biased toward A and T, with A + T content of 58.68%. The mitogenomes of S. grahami exhibited negative AT and positive GC skews. Most of the 13 PCGs had ATN as the start codon, except COX1 start with CGA, and most PCGs ended with the T stop codon. The dN/dS values for most PCGs were lower than 1, suggesting that purifying selection was likely the main driver of mitochondrial PCG evolution. Phylogenetic analyses based on 13 PCGs using BI and ML methods support the classification of genus Spirobolus and Tropostreptus. Glomeridesmus spelaeus is distantly related to the other Diplopoda species.

First mitogenomic characterization of Macromotettixoides (Orthoptera, Tetrigidae), with the descriptions of two new species

Classification of species is commonly based on morphological, molecular, and distribution depending on the taxa. Macromotettixoides Zheng, Wei & Jiang, 2005 and Formosatettix Tinkham, 1937 are both wingless types of Tetrigidae with extremely similar morphological characteristics, and in the current taxonomic system they are placed in two different subfamilies, Metrodorinae and Tetriginae, respectively. It is difficult to clearly identify the species of these two genera by morphological characteristics, and molecular data is often needed to assist identification. Here, the complete mitogenomes of two new species were sequenced and assembled, with that of Macromotettixoidesorthomargina. Molecular data of species of Formosatettix were used to test the monophyly of Macromotettixoides and to re-assess the generic characters, and also to test whether Macromotettixoides belongs to the Asian Metrodorinae or Tetriginae. Furthermore, mitochondrial characteristics were analyzed and the phylogeny of the Tetrigidae reconstructed based on mitochondrial protein-coding genes (PCGs). The results indicated that the two new species were clustered with Macromotettixoides rather than Formosatettix, and the anterior margin of the fastigium and pronotum of the two new species usually had the humeral angle different from that of Formosatettix. Therefore, after integrating morphological and molecular data, the two new species were placed in the genus Macromotettixoides, M.maoershanensissp. nov. and M.brachycornasp. nov. Finally, a phylogenetic reconstruction supported Macromotettixoides being assigned to Tetriginae rather than Metrodorinae, in contrast to the previous classification of this genus.

Complete mitogenome of Lepidopygopsis typus, an evolutionarily-distinct, endangered cyprinid fish from the Western Ghats Biodiversity Hotspot: Phylogenetic relationships and implications for conservation

The mitogenome of Lepidopygopsis typus, an evolutionarily distinct, endangered, cyprinid fish from the Western Ghats Biodiversity Hotspot, was characterized. Total length of the mitogenome was 16,729 bp, and comprised of 13 protein coding, 22 tRNA and two rRNA genes. Thirteen protein coding genes contributed a total nucleotide length of 11,408 bp, which coded for a total of 3794 amino acids and accounting for 68.2 % of the mitogenome. A maximum likelihood phylogenetic tree based on mitogenomes of members of cyprinid subfamilies Torinae, Schizothoracinae and Schizopygopsinae, recovered L. typus in a well-resolved clade of Torinae with members of the Western Ghats endemic genus Hypselobarbus as sister taxa. Selection analysis suggested that the branch of L. typus and Hypselobarbus spp. was under intensified selection with reference to other members of Torinae. There were 19 codons under diversifying selection in L. typus, which could be the result of positive selection for adapting to high-altitude, upstream tributaries of Periyar River, where the species has a restricted distribution. With respect to Hypselobarbus spp., L. typus was under relaxed selection with about 68 % of codons experiencing neutral evolution. Restricted distribution, low population size and relaxation in selection can likely trigger extinction in L. typus, and therefore urgent conservation and monitoring plans are required to secure the future of this evolutionary distinct, and globally endangered species.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

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.

microRNAs in Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma) possibly involved in the mitochondrial function

Mitochondria are essential organelles for maintaining vital cellular functions, and microRNAs (miRNAs) regulate gene expression posttranscriptionally. miRNAs exhibit tissue and time-specific patterns in mitochondria and specifically mitochondrial miRNAs (mitomiRs) can regulate the mRNA expression both originating from mitochondrial and nuclear transcription which affect mitochondrial metabolic activity and cell homeostasis. In this study, miRNAs of two insect species, Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma), were investigated for the first time. The known and possible novel miRNAs were predicted and characterized and their potential effects on mitochondrial transcription were investigated in these insect species using deep sequencing. The previously reported mitomiRs were also investigated and housekeeping miRNAs were characterized. miRNAs that are involved in mitochondrial processes such as apoptosis and signaling and that affect genes encoding the subunits of OXPHOS complexes have been identified in each species. Here, 81 and 161 novel mature miRNA candidates were bioinformatically predicted and 9 and 24 of those were aligned with reference mitogenomes of S. parreyssi and L. saccharina, respectively. As a result of RNAHybrid analysis, 51 and 69 potential targets of miRNAs were found in the mitogenome of S. parreyssi and L. saccharina, respectively. cox1 gene was the most targeted gene and cytB, rrnS, and rrnL genes were highly targeted in both of the species by novel miRNAs, hypothetically. We speculate that these novel miRNAs, originating from or targeting mitochondria, influence on rRNA genes or positively selected mitochondrial protein-coding genes. These findings may provide a new perspective in evaluating miRNAs for maintaining mitochondrial function and transcription.

Genetic and Haplotype Diversity of Manila Clam Ruditapes philippinarum in Different Regions of China Based on Three Molecular Markers

China has the largest production yield of Manila clam Ruditapes philippinarum in the world. Most of the clam seeds for aquaculture are mainly derived from artificial breeding in southern China, likely resulting in the loss of genetic variation and inbreeding depression. To understand the genetic and haplotype diversity of R. philippinarum, 14 clam populations sampled from different regions of China were analyzed by three molecular markers, including COI, 16SrRNA and ITS. Based on the results of the COI and ITS genes, the 14 populations showed a moderate to high level of genetic diversity, with an average haplotype diversity of 0.9242 and nucleotide diversity of 0.05248. AMOVA showed that there was significant genetic differentiation among all populations (mean FST of the total population was 0.4534). Pairwise FST analysis showed that genetic differentiation reached significant levels between Laizhou and other populations. Two Laizhou populations showed great divergence from other populations, forming an independent branch in the phylogenetic tree. The shared haplotypes Hap_2 and Hap_4 of COI appeared most frequently in most clam populations. In contrast, 16SrRNA analysis of the clam populations revealed the dominated haplotype Hap_2, accounting for 70% of the total number of individuals. The haplotype diversity of the Laizhou population (Laizhou shell-wide (KK) and Laizhou dock (LZMT)) was relatively higher than other populations, showing multiple unique haplotypes (e.g., Hap_40, Hap_41 and Hap_42). These findings of genetic and haplotype diversity of clam populations provide guiding information for genetic resource conservation and genetic improvement of the commercially important R. philippinarum.

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

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

Chromosome-level, nanopore-only genome and allele-specific DNA methylation of Pallas's cat, Otocolobus manul

Pallas's cat, or the manul cat (Otocolobus manul), is a small felid native to the grasslands and steppes of central Asia. Population strongholds in Mongolia and China face growing challenges from climate change, habitat fragmentation, poaching, and other sources. These threats, combined with O. manul's zoo collection popularity and value in evolutionary biology, necessitate improvement of species genomic resources. We used standalone nanopore sequencing to assemble a 2.5 Gb, 61-contig nuclear assembly and 17097 bp mitogenome for O. manul. The primary nuclear assembly had 56× sequencing coverage, a contig N50 of 118 Mb, and a 94.7% BUSCO completeness score for Carnivora-specific genes. High genome collinearity within Felidae permitted alignment-based scaffolding onto the fishing cat (Prionailurus viverrinus) reference genome. Manul contigs spanned all 19 felid chromosomes with an inferred total gap length of less than 400 kilobases. Modified basecalling and variant phasing produced an alternate pseudohaplotype assembly and allele-specific DNA methylation calls; 61 differentially methylated regions were identified between haplotypes. Nearest features included classical imprinted genes, non-coding RNAs, and putative novel imprinted loci. The assembled mitogenome successfully resolved existing discordance between Felinae nuclear and mtDNA phylogenies. All assembly drafts were generated from 158 Gb of sequence using seven minION flow cells.

Mitochondrial Genomes Assembled from Non-Invasive eDNA Metagenomic Scat Samples in Critically Endangered Mammals

The abundance of many large-bodied vertebrates, both in marine and terrestrial environments, has declined substantially due to global and regional climate stressors that define the Anthropocene. The development of genetic tools that can serve to monitor population’s health non-intrusively and inform strategies for the recovery of these species is crucial. In this study, we formally evaluate whether whole mitochondrial genomes can be assembled from environmental DNA (eDNA) metagenomics scat samples. Mitogenomes of four different large vertebrates, the panda bear (Ailuropoda melanoleuca), the moon bear (Ursus thibetanus), the Java pangolin (Manis javanica), and the the North Atlantic right whale (Eubalaena glacialis) were assembled and circularized using the pipeline GetOrganelle with a coverage ranging from 12x to 480x in 14 out of 18 different eDNA samples. Partial mitochondrial genomes were retrieved from three other eDNA samples. The complete mitochondrial genomes of the studied species were AT-rich and comprised 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a putative D-loop/control region. Synteny observed in all assembled mitogenomes was identical to that reported for specimens of the same and other closely related species. This study demonstrates that it is possible to assemble accurate whole mitochondrial chromosomes from eDNA samples (scats) using forthright bench and bioinformatics workflows. The retrieval of mitochondrial genomes from eDNA samples represents a tool to support bioprospecting, bio-monitoring, and other non-intrusive conservation strategies in species considered ‘vulnerable’, ‘endangered’, and/or ‘critically endangered’ by the IUCN Red List of Threatened Species.

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 novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Characterization of Two New Apodemus Mitogenomes (Rodentia: Muridae) and Mitochondrial Phylogeny of Muridae

Apodemus is the most common small rodent species in the Palearctic realm and an ideal species for biogeographical research and understanding environmental changes. Elucidating phylogenetic relationships will help us better understand species adaptation and genetic evolution. Due to its stable structure, maternal inheritance, and rapid evolution, the mitogenome has become a hot spot for taxonomic and evolutionary studies. In this research, we determined the mitochondrial genome of Apodemus agrarius ningpoensis and Apodemus draco draco and studied the phylogeny of Muridae using ML and BI trees based on all known complete mitogenomes. The mitochondrial genome of Apodemus agrarius ningpoensis was 16,262 bp, whereas that of Apodemus draco draco was 16,222 bp, and both encoded 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of base composition showed a clear A-T preference. All tRNAs except tRNASer and tRNALys formed a typical trilobal structure. All protein-coding genes contained T- and TAA as stop codons. Phylogeny analysis revealed two main branches in the Muridae family. Apodemus agrarius ningpoensis formed sister species with Apodemus chevrieri, whereas Apodemus draco draco with Apodemus latronum. Our findings provide theoretical basis for future studies focusing on the mitogenome evolution of Apodemus.

Argopistes sexvittatus and Argopistes capensis (Chrysomelidae: Alticini): Mitogenomics and Phylogeny of Two Flea Beetles Affecting Olive Trees

The genus Argopistes (Chrysomelidae: Alticini) is the only group of flea beetles specialized in plant hosts in the family Oleaceae. In southern Africa, Argopistes are often found feeding on African Wild Olive (Olea europaea subsp. cuspidata) and European cultivated olive (O. e. subsp. europaea), and heavy infestations can be devastating to mature trees and compromise the development of young trees. Despite their negative agricultural impact, African Argopistes are an understudied group for which no genetic data were available. We assessed the species diversity of olive flea beetles in the Western Cape province of South Africa, the largest olive-producing region in sub-Saharan Africa, by collecting adult specimens on wild and cultivated olive trees between 2015 and 2017. Argopistes sexvittatus Bryant, 1922 (n = 289) dominated at all sampling sites, and Argopistes capensis Bryant, 1944 (n = 2) was found only once. Argopistes oleae Bryant, 1922, a third species previously reported in the region, was not found. The complete mitogenomes of one A. capensis and two A. sexvittatus (striped and black morphotypes) individuals were sequenced for phylogenetic reconstruction in the context of other 64 species. The two olive flea beetle species form a monophyletic clade with other Argopistes, supporting the hypothesis that the exclusive feeding habit on Oleaceae is an evolutionary adaptation in this genus.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

The complete mitogenome of the Asian lung fluke Paragonimus skrjabini miyazakii and its implications for the family Paragonimidae (Trematoda: Platyhelminthes)

The complete circular mitogenome of Paragonimus skrjabini miyazakii (Platyhelminthes: Paragonimidae) from Japan, obtained by PacBio long-read sequencing, was 17 591 bp and contained 12 protein-coding genes (PCGs), 2 mitoribosomal RNA and 22 transfer RNA genes. The atp8 gene was absent, and there was a 40 bp overlap between nad4L and nad4. The long non-coding region (4.3 kb) included distinct types of long and short repeat units. The pattern of base usage for PCGs and the mtDNA coding region overall in Asian and American Paragonimus species (P. s. miyazakii, P. heterotremus, P. ohirai and P. kellicotti) and the Indian form of P. westermani was T > G > A > C. On the other hand, East-Asian P. westermani used T > G > C > A. Five Asian and American Paragonimus species and P. westermani had TTT/Phe, TTG/Leu and GTT/Val as the most frequently used codons, whereas the least-used codons were different in each species and between regional forms of P. westermani. The phylogenetic tree reconstructed from a concatenated alignment of amino acids of 12 PCGs from 36 strains/26 species/5 families of trematodes confirmed that the Paragonimidae is monophyletic, with 100% nodal support. Paragonimus skrjabini miyazakii was resolved as a sister to P. heterotremus. The P. westermani clade was clearly separate from remaining congeners. The latter clade was comprised of 2 subclades, one of the East-Asian and the other of the Indian Type 1 samples. Additional mitogenomes in the Paragonimidae are needed for genomic characterization and are useful for diagnostics, identification and genetic/ phylogenetic/ epidemiological/ evolutionary studies of the Paragonimidae.

Complete Mitogenome Analysis of Five Leafhopper Species of Idiocerini (Hemiptera: Cicadellidae)

Insect mitochondrial genomes (mitogenomes) are of great interest in exploring molecular evolution, phylogenetics, and biogeography. So far, only 12 mitogenomes of the leafhopper tribe Idiocerini have been released in GenBank, although the tribe comprises 488 known species including some agricultural, forestry, and horticultural pests. In order to compare and analyze the mitochondrial genome structure of Idiocerini and even the selective pressure of 13 protein-coding genes (PCGs) of the family Cicadellidae, the complete mitogenomes of five species including Nabicerus dentimus, Sahlbergotettix salicicola, Podulmorinus opacus, Podulmorinus consimilis, and a new species of a new genus were determined by next-generation sequencing. The size of the newly determined mitogenomes ranged from 14,733 bp to 15,044 bp, comprising the standard set of 13 PCGs, 22 transfer RNA genes, two ribosomal RNA genes, and a long non-coding control region (CR). The extent of purifying selection presented different pictures in the tribe and the family. The less pronounced genes (0.5 < dN/dS < 1) were nad5 and nad4l in Idiocerin, whereas in the family Cicadellidae including the sequences of Idiocerin, nad1-nad6 and cox1 genes were less pronounced. The codon encoding leucine was the most common in all species, and the codon encoding serine 1 was the most common in all species except for P. opacus. Interestingly, in P. opacus, another of the most common codons is that encoding serine 2. Among the 17 examined species of the Idiocerini, 14 species contained the tandem repeats, and 11 species of them contained the motif "TTATA". These findings will promote research on the structure and evolution of the mitochondrial genome and highlight the need for more mitogenomes in Cicadellidae.

Mitogenomic Comparison of the Mole Crickets Gryllotalpidae with the Phylogenetic Implications (Orthoptera: Ensifera)

Owing to limited molecular data, the phylogenetic position of the family Gryllotalpidae is still controversial in the infraorder Gryllidea. Mitochondrial genome (mitogenome) plays a crucial role in reconstructing phylogenetic relationships and revealing the molecular evolution of insects. However, only four mitogenomes have been reported in Gryllotalpidae to date. Herein, we obtained the first mitogenomes of Gryllotalpa henana Cai & Niu, 1998 and the Chinese G. orientalis Burmeister, 1838, made a detailed comparison of all mitogenomes available in Gryllotalpidae and reconstructed the phylogeny of Gryllidea based on mitogenomes using Bayesian inference (BI) and maximum likelihood (ML) methods. The results show that the complete mitogenome sequences of G. henana (15,504 bp) and G. orientalis (15,497 bp) are conserved, both exhibiting the double-stranded circular structure, typical gene content and the ancestral insect gene arrangement. The complete mitogenome of G.henana exhibits the lowest average AT content ever detected in Gryllotalpidae, and even Gryllidea. The gene nad2 of both species has atypical initiation codon GTG. All tRNAs exhibit typical clover-leaf structure, except for trnS1 lacking the dihydrouridine (DHU) arm. A potential stem-loop structure, containing a (T)_n(TC)₂(T)_n sequence, is detected in the control region of all gryllotalpids investigated and is likely related to the replication initiation of the minority strand. The phylogenetic analyses recover the six families of Gryllidea as Gryllotalpidae + (Myrmecophilidae + (Mogoplistidae + (Trigonidiidae + (Phalangopsidae + Gryllidae)))), similar to the trees based on transcriptomic and mitogenomic data. However, the trees are slightly different from the multilocus phylogenies, which show the sister-group relationship of Gryllotalpidae and Myrmecophilidae. The contradictions between mitogenomic and multilocus trees are briefly discussed.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

The Phylogenetic Position of the Enigmatic, Polypodium hydriforme (Cnidaria, Polypodiozoa): Insights from Mitochondrial Genomes

Polypodium hydriforme is an enigmatic parasite that belongs to the phylum Cnidaria. Its taxonomic position has been debated: whereas it was previously suggested to be part of Medusozoa, recent phylogenomic analyses based on nuclear genes support the view that P. hydriforme and Myxozoa form a clade called Endocnidozoa. Medusozoans have linear mitochondrial (mt) chromosomes, whereas myxozoans, as most metazoan species, have circular chromosomes. In this work, we determined the structure of the mt genome of P. hydriforme, using Illumina and Oxford Nanopore Technologies reads, and showed that it is circular. This suggests that P. hydriforme is not nested within Medusozoa, as this would entail linearization followed by recirculation. Instead, our results support the view that P. hydriforme is a sister clade to Myxozoa, and mt linearization in the lineage leading to medusozoans occurred after the divergence of Myxozoa + P. hydriforme. Detailed analyses of the assembled P. hydriforme mt genome show that: (1) it is encoded on a single circular chromosome with an estimated size of ∼93,000 base pairs, making it one of the largest metazoan mt genomes; (2) around 78% of the genome encompasses a noncoding region composed of several repeat types; (3) similar to Myxozoa, no mt tRNAs were identified; (4) the codon TGA is a stop codon and does not encode for tryptophan as in other cnidarians; (5) similar to myxozoan mt genomes, it is extremely fast evolving.

The complete mitochondrial genome of Pseudopleuronectes herzensteini

Pseudopleuronectes herzensteini belonging to Pleuronectiformes (family Pleuronectidae) is important in the fishery industry. However, the molecular biology of this valuable fish has hardly been reported. Thus, here we report the complete mitochondrial genome of P. herzensteini. The mitochondrial DNA (mtDNA) of P. herzensteini is 16,719 bp long and contains 13 mitochondrial protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a putative control region between tRNA-P and tRNA-F distinguished by a single short noncoding region. Phylogenetic analysis using PCGs confirmed that this mtDNA sequence belongs to the family Pleuronectidae. This is the first study reporting the complete mitochondrial genome sequence of P. herzensteini.

Genetic diversity of Diaphorina citri (Hemiptera: Liviidae) unravels phylogeographic structure and invasion history of eastern African populations

The Asian citrus psyllid (Diaphorina citri Kuwayama) is a key pest of Citrus sp. worldwide, as it acts as a vector for Candidatus Liberibacter asiaticus, the bacterial pathogen that causes citrus Huanglongbing. Diaphorina citri has been reported in Kenya, Tanzania, and more recently in Ethiopia. This study assessed the genetic diversity and phylogeographic structure of the pest to gain insights into the potential sources of its introduction into Africa. Population structure and differentiation of D. citri populations from China, Ethiopia, Kenya, Tanzania, and the USA were assessed using 10 microsatellite loci. Additionally, five new complete mitogenomes of D. citri collected in China, Ethiopia, Kenya, Tanzania, and the USA were analyzed in the context of publicly available sequences. Genotype data grouped the D. citri populations from Kenya and Tanzania in one cluster, and those from Ethiopia formed a separate cluster. The two genetic clusters inferred from genotype data were congruent with mitochondrial sequence data. The mitogenomes from Kenya/Tanzania/China had 99.0% similarity, and the Ethiopia/USA had 99.9% similarity. In conclusion, D. citri populations in eastern Africa have different sources, as the Kenyan and Tanzanian populations probably originated from southeastern Asia, while the Ethiopian population most probably originated from the Americas.

Endangered Schizothoracin Fish in the Tarim River Basin Are Threatened by Introgressive Hybridization

Simple Summary

Interspecific hybridization and introgression may threaten the survival of endangered species. Big-head Schizothoracin (Aspiorhynchus laticeps) and Tarim Schizothoracin (Schizothorax biddulphi) are two recognized types of fish of Schizothoracinae which belong to the family Cyprinidae. Big-head Schizothoracin and Tarim Schizothoracin are sympatrically distributed in the Tarim River basin in Xinjiang, China, and have been listed as first-class and second-class national key protected animals, respectively. Based on morphological characteristics, Schizothorax esocinus (another Schizothoracin that occurs in the Tarim River basin) is speculated to be hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin, but there is no direct genetic evidence for this point. In this study, the hybridization status of Schizothorax esocinus was confirmed by comparing genetic constitutions of mitochondrial DNA (mtDNA) and inter transcribed spacer (ITS2) encoded by the nuclear genome. Extensive hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin was detected. Since both Big-head Schizothoracin and Tarim Schizothoracin are critically endangered fishes, risk management due to hybridization is recommended to be considered as a part of current conservation programs.

Abstract

Big-head Schizothoracin (Aspiorhynchus laticeps) and Tarim Schizothoracin (Schizothorax biddulphi) are locally sympatric in the Tarim River Basin. Although another Schizothoracin (Schizothorax esocinus) in Tarim River basin has been speculated to be hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin, there was no genetic evidence. Previous studies on the genetics and evolution of Schizothoracins in Xinjiang Province were mostly based on mitochondrial DNA (mtDNA), whose characteristics of maternal inheritance made it hard to answer the question of whether there was hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin. In this study, cytochrome b (cytb) gene of mtDNA and internal transcribed spacer 2 (ITS2) that is encoded by the nuclear genome were genotyped within the entire samples at the same time. Our results confirmed that Schizothorax esocinus was the hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin. The heterozygous ITS2 genotypes and/or Aspiorhynchus laticeps-like mtDNA were also detected in a subset of samples that should have been identified as pure Schizothorax biddulphi based on morphology. The ITS2 is characterized by multi-copy, concert evolution, and biparental inheritance. Thus, by comparing with mtDNA data, broad-scale bidirectional hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin were revealed. Although interspecific hybridization may play a positive role in ecology and evolution, interspecific hybrids could threaten their parental species by the swamping of genetics and demography. As both parents of hybridization are critically endangered fishes, in this case, it is urgently necessary to strengthen the scientific assessment of the risks of the hybrids and the control of the hybridization and introgression between Aspiorhynchus laticeps and Schizothorax biddulphi in the Tarim River Basin.

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

Simple Summary

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

Abstract

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

Utilisation of Oxford Nanopore sequencing to generate six complete gastropod mitochondrial genomes as part of a biodiversity curriculum

High-throughput sequencing has enabled genome skimming approaches to produce complete mitochondrial genomes (mitogenomes) for species identification and phylogenomics purposes. In particular, the portable sequencing device from Oxford Nanopore Technologies (ONT) has the potential to facilitate hands-on training from sampling to sequencing and interpretation of mitogenomes. In this study, we present the results from sampling and sequencing of six gastropod mitogenomes (Aplysia argus, Cellana orientalis, Cellana toreuma, Conus ebraeus, Conus miles and Tylothais aculeata) from a graduate level biodiversity course. The students were able to produce mitogenomes from sampling to annotation using existing protocols and programs. Approximately 4 Gb of sequence was produced from 16 Flongle and one MinION flow cells, averaging 235 Mb and N50 = 4.4 kb per flow cell. Five of the six 14.1-18 kb mitogenomes were circlised containing all 13 core protein coding genes. Additional Illumina sequencing revealed that the ONT assemblies spanned over highly AT rich sequences in the control region that were otherwise missing in Illumina-assembled mitogenomes, but still contained a base error of one every 70.8-346.7 bp under the fast mode basecalling with the majority occurring at homopolymer regions. Our findings suggest that the portable MinION device can be used to rapidly produce low-cost mitogenomes onsite and tailored to genomics-based training in biodiversity research.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Sorting Signed Permutations by Inverse Tandem Duplication Random Losses

Gene order evolution of unichromosomal genomes, for example mitochondrial genomes, has been modelled mostly by four major types of genome rearrangements: inversions, transpositions, inverse transpositions, and tandem duplication random losses. Generalizing models that include all those rearrangements while admitting computational tractability are rare. In this paper, we study such a rearrangement model, namely the inverse tandem duplication random loss (iTDRL) model, where an iTDRL duplicates and inverts a continuous segment of a gene order followed by the random loss of one of the redundant copies of each gene. The iTDRL rearrangement has currently been proposed by several authors suggesting it to be a possible mechanisms of mitochondrial gene order evolution. We initiate the algorithmic study of this new model of genome rearrangement by proving that a shortest rearrangement scenario that transforms one given gene order into another given gene order can be obtained in quasilinear time. Furthermore, we show that the length of such a scenario, i.e., the minimum number of iTDRLs in the transformation, can be computed in linear time.

Identification of full-length circular nucleic acids using long-read sequencing technologies

Unlike the traditional perception in genomic DNA or linear RNA, circular nucleic acids are a class of functional biomolecules with a circular configuration and are often observed in nature. These circular molecules encompass the full spectrum of size and play an important role in organisms, making circular nucleic acids research worthy. Due to the low abundance of most types of circular nucleic acids and the disadvantages of short-read sequencing platforms, accurate and full-length circular nucleic acid sequencing and identification is difficult. In this review, we have provided insights into full-length circular nucleic acid detection methods using long-read sequencing technologies, with a focus on the experimental and bioinformatics strategies to obtain accurate sequences.

Species Diversity and Phylogenetic Relationships of Olive Lace Bugs (Hemiptera: Tingidae) Found in South Africa

Olive lace bugs (Hemiptera: Tingidae) are small sap-sucking insects that feed on wild and cultivated Olea europaea. The diversity of olive lace bug species in South Africa, the most important olive producer on the continent, has been incompletely surveyed. Adult specimens were collected in the Western Cape province for morphological and DNA-based species identification, and sequencing of complete mitogenomes. Cysteochila lineata, Plerochila australis, Neoplerochila paliatseasi and Neoplerochila sp. were found at 12 sites. Intra- and interspecific genetic divergences and phylogenetic clustering in 30 species in 18 genera of Tingidae using new and publicly available DNA barcodes showed high levels of congruity between taxonomic and genetic data. The phylogenetic position of the four species found in South Africa was inferred using new and available mitogenomes of Tingidae. Notably, olive lace bugs formed a cluster of closely related species. However, Cysteochila was non-monophyletic as C. lineata was recovered as a sister species to P. australis whereas Cysteochila chiniana, the other representative of the genus, was grouped with Trachypeplus jacobsoni and Tingis cardui in a different cluster. This result suggests that feeding on O. europaea may have a common origin in Tingidae and warrants future research on potential evolutionary adaptations of olive lace bugs to this plant host.

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Mitochondrial genomes of five Hyphessobrycon tetras and their phylogenetic implications

To date, the taxonomic status and phylogenetic affinities within Hyphessobrycon, even among other genera in Characidae, remain unclear. Here, we determined five new mitochondrial genomes (mitogenomes) of Hyphessobrycon species (H. elachys, H. flammeus, H. pulchripinnis, H. roseus, and H. sweglesi). The mitogenomes were all classical circular structures, with lengths ranging from 16,008 to 17,224 bp. The type of constitutive genes and direction of the coding strand that appeared in the mitogenomes were identical to those of other species in Characidae. The highest value of the Ka/Ks ratio within 13 protein-coding genes (PCGs) was found in ND2 with 0.83, suggesting that they were subject to purifying selection in the Hyphessobrycon genus. Comparison of the control region sequences among seven Hyphessobrycon fish revealed that repeat units differ in length and copy number across different species, which led to sharp differences in mitogenome sizes. Phylogenetic trees based on the 13 PCGs did not support taxonomic relationships, as the Hyphessobrycon fish mixed with those from other genera. These data were combined to explore higher level relationships within Characidae and could aid in the understanding of the evolution of this group.

Mitogenomes of Accipitriformes and Cathartiformes Were Subjected to Ancestral and Recent Duplications Followed by Gradual Degeneration

The rearrangement of 37 genes with one control region, firstly identified in Gallus gallus mitogenome, is believed to be ancestral for all Aves. However, mitogenomic sequences obtained in recent years revealed that many avian mitogenomes contain duplicated regions that were omitted in previous genomic versions. Their evolution and mechanism of duplication are still poorly understood. The order of Accipitriformes is especially interesting in this context because its representatives contain a duplicated control region in various stages of degeneration. Therefore, we applied an appropriate PCR strategy to look for duplications within the mitogenomes of the early diverged species Sagittarius serpentarius and Cathartiformes, which is a sister order to Accipitriformes. The analyses revealed the same duplicated gene order in all examined taxa and the common ancestor of these groups. The duplicated regions were subjected to gradual degeneration and homogenization during concerted evolution. The latter process occurred recently in the species of Cathartiformes as well as in the early diverged lineages of Accipitriformes, that is, Sagittarius serpentarius and Pandion haliaetus. However, in other lineages, that is, Pernis ptilorhynchus, as well as representatives of Aegypiinae, Aquilinae, and five related subfamilies of Accipitriformes (Accipitrinae, Circinae, Buteoninae, Haliaeetinae, and Milvinae), the duplications were evolving independently for at least 14-47 Myr. Different portions of control regions in Cathartiformes showed conflicting phylogenetic signals indicating that some sections of these regions were homogenized at a frequency higher than the rate of speciation, whereas others have still evolved separately.

Taxonomy classification using genomic footprint of mitochondrial sequences

BACKGROUND

Advancement in the sequencing technology yields a huge number of genomes of a multitude of organisms in our planet. One of the fundamental tasks for processing and analyzing these sequences is to organize them in the existing taxonomic orders. <P> Method: Recently we proposed a novel approach, GenFooT, of taxonomy classification using the concept of genomic footprint (GFP). The technique is further refined and enhanced in this work leading to improved accuracies in the task of taxonomic classification on various benchmark datasets. GenFooT maps a genome sequence in a 2D coordinate space and extracts features from that representation. It uses two hyper-parameters, namely block size and number of fragments of genomic sequence while computing the feature. In this work, we propose an analysis for choosing values of those parameters adaptively from the sequences. The enhanced version of GenFooT is named GenFooT2. <P> Results and Conclusion: We have experimented GenFooT2 on ten different biological datasets of genomic sequences of various organisms belonging to different taxonomy ranks. Our experimental results indicate more than 3% improved classification performance of the proposed features with Logistic regression classifier than the GenFooT. We also performed the statistical test to compare the performance of GenFooT2 with the state-of-the-art methods including our previous method GenFooT. The experimental results as well as the statistical test exhibit that the performance of the proposed GenFooT2 is significantly better.

Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

Complete mitochondrial genome of Black Drongo Dicrurus macrocercus (Passeriformes: Dicruridae) and phylogenetic analyses

The Black Drongo Dicrurus macrocercus is a bird belonging to the group Passeriformes, distributed almost all over the country in China. The conservation status of this species is Least Concern (LC) in IUCN. In this study, the complete mitogenome of D. macrocercus was determined. The mitochondrial DNA is packaged in a compact 17,056 based pair (bp) circular molecule with A + T content of 57.04%. It contains 37 typical mitochondrial genes, including 13 protein-coding genes, 2 rRNAs and 22 tRNAs, and 1 non-coding regions. We reconstructed a phylogenetic tree based on mitogenome sequences of 15 Dicruridae species and one outgroup. In trees, D. macrocercus was clustered as an independent clade with high support value (100). The new mitogenome data would provide useful information for application in conservation genetics and further clarify the phylogenetic evolution of this species.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

Complete mitochondrial DNA sequence of yellow-legged buttonquail (Turnix tanki)

The Yellow-legged Buttonquail Turnix tanki is a species of the genus Turnix, which belongs to the order Charadriiformes. It is distributed across almost all of China. The International Union for Conservation of Nature has assessed the bird's conservation status as 'Least Concern (LC).' We sequenced the complete mitogenome of T. tanki and examined its phylogenetic relationship with other charadriiformes species. The mitochondrial DNA is packaged in a compact 17,620 base pair circular molecule with A + T content of 57.90%. It contains 37 typical mitochondrial genes, including 13 protein-coding genes, two rRNAs and 22 tRNAs, and two non-coding regions. We reconstructed a phylogenetic tree based on mitogenome sequences of five Turnicidae species and one outgroup. Phylogenetic analysis indicated that T. tanki is a sister to T. suscitator.

How are the mitochondrial genomes reorganized in Hexapoda? Differential evolution and the first report of convergences within Hexapoda

The evolution of the Hexapoda mitochondrial genome has been the focus of several genetic and evolutionary studies over the last decades. However, they have concentrated on certain taxonomic orders of economic or health importance. The recent increase of mitochondrial genomes sequencing of diverse taxonomic orders generates an important opportunity to clarify the evolution of this group of organisms. However, there is no comparative study that investigates the evolution of the Hexapoda mitochondrial genome. In order to verify the level of rearrangement and the mitochondrial genome evolution, we performed a comparative genomic analysis of the Hexapoda mitochondrial genome available in the NCBI database. Using a combination of bioinformatics methods to carefully examine the mitochondrial gene rearrangements in 1198 Hexapoda species belonging to 32 taxonomic orders, we determined that there is a great variation in the rate of rearrangement by gene and by taxonomic order. A higher rate of genetic reassortment is observed in Phthiraptera, Thysanoptera, Protura, and Hymenoptera; compared to other taxonomic orders. Twenty-four events of convergence in the genetic order between different taxonomic orders were determined, most of them not previously reported; which proves the great evolutionary dynamics within Hexapoda.

Genetic variability, cryptic species and phylogenetic relationship of six cyathostomin species based on mitochondrial and nuclear sequences

Cyathostomins are important intestinal nematode parasites of equines and include 50 accepted species. Their taxonomy has been frequently revised and the presence of cryptic species suggested. Furthermore, usually molecular- and morphology-based phylogenetic analyses give divergent results. In this study, the nucleotide sequences of the nuclear second internal transcribed spacer (ITS-2) and the mitochondrial partial cytochrome c oxidase subunit I (COI) were determined for adults of six cyathostomin species (Coronocyclus coronatus, Coronocyclus labiatus, Cylicocyclus nassatus, Cylicostephanus calicatus, Cylicostephanus longibursatus, Cylicostephanus minutus) collected from different equine species within two geographic regions. Maximum likelihood trees were calculated for ITS-2, COI, and concatenated data. No obvious differentiation was observed between geographic regions or equine host species. As previously reported, Coronocyclus coronatus and Cylicostephanus calicatus revealed a close relationship. Cryptic species were detected in Cylicostephanus minutus and Cylicostephanus calicatus. Cylicocyclus nassatus and Coronocyclus labiatus showed diverse mitochondrial and nuclear haplotypes occurring in different combinations, while Cylicostephanus longibursatus was comparatively homogenous. In conclusion, a combined analysis of nuclear and mitochondrial haplotypes improved resolution of the phylogeny and should be applied to the remaining cyathostomin species and across additional equine host species and geographic regions.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Mitochondrial genome architecture and phylogenetic relationships of Odontesthes argentinensis within Atherinomorpha

Silversides are a widely distributed group across South America, with several species occupying marine, freshwater and estuarine environments. Several authors suggest main transitions among these environments took place during Pleistocene, and were accompanied with rapid speciation events. This scenario produced very limited genetic and morphological differentiation among the species. However, most of these surveys have an incomplete coverage of the intraspecific genetic diversity of the taxa studied. In this work, we reconstructed six mitochondrial genomes of O. argentinensis using transcriptomic data, and used them-in combination with several nuclear markers retrieved from the same transcriptomes-to explore the effect of additional coverage of intraspecific diversity of this species in phylogenetic reconstructions. Unlike previous works, phylogenetic analyses failed to identify O. argentinensis as a monophyletic group in relation with closely related taxa. Our results suggest that several species of the genus, especially those related to O. argentinensis, need further taxonomic revision.