Gene family expansions and contractions are associated with host range in plant pathogens of the genus Colletotrichum

Guanosine-specific single-stranded ribonuclease effectors of a phytopathogenic fungus potentiate host immune responses

Plants activate immunity upon recognition of pathogen-associated molecular patterns. Although phytopathogens have evolved a set of effector proteins to counteract plant immunity, some effectors are perceived by hosts and induce immune responses. Here, we show that two secreted ribonuclease effectors, SRN1 and SRN2, encoded in a phytopathogenic fungus, Colletotrichum orbiculare, induce cell death in a signal peptide- and catalytic residue-dependent manner, when transiently expressed in Nicotiana benthamiana. The pervasive presence of SRN genes across Colletotrichum species suggested the conserved roles. Using a transient gene expression system in cucumber (Cucumis sativus), an original host of C. orbiculare, we show that SRN1 and SRN2 potentiate host pattern-triggered immunity responses. Consistent with this, C. orbiculare SRN1 and SRN2 deletion mutants exhibited increased virulence on the host. In vitro analysis revealed that SRN1 specifically cleaves single-stranded RNAs at guanosine, leaving a 3'-end phosphate. Importantly, the potentiation of C. sativus responses by SRN1 and SRN2, present in the apoplast, depends on ribonuclease catalytic residues. We propose that the pathogen-derived apoplastic guanosine-specific single-stranded endoribonucleases lead to immunity potentiation in plants.

Comparative genomics points to tandem duplications of SAD gene clusters as drivers of increased α-linolenic (ω-3) content in S. hispanica seeds

Salvia hispanica L. (chia) is a source of abundant ω-3 polyunsaturated fatty acids (ω-3-PUFAs) that are highly beneficial to human health. The genomic basis for this accrued ω-3-PUFA content in this emerging crop was investigated through the assembly and comparative analysis of a chromosome-level reference genome for S. hispanica. The highly contiguous 321.5-Mbp genome assembly covering all six chromosomes enabled the identification of 32,922 protein-coding genes. Two whole-genome duplications (WGD) events were identified in the S. hispanica lineage. However, these WGD events could not be linked to the high α-linolenic acid (ALA, ω-3) accumulation in S. hispanica seeds based on phylogenomics. Instead, our analysis supports the hypothesis that evolutionary expansion through tandem duplications of specific lipid gene families, particularly the stearoyl-acyl carrier protein desaturase (ShSAD) gene family, is the main driver of the abundance of ω-3-PUFAs in S. hispanica seeds. The insights gained from the genomic analysis of S. hispanica will help establish a molecular breeding target that can be leveraged through genome editing techniques to increase ω-3 content in oil crops.

Evolutionary history of the cytochrome P450s from Colletotrichum species and prediction of their putative functional roles during host-pathogen interactions

The genomes of species belonging to the genus Colletotrichum harbor a substantial number of cytochrome P450 monooxygenases (CYPs) encoded by a broad diversity of gene families. However, the biological role of their CYP complement (CYPome) has not been elucidated. Here, we investigated the putative evolutionary scenarios that occurred during the evolution of the CYPome belonging to the Colletotrichum Graminicola species complex (s.c.) and their biological implications. The study revealed that most of the CYPome gene families belonging to the Graminicola s.c. experienced gene contractions. The reductive evolution resulted in species restricted CYPs are predominant in each CYPome of members from the Graminicola s.c., whereas only 18 families are absolutely conserved among these species. However, members of CYP families displayed a notably different phylogenetic relationship at the tertiary structure level, suggesting a putative convergent evolution scenario. Most of the CYP enzymes of the Graminicola s.c. share redundant functions in secondary metabolite biosynthesis and xenobiotic metabolism. Hence, this current work suggests that the presence of a broad CYPome in the genus Colletotrichum plays a critical role in the optimization of the colonization capability and virulence.

Genome evolution and transcriptome plasticity is associated with adaptation to monocot and dicot plants in Colletotrichum fungi

BACKGROUND

Colletotrichum fungi infect a wide diversity of monocot and dicot hosts, causing diseases on almost all economically important plants worldwide. Colletotrichum is also a suitable model for studying gene family evolution on a fine scale to uncover events in the genome associated with biological changes.

RESULTS

Here we present the genome sequences of 30 Colletotrichum species covering the diversity within the genus. Evolutionary analyses revealed that the Colletotrichum ancestor diverged in the late Cretaceous in parallel with the diversification of flowering plants. We provide evidence of independent host jumps from dicots to monocots during the evolution of Colletotrichum, coinciding with a progressive shrinking of the plant cell wall degradative arsenal and expansions in lineage-specific gene families. Comparative transcriptomics of 4 species adapted to different hosts revealed similarity in gene content but high diversity in the modulation of their transcription profiles on different plant substrates. Combining genomics and transcriptomics, we identified a set of core genes such as specific transcription factors, putatively involved in plant cell wall degradation.

CONCLUSIONS

These results indicate that the ancestral Colletotrichum were associated with dicot plants and certain branches progressively adapted to different monocot hosts, reshaping the gene content and its regulation.

Spirochaete genome identified in red abalone sample represents a novel genus Candidatus Haliotispira gen. nov. within the order Spirochaetales

A fully assembled spirochaete genome was identified as a contaminating scaffold in our red abalone (Haliotis rufescens) genome assembly. In this paper, we describe the analysis of this bacterial genome. The assembled spirochaete genome is 3.25 Mb in size with 48.5 mol% G+C content. The proteomes of 38 species were compared with the spirochaete genome and it was discovered to form an independent branch within the family Spirochaetaceae on the phylogenetic tree. The comparison of 16S rRNA sequences and average nucleotide identity scores between the spirochaete genome with known species of different families in Spirochaetia indicate that it is an unknown species. Further, the percentage of conserved proteins compared to neighbouring taxa confirm that it does not belong to a known genus within Spirochaetaceae. We propose the name Candidatus Haliotispira prima gen. nov., sp. nov. based on its taxonomic placement and origin. We also tested for the presence of this species in different species of abalone and found that it is also present in white abalone (Haliotis sorenseni). In addition, we highlight the need for better classification of taxa within the class Spirochaetia.

Gene Family Expansion during the Adaptation of Colletotrichum gloeosporioides to Woody Plants

Gene gains/losses during evolution are critical for the adaptation of organisms to new environments or hosts. However, it remains unknown whether gene family expansions facilitated the adaptation of phytopathogenic fungi to woody plants. In this study, we compared the newly sequenced genome of the Colletotrichum gloeosporioides strain CFCC80308 with the genomes of two other C. gloeosporioides strains, Cg-14 and Lc-1, isolated from Persea americana and Liriodendron leaves, respectively. The genes in the expanded families, which were associated with plant surface signal recognition, encoded various proteins, including glycosyde hydrolases (GHs) and cytochrome P450. Interestingly, there was a substantial increase in the number of GH family genes in CFCC80308. Specifically, there were 368 enriched genes in the GH families (e.g., GH1, GH3, GH10, GH12, GH15, GH16, GH17, GH18, GH25, GH32, GH53, GH61, GH76, and GH81); the expression levels of these genes were highly up-regulated during the infection of poplar trees. Additionally, the GH17 family was larger in CFCC80308 than in C. gloeosporioides strains Cg-14 and Lc-1. Furthermore, the expansion of the MP65-encoding gene family during the adaptation of Colletotrichum species to woody plants was consistent with the importance of gene gains/losses for the adaptation of organisms to their environments. This study has clarified how C. gloeosporioides adapted to woody plants during evolution.

Rethinking underutilized cereal crops: pan-omics integration and green system biology

MAIN CONCLUSION

