The assembled Banana dihaploid mitochondrial genome is compact with a high number of gene copies

Banana being a major food crop all around the world, attracts various research interests in crop improvement. In banana, complete genome sequences of Musa accuminata and Musa balbisiana are available. However, the mitochondrial genome is not sequenced or assembled. Mitochondrial (mt) genes play an important role in flower and seed development and in Cytoplasmic Male Sterility. Unraveling banana mt genome architecture will be a foundation for understanding inheritance of traits and their evolution. In this study, the complete banana mt genome is assembled from the whole genome sequence data of the Musa acuminata subsp. malaccensis DH-Pahang. The mt genome sequence acquired by this approach was 409574 bp and it contains, 54 genes coding for 25 respiratory complex proteins 15 ribosomal proteins, 12 tRNA genes and two ribosomal RNA gene. Except atpB, rps11 and rps19 other genes are in multiple copies. The copy number is 12 in tRNA genes. In addition, nearly 25% tandem repeats are also present in it. These mt proteins are identical to the mt proteins present in the other members of AA genome and share 98% sequence similarity with M. balbisiana. The C to U RNA editing is profoundly higher (87 vs 13%) in transcripts of M. balbisiana (BB) compared to M. accuminata (AA). The banana AA mitochondrial genome is tightly packed with 233 genes, with less rearrangements and just 5.3% chloroplast DNA in it. The maintenance of high copy number of functional mt genes suggest that they have a crucial role in the evolution of banana.

Fragmentation in mitochondrial genomes in relation to elevated sequence divergence and extreme rearrangements

Background A single circular mitochondrial (mt) genome is a common feature across most metazoans. The mt-genome includes protein-coding genes involved in oxidative phosphorylation, as well as RNAs necessary for translation of mt-RNAs, whose order and number are highly conserved across animal clades, with few known exceptions of alternative mt-gene order or mt-genome architectures. One such exception consists of the fragmented mitochondrial genome, a type of genome architecture where mt-genes are split across two or more mt-chromosomes. However, the origins of mt-genome fragmentation and its effects on mt-genome evolution are unknown. Here, we investigate these origin and potential mechanisms underlying mt-genome fragmentation, focusing on a genus of booklice, Liposcelis, which exhibits elevated sequence divergence, frequent rearrangement of mt-gene order, and fragmentation of the mt genome, and compare them to other Metazoan clades. Results We found this genus Liposcelis exhibits very low conservation of mt-gene order across species, relative to other metazoans. Levels of gene order rearrangement were, however, unrelated to whether or not mt-genomes were fragmented or intact, suggesting mitochondrial genome fragmentation is not affecting mt-gene order directly. We further investigated possible mechanisms underpinning these patterns and revealed very high conservation of non-coding sequences at the edges of multiple recombination regions across populations of one particular Liposcelis species, supportive of a hypothesis that mt-fragmentation arises from recombination errors between mt-genome copies. We propose these errors may arise as a consequence of a heightened mutation rate in clades exhibiting mt-fragmentation. Consistent with this, we observed a striking pattern across three Metazoan phyla (Arthropoda, Nematoda, Cnidaria) characterised by members exhibiting high levels of mt-gene order rearrangement and cases of mt-fragmentation, whereby the mt-genomes of species more closely related to species with fragmented mt-genomes diverge more rapidly despite experiencing strong purifying selection. Conclusions We showed that contrary to expectations, mt-genome fragmentation is not correlated with the increase in mt-genome rearrangements. Furthermore, we present evidence that fragmentation of the mt-genome may be part of a general relaxation of a natural selection on the mt-genome, thus providing new insights into the origins of mt-genome fragmentation and evolution.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

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

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

Hookworm infection in central China: morphological and molecular diagnosis

Background Necator americanus is one of the major etiological agents of human ancylostomiasis. Historically, the epidemiology of ancylostomiasis in Henan Province of central China and the molecular characteristics of N. americanus have been poorly understood. Methods In this study, we report a case of ancylostomiasis in Zhengzhou city of Henan Province. We also review the epidemiology of ancylostomiasis in Henan Province from 1949 to 2020. In addition, the complete mitochondrial (mt) genome of one clinical isolate is fully characterized using Illumina sequencing. All available mt genomes of hookworms in GenBank were included to reconstruct the phylogeny using both maximum likelihood (ML) and Bayesian inference (BI) methods. Results A total of three worms were collected from the patient. These worms were identified as N. americanus based on morphological characteristics as well as confirmed by genotyping with the barcoding gene cox1. Although ancylostomiasis cases have dropped substantially in recent years, hookworm infection is still a public health problem in underdeveloped areas and remote rural areas in Henan Province. The mt genome features of the N. americanus contained 12 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a major non-coding region. The nad1 gene showed high sequence variability among isolates, which is worth considering for future genetic studies of N. americanus. Phylogenetic analyses support the monophyly of hookworm isolates from different hosts and distinct geographical locations. Conclusions The mt genome of N. americanus presented here will serve as a useful data set for studying population genetics and phylogenetic relationships of hookworms. Positive measures for preventing and controlling ancylostomiasis are required by both health services and individuals in Henan Province. Graphical abstract

