The mitochondrial DNAs of Agaricus: heterogeneity in A. bitorquis and homogeneity in A. brunnescens

Insight into the evolutionary and domesticated history of the most widely cultivated mushroom Agaricus bisporus via mitogenome sequences of 361 global strains

Agaricus bisporus is the most widely cultivated edible mushroom in the world with a only around three hundred years known history of cultivation. Therefore, it represents an ideal organism not only to investigate the natural evolutionary history but also the understanding on the evolution going back to the early era of domestication. In this study, we generated the mitochondrial genome sequences of 352 A. bisporus strains and 9 strains from 4 closely related species around the world. The population mitogenomic study revealed all A. bisporus strains can be divided into seven clades, and all domesticated cultivars present only in two of those clades. The molecular dating analysis showed this species origin in Europe on 4.6 Ma and we proposed the main dispersal routes. The detailed mitogenome structure studies showed that the insertion of the plasmid-derived dpo gene caused a long fragment (MIR) inversion, and the distributions of the fragments of dpo gene were strictly in correspondence with these seven clades. Our studies also showed A. bisporus population contains 30 intron distribution patterns (IDPs), while all cultivars contain only two IDPs, which clearly exhibit intron loss compared to the others. Either the loss occurred before or after domestication, that could suggest that the change facilitates their adaptation to the cultivated environment.

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.

Molecular evolution of mitochondrial ribosomal DNA in the fungal genus Tricholoma: barcoding implications

The molecular evolution of the V6 and V9 domains of the mitochondrial SSU-rDNA was investigated to evaluate the use of these sequences for DNA barcodes in the Basidiomycota division. The PCR products from 27 isolates belonging to 11 Tricholoma species were sequenced. Both domains in the isolates belonging to the same species had identical sequences. All the species possess distinctive V9 sequences due to point mutations and insertion/deletion events. Secondary structures revealed that the insertion-deletion events occurred in regions not directly involved in the maintenance of the standard SSU-rRNA structure. The inserted sequences possess conserved motifs that enable their alignment among phylogenetically distant species. Hence, the V9 domain by displaying identical sequences within species, an adequate divergence level, easy amplification, and alignment represents an alternative molecular marker for the Basidiomycota division and opens the way for this sequence to be used as specific molecular markers of the fungal kingdom.

Mapping and characterization of polymorphism in mtDNA of Cryphonectria parasitica: evidence of the presence of an optional intron

The mitochondrial DNA (mtDNA) of the filamentous ascomycete Cryphonectria parasitica is large and polymorphic so, to better understand the nature of the polymorphisms within populations, a small collection of Italian strains of the fungus was examined. Known mtDNA polymorphisms were mapped and found to cluster in four regions of the mtDNA molecule, particularly in the RFLP region 2 where five different mtDNA haplotypes out of 13 strains were identified. This region included an area of 8.4kbp which was entirely sequenced in strain Ep155 showing the presence of two introns. An internal 3.2kbp portion was sequenced also in six additional strains. Sequence comparison of the C. parasitica mitochondrial intronic ORFs revealed similarities to known endonucleases such as those of Podospora anserina and Neurospora crassa. DNA sequence analysis showed that three polymorphisms of this mtDNA region within this population of 12 strains were due to the optional presence in the ND5 gene of an intron and of an intervening sequence within the intron. Evidence was also found within this population of mixed mitochondrial types within a single strain.

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.

The plant mitochondrial genome: homologous recombination as a mechanism for generating heterogeneity

The mitochondrial genomes of higher plants are among the largest and most complex organelle genomes described. They are generally multicircular or partly linear; in some species, extrachromosomal plasmids are present. It is proposed that inter- and intramolecular homologous recombination can account for the diversity of the observed genome organizations. The ability of mitochondria to fuse establishes a panmictic mitochondrial DNA population which is in recombinational equilibrium. It is suggested that this suppresses the base mutation rate, and unequal partitioning of the cytoplasm during cell division can lead to the rapid evolution of mitochondrial genome structure. This contrasts with the observed rates of base-sequence and genome evolution in chloroplasts. This difference can be accounted for solely by the inability of chloroplasts to fuse, thereby preventing chloroplast genome panmixis.

Plasmid RNA polymerase-like mitochondrial sequences in Agaricus bitorquis