Due to harsh lifestyle changes, in the present era, nutritional security is needed along with food security so it is necessary to include underutilized cereal crops (UCCs) in our daily diet to counteract the rising danger of human metabolic illness. We can attain both the goal of zero hunger and nutritional security by developing improved UCCs using advanced pan-omics (genomics, transcriptomics, proteomics, metabolomics, nutrigenomics, phenomics and ionomics) practices. Plant sciences research progressed profoundly since the last few decades with the introduction of advanced technologies and approaches, addressing issues of food demand of the growing population, nutritional security challenges and climate change. However, throughout the expansion and popularization of commonly consumed major cereal crops such as wheat and rice, other cereal crops such as millet, rye, sorghum, and others were impeded, despite their potential medicinal and nutraceutical qualities. Undoubtedly neglected underutilized cereal crops (UCCs) also have the capability to withstand diverse climate change. To relieve the burden of major crops, it is necessary to introduce the new crops in our diet in the way of UCCs. Introgression of agronomically and nutritionally important traits by pan-omics approaches in UCCs could be a defining moment for the population's well-being on the globe. This review discusses the importance of underutilized cereal crops, as well as the application of contemporary omics techniques and advanced bioinformatics tools that could open up new avenues for future study and be valuable assets in the development and usage of UCCs in the perspective of green system biology. The increased and improved use of UCCs is dependent on number of factors that necessitate a concerted research effort in agricultural sciences. The emergence of functional genomics with molecular genetics might gear toward the reawakening of interest in underutilized cereals crops. The need of this era is to focus on potential UCCs in advanced agriculture and breeding programmes. Hence, targeting the UCCs, might provide a bright future for better health and scientific rationale for its use.

Phylogeographic and Phylogenomic Structure of the Quarantine Plant Pathogen Colletotrichum liriopes, Including New Reports in the United States

Global agricultural trade has accelerated the emergence and re-emergence of new plant pathogens. In the United States, the fungal pathogen Colletotrichum liriopes is still considered a foreign quarantine pathogen that affects ornamental plants (i.e., Liriope spp.). Even though this species has been reported in East Asia on various asparagaceous hosts, its first and only report in the United States was in 2018. However, that study used only ITS nrDNA for identification, and no available culture or voucher specimen was maintained. The main objective of the present study was to determine the geographic and host distribution of specimens identified as C. liriopes. To accomplish this, new and existing isolates, sequences, and genomes obtained from various hosts and geographic locations (i.e., China, Colombia, Mexico, and the United States) were compared with the ex-type of C. liriopes. Multilocus phylogenetic (ITS, Tub2, GAPDH, CHS-1, and HIS3), phylogenomic, and splits tree analyses revealed that all the studied isolates/sequences form a well-supported clade with little intraspecific variation. Morphological characterizations support these findings. The minimum spanning network, low nucleotide diversity, and negative Tajima's D from both multilocus and genomic data suggest that there was a recent movement/invasion of a few East Asian genotypes to other countries where the ornamental plants are produced (e.g., South America) and subsequently to the importing countries, such as the United States. The study reveals that the geographic and host distribution of C. liriopes sensu stricto is expanded to the United States (i.e., at least Maryland, Mississippi, and Tennessee) and on various hosts in addition to Asparagaceae and Orchidaceae. The present study produces fundamental knowledge that can be used in efforts to reduce costs or losses from agricultural trade and to expand our understanding of pathogen movement.

Whole genome resequencing reveal patterns of genetic variation within Colletotrichum acutatum species complex from rubber trees in China

The Colletotrichum acutatum species complex possesses a diverse number of important traits, such as a wide host range and host preference, different modes of reproduction, and different strategies of host infection. Research using comparative genomics has attempted to find correlations between these traits. Here, we used multi-locus techniques and gene genealogical concordance analysis to investigate the phylogenetic relationships and taxonomic status of the Colletotrichum acutatum species complex using field isolates obtained from rubber trees. The results revealed that the dominant species was C. australisinense, followed by C. bannaense, while strain YNJH17109 was identified as C. laticiphilum. The taxonomic status of strains YNLC510 and YNLC511 was undetermined. Using whole-genome single nucleotide polymorphism data to analyze population structure, 18 strains of C. australisinense were subsequently divided into four populations, one of which was derived from an admixture of two populations. In addition, the strains LD1687, GD1628, and YNLC516, did not belong to any populations, and were considered to be admixtures of two or more populations. A split decomposition network analysis also provided evidence for genetic recombination within Colletotrichum acutatum species complex from rubber trees in China. Overall, a weak phylogeographic sub-structure was observed. Analysis also revealed significant differences in morphological characters and levels of virulence between populations.

Structural insights into latency of the metallopeptidase ulilysin (lysargiNase) and its unexpected inhibition by a sulfonyl-fluoride inhibitor of serine peptidases

Peptidases are regulated by latency and inhibitors, as well as compatibilization and cofactors. Ulilysin from Methanosarcina acetivorans, also called lysargiNase, is an archaeal metallopeptidase (MP) that is biosynthesized as a zymogen with a 60-residue N-terminal prosegment (PS). In the presence of calcium, it self-activates to yield the mature enzyme, which specifically cleaves before basic residues and thus complements trypsin in proteomics workflows. Here, we obtained a low-resolution crystal structure of proulilysin, in which 28 protomers arranged as 14 dimers form a continuous double helix of 544 Å pitch that parallels cell axis b of the crystal. The PS includes two α-helices and obstructs the active-site cleft of the catalytic domain (CD) by traversing it in the opposite orientation of a substrate, and a cysteine blocks the catalytic zinc according to a "cysteine-switch mechanism". Moreover, the PS interacts through its first helix with an "S-loop" of the CD, which acts as an "activation segment" that lacks one of two essential calcium cations. Upon PS removal during maturation, the S-loop adopts its competent conformation and binds the second calcium ion. Next, we found that in addition to general MP inhibitors, ulilysin was competitively and reversibly inhibited by 4-(2-aminoethyl)benzenesulfonyl fluoride (AEBSF; K_i = 4 μM). This is a compound that normally forms an irreversible covalent complex with serine peptidases but does not inhibit MPs. A high-resolution crystal structure of the complex revealed that the inhibitor penetrates the specificity pocket of ulilysin. A primary amine of the inhibitor salt-bridges an aspartate at the pocket bottom, thus mimicking the basic side chain of substrates. In contrast, the sulfonyl fluoride warhead is not involved and the catalytic zinc ion is freely accessible. Thus, the usage of inhibitor cocktails of peptidases, which typically contain AEBSF at ∼25-fold higher concentrations than the determined K_i, should be avoided when working with ulilysin. Finally, the structure of the complex, which occurred as a crystallographic dimer recurring in previous mature ulilysin structures, unveiled an N-terminal product fragment that delineated the non-primed side of the cleft. These results complement prior structures of ulilysin with primed-side product fragments and inhibitors.

Population Genomics Provide Insights into the Global Genetic Structure of Colletotrichum graminicola, the Causal Agent of Maize Anthracnose

Understanding the genetic diversity and mechanisms underlying genetic variation in pathogen populations is crucial to the development of effective control strategies. We investigated the genetic diversity and reproductive biology of Colletotrichum graminicola isolates which infect maize by sequencing the genomes of 108 isolates collected from 14 countries using restriction site-associated DNA sequencing (RAD-seq) and whole-genome sequencing (WGS). Clustering analyses based on single-nucleotide polymorphisms revealed three genetic groups delimited by continental origin, compatible with short-dispersal of the pathogen and geographic subdivision. Intra- and intercontinental migration was observed between Europe and South America, likely associated with the movement of contaminated germplasm. Low clonality, evidence of genetic recombination, and high phenotypic diversity were detected. We show evidence that, although it is rare (possibly due to losses of sexual reproduction- and meiosis-associated genes) C. graminicola can undergo sexual recombination. Our results support the hypotheses that intra- and intercontinental pathogen migration and genetic recombination have great impacts on the C. graminicola population structure. IMPORTANCE Plant pathogens cause significant reductions in yield and crop quality and cause enormous economic losses worldwide. Reducing these losses provides an obvious strategy to increase food production without further degrading natural ecosystems; however, this requires knowledge of the biology and evolution of the pathogens in agroecosystems. We employed a population genomics approach to investigate the genetic diversity and reproductive biology of the maize anthracnose pathogen (Colletotrichum graminicola) in 14 countries. We found that the populations are correlated with their geographical origin and that migration between countries is ongoing, possibly caused by the movement of infected plant material. This result has direct implications for disease management because migration can cause the movement of more virulent and/or fungicide-resistant genotypes. We conclude that genetic recombination is frequent (in contrast to the traditional view of C. graminicola being mainly asexual), which strongly impacts control measures and breeding programs aimed at controlling this disease.

