The mitochondrial genome of the basidiomycete Agrocybe aegerita: molecular cloning, physical mapping and gene location

Large inverted repeats identified by intra-specific comparison of mitochondrial genomes provide insights into the evolution of Agrocybe aegerita

Genomic structure and content of Agrocybe aegerita mitochondrial DNA contain essential information regarding the evolution of this gourmet mushroom. In this study, eight isolates of A. aegerita were sequenced and assembled into complete mitochondrial genomes. The mtDNA of the isolate Ag0067 contained two genotypes, both of which were quadripartite architecture consisting of two identical inverted repeats, separated by a small single-copy region and a large single-copy region. The only difference was opposite directions of the small single-copy region. The mtDNAs ranged from 116,329 bp to 134,035 bp, harboring two large identical inverted repeats. Genes of plasmid-origin were present in regions flanked by inverted repeat ID2. Most of the core genes evolved at a relatively low rate, whereas five tRNA genes located in corresponding regions of Ag0002:1-14000 and Ag0002:50001-61000 showed higher diversity. A long fragment inversion (10 Kb) was suggested to have occurred during the differentiation of two main clades, leading to two different gene orders. The number and distribution of the introns varied greatly among the A. aegerita mtDNAs. Fast invasion of short insertions likely resulted in the diversity of introns as well as other non-coding regions, increasing the variation of the mtDNAs. We raised a model about the evolution of the large repeats to explain the unusual features of A. aegerita mtDNAs. This study constructed quadripartite architecture of A. aegerita mtDNAs analogous to chloroplast DNA, proposed an interconversion model of the divergent mitochondrial genotypes with large inverted repeats. The findings could increase our knowledge of fungal evolution.

Karyotype analysis, genome organization, and stable genetic transformation of the root colonizing fungus Piriformospora indica

Piriformospora indica (Basidiomycota, Sebacinales) is a root colonizing fungus which is able to increase biomass and yield of crop plants and to induce local and systemic resistance to fungal diseases and tolerance to abiotic stress. A prerequisite for the elucidation of the mode of action of this novel kind of symbiosis is knowledge of the genome organization as well as the development of tools to study and modify gene functions. Here we provide data on the karyotype and genetic transformation strategies. The fungus was shown to possess at least six chromosomes and a genome size of about 15.4-24Mb. Sequences of the genes encoding the elongation factor 1-alpha (TEF) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were used for genome size estimation through real-time PCR analysis. Chromosomal location investigated by Southern blot and expression analysis suggested that TEF and GAPDH are single-copy genes with strong and constitutive promoters. A genetic transformation system was established using a fragment of the TEF promoter region for construction of vectors carrying the selectable marker hygromycin B phosphotransferase. Results demonstrate that P. indica can be stably transformed by random genomic integration of foreign DNA and that it posses a relative small genome as compared to other members of the Basidiomycota.

Organization of mitochondrial DNA in the basidiomycetous Dioszegia hungarica (Cryptococcus hungaricus) species

The organization of mitochondrial DNA was investigated in the six collection strains of the basidiomycetous yeast Dioszegia hungarica (Cryptococcus hungaricus) isolated so far. Physical and partial functional maps were constructed. Two strains (CBS 6324 and 6576) were identical while three others (CBS 4214, 5124, 6953) differed not only in the distribution of restriction sites but in gene order as well. Results confirm the hypothesis that these five strains are representatives of different closely related species. The sixth strain CBS 6569 revealed a unique mitochondrial genome organization. Its mtDNA separated into eight bands on agarose gel without enzymatic digestion. These molecules carried mitochondrial genes, and RFLP analysis of the four largest molecules using frequently-cutting restriction enzymes (KpnI and SmaI) showed them to have strongly homologous sequences. This unique mtDNA organization was also observed in a strain of Cystofilobasidium capitatum, providing evidence that CBS 6569 belongs to the Cystofilobasidium clade.

Duplication of a truncated paralog of the family B DNA polymerase gene Aa-polB in the Agrocybe aegerita mitochondrial genome

