Diploid Nuclei Occur throughout the Life Cycles of Pucciniales Fungi.

Within Eukaryotes, fungi are the typical representatives of haplontic life cycles. Basidiomycota fungi are dikaryotic in extensive parts of their life cycle, but diploid nuclei are known to form only in basidia. Among Basidiomycota, the Pucciniales are notorious for presenting the most complex life cycles, with high host specialization, and for their expanded genomes. Using cytogenomic (flow cytometry and cell sorting on propidium iodide-stained nuclei) and cytogenetic (FISH with rDNA probe) approaches, we report the widespread occurrence of replicating haploid and diploid nuclei (i.e., 1C, 2C and a small proportion of 4C nuclei) in diverse life cycle stages (pycnial, aecial, uredinial, and telial) of all 35 Pucciniales species analyzed, but not in sister taxa. These results suggest that the Pucciniales life cycle is distinct from any cycle known, i.e., neither haplontic, diplontic nor haplodiplontic, corroborating patchy and disregarded previous evidence. However, the biological basis and significance of this phenomenon remain undisclosed. IMPORTANCE Within Eukaryotes, fungi are the typical representatives of haplontic life cycles, contrasting with plants and animals. As such, fungi thus contain haploid nuclei throughout their life cycles, with sexual reproduction generating a single diploid cell upon karyogamy that immediately undergoes meiosis, thus resuming the haploid cycle. In this work, using cytogenetic and cytogenomic tools, we demonstrate that a vast group of fungi presents diploid nuclei throughout their life cycles, along with haploid nuclei, and that both types of nuclei replicate. Moreover, haploid nuclei are absent from urediniospores. The phenomenon appears to be transversal to the organisms in the order Pucciniales (rust fungi) and it does not occur in neighboring taxa, but a biological explanation or function for it remains elusive.

Pathological, Morphological, Cytogenomic, Biochemical and Molecular Data Support the Distinction between Colletotrichum cigarro comb. et stat. nov. and Colletotrichum kahawae.

The genus Colletotrichum has witnessed tremendous variations over the years in the number of species recognized, ranging from 11 to several hundreds. Host-specific fungal species, once the rule, are now the exception, with polyphagous behavior regarded as normal in this genus. The species Colletotrichum kahawae was created to accommodate the pathogens that have the unique ability to infect green developing coffee berries causing the devastating Coffee Berry Disease in Africa, but its close phylogenetic relationship to a polyphagous group of fungi in the C. gloeosporioides species complex led some researchers to regard these pathogens as members of a wider species. In this work we combine pathological, morphological, cytogenomic, biochemical, and molecular data of a comprehensive set of phylogenetically-related isolates to show that the Coffee Berry Disease pathogen forms a separate species, C. kahawae, and also to assign the closely related fungi, previously in C. kahawae subsp. cigarro, to a new species, C. cigarro comb. et stat. nov. This taxonomic clarification provides an opportunity to link phylogeny and functional biology, and additionally enables a much-needed tool for plant pathology and agronomy, associating exclusively C. kahawae to the Coffee Berry Disease pathogen.

Dietary Antioxidants in Coffee Leaves: Impact of Botanical Origin and Maturity on Chlorogenic Acids and Xanthones.

Natural polyphenols are important dietary antioxidants that significantly benefit human health. Coffee and tea have been shown to largely contribute to the dietary intake of these antioxidants in several populations. More recently, the use of coffee leaves to produce tea has become a potential commercial target, therefore prompting studies on the quantification of polyphenols in coffee leaves. In this study a variety of coffee leaf species, at different development stages, were analyzed using ultra-high pressure liquid chromatography. The results demonstrate that both the botanical origin of the samples and their maturity influence significantly the concentration of the antioxidants; for total chlorogenic acids a two-fold difference was found between different species and up to a three-fold variation was observed between young and mature leaves. Furthermore, the range of concentrations of chlorogenic acids in young leaves (35.7-80.8 mg/g of dry matter) were found to be comparable to the one reported for green coffee beans. The results provide important data from which potential new commercial products can be developed.

A first insight into the involvement of phytohormones pathways in coffee resistance and susceptibility to Colletotrichum kahawae.

Understanding the molecular mechanisms underlying coffee-pathogen interactions are of key importance to aid disease resistance breeding efforts. In this work the expression of genes involved in salicylic acid (SA), jasmonic acid (JA) and ethylene (ET) pathways were studied in hypocotyls of two coffee varieties challenged with the hemibiotrophic fungus Colletotrichum kahawae, the causal agent of Coffee Berry Disease. Based on a cytological analysis, key time-points of the infection process were selected and qPCR was used to evaluate the expression of phytohormones biosynthesis, reception and responsive-related genes. The resistance to C. kahawae was characterized by restricted fungal growth associated with early accumulation of phenolic compounds in the cell walls and cytoplasmic contents, and deployment of hypersensitive reaction. Similar responses were detected in the susceptible variety, but in a significantly lower percentage of infection sites and with no apparent effect on disease development. Gene expression analysis suggests a more relevant involvement of JA and ET phytohormones than SA in this pathosystem. An earlier and stronger activation of the JA pathway observed in the resistant variety, when compared with the susceptible one, seems to be responsible for the successful activation of defense responses and inhibition of fungal growth. For the ET pathway, the down or non-regulation of ET receptors in the resistant variety, together with a moderate expression of the responsive-related gene ERF1, indicates that this phytohormone may be related with other functions besides the resistance response. However, in the susceptible variety, the stronger activation of ERF1 gene at the beginning of the necrotrophic phase, suggests the involvement of ET in tissue senescence. As far as we know, this is the first attempt to unveil the role of phytohormones in coffee-C. kahawae interactions, thus contributing to deepen our understanding on the complex mechanisms of plant signaling and defense.

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 Validation of Conventional and RNA-Seq Data-Derived Reference Genes for qPCR Expression Studies of Colletotrichum kahawae.