Characterization of the microbiota dynamics associated with Moniliophthora roreri, causal agent of cocoa frosty pod rot disease, reveals new viral species

Introduction

Theobroma cacao, the cocoa tree, is a target for pathogens, such as fungi from the genera Phytophthora, Moniliophthora, Colletotrichum, Ceratocystis, among others. Some cacao pathogens are restricted to specific regions of the world, such as the Cacao swollen shoot virus (CSSV) in West African countries, while others are expanding geographically, such as Moniliophthora roreri in the Americas. M. roreri is one of the most threatening cacao pathogens since it directly attacks the cacao pods driving a significant reduction in production, and therefore economic losses. Despite its importance, the knowledge about the microenvironment of this pathogen and the cocoa pods is still poorly characterized.

Methods

Herein we performed RNA sequencing of spores in differential stages of culture in a medium supplemented with cacao pod extract and mycelium collected of the susceptible variety ICT 7121 naturally infected by the pathogen to evaluate the diversity and transcriptional activity of microorganisms associated with the in vitro sporulation of M. roreri.

Results

Our data revealed a great variety of fungi and bacteria associated with M. roreri, with an exceptional diversity of individuals from the genus Trichoderma sp. Interestingly, the dynamics of microorganisms from different kingdoms varied proportionally, suggesting they are somehow affected by M. roreri culture time. We also identified three sequences similar to viral genomes from the Narnaviridae family, posteriorly confirmed by phylogenetic analysis as members of the genus Narnavirus. Screening of M. roreri public datasets indicated the virus sequences circulating in samples from Ecuador, suggesting a wide spread of these elements. Of note, we did not identify traces of the viral sequences in the M. roreri genome or DNA sequencing, restricting the possibility of these sequences representing endogenized elements.

Discussion

To the best of our knowledge, this is the first report of viruses infecting the fungus of the genus Moniliophthora and only the third description of viruses that are able to parasite elements from the Marasmiaceae family.

Metabolic modelling of the human gut microbiome in type 2 diabetes patients in response to metformin treatment

The human gut microbiome has been associated with several metabolic disorders including type 2 diabetes mellitus. Understanding metabolic changes in the gut microbiome is important to elucidate the role of gut bacteria in regulating host metabolism. Here, we used available metagenomics data from a metformin study, together with genome-scale metabolic modelling of the key bacteria in individual and community-level to investigate the mechanistic role of the gut microbiome in response to metformin. Individual modelling predicted that species that are increased after metformin treatment have higher growth rates in comparison to species that are decreased after metformin treatment. Gut microbial enrichment analysis showed prior to metformin treatment pathways related to the hypoglycemic effect were enriched. Our observations highlight how the key bacterial species after metformin treatment have commensal and competing behavior, and how their cellular metabolism changes due to different nutritional environment. Integrating different diets showed there were specific microbial alterations between different diets. These results show the importance of the nutritional environment and how dietary guidelines may improve drug efficiency through the gut microbiota.

Colletotrichum truncatum Causing Anthracnose of Tomato (Solanum lycopersicum L.) in Malaysia

Tomato (Solanum lycopersicum L.) is a popular nutritious vegetable crop grown in Malaysia and other parts of the world. However, fungal diseases such as anthracnose pose significant threats to tomato production by reducing the fruit quality and food value of tomato, resulting in lower market prices of the crop globally. In the present study, the etiology of tomato anthracnose was investigated in commercial tomato farms in Sabah, Malaysia. A total of 22 fungal isolates were obtained from anthracnosed tomato fruits and identified as Colletotrichum species, using morphological characteristics. The phylogenetic relationships of multiple gene sequence alignments such as internal transcribed spacer (ITS), β-tubulin (tub2), glyceraldehyde 3-phosphate dehydrogenase (gapdh), actin (act), and calmodulin (cal), were adopted to accurately identify the Colletotrichum species as C. truncatum. The results of pathogenicity tests revealed that all C. truncatum isolates caused anthracnose disease symptoms on inoculated tomato fruits. To our knowledge, the present study is the first report of tomato anthracnose caused by C. truncatum in Malaysia. The findings of this study will be helpful in disease monitoring, and the development of strategies for effective control of anthracnose on tomato fruits.

Chromosome-level analysis of the Colletotrichum graminicola genome reveals the unique characteristics of core and minichromosomes

The fungal pathogen Colletotrichum graminicola causes the anthracnose of maize (Zea mays) and is responsible for significant yield losses worldwide. The genome of C. graminicola was sequenced in 2012 using Sanger sequencing, 454 pyrosequencing, and an optical map to obtain an assembly of 13 pseudochromosomes. We re-sequenced the genome using a combination of short-read (Illumina) and long-read (PacBio) technologies to obtain a chromosome-level assembly. The new version of the genome sequence has 13 chromosomes with a total length of 57.43 Mb. We detected 66 (23.62 Mb) structural rearrangements in the new assembly with respect to the previous version, consisting of 61 (21.98 Mb) translocations, 1 (1.41 Mb) inversion, and 4 (221 Kb) duplications. We annotated the genome and obtained 15,118 predicted genes and 3,614 new gene models compared to the previous version of the assembly. We show that 25.88% of the new assembly is composed of repetitive DNA elements (13.68% more than the previous assembly version), which are mostly found in gene-sparse regions. We describe genomic compartmentalization consisting of repeat-rich and gene-poor regions vs. repeat-poor and gene-rich regions. A total of 1,140 secreted proteins were found mainly in repeat-rich regions. We also found that ~75% of the three smallest chromosomes (minichromosomes, between 730 and 551 Kb) are strongly affected by repeat-induced point mutation (RIP) compared with 28% of the larger chromosomes. The gene content of the minichromosomes (MCs) comprises 121 genes, of which 83.6% are hypothetical proteins with no predicted function, while the mean percentage of Chr1-Chr10 is 36.5%. No predicted secreted proteins are present in the MCs. Interestingly, only 2% of the genes in Chr11 have homologs in other strains of C. graminicola, while Chr12 and 13 have 58 and 57%, respectively, raising the question as to whether Chrs12 and 13 are dispensable. The core chromosomes (Chr1-Chr10) are very different with respect to the MCs (Chr11-Chr13) in terms of the content and sequence features. We hypothesize that the higher density of repetitive elements and RIPs in the MCs may be linked to the adaptation and/or host co-evolution of this pathogenic fungus.

Functional analysis of the Nep1-like proteins from Plasmopara viticola

Necrosis and ethylene-inducing peptide 1 (Nep1) -like proteins (NLP) are secreted by multiple taxonomically unrelated plant pathogens (bacteria, fungi, and oomycete) and are best known for inducing cell death and immune responses in dicotyledonous plants. A group of putative NLP genes from obligate biotrophic oomycete Plasmopara viticola were predicted by RNA-Seq in our previous study, but their activity has not been established. Therefore, we analyzed the P. viticola NLP (PvNLP) family and identified seven PvNLP genes. They all belong to type 1 NLP genes and form a P. viticola-specific cluster when compared with other pathogen NLP genes. The expression of PvNLPs was induced during early infection process and the expression patterns could be categorized into two groups. Agrobacterium tumefaciens-mediated transient expression assays revealed that only PvNLP7 was cytotoxic and could induce Phytophthora capsici resistance in Nicotiana benthamiana. Functional analysis showed that PvNLP4, PvNLP5, PvNLP7, and PvNLP10 significantly improved disease resistance of Arabidopsis thaliana to Hyaloperonospora arabidopsidis. Moreover, the four genes caused an inhibition of plant growth which is typically associated with enhanced immunity when over-expressed in Arabidopsis. Further research found that PvNLP7 could activate the expression of defense-related genes and its conserved NPP1 domain was critical for cell death- and immunity-inducing activity. This record of NLP genes from P. viticola showed a functional diversification, laying a foundation for further study on pathogenic mechanism of the devastating pathogen.

