Mitochondrial genome evolution and the origin of eukaryotes

Comparative analyses of Pleurotus pulmonarius mitochondrial genomes reveal two major lineages of mini oyster mushroom cultivars

Pleurotus pulmonarius, commonly known as the mini oyster mushroom, is highly esteemed for its crisp texture and umami flavor. Limited genetic diversity among P. pulmonarius cultivars raises concerns regarding its sustainable industrial production. To delve into the maternal genetic diversity of the principal P. pulmonarius cultivars, 36 cultivars and five wild isolates were subjected to de novo sequencing and assembly to generate high-quality mitogenome sequences. The P. pulmonarius mitogenomes had lengths ranging from 69,096 to 72,905 base pairs. The mitogenome sizes of P. pulmonarius and those of other mushroom species in the Pleurotus genus showed a significant positive correlation with the counts of LAGLIDAG and GIY-YIG homing endonucleases encoded by intronic open reading frames. A comparison of gene arrangements revealed an inversion of a fragment containing atp9-nad3-nad2 between P. pulmonarius and P. ostreatus. The mitogenomes of P. pulmonarius were clustered into three distinct clades, two of which were crowded with commercial cultivars. Clade I, all of which possess an inserted dpo gene, shared a maternal origin linked to an ancestral cultivar from Taiwan. Primers were designed to target the dpo gene, potentially safeguarding intellectual property rights. The wild isolates in Clade III exhibited more divergent mitogenomes, rendering them valuable for breeding.

A gene-rich mitochondrion with a unique ancestral protein transport system

Mitochondria originated from an ancient endosymbiosis involving an alphaproteobacterium.1,2,3 Over time, these organelles reduced their gene content massively, with most genes being transferred to the host nucleus before the last eukaryotic common ancestor (LECA).4 This process has yielded varying gene compositions in modern mitogenomes, including the complete loss of this organellar genome in some extreme cases.5,6,7,8,9,10,11,12,13,14 At the other end of the spectrum, jakobids harbor the most gene-rich mitogenomes, encoding 60-66 proteins.8 Here, we introduce the mitogenome of Mantamonas sphyraenae, a protist from the deep-branching CRuMs supergroup.15,16 Remarkably, it boasts the most gene-rich mitogenome outside of jakobids, by housing 91 genes, including 62 protein-coding ones. These include rare homologs of the four subunits of the bacterial-type cytochrome c maturation system I (CcmA, CcmB, CcmC, and CcmF) alongside a unique ribosomal protein S6. During the early evolution of mitochondria, gene transfer from the proto-mitochondrial endosymbiont to the nucleus became possible thanks to systems facilitating the transport of proteins synthesized in the host cytoplasm back to the mitochondrion. In addition to the universally found eukaryotic protein import systems, jakobid mitogenomes were reported to uniquely encode the SecY transmembrane protein of the Sec general secretory pathway, whose evolutionary origin was however unclear. The Mantamonas mitogenome not only encodes SecY but also SecA, SecE, and SecG, making it the sole eukaryote known to house a complete mitochondrial Sec translocation system. Furthermore, our phylogenetic and comparative genomic analyses provide compelling evidence for the alphaproteobacterial origin of this system, establishing its presence in LECA.

The first complete mitochondrial genome of Grossulariaceae: Molecular features, structure recombination, and genetic evolution

BACKGROUND

Mitochondria play crucial roles in the growth, development, and adaptation of plants. Blackcurrant (Ribes nigrum L.) stands out as a significant berry species due to its rich nutritional profile, medicinal properties, and health benefits. Despite its importance, the mitochondrial genome of blackcurrant remains unassembled.

RESULTS

This study presents the first assembly of the mitochondrial genome of R. nigrum in the Grossulariaceae family. The genome spans 450,227 base pairs (bp) and encompasses 39 protein-coding genes (PCGs), 19 transfer RNAs (tRNAs), and three ribosomal RNAs (rRNAs). Protein-coding regions constitute 8.88% of the entire genome. Additionally, we identified 180 simple sequence repeats, 12 tandem repeats, and 432 pairs of dispersed repeats. Notably, the dispersed sequence R1 (cotig3, 1,129 bp) mediated genome recombination, resulting in the formation of two major conformations, namely master and double circles. Furthermore, we identified 731 C-to-U RNA editing sites within the PCGs. Among these, cox1-2, nad1-2, and nad4L-2 were associated with the creation of start codons, whereas atp6-718 and rps10-391 were linked to termination codons. We also detected fourteen plastome fragments within the mitogenome, constituting 1.11% of the total length. Phylogenetic analysis suggests that R. nigrum might have undergone multiple genomic reorganization and/or gene transfer events, resulting in the loss of two PCGs (rps2 and rps11) during its evolutionary history.

CONCLUSIONS

This investigation unveils the molecular characteristics of the R. nigrum mitogenome, shedding light on its evolutionary trajectory and phylogenetic implications. Furthermore, it serves as a valuable reference for evolutionary research and germplasm identification within the genus.

Distinguishing Sanghuangporus from sanghuang-related fungi: a comparative and phylogenetic analysis based on mitogenomes

The Chinese medicinal fungi "Sanghuang" have been long recognized for their significant and valued medicinal properties, as documented in ancient medical literature. However, in traditional folk medicine, various macrofungi sharing similar appearance, habitat, and therapeutic effects with Sanghuang were erroneously used. These Sanghuang-like fungi mainly belong to the Porodaedalea, Phellinus, and Inonotus genera within the Hymenochaetaceae family. Despite the establishment of the Sanghuangporus genus and the identification of multiple species, the emerging taxonomic references based on morphological, ITS, and mycelial structural features have been inadequate to differentiate Sanghuangporus and Sanghuang-like fungi. To address this limitation, this study presents the first comparative and phylogenetic analysis of Sanghuang-related fungi based on mitogenomes. Our results show that Sanghuangporus species show marked convergence in mitochondrial genomic features and form a distinct monophyletic group based on phylogenetic analyses of five datasets. These results not only deepen our understanding of Sanghuang-like fungi but also offer novel insights into their mitochondrial composition and phylogeny, thereby providing new research tools for distinguishing members of the Sanghuangporus genus. KEY POINTS: • Sanghuangporus, Inonotus, and Porodaedalea are monophyly in sanghuang-like species. • Mitogenome-based analysis exhibits high resolution in sanghuang-like genus. • The mitogenomes provide strong evidence for reclassifying Phellinus gilvus S12 as Sanghuangporus vaninii.

