Draft genome sequencing and secretome analysis of fungal phytopathogen Ascochyta rabiei provides insight into the necrotrophic effector repertoire

Secretome Analysis of Thermothelomyces thermophilus LMBC 162 Cultivated with Tamarindus indica Seeds Reveals CAZymes for Degradation of Lignocellulosic Biomass

The analysis of the secretome allows us to identify the proteins, especially carbohydrate-active enzymes (CAZymes), secreted by different microorganisms cultivated under different conditions. The CAZymes are divided into five classes containing different protein families. Thermothelomyces thermophilus is a thermophilic ascomycete, a source of many glycoside hydrolases and oxidative enzymes that aid in the breakdown of lignocellulosic materials. The secretome analysis of T. thermophilus LMBC 162 cultivated with submerged fermentation using tamarind seeds as a carbon source revealed 79 proteins distributed between the five diverse classes of CAZymes: 5.55% auxiliary activity (AAs); 2.58% carbohydrate esterases (CEs); 20.58% polysaccharide lyases (PLs); and 71.29% glycoside hydrolases (GHs). In the identified GH families, 54.97% are cellulolytic, 16.27% are hemicellulolytic, and 0.05 are classified as other. Furthermore, 48.74% of CAZymes have carbohydrate-binding modules (CBMs). Observing the relative abundance, it is possible to state that only thirteen proteins comprise 92.19% of the identified proteins secreted and are probably the main proteins responsible for the efficient degradation of the bulk of the biomass: cellulose, hemicellulose, and pectin.

Genetic Analysis of Partially Resistant and Susceptible Chickpea Cultivars in Response to Ascochyta rabiei Infection

The molecular mechanism involved in chickpea (Cicer arietinum L.) resistance to the necrotrophic fungal pathogen Ascochyta rabiei is not well documented. A. rabiei infection can cause severe damage in chickpea, resulting in significant economic losses. Understanding the resistance mechanism against ascochyta blight can help to define strategies to develop resistant cultivars. In this study, differentially expressed genes from two partially resistant cultivars (CDC Corinne and CDC Luna) and a susceptible cultivar (ICCV 96029) to ascochyta blight were identified in the early stages (24, 48 and 72 h) of A. rabiei infection using RNA-seq. Altogether, 3073 genes were differentially expressed in response to A. rabiei infection across different time points and cultivars. A larger number of differentially expressed genes (DEGs) were found in CDC Corinne and CDC Luna than in ICCV 96029. Various transcription factors including ERF, WRKY, bHLH and MYB were differentially expressed in response to A. rabiei infection. Genes involved in pathogen detection and immune signalings such as receptor-like kinases (RLKs), Leucine-Rich Repeat (LRR)-RLKs, and genes associated with the post-infection defence response were differentially expressed among the cultivars. GO functional enrichment and pathway analysis of the DEGs suggested that the biological processes such as metabolic process, response to stimulus and catalytic activity were overrepresented in both resistant and susceptible chickpea cultivars. The expression patterns of eight randomly selected genes revealed by RNA-seq were confirmed by quantitative PCR (qPCR) analysis. The results provide insights into the complex molecular mechanism of the chickpea defence in response to the A. rabiei infection.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

The nuclear effector ArPEC25 from the necrotrophic fungus Ascochyta rabiei targets the chickpea transcription factor CaβLIM1a and negatively modulates lignin biosynthesis, increasing host susceptibility

Fungal pathogens deploy a barrage of secreted effectors to subvert host immunity, often by evading, disrupting, or altering key components of transcription, defense signaling, and metabolic pathways. However, the underlying mechanisms of effectors and their host targets are largely unexplored in necrotrophic fungal pathogens. Here, we describe the effector protein Ascochyta rabiei PEXEL-like Effector Candidate 25 (ArPEC25), which is secreted by the necrotroph A. rabiei, the causal agent of Ascochyta blight disease in chickpea (Cicer arietinum), and is indispensable for virulence. After entering host cells, ArPEC25 localizes to the nucleus and targets the host LIM transcription factor CaβLIM1a. CaβLIM1a is a transcriptional regulator of CaPAL1, which encodes phenylalanine ammonia lyase (PAL), the regulatory, gatekeeping enzyme of the phenylpropanoid pathway. ArPEC25 inhibits the transactivation of CaβLIM1a by interfering with its DNA-binding ability, resulting in negative regulation of the phenylpropanoid pathway and decreased levels of intermediates of lignin biosynthesis, thereby suppressing lignin production. Our findings illustrate the role of fungal effectors in enhancing virulence by targeting a key defense pathway that leads to the biosynthesis of various secondary metabolites and antifungal compounds. This study provides a template for the study of less explored necrotrophic effectors and their host target functions.

Whole genome sequencing of Wilsonomyces carpophilus, an incitant of shot hole disease in stone fruits: insights into secreted proteins of a necrotrophic fungal repository

BACKGROUND

Shot hole is one of the important fungal diseases in stone fruits viz., peach, plum, apricot and cherry caused by Wilsonomyces carpophilus and almond among nut crops. Fungicides significantly decrease the disease. Pathogenicity studies proved a wide host range of the pathogen infecting all stone fruits and almond among the nut crops, however, the mechanism underlying host-pathogen interaction is still unknown. Molecular detection of the pathogen using polymerase chain reaction (PCR) based simple sequence repeat (SSR) markers is also unknown due to the unavailability of the pathogen genome.

METHODS AND RESULTS