Colletotrichum kahawae is an emergent fungal pathogen causing severe epidemics of Coffee Berry Disease on Arabica coffee crops in Africa. Currently, the molecular mechanisms underlying the Coffea arabica-C. kahawae interaction are still poorly understood, as well as the differences in pathogen aggressiveness, which makes the development of functional studies for this pathosystem a crucial step. Quantitative real time PCR (qPCR) has been one of the most promising approaches to perform gene expression analyses. However, proper data normalization with suitable reference genes is an absolute requirement. In this study, a set of 8 candidate reference genes were selected based on two different approaches (literature and Illumina RNA-seq datasets) to assess the best normalization factor for qPCR expression analysis of C. kahawae samples. The gene expression stability of candidate reference genes was evaluated for four isolates of C. kahawae bearing different aggressiveness patterns (Ang29, Ang67, Zim12 and Que2), at different stages of fungal development and key time points of the plant-fungus interaction process. Gene expression stability was assessed using the pairwise method incorporated in geNorm and the model-based method used by NormFinder software. For C. arabica-C. kahawae interaction samples, the best normalization factor included the combination of PP1, Act and ck34620 genes, while for C. kahawae samples the combination of PP1, Act and ck20430 revealed to be the most appropriate choice. These results suggest that RNA-seq analyses can provide alternative sources of reference genes in addition to classical reference genes. The analysis of expression profiles of bifunctional catalase-peroxidase (cat2) and trihydroxynaphthalene reductase (thr1) genes further enabled the validation of the selected reference genes. This study provides, for the first time, the tools required to conduct accurate qPCR studies in C. kahawae considering its aggressiveness pattern, developmental stage and host interaction.

Genomic Patterns of Positive Selection at the Origin of Rust Fungi.

Understanding the origin and evolution of pathogenicity and biotrophic life-style of rust fungi has remained a conundrum for decades. Research on the molecular mechanisms responsible for rust fungi evolution has been hampered by their biotrophic life-style until the sequencing of some rust fungi genomes. With the availability of multiple whole genomes and EST data for this group, it is now possible to employ genome-wide surveys and investigate how natural selection shaped their evolution. In this work, we employed a phylogenomics approach to search for positive selection and genes undergoing accelerated evolution at the origin of rust fungi on an assembly of single copy genes conserved across a broad range of basidiomycetes. Up to 985 genes were screened for positive selection on the phylogenetic branch leading to rusts, revealing a pervasive signal of positive selection throughout the data set with the proportion of positively selected genes ranging between 19.6-33.3%. Additionally, 30 genes were found to be under accelerated evolution at the origin of rust fungi, probably due to a mixture of positive selection and relaxation of purifying selection. Functional annotation of the positively selected genes revealed an enrichment in genes involved in the biosynthesis of secondary metabolites and several metabolism and transporter classes.

A method for obtaining RNA from Hemileia vastatrix appressoria produced in planta, suitable for transcriptomic analyses.

Appressoria are the first infection structures developed by rust fungi and require specific topographic signals from the host for their differentiation. The ease in obtaining appressoria in vitro for these biotrophic fungi led to studies concerning gene expression and gene discovery at appressorial level, avoiding the need to distinguish plant and fungal transcripts. However, in some pathosystems, it was observed that gene expression in appressoria seems to be influenced by host-derived signals, suggesting that transcriptomic analyses performed from in planta differentiated appressoria would be potentially more informative than those from in vitro differentiated appressoria. Nevertheless analysing appressorial RNA obtained from in planta samples is often hampered by an excessive dilution of fungal RNA within plant RNA, besides uncertainty regarding the fungal or plant origin of RNA from highly conserved genes. To circumvent these difficulties, we have recovered Hemileia vastatrix appressoria from Arabica coffee leaf surface using a film of nitrocellulose dissolved in butyl and ethyl acetates (nail polish), and extracted fungal RNA from the polish peel. RNA thus obtained is of good quality and usable for cDNA synthesis and transcriptomic (quantitative PCR) studies. This method could provide the means to investigate specific host-induced appressoria-related fungal pathogenicity factors.

Genome size analyses of Pucciniales reveal the largest fungal genomes.

