Evolution of Fusarium tricinctum and Fusarium avenaceum mitochondrial genomes is driven by mobility of introns and of a new type of palindromic microsatellite repeats.

BACKGROUND: Increased contamination of European and Asian wheat and barley crops with "emerging" mycotoxins such as enniatins or beauvericin, produced by Fusarium avenaceum and Fusarium tricinctum, suggest that these phylogenetically close species could be involved in future food-safety crises. RESULTS: The mitochondrial genomes of F. tricinctum strain INRA104 and F. avenaceum strain FaLH27 have been annotated. A comparative analysis was carried out then extended to a set of 25 wild strains. Results show that they constitute two distinct species, easily distinguished by their mitochondrial sequences. The mitochondrial genetic variability is mainly located within the intergenic regions. Marks of variations show they have evolved (i) by Single Nucleotide Polymorphisms (SNPs), (ii) by length variations mediated by insertion/deletion sequences (Indels), and (iii) by length mutations generated by DNA sliding events occurring in mononucleotide (A)n or (T)n microsatellite type sequences arranged in a peculiar palindromic organization. The optionality of these palindromes between both species argues for their mobility. The presence of Indels and SNPs in palindrome neighbouring regions suggests their involvement in these observed variations. Moreover, the intraspecific and interspecific variations in the presence/absence of group I introns suggest a high mobility, resulting from several events of gain and loss during short evolution periods. Phylogenetic analyses of intron orthologous sequences suggest that most introns could have originated from lateral transfers from phylogenetically close or distant species belonging to various Ascomycota genera and even to the Basidiomycota fungal division. CONCLUSIONS: Mitochondrial genome evolution between F. tricinctum and F. avenaceum is mostly driven by two types of mobile genetic elements, implicated in genome polymorphism. The first one is represented by group I introns. Indeed, both genomes harbour optional (inter- or intra-specifically) group I introns, all carrying putatively functional hegs, arguing for a high mobility of these introns during short evolution periods. The gain events were shown to involve, for most of them, lateral transfers between phylogenetically distant species. This study has also revealed a new type of mobile genetic element constituted by a palindromic arrangement of (A) n and (T) n microsatellite sequences whose presence was related to occurrence of SNPs and Indels in the neighbouring regions.

The intraspecific variability of mitochondrial genes of Agaricus bisporus reveals an extensive group I intron mobility combined with low nucleotide substitution rates.

Intraspecific mitochondrial variability was studied in ten strains of A. bisporus var. bisporus, in a strain representative of A. bisporus var. eurotetrasporus and in a strain of the closely related species Agaricus devoniensis. In A. bisporus, the cox1 gene is the richest in group I introns harboring homing endonuclease genes (heg). This study led to identify group I introns as the main source of cox1 gene polymorphism. Among the studied introns, two groups were distinguished according to the heg they contained. One group harbored heg maintained putatively functional. The other group was composed of eroded heg sequences that appeared to evolve toward their elimination. Low nucleotide substitution rates were found in both types of intronic sequences. This feature was also shared by all types of studied mitochondrial sequences, not only intronic but also genic and intergenic ones, when compared with nuclear sequences. Hence, the intraspecific evolution of A. bisporus mitochondrial genome appears characterized by both an important mobility (presence/absence) of large group I introns and by low nt substitution rates. This stringent conservation of mitochondrial sequences, when compared with their nuclear counterparts, appears irrespective of their apparent functionality and contrasts to what is widely accepted in fungal sequence evolution. This strengthens the usefulness of mtDNA sequences to get clues on intraspecific evolution.

The 135 kbp mitochondrial genome of Agaricus bisporus is the largest known eukaryotic reservoir of group I introns and plasmid-related sequences.