We examined the morphology, pathology and genomics of the Wilsonomyces carpophilus. Whole genome sequencing of the W. carpophilus was carried out by Illumina HiSeq and PacBio high throughput sequencing plate-forms through hybrid assembly. Constant selection pressure alters the molecular mechanism of the pathogen causing disease. The studies revealed that the necrotrophs are more lethal with a complex pathogenicity mechanism and little-understood effector repositories. The different isolates of necrotrophic fungus W. carpophilus causing shot hole in stone fruits namely peach, plum, apricot and cherry, and almonds among the nut crops showed a significant variation in their morphology, however, the probability value (p = 0.29) suggests in-significant difference in the pathogenicity. Here, we reported draft genome of W. carpophilus of size 29.9 Mb (Accession number: PRJNA791904). A total of 10,901 protein-coding genes were predicted, including heterokaryon incompatibility genes, cytochrome-p450 genes, kinases, sugar transporters among others. We found 2851 simple sequence repeats (SSRs), tRNAs, rRNAs and pseudogenes in the genome. The most prominent proteins showing necrotrophic lifestyle of the pathogen were hydrolases, polysaccharide-degrading enzymes, esterolytic, lipolytic, and proteolytic enzymes accounted for 225 released proteins. Among the 223 fungal species, top-hit species distribution revealed the majority of hits against the Pyrenochaeta species followed by Ascochyta rabiei and Alternaria alternata.

CONCLUSION

Draft genome of W. carpophilus is 29.9 Mb based on Illumina HiSeq and PacBio hybrid assembly. The necrotrophs are more lethal with a complex pathogenicity mechanism. A significant variation in morphology was observed in different pathogen isolates. A total of 10,901 protein-coding genes were predicted in the pathogen genome including heterokaryon incompatibility, cytochrome-p450 genes, kinases and sugar transporters. We found 2851 SSRs, tRNAs, rRNAs and pseudogenes, and prominent proteins showing necrotrophic lifestyle such as hydrolases, polysaccharide-degrading enzymes, esterolytic, lipolytic and proteolytic enzymes. The top-hit species distribution were against the Pyrenochaeta spp. followed by Ascochyta rabiei.

Ascochyta Blight in Chickpea: An Update

Chickpea (Cicer arietinum L.), one of the most cultivated legumes worldwide, is crucial for the economy of several countries and a valuable source of nutrients. Yields may be severely affected by Ascochyta blight, a disease caused by the fungus Ascochyta rabiei. Molecular and pathological studies have not yet managed to establish its pathogenesis, since it is highly variable. Similarly, much remains to be elucidated about plant defense mechanisms against the pathogen. Further knowledge of these two aspects is fundamental for the development of tools and strategies to protect the crop. This review summarizes up-to-date information on the disease's pathogenesis, symptomatology, and geographical distribution, as well as on the environmental factors that favor infection, host defense mechanisms, and resistant chickpea genotypes. It also outlines existing practices for integrated blight management.

BAR domain is essential for early endosomal trafficking and dynamics in Ascochyta rabiei

Ascochyta blight disease is a devastating disease caused by the fungal pathogen Ascochyta rabiei that threatens chickpea production around the globe. Endocytic mechanism has a significant role in fungal growth and virulence. The underlying biology of biogenesis of central component of endocytosis viz Rab5 vesicles, is not completely understood. The involvement of F-BAR domain containing protein (ArF-BAR) in various cellular processes that collectively make ArF-BAR as an important virulence determinant. Here, we report that ArF-BAR is involved in biogenesis and motility of early endosome. In the absence of ArF-BAR gene (Δarf-bar), fungal mutants exhibited reduced number of EGFP coated ArRab5 vesicles, along with the considerable reduction in their dynamics. Here, we show that ArF-BAR interacts with clathrin light chain (ArCLC), specifically with its F-BAR domain. These findings suggests the novel role of ArF-BAR in biogenesis and dynamics of early endosome. Additionally, ArF-BAR is involved in clathrin-mediated mechanism of endocytosis which is required for host infection and disease development. Identification of this pathway offers new impending targets for disease intervention in plants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03451-5.

A Global Transcriptome and Co-expression Analysis Reveals Robust Host Defense Pathway Reprogramming and Identifies Key Regulators of Early Phases of Cicer-Ascochyta Interactions

Ascochyta blight (AB) caused by the filamentous fungus Ascochyta rabiei is a major threat to global chickpea production. The mechanisms underlying chickpea response to A. rabiei remain elusive to date. Here, we investigated the comparative transcriptional dynamics of AB-resistant and -susceptible chickpea genotypes upon A. rabiei infection, to understand the early host defense response. Our findings revealed that AB-resistant plants underwent rapid and extensive transcriptional reprogramming compared with a susceptible host. At the early stage (24 h postinoculation [hpi]), mainly cell-wall remodeling and secondary metabolite pathways were highly activated, while differentially expressed genes related to signaling components, such as protein kinases, transcription factors, and hormonal pathways, show a remarkable upsurge at 72 hpi, especially in the resistant genotype. Notably, our data suggest an imperative role of jasmonic acid, ethylene, and abscisic acid signaling in providing immunity against A. rabiei. Furthermore, gene co-expression networks and modules corroborated the importance of cell-wall remodeling, signal transduction, and phytohormone pathways. Hub genes such as MYB14, PRE6, and MADS-SOC1 discovered in these modules might be the master regulators governing chickpea immunity. Overall, we not only provide novel insights for comprehensive understanding of immune signaling components mediating AB resistance and susceptibility at early Cicer-Ascochyta interactions but, also, offer a valuable resource for developing AB-resistant chickpea. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Breeding and Genomics Interventions for Developing Ascochyta Blight Resistant Grain Legumes

