Variable Number Tandem Repeats in the Mitochondrial DNA of Lentinula edodes

Structural Analysis of the A Mating Type Locus and Development of the Mating Type Marker of Agaricus bisporus var. bisporus

Karyotyping in Agaricus bisporus is crucial for both the isolation of homokaryotic strains and the confirmation of dikaryon establishment. For the verification of the karyotype, the A mating type loci of two homokaryotic strains, H39 and H97, were analyzed through comparative sequence analysis. The two loci showed major differences in two sequence regions designated as Region 1 and Region 2. H97 had a putative DNA transposon in Region 1 that had target site duplications (TSDs), terminal inverted repeats (TIRs), and a loop sequence, in contrast to H39, which only had the insertional target sequence. Homologous sequences of the transposon were discovered in the two different chromosomes of H97 and in one of H39, all of which have different TSDs but share high sequence homology in TIR. Region 2 shared three consensus sequences between H97 and H39. However, it was only from H97 that a large insertional sequence of unknown origin was discovered between the first and second consensus sequences. The difference in length in Region 1, employed for the verification of the A mating type, resulted in the successful verification of mating types in the heterokaryotic and homokaryotic strains. This length difference enables the discrimination between homo- and heterokaryotic spores by PCR. The present study suggests that the A mating type locus in A. bisporus H97 has evolved through transposon insertion, allowing the discrimination of the mating type, and thus the nuclear type, between A. bisporus H97 and H39.

Comparative structural analysis on the mitochondrial DNAs from various strains of Lentinula edodes

The evolution of mitochondria through variations in mitochondrial DNA (mtDNA) is one of the intriguing questions in eukaryotic cells. In order to assess the causes of the variations in mitochondria, the mtDNAs of the 21 strains of Lentinula edodes were assembled for this study, and analyzed together with four published mtDNA sequences. The mtDNAs were within the sizes of 117 kb ~ 122 kb. The gene number was observed consistent except for two mtDNAs, which carry a duplicated trnG1-trnG2 unit or a putative gene deletion. The size variation was largely attributed to the number of introns, repeated sequences, transposable elements (TEs), and plasmid-related sequences. Intron loss and gain were found from cox1, rnl, and rns of three mtDNAs. Loss of two introns in cox1 of KY217797.1 reduced its size by 2.7 kb, making it the smallest cox1 gene (8.4 kb) among the cox1s of the 25 mtDNAs, whereas gain of a Group II intron (2.65 kb) and loss of a Group I intron (1.7 kb) in cox1 of MF774813.1 resulted in the longest cox1 (12 kb). In rnl of L. edodes, we discovered four intron insertion consensus sequences which were unique to basidiomycetes but not ascomycetes. Differential incorporation of introns was the primary cause of the rnl size polymorphism. Homing endonucleases (HEGs) were suggestively involved in the mobilization of the introns because all of the introns have HEG genes of the LAGRIDADG or GIY-YIG families with the conserved HEG cleavage sites. TEs contributed to 11.04% of the mtDNA size in average, of which 7.08% was LTR-retrotransposon and 3.96% was DNA transposon, whereas the repeated sequences covered 4.6% of the mtDNA. The repeat numbers were variable in a strain-dependent manner. Both the TEs and repeated sequences were mostly found in the intronic and intergenic regions. Lastly, two major deletions were found in the plasmid-related sequence regions (pol2-pol3 and pol1-atp8) in the five mtDNAs. Particularly, the 6.8 kb-long deletion at pol2-pol3 region made MF774813.1 the shortest mtDNA of all. Our results demonstrate that mtDNA is a dynamic molecule that persistently evolves over a short period of time by insertion/deletion and repetition of DNA segments at the strain level.

Mitochondrial Effects on the Physiological Characteristics of Lentinula edodes

