Intronic GIY-YIG endonuclease gene in the mitochondrial genome of Podospora curvicolla: evidence for mobility

Comparative analysis of mitochondrial genomes in Ceratocystis fimbriata complex across diverse hosts

The decline ofAcacia mangiumWilld. in Malaysia, especially in Sabah since 2010, is primarily due to Ceratocystiswilt and canker disease (CWCD) caused by theCeratocystis fimbriataEllis & Halst. complex. This study was aimed to investigate the mitochondrial genome architecture of two differentC. fimbriatacomplex isolates from Malaysia: one fromA. mangiumin Pahang (FRIM1162) and another fromEucalyptus pellitain Sarawak (FRIM1441). This research employed Next-Generation Sequencing (NGS) to contrast genomes from diverse hosts with nine additional mitochondrial sequences, identifying significant genetic diversity and mutational hotspots in the mitochondrial genome alignment. The mitochondrial genome-based phylogenetic analysis revealed a significant genetic relationship between the studied isolates and theC. fimbriatacomplex in the South American Subclade, indicating that theC. fimbriatacomplex discovered in Malaysia isC. manginecans. The comparative mitochondrial genome demonstrates the adaptability of the complex due to mobile genetic components and genomic rearrangements in the studiedfungal isolates. This research enhances our knowledge of the genetic diversity and evolutionary patterns within theC. fimbriatacomplex, aiding in a deeper understanding of fungal disease development and host adaption processes. The acquired insights are crucial for creating specific management strategies for CWCD, improving the overall understanding of fungal disease evolution and control.

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.

Reevaluation of the Toxoplasma gondii and Neospora caninum genomes reveals misassembly, karyotype differences, and chromosomal rearrangements

Neosporacaninum primarily infects cattle, causing abortions, with an estimated impact of a billion dollars on the worldwide economy annually. However, the study of its biology has been unheeded by the established paradigm that it is virtually identical to its close relative, the widely studied human pathogen Toxoplasma gondii. By revisiting the genome sequence, assembly, and annotation using third-generation sequencing technologies, here we show that the N. caninum genome was originally incorrectly assembled under the presumption of synteny with T. gondii. We show that major chromosomal rearrangements have occurred between these species. Importantly, we show that chromosomes originally named Chr VIIb and VIII are indeed fused, reducing the karyotype of both N. caninum and T. gondii to 13 chromosomes. We reannotate the N. caninum genome, revealing more than 500 new genes. We sequence and annotate the nonphotosynthetic plastid and mitochondrial genomes and show that although apicoplast genomes are virtually identical, high levels of gene fragmentation and reshuffling exist between species and strains. Our results correct assembly artifacts that are currently widely distributed in the genome database of N. caninum and T. gondii and, more importantly, highlight the mitochondria as a previously oversighted source of variability and pave the way for a change in the paradigm of synteny, encouraging rethinking the genome as basis of the comparative unique biology of these pathogens.

Mitochondrial evolution in the entomopathogenic fungal genus Beauveria

Species in the fungal genus Beauveria are pathogens of invertebrates and have been commonly used as the active agent in biopesticides. After many decades with few species described, recent molecular approaches to classification have led to over 25 species now delimited. Little attention has been given to the mitochondrial genomes of Beauveria but better understanding may led to insights into the nature of species and evolution in this important genus. In this study, we sequenced the mitochondrial genomes of four new strains belonging to Beauveria bassiana, Beauveria caledonica and Beauveria malawiensis, and compared them to existing mitochondrial sequences of related fungi. The mitochondrial genomes of Beauveria ranged widely from 28,806 to 44,135 base pairs, with intron insertions accounting for most size variation and up to 39% (B. malawiensis) of the mitochondrial length due to introns in genes. Gene order of the common mitochondrial genes did not vary among the Beauveria sequences, but variation was observed in the number of transfer ribonucleic acid genes. Although phylogenetic analysis using whole mitochondrial genomes showed, unsurprisingly, that B. bassiana isolates were the most closely related to each other, mitochondrial codon usage suggested that some B. bassiana isolates were more similar to B. malawiensis and B. caledonica than the other B. bassiana isolates analyzed.