At 135,005 nt, the mitochondrial genome in Agaricus bisporus represents the largest fungal mitochondrial genome sequenced to date. Its large size is mainly due to the presence of mobile genetic elements, including a total of 43 group I introns, three group II introns, and five DNA fragments that show sequence similarity to linear invertron-like plasmids. The introns are distributed in eight of the 15 protein coding genes. These introns contain a total of 61,092 nt (∼45.3% of the whole mitochondrial genome) and include representatives of most of the group I introns so far found in mitochondrial genomes of Basidiomycota. The plasmid-like sequences include 6730 nt total representing 5.0% of the genome. These sequences showed high-level similarities to two different mitochondrial plasmids reported for basidiomycete mushrooms: the autonomously replicating pEM in Agaricus bitorquis and the integrated linear plasmid sequences in Agrocybe aegerita and Moniliophthora perniciosa. Moreover, the plasmid-related sequences are located within or adjacent to two large (4559 nt) inverted repeats containing also two sets of mitochondrial tRNA genes. Our analyses are consistent with the hypothesis that horizontal DNA transfer has played a significant role in the evolution of the A. bisporus mitochondrial genome.

Experimental Outcrossing in Agaricus bisporus Revealed a Major and Unexpected Involvement of Airborne Mycelium Fragments.

Agaricus bisporus var. bisporus, the button mushroom, has a predominantly pseudohomothallic life cycle. Most of its spores are heterokaryotic and give rise to fertile heterokaryons. However, previous studies have suggested that outcrossing should not be rare in wild populations. In order to discover how outcrossing occurs, we experimentally favored it between aerial propagules of a fruiting donor mycelium and a delayed receiver mycelium that only invaded culture trays. Two donor/receiver pairs were studied, and potentially hybrid basidiomata collected on the receiver trays were analyzed with a mitochondrial marker, two unlinked nuclear CAPS markers, then haplotype markers based on DNA sequences obtained after PCR cloning of the rDNA ITS region and the fruk gene. For one of the two pairs, most basidiomata were hybrids between the donor and the receiver. Genotyping of the hybrids revealed only two genotypes consistent with outcrossing involving airborne mycelium fragments rather than basidiospores. The resident receiver heterokaryon that provided its mitochondria to the hybrid basidiomata is suspected to have had a trophic contribution to their growth and successful fruiting. The high level of heterozygosity and the cultivar introgression previously revealed in wild populations of this pseudohomothallic species may result from outcrossing involving airborne pieces of mycelium.

Genome Sequence of the Emerging Mycotoxin-Producing Filamentous Fungus Fusarium tricinctum Strain INRA104.

The genome of the phytopathogenic fungus Fusarium tricinctum strain INRA104 was sequenced at a fold-coverage of more than 500×. This led to 23 scaffolds, including one scaffold for the mitochondrial genome, for a total genome size of 42.8 Mb, with an average GC content of 45% and 13,387 predicted genes.

Unusual accumulation of polymorphic microsatellite loci in a specific region of the mitochondrial genome of two mushroom-forming Agrocybe species.

The cob/tRNA(Tyr) mitochondrial regions of Agrocybe aegerita and of the related species Agrocybe chaxingu display an unusual clustering of four microsatellite loci constituted by motifs of one to six nucleotides whose number of repeats varied from three to 18. In A. chaxingu, these microsatellite loci are followed by a small region bearing one additional microsatelite and one minisatellite locus constituted by an octanucleotide motif repeated 13-18 times. In A. aegerita, this latter region is deleted. This is the first evidence of such an accumulation of microsatellites in mitochondrial genomes. The analyses of the microsatellite loci in 11 A. aegerita and in four A. chaxingu wild strains have shown extensive intraspecific and interspecific variations in the number of tandem repeats (VNTRs), suggesting that these loci could represent powerful molecular markers for strain fingerprinting. Up to 23 different alleles were present in the 15 Agrocybe studied strains, allowing the definition of 12 different haplotypes.

A CAPS test allowing a rapid distinction of Penicillium expansum among fungal species collected on grape berries, inferred from the sequence and secondary structure of the mitochondrial SSU-rRNA.

