454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction

Rust Fungi on Medicinal Plants in Guizhou Province with Descriptions of Three New Species

During the research on rust fungi in medicinal plants of Guizhou Province, China, a total of 9 rust fungal species were introduced, including 3 new species (Hamaspora rubi-alceifolii, Nyssopsora altissima, and Phragmidium cymosum), as well as 6 known species (Melampsora laricis-populina, Melampsoridium carpini, Neophysopella ampelopsidis, Nyssopsora koelrezidis, P. rosae-roxburghii, P. tormentillae). Notably, N. ampelopsidis and P. tormentillae were discovered for the first time in China, while M. laricis-populina, Me. carpini, and Ny. koelreuteriae were first documented in Guizhou Province. Morphological observation and molecular phylogenetic analyses of these species with similar taxa were compared to confirm their taxonomic identities, and taxonomic descriptions, illustrations and host species of those rust fungi on medicinal plant are provided.

A Remarkable Expansion of Oligopeptide Transporter Genes in Rust Fungi (Pucciniales) Suggests a Specialization in Nutrient Acquisition for Obligate Biotrophy

Nutrient acquisition by rust fungi during their biotrophic growth has been assigned to a few transporters expressed in haustorial infection structures. We performed a comparative genomic analysis of all transporter genes (hereafter termed transportome) classified according to the Transporter Classification Database, focusing specifically on rust fungi (order Pucciniales) versus other species in the Dikarya. We also surveyed expression of transporter genes in the poplar rust fungus for which transcriptomics data are available across the whole life cycle. Despite a significant increase in gene number, rust fungi presented a reduced transportome compared with most fungi in the Dikarya. However, a few transporter families in the subclass Porters showed significant expansions. Notably, three metal transport-related families involved in the import, export, and sequestration of metals were expanded in Pucciniales and expressed at various stages of the rust life cycle, suggesting a tight regulation of metal homeostasis. The most remarkable gene expansion in the Pucciniales was observed for the oligopeptide transporter (OPT) family, with 25 genes on average compared with seven to 14 genes in the other surveyed taxonomical ranks. A phylogenetic analysis showed several specific expansion events at the root of the order Pucciniales with subsequent expansions in rust taxonomical families. The OPT genes showed dynamic expression patterns along the rust life cycle and more particularly during infection of the poplar host tree, suggesting a possible specialization for the acquisition of nitrogen and sulfur through the transport of oligopeptides from the host during biotrophic growth.

Identification and verification of seed development related miRNAs in kernel almond by small RNA sequencing and qPCR

Many studies have investigated the role of miRNAs on the yield of various plants, but so far, no report is available on the identification and role of miRNAs in fruit and seed development of almonds. In this study, preliminary analysis by high-throughput sequencing of short RNAs of kernels from the crosses between almond cultivars 'Sefid' × 'Mamaee' (with small and large kernels, respectively) and 'Sefid' × 'P. orientalis' (with small kernels) showed that the expressions of several miRNAs such as Pdu-miR395a-3p, Pdu-miR8123-5p, Pdu-miR482f, Pdu-miR6285, and Pdu-miR396a were significantly different. These miRNAs targeted genes encoding different proteins such as NYFB-3, SPX1, PGSIP3 (GUX2), GH3.9, and BEN1. The result of RT-qPCR revealed that the expression of these genes showed significant differences between the crosses and developmental stages of the seeds, suggesting that these genes might be involved in controlling kernel size because the presence of these miRNAs had a negative effect on their target genes. Pollen source can influence kernel size by affecting hormonal signaling and metabolic pathways through related miRNAs, a phenomenon known as xenia.

Genome-Wide Analysis of Four Pathotypes of Wheat Rust Pathogen (Puccinia graminis) Reveals Structural Variations and Diversifying Selection

Diseases caused by Puccinia graminis are some of the most devastating diseases of wheat. Extensive genomic understanding of the pathogen has proven helpful not only in understanding host- pathogen interaction but also in finding appropriate control measures. In the present study, whole-genome sequencing of four diverse P. graminis pathotypes was performed to understand the genetic variation and evolution. An average of 63.5 Gb of data per pathotype with about 100× average genomic coverage was achieved with 100-base paired-end sequencing performed with Illumina Hiseq 1000. Genome structural annotations collectively predicted 9273 functional proteins including ~583 extracellular secreted proteins. Approximately 7.4% of the genes showed similarity with the PHI database which is suggestive of their significance in pathogenesis. Genome-wide analysis demonstrated pathotype 117-6 as likely distinct and descended through a different lineage. The 3-6% more SNPs in the regulatory regions and 154 genes under positive selection with their orthologs and under negative selection in the other three pathotypes further supported pathotype 117-6 to be highly diverse in nature. The genomic information generated in the present study could serve as an important source for comparative genomic studies across the genus Puccinia and lead to better rust management in wheat.

Host Adaptation and Virulence in Heteroecious Rust Fungi

Rust fungi (Pucciniales, Basidiomycota) are obligate biotrophic pathogens that cause rust diseases in plants, inflicting severe damage to agricultural crops. Pucciniales possess the most complex life cycles known in fungi. These include an alternation of generations, the development of up to five different sporulating stages, and, for many species, the requirement of infecting two unrelated host plants during different parts of their life cycle, termed heteroecism. These fungi have been extensively studied in the past century through microscopy and inoculation studies, providing precise descriptions of their infection processes, although the molecular mechanisms underlying their unique biology are poorly understood. In this review, we cover recent genomic and life cycle transcriptomic studies in several heteroecious rust species, which provide insights into the genetic tool kits associated with host adaptation and virulence, opening new avenues for unraveling their unique evolution.

The concentration of polyphenolic compounds and trace elements in the Coffea arabica leaves: Potential chemometric pattern recognition of coffee leaf rust resistance

Coffee (Coffea arabica L.) is an important commodity, involving about 500 million people from the cultivation of the coffee trees to final consumption of infusions of the ground roasted coffee beans. In contrast to a considerable amount of research performed on green coffee beans, there are relatively few studies regarding the chemical constituents of coffee leaves. Hemileia vastatrix is a parasite, specific to coffee plants and causes coffee leaf rust, which is a very destructive disease. Some coffee plants have natural resistance which is mainly linked to a gene and specific host resistance response. An increase in flavonoid production may be related to fungal disease resistance, with the levels and flavonoid types being an early physiological response to rust infection. Trace inorganic elements can be related to many roles in the defense response of higher plants and can be used as a biomarker for some diseases. To address this, coffee leaves from 16 different cultivars of Coffea arabica were harvested from Minas Gerais, Brazil (susceptible and resistant to rust) and their polyphenolic compounds were extracted using the QuEChERS technique and quantitated by HPLC-ESI-MS. The same leaves were decomposed using an acid mixture in a block digester and the content of Al, Cu, Mg, Mn, Ni, Sn and Zn was quantitated by ICP-OES. Principal component analysis (PCA) was applied and we could establish a relation between polyphenolic and trace element concentration in the leaves with resistance to rust infection. On this basis in this preliminary study we were able to separate the resistant from the susceptible cultivars. The main compounds responsible for this differentiation were the content of chlorogenic acid and magnesium in the leaves. The content of polyphenolic compounds was lower in susceptible cultivars and a diametric effect was observed between Mn and Mg concentrations. This study shows potential for the discrimination of resistant and susceptible coffee trees based on the analyses of both trace element and polyphenolic concentration.

Progress on Molecular Genetics and Manipulation of Rust Fungi

Among the thousands of rust species described, many are known for their devastating effects on their hosts, which include major agriculture crops and trees. Hence, for over a century, these basidiomycete pathogenic fungi have been researched and experimented with. However, due to their biotrophic nature, they are challenging organisms to work with and, needing their hosts for propagation, represent pathosystems that are not easily experimentally accessible. Indeed, efforts to perform genetics have been few and far apart for the rust fungi, though one study performed in the 1940s was famously instrumental in formulating the gene-for-gene hypothesis describing pathogen-host interactions. By taking full advantage of the molecular genetic tools developed in the 1980s, research on many plant pathogenic microbes thrived, yet similar work on the rusts remained very challenging though not without some successes. However, the genomics era brought real breakthrough research for the biotrophic fungi and with innovative experimentation and the use of heterologous systems, molecular genetic analyses over the last 2 decades have significantly advanced our insight into the function of many rust fungus genes and their role in the interaction with their hosts. This has allowed optimizing efforts for resistance breeding and the design and testing of various novel strategies to reduce the devastating diseases they cause.