In the mating of filamentous basidiomycetes, dikaryotic mycelia are generated through the reciprocal movement of nuclei to a monokaryotic cytoplasm where a nucleus of compatible mating type resides, resulting in the establishment of two different dikaryotic strains having the same nuclei but different mitochondria. To better understand the role of mitochondria in mushrooms, we created four sets of dikaryotic strains of Lentinula edodes, including B2 × E13 (B2 side) and B2 × E13 (E13 side), B5 × E13 (B5 side) and B5 × E13 (E13 side), E8 × H3 (E8 side) and E8 × H3 (H3 side), and K3 × H3 (K3 side) and K3 × H3 (H3 side). The karyotypes and mitochondrial types of the dikaryotic strains were successfully identified by the A mating type markers and the mitochondrial variable length tandem repeat markers, respectively. Comparative analyses of the dikaryotic strains on the mycelial growth, substrate browning, fruiting characteristics, and mitochondrial gene expression revealed that certain mitochondria are more effective in the mycelial growth and the production of fruiting body, possibly through the activated energy metabolism. Our findings indicate that mitochondria affect the physiology of dikaryotic strains having the same nuclear information and therefore a selection strategy aimed at mitochondrial function is needed in the development of new mushroom strain.

Mitochondrial Genomics of Six Cacao Pathogens From the Basidiomycete Family Marasmiaceae

Thread blight disease has recently been described as an emerging disease on cacao (Theobroma cacao) in Ghana. In Ghana, thread blight disease is caused by multiple species of the Marasmiaceae family: Marasmius tenuissimus, M. crinis-equi, M. palmivorus, and Marasmiellus scandens. Interestingly, two additional members of the Marasmiaceae; Moniliophthora roreri (frosty pod rot) and Moniliophthora perniciosa (witches’ broom disease), are major pathogens of cacao in the Western hemisphere. It is important to accurately characterize the genetic relationships among these economically important species in support of their disease management. We used data from Illumina NGS-based genome sequencing efforts to study the mitochondrial genomes (mitogenomes) of the four cacao thread blight associated pathogens from Ghana and compared them with published mitogenomes of Mon. roreri and Mon. perniciosa. There is a remarkable interspecies variation in mitogenome size within the six cacao-associated Marasmiaceae species, ranging from 43,121 to 109,103 bp. The differences in genome lengths are primarily due to the number and lengths of introns, differences in intergenic space, and differences in the size and numbers of unidentified ORFs (uORF). Among seven M. tenuissimus mitogenomes sequenced, there is variation in size and sequence pointing to divergent evolution patterns within the species. The intronic regions show a high degree of sequence variation compared to the conserved sequences of the 14 core genes. The intronic ORFs identified, regardless of species, encode GIY-YIG or LAGLIDADG domain-containing homing endonuclease genes. Phylogenetic relationships using the 14 core proteins largely mimic the phylogenetic relationships observed in gene order patterns, grouping M. tenuissimus with M. crinis-equi, and M. palmivorus with Mon. roreri and Mon. perniciosa, leaving Mar. scandens as an outlier. The results from this study provide evidence of independent expansion/contraction events and sequence diversification in each species and establish a foundation for further exploration of the evolutionary trajectory of the fungi in Marasmiaceae family.

Molecular differences in mitochondrial DNA genomes of dogs with malignant mammary tumours

The aim of this study was to determine molecular defects in mitochondrial DNA (mtDNA) with the use of large-scale genome analysis in malignant canine mammary gland tumours and indicate whether these changes were linked with the carcinogenesis process. With the use of the NGS technology, we sequenced 27 samples of mtDNA isolated from blood and tumours obtained from 13 dogs with mammary gland tumours. The total number of mutations and polymorphisms in the analysed mitochondrial genomes was 557. We identified 383 single nucleotide polymorphisms (SNP), 32 indels (or length polymorphisms), 4 mutations, 137 heteroplasmic positions and 1 indel mutation. The highest variability (132 changes) was observed in the variable number of tandem repeats (VNTR) region. The heteroplasmy rate in VNTR varied among individuals and even between two tumours in one organism. Our previous study resulted in determination of a probable CpG island in this region, thus it is not excluded that these changes might alter mtDNA methylation. Only the ATP8 gene was not affected by any polymorphisms or mutations, whereas the COX1 gene had the highest number of polymorphisms from all protein-coding genes. One change m.13594G>A was detected in a region spanning two genes: ND5 and ND6, from which a deleterious effect was observed for the ND5 protein. Molecular changes were frequently observed in the TΨC loop, which is thought to interact with ribosomal RNA.