Mechanisms and Future Research Perspectives on Mitochondrial Diseases Associated with Isoleucyl-tRNA Synthetase Gene Mutations

Aminoacyl-tRNA synthetases are essential enzymes for the accurate translation of genetic information. IARS1 and IARS2 are isoleucyl-tRNA synthetases functioning in the cytoplasm and mitochondria, respectively, with genetic mutations in these enzymes causing diverse clinical phenotypes in specific organs and tissues. Mutations in IARS1 and IARS2 have recently been linked to mitochondrial diseases. This review aims to explore the relationship between IARS1 and IARS2 and these diseases, providing a comprehensive overview of their association with mitochondrial diseases. Mutations in IARS1 cause weak calf syndrome in cattle and mitochondrial diseases in humans, leading to growth retardation and liver dysfunction. Mutations in IARS2 are associated with Leigh syndrome, craniosynostosis and abnormal genitalia syndrome. Future research is expected to involve genetic analysis of a larger number of patients, identifying new mutations in IARS1 and IARS2, and elucidating their impact on mitochondrial function. Additionally, genetically modified mice and the corresponding phenotypic analysis will serve as powerful tools for understanding the functions of these gene products and unraveling disease mechanisms. This will likely promote the development of new therapies and preventive measures.

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BACKGROUND

Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.

RESULTS

We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.

CONCLUSIONS

The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Friend turned foe: selfish behavior of a spontaneously arising mitochondrial deletion in an experimentally evolved Caenorhabditis elegans population

Selfish mitochondrial DNA (mtDNA) mutations are variants that can proliferate within cells and enjoy a replication or transmission bias without fitness benefits for the host. mtDNA deletions in Caenorhabditis elegans can reach high heteroplasmic frequencies despite significantly reducing fitness, illustrating how new mtDNA variants can give rise to genetic conflict between different levels of selection and between the nuclear and mitochondrial genomes. During a mutation accumulation experiment in C. elegans, a 1,034-bp deletion originated spontaneously and reached an 81.7% frequency within an experimental evolution line. This heteroplasmic mtDNA deletion, designated as meuDf1, eliminated portions of 2 protein-coding genes (coxIII and nd4) and tRNA-thr in entirety. mtDNA copy number in meuDf1 heteroplasmic individuals was 35% higher than in individuals with wild-type mitochondria. After backcrossing into a common genetic background, the meuDf1 mitotype was associated with reduction in several fitness traits and independent competition experiments found a 40% reduction in composite fitness. Experiments that relaxed individual selection by single individual bottlenecks demonstrated that the deletion-bearing mtDNA possessed a strong transmission bias, thereby qualifying it as a novel selfish mitotype.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

An overview on the antimalarial activity of 1,2,4-trioxanes, 1,2,4-trioxolanes and 1,2,4,5-tetraoxanes

The demand for novel, fast-acting, and effective antimalarial medications is increasing exponentially. Multidrug resistant forms of malarial parasites, which are rapidly spreading, pose a serious threat to global health. Drug resistance has been addressed using a variety of strategies, such as targeted therapies, the hybrid drug idea, the development of advanced analogues of pre-existing drugs, and the hybrid model of resistant strains control mechanisms. Additionally, the demand for discovering new potent drugs grows due to the prolonged life cycle of conventional therapy brought on by the emergence of resistant strains and ongoing changes in existing therapies. The 1,2,4-trioxane ring system in artemisinin (ART) is the most significant endoperoxide structural scaffold and is thought to be the key pharmacophoric moiety required for the pharmacodynamic potential of endoperoxide-based antimalarials. Several derivatives of artemisinin have also been found as potential treatments for multidrug-resistant strain in this area. Many 1,2,4-trioxanes, 1,2,4-trioxolanes, and 1,2,4,5-tetraoxanes derivatives have been synthesised as a result, and many of these have shown promise antimalarial activity both in vivo and in vitro against Plasmodium parasites. As a consequence, efforts to develop a functionally straight-forward, less expensive, and vastly more effective synthetic pathway to trioxanes continue. This study aims to give a thorough examination of the biological properties and mode of action of endoperoxide compounds derived from 1,2,4-trioxane-based functional scaffolds. The present system of 1,2,4-trioxane, 1,2,4-trioxolane, and 1,2,4,5-tetraoxane compounds and dimers with potentially antimalarial activity will be highlighted in this systematic review (January 1963-December 2022).

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Phenoptosis and the Various Types of Natural Selection

In the first description of evolution, the fundamental mechanism is the natural selection favoring the individuals best suited for survival and reproduction (selection at the individual level or classical Darwinian selection). However, this is a very reductive description of natural selection that does not consider or explain a long series of known phenomena, including those in which an individual sacrifices or jeopardizes his life on the basis of genetically determined mechanisms (i.e., phenoptosis). In fact, in addition to (i) selection at the individual level, it is essential to consider other types of natural selection such as those concerning: (ii) kin selection and some related forms of group selection; (iii) the interactions between the innumerable species that constitute a holobiont; (iv) the origin of the eukaryotic cell from prokaryotic organisms; (v) the origin of multicellular eukaryotic organisms from unicellular organisms; (vi) eusociality (e.g., in many species of ants, bees, termites); (vii) selection at the level of single genes, or groups of genes; (viii) the interactions between individuals (or more precisely their holobionts) of the innumerable species that make up an ecosystem. These forms of natural selection, which are all effects and not violations of the classical Darwinian selection, also show how concepts as life, species, individual, and phenoptosis are somewhat not entirely defined and somehow arbitrary. Furthermore, the idea of organisms selected on the basis of their survival and reproduction capabilities is intertwined with that of organisms also selected on the basis of their ability to cooperate and interact, even by losing their lives or their distinct identities.

Inflammatory signaling in NASH driven by hepatocyte mitochondrial dysfunctions

Liver steatosis, inflammation, and variable degrees of fibrosis are the pathological manifestations of nonalcoholic steatohepatitis (NASH), an aggressive presentation of the most prevalent chronic liver disease in the Western world known as nonalcoholic fatty liver (NAFL). Mitochondrial hepatocyte dysfunction is a primary event that triggers inflammation, affecting Kupffer and hepatic stellate cell behaviour. Here, we consider the role of impaired mitochondrial function caused by lipotoxicity during oxidative stress in hepatocytes. Dysfunction in oxidative phosphorylation and mitochondrial ROS production cause the release of damage-associated molecular patterns from dying hepatocytes, leading to activation of innate immunity and trans-differentiation of hepatic stellate cells, thereby driving fibrosis in NASH.