rps3 as a Candidate Mitochondrial Gene for the Molecular Identification of Species from the Colletotrichum acutatum Species Complex

Colletotrichum species form one of the most economically significant groups of pathogenic fungi and lead to significant losses in the production of major crops-in particular, fruits, vegetables, ornamental plants, shrubs, and trees. Members of the genus Colletotrichum cause anthracnose disease in many plants. Due to their considerable variation, these fungi have been widely investigated in genetic studies as model organisms. Here, we report the complete mitochondrial genome sequences of four Colletotrichum species (C. fioriniae, C. lupini, C. salicis, and C. tamarilloi). The reported circular mitogenomes range from 30,020 (C. fioriniae) to 36,554 bp (C. lupini) in size and have identical sets of genes, including 15 protein-coding genes, two ribosomal RNA genes, and 29 tRNA genes. All four mitogenomes are characterized by a rather poor repetitive sequence content with only forward repeat representatives and a low number of microsatellites. The topology of the phylogenetic tree reflects the systematic positions of the studied species, with representatives of each Colletotrichum species complex gathered in one clade. A comparative analysis reveals consistency in the gene composition and order of Colletotrichum mitogenomes, although some highly divergent regions are also identified, like the rps3 gene which appears as a source of potential diagnostic markers for all studied Colletotrichum species.

Annotation and analysis of the mitochondrial genome of Coniothyrium glycines, causal agent of red leaf blotch of soybean, reveals an abundance of homing endonucleases

Coniothyrium glycines, the causal agent of soybean red leaf blotch, is a USDA APHIS-listed Plant Pathogen Select Agent and potential threat to US agriculture. Sequencing of the C. glycines mt genome revealed a circular 98,533-bp molecule with a mean GC content of 29.01%. It contains twelve of the mitochondrial genes typically involved in oxidative phosphorylation (atp6, cob, cox1-3, nad1-6, and nad4L), one for a ribosomal protein (rps3), four for hypothetical proteins, one for each of the small and large subunit ribosomal RNAs (rns and rnl) and a set of 30 tRNAs. Genes were encoded on both DNA strands with cox1 and cox2 occurring as adjacent genes having no intergenic spacers. Likewise, nad2 and nad3 are adjacent with no intergenic spacers and nad5 is immediately followed by nad4L with an overlap of one base. Thirty-two introns, comprising 54.1% of the total mt genome, were identified within eight protein-coding genes and the rnl. Eighteen of the introns contained putative intronic ORFs with either LAGLIDADG or GIY-YIG homing endonuclease motifs, and an additional eleven introns showed evidence of truncated or degenerate endonuclease motifs. One intron possessed a degenerate N-acetyl-transferase domain. C. glycines shares some conservation of gene order with other members of the Pleosporales, most notably nad6-rnl-atp6 and associated conserved tRNA clusters. Phylogenetic analysis of the twelve shared protein coding genes agrees with commonly accepted fungal taxonomy. C. glycines represents the second largest mt genome from a member of the Pleosporales sequenced to date. This research provides the first genomic information on C. glycines, which may provide targets for rapid diagnostic assays and population studies.

Whole genome sequencing of Rhodotorula mucilaginosa isolated from the chewing stick (Distemonanthus benthamianus): insights into Rhodotorula phylogeny, mitogenome dynamics and carotenoid biosynthesis

