Rust expression browser: an open source database for simultaneous analysis of host and pathogen gene expression profiles with expVIP

New Approaches to Plant Pathogen Detection and Disease Diagnosis

Detecting plant pathogens and diagnosing diseases are critical components of successful pest management. These key areas have undergone significant advancements driven by breakthroughs in molecular biology and remote sensing technologies within the realm of precision agriculture. Notably, nucleic acid amplification techniques, with recent emphasis on sequencing procedures, particularly next-generation sequencing, have enabled improved DNA or RNA amplification detection protocols that now enable previously unthinkable strategies aimed at dissecting plant microbiota, including the disease-causing components. Simultaneously, the domain of remote sensing has seen the emergence of cutting-edge imaging sensor technologies and the integration of powerful computational tools, such as machine learning. These innovations enable spectral analysis of foliar symptoms and specific pathogen-induced alterations, making imaging spectroscopy and thermal imaging fundamental tools for large-scale disease surveillance and monitoring. These technologies contribute significantly to understanding the temporal and spatial dynamics of plant diseases.

A comprehensive workflow for optimizing RNA-seq data analysis

BACKGROUND

Current RNA-seq analysis software for RNA-seq data tends to use similar parameters across different species without considering species-specific differences. However, the suitability and accuracy of these tools may vary when analyzing data from different species, such as humans, animals, plants, fungi, and bacteria. For most laboratory researchers lacking a background in information science, determining how to construct an analysis workflow that meets their specific needs from the array of complex analytical tools available poses a significant challenge.

RESULTS

By utilizing RNA-seq data from plants, animals, and fungi, it was observed that different analytical tools demonstrate some variations in performance when applied to different species. A comprehensive experiment was conducted specifically for analyzing plant pathogenic fungal data, focusing on differential gene analysis as the ultimate goal. In this study, 288 pipelines using different tools were applied to analyze five fungal RNA-seq datasets, and the performance of their results was evaluated based on simulation. This led to the establishment of a relatively universal and superior fungal RNA-seq analysis pipeline that can serve as a reference, and certain standards for selecting analysis tools were derived for reference. Additionally, we compared various tools for alternative splicing analysis. The results based on simulated data indicated that rMATS remained the optimal choice, although consideration could be given to supplementing with tools such as SpliceWiz.

CONCLUSION

The experimental results demonstrate that, in comparison to the default software parameter configurations, the analysis combination results after tuning can provide more accurate biological insights. It is beneficial to carefully select suitable analysis software based on the data, rather than indiscriminately choosing tools, in order to achieve high-quality analysis results more efficiently.

Co-expression network analysis and identification of core genes in the interaction between wheat and Puccinia striiformis f. sp. tritici

The epidemic of stripe rust, caused by the pathogen Puccinia striiformis f. sp. tritici (Pst), would reduce wheat (Triticum aestivum) yields seriously. Traditional experimental methods are difficult to discover the interaction between wheat and Pst. Multi-omics data analysis provides a new idea for efficiently mining the interactions between host and pathogen. We used 140 wheat-Pst RNA-Seq data to screen for differentially expressed genes (DEGs) between low susceptibility and high susceptibility samples, and carried out Gene Ontology (GO) enrichment analysis. Based on this, we constructed a gene co-expression network, identified the core genes and interacted gene pairs from the conservative modules. Finally, we checked the distribution of Nucleotide-binding and leucine-rich repeat (NLR) genes in the co-expression network and drew the wheat NLR gene co-expression network. In order to provide accessible information for related researchers, we built a web-based visualization platform to display the data. Based on the analysis, we found that resistance-related genes such as TaPR1, TaWRKY18 and HSP70 were highly expressed in the network. They were likely to be involved in the biological processes of Pst infecting wheat. This study can assist scholars in conducting studies on the pathogenesis and help to advance the investigation of wheat-Pst interaction patterns.

Transcriptome-Wide N6-Methyladenosine (m6A) Methylation Analyses in a Compatible Wheat-Puccinia striiformis f. sp. tritici Interaction