Mitochondrial Cristae Morphology Reflecting Metabolism, Superoxide Formation, Redox Homeostasis, and Pathology

Significance: Mitochondrial (mt) reticulum network in the cell possesses amazing ultramorphology of parallel lamellar cristae, formed by the invaginated inner mitochondrial membrane. Its non-invaginated part, the inner boundary membrane (IBM) forms a cylindrical sandwich with the outer mitochondrial membrane (OMM). Crista membranes (CMs) meet IBM at crista junctions (CJs) of mt cristae organizing system (MICOS) complexes connected to OMM sorting and assembly machinery (SAM). Cristae dimensions, shape, and CJs have characteristic patterns for different metabolic regimes, physiological and pathological situations. Recent Advances: Cristae-shaping proteins were characterized, namely rows of ATP-synthase dimers forming the crista lamella edges, MICOS subunits, optic atrophy 1 (OPA1) isoforms and mitochondrial genome maintenance 1 (MGM1) filaments, prohibitins, and others. Detailed cristae ultramorphology changes were imaged by focused-ion beam/scanning electron microscopy. Dynamics of crista lamellae and mobile CJs were demonstrated by nanoscopy in living cells. With tBID-induced apoptosis a single entirely fused cristae reticulum was observed in a mitochondrial spheroid. Critical Issues: The mobility and composition of MICOS, OPA1, and ATP-synthase dimeric rows regulated by post-translational modifications might be exclusively responsible for cristae morphology changes, but ion fluxes across CM and resulting osmotic forces might be also involved. Inevitably, cristae ultramorphology should reflect also mitochondrial redox homeostasis, but details are unknown. Disordered cristae typically reflect higher superoxide formation. Future Directions: To link redox homeostasis to cristae ultramorphology and define markers, recent progress will help in uncovering mechanisms involved in proton-coupled electron transfer via the respiratory chain and in regulation of cristae architecture, leading to structural determination of superoxide formation sites and cristae ultramorphology changes in diseases. Antioxid. Redox Signal. 39, 635-683.

Analysis of the Tomato mTERF Gene Family and Study of the Stress Resistance Function of SLmTERF-13

Mitochondrial transcription termination factor (mTERF) is a DNA-binding protein that is encoded by nuclear genes, ultimately functions in mitochondria and can affect gene expression. By combining with mitochondrial nucleic acids, mTERF regulates the replication, transcription and translation of mitochondrial genes and plays an important role in the response of plants to abiotic stress. However, there are few studies on mTERF genes in tomato, which limits the in-depth study and utilization of mTERF family genes in tomato stress resistance regulation. In this study, a total of 28 mTERF gene family members were obtained through genome-wide mining and identification of the tomato mTERF gene family. Bioinformatics analysis showed that all members of the family contained environmental stress or hormone response elements. Gene expression pattern analysis showed that the selected genes had different responses to drought, high salt and low temperature stress. Most of the genes played key roles under drought and salt stress, and the response patterns were more similar. The VIGS method was used to silence the SLmTERF13 gene, which was significantly upregulated under drought and salt stress, and it was found that the resistance ability of silenced plants was decreased under both kinds of stress, indicating that the SLmTERF13 gene was involved in the regulation of the tomato abiotic stress response. These results provide important insights for further evolutionary studies and contribute to a better understanding of the role of the mTERF genes in tomato growth and development and abiotic stress response, which will ultimately play a role in future studies of tomato gene function.

Organelle Genetics in Plants 2.0

Most of the DNA of eukaryotes is located in the nucleus [...].

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

Phenotypic and Genotypic Identification of Dermatophytes from Mexico and Central American Countries

Dermatophytes are fungi included in the genera Trichophyton, Microsporum, Epidermophyton, Nannizzia, Paraphyton, Lophophyton, and Arthroderma. Molecular techniques have contributed to faster and more precise identification, allowing significant advances in phylogenetic studies. This work aimed to identify clinical isolates of dermatophytes through phenotypic (macro- and micromorphology and conidia size) and genotypic methods (sequences of ITS regions, genes of β tubulin (BT2), and elongation factor α (Tef-1α)) and determine the phylogenetic relationships between isolates. Ninety-four dermatophyte isolates from Costa Rica, Guatemala, Honduras, Mexico, and the Dominican Republic were studied. The isolates presented macro- and micromorphology and conidia size described for the genera Trichophyton, Microsporum, and Epidermophyton. Genotypic analysis classified the isolates into the genera Trichophyton (63.8%), Nannizzia (25.5%), Arthroderma (9.6%), and Epidermophyton (1.1%). The most frequent species were T. rubrum (26 isolates, 27.6%), T. interdigitale (26 isolates, 27.6%), and N. incurvata (11 isolates, 11.7%), N. gypsea and A. otae (nine isolates, 9.6%), among others. The genotypic methods clarified the taxonomic status of closely related species. For instance, the ITS and BT2 markers of T. rubrum/T. violaceum did not differ but the Tef-1α gene did. On the other hand, the three markers differed in T. equinum/T. tonsurans. Therefore, the ITS, BT2, and Tef-1α genes are useful for typing in phylogenetic analyses of dermatophytes, with Tef-1α being the most informative locus. It should be noted that isolate MM-474 was identified as T. tonsurans when using ITS and Tef-1α, but when using BT2, it was identified as T. rubrum. On the other hand, no significant difference was found when comparing the methods for constructing phylogenies, as the topologies were similar.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

How mitochondria showcase evolutionary mechanisms and the importance of oxygen

Darwinian evolution can be simply stated: natural selection of inherited variations increasing differential reproduction. However, formulated thus, links with biochemistry, cell biology, ecology, and population dynamics remain unclear. To understand interactive contributions of chance and selection, higher levels of biological organization (e.g., endosymbiosis), complexities of competing selection forces, and emerging biological novelties (such as eukaryotes or meiotic sex), we must analyze actual examples. Focusing on mitochondria, I will illuminate how biology makes sense of life's evolution, and the concepts involved. First, looking at the bacterium - mitochondrion transition: merging with an archaeon, it lost its independence, but played a decisive role in eukaryogenesis, as an extremely efficient aerobic ATP generator and internal ROS source. Second, surveying later mitochondrion adaptations and diversifications illustrates concepts such as constructive neutral evolution, dynamic interactions between endosymbionts and hosts, the contingency of life histories, and metabolic reprogramming. Without oxygen, mitochondria disappear; with (intermittent) oxygen diversification occurs in highly complex ways, especially upon (temporary) phototrophic substrate supply. These expositions show the Darwinian model to be a highly fruitful paradigm.