Tandem metalloenzymes gate plant cell entry by pathogenic fungi

Global food security is endangered by fungal phytopathogens causing devastating crop production losses. Many of these pathogens use specialized appressoria cells to puncture plant cuticles. Here, we unveil a pair of alcohol oxidase-peroxidase enzymes to be essential for pathogenicity. Using Colletotrichum orbiculare, we show that the enzyme pair is cosecreted by the fungus early during plant penetration and that single and double mutants have impaired penetration ability. Molecular modeling, biochemical, and biophysical approaches revealed a fine-tuned interplay between these metalloenzymes, which oxidize plant cuticular long-chain alcohols into aldehydes. We show that the enzyme pair is involved in transcriptional regulation of genes necessary for host penetration. The identification of these infection-specific metalloenzymes opens new avenues on the role of wax-derived compounds and the design of oxidase-specific inhibitors for crop protection.

Comparative genomic analysis reveals contraction of gene families with putative roles in pathogenesis in the fungal boxwood pathogens Calonectria henricotiae and C. pseudonaviculata

Background

Boxwood blight disease caused by Calonectria henricotiae and C. pseudonaviculata is of ecological and economic significance in cultivated and native ecosystems worldwide. Prior research has focused on understanding the population genetic and genomic diversity of C. henricotiae and C. pseudonaviculata, but gene family evolution in the context of host adaptation, plant pathogenesis, and trophic lifestyle is poorly understood. This study applied bioinformatic and phylogenetic methods to examine gene family evolution in C. henricotiae, C. pseudonaviculata and 22 related fungi in the Nectriaceae that vary in pathogenic and saprobic (apathogenic) lifestyles.

Results

A total of 19,750 gene families were identified in the 24 genomes, of which 422 were rapidly evolving. Among the six Calonectria species, C. henricotiae and C. pseudonaviculata were the only species to experience high levels of rapid contraction of pathogenesis-related gene families (89% and 78%, respectively). In contrast, saprobic species Calonectria multiphialidica and C. naviculata, two of the closest known relatives of C. henricotiae and C. pseudonaviculata, showed rapid expansion of pathogenesis-related gene families.

Conclusions

Our results provide novel insight into gene family evolution within C. henricotiae and C. pseudonaviculata and suggest gene family contraction may have contributed to limited host-range expansion of these pathogens within the plant family Buxaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02035-4.

Comparative Genomics Reveal the Utilization Ability of Variable Carbohydrates as Key Genetic Features of Listeria Pathogens in Their Pathogenic Lifestyles

BACKGROUND

L. monocytogenes and L. ivanovii, the only two pathogens of Listeria, can survive in various environments, having different pathogenic characteristics. However, the genetic basis of their excellent adaptability and differences in pathogenicity has still not been completely elucidated.

METHODS

We performed a comparative genomic analysis based on 275 L. monocytogenes, 10 L. ivanovii, and 22 non-pathogenic Listeria strains.

RESULTS

Core/pan-genome analysis revealed that 975 gene families were conserved in all the studied strains. Additionally, 204, 242, and 756 gene families existed uniquely in L. monocytogenes, L. ivanovii, and both, respectively. Functional annotation partially verified that these unique gene families were closely related to their adaptability and pathogenicity. Moreover, the protein-protein interaction (PPI) network analysis of these unique gene sets showed that plenty of carbohydrate transport systems and energy metabolism enzymes were clustered in the networks. Interestingly, ethanolamine-metabolic-process-related proteins were significantly enriched in the PPI network of the unique genes of the Listeria pathogens, which can be understood as a determining factor of their pathogenicity.

CONCLUSIONS

The utilization capacity of multiple carbon sources of Listeria pathogens, especially ethanolamine, is the key genetic basis for their ability to adapt to various environments and pathogenic lifestyles.

Pan-Genome Analysis of Campylobacter: Insights on the Genomic Diversity and Virulence Profile

The genus Campylobacter contains pathogens that cause bacterial gastroenteritis in humans and animals. Despite large-scale sequencing efforts to raise clinical awareness of Campylobacter, little is known about the diversity and functions of virulence factors. Here, we constructed the pan-genome of Campylobacter using 39 representative genomes, elucidating their genetic diversity, evolutionary characteristics, and virulence and resistance profiles. The Campylobacter pan-genome was open and showed extensive genome variability, with high levels of gene expansion and contraction as the organism evolved. These Campylobacter members had diverse virulence gene content, and six potential core virulence genes (porA, PEB4, cheY, htrB, Cj1135, and kpsF) have been identified. The conserved mechanisms for Campylobacter pathogenicity were related to adherence, motility, and immune modulation. We emphasized the relative importance of variable virulence genes. Many virulence genes have experienced expansion or contraction in specific lineages, which may be one of the factors causing differences in the content of virulence genes. Additionally, these Campylobacter genomes have a high prevalence of the cmeA and cmeC genes, which are linked to the CmeABC pump and contribute to multidrug resistance. The genomic variations, core and variable virulence factors, and resistance genes of Campylobacter characterized in this study would contribute to a better understanding of the virulence of Campylobacter and more effective use of candidates for drug development and prevention of Campylobacter infections. IMPORTANCE Pathogenic members of the genus Campylobacter are recognized as one of the major causative agents of human bacterial gastroenteritis. This study revealed the pan-genome of 39 Campylobacter species, provided the most updated reconstruction of the global virulence gene pool of 39 Campylobacter species, and identified species-related virulence differences. This study highlighted the basic conserved functionality and specificity of pathogenicity that are crucial to infection, which was critical for improving the diagnosis and prevention of Campylobacter infections.

Integrative functional analysis uncovers metabolic differences between Candida species

Candida species are a dominant constituent of the human mycobiome and associated with the development of several diseases. Understanding the Candida species metabolism could provide key insights into their ability to cause pathogenesis. Here, we have developed the BioFung database, providing an efficient annotation of protein-encoding genes. Along, with BioFung, using carbohydrate-active enzyme (CAZymes) analysis, we have uncovered core and accessory features across Candida species demonstrating plasticity, adaption to the environment and acquired features. We show a greater importance of amino acid metabolism, as functional analysis revealed that all Candida species can employ amino acid metabolism. However, metabolomics revealed that only a specific cluster of species (AGAu species—C. albicans, C. glabrata and C. auris) utilised amino acid metabolism including arginine, cysteine, and methionine metabolism potentially improving their competitive fitness in pathogenesis. We further identified critical metabolic pathways in the AGAu cluster with biomarkers and anti-fungal target potential in the CAZyme profile, polyamine, choline and fatty acid biosynthesis pathways. This study, combining genomic analysis, and validation with gene expression and metabolomics, highlights the metabolic diversity with AGAu species that underlies their remarkable ability to dominate they mycobiome and cause disease.

Metabolic differences between Candida species are uncovered using the BioFung database alongside genomic and metabolic analysis.

Characterization and necrosis-inducing activity of necrosis- and ethylene-inducing peptide 1-like proteins from Colletotrichum australisinense, the causative agent of rubber tree anthracnose