Grain legumes are a key food source for ensuring global food security and sustaining agriculture. However, grain legume production is challenged by growing disease incidence due to global climate change. Ascochyta blight (AB) is a major disease, causing substantial yield losses in grain legumes worldwide. Harnessing the untapped reserve of global grain legume germplasm, landraces, and crop wild relatives (CWRs) could help minimize yield losses caused by AB infection in grain legumes. Several genetic determinants controlling AB resistance in various grain legumes have been identified following classical genetic and conventional breeding approaches. However, the advent of molecular markers, biparental quantitative trait loci (QTL) mapping, genome-wide association studies, genomic resources developed from various genome sequence assemblies, and whole-genome resequencing of global germplasm has revealed AB-resistant gene(s)/QTL/genomic regions/haplotypes on various linkage groups. These genomics resources allow plant breeders to embrace genomics-assisted selection for developing/transferring AB-resistant genomic regions to elite cultivars with great precision. Likewise, advances in functional genomics, especially transcriptomics and proteomics, have assisted in discovering possible candidate gene(s) and proteins and the underlying molecular mechanisms of AB resistance in various grain legumes. We discuss how emerging cutting-edge next-generation breeding tools, such as rapid generation advancement, field-based high-throughput phenotyping tools, genomic selection, and CRISPR/Cas9, could be used for fast-tracking AB-resistant grain legumes to meet the increasing demand for grain legume-based protein diets and thus ensuring global food security.

Gcorn fungi: A Web Tool for Detecting Biases between Gene Evolution and Speciation in Fungi

(1) Background: Fungi contain several millions of species, and the diversification of fungal genes has been achieved by speciation, gene duplication, and horizontal gene transfer. Although several databases provide information on orthologous and paralogous events, these databases show no information on biases between gene mutation and speciation. Here, we designed the Gcorn fungi database to better understand such biases. (2) Methods: Amino acid sequences of fungal genes in 249 species, which contain 2,345,743 sequences, were used for this database. Homologous genes were grouped at various thresholds of the homology index, which was based on the percentages of gene mutations. By grouping genes that showed highly similar homology indices to each other, we showed functional and evolutionary traits in the phylogenetic tree depicted for the gene of interest. (3) Results: Gcorn fungi provides well-summarized information on the evolution of a gene lineage and on the biases between gene evolution and speciation, which are quantitatively identified by the Robinson–Foulds metric. The database helps users visualize these traits using various depictions. (4) Conclusions: Gcorn fungi is an open access database that provides a variety of information with which to understand gene function and evolution.

Candidate Effectors of Plasmodiophora brassicae Pathotype 5X During Infection of Two Brassica napus Genotypes

Clubroot, caused by Plasmodiophora brassicae, is one of the most important diseases of canola (Brassica napus) in Canada. Disease management relies heavily on planting clubroot resistant (CR) cultivars, but in recent years, new resistance-breaking pathotypes of P. brassicae have emerged. Current efforts against the disease are concentrated in developing host resistance using traditional genetic breeding, omics and molecular biology. However, because of its obligate biotrophic nature, limited resources have been dedicated to investigating molecular mechanisms of pathogenic infection. We previously performed a transcriptomic study with the cultivar resistance-breaking pathotype 5X on two B. napus hosts presenting contrasting resistance/susceptibility, where we evaluated the mechanisms of host response. Since cultivar-pathotype interactions are very specific, and pathotype 5X is one of the most relevant resistance-breaking pathotypes in Canada, in this study, we analyze the expression of genes encoding putative secreted proteins from this pathotype, predicted using a bioinformatics pipeline, protein modeling and orthologous comparisons with effectors from other pathosystems. While host responses were found to differ markedly in our previous study, many common effectors are found in the pathogen while infecting both hosts, and the gene response among biological pathogen replicates seems more consistent in the effectors associated with the susceptible interaction, especially at 21 days after inoculation. The predicted effectors indicate the predominance of proteins with interacting domains (e.g., ankyrin), and genes bearing kinase and NUDIX domains, but also proteins with protective action against reactive oxygen species from the host. Many of these genes confirm previous predictions from other clubroot studies. A benzoic acid/SA methyltransferase (BSMT), which methylates SA to render it inactive, showed high levels of expression in the interactions with both hosts. Interestingly, our data indicate that E3 ubiquitin proteasome elements are also potentially involved in pathogenesis. Finally, a gene with similarity to indole-3-acetaldehyde dehydrogenase is a promising candidate effector because of its involvement in indole acetic acid synthesis, since auxin is one of the major players in clubroot development.

Increased abundance of secreted hydrolytic enzymes and secondary metabolite gene clusters define the genomes of latent plant pathogens in the Botryosphaeriaceae

BACKGROUND

The Botryosphaeriaceae are important plant pathogens, but also have the ability to establish asymptomatic infections that persist for extended periods in a latent state. In this study, we used comparative genome analyses to shed light on the genetic basis of the interactions of these fungi with their plant hosts. For this purpose, we characterised secreted hydrolytic enzymes, secondary metabolite biosynthetic gene clusters and general trends in genomic architecture using all available Botryosphaeriaceae genomes, and selected Dothideomycetes genomes.

RESULTS

The Botryosphaeriaceae genomes were rich in carbohydrate-active enzymes (CAZymes), proteases, lipases and secondary metabolic biosynthetic gene clusters (BGCs) compared to other Dothideomycete genomes. The genomes of Botryosphaeria, Macrophomina, Lasiodiplodia and Neofusicoccum, in particular, had gene expansions of the major constituents of the secretome, notably CAZymes involved in plant cell wall degradation. The Botryosphaeriaceae genomes were shown to have moderate to high GC contents and most had low levels of repetitive DNA. The genomes were not compartmentalized based on gene and repeat densities, but genes of secreted enzymes were slightly more abundant in gene-sparse regions.

CONCLUSION

The abundance of secreted hydrolytic enzymes and secondary metabolite BGCs in the genomes of Botryosphaeria, Macrophomina, Lasiodiplodia, and Neofusicoccum were similar to those in necrotrophic plant pathogens and some endophytes of woody plants. The results provide a foundation for comparative genomic analyses and hypotheses to explore the mechanisms underlying Botryosphaeriaceae host-plant interactions.

Genome sequences and in silico effector mining of Corynespora cassiicola CC_29 and Corynespora olivacea CBS 114450