Analysis of synonymous codon usage patterns in mitochondrial genomes of nine Amanita species

Introduction

Codon basis is a common and complex natural phenomenon observed in many kinds of organisms.

Methods

In the present study, we analyzed the base bias of 12 mitochondrial core protein-coding genes (PCGs) shared by nine Amanita species.

Results

The results showed that the codons of all Amanita species tended to end in A/T, demonstrating the preference of mitochondrial codons of Amanita species for a preference for this codon. In addition, we detected the correlation between codon base composition and the codon adaptation index (CAI), codon bias index (CBI), and frequency of optimal codons (FOP) indices, indicating the influence of base composition on codon bias. The average effective number of codons (ENC) of mitochondrial core PCGs of Amanita is 30.81, which is <35, demonstrating the strong codon preference of mitochondrial core PCGs of Amanita. The neutrality plot analysis and PR2-Bias plot analysis further demonstrated that natural selection plays an important role in Amanita codon bias. In addition, we obtained 5–10 optimal codons (ΔRSCU > 0.08 and RSCU > 1) in nine Amanita species, and GCA and AUU were the most widely used optimal codons. Based on the combined mitochondrial sequence and RSCU value, we deduced the genetic relationship between different Amanita species and found large variations between them.

Discussion

This study promoted the understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Supply-driven evolution: Mutation bias and trait-fitness distributions can drive macro-evolutionary dynamics

Many well-documented macro-evolutionary phenomena still challenge current evolutionary theory. Examples include long-term evolutionary trends, major transitions in evolution, conservation of certain biological features such as hox genes, and the episodic creation of new taxa. Here, we present a framework that may explain these phenomena. We do so by introducing a probabilistic relationship between trait value and reproductive fitness. This integration allows mutation bias to become a robust driver of long-term evolutionary trends against environmental bias, in a way that is consistent with all current evolutionary theories. In cases where mutation bias is strong, such as when detrimental mutations are more common than beneficial mutations, a regime called “supply-driven” evolution can arise. This regime can explain the irreversible persistence of higher structural hierarchies, which happens in the major transitions in evolution. We further generalize this result in the long-term dynamics of phenotype spaces. We show how mutations that open new phenotype spaces can become frozen in time. At the same time, new possibilities may be observed as a burst in the creation of new taxa.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown.

Methods

In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains.

Results

The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species.

Discussion

This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

Comparison of detection methods and genome quality when quantifying nuclear mitochondrial insertions in vertebrate genomes

The integration of mitochondrial genome fragments into the nuclear genome is well documented, and the transfer of these mitochondrial nuclear pseudogenes (numts) is thought to be an ongoing evolutionary process. With the increasing number of eukaryotic genomes available, genome-wide distributions of numts are often surveyed. However, inconsistencies in genome quality can reduce the accuracy of numt estimates, and methods used for identification can be complicated by the diverse sizes and ages of numts. Numts have been previously characterized in rodent genomes and it was postulated that they might be more prevalent in a group of voles with rapidly evolving karyotypes. Here, we examine 37 rodent genomes, and an additional 26 vertebrate genomes, while also considering numt detection methods. We identify numts using DNA:DNA and protein:translated-DNA similarity searches and compare numt distributions among rodent and vertebrate taxa to assess whether some groups are more susceptible to transfer. A combination of protein sequence comparisons (protein:translated-DNA) and BLASTN genomic DNA searches detect 50% more numts than genomic DNA:DNA searches alone. In addition, higher-quality RefSeq genomes produce lower estimates of numts than GenBank genomes, suggesting that lower quality genome assemblies can overestimate numts abundance. Phylogenetic analysis shows that mitochondrial transfers are not associated with karyotypic diversity among rodents. Surprisingly, we did not find a strong correlation between numt counts and genome size. Estimates using DNA: DNA analyses can underestimate the amount of mitochondrial DNA that is transferred to the nucleus.

The complete mitochondrial genome of Isochrysis galbana harbors a unique repeat structure and a specific trans-spliced cox1 gene

The haptophyte Isochrysis galbana is considered as a promising source for food supplements due to its rich fucoxanthin and polyunsaturated fatty acids content. Here, the I. galbana mitochondrial genome (mitogenome) was sequenced using a combination of Illumina and PacBio sequencing platforms. This 39,258 bp circular mitogenome has a total of 46 genes, including 20 protein-coding genes, 24 tRNA genes and two rRNA genes. A large block of repeats (~12.7 kb) was segregated in one region of the mitogenome, accounting for almost one third of the total size. A trans-spliced gene cox1 was first identified in I. galbana mitogenome and was verified by RNA-seq and DNA-seq data. The massive expansion of tandem repeat size and cis- to trans-splicing shift could be explained by the high mitogenome rearrangement rates in haptophytes. Strict SNP calling based on deep transcriptome sequencing data suggested the lack of RNA editing in both organelles in this species, consistent with previous studies in other algal lineages. To gain insight into haptophyte mitogenome evolution, a comparative analysis of mitogenomes within haptophytes and among eight main algal lineages was performed. A core gene set of 15 energy and metabolism genes is present in haptophyte mitogenomes, consisting of 1 cob, 3 cox, 7 nad, 2 atp and 2 ribosomal genes. Gene content and order was poorly conserved in this lineage. Haptophyte mitogenomes have lost many functional genes found in many other eukaryotes including rps/rpl, sdh, tat, secY genes, which make it contain the smallest gene set among all algal taxa. All these implied the rapid-evolving and more recently evolved mitogenomes of haptophytes compared to other algal lineages. The phylogenetic tree constructed by cox1 genes of 204 algal mitogenomes yielded well-resolved internal relationships, providing new evidence for red-lineages that contained plastids of red algal secondary endosymbiotic origin. This newly assembled mitogenome will add to our knowledge of general trends in algal mitogenome evolution within haptophytes and among different algal taxa.

Comparative analyses of Flammulina filiformis mitochondrial genomes reveal high length polymorphism in intergenic regions and multiple intron gain/loss in cox1

