Transcription Factor Repertoire of Necrotrophic Fungal Phytopathogen Ascochyta rabiei: Predominance of MYB Transcription Factors As Potential Regulators of Secretome

Systematic identification of CAZymes and transcription factors in the hypercellulolytic fungus Penicillium funiculosum NCIM1228 involved in lignocellulosic biomass degradation

BACKGROUND

Penicillium funiculosum NCIM1228 is a filamentous fungus that was identified in our laboratory to have high cellulolytic activity. Analysis of its secretome suggested that it responds to different carbon substrates by secreting specific enzymes capable of digesting those substrates. This phenomenon indicated the presence of a regulatory system guiding the expression of these hydrolyzing enzymes. Since transcription factors (TFs) are the key players in regulating the expression of enzymes, this study aimed first to identify the complete repertoire of Carbohydrate Active Enzymes (CAZymes) and TFs coded in its genome. The regulation of CAZymes was then analysed by studying the expression pattern of these CAZymes and TFs in different carbon substrates-Avicel (cellulosic substrate), wheat bran (WB; hemicellulosic substrate), Avicel + wheat bran, pre-treated wheat straw (a potential substrate for lignocellulosic ethanol), and glucose (control).

RESULTS

The P. funiculosum NCIM1228 genome was sequenced, and 10,739 genes were identified in its genome. These genes included a total of 298 CAZymes and 451 TF coding genes. A distinct expression pattern of the CAZymes was observed in different carbon substrates tested. Core cellulose hydrolyzing enzymes were highly expressed in the presence of Avicel, while pre-treated wheat straw and Avicel + wheat bran induced a mixture of CAZymes because of their heterogeneous nature. Wheat bran mainly induced hemicellulases, and the least number of CAZymes were expressed in glucose. TFs also exhibited distinct expression patterns in each of the carbon substrates. Though most of these TFs have not been functionally characterized before, homologs of NosA, Fcr1, and ATF21, which have been known to be involved in fruiting body development, protein secretion and stress response, were identified.

CONCLUSIONS

Overall, the P. funiculosum NCIM1228 genome was sequenced, and the CAZymes and TFs present in its genome were annotated. The expression of the CAZymes and TFs in response to various polymeric sugars present in the lignocellulosic biomass was identified. This work thus provides a comprehensive mapping of transcription factors (TFs) involved in regulating the production of biomass hydrolyzing enzymes.

The SbbHLH041-SbEXPA11 Module Enhances Cadmium Accumulation and Rescues Biomass by Increasing Photosynthetic Efficiency in Sorghum

In plants, expansin genes are responsive to heavy metal exposure. To study the bioremediary potential of this important gene family, we discovered a root-expressed expansin gene in sorghum, SbEXPA11, which is notably upregulated following cadmium (Cd) exposure. However, the mechanism underlying the Cd detoxification and accumulation mediated by SbEXPA11 in sorghum remains unclear. We overexpressed SbEXPA11 in sorghum and compared wild-type (WT) and SbEXPA11-overexpressing transgenic sorghum in terms of Cd accumulation and physiological indices following Cd. Compared with the WT, we found that SbEXPA11 mediates Cd tolerance by exerting reactive oxygen species (ROS)-scavenging effects through upregulating the expression of antioxidant enzymes. Moreover, the overexpression of SbEXPA11 rescued biomass production by increasing the photosynthetic efficiency of transgenic plants. In the pot experiment with a dosage of 10 mg/kg Cd, transgenic sorghum plants demonstrated higher efficacy in reducing the Cd content of the soil (8.62 mg/kg) compared to WT sorghum plants (9.51 mg/kg). Subsequent analysis revealed that the SbbHLH041 transcription factor has the ability to induce SbEXPA11 expression through interacting with the E-box located within the SbEXPA11 promoter. These findings suggest that the SbbHLH041-SbEXPA11 cascade module may be beneficial for the development of phytoremediary sorghum varieties.

Ascochyta rabiei: A threat to global chickpea production

The necrotrophic fungus Ascochyta rabiei causes Ascochyta blight (AB) disease in chickpea. A. rabiei infects all aerial parts of the plant, which results in severe yield loss. At present, AB disease occurs in most chickpea-growing countries. Globally increased incidences of A. rabiei infection and the emergence of new aggressive isolates directed the interest of researchers toward understanding the evolution of pathogenic determinants in this fungus. In this review, we summarize the molecular and genetic studies of the pathogen along with approaches that are helping in combating the disease. Possible areas of future research are also suggested.

TAXONOMY

kingdom Mycota, phylum Ascomycota, class Dothideomycetes, subclass Coelomycetes, order Pleosporales, family Didymellaceae, genus Ascochyta, species rabiei.

PRIMARY HOST

A. rabiei survives primarily on Cicer species.

DISEASE SYMPTOMS