Transcriptome Analysis of Apple Leaves Infected by the Rust Fungus Gymnosporangium yamadae at Two Sporulation Stages

Apple rust disease caused by Gymnosporangium yamadae is one of the major threats to apple orchards. In this study, dual RNA-seq analysis was conducted to simultaneously monitor gene expression profiles of G. yamadae and infected apple leaves during the formation of rust spermogonia and aecia. The molecular mechanisms underlying this compatible interaction at 10 and 30 days postinoculation (dpi) indicate a significant reaction from the host plant and comprise detoxication pathways at the earliest stage and the induction of secondary metabolism pathways at 30 dpi. Such host reactions have been previously reported in other rust pathosystems and may represent a general reaction to rust infection. G. yamadae transcript profiling indicates a conserved genetic program in spermogonia and aecia that is shared with other rust fungi, whereas secretome prediction reveals the presence of specific secreted candidate effector proteins expressed during apple infection. Unexpectedly, the survey of fungal unigenes in the transcriptome assemblies of inoculated and mock-inoculated apple leaves reveals that G. yamadae infection may modify the fungal community composition in the apple phyllosphere at 30 dpi. Collectively, our results provide novel insights into the compatible apple-G. yamadae interaction and advance the knowledge of this heteroecious demicyclic rust fungus.

Prehaustorial local resistance to coffee leaf rust in a Mexican cultivar involves expression of salicylic acid-responsive genes

Background

In Mexico, coffee leaf rust (CLR) is the main disease that affects the Arabica coffee crop. In this study, the local response of two Mexican cultivars of Coffea arabica (Oro Azteca and Garnica) in the early stages of Hemileia vastatrix infection was evaluated.

Methods

We quantified the development of fungal structures in locally-infected leaf disks from both cultivars, using qRT-PCR to measure the relative expression of two pathogenesis recognition genes (CaNDR1 and CaNBS-LRR) and three genes associated with the salicylic acid (SA)-related pathway (CaNPR1, CaPR1, and CaPR5).

Results

Resistance of the cv. Oro Azteca was significantly higher than that of the cv. Garnica, with 8.2% and 53.3% haustorial detection, respectively. In addition, the non-race specific disease resistance gene (CaNDR1), a key gene for the pathogen recognition, as well as the genes associated with SA, CaNPR1, CaPR1, and CaPR5, presented an increased expression in response to infection by H. vastatrix in cv. Oro Azteca if comparing with cv. Garnica. Our results suggest that Oro Azteca’s defense mechanisms could involve early recognition of CLR by NDR1 and the subsequent activation of the SA signaling pathway.

Transcriptional patterns of Coffea arabica L. nitrate reductase, glutamine and asparagine synthetase genes are modulated under nitrogen suppression and coffee leaf rust

This study evaluated the transcriptional profile of genes related to nitrogen (N) assimilation in coffee plants susceptible and resistant to rust fungi under N sufficiency and N suppression. For this purpose, we inoculated young coffee leaves with Hemileia vastatrix uredospores and collected them at 0, 12, 24 and 48 hours post-inoculation (HPI) to evaluate the relative expressions of genes encoding cytosolic glutamine synthetase (CaGS₁), plastid glutamine synthetase (CaGS₂), nitrate reductase (CaNR), and asparagine synthetase (CaAS). The genes exhibited distinct patterns of transcriptional modulation for the different genotypes and N nutritional regimes. The resistant genotype (I59) presented high levels of transcription in response to pathogen inoculation for CaNR and CaGS₁ genes, evaluated under N sufficiency in the initial moments of infection (12 HPI). The gene CaGS₁ also showed a peak at 48 HPI. The susceptible genotype (CV99) showed increased transcript rates of CaNR at 12 and 24 HPI in response to rust inoculation. The transcriptional patterns observed for CV99, under N suppression, were high levels for CaAS and CaGS₂ at all post-inoculation times in response to coffee leaf rust disease. In addition, CaGS₁ was up-regulated at 48 HPI for CV99. Cultivar I59 showed high transcript levels at 12 HPI for CaAS and peaks at 24 and 48 HPI for CaGS₂ in inoculated samples. Consequently, total chlorophyl concentration was influenced by N suppression and by rust infection. Regarding enzyme activities in vitro for glutamine synthetase and CaNR, there was an increase in infected coffee leaves (I59) and under N sufficiency. Moreover, CV99 was modulated in both N nutritional regimes for GS activity in response to rust. Our results indicate that N transport genes trigger a differential modulation between genotypes through the action of rust disease.

Flax rust infection transcriptomics reveals a transcriptional profile that may be indicative for rust Avr genes

Secreted effectors of fungal pathogens are essential elements for disease development. However, lack of sequence conservation among identified effectors has long been a problem for predicting effector complements in fungi. Here we have explored the expression characteristics of avirulence (Avr) genes and candidate effectors of the flax rust fungus, Melampsora lini. We performed transcriptome sequencing and real-time quantitative PCR (qPCR) on RNA extracted from ungerminated spores, germinated spores, isolated haustoria and flax seedlings inoculated with M. lini isolate CH5 during plant infection. Genes encoding two categories of M. lini proteins, namely Avr proteins and plant cell wall degrading enzymes (CWDEs), were investigated in detail. Analysis of the expression profiles of 623 genes encoding predicted secreted proteins in the M. lini transcriptome shows that the six known Avr genes (i.e. AvrM (avrM), AvrM14, AvrL2, AvrL567, AvrP123 (AvrP) and AvrP4) fall within a group of 64 similarly expressed genes that are induced in planta and show a peak of expression early in infection with a subsequent decline towards sporulation. Other genes within this group include two paralogues of AvrL2, an AvrL567 virulence allele, and a number of genes encoding putative effector proteins. By contrast, M. lini genes encoding CWDEs fall into different expression clusters with their distribution often unrelated to their catalytic activity or substrate targets. These results suggest that synthesis of M. lini Avr proteins may be regulated in a coordinated fashion and that the expression profiling-based analysis has significant predictive power for the identification of candidate Avr genes.

Comparative transcriptomics of Gymnosporangium spp. teliospores reveals a conserved genetic program at this specific stage of the rust fungal life cycle

BACKGROUND

Gymnosporangium spp. are fungal plant pathogens causing rust disease and most of them are known to infect two different host plants (heteroecious) with four spore stages (demicyclic). In the present study, we sequenced the transcriptome of G. japonicum teliospores on its host plant Juniperus chinensis and we performed comparison to the transcriptomes of G. yamadae and G. asiaticum at the same life stage, that happens in the same host but on different organs.

RESULTS

Functional annotation for the three Gymnosporangium species showed the expression of a conserved genetic program with the top abundant cellular categories corresponding to energy, translation and signal transduction processes, indicating that this life stage is particularly active. Moreover, the survey of predicted secretomes in the three Gymnosporangium transcriptomes revealed shared and specific genes encoding carbohydrate active enzymes and secreted proteins of unknown function that could represent candidate pathogenesis effectors. A transcript encoding a hemicellulase of the glycoside hydrolase 26 family, previously identified in other rust fungi, was particularly highly expressed suggesting a general role in rust fungi. The comparison between the transcriptomes of the three Gymnosporangium spp. and selected Pucciniales species in different taxonomical families allowed to identify lineage-specific protein families that may relate to the biology of teliospores in rust fungi. Among clustered gene families, 205, 200 and 152 proteins were specifically identified in G. japonicum, G. yamadae and G. asiaticum, respectively, including candidate effectors expressed in teliospores.

CONCLUSIONS

This comprehensive comparative transcriptomics study of three Gymnosporangium spp. identified gene functions and metabolic pathways particularly expressed in teliospores, a stage of the life cycle that is mostly overlooked in rust fungi. Secreted protein encoding transcripts expressed in teliospores may reveal new candidate effectors related to pathogenesis. Although this spore stage is not involved in host plant infection but in the production of basidiospores infecting plants in the Amygdaloideae, we speculate that candidate effectors may be expressed as early as the teliospore stage for preparing further infection by basidiospores.

Advances in understanding obligate biotrophy in rust fungi