Frequent tRNA gene translocation towards the boundaries with control regions contributes to the highly dynamic mitochondrial genome organization of the parasitic lice of mammals

Background The typical single-chromosome mitochondrial (mt) genome of animals has fragmented into multiple minichromosomes in the lineage Mitodivisia, which contains most of the parasitic lice of eutherian mammals. These parasitic lice differ from each other even among congeneric species in mt karyotype, i.e. the number of minichromosomes, and the gene content and gene order in each minichromosome, which is in stark contrast to the extremely conserved single-chromosome mt genomes across most animal lineages. How fragmented mt genomes evolved is still poorly understood. We use Polyplax sucking lice as a model to investigate how tRNA gene translocation shapes the dynamic mt karyotypes. Results We sequenced the full mt genome of the Asian grey shrew louse, Polyplax reclinata. We then inferred the ancestral mt karyotype for Polyplax lice and compared it with the mt karyotypes of the three Polyplax species sequenced to date. We found that tRNA genes were entirely responsible for mt karyotype variation among these three species of Polyplax lice. Furthermore, tRNA gene translocation observed in Polyplax lice was only between different types of minichromosomes and towards the boundaries with the control region. A similar pattern of tRNA gene translocation can also been seen in other sucking lice with fragmented mt genomes. Conclusions We conclude that inter-minichromosomal tRNA gene translocation orientated towards the boundaries with the control region is a major contributing factor to the highly dynamic mitochondrial genome organization in the parasitic lice of mammals.

Insights into molecular structure, genome evolution and phylogenetic implication through mitochondrial genome sequence of Gleditsia sinensis

Gleditsia sinensis is an endemic species widely distributed in China with high economic and medicinal value. To explore the genomic evolution and phylogenetic relationships of G. sinensis, the complete mitochondrial (mt) genome of G. sinensis was sequenced and assembled, which was firstly reported in Gleditsia. The mt genome was circular and 594,121 bp in length, including 37 protein-coding genes (PCGs), 19 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. The overall base composition of the G. sinensis mt genome was 27.4% for A, 27.4% for T, 22.6% for G, 22.7% for C. The comparative analysis of PCGs in Fabaceae species showed that most of the ribosomal protein genes and succinate dehydrogenase genes were lost. In addition, we found that the rps4 gene was only lost in G. sinensis, whereas it was retained in other Fabaceae species. The phylogenetic analysis based on shared PCGs of 24 species (22 Fabaceae and 2 Solanaceae) showed that G. sinensis is evolutionarily closer to Senna species. In general, this research will provide valuable information for the evolution of G. sinensis and provide insight into the phylogenetic relationships within the family Fabaceae.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts contain diverse clades within the kingdom of Fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 141 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions' diversity, genetic code alterations and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout Saccharomycotina evolution. Overall, this work shows that phylogenetic studies that employ the mt genome of yeasts highlight major evolutionary events.

Precise annotation of Drosophila mitochondrial genomes leads to insights into AT-rich regions

In the present study, we performed precise annotation of Drosophila melanogaster, D. simulans, D. grimshawi, Bactrocera oleae mitochondrial (mt) genomes by pan RNA-seq analysis. Our new annotations corrected or modified some of the previous annotations and two important findings were reported for the first time, including the discovery of the conserved polyA(+) and polyA(-) motifs in the control regions (CRs) of insect mt genomes and the adding of CCAs to the 3' ends of two antisense tRNAs in D. melanogaster mt genome. Using PacBio cDNA-seq data from D. simulans, we precisely annotated the Transcription Initiation Sites (TISs) of the mt Heavy and Light strands in Drosophila mt genomes and reported that the polyA(+) and polyA(-) motifs in the CRs are associated with TISs. The discovery of the conserved polyA(+) and polyA(-) motifs provides insights into many polyA and polyT sequences in CRs of insect mt genomes, leading to reveal the mt transcription and its regulation in invertebrates. In addition, we provided a high-quality, well-curated and precisely annotated D. simulans mt genome (GenBank: MN611461), which should be included into the NCBI RefSeq database to replace the current reference genome NC_005781.