N6-methyladenosine (m6A) is a prevalent internal modification in eukaryotic mRNA, tRNA, miRNA, and long non-coding RNA. It is also known for its role in plant responses to biotic and abiotic stresses. However, a comprehensive m6A transcriptome-wide map for Puccinia striiformis f. sp. tritici (Pst) infections in wheat (Triticum aestivum) is currently unavailable. Our study is the first to profile m6A modifications in wheat infected with a virulent Pst race. Analysis of RNA-seq and MeRIP-seq data revealed that the majority of differentially expressed genes are up-regulated and hyper-methylated. Some of these genes are enriched in the plant-pathogen interaction pathway. Notably, genes related to photosynthesis showed significant down-regulation and hypo-methylation, suggesting a potential mechanism facilitating successful Pst invasion by impairing photosynthetic function. The crucial genes, epitomizing the core molecular constituents that fortify plants against pathogenic assaults, were detected with varying expression and methylation levels, together with a newly identified methylation motif. Additionally, m6A regulator genes were also influenced by m6A modification, and their expression patterns varied at different time points of post-inoculation, with lower expression at early stages of infection. This study provides insights into the role of m6A modification regulation in wheat's response to Pst infection, establishing a foundation for understanding the potential function of m6A RNA methylation in plant resistance or susceptibility to pathogens.

Genome-Wide Characterization of Berberine Bridge Enzyme Gene Family in Wheat (Triticum aestivum L.) and the Positive Regulatory Role of TaBBE64 in Response to Wheat Stripe Rust

Berberine bridge enzymes (BBEs), functioning as carbonate oxidases, enhance disease resistance in Arabidopsis and tobacco. However, the understanding of BBEs' role in monocots against pathogens remains limited. This study identified 81 TaBBEs with FAD_binding_4 and BBE domains. Phylogenetic analysis revealed a separation of the BBE gene family between monocots and dicots. Notably, RNA-seq showed TaBBE64's significant induction by both pathogen-associated molecular pattern treatment and Puccinia striiformis f. sp. tritici (Pst) infection at early stages. Subcellular localization revealed TaBBE64 at the cytoplasmic membrane. Knocking down TaBBE64 compromised wheat's Pst resistance, reducing reactive oxygen species and promoting fungal growth, confirming its positive role. Molecular docking and enzyme activity assays confirmed TaBBE64's glucose oxidation to produce H2O2. Since Pst relies on living wheat cells for carbohydrate absorption, TaBBE64's promotion of glucose oxidation limits fungal growth and resists pathogen infection. This study sheds light on BBEs' role in wheat resistance against biotrophic fungi.

Deciphering key genomic regions controlling flag leaf size in wheat via integration of meta-QTL and in silico transcriptome assessment

BACKGROUND

Grain yield is a complex and polygenic trait influenced by the photosynthetic source-sink relationship in wheat. The top three leaves, especially the flag leaf, are considered the major sources of photo-assimilates accumulated in the grain. Determination of significant genomic regions and candidate genes affecting flag leaf size can be used in breeding for grain yield improvement.

RESULTS

With the final purpose of understanding key genomic regions for flag leaf size, a meta-analysis of 521 initial quantitative trait loci (QTLs) from 31 independent QTL mapping studies over the past decades was performed, where 333 loci eventually were refined into 64 meta-QTLs (MQTLs). The average confidence interval (CI) of these MQTLs was 5.28 times less than that of the initial QTLs. Thirty-three MQTLs overlapped the marker trait associations (MTAs) previously reported in genome-wide association studies (GWAS) for flag leaf traits in wheat. A total of 2262 candidate genes for flag leaf size, which were involved in the peroxisome, basal transcription factor, and tyrosine metabolism pathways were identified in MQTL regions by the in silico transcriptome assessment. Of these, the expression analysis of the available genes revealed that 134 genes with > 2 transcripts per million (TPM) were highly and specifically expressed in the leaf. These candidate genes could be critical to affect flag leaf size in wheat.

CONCLUSIONS

The findings will make further insight into the genetic determinants of flag leaf size and provide some reliable MQTLs and putative candidate genes for the genetic improvement of flag leaf size in wheat.

The Botrytis cinerea Gene Expression Browser