The Agrocybe aegerita mitochondrial genome contains a truncated family B DNA polymerase gene (Aa-polB P1) whose nucleotide sequence is 86% identical to the previously described and potentially functional Aa-polB gene. A tRNA(Met) gene occurs at the 3' end of the Aa-polB P1 gene. The Aa-polB P1 gene could result from reverse transcription of an Aa-polB mRNA primed by a tRNA(Met) followed by the integration of the cDNA after recombination at the mitochondrial tRNA locus. Two naturally occurring alleles of Aa-polB P1 carry one or two copies of the disrupted sequence. In strains with two copies of Aa-polB P1, these copies are inverted relative to one another and separated by a short sequence carrying the tRNA(Met) gene. Both A. aegerita mitochondrial family B DNA polymerases were found to be related to other family B DNA polymerases (36 to 53% amino acid similarity), including the three enzymes of the archaebacterium Sulfolobus solfataricus. If mitochondria originated from a fusion between a Clostridium-like eubacterium and a Sulfolobus-like archaebacterium, then the A. aegerita family B DNA polymerase genes could be remnants of the archaebacterial genes.

Sequence and secondary structure of the mitochondrial small-subunit rRNA V4, V6, and V9 domains reveal highly species-specific variations within the genus Agrocybe

A comparative study of variable domains V4, V6, and V9 of the mitochondrial small-subunit (SSU) rRNA was carried out with the genus Agrocybe by PCR amplification of 42 wild isolates belonging to 10 species, Agrocybe aegerita, Agrocybe dura, Agrocybe chaxingu, Agrocybe erebia, Agrocybe firma, Agrocybe praecox, Agrocybe paludosa, Agrocybe pediades, Agrocybe alnetorum, and Agrocybe vervacti. Sequencing of the PCR products showed that the three domains in the isolates belonging to the same species were the same length and had the same sequence, while variations were found among the 10 species. Alignment of the sequences showed that nucleotide motifs encountered in the smallest sequence of each variable domain were also found in the largest sequence, indicating that the sequences evolved by insertion-deletion events. Determination of the secondary structure of each domain revealed that the insertion-deletion events commonly occurred in regions not directly involved in the secondary structure (i.e., the loops). Moreover, conserved sequences ranging from 4 to 25 nucleotides long were found at the beginning and end of each domain and could constitute genus-specific sequences. Comparisons of the V4, V6, and V9 secondary structures resulted in identification of the following four groups: (i) group I, which was characterized by the presence of additional P23-1 and P23-3 helices in the V4 domain and the lack of the P49-1 helix in V9 and included A. aegerita, A. chaxingu, and A. erebia; (ii) group II, which had the P23-3 helix in V4 and the P49-1 helix in V9 and included A. pediades; (iii) group III, which did not have additional helices in V4, had the P49-1 helix in V9 and included A. paludosa, A. firma, A. alnetorum, and A. praecox; and (iv) group IV, which lacked both the V4 additional helices and the P49-1 helix in V9 and included A. vervacti and A. dura. This grouping of species was supported by the structure of a consensus tree based on the variable domain sequences. The conservation of the sequences of the V4, V6, and V9 domains of the mitochondrial SSU rRNA within species and the high degree of interspecific variation found in the Agrocybe species studied open the way for these sequences to be used as specific molecular markers of the Basidiomycota.

Molecular analysis of the split cox1 gene from the Basidiomycota Agrocybe aegerita: relationship of its introns with homologous Ascomycota introns and divergence levels from common ancestral copies

The Basidiomycota Agrocybe aegerita (Aa) mitochondrial cox1 gene (6790 nucleotides), encoding a protein of 527aa (58377Da), is split by four large subgroup IB introns possessing site-specific endonucleases assumed to be involved in intron mobility. When compared to other fungal COX1 proteins, the Aa protein is closely related to the COX1 one of the Basidiomycota Schizophyllum commune (Sc). This clade reveals a relationship with the studied Ascomycota ones, with the exception of Schizosaccharomyces pombe (Sp) which ranges in an out-group position compared with both higher fungi divisions. When comparison is extended to other kingdoms, fungal COX1 sequences are found to be more related to algae and plant ones (more than 57.5% aa similarity) than to animal sequences (53.6% aa similarity), contrasting with the previously established close relationship between fungi and animals, based on comparisons of nuclear genes. The four Aa cox1 introns are homologous to Ascomycota or algae cox1 introns sharing the same location within the exonic sequences. The percentages of identity of the intronic nucleotide sequences suggest a possible acquisition by lateral transfers of ancestral copies or of their derived sequences. These identities extend over the whole intronic sequences, arguing in favor of a transfer of the complete intron rather than a transfer limited to the encoded ORF. The intron i4 shares 74% of identity, at the nucleotidic level, with the Podospora anserina (Pa) intron i14, and up to 90.5% of aa similarity between the encoded proteins, i.e. the highest values reported to date between introns of two phylogenetically distant species. This low divergence argues for a recent lateral transfer between the two species. On the contrary, the low sequence identities (below 36%) observed between Aa i1 and the homologous Sp i1 or Prototheca wickeramii (Pw) i1 suggest a long evolution time after the separation of these sequences. The introns i2 and i3 possessed intermediate percentages of identity with their homologous Ascomycota introns. This is the first report of the complete nucleotide sequence and molecular organization of a mitochondrial cox1 gene of any member of the Basidiomycota division.