A linear mitochondrial plasmid, pEM, found in certain isolates of the basidiomycete Agaricus bitorquis, potentially encodes virus-like DNA and RNA polymerases. Mitochondrial DNA from Agaricus bisporus that hybridizes to an internal region of pEM contains a fragmented and potentially non-functional version of the carboxy terminal end of the plasmid RNA polymerase. In this study, we present the sequence of the corresponding region of mitochondrial DNA from A. bitorquis. This sequence contained the same region of the plasmid RNA polymerase gene as was reported for the mitochondrial DNA of A. bisporus, and the level of similarity between the A. bisporus and A. bitorquis mitochondrial sequences was much higher than the level of similarity between either mitochondrial sequence and the plasmid. We propose that this plasmid RNA polymerase-like sequence was present in the Agaricus mitochondrial genome before the divergence of A. bisporus and A. bitorquis, and thus is unlikely to be a recent derivative of the plasmid pEM.

Isolation of a species-specific mitochondrial DNA sequence for identification of Tilletia indica, the Karnal bunt of wheat fungus

Mitochondrial DNA (mtDNA) from five isolates of Tilletia indica was isolated and digested with several restriction enzymes. A 2.3-kb EcoRI fragment was chosen, cloned, and shown to hybridize with total DNA restricted with EcoRI from T. indica and not from a morphologically similar smut fungus, Tilletia barclayana. The clone was partially sequenced, and primers were designed and tested under high-stringency conditions in PCR assays. The primer pair Ti1/Ti4 amplified a 2.3-kb fragment from total DNA of 17 T. indica isolates from India, Pakistan, and Mexico. DNA from 25 isolates of other smut fungi (T. barclayana, Tilletia foetida, Tilletia caries, Tilletia fusca, and Tilletia controversa) did not produce any bands, as detected by ethidium bromide-stained agarose gels and Southern hybridizations. The sensitivity of the assay was determined and increased by using a single nested primer in a second round of amplification, so that 1 pg of total mycelial DNA could be detected. The results indicated that the primers which originated from a cloned mtDNA sequence can be used to differentiate T. indica from other Tilletia species and have the potential to identify teliospores contaminating wheat seeds.

The utility of mitochondrial DNA restriction analysis in the classification of strains of Chrysosporium (hyphomycetes)

The taxonomy of the fungal genus Chrysosporium is mainly based on morphological features. In our current studies we have found several Chrysosporium species which showed intermediate morphological characteristics between several species. For this reason, we have carried out an analysis of the mitochondrial DNA restriction fragments of these strains that have permit us to classify each isolate strain in a species.

Physical map and gene organization of the mitochondrial genome from the unicellular green alga Platymonas (Tetraselmis) subcordiformis (Prasinophyceae)

The entire mitochondrial genome (mt genome) of the unicellular green alga Platymonas subcordiformis (synonym Tetraselmis subcordiformis; Prasinophyceae) was cloned and a physical map for the four restriction enzymes Hind III, Eco RI, Bg/II and Xba I was constructed. The mt genome of P. subcordiformis is a 42.8 kb circular molecule, coding for at least 23 genes. Hybridization and sequence analysis revealed the presence of a ca. 1.5 kb inverted repeat on the mt genome of P. subcordiformis. Phylogenetic analyses based on sequences of several coxI genes were carried out. Our data indicate that mitochondria from P. subcordiformis and from land plants form a natural, monophyletic group.

Indigenous and introduced populations ofAgaricus bisporus, the cultivated button mushroom, in eastern and western Canada: implications for population biology, resource management, and conservation of genetic diversity

Agaricus bisporus is known from field collections in several parts of southern Canada and the border states of Washington and Idaho. In Ontario, the species is associated with urban horticultural sites, agricultural areas, and thoroughfares. In British Columbia, Washington, and Idaho, the species is less well documented but occurs in habitats similar to those in Ontario. All studied isolates from these two regions were genotypically similar to those collected in Europe and also to a representative sample of cultivar isolates believed to be of European origin. In contrast, a population from forests of Picea in the Rocky Mountains of Alberta had genotypes that are very different from European and cultivar groups, including the putatively introduced populations in other parts of Canada. To a lesser extent, the Albertan individuals also differed from three other isolated natural populations, two in California and one in Israel. Two of 35 isolates from Alberta had nuclear genotypes that were partially or entirely consistent with European ancestry, indicating that foreign (probably cultivar) germ plasm has become established in the native habitat. Thirty of 35 isolates from the Alberta field sample had one of six mitochondrial DNA (mtDNA) haplotypes known only from Alberta; the remaining five had an mtDNA type known from cultivar strains. However, nuclear genotypic similarity to other Albertan or European-cultivar isolates was not significantly correlated with presence or absence of this latter mtDNA type. The genetic diversity represented by the previously unknown Albertan population, and to a lesser extent by the heirloom varietal lineages and their derivatives that have become established in Ontario, British Columbia, Washington, and Idaho, is an important resource for breeders of this valuable fungal crop species. Threats to resource populations of this species, and possible responses, are discussed. Key words: Agaricus bisporus, mushroom population biology, resource management, germ plasm conservation, genetic diversity, microphylogeny.

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.