The placement of Corynespora olivacea within the large genus Corynespora (Pleosporales) is controversial, because the species is distantly related to other congeners, including the type species C. cassiicola. Corynespora cassiicola is a polyphagous, cosmopolitan plant pathogen. Successful colonization of plant tissues requires the pathogen's effector repertoire to modulate host cell physiology and facilitate the infection process. We sequenced and performed functional annotations on the genomes of C. cassiicola CC_29 (genome size about 44.8 Mb; isolated from soybean leaves) and C. olivacea CBS 114450 (32.3 Mb). Our phylogenomic approach showed that C. cassiicola is distantly related to C. olivacea, which clustered among the Massarinaceae family members, supporting a hypothesis that C. olivacea was originally misclassified. The predicted sizes for the proteome and secretome of C. cassiicola (18,487 and 1327, respectively) were larger than those of C. olivacea (13,501 and 920; respectively). Corynespora cassiicola had a richer repertoire of effector proteins (CAZymes, proteases, lipases, and effectors) and genes associated with secondary metabolism than did C. olivacea.

De novo genome sequencing of mycoparasite Mycogone perniciosa strain MgR1 sheds new light on its biological complexity

Mycogone perniciosa is a mycoparasite causing Wet Bubble Diseases (WBD) of Agaricus bisporus. In the present study, the whole genome of M. perniciosa strain MgR1 was sequenced using Illumina NextSeq500 platform. This sequencing generated 8.03 Gb of high-quality data and a draft genome of 39 Mb was obtained through a de novo assembly of the high-quality reads. The draft genome resulted into prediction of 9276 genes from the 1597 scaffolds. NCBI-based homology analysis revealed the identification of 8660 genes. Notably, non-redundant protein database analysis of the M. perniciosa strain MgR1 revealed its close relation with the Trichoderma arundinaceum. Moreover, ITS-based phylogenetic analysis showed the highest similarity of M. perniciosa strain MgR1 with Hypomyces perniciosus strain CBS 322.22 and Mycogone perniciosa strain PPRI 5784. Annotation of the 3917 genes of M. perniciosa strain MgR1 grouped in three major categories viz. biological process (2583 genes), cellular component (2013 genes), and molecular function (2919 genes). UniGene analysis identified 2967 unique genes in M. perniciosa strain MgR1. In addition, prediction of the secretory and pathogenicity-related genes based on the fungal database indicates that 1512 genes (16% of predicted genes) encode for secretory proteins. Moreover, out of 9276 genes, 1296 genes were identified as pathogenesis-related proteins matching with 51 fungal and bacterial genera. Overall, the key pathogenic genes such as lysine M protein domain genes, G protein, hydrophobins, and cytochrome P450 were also observed. The draft genome of MgR1 provides an understanding of pathogenesis of WBD in A. bisporus and could be utilized to develop novel management strategies.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Modulation of fungal virulence through CRZ1 regulated F-BAR-dependent actin remodeling and endocytosis in chickpea infecting phytopathogen Ascochyta rabiei

Polarized hyphal growth of filamentous pathogenic fungi is an essential event for host penetration and colonization. The long-range early endosomal trafficking during hyphal growth is crucial for nutrient uptake, sensing of host-specific cues, and regulation of effector production. Bin1/Amphiphysin/Rvs167 (BAR) domain-containing proteins mediate fundamental cellular processes, including membrane remodeling and endocytosis. Here, we identified a F-BAR domain protein (ArF-BAR) in the necrotrophic fungus Ascochyta rabiei and demonstrate its involvement in endosome-dependent fungal virulence on the host plant Cicer arietinum. We show that ArF-BAR regulates endocytosis at the hyphal tip, localizes to the early endosomes, and is involved in actin dynamics. Functional studies involving gene knockout and complementation experiments reveal that ArF-BAR is necessary for virulence. The loss-of-function of ArF-BAR gene results in delayed formation of apical septum in fungal cells near growing hyphal tip that is crucial for host penetration, and impaired secretion of a candidate effector having secretory signal peptide for translocation across the endoplasmic reticulum membrane. The mRNA transcripts of ArF-BAR were induced in response to oxidative stress and infection. We also show that ArF-BAR is able to tubulate synthetic liposomes, suggesting the functional role of F-BAR domain in membrane tubule formation in vivo. Further, our studies identified a stress-induced transcription factor, ArCRZ1 (Calcineurin-responsive zinc finger 1), as key transcriptional regulator of ArF-BAR expression. We propose a model in which ArCRZ1 functions upstream of ArF-BAR to regulate A. rabiei virulence through a mechanism that involves endocytosis, effector secretion, and actin cytoskeleton regulation.

First Draft Genome Sequence Resource of a Strain of Pseudocercospora fijiensis Isolated in North America

Black Sigatoka disease, caused by the fungus Pseudocercospora fijiensis, is one of the most devastating diseases of banana around the world. Fungicide applications are the primary tool used to manage black Sigatoka, but fungicide resistance in P. fijiensis, as in other fungal pathogens, is one of the major limitations in the efficient management and prevention of this disease. In the current study, we present the draft genome of P. fijiensis strain IIL-20, the first genomic sequence published from a strain of this fungus isolated in North America. Bioinformatic analysis showed putative genes involved in fungus virulence and fungicide resistance. These findings may lead us to a better understanding of the molecular pathogenesis of this fungal pathogen and also to the discovery of the mechanisms conferring fungicide resistance.

Genomic analysis of facultatively oligotrophic haloarchaea of the genera Halarchaeum, Halorubrum, and Halolamina, isolated from solar salt