Genetic evidence for nonrandom sorting of mitochondria in the basidiomycete Agrocybe aegerita

We studied mitochondrial transmission in the homobasidiomycete Agrocybe aegerita during plasmogamy, vegetative growth, and basidiocarp differentiation. Plasmogamy between homokaryons from progeny of three wild-type strains resulted in bidirectional nuclear migration, and the dikaryotization speed was dependent on the nuclear genotype of the recipient homokaryon. Little mitochondrial migration accompanied the nuclear migration. A total of 75% of the dikaryons from the fusion lines had both parental mitochondrial haplotypes (mixed dikaryons), and 25% had only a single haplotype (homoplasmic dikaryons); with some matings, there was a strong bias in favor of one parental haplotype. We demonstrated the heteroplasmic nature of mixed dikaryons by (i) isolating and subculturing apical cells in micromanipulation experiments and (ii) identifying recombinant mitochondrial genomes. This heteroplasmy is consistent with the previously reported suggestion that there is recombination between mitochondrial alleles in A. aegerita. Conversion of heteroplasmons into homoplasmons occurred (i) during long-term storage, (ii) in mycelia regenerated from isolated apical cells, and (iii) during basidiocarp differentiation. Homokaryons that readily accepted foreign nuclei were the most efficient homokaryons in maintaining their mitochondrial haplotype during plasmogamy, long-term storage, and basidiocarp differentiation. This suggests that the mechanism responsible for the nonrandom retention or elimination of a given haplotype may be related to the nuclear genotype or the mitochondrial haplotype or both.

DNA sequence and secondary structure of the mitochondrial small subunit ribosomal RNA coding region including a group-IC2 intron from the cultivated basidiomycete Agrocybe aegerita

Due to their structural complexity and their evolutionary dimension, rRNAs are the most investigated nucleic acids in prokaryotes, eukaryotes and organelles. However, no complete sequence of a mitochondrial small subunit (SSU) rRNA was available in the basidiomycotina subdivision. The mitochondrial gene encoding the SSU rRNA of the cultivated basidiomycete Agrocybe aegerita was cloned and its complete nucleotide sequence achieved; the 5'- and 3'-ends were localized by nuclease SI mapping, leading to a size of 3277 nt. The secondary structure of the SSU rRNA (1906 nt in size) possessed all the helices and loops of the prokaryotic model; a unique modification was found in a conserved nucleotide predicted by the model: the nt 487 was A instead of C. The same modification, has been found in all the partial basidiomycete mitochondrial sequences available in databases. The Agrocybe aegerita SSU rRNA was characterized by large and unusual extensions leading to additional helices in the variable domains V4, V6 and V9, which were the longest of the known prokaryotic or mitochondrial SSU rRNAs. Nucleotide sequence analysis indicated a 1371-bp intron, belonging to subgroup-IC2, located in a conserved loop in the 3'-part of the SSU rRNA. This intron, which is the second example reported in a fungal mitochondrial SSU rDNA, encoded a putative protein (407 aa) sharing homologies with endonucleases involved in group-I intron mobility. This report constitutes the first complete mitochondrial SSU rRNA sequence and secondary structure of any member of the basidiomycotina subdivision.

Wide Distribution of Mitochondrial Genome Rearrangements in Wild Strains of the Cultivated Basidiomycete Agrocybe aegerita