Contents Summary 1190 I. Introduction 1190 II. Rust fungi: a diverse and serious threat to agriculture 1191 III. The different facets of rust life cycles and unresolved questions about their evolution 1191 IV. The biology of rust infection 1192 V. Rusts in the genomics era: the ever-expanding list of candidate effector genes 1195 VI. Functional characterization of rust effectors 1197 VII. Putting rusts to sleep: Pucciniales research outlooks 1201 Acknowledgements 1202 References 1202 SUMMARY: Rust fungi (Pucciniales) are the largest group of plant pathogens and represent one of the most devastating threats to agricultural crops worldwide. Despite the economic importance of these highly specialized pathogens, many aspects of their biology remain obscure, largely because rust fungi are obligate biotrophs. The rise of genomics and advances in high-throughput sequencing technology have presented new options for identifying candidate effector genes involved in pathogenicity mechanisms of rust fungi. Transcriptome analysis and integrated bioinformatics tools have led to the identification of key genetic determinants of host susceptibility to infection by rusts. Thousands of genes encoding secreted proteins highly expressed during host infection have been reported for different rust species, which represents significant potential towards understanding rust effector function. Recent high-throughput in planta expression screen approaches (effectoromics) have pushed the field ahead even further towards predicting high-priority effectors and identifying avirulence genes. These new insights into rust effector biology promise to inform future research and spur the development of effective and sustainable strategies for managing rust diseases.

Genome sequencing and transcript analysis of Hemileia vastatrix reveal expression dynamics of candidate effectors dependent on host compatibility

Coffee leaf rust caused by the fungus Hemileia vastatrix is one of the most important leaf diseases of coffee plantations worldwide. Current knowledge of the H. vastatrix genome is limited and only a small fraction of the total fungal secretome has been identified. In order to obtain a more comprehensive understanding of its secretome, we aimed to sequence and assemble the entire H. vastatrix genome using two next-generation sequencing platforms and a hybrid assembly strategy. This resulted in a 547 Mb genome of H. vastatrix race XXXIII (Hv33), with 13,364 predicted genes that encode 13,034 putative proteins with transcriptomic support. Based on this proteome, 615 proteins contain putative secretion peptides, and lack transmembrane domains. From this putative secretome, 111 proteins were identified as candidate effectors (EHv33) unique to H. vastatrix, and a subset consisting of 17 EHv33 genes was selected for a temporal gene expression analysis during infection. Five genes were significantly induced early during an incompatible interaction, indicating their potential role as pre-haustorial effectors possibly recognized by the resistant coffee genotype. Another nine genes were significantly induced after haustorium formation in the compatible interaction. Overall, we suggest that this fungus is able to selectively mount its survival strategy with effectors that depend on the host genotype involved in the infection process.

Alterations in plant sugar metabolism: signatory of pathogen attack

This review summarizes the current understanding, future challenges and ongoing quest on sugar metabolic alterations that influence the outcome of plant-pathogen interactions. Intricate cellular and molecular events occur during plant-pathogen interactions. They cause major metabolic perturbations in the host and alterations in sugar metabolism play a pivotal role in governing the outcome of various kinds of plant-pathogen interactions. Sugar metabolizing enzymes and transporters of both host and pathogen origin get differentially regulated during the interactions. Both plant and pathogen compete for utilizing the host sugar metabolic machinery and in turn promote resistant or susceptible responses. However, the kind of sugar metabolism alteration that is beneficial for the host or pathogen is yet to be properly understood. Recently developed tools and methodologies are facilitating research to understand the intricate dynamics of sugar metabolism during the interactions. The present review elaborates current understanding, future challenges and ongoing quest on sugar metabolism, mobilization and regulation during various plant-pathogen interactions.

An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: insights on nitrogen starvation responses

Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.

High throughput transcriptome analysis of coffee reveals prehaustorial resistance in response to Hemileia vastatrix infection

We provide a transcriptional profile of coffee rust interaction and identified putative up regulated resistant genes Coffee rust disease, caused by the fungus Hemileia vastatrix, is one of the major diseases in coffee throughout the world. The use of resistant cultivars is considered to be the most effective control strategy for this disease. To identify candidate genes related to different mechanism defense in coffee, we present a time-course comparative gene expression profile of Caturra (susceptible) and Híbrido de Timor (HdT, resistant) in response to H. vastatrix race XXXIII infection. The main objectives were to obtain a global overview of transcriptome in both interaction, compatible and incompatible, and, specially, analyze up-regulated HdT specific genes with inducible resistant and defense signaling pathways. Using both Coffea canephora as a reference genome and de novo assembly, we obtained 43,159 transcripts. At early infection events (12 and 24 h after infection), HdT responded to the attack of H. vastatrix with a larger number of up-regulated genes than Caturra, which was related to prehaustorial resistance. The genes found in HdT at early hours were involved in receptor-like kinases, response ion fluxes, production of reactive oxygen species, protein phosphorylation, ethylene biosynthesis and callose deposition. We selected 13 up-regulated HdT-exclusive genes to validate by real-time qPCR, which most of them confirmed their higher expression in HdT than in Caturra at early stage of infection. These genes have the potential to assist the development of new coffee rust control strategies. Collectively, our results provide understanding of expression profiles in coffee-H. vastatrix interaction over a time course in susceptible and resistant coffee plants.

The coffee leaf rust pathogen Hemileia vastatrix: one and a half centuries around the tropics

TAXONOMY AND HISTORY

Hemileia vastatrix Berk. and Broome (Basidiomycota, Pucciniales) was described in 1869 in eastern Africa and Ceylon as the agent of coffee leaf rust and has spread to all coffee cultivation areas worldwide. Major disease outbreaks in Asia, Africa and America caused and continue to cause severe yield losses, making this the most important disease of Arabica coffee, a cash crop for many tropical and sub-tropical countries.

LIFE CYCLE AND DISEASE SYMPTOMS

Hemileia vastatrix is a hemicyclic fungus with the urediniosporic life cycle as its most important (if not only) source of inoculum. Chlorotic spots are the first macroscopic symptoms, preceding the differentiation of suprastomatal, bouquet-shaped, orange-coloured uredinia. The disease can cause yield losses of up to 35% and have a polyetic epidemiological impact on subsequent years.

DISEASE CONTROL

Although the use of fungicides is one of the preferred immediate control measures, the use of resistant cultivars is considered to be the most effective and durable disease control strategy. The discovery of 'Híbrido de Timor' provided sources of resistance that have been used in several breeding programmes and that have been proven to be effective and durable, as some have been in use for more than 30 years.

GENETIC DIVERSITY AND MOLECULAR PATHOGENICITY

Although exhibiting limited genetic polymorphism, the very large genome of H. vastatrix (c. 797 Mbp) conceals great pathological diversity, with more than 50 physiological races. Gene expression studies have revealed a very precocious activation of signalling pathways and production of putative effectors, suggesting that the plant-fungus dialogue starts as early as at the germ tube stage, and have provided clues for the identification of avr genes.

Comparative transcriptome analysis and identification of candidate effectors in two related rust species (Gymnosporangium yamadae and Gymnosporangium asiaticum)

BACKGROUND

Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection.

RESULTS

The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya.

CONCLUSION

This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus.

The Hemileia vastatrix effector HvEC-016 suppresses bacterial blight symptoms in coffee genotypes with the SH 1 rust resistance gene

A number of genes that confer resistance to coffee leaf rust (S_H 1-S_H 9) have been identified within the genus Coffea, but despite many years of research on this pathosystem, the complementary avirulence genes of Hemileia vastatrix have not been reported. After identification of H. vastatrix effector candidate genes (HvECs) expressed at different stages of its lifecycle, we established an assay to characterize HvEC proteins by delivering them into coffee cells via the type-three secretion system (T3SS) of Pseudomonas syringae pv. garcae (Psgc). Employing a calmodulin-dependent adenylate cyclase assay, we demonstrate that Psgc recognizes a heterologous P. syringae T3SS secretion signal which enables us to translocate HvECs into the cytoplasm of coffee cells. Using this Psgc-adapted effector detector vector (EDV) system, we found that HvEC-016 suppresses the growth of Psgc on coffee genotypes with the S_H 1 resistance gene. Suppression of bacterial blight symptoms in S_H 1 plants was associated with reduced bacterial multiplication. By contrast, HvEC-016 enhanced bacterial multiplication in S_H 1-lacking plants. Our findings suggest that HvEC-016 may be recognized by the plant immune system in a S_H 1-dependent manner. Thus, our experimental approach is an effective tool for the characterization of effector/avirulence proteins of this important pathogen.