The golden-needle mushroom Flammulina filiformis is one of the bulk mushroom products in the world. This study obtained complete mitogenomes of 44 wild isolates collected from nine provinces and two artificially bred cultivars of F. filiformis, together with three Flammulina rossica isolates and one Flammulina fennae isolate for comparison. The mitogenome of F. filiformis ranged from 83,540 bp to 90,938 bp, consisting of 14 conserved protein-coding genes (PCGs), two rRNA genes, and 25 tRNA genes. To the best of our knowledge, it contained the highest proportion of intergenic regions compared to the other known Basidiomycota mitogenomes. Introns and intergenic regions were two major contributing factors to the total size of the F. filiformis mitogenome. The conserved PCG cox3 is located in an intron of another conserved PCG, nad5. This is a unique phenomenon in all known fungal mitogenomes. Gain/loss of introns was observed in cox1, nad5, and rnl. Length polymorphism was widely observed in intergenic regions. Accordingly, primers were designed as useful markers for rapid identification of F. filiformis isolates with differentiated mitogenomes. Our findings provide a basis for further studies related to variety identification and population genetics of this economically important mushroom.

Dynamic features of human mitochondrial DNA maintenance and transcription

Mitochondria are the primary sites for cellular energy production and are required for many essential cellular processes. Mitochondrial DNA (mtDNA) is a 16.6 kb circular DNA molecule that encodes only 13 gene products of the approximately 90 different proteins of the respiratory chain complexes and an estimated 1,200 mitochondrial proteins. MtDNA is, however, crucial for organismal development, normal function, and survival. MtDNA maintenance requires mitochondrially targeted nuclear DNA repair enzymes, a mtDNA replisome that is unique to mitochondria, and systems that control mitochondrial morphology and quality control. Here, we provide an overview of the current literature on mtDNA repair and transcription machineries and discuss how dynamic functional interactions between the components of these systems regulate mtDNA maintenance and transcription. A profound understanding of the molecular mechanisms that control mtDNA maintenance and transcription is important as loss of mtDNA integrity is implicated in normal process of aging, inflammation, and the etiology and pathogenesis of a number of diseases.

Variation in mitochondrial DNA affects locomotor activity and sleep in Drosophila melanogaster

Mitochondria are organelles that produce cellular energy in the form of ATP through oxidative phosphorylation, and this primary function is conserved among many taxa. Locomotion is a trait that is highly reliant on metabolic function and expected to be greatly affected by disruptions to mitochondrial performance. To this end, we aimed to examine how activity and sleep vary between Drosophila melanogaster strains with different geographic origins, how these patterns are affected by mitochondrial DNA (mtDNA) variation, and how breaking up co-evolved mito-nuclear gene combinations affect the studied activity traits. Our results demonstrate that Drosophila strains from different locations differ in sleep and activity, and that females are generally more active than males. By comparing activity and sleep of mtDNA variants introgressed onto a common nuclear background in cytoplasmic hybrid (cybrid) strains, we were able to quantify the among-line variance attributable to mitochondrial DNA, and we establish that mtDNA variation affects both activity and sleep, in a sex-specific manner. Altogether our study highlights the important role that mitochondrial genome variation plays on organismal physiology and behaviour.

Analysis of Mitochondria by Single-Organelle Resolution

Cellular organelles are highly specialized compartments with distinct functions. With the increasing resolution of detection methods, it is becoming clearer that same organelles may have different functions or properties not only within different cell populations of a tissue but also within the same cell. Dysfunction or altered function affects the organelle itself and may also lead to malignancies or undesirable cell death. To understand cellular function or dysfunction, it is therefore necessary to analyze cellular components at the single-organelle level. Here, we review the recent advances in analyzing cellular function at single-organelle resolution using high-parameter flow cytometry or multicolor confocal microscopy. We focus on the analysis of mitochondria, as they are organelles at the crossroads of various cellular signaling pathways and functions. However, most of the applied methods/technologies are transferable to any other organelle, such as the endoplasmic reticulum, lysosomes, or peroxisomes.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to 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 155 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 the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

Refining Mitochondrial Intron Classification With ERPIN: Identification Based on Conservation of Sequence Plus Secondary Structure Motifs

Mitochondrial genomes—in particular those of fungi—often encode genes with a large number of Group I and Group II introns that are conserved at both the sequence and the RNA structure level. They provide a rich resource for the investigation of intron and gene structure, self- and protein-guided splicing mechanisms, and intron evolution. Yet, the degree of sequence conservation of introns is limited, and the primary sequence differs considerably among the distinct intron sub-groups. It makes intron identification, classification, structural modeling, and the inference of gene models a most challenging and error-prone task—frequently passed on to an “expert” for manual intervention. To reduce the need for manual curation of intron structures and mitochondrial gene models, computational methods using ERPIN sequence profiles were initially developed in 2007. Here we present a refinement of search models and alignments using the now abundant publicly available fungal mtDNA sequences. In addition, we have tested in how far members of the originally proposed sub-groups are clearly distinguished and validated by our computational approach. We confirm clearly distinct mitochondrial Group I sub-groups IA1, IA3, IB3, IC1, IC2, and ID. Yet, IB1, IB2, and IB4 ERPIN models are overlapping substantially in predictions, and are therefore combined and reported as IB. We have further explored the conversion of our ERPIN profiles into covariance models (CM). Current limitations and prospects of the CM approach will be discussed.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the pathogenic fungus Ilyonectria destructans

Ilyonectria destructans is a pathogenic fungus causing root rot and other symptoms on trees and many crops. This paper analyses the mitochondrial genome of I. destructans and compares it with other published Nectriaceae mitogenomes. The I. destructans mitogenome appears as a circular DNA molecule of 42,895 bp and an overall GC content of 28.23%. It contains 28 protein-coding genes (15 core protein genes and 13 free-standing ORFs), two rRNAs and 27 tRNAs. The gene content and order were found to be conserved in the mitogenome of I. destructans and other Nectriaceae, although the genome size varies because of the variation in the number and length of intergenic regions and introns. For most core protein-coding genes in Nectriaceae species, Ka/Ks < 1 indicates purifying selection. Among some Nectriaceae representatives, only the rps3 gene was found under positive selection. Phylogenetic analyses based on nucleotide sequences of 15 protein-coding genes divided 45 Hypocreales species into six major clades matching the families Bionectriaceae, Cordycipitaceae, Clavicipitaceae, Ophiocordycipitaceae, Hypocreaceae and Nectriaceae. I. destructans appeared as a sister species to unidentified Ilyonectia sp., closely related to C. ilicicola, N. cinnabarina and a clad of ten Fusarium species and G. moniliformis. The complete mitogenome of I. destructans reported in the current paper will facilitate the study of epidemiology, biology, genetic diversity of the species and the evolution of family Nectriace and the Hypocreales order.