Extremely halophilic archaea (haloarchaea) belonging to the phylum Euryarchaeota have been found in high-salinity environments. In this study, Halarchaeum sp. CBA1220, Halorubrum sp. CBA1229, and Halolamina sp. CBA1230, which are facultatively oligotrophic haloarchaea, were isolated from solar salt by culture under oligotrophic culture conditions. The complete genomes of strains CBA1220, CBA1229, and CBA1230 were sequenced and were found to contain 3,175,875, 3,582,278, and 3,465,332 bp, with a G + C content of 68.25, 67.66, and 66.75 mol %, respectively. In total, 60, 36, and 33 carbohydrate-active enzyme genes were determined in the respective strains. The strains harbored various genes encoding stress-tolerance proteins, including universal stress proteins, cold-shock proteins, and rubrerythrin and rubrerythrin-related proteins. The genome data produced in this study will facilitate further research to improve our understanding of other halophilic strains and promote their industrial application.

Genome Sequence Data of the Soybean Pathogen Stagonosporopsis vannaccii: A Resource for Studies on Didymellaceae Evolution

The genus Stagonosporopsis is classified within the Didymellaceae family and has around 40 associated species. Among them, several species are important plant pathogens responsible for significant losses in economically important crops worldwide. Stagonosporopsis vannaccii is a newly described species pathogenic to soybean. Here, we present the draft whole-genome sequence, gene prediction, and annotation of S. vannaccii isolate LFN0148 (also known as IMI 507030). To our knowledge, this is the first genome sequenced of this species and represents a new useful source for future research on fungal comparative genomics studies.

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.

Transcript profiling reveals potential regulators for oxidative stress response of a necrotrophic chickpea pathogen Ascochyta rabiei

Necrotrophic pathogens experience host-generated oxidative stress during pathogenesis. They overcome such hostile environment by intricate mechanisms which are largely understudied. In this article, reference-based transcriptome analysis of a devastating Ascochyta Blight (AB) disease causing chickpea pathogen Ascochyta rabiei was explored to get insights into survival mechanisms under oxidative stress. Here, expression profiling of mock-treated and menadione-treated fungus was carried out by RNA-Seq approach. A significant number of genes in response to oxidative stress were overrepresented, suggestive of a robust and coordinated defense system of A. rabiei. A total 73 differentially expressed genes were filtered out from both the transcriptomes, among them 64 were up-regulated and 9 were found down-regulated. The gene ontology and KEGG mapping were conducted to comprehend the possible regulatory roles of differentially expressed genes in metabolic networks and biosynthetic pathways. Transcript profiling, KEGG pathway and gene ontology-based enrichment analysis revealed 12 (16.43%) stress responsive factors, 25 (34.24%) virulence associated genes, 10 (13.69%) putative effectors and 28 (38.35%) important interacting proteins associated with various metabolic pathways. In addition, genes with differential expression were further explored for underlying putative pathogenicity factors. We identified five genes ST47_g10291, ST47_g9396, ST47_g10294, ST47_g4395, and ST47_g7191 that were common to stress and fungal pathogenicity. The factors recognized in this work can be used to establish molecular tools to explain the regulatory gene networks engaged in stress response of fungal pathogens and disease management.

Whole-genome and time-course dual RNA-Seq analyses reveal chronic pathogenicity-related gene dynamics in the ginseng rusty root rot pathogen Ilyonectria robusta

Ilyonectria robusta causes rusty root rot, the most devastating chronic disease of ginseng. Here, we for the first time report the high-quality genome of the I. robusta strain CD-56. Time-course (36 h, 72 h, and 144 h) dual RNA-Seq analysis of the infection process was performed, and many genes, including candidate effectors, were found to be associated with the progression and success of infection. The gene expression profile of CD-56 showed a trend of initial inhibition and then gradually returned to a profile similar to that of the control. Analyses of the gene expression patterns and functions of pathogenicity-related genes, especially candidate effector genes, indicated that the stress response changed to an adaptive response during the infection process. For ginseng, gene expression patterns were highly related to physiological conditions. Specifically, the results showed that ginseng defenses were activated by CD-56 infection and persisted for at least 144 h thereafter but that the mechanisms invoked were not effective in preventing CD-56 growth. Moreover, CD-56 did not appear to fully suppress plant defenses, even in late stages after infection. Our results provide new insight into the chronic pathogenesis of CD-56 and the comprehensive and complex inducible defense responses of ginseng root to I. robusta infection.

Genome-scale data resolve ancestral rock-inhabiting lifestyle in Dothideomycetes (Ascomycota)

Dothideomycetes is the most diverse fungal class in Ascomycota and includes species with a wide range of lifestyles. Previous multilocus studies have investigated the taxonomic and evolutionary relationships of these taxa but often failed to resolve early diverging nodes and frequently generated inconsistent placements of some clades. Here, we use a phylogenomic approach to resolve relationships in Dothideomycetes, focusing on two genera of melanized, extremotolerant rock-inhabiting fungi, Lichenothelia and Saxomyces, that have been suggested to be early diverging lineages. We assembled phylogenomic datasets from newly sequenced (4) and previously available genomes (238) of 242 taxa. We explored the influence of tree inference methods, supermatrix vs. coalescent-based species tree, and the impact of varying amounts of genomic data. Overall, our phylogenetic reconstructions provide consistent and well-supported topologies for Dothideomycetes, recovering Lichenothelia and Saxomyces among the earliest diverging lineages in the class. In addition, many of the major lineages within Dothideomycetes are recovered as monophyletic, and the phylogenomic approach implemented strongly supports their relationships. Ancestral character state reconstruction suggest that the rock-inhabiting lifestyle is ancestral within the class.

De novo genome sequencing and secretome analysis of Tilletia indica inciting Karnal bunt of wheat provides pathogenesis-related genes