Penicillium expansum is a fungal species highly damageable for the postharvest conservation of numerous fruits. In vineyards, this fungus is sometimes isolated from grape berries where its presence may lead to the production of geosmin, a powerful earthy odorant, which can impair grapes and wines aromas. However, the discrimination of P. expansum from related fungi is difficult because it is based on ambiguous phenotypic characters and/or expensive and time-consuming molecular tests. In this context, the complete sequences and secondary structures of Penicillium expansum and Penicillium thomii mitochondrial SSU-rRNAs were achieved and compared with those of two other phylogenetically related Ascomycota: Penicillium chrysogenum and Emericella nidulans. The comparison has shown a high conservation in size and sequence of the core and of the variable domains (more than 80% of nt identity) of the four SSU-rRNAs, arguing for a close phylogenetic relationship between these four species of the Trichocomaceae family. Large (from 10 to 18 nt) inserted/deleted (indel) sequences were evidenced in the V1, V5 and V6 variable domains. The size variations (10 to 18 nt) of the V1 indel sequence allowed the distinction of the four species; the V5 indel (15 nt) was specifically recovered in E. nidulans; the V6 indel (16 nt), shared by the three Penicillium species, was lacking in E. nidulans. A couple of conserved primers (UI/R2) were defined to generate a PCR product containing the V1 to V5 variable domains. This product contained the two regions of the four SSU-rRNAs showing the highest rates of nt substitutions, namely the V2 variable domain and, surprisingly, a helix (H17) of the core. The H17 sequence was shown to specifically possess in P. expansum a recognition site for the ClaI restriction endonuclease. Hence, this enzyme generates a digestion pattern of the PCR product with two bands (350 bp+500 bp), specific to P. expansum and easily separable by agarose gel electrophoresis. This leads to a CAPS test usable for P. expansum discrimination among grape berries isolated filamentous fungi. The CAPS test was validated by a comparative analysis involving 29 strains belonging to 17 species currently isolated from grape berries in the Bordeaux vineyards.

A molecular contribution to the assessment of the Tricholoma equestre species complex.

In recent years, interest in the Tricholoma equestre species complex has increased because of several cases of severe and sometimes fatal rhabdomyolysis reported in France and Poland. These occurred after repeated consumption of large portions of T. equestre sporophores during consecutive meals, despite the fact that this species is renowned as a tasty edible wild mushroom. The T. equestre species complex includes three ectomycorrhizal species Tricholoma flavovirens (Pers.) S. Lundell, Tricholoma auratum (Paulet) Gillet, and T. equestre (L.) P. Kummer. All these species produce sporophores with intense yellow gills but are difficult to distinguish by morphological analyses at both the macroscopic and microscopic levels. In T. equestre, two additional varieties are recognized: T. equestre var. populinum (Christensen & Noordeloos) associated with Populus sp. and/or Betula sp. trees and sometimes recognized as Tricholoma frondosae (Kalamees & Shchukin) and T. equestre var. pallidifolia characterized by pale to white gills, frequently recognized as Tricholoma joachimii (Bon & Riva). To explore the taxonomic (species delimitation), ecological, and geographical extent and limits of the T. equestre species complex, we have carried out a molecular comparison of worldwide strains belonging to this complex by using sequences of two molecular markers: the internal transcript spacer (ITS)1/5.8S/ITS2 region of the nuclear ribosomal unit and the 5' part of the mitochondrial cox1 gene. Phylogenetic analyses support the placement of European T. equestre, T. flavovirens, and T. auratum strains as representatives of a single species. This species appears associated with various conifers trees, depending on the geographic origin (Pinus pinaster for T. auratum, Pinus sylvestris or Abies alba for T. equestre and T. flavovirens). However, in the context of a single T. equestre species, the geographical location could lead to the characterization of sub-species or varieties, as suggested by the gathering of the four Asian (Japanese) T. auratum strains in a strongly supported distinct phylogenetic clade. Moreover, our analysis strongly argues for considering T. joachimii and the synonymised T. equestre var. pallidifolia as two representatives of a different species not belonging to the T. equestre group. This species would be phylogenetically related to the Tricholoma columbetta species with which they share white gills. Similarly, the phylogenetic analysis of the molecular data and the lack of gene flow between the strains associated with broad-leaved trees and those of the T. equestre complex, rather argues for two distinct species depending on the ecological niche: T. frondosae under broad-leaved trees and T. equestre under conifers.

Performance of the COX1 gene as a marker for the study of metabolically active Pezizomycotina and Agaricomycetes fungal communities from the analysis of soil RNA.