Colletotrichum australisinense, a member of the Colletotrichum acutatum species complex, is an important pathogen causing rubber tree anthracnose. Genome-wide comparative analysis showed this species complex contains more genes encoding necrosis- and ethylene-inducing peptide 1-like proteins (NLPs) than other Colletotrichum species complexes, but little is known about their necrosis-inducing roles in host. The aim of this study was to analyze NLPs number and type in C. australisinense, and characterize their necrosis-inducing activity in host or non-host. According to phylogenetic relationship, conserved the cysteine residues and the heptapeptide motif (GHRHDWE), 11 NLPs were identified and classified into three types. Five of the eleven NLPs were evaluated for necrosis-inducing activity. CaNLP4 (type 1) could not induce necrosis in host or non-host plants. By contrast, both CaNLP5 and CaNLP9 (type 1) induced necrosis in host and non-host plants, and necrosis-inducing activity was strongest for CaNLP9. CaNLP10 (type 2) and CaNLP11 (type 3) induced necrosis in host but not non-host plants. Substitution of key amino acid residues essential for necrosis induction activity led to loss of CaNLP4 activity. Structural characterization of CaNLP5 and CaNLP9 may explain differences in necrosis-inducing activity. We evaluated the expression of genes coding CaNLP by reverse transcription polymerase chain reaction (RT-PCR) and quantitative real-time PCR (qRT-PCR) at different time-points after pathogen infection. It was found that genes encoding CaNLPs with different activities exhibited significantly different expression patterns. The results demonstrate that CaNLPs are functionally and spatially distinct, and may play different but important roles in C. australisinense pathogenesis.

Genome Analyses of Two Blueberry Pathogens: Diaporthe amygdali CAA958 and Diaporthe eres CBS 160.32

The genus Diaporthe includes pathogenic species distributed worldwide and affecting a wide variety of hosts. Diaporthe amygdali and Diaporthe eres have been found to cause cankers, dieback, or twig blights on economically important crops such as soybean, almond, grapevine, and blueberry. Despite their importance as plant pathogens, the strategies of species of Diaporthe to infect host plants are poorly explored. To provide a genomic basis of pathogenicity, the genomes of D. amygdali CAA958 and D. eres CBS 160.32 were sequenced and analyzed. Cellular transporters involved in the transport of toxins, ions, sugars, effectors, and genes implicated in pathogenicity were detected in both genomes. Hydrolases and oxidoreductases were the most prevalent carbohydrate-active enzymes (CAZymes). However, analyses of the secreted proteins revealed that the secretome of D. eres CBS 160.32 is represented by 5.4% of CAZymes, whereas the secreted CAZymes repertoire of D. amygdali CAA958 represents 29.1% of all secretomes. Biosynthetic gene clusters (BGCs) encoding compounds related to phytotoxins and mycotoxins were detected in D. eres and D. amygdali genomes. The core gene clusters of the phytotoxin Fusicoccin A in D. amygdali are reported here through a genome-scale assembly. Comparative analyses of the genomes from 11 Diaporthe species revealed an average of 874 CAZymes, 101 secondary metabolite BGCs, 1640 secreted proteins per species, and genome sizes ranging from 51.5 to 63.6 Mbp. This study offers insights into the overall features and characteristics of Diaporthe genomes. Our findings enrich the knowledge about D. eres and D. amygdali, which will facilitate further research into the pathogenicity mechanisms of these species.

Comparative genomic analysis reveals cellulase plays an important role in the pathogenicity of Setosphaeria turcica f. sp. zeae

Setosphaeria turcica f. sp. zeae and S. turcica f. sp. sorghi, the two formae speciales of S. turcica, cause northern leaf blight disease of corn and sorghum, respectively, and often cause serious economic losses. They have obvious physiological differentiation and show complete host specificity. Host specificity is often closely related to pathogen virulence factors, including secreted protein effectors and secondary metabolites. Genomic sequencing can provide more information for understanding the virulence mechanisms of pathogens. However, the complete genomic sequence of S. turcica f. sp. sorghi has not yet been reported, and no comparative genomic information is available for the two formae speciales. In this study, S. turcica f. sp. zeae was predicted to have fewer secreted proteins, pathogen-host interaction (PHI) genes and carbohydrate-active enzymes (CAZys) than S. turcica f. sp. sorghi. Fifteen and 20 polyketide synthase (PKS) genes were identified in S. turcica f. sp. zeae and S. turcica f. sp. sorghi, respectively, which maintained high homology. There were eight functionally annotated effector protein-encoding genes specifically in S. turcica f. sp. zeae, among which the encoding gene StCEL2 of endo-1, 4-β-D-glucanase, an important component of cellulase, was significantly up-regulated during the interaction process. Finally, gluconolactone inhibited cellulase activity and decreased infection rate and pathogenicity, which indicates that cellulase is essential for maintaining virulence. These findings demonstrate that cellulase plays an important role in the pathogenicity of S. turcica f. sp. zeae. Our results also provide a theoretical basis for future research on the molecular mechanisms underlying the pathogenicity of the two formae speciales and for identifying any associated genes.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

Proteinaceous Effector Discovery and Characterization in Plant Pathogenic Colletotrichum Fungi

Anthracnose caused by plant pathogenic Colletotrichum fungi results in large economic losses in field crop production worldwide. To aid the establishment of plant host infection, Colletotrichum pathogens secrete numerous effector proteins either in apoplastic space or inside of host cells for effective colonization. Understanding these effector repertoires is critical for developing new strategies for resistance breeding and disease management. With the advance of genomics and bioinformatics tools, a large repertoire of putative effectors has been identified in Colletotrichum genomes, and the biological functions and molecular mechanisms of some studied effectors have been summarized. Here, we review recent advances in genomic identification, understanding of evolutional characteristics, transcriptional profiling, and functional characterization of Colletotrichum effectors. We also offer a perspective on future research.

Phylogenetic and ecological reevaluation of the order Onygenales

The order Onygenales is classified in the class Eurotiomycetes of the subphylum Pezizomycotina. Families in this order have classically been isolated from soil and dung, and two lineages contain causative agents of superficial, cutaneous and systemic infections in mammals. The ecology and habitat choices of the species are driven mainly by the keratin and cellulose degradation abilities. The present study aimed to investigate whether the ecological trends of the members of Onygenales can be interpreted in an evolutionary sense, linking phylogenetic parameters with habitat preferences, to achieve polyphasic definitions of the main taxonomic groups. Evolutionary processes were estimated by multiple gene genealogies and divergence time analysis. Previously described families, namely, Arthrodermataceae, Ajellomycetaceae, Ascosphaeraceae, Eremascaceae, Gymnoascaceae, Onygenaceae and Spiromastigoidaceae, were accepted in Onygenales, and two new families, Malbrancheaceae and Neogymnomycetaceae, were introduced. A number of species could not be assigned to any of the defined families. Our study provides a revised overview of the main lines of taxonomy of Onygenales, supported by multilocus analyses of ITS, LSU, TUB, TEF1, TEF3, RPB1, RPB2, and ribosomal protein 60S L10 (L1) (RP60S) sequences, combined with available data on ecology, physiology, morphology, and genomics.

Comparative Genomics of Three Colletotrichum scovillei Strains and Genetic Analysis Revealed Genes Involved in Fungal Growth and Virulence on Chili Pepper

Colletotrichum scovillei causes anthracnose of chili pepper in many countries. Three strains of this pathogen, Coll-524, Coll-153, and Coll-365, show varied virulence on chili pepper. Among the three strains, Coll-365 showed significant defects in growth and virulence. To decipher the genetic variations among these strains and identify genes contributing to growth and virulence, comparative genomic analysis and gene transformation to show gene function were applied in this study. Compared to Coll-524, Coll-153, and Coll-365 had numerous gene losses including 32 candidate effector genes that are mainly exist in acutatum species complex. A cluster of 14 genes in a 34-kb genomic fragment was lost in Coll-365. Through gene transformation, three genes in the 34-kb fragment were identified to have functions in growth and/or virulence of C. scovillei. CsPLAA encoding a phospholipase A2-activating protein enhanced the growth of Coll-365. A combination of CsPLAA with one transcription factor CsBZTF and one C6 zinc finger domain-containing protein CsCZCP was found to enhance the pathogenicity of Coll-365. Introduction of CsGIP, which encodes a hypothetical protein, into Coll-365 caused a reduction in the germination rate of Coll-365. In conclusion, the highest virulent strain Coll-524 had more genes and encoded more pathogenicity related proteins and transposable elements than the other two strains, which may contribute to the high virulence of Coll-524. In addition, the absence of the 34-kb fragment plays a critical role in the defects of growth and virulence of strain Coll-365.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.