RFLPs in mitochondrial and nuclear DNA indicate low levels of genetic diversity in the oak wilt pathogen Ceratocystis fagacearum

Genetic diversity in the oak wilt pathogen Ceratocystis fagacearum was assessed using restriction fragment length polymorphisms (RFLPs) of the mitochondrial DNA (mtDNA) and anonymous RFLP loci in the nuclear DNA (nuDNA). No genetic variation was detected in the mtDNA among 27 isolates sampled from a broad geographical area. Southern hybridization to 100 anonymous, random, nuDNA probes detected a low level of variation among nine of the isolates. Only 35 out of 437 probe-enzyme combinations detected RFLPs. Most of the RFLPs appeared to result from insertions and deletions of less than 200 bp. A composite multilocus haplotype based on hybridization to six anonymous probes could differentiate each of the nine isolates tested, suggesting that these probes may be useful for further studies of the population biology and epidemiology of this pathogen. Hypotheses are presented to account for the low level of genetic variation.

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.

The sugar beet mitochondrial genome contains an ORF sharing sequence homology with the gene for the 30 kDa subunit of bovine mitochondrial complex I

From a sugar beet mitochondrial DNA library, we have isolated an open reading frame (ORF192) showing extensive homology to the gene for the 30 kDa subunit of the bovine mitochondrial complex I (NADH: ubiquinone reductase). The ORF192 was found to be actively transcribed to give an RNA of approximately 1.0 kb. We have designated this gene nad9. Transcripts from the nad9 locus are edited by five C to U transitions, increasing similarity with the amino acid sequence of the corresponding bovine and Neurospora crassa polypeptides. Southern blot hybridization also indicates that nad9 is present in the mitochondrial genomes of a variety of higher plant species.

Mitochondrial DNA of Schizophyllum commune: restriction map, genetic map, and mode of inheritance

Mitochondrial DNA (mtDNA) found in the basidiomycete Schizophyllum commune (strain 4-40) is a circular molecule 49.75 kbp in length. A physical map containing 61 restriction sites revealed no repeat structures. Cloned genes from Neurospora crassa, Aspergillus nidulans, and Saccharomyces cerevisiae were used in Southern hybridizations to locate nine mitochondrial genes, including a possible pseudogene of ATPase 9, on the restriction map. A probe from a functional ATPase 9 gene identified homologous fragments only in the nuclear genome of S. commune. Restriction fragment length polymorphisms (RFLPs) between mtDNA isolated from different strains of S. commune were used to show that mitochondria do not migrate with nuclei during dikaryosis.

Investigation of Mitochondrial Transmission in Selected Matings between Homokaryons from Commercial and Wild-Collected Isolates of Agaricus bisporus (= Agaricus brunnescens )

Ten heterokaryons of Agaricus bisporus (= Agaricus brunnescens ) were shown to carry four different mitochondrial (mt) genotypes by analysis of mt restriction fragment length polymorphisms (RFLPs). Fifteen homokaryons derived from these strains were used to investigate mt inheritance in A. bisporus. One hundred eighty-nine pairings were performed in 25 different combinations. Pairings in 15 different combinations produced heterokaryons on the basis of nuclear RFLP analyses and/or fruiting trials. The mt genotype of each new intraspecies hybrid was examined by using mt RFLPs as genetic markers. Our results suggest the following. (i) Recombination between the mt genomes was not a common event. (ii) From most individual pairings, all heterokaryons carried the same mt genotype. (iii) Heterokaryons carrying either of the two possible mt genotypes were observed in certain crosses after modification of the pairing procedure. A biparental transmission pattern was demonstrated for some crosses, but there appears to be a preference for one of the mt genotypes to predominate in any specific pairing.

Involvement of a large inverted repeated sequence in a recombinational rearrangement of the mitochondrial genome of the higher fungus Agrocybe aegerita

Southern hybridization of the total DNA of Agrocybe aegerita with cloned mitochondrial (mt) probes revealed a sequence homology between two distant mitochondrial restriction fragments. From the mtDNA restriction map and the distribution of restriction sites on the cross-hybridizing mitochondrial fragments, two copies of a large inverted repeated sequence (IR) of 3 kbp were located on the mitochondrial genome. These IR sequences divided the 80 kbp mtDNA into two single-copy regions of 24 kbp (SSC) and 50 kbp (LSC). For the first time in higher fungi, this IR sequence has been shown to be involved in an intramolecular homologous recombinational event. Such a rearrangement led to an inversion of the orientation of the two unique-copy regions, without any change in mtDNA complexity. The location of the recombinational event was compared with previously reported plant and fungal mitochondrial rearrangements and the potential role of the IR sequence was discussed.