Production of small cysteine-rich effector proteins in Escherichia coli for structural and functional studies

Although the lifestyles and infection strategies of plant pathogens are diverse, a prevailing feature is the use of an arsenal of secreted proteins, known as effectors, which aid in microbial infection. In the case of eukaryotic filamentous pathogens, such as fungi and oomycetes, effector proteins are typically dissimilar, at the protein sequence level, to known protein families and functional domains. Consequently, we currently have a limited understanding of how fungal and oomycete effectors promote disease. Protein biochemistry and structural biology are two methods that can contribute greatly to the understanding of protein function. Both techniques are dependent on obtaining proteins that are pure and functional, and generally require the use of heterologous recombinant protein expression systems. Here, we present a general scheme and methodology for the production and characterization of small cysteine-rich (SCR) effectors utilizing Escherichia coli expression systems. Using this approach, we successfully produced cysteine-rich effectors derived from the biotrophic fungal pathogen Melampsora lini and the necrotrophic fungal pathogen Parastagonospora nodorum. Access to functional recombinant proteins facilitated crystallization and functional experiments. These results are discussed in the context of a general workflow that may serve as a template for others interested in understanding the function of SCR effector(s) from their plant pathogen(s) of interest.

Dual RNA-seq reveals Meloidogyne graminicola transcriptome and candidate effectors during the interaction with rice plants

Root-knot nematodes secrete proteinaceous effectors into plant tissues to facilitate infection by suppressing host defences and reprogramming the host metabolism to their benefit. Meloidogyne graminicola is a major pest of rice (Oryza sativa) in Asia and Latin America, causing important crop losses. The goal of this study was to identify M. graminicola pathogenicity genes expressed during the plant-nematode interaction. Using the dual RNA-sequencing (RNA-seq) strategy, we generated transcriptomic data of M. graminicola samples covering the pre-parasitic J2 stage and five parasitic stages in rice plants, from the parasitic J2 to the adult female. In the absence of a reference genome, a de novo M. graminicola transcriptome of 66 396 contigs was obtained from those reads that were not mapped on the rice genome. Gene expression profiling across the M. graminicola life cycle revealed key genes involved in nematode development and provided insights into the genes putatively associated with parasitism. The development of a 'secreted protein prediction' pipeline revealed a typical set of proteins secreted by nematodes, as well as a large number of cysteine-rich proteins and putative nuclear proteins. Combined with expression data, this pipeline enabled the identification of 15 putative effector genes, including two homologues of well-characterized effectors from cyst nematodes (CLE-like and VAP1) and a metallothionein. The localization of gene expression was assessed by in situ hybridization for a subset of candidates. All of these data represent important molecular resources for the elucidation of M. graminicola biology and for the selection of potential targets for the development of novel control strategies for this nematode species.

Identification of candidate genes for drought tolerance in coffee by high-throughput sequencing in the shoot apex of different Coffea arabica cultivars

BACKGROUND

Drought is a widespread limiting factor in coffee plants. It affects plant development, fruit production, bean development and consequently beverage quality. Genetic diversity for drought tolerance exists within the coffee genus. However, the molecular mechanisms underlying the adaptation of coffee plants to drought are largely unknown. In this study, we compared the molecular responses to drought in two commercial cultivars (IAPAR59, drought-tolerant and Rubi, drought-susceptible) of Coffea arabica grown in the field under control (irrigation) and drought conditions using the pyrosequencing of RNA extracted from shoot apices and analysing the expression of 38 candidate genes.

RESULTS

Pyrosequencing from shoot apices generated a total of 34.7 Mbp and 535,544 reads enabling the identification of 43,087 clusters (41,512 contigs and 1,575 singletons). These data included 17,719 clusters (16,238 contigs and 1,575 singletons) exclusively from 454 sequencing reads, along with 25,368 hybrid clusters assembled with 454 sequences. The comparison of DNA libraries identified new candidate genes (n = 20) presenting differential expression between IAPAR59 and Rubi and/or drought conditions. Their expression was monitored in plagiotropic buds, together with those of other (n = 18) candidates genes. Under drought conditions, up-regulated expression was observed in IAPAR59 but not in Rubi for CaSTK1 (protein kinase), CaSAMT1 (SAM-dependent methyltransferase), CaSLP1 (plant development) and CaMAS1 (ABA biosynthesis). Interestingly, the expression of lipid-transfer protein (nsLTP) genes was also highly up-regulated under drought conditions in IAPAR59. This may have been related to the thicker cuticle observed on the abaxial leaf surface in IAPAR59 compared to Rubi.

CONCLUSIONS

The full transcriptome assembly of C. arabica, followed by functional annotation, enabled us to identify differentially expressed genes related to drought conditions. Using these data, candidate genes were selected and their differential expression profiles were confirmed by qPCR experiments in plagiotropic buds of IAPAR59 and Rubi under drought conditions. As regards the genes up-regulated under drought conditions, specifically in the drought-tolerant IAPAR59, several corresponded to orphan genes but also to genes coding proteins involved in signal transduction pathways, as well as ABA and lipid metabolism, for example. The identification of these genes should help advance our understanding of the genetic determinism of drought tolerance in coffee.

Identification of Phakopsora pachyrhizi Candidate Effectors with Virulence Activity in a Distantly Related Pathosystem

Phakopsora pachyrhizi is the causal agent of Asian Soybean Rust, a disease that causes enormous economic losses, most markedly in South America. P. pachyrhizi is a biotrophic pathogen that utilizes specialized feeding structures called haustoria to colonize its hosts. In rusts and other filamentous plant pathogens, haustoria have been shown to secrete effector proteins into their hosts to permit successful completion of their life cycle. We have constructed a cDNA library from P. pachyrhizi haustoria using paramagnetic bead-based methodology and have identified 35 P. pachyrhizi candidate effector (CE) genes from this library which are described here. In addition, we quantified the transcript expression pattern of six of these genes and show that two of these CEs are able to greatly increase the susceptibility of Nicotiana benthamiana to Phytophthora infestans. This strongly suggests that these genes play an important role in P. pachyrhizi virulence on its hosts.

Transcriptional dynamics of Phytophthora infestans during sequential stages of hemibiotrophic infection of tomato

Hemibiotrophic plant pathogens, such as the oomycete Phytophthora infestans, employ a biphasic infection strategy, initially behaving as biotrophs, where minimal symptoms are exhibited by the plant, and subsequently as necrotrophs, feeding on dead plant tissue. The regulation of this transition and the breadth of molecular mechanisms that modulate plant defences are not well understood, although effector proteins secreted by the pathogen are thought to play a key role. We examined the transcriptional dynamics of P. infestans in a compatible interaction with its host tomato (Solanum lycopersicum) at three infection stages: biotrophy; the transition from biotrophy to necrotrophy; and necrotrophy. The expression data suggest a tight temporal regulation of many pathways associated with the suppression of plant defence mechanisms and pathogenicity, including the induction of putative cytoplasmic and apoplastic effectors. Twelve of these were experimentally evaluated to determine their ability to suppress necrosis caused by the P. infestans necrosis-inducing protein PiNPP1.1 in Nicotiana benthamiana. Four effectors suppressed necrosis, suggesting that they might prolong the biotrophic phase. This study suggests that a complex regulation of effector expression modulates the outcome of the interaction.

Effector-Mining in the Poplar Rust Fungus Melampsora larici-populina Secretome

The poplar leaf rust fungus, Melampsora larici-populina has been established as a tree-microbe interaction model. Understanding the molecular mechanisms controlling infection by pathogens appears essential for durable management of tree plantations. In biotrophic plant-parasites, effectors are known to condition host cell colonization. Thus, investigation of candidate secreted effector proteins (CSEPs) is a major goal in the poplar-poplar rust interaction. Unlike oomycetes, fungal effectors do not share conserved motifs and candidate prediction relies on a set of a priori criteria established from reported bona fide effectors. Secretome prediction, genome-wide analysis of gene families and transcriptomics of M. larici-populina have led to catalogs of more than a thousand secreted proteins. Automatized effector-mining pipelines hold great promise for rapid and systematic identification and prioritization of CSEPs for functional characterization. In this review, we report on and discuss the current status of the poplar rust fungus secretome and prediction of candidate effectors from this species.