From Genome Variation to Molecular Mechanisms: What we Have Learned From Yeast Mitochondrial Genomes?

Analysis of genome variation provides insights into mechanisms in genome evolution. This is increasingly appreciated with the rapid growth of genomic data. Mitochondrial genomes (mitogenomes) are well known to vary substantially in many genomic aspects, such as genome size, sequence context, nucleotide base composition and substitution rate. Such substantial variation makes mitogenomes an excellent model system to study the mechanisms dictating mitogenome variation. Recent sequencing efforts have not only covered a rich number of yeast species but also generated genomes from abundant strains within the same species. The rich yeast genomic data have enabled detailed investigation from genome variation into molecular mechanisms in genome evolution. This mini-review highlights some recent progresses in yeast mitogenome studies.

Global Characterization of Fungal Mitogenomes: New Insights on Genomic Diversity and Dynamism of Coding Genes and Accessory Elements

Fungi comprise a great diversity of species with distinct ecological functions and lifestyles. Similar to other eukaryotes, fungi rely on interactions with prokaryotes and one of the most important symbiotic events was the acquisition of mitochondria. Mitochondria are organelles found in eukaryotic cells whose main function is to generate energy through aerobic respiration. Mitogenomes (mtDNAs) are double-stranded circular or linear DNA from mitochondria that may contain core genes and accessory elements that can be replicated, transcribed, and independently translated from the nuclear genome. Despite their importance, investigative studies on the diversity of fungal mitogenomes are scarce. Herein, we have evaluated 788 curated fungal mitogenomes available at NCBI database to assess discrepancies and similarities among them and to better understand the mechanisms involved in fungal mtDNAs variability. From a total of 12 fungal phyla, four do not have any representative with available mitogenomes, which highlights the underrepresentation of some groups in the current available data. We selected representative and non-redundant mitogenomes based on the threshold of 90% similarity, eliminating 81 mtDNAs. Comparative analyses revealed considerable size variability of mtDNAs with a difference of up to 260 kb in length. Furthermore, variation in mitogenome length and genomic composition are generally related to the number and length of accessory elements (introns, HEGs, and uORFs). We identified an overall average of 8.0 (0–39) introns, 8.0 (0–100) HEGs, and 8.2 (0–102) uORFs per genome, with high variation among phyla. Even though the length of the core protein-coding genes is considerably conserved, approximately 36.3% of the mitogenomes evaluated have at least one of the 14 core coding genes absent. Also, our results revealed that there is not even a single gene shared among all mitogenomes. Other unusual genes in mitogenomes were also detected in many mitogenomes, such as dpo and rpo, and displayed diverse evolutionary histories. Altogether, the results presented in this study suggest that fungal mitogenomes are diverse, contain accessory elements and are absent of a conserved gene that can be used for the taxonomic classification of the Kingdom Fungi.

Characterization of the mitochondrial genomes of three powdery mildew pathogens reveals remarkable variation in size and nucleotide composition

Powdery mildews comprise a large group of economically important phytopathogenic fungi. However, limited information exists on their mitochondrial genomes. Here, we assembled and compared the mitochondrial genomes of the powdery mildew pathogens Blumeria graminis f. sp. tritici, Erysiphe pisi, and Golovinomyces cichoracearum. Included in the comparative analysis was also the mitochondrial genome of Erysiphe necator that was previously analysed. The mitochondrial genomes of the four Erysiphales exhibit a similar gene content and organization but a large variation in size, with sizes ranging from 109800 bp in B. graminis f. sp. tritici to 332165 bp in G. cichoracearum, which is the largest mitochondrial genome of a fungal pathogen reported to date. Further comparative analysis revealed an unusual bimodal GC distribution in the mitochondrial genomes of B. graminis f. sp. tritici and G. cichoracearum that was not previously observed in fungi. The cytochrome b (cob) genes of E. necator, E. pisi, and G. cichoracearum were also exceptionally rich in introns, which in turn harboured rare open reading frames encoding reverse transcriptases that were likely acquired horizontally. Golovinomyces cichoracearum had also the longest cob gene (45 kb) among 703 fungal cob genes analysed. Collectively, these results provide novel insights into the organization of mitochondrial genomes of powdery mildew pathogens and represent valuable resources for population genetics and evolutionary studies.

Complete mitochondrial genome sequencing of Oxycarenus laetus (Hemiptera: Lygaeidae) from two geographically distinct regions of India

Oxycarenus laetus is a seed-sap sucking pest affecting a variety of crops, including cotton plants. Rising incidence and pesticide resistance by O. laetus have been reported from India and neighbouring countries. In this study, O. laetus samples were collected from Bhatinda and Coimbatore (India). Pure mtDNA was isolated and sequenced using Illumina MiSeq. Both the samples were found to be identical species (99.9%), and the complete genome was circular (15,672 bp), consisting of 13 PCGs, 2 rRNA, 23 tRNA genes, and a 962 bp control region. The mitogenome is 74.1% AT-rich, 0.11 AT, and - 0.19 GC skewed. All the genes had ATN as the start codon except cox1 (TTG), and an additional trnT was predicted. Nearly all tRNAs folded into the clover-leaf structure, except trnS1 and trnV. The intergenic space between trnH and nad4, considered as a synapomorphy of Lygaeoidea, was displaced. Two 5 bp motifs AATGA and ACCTA, two tandem repeats, and a few microsatellite sequences, were also found. The phylogenetic tree was constructed using 36 mitogenomes from 7 super-families of Hemiptera by employing rigorous bootstrapping and ML. Ours is the first study to sequence the complete mitogenome of O. laetus or any Oxycarenus species. The findings from this study would further help in the evolutionary studies of Lygaeidae.

Melatonin and Pathological Cell Interactions: Mitochondrial Glucose Processing in Cancer Cells