Tilletia indica is an internationally quarantined fungal pathogen causing Karnal bunt of wheat. The present study carried out that the whole genome of T. indica was sequenced and identified transposable elements, pathogenicity-related genes using a comparative genomics approach. The T. indica genome assembly size of 33.7 MB was generated using Illumina and Pac Bio platforms with GC content of 55.0%. A total of 1737 scaffolds were obtained with N₅₀ of 58,667 bp. The ab initio gene prediction was performed using Ustilago maydis as the reference species. A total number of 10,113 genes were predicted with an average gene size of 1945 bp out of which functionally annotated genes were 7262. A total number of 3216 protein-coding genes were assigned in different categories. Out of a total number of 1877 transposable elements, gypsy had the highest count (573). Total 5772 simple sequence repeats were identified in the genome assembly, and the most abundant simple sequence repeat type was trinucleotide having 42% of total SSRs. The comparative genome analysis suggested 3751 proteins of T. indica had orthologs in five fungi, whereas 126 proteins were unique to T. indica. Secretome analysis revealed the presence of 1014 secretory proteins and few carbohydrate-active enzymes in the genome. Some putative candidate pathogenicity-related genes were identified in the genome. The whole genome of T. indica will provide a window to understand the pathogenesis mechanism, fungal life cycle, survival of teliospores, and novel strategies for management of Karnal bunt disease of wheat.

Secretome Analysis and In Planta Expression of Salivary Proteins Identify Candidate Effectors from the Brown Planthopper Nilaparvata lugens

The brown planthopper (BPH), Nilaparvata lugens (Stål), is a phloem sap-feeding insect. During feeding on rice plants, BPH secretes salivary proteins with potential effector functions, which may play a critical role in the plant-insect interactions. However, a limited number of BPH effector proteins have been identified to date. Here, we sequenced the salivary gland transcriptomes of five BPH populations and subsequently established a N. lugens secretome consisting of 1,140 protein-encoding genes. Secretome analysis revealed the presence of both conserved and rapidly evolving salivary proteins. A screen for potential effectors that elicit responses in the plant was performed via the transient expression analysis of 64 BPH salivary proteins in Nicotiana benthamiana leaves and rice protoplasts. The salivary proteins Nl12, Nl16, Nl28, and Nl43 induced cell death, whereas Nl40 induced chlorosis and Nl32 induced a dwarf phenotype in N. benthamiana, indicating effector properties of these proteins. Ectopic expression of the six salivary proteins in N. benthamiana upregulated expression of defense-related genes and callose deposition. Tissue expression analysis showed a higher expression level of the six candidate effectors in salivary glands than in other tissues. Subcellular localization and analysis of the domain required for cell death showed a diverse structure of the six effectors. Nl28, Nl40, and Nl43 are N. lugens specific; in contrast, Nl12, Nl16, and Nl32 are conserved among insects. The Nl40 family has numerous isoforms produced by alternative splicing, exemplifying rapid evolution and expansion of effector proteins in the BPH. Our results suggest a potential large effector repertoire in BPH and a higher level of effector conservation exist in BPH compared with that in plant pathogens.

Carbon sources and XlnR-dependent transcriptional landscape of CAZymes in the industrial fungus Talaromyces versatilis: when exception seems to be the rule

Background

Research on filamentous fungi emphasized the remarkable redundancy in genes encoding hydrolytic enzymes, the similarities but also the large differences in their expression, especially through the role of the XlnR/XYR1 transcriptional activator. The purpose of this study was to evaluate the specificities of the industrial fungus Talaromyces versatilis, getting clues into the role of XlnR and the importance of glucose repression at the transcriptional level, to provide further levers for cocktail production.

Results

By studying a set of 62 redundant genes representative of several categories of enzymes, our results underlined the huge plasticity of transcriptional responses when changing nutritional status. As a general trend, the more heterogeneous the substrate, the more efficient to trigger activation. Genetic modifications of xlnR led to significant reorganisation of transcriptional patterns. Just a minimal set of genes actually fitted in a simplistic model of regulation by a transcriptional activator, and this under specific substrates. On the contrary, the diversity of xlnR⁺ versus ΔxlnR responses illustrated the existence of complex and unpredicted patterns of co-regulated genes that were highly dependent on the culture condition, even between genes that encode members of a functional category of enzymes. They notably revealed a dual, substrate-dependant repressor-activator role of XlnR, with counter-intuitive transcripts regulations that targeted specific genes. About glucose, it appeared as a formal repressive sugar as we observed a massive repression of most genes upon glucose addition to the mycelium grown on wheat straw. However, we also noticed a positive role of this sugar on the basal expression of a few genes, (notably those encoding cellulases), showing again the strong dependence of these regulatory mechanisms upon promoter and nutritional contexts.

Conclusions

The diversity of transcriptional patterns appeared to be the rule, while common and stable behaviour, both within gene families and with fungal literature, the exception. The setup of a new biotechnological process to reach optimized, if not customized expression patterns of enzymes, hence appeared tricky just relying on published data that can lead, in the best scenario, to approximate trends. We instead encourage preliminary experimental assays, carried out in the context of interest to reassess gene responses, as a mandatory step before thinking in (genetic) strategies for the improvement of enzyme production in fungi.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1062-8) contains supplementary material, which is available to authorized users.

Powdery Mildews Are Characterized by Contracted Carbohydrate Metabolism and Diverse Effectors to Adapt to Obligate Biotrophic Lifestyle

Powdery mildew is a widespread plant disease caused by obligate biotrophic fungal pathogens involving species-specific interactions between host and parasite. To gain genomic insights into the underlying obligate biotrophic mechanisms, we analyzed 15 microbial genomes covering powdery and downy mildews and rusts. We observed a genome-wide, massive contraction of multiple gene families in powdery mildews, such as enzymes in the carbohydrate metabolism pathway, when compared with ascomycete phytopathogens, while the fatty acid metabolism pathway maintained its integrity. We also observed significant differences in candidate secreted effector protein (CSEP) families between monocot and dicot powdery mildews, perhaps due to different selection forces. While CSEPs in monocot mildews are likely subject to positive selection causing rapid expansion, CSEP families in dicot mildews are shrinking under strong purifying selection. Our results not only illustrate obligate biotrophic mechanisms of powdery mildews driven by gene family evolution in nutrient metabolism, but also demonstrate how the divergence of CSEPs between monocot and dicot lineages might contribute to species-specific adaption.