Genes related to antioxidant metabolism are involved in Methylobacterium mesophilicum-soybean interaction

The genus Methylobacterium is composed of pink-pigmented methylotrophic bacterial species that are widespread in natural environments, such as soils, stream water and plants. When in association with plants, this genus colonizes the host plant epiphytically and/or endophytically. This association is known to promote plant growth, induce plant systemic resistance and inhibit plant infection by phytopathogens. In the present study, we focused on evaluating the colonization of soybean seedling-roots by Methylobacterium mesophilicum strain SR1.6/6. We focused on the identification of the key genes involved in the initial step of soybean colonization by methylotrophic bacteria, which includes the plant exudate recognition and adaptation by planktonic bacteria. Visualization by scanning electron microscopy revealed that M. mesophilicum SR1.6/6 colonizes soybean roots surface effectively at 48 h after inoculation, suggesting a mechanism for root recognition and adaptation before this period. The colonization proceeds by the development of a mature biofilm on roots at 96 h after inoculation. Transcriptomic analysis of the planktonic bacteria (with plant) revealed the expression of several genes involved in membrane transport, thus confirming an initial metabolic activation of bacterial responses when in the presence of plant root exudates. Moreover, antioxidant genes were mostly expressed during the interaction with the plant exudates. Further evaluation of stress- and methylotrophic-related genes expression by qPCR showed that glutathione peroxidase and glutathione synthetase genes were up-regulated during the Methylobacterium-soybean interaction. These findings support that glutathione (GSH) is potentially a key molecule involved in cellular detoxification during plant root colonization. In addition to methylotrophic metabolism, antioxidant genes, mainly glutathione-related genes, play a key role during soybean exudate recognition and adaptation, the first step in bacterial colonization.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors

Background

The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq.

Results

cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses.

Conclusion

This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1861-1) contains supplementary material, which is available to authorized users.

Proteomic analysis of apoplastic fluid of Coffea arabica leaves highlights novel biomarkers for resistance against Hemileia vastatrix

A proteomic analysis of the apoplastic fluid (APF) of coffee leaves was conducted to investigate the cellular processes associated with incompatible (resistant) and compatible (susceptible) Coffea arabica-Hemileia vastatrix interactions, during the 24-96 hai period. The APF proteins were extracted by leaf vacuum infiltration and protein profiles were obtained by 2-DE. The comparative analysis of the gels revealed 210 polypeptide spots whose volume changed in abundance between samples (control, resistant and susceptible) during the 24-96 hai period. The proteins identified were involved mainly in protein degradation, cell wall metabolism and stress/defense responses, most of them being hydrolases (around 70%), particularly sugar hydrolases and peptidases/proteases. The changes in the APF proteome along the infection process revealed two distinct phases of defense responses, an initial/basal one (24-48 hai) and a late/specific one (72-96 hai). Compared to susceptibility, resistance was associated with a higher number of proteins, which was more evident in the late/specific phase. Proteins involved in the resistance response were mainly, glycohydrolases of the cell wall, serine proteases and pathogen related-like proteins (PR-proteins), suggesting that some of these proteins could be putative candidates for resistant markers of coffee to H. vastatrix. Antibodies were produced against chitinase, pectin methylesterase, serine carboxypeptidase, reticuline oxidase and subtilase and by an immunodetection assay it was observed an increase of these proteins in the resistant sample. With this methodology we have identified proteins that are candidate markers of resistance and that will be useful in coffee breeding programs to assist in the selection of cultivars with resistance to H. vastatrix.

Host-microbe and microbe-microbe interactions in the evolution of obligate plant parasitism

Research on obligate biotrophic plant parasites, which reproduce only on living hosts, has revealed a broad diversity of filamentous microbes that have independently acquired complex morphological structures, such as haustoria. Genome studies have also demonstrated a concerted loss of genes for metabolism and lytic enzymes, and gain of diversity of genes coding for effectors involved in host defense suppression. So far, these traits converge in all known obligate biotrophic parasites, but unexpected genome plasticity remains. This plasticity is manifested as transposable element (TE)-driven increases in genome size, observed to be associated with the diversification of virulence genes under selection pressure. Genome expansion could result from the governing of the pathogen response to ecological selection pressures, such as host or nutrient availability, or to microbial interactions, such as competition, hyperparasitism and beneficial cooperations. Expansion is balanced by alternating sexual and asexual cycles, as well as selfing and outcrossing, which operate to control transposon activity in populations. In turn, the prevalence of these balancing mechanisms seems to be correlated with external biotic factors, suggesting a complex, interconnected evolutionary network in host-pathogen-microbe interactions. Therefore, the next phase of obligate biotrophic pathogen research will need to uncover how this network, including multitrophic interactions, shapes the evolution and diversity of pathogens.

Analysis of gene expression changes in peach leaves in response to Plum pox virus infection using RNA-Seq