In temperate forest soils, filamentous ectomycorrhizal and saprotrophic fungi affiliated to the Agaricomycetes and Pezizomycotina contribute to key biological processes. The diversity of soil fungal communities is usually estimated by studying molecular markers such as nuclear ribosomal gene regions amplified from soil-extracted DNA. However, this approach only reveals the presence of the corresponding genomic DNA in the soil sample and may not reflect the diversity of the metabolically active species. To circumvent this problem, we investigated the performance of the mitochondrial cytochrome c oxidase 1 (COX1)-encoding gene as a fungal molecular marker for environmental RNA-based studies. We designed PCR primers to specifically amplify Agaricomycetes and Pezizomycotina COX1 partial sequences and amplified them from both soil DNA and reverse-transcribed soil RNA. As a control, we also amplified the nuclear internal transcribed spacer ribosomal region from soil DNA. Fungal COX1 sequences were readily amplified from soil-extracted nucleic acids and were not significantly contaminated by nontarget sequences. We show that the relative abundance of fungal taxonomic groups differed between the different sequence data sets, with for example ascomycete COX1 sequences being more abundant among sequences amplified from soil DNA than from soil cDNAs.

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.

The mitochondrial apocytochrome b genes of two Agrocybe species suggest lateral transfers of group I homing introns among phylogenetically distant fungi.

The Agrocybe chaxingu and Agrocybe aegerita mitochondrial apocytochrome b coding sequences are highly similar (97% of nt identity), but have highly different sizes (2312 and 4867nt, respectively), due to the presence of three large group IB introns: two (iAae1 and iAae2) in A. aegerita, one (iAch1) in A. chaxingu. All these introns encode a homing endonuclease (HE) similar to those described in introns of mitochondrial genes (cob, cox1, and nad5) from various organisms. Phylogenetic trees were built with these HE sequences. From these trees, the Agrocybe coding introns argue for recent lateral transfers, i.e., occurring after the separation of the two Agrocybe species, involving phylogenetically distant fungi such as members of the Ascomycota phylum (for iAch1 and iAae2) and, for the first time to our knowledge, a member of the Chytridiomycota phylum (for iAae1). The grouping of the HE gene (HEG) sequences according to the mitochondrial gene (cob, cox1, and nad5) where they are inserted, suggests modifications of the interactions between the HE and the recognized sequences, leading to new target genes. The largest distribution of the iAch1 HE, shared by several cob and cox1 mitochondrial genes from Ascomycota, Basidiomycota, and Chytridiomycota phyla, suggests a higher target flexibility of this HE, perhaps related to the presence of two different LAGLIDADG motifs in the catalytic site of the enzyme.

Secondary structure and molecular evolution of the mitochondrial small subunit ribosomal RNA in Agaricales (Euagarics clade, Homobasidiomycota).

The complete sequences and secondary structures of the mitochondrial small subunit (SSU) ribosomal RNAs of both mostly cultivated mushrooms Agaricus bisporus (1930 nt) and Lentinula edodes (2164 nt) were achieved. These secondary structures and that of Schizophyllum commune (1872 nt) were compared to that previously established for Agrocybe aegerita. The four structures are near the model established for Archae, Bacteria, plastids, and mitochondria; particularly the helices 23 and 37, described as specific to bacteria, are present. Within the four Agaricales (Homobasidiomycota), the SSU-rRNA "core" is conserved in size (966 to 1009 nt) with the exception of an unusual extension of 40 nt in the H17 helix of S. commune. The four core sequences possess 76% of conserved positions and a cluster of C in their 3' end, which could constitute a signal involved in the RNA maturation process. Among the nine putative variable domains, three (V3, V5, V7) do not show significant length variations and possess similar percentages of conserved positions (69%) than the core. The other six variable domains show important length variations, due to independent large size inserted/deleted sequences, and higher rates of nucleotide substitutions than the core (only 31% of conserved positions between the four species). Interestingly, the inserted/deleted sequences are located in few preferential sites (hot spots for insertion/deletion) where they seem to arise or disappear haphazardly during evolution. These sites are located on the surface of the tertiary structure of the 30S ribosomal subunit, at the beginning of hairpin loops; the insertions lead to a lengthening of existing hairpins or to branching loops bearing up to five additional helices.