A. rabiei infects aboveground parts of the plant including leaves, petioles, stems, pods, and seeds. The disease symptoms first appear as watersoaked lesions on the leaves and stems, which turn brown or dark brown. Early symptoms include small circular necrotic lesions visible on the leaves and oval brown lesions on the stem. At later stages of infection, the lesions may girdle the stem and the region above the girdle falls off. The disease severity increases at the reproductive stage and rounded lesions with concentric rings, due to asexual structures called pycnidia, appear on leaves, stems, and pods. The infected pod becomes blighted and often results in shrivelled and infected seeds.

DISEASE MANAGEMENT STRATEGIES

Crop failures may be avoided by judicious practices of integrated disease management based on the use of resistant or tolerant cultivars and growing chickpea in areas where conditions are least favourable for AB disease development. Use of healthy seeds free of A. rabiei, seed treatments with fungicides, and proper destruction of diseased stubbles can also reduce the fungal inoculum load. Crop rotation with nonhost crops is critical for controlling the disease. Planting moderately resistant cultivars and prudent application of fungicides is also a way to combat AB disease. However, the scarcity of AB-resistant accessions and the continuous evolution of the pathogen challenges the disease management process.

USEFUL WEBSITES

https://www.ndsu.edu/pubweb/pulse-info/resourcespdf/Ascochyta%20blight%20of%20chickpea.pdf https://saskpulse.com/files/newsletters/180531_ascochyta_in_chickpeas-compressed.pdf http://www.pulseaus.com.au/growing-pulses/bmp/chickpea/ascochyta-blight http://agriculture.vic.gov.au/agriculture/pests-diseases-and-weeds/plant-diseases/grains-pulses-and-cereals/ascochyta-blight-of-chickpea http://www.croppro.com.au/crop_disease_manual/ch05s02.php https://www.northernpulse.com/uploads/resources/722/handout-chickpeaascochyta-nov13-2011.pdf http://oar.icrisat.org/184/1/24_2010_IB_no_82_Host_Plant https://www.crop.bayer.com.au/find-crop-solutions/by-pest/diseases/ascochyta-blight.

An MYB Transcription Factor Modulates Panax notoginseng Resistance Against the Root Rot Pathogen Fusarium solani by Regulating the Jasmonate Acid Signaling Pathway and Photosynthesis

Root rot of Panax notoginseng, a precious Chinese medicinal plant, seriously impacts its sustainable production. However, the molecular regulatory mechanisms employed by P. notoginseng against root rot pathogens, including Fusarium solani, are still unclear. In this study, the PnMYB2 gene was isolated, and its expression was affected by independent treatments with four signaling molecules (methyl jasmonate, ethephon, salicylic acid, and hydrogen peroxide) as assessed by quantitative real-time PCR. Moreover, the PnMYB2 expression level was induced by F. solani infection. The PnMYB2 protein localized to the nucleus and may function as a transcription factor. When overexpressed in transgenic tobacco, the PnMYB2 gene conferred resistance to F. solani. Jasmonic acid (JA) metabolism and disease resistance-related genes were induced in the transgenic tobacco, and the JA content significantly increased compared with in the wild type. Additionally, transcriptome sequencing, Kyoto Encyclopedia of Genes and Genomes annotation enrichment, and metabolic pathway analyses of the differentially expressed genes in the transgenic tobacco revealed that JA metabolic, photosynthetic, and defense response-related pathways were activated. In summary, PnMYB2 is an important transcription factor in the defense responses of P. notoginseng against root rot pathogens that acts by regulating JA signaling, photosynthesis, and disease-resistance genes.

Comparative Transcriptome Analysis of Onion in Response to Infection by Alternaria porri (Ellis) Cifferi

Purple blotch (PB) is one of the most destructive foliar diseases of onion and other alliums, caused by a necrotrophic fungal pathogen Alternaria porri. There are no reports on the molecular response of onion to PB infection. To elucidate the response of onion to A. porri infection, we consequently carried out an RNAseq analysis of the resistant (Arka Kalyan; AK) and susceptible (Agrifound rose; AFR) genotype after an artificial infection. Through differential expression analyses between control and pathogen-treated plants, we identified 8,064 upregulated and 248 downregulated genes in AFR, while 832 upregulated and 564 downregulated genes were identified in AK. A further significant reprogramming in the gene expression profile was also demonstrated by a functional annotation analysis. Gene ontology (GO) terms, which are particularly involved in defense responses and signaling, are overrepresented in current analyses such as "oxidoreductase activity," "chitin catabolic processes," and "defense response." Several key plant defense genes were differentially expressed on A. porri infection, which includes pathogenesis-related (PR) proteins, receptor-like kinases, phytohormone signaling, cell-wall integrity, cytochrome P450 monooxygenases, and transcription factors. Some of the genes were exclusively overexpressed in resistant genotype, namely, GABA transporter1, ankyrin repeat domain-containing protein, xyloglucan endotransglucosylase/hydrolase, and PR-5 (thaumatin-like). Antioxidant enzyme activities were observed to be increased after infection in both genotypes but higher activity was found in the resistant genotype, AK. This is the first report of transcriptome profiling in onion in response to PB infection and will serve as a resource for future studies to elucidate the molecular mechanism of onion-A. porri interaction and to improve PB resistance in onions.

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.

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