Melatonin is synthesized in the pineal gland at night. Since melatonin is produced in the mitochondria of all other cells in a non-circadian manner, the amount synthesized by the pineal gland is less than 5% of the total. Melatonin produced in mitochondria influences glucose metabolism in all cells. Many pathological cells adopt aerobic glycolysis (Warburg effect) in which pyruvate is excluded from the mitochondria and remains in the cytosol where it is metabolized to lactate. The entrance of pyruvate into the mitochondria of healthy cells allows it to be irreversibly decarboxylated by pyruvate dehydrogenase (PDH) to acetyl coenzyme A (acetyl-CoA). The exclusion of pyruvate from the mitochondria in pathological cells prevents the generation of acetyl-CoA from pyruvate. This is relevant to mitochondrial melatonin production, as acetyl-CoA is a required co-substrate/co-factor for melatonin synthesis. When PDH is inhibited during aerobic glycolysis or during intracellular hypoxia, the deficiency of acetyl-CoA likely prevents mitochondrial melatonin synthesis. When cells experiencing aerobic glycolysis or hypoxia with a diminished level of acetyl-CoA are supplemented with melatonin or receive it from another endogenous source (pineal-derived), pathological cells convert to a more normal phenotype and support the transport of pyruvate into the mitochondria, thereby re-establishing a healthier mitochondrial metabolic physiology.

The Gene Rearrangement, Loss, Transfer, and Deep Intronic Variation in Mitochondrial Genomes of Conidiobolus

The genus Conidiobolus s.s. was newly delimited from Conidiobolus s.l. In order to gain insight into its mitochondrial genetic background, this study sequenced six mitochondrial genomes of the genus Conidiobolus s.s. These mitogenomes were all composed of circular DNA molecules, ranging from 29,253 to 48,417 bp in size and from 26.61 to 27.90% in GC content. The order and direction for 14 core protein-coding genes (PCGs) were identical, except for the atp8 gene lost in Conidiobolus chlamydosporus, Conidiobolus polyspermus, and Conidiobolus polytocus, and rearranged in the other Conidiobolus s.s. species. Besides, the atp8 gene split the cox1 gene in Conidiobolus taihushanensis. Phylogenomic analysis based on the 14 core PCGs confirmed that all Conidiobolus s.s. species formed a monophyly in the Entomophthoromycotina lineage. The number and length of introns were the main factors contributing to mitogenomic size, and deep variations and potential transfer were detected in introns. In addition, gene transfer occurred between the mitochondrial and nuclear genomes. This study promoted the understanding of the evolution and phylogeny of the Conidiobolus s.s. genus.

The Elusive Mitochondrial Genomes of Apicomplexa: Where Are We Now?

Mitochondria are vital organelles of eukaryotic cells, participating in key metabolic pathways such as cellular respiration, thermogenesis, maintenance of cellular redox potential, calcium homeostasis, cell signaling, and cell death. The phylum Apicomplexa is entirely composed of obligate intracellular parasites, causing a plethora of severe diseases in humans, wild and domestic animals. These pathogens include the causative agents of malaria, cryptosporidiosis, neosporosis, East Coast fever and toxoplasmosis, among others. The mitochondria in Apicomplexa has been put forward as a promising source of undiscovered drug targets, and it has been validated as the target of atovaquone, a drug currently used in the clinic to counter malaria. Apicomplexans present a single tubular mitochondria that varies widely both in structure and in genomic content across the phylum. The organelle is characterized by massive gene migrations to the nucleus, sequence rearrangements and drastic functional reductions in some species. Recent third generation sequencing studies have reignited an interest for elucidating the extensive diversity displayed by the mitochondrial genomes of apicomplexans and their intriguing genomic features. The underlying mechanisms of gene transcription and translation are also ill-understood. In this review, we present the state of the art on mitochondrial genome structure, composition and organization in the apicomplexan phylum revisiting topological and biochemical information gathered through classical techniques. We contextualize this in light of the genomic insight gained by second and, more recently, third generation sequencing technologies. We discuss the mitochondrial genomic and mechanistic features found in evolutionarily related alveolates, and discuss the common and distinct origins of the apicomplexan mitochondria peculiarities.

Role of Mitochondria in Radiation Responses: Epigenetic, Metabolic, and Signaling Impacts

Until recently, radiation effects have been considered to be mainly due to nuclear DNA damage and their management by repair mechanisms. However, molecular biology studies reveal that the outcomes of exposures to ionizing radiation (IR) highly depend on activation and regulation through other molecular components of organelles that determine cell survival and proliferation capacities. As typical epigenetic-regulated organelles and central power stations of cells, mitochondria play an important pivotal role in those responses. They direct cellular metabolism, energy supply and homeostasis as well as radiation-induced signaling, cell death, and immunological responses. This review is focused on how energy, dose and quality of IR affect mitochondria-dependent epigenetic and functional control at the cellular and tissue level. Low-dose radiation effects on mitochondria appear to be associated with epigenetic and non-targeted effects involved in genomic instability and adaptive responses, whereas high-dose radiation effects (>1 Gy) concern therapeutic effects of radiation and long-term outcomes involving mitochondria-mediated innate and adaptive immune responses. Both effects depend on radiation quality. For example, the increased efficacy of high linear energy transfer particle radiotherapy, e.g., C-ion radiotherapy, relies on the reduction of anastasis, enhanced mitochondria-mediated apoptosis and immunogenic (antitumor) responses.

Mitochondrial phylogenetic and diversity analysis in Azi-Kheli buffalo

NOVELTY STATEMENT

The present study was conducted for the first time in Pakistan to investigate Cytochrome C Oxidase Subunit 1 (CO1) gene and full-length Displacement Loop (D-loop) region of mitochondrial DNA in Azi-Kheli buffalo breed native to northern hilly areas of Khyber Pakhtunkhwa Province of Pakistan. The present study was designed to investigate phylogeny and diversity in Azi-Kheli buffalo, through two mitochondrial DNA regions, i.e., Cytochrome C Oxidase Subunit-I (CO1) and Displacement Loop (D-loop) region. Thirty (30) blood samples were taken from Azi-Kheli pure breed animals from original breeding tract, i.e., Khwazakhela, Swat. Polymerase chain reactions using gene-specific primers were carried out for amplifying 709-bp region of CO1 gene and 1159-bp region of D-Loop for identification, phylogeny, and diversity in Azi-Kheli buffalo, respectively. The sequences of CO1 gene revealed four (04) haplotypes, whereas D-loop sequences revealed five (05) haplotypes. Mean interspecific diversity with related species was 2.56%, and mean intraspecific diversity within Azi-Kheli buffalo was 0.25%, estimated via Kimura-2 parameter. Phylogenetic tree (maximum likelihood) revealed clustering of Azi-Kheli haplotypes with river buffalo and is distinct from swamp buffalo and other related species of genus Bubalus. Mean haplotype and nucleotide diversity of D-loop were Hd = 0.9601 ± SD = 0.096 and π = 0.01208 ± SD = 0.00182, respectively. Phylogenetic tree (neighbor-joining) revealed two main clades, i.e., river buffalo and swamp buffalo clade. The haplotypes of Azi-Kheli clustered with haplotypes of different river buffalo breeds at different positions. The current study suggests that Azi-Kheli has common origin with other river buffalo breeds; hence, it is river buffalo which harbors high genetic diversity.

