Genome-Wide Characterization and Comparative Analysis of MYB Transcription Factors in Ganoderma Species

Genome-wide identification and expression analysis revealed key transcription factors as potential regulators of high-temperature adaptation of Coriolopsis trogii

Transcription factors (TFs) play a crucial role in gene expression, and studying them can lay the foundation for future research on the functional characterization of TFs involved in various biological processes. In this study, we conducted a genome-wide identification and analysis of TFs in the thermotolerant basidiomycete fungus, Coriolopsis trogii. The TF repertoire of C. trogii consisted of 350 TFs, with C2H2 and Zn2C6 being the largest TF families. When the mycelia of C. trogii were cultured on PDA and transferred from 25 to 35 °C, 14 TFs were up-regulated and 14 TFs were down-regulated. By analyzing RNA-seq data from mycelia cultured at different temperatures and under different carbon sources, we identified 22 TFs that were differentially expressed in more than three comparisons. Co-expression analysis revealed that seven differentially expressed TFs, including four Zn2C6s, one Hap4_Hap_bind, one HMG_box, and one Zinc_knuckle, showed significant correlation with 729 targeted genes. Overall, this study provides a comprehensive characterization of the TF family and systematically screens TFs involved in the high-temperature adaptation of C. trogii, laying the groundwork for further research into the specific roles of TFs in the heat tolerance mechanisms of filamentous fungi.

The regulatory pathway of transcription factor MYB36 from Trichoderma asperellum Tas653 resistant to poplar leaf blight pathogen Alternaria alternata Aal004