Rust fungi (Basidiomycota, Pucciniales) are biotrophic plant pathogens which exhibit diverse complexities in their life cycles and host ranges. The completion of genome sequencing of a few rust fungi has revealed the occurrence of large genomes. Sequencing efforts for other rust fungi have been hampered by uncertainty concerning their genome sizes. Flow cytometry was recently applied to estimate the genome size of a few rust fungi, and confirmed the occurrence of large genomes in this order (averaging 225.3 Mbp, while the average for Basidiomycota was 49.9 Mbp and was 37.7 Mbp for all fungi). In this work, we have used an innovative and simple approach to simultaneously isolate nuclei from the rust and its host plant in order to estimate the genome size of 30 rust species by flow cytometry. Genome sizes varied over 10-fold, from 70 to 893 Mbp, with an average genome size value of 380.2 Mbp. Compared to the genome sizes of over 1800 fungi, Gymnosporangium confusum possesses the largest fungal genome ever reported (893.2 Mbp). Moreover, even the smallest rust genome determined in this study is larger than the vast majority of fungal genomes (94%). The average genome size of the Pucciniales is now of 305.5 Mbp, while the average Basidiomycota genome size has shifted to 70.4 Mbp and the average for all fungi reached 44.2 Mbp. Despite the fact that no correlation could be drawn between the genome sizes, the phylogenomics or the life cycle of rust fungi, it is interesting to note that rusts with Fabaceae hosts present genomes clearly larger than those with Poaceae hosts. Although this study comprises only a small fraction of the more than 7000 rust species described, it seems already evident that the Pucciniales represent a group where genome size expansion could be a common characteristic. This is in sharp contrast to sister taxa, placing this order in a relevant position in fungal genomics research.

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.

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.

Coffee (Coffea arabica L.), one of the key export and cash crops in tropical and subtropical countries, suffers severe losses from the rust fungus Hemileia vastatrix. The transcriptome of H. vastatrix was analysed during a compatible interaction with coffee to obtain an exhaustive repertoire of the genes expressed during infection and to identify potential effector genes. Large-scale sequencing (454-GS-FLEX Titanium) of mixed coffee and rust cDNAs obtained from 21-day rust-infected leaves generated 352 146 sequences which assembled into 22 774 contigs. In the absence of any reference genomic sequences for Coffea or Hemileia, specific trinucleotide frequencies within expressed sequence tags (ESTs) and blast homology against a set of dicots and basidiomycete genomes were used to distinguish pathogen from plant sequences. About 30% (6763) of the contigs were assigned to H. vastatrix and 61% (13 951) to C. arabica. The majority (60%) of the rust sequences did not show homology to any genomic database, indicating that they were potential novel fungal genes. In silico analyses of the 6763 H. vastatrix contigs predicted 382 secreted proteins and identified homologues of the flax rust haustorially expressed secreted proteins (HESPs) and bean rust transferred protein 1 (RTP1). These rust candidate effectors showed conserved amino-acid domains and conserved patterns of cysteine positions suggestive of conserved functions during infection of host plants. Quantitative reverse transcription-polymerase chain reaction profiling of selected rust genes revealed dynamic expression patterns during the time course of infection of coffee leaves. This study provides the first valuable genomic resource for the agriculturally important plant pathogen H. vastatrix and the first comprehensive C. arabica EST dataset.

Validation of RT-qPCR reference genes for in planta expression studies in Hemileia vastatrix, the causal agent of coffee leaf rust.

Hemileia vastatrix is a biotrophic fungus, causing coffee leaf rust in all coffee growing countries, leading to serious social and economic problems. Gene expression studies may have a key role unravelling the transcriptomics of this pathogen during interaction with the plant host. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is currently the golden standard for gene expression analysis, although an accurate normalisation is essential for adequate conclusions. Reference genes are often used for this purpose, but the stability of their expression levels requires validation under experimental conditions. Moreover, pathogenic fungi undergo important biomass variations along their infection process in planta, which raises the need for an adequate method to further normalise the proportion of fungal cDNA in the total plant and fungus cDNA pool. In this work, the expression profiles of seven reference genes [glyceraldehyde-3-phosphate dehydrogenase (GADPH), elongation factor (EF-1), Beta tubulin (ß-tubulin), cytochrome c oxidase subunit III (Cyt III), cytochrome b (Cyt b), Hv00099, and 40S ribosomal protein (40S_Rib)] were analysed across 28 samples, obtained in vitro (germinated uredospores and appressoria) and in planta (post-penetration fungal growth phases). Gene stability was assessed using the statistical algorithms incorporated in geNorm and NormFinder tools. Cyt b, 40S_Rib, and Hv00099 were the most stable genes for the in vitro dataset, while 40S_Rib, GADPH, and Cyt III were the most stable in planta. For the combined datasets (in vitro and in planta), 40S_Rib, GADPH, and Hv00099 were selected as the most stable. Subsequent expression analysis for a gene encoding an alpha subunit of a heterotrimeric G-protein showed that the reference genes selected for the combined dataset do not differ significantly from those selected specifically for the in vitro and in planta datasets. Our study provides tools for correct validation of reference genes in obligate biotrophic plant pathogens, as well as the basis for RT-qPCR studies in H. vastatrix.