In industry, the yeast Rhodotorula mucilaginosa is commonly used for the production of carotenoids. The production of carotenoids is important because they are used as natural colorants in food and some carotenoids are precursors of retinol (vitamin A). However, the identification and molecular characterization of the carotenoid pathway/s in species belonging to the genus Rhodotorula is scarce due to the lack of genomic information thus potentially impeding effective metabolic engineering of these yeast strains for improved carotenoid production. In this study, we report the isolation, identification, characterization and the whole nuclear genome and mitogenome sequence of the endophyte R. mucilaginosa RIT389 isolated from Distemonanthus benthamianus, a plant known for its anti-fungal and antibacterial properties and commonly used as chewing sticks. The assembled genome of R. mucilaginosa RIT389 is 19 Mbp in length with an estimated genomic heterozygosity of 9.29%. Whole genome phylogeny supports the species designation of strain RIT389 within the genus in addition to supporting the monophyly of the currently sequenced Rhodotorula species. Further, we report for the first time, the recovery of the complete mitochondrial genome of R. mucilaginosa using the genome skimming approach. The assembled mitogenome is at least 7,000 bases larger than that of Rhodotorula taiwanensis which is largely attributed to the presence of large intronic regions containing open reading frames coding for homing endonuclease from the LAGLIDADG and GIY-YIG families. Furthermore, genomic regions containing the key genes for carotenoid production were identified in R. mucilaginosa RIT389, revealing differences in gene synteny that may play a role in the regulation of the biotechnologically important carotenoid synthesis pathways in yeasts.

The Agaricus bisporus cox1 gene: the longest mitochondrial gene and the largest reservoir of mitochondrial group i introns

In eukaryotes, introns are located in nuclear and organelle genes from several kingdoms. Large introns (up to 5 kbp) are frequent in mitochondrial genomes of plant and fungi but scarce in Metazoa, even if these organisms are grouped with fungi among the Opisthokonts. Mitochondrial introns are classified in two groups (I and II) according to their RNA secondary structure involved in the intron self-splicing mechanism. Most of these mitochondrial group I introns carry a “Homing Endonuclease Gene” (heg) encoding a DNA endonuclease acting in transfer and site-specific integration (“homing”) and allowing intron spreading and gain after lateral transfer even between species from different kingdoms. Opposed to this gain mechanism, is another which implies that introns, which would have been abundant in the ancestral genes, would mainly evolve by loss. The importance of both mechanisms (loss and gain) is matter of debate. Here we report the sequence of the cox1 gene of the button mushroom Agaricus bisporus, the most widely cultivated mushroom in the world. This gene is both the longest mitochondrial gene (29,902 nt) and the largest group I intron reservoir reported to date with 18 group I and 1 group II. An exhaustive analysis of the group I introns available in cox1 genes shows that they are mobile genetic elements whose numerous events of loss and gain by lateral transfer combine to explain their wide and patchy distribution extending over several kingdoms. An overview of intron distribution, together with the high frequency of eroded heg, suggests that they are evolving towards loss. In this landscape of eroded and lost intron sequences, the A. bisporus cox1 gene exhibits a peculiar dynamics of intron keeping and catching, leading to the largest collection of mitochondrial group I introns reported to date in a Eukaryote.

Hpy188I-DNA pre- and post-cleavage complexes--snapshots of the GIY-YIG nuclease mediated catalysis

The GIY-YIG nuclease domain is present in all kingdoms of life and has diverse functions. It is found in the eukaryotic flap endonuclease and Holliday junction resolvase Slx1–Slx4, the prokaryotic nucleotide excision repair proteins UvrC and Cho, and in proteins of ‘selfish’ genetic elements. Here we present the structures of the ternary pre- and post-cleavage complexes of the type II GIY-YIG restriction endonuclease Hpy188I with DNA and a surrogate or catalytic metal ion, respectively. Our structures suggest that GIY-YIG nucleases catalyze DNA hydrolysis by a single substitution reaction. They are consistent with a previous proposal that a tyrosine residue (which we expect to occur in its phenolate form) acts as a general base for the attacking water molecule. In contrast to the earlier proposal, our data identify the general base with the GIY and not the YIG tyrosine. A conserved glutamate residue (Glu149 provided in trans in Hpy188I) anchors a single metal cation in the active site. This metal ion contacts the phosphate proS oxygen atom and the leaving group 3′-oxygen atom, presumably to facilitate its departure. Taken together, our data reveal striking analogy in the absence of homology between GIY-YIG and ββα-Me nucleases.

Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs

Homing endonucleases typically contain one of four conserved catalytic motifs, and other elements that confer tight DNA binding. I-CreII, which catalyzes homing of the Cr.psbA4 intron, is unusual in containing two potential catalytic motifs, H-N-H and GIY-YIG. Previously, we showed that cleavage by I-CreII leaves ends (2-nt 3' overhangs) that are characteristic of GIY-YIG endonucleases, yet it has a relaxed metal requirement like H-N-H enzymes. Here we show that I-CreII can bind DNA without an added metal ion, and that it binds as a monomer, akin to GIY-YIG enzymes. Moreover, cleavage of supercoiled DNA, and estimates of strand-specific cleavage rates, suggest that I-CreII uses a sequential cleavage mechanism. Alanine substitution of a number of residues in the GIY-YIG motif, however, did not block cleavage activity, although DNA binding was substantially reduced in several variants. Substitution of conserved histidines in the H-N-H motif resulted in variants that did not promote DNA cleavage, but retained high-affinity DNA binding-thus identifying it as the catalytic motif. Unlike the non-specific H-N-H colicins, however; substitution of the conserved asparagine substantially reduced DNA binding (though not the ability to promote cleavage). These results indicate that, in I-CreII, two catalytic motifs have evolved to play important roles in specific DNA binding. The data also indicate that only the H-N-H motif has retained catalytic ability.

Insertion of a homing endonuclease creates a genes-in-pieces ribonucleotide reductase that retains function

In bacterial and phage genomes, coding regions are sometimes interrupted by self-splicing introns or inteins, which can encode mobility-promoting homing endonucleases. Homing endonuclease genes are also found free-standing (not intron- or intein-encoded) in phage genomes where they are inserted in intergenic regions. One example is the HNH family endonuclease, mobE, inserted between the large (nrdA) and small (nrdB) subunit genes of aerobic ribonucleotide reductase (RNR) of T-even phages T4, RB2, RB3, RB15, and LZ7. Here, we describe an insertion of mobE into the nrdA gene of Aeromonas hydrophila phage Aeh1. The insertion creates a unique genes-in-pieces arrangement, where nrdA is split into two independent genes, nrdA-a and nrdA-b, each encoding cysteine residues that correspond to the active-site residues of uninterrupted NrdA proteins. Remarkably, the mobE insertion does not inactivate NrdA function, although the insertion is not a self-splicing intron or intein. We copurified the NrdA-a, NrdA-b, and NrdB proteins as complex from Aeh1-infected cells and also showed that a reconstituted complex has RNR activity. Class I RNR activity in phage Aeh1 is thus assembled from separate proteins that interact to form a composite active site, demonstrating that the mobE insertion is phenotypically neutral in that its presence as an intervening sequence does not disrupt the function of the surrounding gene.

Interpretation of mtDNA RFLP variability among Aspergillus tubingensis isolates

Aspergillus tubingensis isolates collected from distant geographic areas were earlier classified into six groups on the basis of the mtDNA RFLP variability they exhibited (mtDNA types 2a-2f). In the present work, we investigated the reason for the intraspecific mtDNA variability and we describe here how this fungus, with a relatively small mitochondrial genome, can display intraspecific polymorphism due to intron acquisition and also sporadic point mutations affecting the recognition motifs of the restriction enzymes employed in the RFLP analysis. Three different LAGLI-DADG type group I introns were identified in the cox1 gene amongst the six mtDNA RFLP types. MtDNAs of types 2b and 2d contain all of the three introns, mtDNA of type 2f carries only one, and the other mtDNA types contain two introns each. Comparative analysis showed that the first and second introns of mtDNAs of types 2b and 2d are well distributed among fungi, indicating their active horizontal transfer capacity. The third intron occurs rarely among fungi and is restricted to a limited number of fungal species, namely to A. tubingensis and the yeast Candida stellata. It is interesting that this intron is present in a small mitochondrial genome such as that of A. tubingensis and, considering its rarity, its presence amongst black Aspergillus isolates is recommended to be considered as a tool to establish taxonomical unit(s) or to track down evolutionary divergence of closely related taxonomical units.

Expression, purification, and biochemical characterization of the intron-encoded endonuclease, I-CreII

The ORF of the Cr.psbA4 intron of Chlamydomonas reinhardtii mediates efficient intron homing, and contains an H-N-H and possibly a GIY-YIG motif. The ORF was over-expressed in Escherichia coli without non-native amino acids, but was mostly insoluble. However, co-over-expression of E. coli chaperonins GroEL/GroES solubilized approximately 50% of the protein, which was purified by ion-exchange and heparin-affinity chromatography. Biochemical characterization showed that the protein is a double-strand-specific endonuclease that cleaves fused psbA exon 4-exon 5 DNA, and was named I-CreII. I-CreII has a relatively relaxed divalent metal ion requirement (Mg(2+), Mn(2+), Ca(2+), and Fe(2+) supported cleavage), is insensitive to salt <350 mM, and is stabilized by DNA. Cleavage of target DNA occurs close (4 nt on the top strand) to the intron-insertion site, and leaves 2-nt 3'-OH overhangs, similar to GIY-YIG endonucleases. The boundaries of the recognition sequence span approximately 30 bp, and encompass the cleavage and intron-insertion sites. Cleavage of heterologous psbA DNAs indicates the enzyme can tolerate multiple, but not all, substitutions in the recognition site. This work will facilitate further study of this novel endonuclease, which may also find use in site-specific manipulation of chloroplast DNA.

Optional mitochondrial introns and evidence for a homing-endonuclease gene in the mtDNA rnl gene in Ophiostoma ulmi s. lat

Strains of Ophiostoma ulmi, O. novo-ulmi subsp. americana, O. novo-ulmi subsp. novo-ulmi and O. himal-ulmi were examined for optional introns/insertions within the following mitochondrial genes: small subunit RNA gene (rns), large ribosomal subunit gene (rnl) and the cytochrome oxidase subunit I gene (coxI). Insertions were noted in the rns and coxI genes in strains of O. ulmi, the less aggressive species, but absent in strains of the more aggressive O. novo-ulmi subsp. americana. Strains of all species examined had a group I intron present in the U11 region of the mitochondrial-rnl gene. In all but two strains of O. novo-ulmi subsp. americana, this rnl-U11 intron was about 1.5 kb in length whereas a 2.6 kb version of this element was present in all strains representing O. ulmi, O. novo-ulmi subsp. novo-ulmi, and Ophiostoma himal-ulmi. Irrespective of size, this intron based on RNA folds is a class IA1 group I intron and it encodes a putative ORF for the rps3 ribosomal protein. The size variation of the rnl-U11 intron was examined in detail for two strains of O. novo-ulmi subsp. americana and sequence data suggests the presence of a complex ORF within the 2.6 kb version of this intron; here a homing endonuclease-like gene has been inserted in frame and fused to the carboxyl-terminus of the putative rps3 coding region. The mitochondrial optional introns/insertions in combination with nuclear markers might be useful in distinguishing among the various species and subspecies of the O. ulmi s. lat. complex.

Mitochondrial DNA organisation of the mtDNA type 2b of Aspergillus tubingensis compared to the Aspergillus niger mtDNA type 1a

A physical and functional map of Aspergillus tubingensis mtDNA type 2b was constructed and compared to Aspergillus niger typeI a mtDNA. The gene content and order, as well as the patterns of restriction sites, were similar. Two unidentified ORFs and several repeat elements were found in the region between the ndhl and ndh4 genes on both mtDNAs. The sizes of the A. niger and A. tubingensis mtDNAs were 31.23 and 33.06 kb. respectively, the difference was principally attributed to the altered intron content of their coxl genes.

Interpretation of intraspecific variability in mtDNAs of Aspergillus niger strains and rearrangement of their mtDNAs following mitochondrial transmissions

Physical and functional maps of mitochondrial DNAs of Aspergillus niger strains representing different mitochondrial DNA RFLP patterns were constructed and compared. In spite of the high similarity in the organisation of mitochondrial DNAs among examined strains, differences could be easily recognised by applying molecular markers, such as the different intron content of the cox1 genes, the sequence of the intergenic regions between the Met- and His-tRNA genes and downstream of the tRNA-Gly gene. Intraspecific mitochondrial transfers between the heterokaryon incompatible mitochondrial oligomycin-resistant A. niger strain, as the donor, and other A. niger-sensitive strains bearing different RFLP patterns resulted in oligomycin-resistant progeny possessing either rearranged or unchanged donor mitochondrial DNA and recipient nuclei. Since the intergenic marker sequences of mitochondrial DNAs turned out to be identical in the donor and the progeny, it can be assumed that the oligomycin-resistant progeny inherit the mitochondrial DNA of the donor strain; this may either remain unchanged or may be modified by a mobile intron of the cox1 gene of the recipient mitochondria.

Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility

Homing endonucleases confer mobility to their host intervening sequence, either an intron or intein, by catalyzing a highly specific double-strand break in a cognate allele lacking the intervening sequence. These proteins are characterized by their ability to bind long DNA target sites (14-40 bp) and their tolerance of minor sequence changes in these sites. A wealth of biochemical and structural data has been generated for these enzymes over the past few years. Herein we review our current understanding of homing endonucleases, including their diversity and evolution, DNA-binding and catalytic mechanisms, and attempts to engineer them to bind novel DNA substrates.

Related homing endonucleases I-BmoI and I-TevI use different strategies to cleave homologous recognition sites

A typical homing endonuclease initiates mobility of its group I intron by recognizing DNA both upstream and downstream of the intron insertion site of intronless alleles, preventing the endonuclease from binding and cleaving its own intron-containing allele. Here, we describe a GIY-YIG family homing endonuclease, I-BmoI, that possesses an unusual recognition sequence, encompassing 1 base pair upstream but 38 base pairs downstream of the intron insertion site. I-BmoI binds intron-containing and intronless substrates with equal affinity but can nevertheless discriminate between the two for cleavage. I-BmoI is encoded by a group I intron that interrupts the thymidylate synthase (TS) gene (thyA) of Bacillus mojavensis s87-18. This intron resembles one inserted 21 nucleotides further downstream in a homologous TS gene (td) of Escherichia coli phage T4. I-TevI, the T4 td intron-encoded GIY-YIG endonuclease, is very similar to I-BmoI, but each endonuclease gene is inserted within a different position of its respective intron. Remarkably, I-TevI and I-BmoI bind a homologous stretch of TS-encoding DNA and cleave their intronless substrates in very similar positions. Our results suggest that each endonuclease has independently evolved the ability to distinguish intron-containing from intronless alleles while maintaining the same conserved recognition sequence centered on DNA-encoding active site residues of TS.

Intron mobility results in rearrangement in mitochondrial DNAs of heterokaryon incompatible Aspergillus japonicus strains after protoplast fusion

Mitochondrial transmission was carried out under selective conditions between incompatible Aspergillus japonicus strains always using an oligomycin-resistant mitochondrial donor and selecting for recipient nuclei and oligomycin-resistant mitochondria. All attempted intraspecific mitochondrial transmissions were successful, but the transmission between closely related A. japonicus and A. aculeatus failed. Under selection pressure, resistant progeny harbor the mitochondrial DNA (mtDNA) of the donor strain, which may remain unchanged or may be modified by the introns of the recipient mitochondrial genome. Detailed analysis of a certain strain harboring rearranged mtDNA suggests that the mtDNA profiles of recombinant-like progeny are strongly influenced by the characteristics and mobility of introns of both parental mtDNAs. Both intron loss and intron acquisition play a role in the rearrangement of mtDNA. In certain parental combinations, a particular intron was lost very frequently.