In fungi, MYB transcription factors (TFs) mainly regulate growth, development, and resistance to stress. However, as major disease-resistance TFs, they have rarely been studied in biocontrol fungi. In this study, MYB36 of Trichoderma asperellum Tas653 (Ta) was shown to respond strongly to the stress caused by Alternaria alternata Aa1004. Compared with wild-type Ta (Ta-Wt), the inhibition rate of the MYB36 knockout strain (Ta-Kn) on Aa1004 decreased by 11.06%; the superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities decreased by 82.15 U/g, 0.19 OD470/min/g, and 1631.2 μmol/min/g, respectively. The MYB36 overexpression strain (Ta-Oe) not only enhanced hyperparasitism on Aa1004, caused its hyphae to swell, deform, or even rupture, but also reduced the incidence rate of poplar leaf blight. MYB36 regulates downstream (TFs, detoxification genes, defense genes, and other antifungal-related genes by binding to the cis-acting elements "ACAT" and "ATCG". Zinc finger TFs, as the main antifungal TFs, account for 90% of the total TFs, and Zn37.5 (23.24-) and Zn83.7 (23.18-fold) showed the greatest expression difference when regulated directly by MYB36. The detoxification genes mainly comprised 11 major major facilitator superfamily (MFS) genes, among which MYB36 directly increased the expression levels of three genes by more than 2-3.44-fold. The defense genes mainly encoded cytochrome P450 (P450) and hydrolases. e.g., P45061.3 (2-10.95-), P45060.2 (2-7.07-), and Hyd44.6 (2-2.30-fold). This study revealed the molecular mechanism of MYB36 regulation of the resistance of T. asperellum to A. alternata and provides theoretical guidance for the biocontrol of poplar leaf blight and the anti-disease mechanism of biocontrol fungi.

Genome-wide characterization of regulator of chromosome condensation 1 (RCC1) gene family in Artemisia annua L. revealed a conservation evolutionary pattern

BACKGROUND

Artemisia annua is the major source for artemisinin production. The artemisinin content in A. annua is affected by different types of light especially the UV light. UVR8, a member of RCC1 gene family was found to be the UV-B receptor in plants. The gene structures, evolutionary history and expression profile of UVR8 or RCC1 genes remain undiscovered in A. annua.

RESULTS

Twenty-two RCC1 genes (AaRCC1) were identified in each haplotype genome of two diploid strains of A. annua, LQ-9 and HAN1. Varied gene structures and sequences among paralogs were observed. The divergence of most RCC1 genes occurred at 46.7 - 51 MYA which overlapped with species divergence of core Asteraceae during the Eocene, while no recent novel RCC1 members were found in A. annua genome. The number of RCC1 genes remained stable among eudicots and RCC1 genes underwent purifying selection. The expression profile of AaRCC1 is analogous to that of Arabidopsis thaliana (AtRCC1) when responding to environmental stress.

CONCLUSIONS

This study provided a comprehensive characterization of the AaRCC1 gene family and suggested that RCC1 genes were conserved in gene number, structures, constitution of amino acids and expression profiles among eudicots.

Transcription factors: switches for regulating growth and development in macrofungi

Macrofungi (or mushrooms) act as an extraordinarily important part to human health due to their nutritional and/or medicinal value, but the detailed researches in growth and development mechanisms have yet to be explored further. Transcription factors (TFs) play indispensable roles in signal transduction and affect growth, development, and metabolism of macrofungi. In recent years, increasing research effort has been employed to probe the relationship between the development of macrofungi and TFs. Herein, the present review comprehensively summarized the functional TFs researched in macrofungi, including modulating mycelial growth, fructification, sclerotial formation, sexual reproduction, spore formation, and secondary metabolism. Meanwhile, the possible effect mechanisms of TFs on the growth and development of some macrofungi were also revealed. Specific examples of functional characterizations of TFs in macrofungi (such as Schizophyllum commune and Coprinopsis cinerea) were described to a better comprehension of regulatory effect. Future research prospects in the field of TFs of macrofungi are discussed. We illustrated the functional versatility of the TFs in macrofungi based on specific examples. A systematical realization of the interaction and possible mechanisms between TFs and macrofungi can supply possible solutions to regulate genetic characteristics, which supply novel insights into the regulation of growth, development and metabolism of macrofungi. KEY POINTS: • The functional TFs researched in macrofungi were summarized. • The possible effect mechanisms of TFs in macrofungal were described. • The multiple physiological functions of TFs in macrofungi were discussed.

Function of Transcription Factors PoMYB12, PoMYB15, and PoMYB20 in Heat Stress and Growth of Pleurotus ostreatus

MYB transcription factors (TFs) have been extensively studied in plant abiotic stress responses and growth and development. However, the role of MYB TFs in the heat stress response and growth and development of Pleurotus ostreatus remains unclear. To investigate the function of PoMYB12, PoMYB15, and PoMYB20 TFs in P. ostreatus, mutant strains of PoMYB12, PoMYB15, and PoMYB20 were generated using RNA interference (RNAi) and overexpression (OE) techniques. The results indicated that the mycelia of OE-PoMYB12, OE-PoMYB20, and RNAi-PoMYB15 mutant strains exhibited positive effects under heat stress at 32 °C, 36 °C, and 40 °C. Compared to wild-type strains, the OE-PoMYB12, OE-PoMYB20, and RNAi-PoMYB15 mutant strains promoted the growth and development of P. ostreatus. These mutant strains also facilitated the recovery of growth and development of P. ostreatus after 24 h of 36 °C heat stress. In conclusion, the expression of PoMYB12 and PoMYB20 supports the mycelium's response to heat stress and enhances the growth and development of P. ostreatus, whereas PoMYB15 produces the opposite effect.

Comprehensive analysis of metabolome and transcriptome reveals the mechanism of color formation in different leave of Loropetalum Chinense var. Rubrum

BACKGROUND

Loropetalum chinense var. rubrum (L. chinense var. rubrum) is a precious, coloured-leaf native ornamental plant in the Hunan Province. We found an L. chinense var. rubrum tree with three different leaf colours: GL (green leaf), ML (mosaic leaf), and PL (purple leaf). The mechanism of leaf coloration in this plant is still unclear. Therefore, this study aimed to identify the metabolites and genes involved in determining the colour composition of L. chinense var. rubrum leaves, using phenotypic/anatomic observations, pigment content detection, and comparative metabolomics and transcriptomics.

RESULTS

We observed that the mesophyll cells in PL were purple, while those in GL were green and those in ML were a mix of purple-green. The contents of chlorophyll a, b, carotenoids, and total chlorophyll in PL and ML were significantly lower than those in GL. While the anthocyanin content in PL and ML was significantly higher than that in GL. The metabolomics results showed the differences in the content of cyanidin 3-O-glucoside, delphinidin 3-O-glucoside, cyanidin 3,5-O-diglucoside, pelargonidin, and petunidin 3,5-diglucoside in ML, GL, and PL were significant. Considering that the change trend of anthocyanin content change was consistent with the leaf colour difference, we speculated that these compounds might influence the colour of L. chinense var. rubrum leaves. Using transcriptomics, we finally identified nine differentially expressed structural genes (one ANR (ANR1217); four CYP75As (CYP75A1815, CYP75A2846, CYP75A2909, and CYP75A1716); four UFGTs (UFGT1876, UFGT1649, UFGT1839, and UFGT3273) and nine transcription factors (two MYBs (MYB1057 and MYB1211), one MADS-box (MADS1235), two AP2-likes (AP2-like1779 and AP2-like2234), one bZIP (bZIP3720), two WD40s (WD2173 and WD1867) and one bHLH (bHLH1631) that might be related to flavonoid biosynthesis and then impacted the appearance of colour in L. chinense var. rubrum leaves.

CONCLUSION

This study revealed potential molecular mechanisms associated with leaf coloration in L. chinense var. rubrum by analyzing differential metabolites and genes related to the anthocyanin biosynthesis pathway. It also provided a reference for research on leaf colour variation in other ornamental plants.

Transcription factor FfMYB15 regulates the expression of cellulase gene FfCEL6B during mycelial growth of Flammulina filiformis

Background

Cellulose degradation can determine mycelial growth rate and affect yield during the growth of Flammulina filiformis. The degradation of cellulose requires the joint action of a variety of cellulases, and some cellulase-related genes have been detected in mushrooms. However, little is known about the transcriptional regulatory mechanisms of cellulose degradation.

Results

In this study, FfMYB15 that may regulate the expression of cellulase gene FfCEL6B in F. filiformis was identified. RNA interference (RNAi) showed that FfCEL6B positively regulated mycelial growth. Gene expression analyses indicated that the expression patterns of FfCEL6B and FfMYB15 in mycelia cultured on the 0.9% cellulose medium for different times were similar with a correlation coefficient of 0.953. Subcellular localization and transcriptional activity analyses implied that FfMYB15 was located in the nucleus and was a transcriptional activator. Electrophoretic mobility shift assay (EMSA) and dual-luciferase assays demonstrated that FfMYB15 could bind and activate FfCEL6B promoter by recognizing MYB cis-acting element.

Conclusions

This study indicated that FfCEL6B played an active role in mycelial growth of F. filiformis and was regulated by FfMYB15.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01932-z.

Gllac7 Is Induced by Agricultural and Forestry Residues and Exhibits Allelic Expression Bias in Ganoderma lucidum

Ganoderma lucidum has a wide carbon spectrum, while the expression profile of key genes relevant to carbon metabolism on different carbon sources has been seldom studied. Here, the transcriptomes of G. lucidum mycelia cultured on each of 19 carbon sources were conducted. In comparison with glucose, 16 to 1,006 genes were upregulated and 7 to 1,865 genes were downregulated. Significant gene expression dynamics and induced activity were observed in laccase genes when using agricultural and forestry residues (AFRs) as solo carbon sources. Furthermore, study of laccase gene family in two haploids of G. lucidum GL0102 was conducted. Totally, 15 and 16 laccase genes were identified in GL0102_53 and GL0102_8, respectively, among which 15 pairs were allelic genes. Gene structures were conserved between allelic laccase genes, while sequence variations (most were SNPs) existed. Nine laccase genes rarely expressed on all the tested carbon sources, while the other seven genes showed high expression level on AFRs, especially Gllac2 and Gllac7, which showed 5- to 1,149-fold and 4- to 94-fold upregulation in mycelia cultured for 5 days, respectively. The expression of H53lac7 was consistently higher than that of H8lac7_1 on all the carbon sources except XM, exhibiting a case of allelic expression bias. A total of 47 SNPs and 3 insertions/deletions were observed between promoters of H53lac7 and H8lac7_1, which lead to differences in predicted binding sites of zinc fingers. These results provide scientific data for understanding the gene expression profile and regulatory role on different carbon sources and may support further functional research of laccase.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

Independent Evolution of the MYB Family in Brown Algae

Myeloblastosis (MYB) proteins represent one of the largest families of eukaryotic transcription factors and regulate important processes in growth and development. Studies on MYBs have mainly focused on animals and plants; however, comprehensive analysis across other supergroups such as SAR (stramenopiles, alveolates, and rhizarians) is lacking. This study characterized the structure, evolution, and expression of MYBs in four brown algae, which comprise the biggest multicellular lineage of SAR. Subfamily 1R-MYB comprised heterogeneous proteins, with fewer conserved motifs found outside the MYB domain. Unlike the SHAQKY subgroup of plant 1R-MYB, THAQKY comprised the largest subgroup of brown algal 1R-MYBs. Unlike the expansion of 2R-MYBs in plants, brown algae harbored more 3R-MYBs than 2R-MYBs. At least ten 2R-MYBs, fifteen 3R-MYBs, and one 6R-MYB orthologs existed in the common ancestor of brown algae. Phylogenetic analysis showed that brown algal MYBs had ancient origins and a diverged evolution. They showed strong affinity with stramenopile species, while not with red algae, green algae, or animals, suggesting that brown algal MYBs did not come from the secondary endosymbiosis of red and green plastids. Sequence comparison among all repeats of the three types of MYB subfamilies revealed that the repeat of 1R-MYBs showed higher sequence identity with the R3 of 2R-MYBs and 3R-MYBs, which supports the idea that 1R-MYB was derived from loss of the first and second repeats of the ancestor MYB. Compared with other species of SAR, brown algal MYB proteins exhibited a higher proportion of intrinsic disordered regions, which might contribute to multicellular evolution. Expression analysis showed that many MYB genes are responsive to different stress conditions and developmental stages. The evolution and expression analyses provided a comprehensive analysis of the phylogeny and functions of MYBs in brown algae.

Haploid Genome Analysis Reveals a Tandem Cluster of Four HSP20 Genes Involved in the High-Temperature Adaptation of Coriolopsis trogii

Coriolopsis trogii is a typical thermotolerant basidiomycete fungus, but its thermotolerance mechanisms are currently unknown. In this study, two monokaryons of C. trogii strain Ct001 were assembled: Ct001_29 had a genome assembly size of 38.85 Mb and encoded 13,113 genes, while Ct001_31 was 40.19 Mb in length and encoded 13,309 genes. Comparative intra- and interstrain genomic analysis revealed the rich genetic diversity of C. trogii, which included more than 315,194 single-nucleotide polymorphisms (SNPs), 30,387 insertion/deletions (indels), and 1,460 structural variations. Gene family analysis showed that the expanded families of C. trogii were functionally enriched in lignocellulose degradation activities. Furthermore, a total of 14 allelic pairs of heat shock protein 20 (HSP20) genes were identified in the C. trogii genome. The expression profile obtained from RNA sequencing (RNA-Seq) showed that four tandem-duplicated allelic pairs, HSP20.5 to HSP20.8, had more than 5-fold higher expression at 35°C than at 25°C. In particular, HSP20.5 and HSP20.8 were the most highly expressed HSP20 genes. Allelic expression bias was found for HSP20.5 and HSP20.8; the expression of Ct29HSP20.8 was at least 1.34-fold higher than that of Ct31HSP20.8, and that of Ct31HSP20.5 was at least 1.5-fold higher than that of Ct29HSP20.5. The unique structural and expression profiles of the HSP20 genes revealed by these haplotype-resolved genomes provide insight into the molecular mechanisms of high-temperature adaptation in C. trogii. IMPORTANCE Heat stress is one of the most frequently encountered environmental stresses for most mushroom-forming fungi. Currently available fungal genomes are mostly haploid because high heterozygosity hinders diploid genome assembly. Here, two haplotype genomes of C. trogii, a thermotolerant basidiomycete, were assembled separately. A conserved tandem cluster of four HSP20 genes showing allele-specific expression was found to be closely related to high-temperature adaptation in C. trogii. The obtained haploid genomes and their comparison offer a more thorough understanding of the genetic background of C. trogii. In addition, the responses of HSP20 genes at 35°C, which may contribute to the growth and survival of C. trogii at high temperatures, could inform the selection and breeding of elite strains in the future.