Skeletal Muscle Mitochondrial Physiology in Children With Cerebral Palsy: Considerations for Healthy Aging

Skeletal muscle contractile proteins require a constant supply of energy to produce force needed for movement. Energy (ATP) is primarily produced by mitochondrial organelles, located within and around muscle fibers, by oxidative phosphorylation that couples electron flux through the electron transport chain to create a proton gradient across the inner mitochondrial membrane that is in turn used by the ATP synthase. Mitochondrial networks increase in size by biogenesis to increase mitochondrial abundance and activity in response to endurance exercise, while their function and content reduce with constant inactivity, such as during muscle atrophy. During healthy aging, there is an overall decline in mitochondrial activity and abundance, increase in mitochondrial DNA mutations, potential increase in oxidative stress, and reduction in overall muscular capacity. Many of these alterations can be attenuated by consistent endurance exercise. Children with cerebral palsy (CP) have significantly increased energetics of movement, reduced endurance capacity, and increased perceived effort. Recent work in leg muscles in ambulatory children with CP show a marked reduction in mitochondrial function. Arm muscles show that mitochondrial protein content and mitochondria DNA copy number are lower, suggesting a reduction in mitochondrial abundance, along with a reduction in markers for mitochondrial biogenesis. Gene expression networks are reduced for glycolytic and mitochondrial pathways and share similarities with gene networks with aging and chronic inactivity. Given the importance of mitochondria for energy production and changes with aging, future work needs to assess changes in mitochondria across the lifespan in people with CP and the effect of exercise on promoting metabolic health.

Photosynthetic Systems Suggest an Evolutionary Pathway to Diderms

Bacteria are divided primarily into monoderms (with one cell membrane, and usually Gram-positive, due to a thick peptidoglycan layer) and diderms (with two cell membranes, and mostly Gram-negative, due to a thin peptidoglycan layer sandwiched between the two membranes). Photosynthetic species are spread among the taxonomic groups, some having type I reaction centers (RCI in monoderm phylum Firmicutes; and diderm phyla Acidobacteria and Chlorobi), others with type II reaction centers (RCII in monoderm phylum Chloroflexi; and diderm taxa Gemmatimonadetes, and alpha-, beta-, and gamma-Proteobacteria), and some containing both (RCI and RCII, only in diderm phylum Cyanobacteria). In most bacterial phylograms, photosystem types and diderm taxa are polyphyletic. A more parsimonious arrangement, which is supported by photosystem evolution, as well as additional sets of molecular characters, suggests that endosymbiotic events resulted in the formation of the diderms. In the model presented, monoderms readily form a monophyletic group, while diderms are produced by at least two endosymbiotic events, followed by additional evolutionary changes.

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.

The Coevolution of Fungal Mitochondrial Introns and Their Homing Endonucleases (GIY-YIG and LAGLIDADG)

Fungal mitochondrial (mt) genomes exhibit great diversity in size which is partially attributed to their variable intergenic regions and most importantly to the inclusion of introns within their genes. These introns belong to group I or II, and both of them are self-splicing. The majority of them carry genes encoding homing endonucleases, either LAGLIDADG or GIY-YIG. In this study, it was found that these intronic homing endonucleases genes (HEGs) may originate from mt free-standing open reading frames which can be found nowadays in species belonging to Early Diverging Fungi as “living fossils.” A total of 487 introns carrying HEGs which were located in the publicly available mt genomes of representative species belonging to orders from all fungal phyla was analyzed. Their distribution in the mt genes, their insertion target sequence, and the phylogenetic analyses of the HEGs showed that these introns along with their HEGs form a composite structure in which both selfish elements coevolved. The invasion of the ancestral free-standing HEGs in the introns occurred through a perpetual mechanism, called in this study as “aenaon” hypothesis. It is based on recombination, transpositions, and horizontal gene transfer events throughout evolution. HEGs phylogenetically clustered primarily according to their intron hosts and secondarily to the mt genes carrying the introns and their HEGs. The evolutionary models created revealed an “intron-early” evolution which was enriched by “intron-late” events through many different independent recombinational events which resulted from both vertical and horizontal gene transfers.

Pervasive transcription of the mitochondrial genome in Candida albicans is revealed in mutants lacking the mtEXO RNase complex

The mitochondrial genome of the pathogenic yeast Candida albicans displays a typical organization of several (eight) primary transcription units separated by noncoding regions. Presence of genes encoding Complex I subunits and the stability of its mtDNA sequence make it an attractive model to study organellar genome expression using transcriptomic approaches. The main activity responsible for RNA degradation in mitochondria is a two-component complex (mtEXO) consisting of a 3ʹ-5ʹ exoribonuclease, in yeasts encoded by the DSS1 gene, and a conserved Suv3p helicase. In C. albicans, deletion of either DSS1 or SUV3 gene results in multiple defects in mitochondrial genome expression leading to the loss of respiratory competence. Transcriptomic analysis reveals pervasive transcription in mutants lacking the mtEXO activity, with evidence of the entire genome being transcribed, whereas in wild-type strains no RNAs corresponding to a significant fraction of the noncoding genome can be detected. Antisense (‘mirror’) transcripts, absent from normal mitochondria are also prominent in the mutants. The expression of multiple mature transcripts, particularly those translated from bicistronic mRNAs, as well as those that contain introns is affected in the mutants, resulting in a decreased level of proteins and reduced respiratory complex activity. The phenotype is most severe in the case of Complex IV, where a decrease of mature COX1 mRNA level to ~5% results in a complete loss of activity. These results show that RNA degradation by mtEXO is essential for shaping the mitochondrial transcriptome and is required to maintain the functional demarcation between transcription units and non-coding genome segments.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.