For comprehensive gene expression analyses of the phytopathogenic fungus Botrytis cinerea, which infects a number of plant taxa and is a cause of substantial agricultural losses worldwide, we developed BEB, a web-based B. cinerea gene Expression Browser. This computationally inexpensive web-based application and its associated database contain manually curated RNA-Seq data for B. cinerea. BEB enables expression analyses of genes of interest under different culture conditions by providing publication-ready heatmaps depicting transcript levels, without requiring advanced computational skills. BEB also provides details of each experiment and user-defined gene expression clustering and visualization options. If needed, tables of gene expression values can be downloaded for further exploration, including, for instance, the determination of differentially expressed genes. The BEB implementation is based on open-source computational technologies that can be deployed for other organisms. In this case, the new implementation will be limited only by the number of transcriptomic experiments that are incorporated into the platform. To demonstrate the usability and value of BEB, we analyzed gene expression patterns across different conditions, with a focus on secondary metabolite gene clusters, chromosome-wide gene expression, previously described virulence factors, and reference genes, providing the first comprehensive expression overview of these groups of genes in this relevant fungal phytopathogen. We expect this tool to be broadly useful in B. cinerea research, providing a basis for comparative transcriptomics and candidate gene identification for functional assays.

Puccinia striiformis f. sp. tritici effectors in wheat immune responses

The obligate biotrophic fungus Puccinia striiformis f. sp. tritici, which causes yellow (stripe) rust disease, is among the leading biological agents resulting in tremendous yield losses on global wheat productions per annum. The combatting strategies include, but are not limited to, fungicide applications and the development of resistant cultivars. However, evolutionary pressure drives rapid changes, especially in its "effectorome" repertoire, thus allowing pathogens to evade and breach resistance. The extracellular and intracellular effectors, predominantly secreted proteins, are tactical arsenals aiming for many defense processes of plants. Hence, the identity of the effectors and the molecular mechanisms of the interactions between the effectors and the plant immune system have long been targeted in research. The obligate biotrophic nature of P. striiformis f. sp. tritici and the challenging nature of its host, the wheat, impede research on this topic. Next-generation sequencing and novel prediction algorithms in bioinformatics, which are accompanied by in vitro and in vivo validation approaches, offer a speedy pace for the discovery of new effectors and investigations of their biological functions. Here, we briefly review recent findings exploring the roles of P. striiformis f. sp. tritici effectors together with their cellular/subcellular localizations, host responses, and interactors. The current status and the challenges will be discussed. We hope that the overall work will provide a broader view of where we stand and a reference point to compare and evaluate new findings.

A decade after the first Pucciniales genomes: A bibliometric snapshot of (post) genomics studies in three model rust fungi

Pucciniales (rust fungi) are one of the largest fungal order of plant pathogens. They collectively infect key crops such as wheat and soybean, and threaten global food security. In the early 2010s, the genome sequences of three rust fungi were released: Melampsora larici-populina (the poplar leaf rust fungus), Puccinia graminis f. sp. tritici (the wheat stem rust fungus), and Puccinia striiformis f. sp. triciti (the wheat stripe rust or wheat yellow rust fungus). The availability of those genomes has forwarded rust biology into the post-genomic era, sparking a series of genomics, transcriptomics, in silico, and functional studies. Here, we snapshot the last 10 years of post-genomics studies addressing M. larici-populina, P. graminis f. sp. tritici, and/or P. striiformis f. sp. tritici. This mini-review notably reveals the model species-centered structure of the research community, and highlights the drastic increase of the number of functional studies focused on effectors since 2014, which notably revealed chloroplasts as a central host compartment targeted by rust fungi. This mini-review also discusses genomics-facilitated studies in other rust species, and emerging post-genomic research trends related to fully-phased rust genomes.

GrainGenes: a data-rich repository for small grains genetics and genomics