Numerous studies in plants have shown the vital roles of MYB transcription factors in signal transduction, developmental regulation, biotic/abiotic stress responses and secondary metabolism regulation. However, less is known about the functions of MYBs in Ganoderma. In this study, five medicinal macrofungi of genus Ganoderma were subjected to a genome-wide comparative analysis of MYB genes. A total of 75 MYB genes were identified and classified into four types: 1R-MYBs (52), 2R-MYBs (19), 3R-MYBs (2) and 4R-MYBs (2). Gene structure analysis revealed varying exon numbers (3-14) and intron lengths (7-1058 bp), and noncanonical GC-AG introns were detected in G. lucidum and G. sinense. In a phylogenetic analysis, 69 out of 75 MYB genes were clustered into 15 subgroups, and both single-copy orthologous genes and duplicated genes were identified. The promoters of the MYB genes harbored multiple cis-elements, and specific genes were co-expressed with the G. lucidum MYB genes, indicating the potential roles of these MYB genes in stress response, development and metabolism. This comprehensive and systematic study of MYB family members provides a reference and solid foundation for further functional analysis of MYB genes in Ganoderma species.

Reference Genome Assembly for Australian Ascochyta rabiei Isolate ArME14

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

Assembly, annotation, and comparison of Macrophomina phaseolina isolates from strawberry and other hosts

Background

Macrophomina phaseolina is a fungal plant pathogen with a broad host range, but one genotype was shown to exhibit host preference/specificity on strawberry. This pathogen lacked a high-quality genome assembly and annotation, and little was known about genomic differences among isolates from different hosts.

Results

We used PacBio sequencing and Hi-C scaffolding to provide nearly complete genome assemblies for M. phaseolina isolates representing the strawberry-specific genotype and another genotype recovered from alfalfa. The strawberry isolate had 59 contigs/scaffolds with an N50 of 4.3 Mb. The isolate from alfalfa had an N50 of 5.0 Mb and 14 nuclear contigs with half including telomeres. Both genomes were annotated with MAKER using transcript evidence generated in this study with over 13,000 protein-coding genes predicted. Unique groups of genes for each isolate were identified when compared to closely related fungal species. Structural comparisons between the isolates reveal large-scale rearrangements including chromosomal inversions and translocations. To include isolates representing a range of pathogen genotypes, an additional 30 isolates were sequenced with Illumina, assembled, and compared to the strawberry genotype assembly. Within the limits of comparing Illumina and PacBio assemblies, no conserved structural rearrangements were identified among the isolates from the strawberry genotype compared to those from other hosts, but some candidate genes were identified that were largely present in isolates of the strawberry genotype and absent in other genotypes.

Conclusions

High-quality reference genomes of M. phaseolina have allowed for the identification of structural changes associated with a genotype that has a host preference toward strawberry and will enable future comparative genomics studies. Having more complete assemblies allows for structural rearrangements to be more fully assessed and ensures a greater representation of all the genes. Work with Illumina data from additional isolates suggests that some genes are predominately present in isolates of the strawberry genotype, but additional work is needed to confirm the role of these genes in pathogenesis. Additional work is also needed to complete the scaffolding of smaller contigs identified in the strawberry genotype assembly and to determine if unique genes in the strawberry genotype play a role in pathogenicity.

Genome-Wide Characterization and Expression Analyses of Pleurotus ostreatus MYB Transcription Factors during Developmental Stages and under Heat Stress Based on de novo Sequenced Genome

Pleurotus ostreatus is a commercially grown mushroom species in China. However, studies on the mechanisms of the fruiting body development and stress response of P. ostreatus are still at a primary stage. In this study, we report the entire genome sequence of P. ostreatus CCMSSC03989. Then, we performed comprehensive genome-wide characterization and expression analysis of the MYB transcription factor family during a series of developmental stages and under the condition of heat stress. A 34.76 Mb genome was obtained through next-generation sequencing (NGS) and Bionano optical mapping approaches. The genome has a scaffold N50 of 1.1 Mb and contains 10.11% repeats, and 10,936 gene models were predicted. A total of 20 MYB genes (PoMYB) were identified across the genome, and the full-length open reading frames were isolated. The PoMYBs were classified into 1 repeat (1R), 2R, and 3R-MYB groups according to their MYB domain repeat numbers, and 3R-MYBs possessed relatively more introns than 1R and 2R-MYBs. Based on phylogenetic analysis, the PoMYBs were divided into four groups and showed close relationships with the MYB genes of plants and fungi. RNA-sequencing (RNA-Seq) and quantitative PCR (qPCR) analyses revealed that PoMYB expression showed stage-specific patterns in reproductive stages and could be induced by heat stress. The P. ostreatus draft genome will promote genome-wide analysis, and our study of PoMYBs will promote further functional analysis of MYB genes in mushrooms.