Fungi as mediators linking organisms and ecosystems

Fungi form a major and diverse component of most ecosystems on Earth. They are both micro and macroorganisms with high and varying functional diversity as well as great variation in dispersal modes. With our growing knowledge of microbial biogeography, it has become increasingly clear that fungal assembly patterns and processes differ from other microorganisms such as bacteria, but also from macroorganisms such as plants. The success of fungi as organisms and their influence on the environment lies in their ability to span multiple dimensions of time, space, and biological interactions, that is not rivalled by other organism groups. There is also growing evidence that fungi mediate links between different organisms and ecosystems, with the potential to affect the macroecology and evolution of those organisms. This suggests that fungal interactions are an ecological driving force, interconnecting different levels of biological and ecological organisation of their hosts, competitors, and antagonists with the environment and ecosystem functioning. Here we review these emerging lines of evidence by focusing on the dynamics of fungal interactions with other organism groups across various ecosystems. We conclude that the mediating role of fungi through their complex and dynamic ecological interactions underlie their importance and ubiquity across Earth's ecosystems.

Complete Genome Sequence of the Plant-Pathogenic Fungus Colletotrichum lupini

Colletotrichum is a fungal genus (Ascomycota, Sordariomycetes, Glomerellaceae) that includes many economically important plant pathogens that cause devastating diseases of a wide range of plants. In this work, using a combination of long- and short-read sequencing technologies, we sequenced the genome of Colletotrichum lupini RB221, isolated from white lupin (Lupinus albus) in France during a survey in 2014. The genome was assembled into 11 nuclear chromosomes and a mitochondrial genome with a total assembly size of 63.41 Mb and 36.55 kb, respectively. In total, 18,324 protein-encoding genes have been predicted, of which only 39 are specific to C. lupini. This resource will provide insight into pathogenicity factors and will help provide a better understanding of the evolution and genome structure of this important plant pathogen.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative Analysis of Herbaceous and Woody Cell Wall Digestibility by Pathogenic Fungi

Fungal pathogens have evolved combinations of plant cell-wall-degrading enzymes (PCWDEs) to deconstruct host plant cell walls (PCWs). An understanding of this process is hoped to create a basis for improving plant biomass conversion efficiency into sustainable biofuels and bioproducts. Here, an approach integrating enzyme activity assay, biomass pretreatment, field emission scanning electron microscopy (FESEM), and genomic analysis of PCWDEs were applied to examine digestibility or degradability of selected woody and herbaceous biomass by pathogenic fungi. Preferred hydrolysis of apple tree branch, rapeseed straw, or wheat straw were observed by the apple-tree-specific pathogen Valsa mali, the rapeseed pathogen Sclerotinia sclerotiorum, and the wheat pathogen Rhizoctonia cerealis, respectively. Delignification by peracetic acid (PAA) pretreatment increased PCW digestibility, and the increase was generally more profound with non-host than host PCW substrates. Hemicellulase pretreatment slightly reduced or had no effect on hemicellulose content in the PCW substrates tested; however, the pretreatment significantly changed hydrolytic preferences of the selected pathogens, indicating a role of hemicellulose branching in PCW digestibility. Cellulose organization appears to also impact digestibility of host PCWs, as reflected by differences in cellulose microfibril organization in woody and herbaceous PCWs and variation in cellulose-binding domain organization in cellulases of pathogenic fungi, which is known to influence enzyme access to cellulose. Taken together, this study highlighted the importance of chemical structure of both hemicelluloses and cellulose in host PCW digestibility by fungal pathogens.

Identification of Copper-Containing Oxidoreductases in the Secretomes of Three Colletotrichum Species with a Focus on Copper Radical Oxidases for the Biocatalytic Production of Fatty Aldehydes

Copper radical alcohol oxidases (CRO-AlcOx), which have been recently discovered among fungal phytopathogens, are attractive for the production of fragrant fatty aldehydes. With the initial objective to investigate the secretion of CRO-AlcOx by natural fungal strains, we undertook time course analyses of the secretomes of three Colletotrichum species (C. graminicola, C. tabacum, and C. destructivum) using proteomics. The addition of a copper-manganese-ethanol mixture in the absence of any plant-biomass mimicking compounds to Colletotrichum cultures unexpectedly induced the secretion of up to 400 proteins, 29 to 52% of which were carbohydrate-active enzymes (CAZymes), including a wide diversity of copper-containing oxidoreductases from the auxiliary activities (AA) class (AA1, AA3, AA5, AA7, AA9, AA11, AA12, AA13, and AA16). Under these specific conditions, while a CRO-glyoxal oxidase from the AA5_1 subfamily was among the most abundantly secreted proteins, the targeted AA5_2 CRO-AlcOx were secreted at lower levels, suggesting heterologous expression as a more promising strategy for CRO-AlcOx production and utilization. C. tabacum and C. destructivum CRO-AlcOx were thus expressed in Pichia pastoris, and their preference toward both aromatic and aliphatic primary alcohols was assessed. The CRO-AlcOx from C. destructivum was further investigated in applied settings, revealing a full conversion of C₆ and C₈ alcohols into their corresponding fragrant aldehydes. IMPORTANCE In the context of the industrial shift toward greener processes, the biocatalytic production of aldehydes is of utmost interest owing to their importance for their use as flavor and fragrance ingredients. Copper radical alcohol oxidases (CRO-AlcOx) have the potential to become platform enzymes for the oxidation of alcohols to aldehydes. However, the secretion of CRO-AlcOx by natural fungal strains has never been explored, while the use of crude fungal secretomes is an appealing approach for industrial applications to alleviate various costs pertaining to biocatalyst production. While investigating this primary objective, the secretomics studies revealed unexpected results showing that under the oxidative stress conditions we probed, Colletotrichum species can secrete a broad diversity of copper-containing enzymes (laccases, sugar oxidoreductases, and lytic polysaccharide monooxygenases [LPMOs]) usually assigned to "plant cell wall degradation," despite the absence of any plant-biomass mimicking compound. However, in these conditions, only small amounts of CRO-AlcOx were secreted, pointing out recombinant expression as the most promising path for their biocatalytic application.

Identification and Comparison of Colletotrichum Secreted Effector Candidates Reveal Two Independent Lineages Pathogenic to Soybean

Colletotrichum is one of the most important plant pathogenic genus of fungi due to its scientific and economic impact. A wide range of hosts can be infected by Colletotrichum spp., which causes losses in crops of major importance worldwide, such as soybean. Soybean anthracnose is mainly caused by C. truncatum, but other species have been identified at an increasing rate during the last decade, becoming one of the most important limiting factors to soybean production in several regions. To gain a better understanding of the evolutionary origin of soybean anthracnose, we compared the repertoire of effector candidates of four Colletotrichum species pathogenic to soybean and eight species not pathogenic. Our results show that the four species infecting soybean belong to two lineages and do not share any effector candidates. These results strongly suggest that two Colletotrichum lineages have acquired the capability to infect soybean independently. This study also provides, for each lineage, a set of candidate effectors encoding genes that may have important roles in pathogenicity towards soybean offering a new resource useful for further research on soybean anthracnose management.

FunOrder: A robust and semi-automated method for the identification of essential biosynthetic genes through computational molecular co-evolution