The Landscape of Repetitive Elements in the Refined Genome of Chilli Anthracnose Fungus Colletotrichum truncatum

The ascomycete fungus Colletotrichum truncatum is a major phytopathogen with a broad host range which causes anthracnose disease of chilli. The genome sequencing of this fungus led to the discovery of functional categories of genes that may play important roles in fungal pathogenicity. However, the presence of gaps in C. truncatum draft assembly prevented the accurate prediction of repetitive elements, which are the key players to determine the genome architecture and drive evolution and host adaptation. We re-sequenced its genome using single-molecule real-time (SMRT) sequencing technology to obtain a refined assembly with lesser and smaller gaps and ambiguities. This enabled us to study its genome architecture by characterising the repetitive sequences like transposable elements (TEs) and simple sequence repeats (SSRs), which constituted 4.9 and 0.38% of the assembled genome, respectively. The comparative analysis among different Colletotrichum species revealed the extensive repeat rich regions, dominated by Gypsy superfamily of long terminal repeats (LTRs), and the differential composition of SSRs in their genomes. Our study revealed a recent burst of LTR amplification in C. truncatum, C. higginsianum, and C. scovillei. TEs in C. truncatum were significantly associated with secretome, effectors and genes in secondary metabolism clusters. Some of the TE families in C. truncatum showed cytosine to thymine transitions indicative of repeat-induced point mutation (RIP). C. orbiculare and C. graminicola showed strong signatures of RIP across their genomes and "two-speed" genomes with extensive AT-rich and gene-sparse regions. Comparative genomic analyses of Colletotrichum species provided an insight into the species-specific SSR profiles. The SSRs in the coding and non-coding regions of the genome revealed the composition of trinucleotide repeat motifs in exons with potential to alter the translated protein structure through amino acid repeats. This is the first genome-wide study of TEs and SSRs in C. truncatum and their comparative analysis with six other Colletotrichum species, which would serve as a useful resource for future research to get insights into the potential role of TEs in genome expansion and evolution of Colletotrichum fungi and for development of SSR-based molecular markers for population genomic studies.

mQTL-seq and classical mapping implicates the role of an AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED (AHL) family gene in Ascochyta blight resistance of chickpea

Ascochyta blight (AB) caused by the fungal pathogen Ascochyta rabiei is a serious foliar disease of chickpea (Cicer arietinum L.). Despite many genetic studies on chickpea-Ascochyta interaction, genome-wide scan of chickpea for the identification of AB-associated quantitative trait loci (QTLs) and their gene(s) has not been accomplished. To elucidate narrow QTLs for AB resistance, here, we report the use of multiple QTL-sequencing approach on 2 sets of extreme AB phenotype bulks derived from Cicer intraspecific and interspecific crosses. Two major QTLs, qABR4.1 and qABR4.2, and a minor QTL, qABR4.3, were identified on assembled chickpea pseudomolecule 4. We narrowed qABR4.1 to a "robust region" at 4.568-4.618 Mb through mapping on a larger intraspecific cross-derived population and comparative analysis. Among 4 genes, the CaAHL18 gene showed higher expression under Ascochyta stress in AB resistant parent suggesting that it is the candidate gene under "robust qABR4.1." Dual-luciferase assay with CaAHL18 polymorphic cis-regulatory sequences showed that allelic variation is associated with higher expression. Thus, our findings on chickpea-Ascochyta interaction have narrowed down AB resistance associated QTLs on chickpea physical map. The narrowed QTLs and gene-associated markers will help in biotechnological and breeding programs for chickpea improvement.

Annotation resource of tandem repeat-containing secretory proteins in sixty fungi

Fungal secretory proteins that interact with host plants are regarded as effectors. Because fungal effectors rarely contain conserved sequence features, identification and annotation of fungal effectors from predicted secretory proteins are difficult using outward comparison methods such as BLAST or hidden Markov model. In desire of more sequence features to prioritize research interests of fungal secretory proteins, this study developed a pipeline to identify tandem repeat (TR) domain within putative secretory proteins and tested a hypothesis that at least one type of TR domain in non-orthologous secretory proteins has emerged from convergent evolution for plant pathogenicity. There were 2804 types of TR domains and a total of 2925 TR-containing secretory proteins found from 60 fungi. There was no conserved type of TR domain shared only by plant pathogens, indicating functional divergence for different types of TR domain and TR-containing secretory proteins. The annotation resource of putative fungal TR-containing secretory proteins provides new sequence features that will be useful for the community interested in fungal effector biology.

Lifestyle, gene gain and loss, and transcriptional remodeling cause divergence in the transcriptomes of Phytophthora infestans and Pythium ultimum during potato tuber colonization

BACKGROUND

How pathogen genomes evolve to support distinct lifestyles is not well-understood. The oomycete Phytophthora infestans, the potato blight agent, is a largely biotrophic pathogen that feeds from living host cells, which become necrotic only late in infection. The related oomycete Pythium ultimum grows saprophytically in soil and as a necrotroph in plants, causing massive tissue destruction. To learn what distinguishes their lifestyles, we compared their gene contents and expression patterns in media and a shared host, potato tuber.

RESULTS

Genes related to pathogenesis varied in temporal expression pattern, mRNA level, and family size between the species. A family's aggregate expression during infection was not proportional to size due to transcriptional remodeling and pseudogenization. Ph. infestans had more stage-specific genes, while Py. ultimum tended towards more constitutive expression. Ph. infestans expressed more genes encoding secreted cell wall-degrading enzymes, but other categories such as secreted proteases and ABC transporters had higher transcript levels in Py. ultimum. Species-specific genes were identified including new Pythium genes, perforins, which may disrupt plant membranes. Genome-wide ortholog analyses identified substantial diversified expression, which correlated with sequence divergence. Pseudogenization was associated with gene family expansion, especially in gene clusters.