Differences in gene expression were studied after Plum pox virus (PPV, sharka disease) infection in peach GF305 leaves with and without sharka symptoms using RNA-Seq. For each sample, more than 80% of 100-nucleotide paired-end (PE) Illumina reads were aligned on the peach reference genome. In the symptomatic sample, a significant proportion of reads were mapped to PPV reference genomes (1.04% compared with 0.00002% in non-symptomatic leaves), allowing for the ultra-deep assembly of the complete genome of the PPV isolate used (9775 nucleotides, missing only 11 nucleotides at the 5' genome end). In addition, significant alternative splicing events were detected in 359 genes and 12 990 single nucleotide polymorphisms (SNPs) were identified, 425 of which could be annotated. Gene ontology annotation revealed that the high-ranking mRNA target genes associated with the expression of sharka symptoms are mainly related to the response to biotic stimuli, to lipid and carbohydrate metabolism and to the negative regulation of catalytic activity. A greater number of differentially expressed genes were observed in the early asymptomatic phase of PPV infection in comparison with the symptomatic phase. These early infection events were associated with the induction of genes related to pathogen resistance, such as jasmonic acid, chitinases, cytokinin glucosyl transferases and Lys-M proteins. Once the virus had accumulated, the overexpression of Dicer protein 2a genes suggested a gene silencing plant response that was suppressed by the virus HCPro and P1 proteins. These results illustrate the dynamic nature of the peach-PPV interaction at the transcriptome level and confirm that sharka symptom expression is a complex process that can be understood on the basis of changes in plant gene expression.

Allelic diversity in the transcriptomes of contrasting rust-infected genotypes of Lathyrus sativus, a lasting resource for smart breeding

BACKGROUND

Grass pea (Lathyrus sativus L.) is a valuable resource for potentially durable partial resistance to rust. To gain insight into the resistance mechanism and identify potential resistance genes, we generated the first comprehensive transcriptome assemblies from control and Uromyces pisi inoculated leafs of a susceptible and a partially rust-resistant grass pea genotype by RNA-seq.

RESULTS

134,914 contigs, shared by both libraries, were used to analyse their differential expression in response to rust infection. Functional annotation grouped 60.4% of the contigs present in plant databases (37.8% of total) to 33 main functional categories, being "protein", "RNA", "signalling", "transport" and "stress" the most represented. Transcription profiles revealed considerable differences in regulation of major phytohormone signalling pathways: whereas Salicylic and Abscisic Acid pathways were up-regulated in the resistant genotype, Jasmonate and Ethylene pathways were down-regulated in the susceptible one. As potential Resistance-genes we identified a mildew resistance locus O (MLO)-like gene, and MLO-related transcripts. Also, several pathogenesis-related genes were up-regulated in the resistant and exclusively down regulated in the susceptible genotype. Pathogen effectors identified in both inoculated libraries, as e.g. the rust Rtp1 transcript, may be responsible for the down-regulation of defence-related transcripts. The two genotypes contained 4,892 polymorphic contigs with SNPs unevenly distributed between different functional categories. Protein degradation (29.7%) and signalling receptor kinases (8.2%) were the most diverged, illustrating evolutionary adaptation of grass pea to the host/pathogens arms race.

CONCLUSIONS

The vast array of novel, resistance-related genomic information we present here provides a highly valuable resource for future smart breeding approaches in this hitherto under-researched, valuable legume crop.

High-resolution transcript profiling of the atypical biotrophic interaction between Theobroma cacao and the fungal pathogen Moniliophthora perniciosa

Witches' broom disease (WBD), caused by the hemibiotrophic fungus Moniliophthora perniciosa, is one of the most devastating diseases of Theobroma cacao, the chocolate tree. In contrast to other hemibiotrophic interactions, the WBD biotrophic stage lasts for months and is responsible for the most distinctive symptoms of the disease, which comprise drastic morphological changes in the infected shoots. Here, we used the dual RNA-seq approach to simultaneously assess the transcriptomes of cacao and M. perniciosa during their peculiar biotrophic interaction. Infection with M. perniciosa triggers massive metabolic reprogramming in the diseased tissues. Although apparently vigorous, the infected shoots are energetically expensive structures characterized by the induction of ineffective defense responses and by a clear carbon deprivation signature. Remarkably, the infection culminates in the establishment of a senescence process in the host, which signals the end of the WBD biotrophic stage. We analyzed the pathogen's transcriptome in unprecedented detail and thereby characterized the fungal nutritional and infection strategies during WBD and identified putative virulence effectors. Interestingly, M. perniciosa biotrophic mycelia develop as long-term parasites that orchestrate changes in plant metabolism to increase the availability of soluble nutrients before plant death. Collectively, our results provide unique insight into an intriguing tropical disease and advance our understanding of the development of (hemi)biotrophic plant-pathogen interactions.

Gene expression in mycorrhizal orchid protocorms suggests a friendly plant-fungus relationship

Orchids fully depend on symbiotic interactions with specific soil fungi for seed germination and early development. Germinated seeds give rise to a protocorm, a heterotrophic organ that acquires nutrients, including organic carbon, from the mycorrhizal partner. It has long been debated if this interaction is mutualistic or antagonistic. To investigate the molecular bases of the orchid response to mycorrhizal invasion, we developed a symbiotic in vitro system between Serapias vomeracea, a Mediterranean green meadow orchid, and the rhizoctonia-like fungus Tulasnella calospora. 454 pyrosequencing was used to generate an inventory of plant and fungal genes expressed in mycorrhizal protocorms, and plant genes could be reliably identified with a customized bioinformatic pipeline. A small panel of plant genes was selected and expression was assessed by real-time quantitative PCR in mycorrhizal and non-mycorrhizal protocorm tissues. Among these genes were some markers of mutualistic (e.g. nodulins) as well as antagonistic (e.g. pathogenesis-related and wound/stress-induced) genes. None of the pathogenesis or wound/stress-related genes were significantly up-regulated in mycorrhizal tissues, suggesting that fungal colonization does not trigger strong plant defence responses. In addition, the highest expression fold change in mycorrhizal tissues was found for a nodulin-like gene similar to the plastocyanin domain-containing ENOD55. Another nodulin-like gene significantly more expressed in the symbiotic tissues of mycorrhizal protocorms was similar to a sugar transporter of the SWEET family. Two genes coding for mannose-binding lectins were significantly up-regulated in the presence of the mycorrhizal fungus, but their role in the symbiosis is unclear.

Functional annotation of Cotesia congregata bracovirus: identification of viral genes expressed in parasitized host immune tissues

Bracoviruses (BVs) from the Polydnaviridae family are symbiotic viruses used as biological weapons by parasitoid wasps to manipulate lepidopteran host physiology and induce parasitism success. BV particles are produced by wasp ovaries and injected along with the eggs into the caterpillar host body, where viral gene expression is necessary for wasp development. Recent sequencing of the proviral genome of Cotesia congregata BV (CcBV) identified 222 predicted virulence genes present on 35 proviral segments integrated into the wasp genome. To date, the expressions of only a few selected candidate virulence genes have been studied in the caterpillar host, and we lacked a global vision of viral gene expression. In this study, a large-scale transcriptomic analysis by 454 sequencing of two immune tissues (fat body and hemocytes) of parasitized Manduca sexta caterpillar hosts allowed the detection of expression of 88 CcBV genes expressed 24 h after the onset of parasitism. We linked the expression profiles of these genes to several factors, showing that different regulatory mechanisms control viral gene expression in the host. These factors include the presence of signal peptides in encoded proteins, diversification of promoter regions, and, more surprisingly, gene position on the proviral genome. Indeed, most genes for which expression could be detected are localized in particular proviral regions globally producing higher numbers of circles. Moreover, this polydnavirus (PDV) transcriptomic analysis also reveals that a majority of CcBV genes possess at least one intron and an arthropod transcription start site, consistent with an insect origin of these virulence genes.

IMPORTANCE

Bracoviruses (BVs) are symbiotic polydnaviruses used by parasitoid wasps to manipulate lepidopteran host physiology, ensuring wasp offspring survival. To date, the expressions of only a few selected candidate BV virulence genes have been studied in caterpillar hosts. We performed a large-scale analysis of BV gene expression in two immune tissues of Manduca sexta caterpillars parasitized by Cotesia congregata wasps. Genes for which expression could be detected corresponded to genes localized in particular regions of the viral genome globally producing higher numbers of circles. Our study thus brings an original global vision of viral gene expression and paves the way to the determination of the regulatory mechanisms enabling the expression of BV genes in targeted organisms, such as major insect pests. In addition, we identify sequence features suggesting that most BV virulence genes were acquired from insect genomes.

Analysis of plant microbe interactions in the era of next generation sequencing technologies

Next generation sequencing (NGS) technologies have impressively accelerated research in biological science during the last years by enabling the production of large volumes of sequence data to a drastically lower price per base, compared to traditional sequencing methods. The recent and ongoing developments in the field allow addressing research questions in plant-microbe biology that were not conceivable just a few years ago. The present review provides an overview of NGS technologies and their usefulness for the analysis of microorganisms that live in association with plants. Possible limitations of the different sequencing systems, in particular sources of errors and bias, are critically discussed and methods are disclosed that help to overcome these shortcomings. A focus will be on the application of NGS methods in metagenomic studies, including the analysis of microbial communities by amplicon sequencing, which can be considered as a targeted metagenomic approach. Different applications of NGS technologies are exemplified by selected research articles that address the biology of the plant associated microbiota to demonstrate the worth of the new methods.

The haustorial transcriptomes of Uromyces appendiculatus and Phakopsora pachyrhizi and their candidate effector families

Haustoria of biotrophic rust fungi are responsible for the uptake of nutrients from their hosts and for the production of secreted proteins, known as effectors, which modulate the host immune system. The identification of the transcriptome of haustoria and an understanding of the functions of expressed genes therefore hold essential keys for the elucidation of fungus-plant interactions and the development of novel fungal control strategies. Here, we purified haustoria from infected leaves and used 454 sequencing to examine the haustorial transcriptomes of Phakopsora pachyrhizi and Uromyces appendiculatus, the causal agents of soybean rust and common bean rust, respectively. These pathogens cause extensive yield losses in their respective legume crop hosts. A series of analyses were used to annotate expressed sequences, including transposable elements and viruses, to predict secreted proteins from the assembled sequences and to identify families of candidate effectors. This work provides a foundation for the comparative analysis of haustorial gene expression with further insights into physiology and effector evolution.

The genome sequence and effector complement of the flax rust pathogen Melampsora lini

Rust fungi cause serious yield reductions on crops, including wheat, barley, soybean, coffee, and represent real threats to global food security. Of these fungi, the flax rust pathogen Melampsora lini has been developed most extensively over the past 80 years as a model to understand the molecular mechanisms that underpin pathogenesis. During infection, M. lini secretes virulence effectors to promote disease. The number of these effectors, their function and their degree of conservation across rust fungal species is unknown. To assess this, we sequenced and assembled de novo the genome of M. lini isolate CH5 into 21,130 scaffolds spanning 189 Mbp (scaffold N50 of 31 kbp). Global analysis of the DNA sequence revealed that repetitive elements, primarily retrotransposons, make up at least 45% of the genome. Using ab initio predictions, transcriptome data and homology searches, we identified 16,271 putative protein-coding genes. An analysis pipeline was then implemented to predict the effector complement of M. lini and compare it to that of the poplar rust, wheat stem rust and wheat stripe rust pathogens to identify conserved and species-specific effector candidates. Previous knowledge of four cloned M. lini avirulence effector proteins and two basidiomycete effectors was used to optimize parameters of the effector prediction pipeline. Markov clustering based on sequence similarity was performed to group effector candidates from all four rust pathogens. Clusters containing at least one member from M. lini were further analyzed and prioritized based on features including expression in isolated haustoria and infected leaf tissue and conservation across rust species. Herein, we describe 200 of 940 clusters that ranked highest on our priority list, representing 725 flax rust candidate effectors. Our findings on this important model rust species provide insight into how effectors of rust fungi are conserved across species and how they may act to promote infection on their hosts.

Overview of the functional virulent genome of the coffee leaf rust pathogen Hemileia vastatrix with an emphasis on early stages of infection

Hemileia vastatrix is the causal agent of coffee leaf rust, the most important disease of coffee Arabica. In this work, a 454-pyrosequencing transcriptome analysis of H. vastatrix germinating urediniospores (gU) and appressoria (Ap) was performed and compared to previously published in planta haustoria-rich (H) data. A total of 9234 transcripts were identified and annotated. Ca. 50% of these transcripts showed no significant homology to international databases. Only 784 sequences were shared by the three conditions, and 75% were exclusive of either gU (2146), Ap (1479) or H (3270). Relative transcript abundance and RT-qPCR analyses for a selection of genes indicated a particularly active metabolism, translational activity and production of new structures in the appressoria and intense signaling, transport, secretory activity and cellular multiplication in the germinating urediniospores, suggesting the onset of a plant-fungus dialogue as early as at the germ tube stage. Gene expression related to the production of carbohydrate-active enzymes and accumulation of glycerol in germinating urediniospores and appressoria suggests that combined lytic and physical mechanisms are involved in appressoria-mediated penetration. Besides contributing to the characterization of molecular processes leading to appressoria-mediated infection by rust fungi, these results point toward the identification of new H. vastatrix candidate virulence factors, with 516 genes predicted to encode secreted proteins.

Using transcription of six Puccinia triticina races to identify the effective secretome during infection of wheat

Wheat leaf rust, caused by the basidiomycete Puccinia triticina, can cause yield losses of up to 20% in wheat producing regions. During infection, the fungus forms haustoria that secrete proteins into the plant cell and effect changes in plant transcription, metabolism, and defense. It is hypothesized that new races emerge as a result of overcoming plant resistance via changes in the secreted effector proteins. To understand gene expression during infection and find genetic differences associated with races, RNA from wheat leaves infected with six different rust races, at 6 days post inoculation, was sequenced using Illumina. As P. triticina is an obligate biotroph, RNA from both the host and fungi were present and separated by alignment to the P. triticina genome and a wheat EST reference. A total of 222,571 rust contigs were assembled from 165 million reads. An examination of the resulting contigs revealed 532 predicted secreted proteins among the transcripts. Of these, 456 were found in all races. Fifteen genes were found with amino acid changes, corresponding to putative avirulence effectors potentially recognized by 11 different leaf rust resistance (Lr) genes. Twelve of the potential avirulence effectors have no homology to known genes. One gene had significant similarity to cerato-platanin, a known fungal elicitor, and another showed similarity to fungal tyrosinase, an enzyme involved in melanin synthesis. Temporal expression profiles were developed for these genes by qRT-PCR and show that the genes expression patterns were consistent between races from infection initiation to just prior to spore eruption.

Annotation of a hybrid partial genome of the coffee rust (Hemileia vastatrix) contributes to the gene repertoire catalog of the Pucciniales

Coffee leaf rust caused by the fungus Hemileia vastatrix is the most damaging disease to coffee worldwide. The pathogen has recently appeared in multiple outbreaks in coffee producing countries resulting in significant yield losses and increases in costs related to its control. New races/isolates are constantly emerging as evidenced by the presence of the fungus in plants that were previously resistant. Genomic studies are opening new avenues for the study of the evolution of pathogens, the detailed description of plant-pathogen interactions and the development of molecular techniques for the identification of individual isolates. For this purpose we sequenced 8 different H. vastatrix isolates using NGS technologies and gathered partial genome assemblies due to the large repetitive content in the coffee rust hybrid genome; 74.4% of the assembled contigs harbor repetitive sequences. A hybrid assembly of 333 Mb was built based on the 8 isolates; this assembly was used for subsequent analyses. Analysis of the conserved gene space showed that the hybrid H. vastatrix genome, though highly fragmented, had a satisfactory level of completion with 91.94% of core protein-coding orthologous genes present. RNA-Seq from urediniospores was used to guide the de novo annotation of the H. vastatrix gene complement. In total, 14,445 genes organized in 3921 families were uncovered; a considerable proportion of the predicted proteins (73.8%) were homologous to other Pucciniales species genomes. Several gene families related to the fungal lifestyle were identified, particularly 483 predicted secreted proteins that represent candidate effector genes and will provide interesting hints to decipher virulence in the coffee rust fungus. The genome sequence of Hva will serve as a template to understand the molecular mechanisms used by this fungus to attack the coffee plant, to study the diversity of this species and for the development of molecular markers to distinguish races/isolates.

Effector proteins of rust fungi

Rust fungi include many species that are devastating crop pathogens. To develop resistant plants, a better understanding of rust virulence factors, or effector proteins, is needed. Thus far, only six rust effector proteins have been described: AvrP123, AvrP4, AvrL567, AvrM, RTP1, and PGTAUSPE-10-1. Although some are well established model proteins used to investigate mechanisms of immune receptor activation (avirulence activities) or entry into plant cells, how they work inside host tissues to promote fungal growth remains unknown. The genome sequences of four rust fungi (two Melampsoraceae and two Pucciniaceae) have been analyzed so far. Genome-wide analyses of these species, as well as transcriptomics performed on a broader range of rust fungi, revealed hundreds of small secreted proteins considered as rust candidate secreted effector proteins (CSEPs). The rust community now needs high-throughput approaches (effectoromics) to accelerate effector discovery/characterization and to better understand how they function in planta. However, this task is challenging due to the non-amenability of rust pathosystems (obligate biotrophs infecting crop plants) to traditional molecular genetic approaches mainly due to difficulties in culturing these species in vitro. The use of heterologous approaches should be promoted in the future.

Transcriptome analysis of poplar rust telia reveals overwintering adaptation and tightly coordinated karyogamy and meiosis processes

Most rust fungi have a complex life cycle involving up to five different spore-producing stages. The telial stage that produces melanized overwintering teliospores is one of these and plays a fundamental role for generating genetic diversity as karyogamy and meiosis occur at that stage. Despite the importance of telia for the rust life cycle, almost nothing is known about the fungal genetic programs that are activated in this overwintering structure. In the present study, the transcriptome of telia produced by the poplar rust fungus Melampsora larici-populina has been investigated using whole genome exon oligoarrays and RT-qPCR. Comparative expression profiling at the telial and uredinial stages identifies genes specifically expressed or up-regulated in telia including osmotins/thaumatin-like proteins (TLPs) and aquaporins that may reflect specific adaptation to overwintering as well numerous lytic enzymes acting on plant cell wall, reflecting extensive cell wall remodeling at that stage. The temporal dynamics of karyogamy was followed using combined RT-qPCR and DAPI-staining approaches. This reveals that fusion of nuclei and induction of karyogamy-related genes occur simultaneously between the 25 and 39 days post inoculation time frame. Transcript profiling of conserved meiosis genes indicates a preferential induction right after karyogamy and corroborates that meiosis begins prior to overwintering and is interrupted in Meiosis I (prophase I, diplonema stage) until teliospore germination in early spring.

Identification of candidate effector genes in the transcriptome of the rice root knot nematode Meloidogyne graminicola

Plant-parasitic nematodes secrete so-called effectors into their host plant which are able to suppress the plant's defence responses, alter plant signalling pathways and, in the case of root knot nematodes, induce the formation of giant cells. Putative effectors have been successfully identified by genomics, transcriptomics and proteomics approaches. In this study, we investigated the transcriptome of the rice root knot nematode Meloidogyne graminicola by 454 sequencing of second-stage juveniles as well as mRNA-seq of rice infected tissue. Over 350 000 reads derived from M. graminicola preparasitic juveniles were assembled, annotated and checked for homologues in different databases. From infected rice tissue, 1.4% of all reads generated were identified as being derived from the nematode. Using multiple strategies, several putative effector genes were identified, both pioneer genes and genes corresponding to already known effectors. To check whether these genes could be involved in the interaction with the plant, in situ hybridization was performed on a selection of genes to localize their expression in the nematode. Most were expressed in the gland cells or amphids of the nematode, confirming possible secretion of the proteins and hence a role in infection. Other putative effectors showed a different expression pattern, potentially linked with the excretory/secretory system. This transcriptome study is a good starting point to functionally investigate novel effectors derived from M. graminicola. This will lead to better insights into the interaction between these nematodes and the model plant rice. Moreover, the transcriptome can be used to identify possible target genes for RNA interference (RNAi)-based control strategies. Four genes proved to be interesting targets by showing up to 40% higher mortality relative to the control treatment when soaked in gene-specific small interfering RNAs (siRNAs).

Identification of candidate effector genes in the transcriptome of the rice root knot nematode Meloidogyne graminicola

Plant-parasitic nematodes secrete so-called effectors into their host plant which are able to suppress the plant's defence responses, alter plant signalling pathways and, in the case of root knot nematodes, induce the formation of giant cells. Putative effectors have been successfully identified by genomics, transcriptomics and proteomics approaches. In this study, we investigated the transcriptome of the rice root knot nematode Meloidogyne graminicola by 454 sequencing of second-stage juveniles as well as mRNA-seq of rice infected tissue. Over 350 000 reads derived from M. graminicola preparasitic juveniles were assembled, annotated and checked for homologues in different databases. From infected rice tissue, 1.4% of all reads generated were identified as being derived from the nematode. Using multiple strategies, several putative effector genes were identified, both pioneer genes and genes corresponding to already known effectors. To check whether these genes could be involved in the interaction with the plant, in situ hybridization was performed on a selection of genes to localize their expression in the nematode. Most were expressed in the gland cells or amphids of the nematode, confirming possible secretion of the proteins and hence a role in infection. Other putative effectors showed a different expression pattern, potentially linked with the excretory/secretory system. This transcriptome study is a good starting point to functionally investigate novel effectors derived from M. graminicola. This will lead to better insights into the interaction between these nematodes and the model plant rice. Moreover, the transcriptome can be used to identify possible target genes for RNA interference (RNAi)-based control strategies. Four genes proved to be interesting targets by showing up to 40% higher mortality relative to the control treatment when soaked in gene-specific small interfering RNAs (siRNAs).

Analysis of the leaf transcriptome of Musa acuminata during interaction with Mycosphaerella musicola: gene assembly, annotation and marker development

BACKGROUND

Although banana (Musa sp.) is an important edible crop, contributing towards poverty alleviation and food security, limited transcriptome datasets are available for use in accelerated molecular-based breeding in this genus. 454 GS-FLX Titanium technology was employed to determine the sequence of gene transcripts in genotypes of Musa acuminata ssp. burmannicoides Calcutta 4 and M. acuminata subgroup Cavendish cv. Grande Naine, contrasting in resistance to the fungal pathogen Mycosphaerella musicola, causal organism of Sigatoka leaf spot disease. To enrich for transcripts under biotic stress responses, full length-enriched cDNA libraries were prepared from whole plant leaf materials, both uninfected and artificially challenged with pathogen conidiospores.

RESULTS

The study generated 846,762 high quality sequence reads, with an average length of 334 bp and totalling 283 Mbp. De novo assembly generated 36,384 and 35,269 unigene sequences for M. acuminata Calcutta 4 and Cavendish Grande Naine, respectively. A total of 64.4% of the unigenes were annotated through Basic Local Alignment Search Tool (BLAST) similarity analyses against public databases.Assembled sequences were functionally mapped to Gene Ontology (GO) terms, with unigene functions covering a diverse range of molecular functions, biological processes and cellular components. Genes from a number of defense-related pathways were observed in transcripts from each cDNA library. Over 99% of contig unigenes mapped to exon regions in the reference M. acuminata DH Pahang whole genome sequence. A total of 4068 genic-SSR loci were identified in Calcutta 4 and 4095 in Cavendish Grande Naine. A subset of 95 potential defense-related gene-derived simple sequence repeat (SSR) loci were validated for specific amplification and polymorphism across M. acuminata accessions. Fourteen loci were polymorphic, with alleles per polymorphic locus ranging from 3 to 8 and polymorphism information content ranging from 0.34 to 0.82.

CONCLUSIONS

A large set of unigenes were characterized in this study for both M. acuminata Calcutta 4 and Cavendish Grande Naine, increasing the number of public domain Musa ESTs. This transcriptome is an invaluable resource for furthering our understanding of biological processes elicited during biotic stresses in Musa. Gene-based markers will facilitate molecular breeding strategies, forming the basis of genetic linkage mapping and analysis of quantitative trait loci.

The rust transferred proteins-a new family of effector proteins exhibiting protease inhibitor function

Only few fungal effectors have been described to be delivered into the host cell during obligate biotrophic interactions. RTP1p, from the rust fungi Uromyces fabae and U. striatus, was the first fungal protein for which localization within the host cytoplasm could be demonstrated directly. We investigated the occurrence of RTP1 homologues in rust fungi and examined the structural and biochemical characteristics of the corresponding gene products. The analysis of 28 homologues showed that members of the RTP family are most likely to occur ubiquitously in rust fungi and to be specific to the order Pucciniales. Sequence analyses indicated that the structure of the RTPp effectors is bipartite, consisting of a variable N-terminus and a conserved and structured C-terminus. The characterization of Uf-RTP1p mutants showed that four conserved cysteine residues sustain structural stability. Furthermore, the C-terminal domain exhibits similarities to that of cysteine protease inhibitors, and it was shown that Uf-RTP1p and Us-RTP1p are able to inhibit proteolytic activity in Pichia pastoris culture supernatants. We conclude that the RTP1p homologues constitute a rust fungi-specific family of modular effector proteins comprising an unstructured N-terminal domain and a structured C-terminal domain, which exhibit protease inhibitory activity possibly associated with effector function during biotrophic interactions.

Transcriptome analysis of symptomatic and recovered leaves of geminivirus-infected pepper (Capsicum annuum)

BACKGROUND

Geminiviruses are a large and important family of plant viruses that infect a wide range of crops throughout the world. The Begomovirus genus contains species that are transmitted by whiteflies and are distributed worldwide causing disease on an array of horticultural crops. Symptom remission, in which newly developed leaves of systemically infected plants exhibit a reduction in symptom severity (recovery), has been observed on pepper (Capsicum annuum) plants infected with Pepper golden mosaic virus (PepGMV). Previous studies have shown that transcriptional and post-transcriptional gene silencing mechanisms are involved in the reduction of viral nucleic acid concentration in recovered tissue. In this study, we employed deep transcriptome sequencing methods to assess transcriptional variation in healthy (mock), symptomatic, and recovered pepper leaves following PepGMV infection.

RESULTS

Differential expression analyses of the pepper leaf transcriptome from symptomatic and recovered stages revealed a total of 309 differentially expressed genes between healthy (mock) and symptomatic or recovered tissues. Computational prediction of differential expression was validated using quantitative reverse-transcription PCR confirming the robustness of our bioinformatic methods. Within the set of differentially expressed genes associated with the recovery process were genes involved in defense responses including pathogenesis-related proteins, reactive oxygen species, systemic acquired resistance, jasmonic acid biosynthesis, and ethylene signaling. No major differences were found when compared the differentially expressed genes in symptomatic and recovered tissues. On the other hand, a set of genes with novel roles in defense responses was identified including genes involved in histone modification. This latter result suggested that post-transcriptional and transcriptional gene silencing may be one of the major mechanisms involved in the recovery process. Genes orthologous to the C. annuum proteins involved in the pepper-PepGMV recovery response were identified in both Solanum lycopersicum and Solanum tuberosum suggesting conservation of components of the viral recovery response in the Solanaceae.

CONCLUSION

These data provide a valuable source of information for improving our understanding of the underlying molecular mechanisms by which pepper leaves become symptomless following infection with geminiviruses. The identification of orthologs for the majority of genes differentially expressed in recovered tissues in two major solanaceous crop species provides the basis for future comparative analyses of the viral recovery process across related taxa.