Secondary metabolites (SMs) are a vast group of compounds with different structures and properties that have been utilized as drugs, food additives, dyes, and as monomers for novel plastics. In many cases, the biosynthesis of SMs is catalysed by enzymes whose corresponding genes are co-localized in the genome in biosynthetic gene clusters (BGCs). Notably, BGCs may contain so-called gap genes, that are not involved in the biosynthesis of the SM. Current genome mining tools can identify BGCs, but they have problems with distinguishing essential genes from gap genes. This can and must be done by expensive, laborious, and time-consuming comparative genomic approaches or transcriptome analyses. In this study, we developed a method that allows semi-automated identification of essential genes in a BGC based on co-evolution analysis. To this end, the protein sequences of a BGC are blasted against a suitable proteome database. For each protein, a phylogenetic tree is created. The trees are compared by treeKO to detect co-evolution. The results of this comparison are visualized in different output formats, which are compared visually. Our results suggest that co-evolution is commonly occurring within BGCs, albeit not all, and that especially those genes that encode for enzymes of the biosynthetic pathway are co-evolutionary linked and can be identified with FunOrder. In light of the growing number of genomic data available, this will contribute to the studies of BGCs in native hosts and facilitate heterologous expression in other organisms with the aim of the discovery of novel SMs.

The discovery and description of novel fungal secondary metabolites promises novel antibiotics, pharmaceuticals, and other useful compounds. A way to identify novel secondary metabolites is to express the corresponding genes in a suitable expression host. Consequently, a detailed knowledge or an accurate prediction of these genes is necessary. In fungi, the genes are co-localized in so-called biosynthetic gene clusters. Notably, the clusters may also contain genes that are not necessary for the biosynthesis of the secondary metabolites, so-called gap genes. We developed a method to detect co-evolved genes within the clusters and demonstrated that essential genes are co-evolving and can thus be differentiated from the gap genes. This adds an additional layer of information, which can support researchers with their decisions on which genes to study and express for the discovery of novel secondary metabolites.

The Conserved Colletotrichum spp. Effector Candidate CEC3 Induces Nuclear Expansion and Cell Death in Plants

Plant pathogens secrete proteins, known as effectors, that promote infection by manipulating host cells. Members of the phytopathogenic fungal genus Colletotrichum collectively have a broad host range and generally adopt a hemibiotrophic lifestyle that includes an initial biotrophic phase and a later necrotrophic phase. We hypothesized that Colletotrichum fungi use a set of conserved effectors during infection to support the two phases of their hemibiotrophic lifestyle. This study aimed to examine this hypothesis by identifying and characterizing conserved effectors among Colletotrichum fungi. Comparative genomic analyses using genomes of ascomycete fungi with different lifestyles identified seven effector candidates that are conserved across the genus Colletotrichum. Transient expression assays showed that one of these putative conserved effectors, CEC3, induces nuclear expansion and cell death in Nicotiana benthamiana, suggesting that CEC3 is involved in promoting host cell death during infection. Nuclear expansion and cell death induction were commonly observed in CEC3 homologs from four different Colletotrichum species that vary in host specificity. Thus, CEC3 proteins could represent a novel class of core effectors with functional conservation in the genus Colletotrichum.

Homeobox Transcription Factors Are Required for Fungal Development and the Suppression of Host Defense Mechanisms in the Colletotrichum scovillei-Pepper Pathosystem

Colletotrichum scovillei, an ascomycete phytopathogenic fungus, is the main causal agent of serious yield losses of economic crops worldwide. The fungus causes anthracnose disease on several fruits, including peppers. However, little is known regarding the underlying molecular mechanisms involved in the development of anthracnose caused by this fungus. In an initial step toward understanding the development of anthracnose on pepper fruits, we retrieved 624 transcription factors (TFs) from the whole genome of C. scovillei and comparatively analyzed the entire repertoire of TFs among phytopathogenic fungi. Evolution and proliferation of members of the homeobox-like superfamily, including homeobox (HOX) TFs that regulate the development of eukaryotic organisms, were demonstrated in the genus Colletotrichum. C. scovillei was found to contain 10 HOX TF genes (CsHOX1 to CsHOX10), which were functionally characterized using deletion mutants of each CsHOX gene. Notably, CsHOX1 was identified as a pathogenicity factor required for the suppression of host defense mechanisms, which represents a new role for HOX TFs in pathogenic fungi. CsHOX2 and CsHOX7 were found to play essential roles in conidiation and appressorium development, respectively, in a stage-specific manner in C. scovillei. Our study provides a molecular basis for understanding the mechanisms associated with the development of anthracnose on fruits caused by C. scovillei, which will aid in the development of novel approaches for disease management. IMPORTANCE The ascomycete phytopathogenic fungus, Colletotrichum scovillei, causes serious yield loss on peppers. However, little is known about molecular mechanisms involved in the development of anthracnose caused by this fungus. We analyzed whole-genome sequences of C. scovillei and isolated 624 putative TFs, revealing the existence of 10 homeobox (HOX) transcription factor (TF) genes. We found that CsHOX1 is a pathogenicity factor required for the suppression of host defense mechanism, which represents a new role for HOX TFs in pathogenic fungi. We also found that CsHOX2 and CsHOX7 play essential roles in conidiation and appressorium development, respectively, in a stage-specific manner in C. scovillei. Our study contributes to understanding the mechanisms associated with the development of anthracnose on fruits caused by C. scovillei, which will aid for initiating novel approaches for disease management.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07902-w.

Sequence and Gene Expression Analysis of Recently Identified NLP from Plasmopara viticola

Grapevine downy mildew, evoked by the obligate biotrophic oomycete Plasmopara viticola, is one of the most challenging diseases in viticulture. P. viticola establishes an infection by circumvention of plant immunity, which is achieved by the secretion of effector molecules. One family of potential effectors are the necrosis- and ethylene-inducing peptide 1 (Nep1)-like proteins (NLP). NLP are most abundant in plant pathogenic microorganisms and exist in cytotoxic and non-cyctotoxic forms. Cytotoxic NLP often act as virulence factors and are synthesized in necrotrophic or hemibiotrophic pathogens during the transition from biotrophic to necrotrophic growth. In addition to these cytotoxic NLP, many non-cytotoxic NLP have been identified; their function in biotrophic pathogens is still unknown. In 2020, eight different NLP coding genes were identified in P. viticola and named PvNLP1 to PvNLP8 (Plasmopara viticolaNLP 1–8). In the present study, PvNLP4 to PvNLP8 were characterized by using qPCR analysis and transient expression in the model plant Nicotiana benthamiana. Gene expression analysis showed high PvNLP expression during the early stages of infection. Necrosis-inducing activity of PvNLP was not observed in the nonhost N. benthamiana.

The Histone Acetyltransferase CfGcn5 Regulates Growth, Development, and Pathogenicity in the Anthracnose Fungus Colletotrichum fructicola on the Tea-Oil Tree

The tea-oil tree (Camellia oleifera Abel.) is a commercial edible-oil tree in China, and anthracnose commonly occurs in its plantations, causing great losses annually. We have previously revealed that CfSnf1 is essential for pathogenicity in Colletotrichum fructicola, the major pathogen of anthracnose on the tea-oil tree. Here, we identified CfGcn5 as the homolog of yeast histone acetyltransferase ScGcn5, which cooperates with ScSnf1 to modify histone H3 in Saccharomyces cerevisiae. Targeted gene deletion revealed that CfGcn5 is important in fungi growth, conidiation, and responses to environmental stresses. Pathogenicity assays indicated that CfGcn5 is essential for C. fructicola virulence both in unwounded and wounded tea-oil tree leaves. Further, we found that CfGcn5 is localized to the nucleus and this specific localization is dependent on both NLS region and HAT domain. Moreover, we provided evidence showing that the nuclear localization is essential but not sufficient for the full function of CfGcn5, and the NLS, HAT, and Bromo domains were proven to be important for normal CfGcn5 functions. Taken together, our studies not only illustrate the key functions of CfGcn5 in growth, development, and pathogenicity but also highlight the relationship between its locations with functions in C. fructicola.

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.