CONCLUSION

This first large-scale analysis of transcriptional divergence within oomycetes revealed major shifts in genome composition and expression, including subfunctionalization within gene families. Biotrophy and necrotrophy seem determined by species-specific genes and the varied expression of shared pathogenicity factors, which may be useful targets for crop protection.

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

Transcription factors (TFs) are the key players in gene expression and their study is highly significant for shedding light on the molecular mechanisms and evolutionary history of organisms. During host-pathogen interaction, extensive reprogramming of gene expression facilitated by TFs is likely to occur in both host and pathogen. To date, the knowledge about TF repertoire in filamentous fungi is in infancy. The necrotrophic fungus Ascochyta rabiei, that causes destructive Ascochyta blight (AB) disease of chickpea (Cicer arietinum), demands more comprehensive study for better understanding of Ascochyta-legume pathosystem. In the present study, we performed the genome-wide identification and analysis of TFs in A. rabiei. Taking advantage of A. rabiei genome sequence, we used a bioinformatic approach to predict the TF repertoire of A. rabiei. For identification and classification of A. rabiei TFs, we designed a comprehensive pipeline using a combination of BLAST and InterProScan software. A total of 381 A. rabiei TFs were predicted and divided into 32 fungal specific families of TFs. The gene structure, domain organization and phylogenetic analysis of abundant families of A. rabiei TFs were also carried out. Comparative study of A. rabiei TFs with that of other necrotrophic, biotrophic, hemibiotrophic, symbiotic, and saprotrophic fungi was performed. It suggested presence of both conserved as well as unique features among them. Moreover, cis-acting elements on promoter sequences of earlier predicted A. rabiei secretome were also identified. With the help of published A. rabiei transcriptome data, the differential expression of TF and secretory protein coding genes was analyzed. Furthermore, comprehensive expression analysis of few selected A. rabiei TFs using quantitative real-time polymerase chain reaction revealed variety of expression patterns during host colonization. These genes were expressed in at least one of the time points tested post infection. Overall, this study illustrates the first genome-wide identification and analysis of TF repertoire of A. rabiei. This work would provide the basis for further studies to dissect role of TFs in the molecular mechanisms during A. rabiei-chickpea interactions.

Genomics and Comparative Genomic Analyses Provide Insight into the Taxonomy and Pathogenic Potential of Novel Emmonsia Pathogens

Over the last 50 years, newly described species of Emmonsia-like fungi have been implicated globally as sources of systemic human mycosis (emmonsiosis). Their ability to convert into yeast-like cells capable of replication and extra-pulmonary dissemination during the course of infection differentiates them from classical Emmonsia species. Immunocompromised patients are at highest risk of emmonsiosis and exhibit high mortality rates. In order to investigate the molecular basis for pathogenicity of the newly described Emmonsia species, genomic sequencing and comparative genomic analyses of Emmonsia sp. 5z489, which was isolated from a non-deliberately immunosuppressed diabetic patient in China and represents a novel seventh isolate of Emmonsia-like fungi, was performed. The genome size of 5z489 was 35.5 Mbp in length, which is ~5 Mbp larger than other Emmonsia strains. Further, 9,188 protein genes were predicted in the 5z489 genome and 16% of the assembly was identified as repetitive elements, which is the largest abundance in Emmonsia species. Phylogenetic analyses based on whole genome data classified 5z489 and CAC-2015a, another novel isolate, as members of the genus Emmonsia. Our analyses showed that divergences among Emmonsia occurred much earlier than other genera within the family Ajellomycetaceae, suggesting relatively distant evolutionary relationships among the genus. Through comparisons of Emmonsia species, we discovered significant pathogenicity characteristics within the genus as well as putative virulence factors that may play a role in the infection and pathogenicity of the novel Emmonsia strains. Moreover, our analyses revealed a novel distribution mode of DNA methylation patterns across the genome of 5z489, with >50% of methylated bases located in intergenic regions. These methylation patterns differ considerably from other reported fungi, where most methylation occurs in repetitive loci. It is unclear if this difference is related to physiological adaptations of new Emmonsia, but this question warrants further investigation. Overall, our analyses provide a framework from which to further study the evolutionary dynamics of Emmonsia strains and identity the underlying molecular mechanisms that determine the infectious and pathogenic potency of these fungal pathogens, and also provide insight into potential targets for therapeutic intervention of emmonsiosis and further research.

Genome Analysis of a Zygomycete Fungus Choanephora cucurbitarum Elucidates Necrotrophic Features Including Bacterial Genes Related to Plant Colonization

A zygomycete fungus, Choanephora cucurbitarum is a plant pathogen that causes blossom rot in cucurbits and other plants. Here we report the genome sequence of Choanephora cucurbitarum KUS-F28377 isolated from squash. The assembled genome has a size of 29.1 Mbp and 11,977 protein-coding genes. The genome analysis indicated that C. cucurbitarum may employ a plant pathogenic mechanism similar to that of bacterial plant pathogens. The genome contained 11 genes with a Streptomyces subtilisin inhibitor-like domain, which plays an important role in the defense against plant immunity. This domain has been found only in bacterial genomes. Carbohydrate active enzyme analysis detected 312 CAZymes in this genome where carbohydrate esterase family 6, rarely found in dikaryotic fungal genomes, was comparatively enriched. The comparative genome analysis showed that the genes related to sexual communication such as the biosynthesis of β-carotene and trisporic acid were conserved and diverged during the evolution of zygomycete genomes. Overall, these findings will help us to understand how zygomycetes are associated with plants.