Transmission of mitochondrial DNA in Ustilago violacea

Mitochondrial DNA (mtDNA) restriction fragment length polymorphisms (RFLPs) were used as genetic markers for following mitochondrial transmission in the basidiomycete Ustilago violacea. Yeast-like cells of opposite mating types (a1 and a2) were mated on 2% water agar and were treated with alpha-tocopherol to induce formation of dikaryotic hyphae. Upon depletion of the alpha-tocopherol, the hyphae budded off haploid cells with parental nuclear genotypes. These cells were examined for mitochondrial RFLP phenotype. In progeny expressing the a1 mating type, mitochondria from either parent were observed equally frequently. In progeny with the a2 mating type, mitochondria were almost exclusively (94%) from the a2 parent.

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

The mtDNA of a wild-type strain of Agrocybe aegerita was purified from mitochondria isolated by cellular fractionation. A representative library was constructed in E. coli by molecular cloning at the HindIII restriction site of pBR322. Southern hybridizations between total DNA of the fungal strain and cloned mitochondrial insert probes were used to establish the restriction map of the mtDNA molecule. Its size was assessed at about 80,500 bp. Four structural genes (for Cox 1, Cox 2, Atp 6, and Atp 8) were located on the map by heterologous hybridizations with oligonucleote probes specific for yeast mitochondrial genes. The location of the genes coding for the large and the small RNAs of the mitochondrial ribosome was determined by hybridization with the E. coli rrnB operon. A comparison of A. aegerita mtDNA organization with that of both phylogenetically close and distant fungi is discussed.

Mitochondrial DNAs and plasmids as taxonomic characteristics in Trichoderma viride

Mitochondrial DNA (mtDNA) was purified from 12 isolates of the Trichoderma viride aggregate and found to be, on the average, 32.7 kb in size. Plasmids were present in the mtDNA preparations from 8 of 12 strains of T. viride examined. Plasmids in four of the strains produced ladderlike banding patterns on gels, and these plasmids were studied in detail. The ladderlike patterns were produced by single molecules that were supercoiled to various degrees. Plasmids from two of the strains do not have homology with the mtDNA but do have a limited amount of homology with each other. No phenotype could be associated with the presence of a plasmid. Restriction endonuclease digestion of the mtDNAs produced patterns in which the presence or absence of certain fragments correlated with the classification of the strains into T. viride group I or II. Phenetic cluster analysis and parsimony analysis of the fragment patterns produced groups that corresponded to T. viride groups I and II. The fragment patterns were very diverse, with nearly all strains having a unique pattern. However, two strains of T. viride group I from widely different geographical locations did have identical restriction patterns for all the enzymes used in this study. This result indicates that it may not be possible to use mtDNA restriction patterns alone to identify Trichoderma strains.

DNA polymorphisms in commercial and wild strains of the cultivated mushroom, Agaricus bisporus

DNA from the cultivated mushroom, Agaricus bisporus, was cloned into the bacteriophage lambda vector EMBL3 creating a partial genomic library. Ten random clones from the library were used to probe for restriction fragment length polymorphisms (RFLPs). Six of the ten probes detected polymorphisms and were used to demonstrate variation in wild and cultivated strains of the mushroom. These results suggest that RFLPs could form a basis for genetic finger-printing and subsequent strain protection in A. bisporus. In single spore progeny, RFLPs were used to demonstrate normal meiotic segregation and to differentiate between homokaryons and heterokaryons. RFLPs therefore have great potential in the development of the genetics and breeding of this commercially important species.

Mitochondrial DNAs of Suillus: three fold size change in molecules that share a common gene order

We constructed restriction-site and gene maps for mitochondrial DNAs from seven isolates of five species of Suillus (Boletaceae, Basidiomycotina). Each mitochondrial genome exists as a single circular chromosome, ranging in size from 36 to 121 kb. Comparisons within species and between two closely related species revealed that insertions and deletions are the major form of genome change, whereas most restriction sites are conserved. Among more distantly related species, size and restriction-site differences were too great to allow precise alignments of maps, but small clusters of putatively homologous restriction sites were found. Two mitochondrial gene orders exist in the five species. These orders differ only by the relative positions of the genes for ATPase subunit 9 and the small ribosomal RNA and are interconvertible by a single transposition. One of the two gene arrangements is shared by four species whose mitochondrial DNAs span the entire size range of 36 to 121 kb. The conservation of gene order in molecules that vary over three-fold in size and share few restriction sites demonstrates a low frequency of rearrangements relative to insertions, deletions, and base substitutions.