As one of the US Department of Agriculture-Agricultural Research Service flagship databases, GrainGenes (https://wheat.pw.usda.gov) serves the data and community needs of globally distributed small grains researchers for the genetic improvement of the Triticeae family and Avena species that include wheat, barley, rye and oat. GrainGenes accomplishes its mission by continually enriching its cross-linked data content following the findable, accessible, interoperable and reusable principles, enhancing and maintaining an intuitive web interface, creating tools to enable easy data access and establishing data connections within and between GrainGenes and other biological databases to facilitate knowledge discovery. GrainGenes operates within the biological database community, collaborates with curators and genome sequencing groups and contributes to the AgBioData Consortium and the International Wheat Initiative through the Wheat Information System (WheatIS). Interactive and linked content is paramount for successful biological databases and GrainGenes now has 2917 manually curated gene records, including 289 genes and 254 alleles from the Wheat Gene Catalogue (WGC). There are >4.8 million gene models in 51 genome browser assemblies, 6273 quantitative trait loci and >1.4 million genetic loci on 4756 genetic and physical maps contained within 443 mapping sets, complete with standardized metadata. Most notably, 50 new genome browsers that include outputs from the Wheat and Barley PanGenome projects have been created. We provide an example of an expression quantitative trait loci track on the International Wheat Genome Sequencing Consortium Chinese Spring wheat browser to demonstrate how genome browser tracks can be adapted for different data types. To help users benefit more from its data, GrainGenes created four tutorials available on YouTube. GrainGenes is executing its vision of service by continuously responding to the needs of the global small grains community by creating a centralized, long-term, interconnected data repository. Database URL:https://wheat.pw.usda.gov.

The evolving battle between yellow rust and wheat: implications for global food security

Wheat (Triticum aestivum L.) is a global commodity, and its production is a key component underpinning worldwide food security. Yellow rust, also known as stripe rust, is a wheat disease caused by the fungus Puccinia striiformis Westend f. sp. tritici (Pst), and results in yield losses in most wheat growing areas. Recently, the rapid global spread of genetically diverse sexually derived Pst races, which have now largely replaced the previous clonally propagated slowly evolving endemic populations, has resulted in further challenges for the protection of global wheat yields. However, advances in the application of genomics approaches, in both the host and pathogen, combined with classical genetic approaches, pathogen and disease monitoring, provide resources to help increase the rate of genetic gain for yellow rust resistance via wheat breeding while reducing the carbon footprint of the crop. Here we review key elements in the evolving battle between the pathogen and host, with a focus on solutions to help protect future wheat production from this globally important disease.

Wheat genetic loci conferring resistance to stripe rust in the face of genetically diverse races of the fungus Puccinia striiformis f. sp. tritici

KEY MESSAGE

Analysis of a wheat multi-founder population identified 14 yellow rust resistance QTL. For three of the four most significant QTL, haplotype analysis indicated resistance alleles were rare in European wheat. Stripe rust, or yellow rust (YR), is a major fungal disease of wheat (Triticum aestivum) caused by Puccinia striiformis Westend f. sp. tritici (Pst). Since 2011, the historically clonal European Pst races have been superseded by the rapid incursion of genetically diverse lineages, reducing the resistance of varieties previously showing durable resistance. Identification of sources of genetic resistance to such races is a high priority for wheat breeding. Here we use a wheat eight-founder multi-parent population genotyped with a 90,000 feature single nucleotide polymorphism array to genetically map YR resistance to such new Pst races. Genetic analysis of five field trials at three UK sites identified 14 quantitative trait loci (QTL) conferring resistance. Of these, four highly significant loci were consistently identified across all test environments, located on chromosomes 1A (QYr.niab-1A.1), 2A (QYr.niab-2A.1), 2B (QYr.niab-2B.1) and 2D (QYr.niab-2D.1), together explaining ~ 50% of the phenotypic variation. Analysis of these four QTL in two-way and three-way combinations showed combinations conferred greater resistance than single QTL, and genetic markers were developed that distinguished resistant and susceptible alleles. Haplotype analysis in a collection of wheat varieties found that the haplotypes associated with YR resistance at three of these four major loci were rare (≤ 7%) in European wheat, highlighting their potential utility for future targeted improvement of disease resistance. Notably, the physical interval for QTL QYr.niab-2B.1 contained five nucleotide-binding leucine-rich repeat candidate genes with integrated BED domains, of which two corresponded to the cloned resistance genes Yr7 and Yr5/YrSp.