Comparative genomics reveals dynamic genome evolution in host specialist ectomycorrhizal fungi

While there has been significant progress characterizing the 'symbiotic toolkit' of ectomycorrhizal (ECM) fungi, how host specificity may be encoded into ECM fungal genomes remains poorly understood. We conducted a comparative genomic analysis of ECM fungal host specialists and generalists, focusing on the specialist genus Suillus. Global analyses of genome dynamics across 46 species were assessed, along with targeted analyses of three classes of molecules previously identified as important determinants of host specificity: small secreted proteins (SSPs), secondary metabolites (SMs) and G-protein coupled receptors (GPCRs). Relative to other ECM fungi, including other host specialists, Suillus had highly dynamic genomes including numerous rapidly evolving gene families and many domain expansions and contractions. Targeted analyses supported a role for SMs but not SSPs or GPCRs in Suillus host specificity. Phylogenomic-based ancestral state reconstruction identified Larix as the ancestral host of Suillus, with multiple independent switches between white and red pine hosts. These results suggest that like other defining characteristics of the ECM lifestyle, host specificity is a dynamic process at the genome level. In the case of Suillus, both SMs and pathways involved in the deactivation of reactive oxygen species appear to be strongly associated with enhanced host specificity.

Diversity and Cross-Infection Potential of Colletotrichum Causing Fruit Rots in Mixed-Fruit Orchards in Kentucky

Fungi in the genus Colletotrichum cause apple, blueberry, and strawberry fruit rots, which can result in significant losses. Accurate identification is important because species differ in aggressiveness, fungicide sensitivity, and other factors affecting management. Multiple Colletotrichum species can cause similar symptoms on the same host, and more than one fruit type can be infected by a single Colletotrichum species. Mixed-fruit orchards may facilitate cross-infection, with significant management implications. Colletotrichum isolates from small fruits in Kentucky orchards were characterized and compared with apple isolates via a combination of morphotyping, sequencing of voucher loci and whole genomes, and cross-inoculation assays. Seven morphotypes representing two species complexes (C. acutatum and C. gloeosporioides) were identified. Morphotypes corresponded with phylogenetic species C. fioriniae, C. fructicola, C. nymphaeae, and C. siamense, identified by TUB2 or GAPDH barcodes. Phylogenetic trees built from nine single-gene sequences matched barcoding results with one exception, later determined to belong to an undescribed species. Comparison of single-gene trees with representative whole genome sequences revealed that CHS and ApMat were the most informative for diagnosis of fruit rot species and individual morphotypes within the C. acutatum or C. gloeosporioides complexes, respectively. All blueberry isolates belonged to C. fioriniae, and most strawberry isolates were C. nymphaeae, with a few C. siamense and C. fioriniae also recovered. All three species cause fruit rot on apples in Kentucky. Cross-inoculation assays on detached apple, blueberry, and strawberry fruits showed that all species were pathogenic on all three hosts but with species-specific differences in aggressiveness.

CgNPG1 as a Novel Pathogenic Gene of Colletotrichum gloeosporioides From Hevea brasiliensis in Mycelial Growth, Conidiation, and the Invasive Structures Development

The rubber tree (Hevea brasiliensis) is a tropical perennial crop for the primary source of natural rubber. Colletotrichum gloeosporioides from Hevea brasiliensis (C. gloeosporioides Hb) and Colletotrichum acutatum from Hevea brasiliensis (C. acutatum Hb) are the causal agents of rubber tree anthracnose and lead to serious loss of natural rubber production. Inoculation tests showed that C. gloeosporioides Hb possessed higher pathogenicity than C. acutatum Hb to the rubber tree. Genomic analysis revealed that an unknown gene, named CgNPG1 (a Novel Pathogenic Gene 1), was presented in the genome of C. gloeosporioides Hb but not identified in C. acutatum Hb. CgNPG1 was predicted to encode a small secretory protein without any conserved domain. To investigate the functions of CgNPG1 in C. gloeosporioides Hb and in C. acutatum Hb, the gene deletion and overexpression mutants were generated. The phenotype analysis showed that deletion of CgNPG1 led to changed conidia morphology, decreased mycelial growth, conidiation, conidia germination rate, appressorium formation rate, and pathogenicity of C. gloeosporioides Hb to the rubber tree. Meanwhile, heterogeneous expression of CgNPG1 in C. acutatum Hb significantly changed the conidia morphology and improved the mycelial growth rate, conidiation, conidia germination rate, appressorium formation rate, and the pathogenicity of C. acutatum Hb to the rubber tree. Consistently, CgNPG1 increased the expression level of CaCRZ1 and CaCMK1 in C. acutatum Hb. These data suggested that CgNPG1 contributed to mycelial growth, conidiation, the development of invasive structures, and the pathogenicity of Colletotrichum to the rubber tree, which might be related to the modulation of CaCRZ1 and mitogen-activated protein kinase CMK1. Our results provided new insight into CgNPG1 in regulating growth and pathogenicity of the Colletotrichum spp.

Phylogenomics of Plant-Associated Botryosphaeriaceae Species

The Botryosphaeriaceae is a fungal family that includes many destructive vascular pathogens of woody plants (e.g., Botryosphaeria dieback of grape, Panicle blight of pistachio). Species in the genera Botryosphaeria, Diplodia, Dothiorella, Lasiodiplodia, Neofusicoccum, and Neoscytalidium attack a range of horticultural crops, but they vary in virulence and their abilities to infect their hosts via different infection courts (flowers, green shoots, woody twigs). Isolates of seventeen species, originating from symptomatic apricot, grape, pistachio, and walnut were tested for pathogenicity on grapevine wood after 4 months of incubation in potted plants in the greenhouse. Results revealed significant variation in virulence in terms of the length of the internal wood lesions caused by these seventeen species. Phylogenomic comparisons of the seventeen species of wood-colonizing fungi revealed clade-specific expansion of gene families representing putative virulence factors involved in toxin production and mobilization, wood degradation, and nutrient uptake. Statistical analyses of the evolution of the size of gene families revealed expansions of secondary metabolism and transporter gene families in Lasiodiplodia and of secreted cell wall degrading enzymes (CAZymes) in Botryosphaeria and Neofusicoccum genomes. In contrast, Diplodia, Dothiorella, and Neoscytalidium generally showed a contraction in the number of members of these gene families. Overall, species with expansions of gene families, such as secreted CAZymes, secondary metabolism, and transporters, were the most virulent (i.e., were associated with the largest lesions), based on our pathogenicity tests and published reports. This study represents the first comparative phylogenomic investigation into the evolution of possible virulence factors from diverse, cosmopolitan members of the Botryosphaeriaceae.

Systematic analysis of gut microbiome reveals the role of bacterial folate and homocysteine metabolism in Parkinson's disease

Parkinson's disease (PD) is the most common progressive neurological disorder compromising motor functions. However, nonmotor symptoms, such as gastrointestinal (GI) dysfunction, precede those affecting movement. Evidence of an early involvement of the GI tract and enteric nervous system highlights the need for better understanding of the role of gut microbiota in GI complications in PD. Here, we investigate the gut microbiome of patients with PD using metagenomics and serum metabolomics. We integrate these data using metabolic modeling and construct an integrative correlation network giving insight into key microbial species linked with disease severity, GI dysfunction, and age of patients with PD. Functional analysis reveals an increased microbial capability to degrade mucin and host glycans in PD. Personalized community-level metabolic modeling reveals the microbial contribution to folate deficiency and hyperhomocysteinemia observed in patients with PD. The metabolic modeling approach could be applied to uncover gut microbial metabolic contributions to PD pathophysiology.