We used restriction fragment length polymorphisms to examine mitochondrial genome rearrangements in 36 wild strains of the cultivated basidiomycete Agrocybe aegerita, collected from widely distributed locations in Europe. We identified two polymorphic regions within the mitochondrial DNA which varied independently: one carrying the Cox II coding sequence and the other carrying the Cox I, ATP6, and ATP8 coding sequences. Two types of mutations were responsible for the restriction fragment length polymorphisms that we observed and, accordingly, were involved in the A. aegerita mitochondrial genome evolution: (i) point mutations, which resulted in strain-specific mitochondrial markers, and (ii) length mutations due to genome rearrangements, such as deletions, insertions, or duplications. Within each polymorphic region, the length differences defined only two mitochondrial types, suggesting that these length mutations were not randomly generated but resulted from a precise rearrangement mechanism. For each of the two polymorphic regions, the two molecular types were distributed among the 36 strains without obvious correlation with their geographic origin. On the basis of these two polymorphisms, it is possible to define four mitochondrial haplotypes. The four mitochondrial haplotypes could be the result of intermolecular recombination between allelic forms present in the population long enough to reach linkage equilibrium. All of the 36 dikaryotic strains contained only a single mitochondrial type, confirming the previously described mitochondrial sorting out after cytoplasmic mixing in basidiomycetes.

Mitochondrial DNAs of the fungus Armillaria ostoyae: restriction map and length variation

A restriction-enzyme-site map is presented for the 147-kb mtDNA of North American Armillaria ostoyae. The locations of five structural genes, atp6, atp8, coxI, coxIII, and cob, along with the location and orientation of the large and small ribosomal RNA genes, were determined through Southern hybridizations with cloned genes from other fungal mtDNAs. Based on this map, the variation in mtDNA suggested geographic structure at two different levels. On a large geographic scale, 17 mtDNA types from North America were distinct, with respect to both size and restriction maps, from three mtDNA types from Europe. At the local scale, identical mtDNA types were evident among several different genetic individuals located no more than 1 km apart at a site in Michigan. No mtDNA type occurred more than once among genetic individuals from different regions of North America, although the occurrence of similar mtDNAs in isolates from distant regions suggested that this may occur at a low frequency with large sample sizes. Among the North American mtDNA types, analysis of discrete length variants was inconsistent with the hypothesis that the mtDNA of A. ostoyae evolves as a clonal lineage in which each length mutation represents a unique event. The two remaining hypotheses, that similar mutational events have occurred independently and that genetic exchange and recombination occurs among mtDNAs in natural populations of this species, remain to be tested.

Homologous transformation of the edible basidiomycete Agrocybe aegerita with the URA1 gene: characterization of integrative events and of rearranged free plasmids in transformants

The URA1 gene, encoding dihydroorotate dehydrogenase of the pyrimidine pathway, cloned into pUC18 (pUra1-1) was used to develop an homologous transformation system for the cultivated mushroom Agrocybe aegerita. Protoplasts of a ura1 auxotrophic strain were transformed by electroporation with efficiencies ranging from 1 to 26 transformants per micrograms of DNA. The phenotype of the stable Ura+ transformants suggested a strong nuclear heterogeneity further confirmed by Southern-blot analysis. All transformants acquired extrachromosomal forms derived from pUra1-1. Integration of pUra1-1 into chromosomal DNA occurred for some transformants. Plasmids containing the integrant of pUC18 recombined to different parts of the URA1 gene were rescued from A. aegerita transformants through transformation of E. coli. Their molecular analysis indicated that they represent products of the continuous excision of primary-integrated vector sequences rather than ARS-dependent autoreplicative forms.

Mitochondrial genome of the dimorphic zygomycete Mucor racemosus

Mitochondria were isolated from the dimorphic zygomycete Mucor racemosus by differential centrifugation. DNA from the organelles was purified by cesium chloride-ethidium bromide isopycnic centrifugation. Examination of the mitochondrial DNA by electron microscopy revealed a circular chromosome approximately 63.8 kbp in circumference. The chromosome was digested with restriction endonucleases and the resulting DNA fragments were separated by agarose-gel electrophoresis. Electophoretic mobilities and stoichiometry of the fragments indicated a mixed population of mtDNA molecules each with a size of about 63.4 kbp. Physical maps were constructed from analyses of fragments generated in single and double restriction digests and from the hybridization of fragments to probes for the large and small mitochondrial rRNA genes from Saccharomyces cerevisiae. The Mucor mitochondrial chromosome was found to exist in the form of two flip-flop isomers with inverted repeat sequences encoding both rRNA genes.