Plant Pathogenomics

Background. Plant-pathogenic oomycetes are responsible for economically important losses on crops worldwide. Phytophthora palmivora, a broad-host-range tropical relative of the potato late blight pathogen, causes rotting diseases in many important tropical crops including papaya, cocoa, oil palm, black pepper, rubber, coconut, durian, mango, cassava and citrus. Transcriptomics have helped to identify repertoires of host-translocated microbial effector proteins which counteract defenses and reprogram the host in support of infection. As such, these studies have helped understanding of how pathogens cause diseases. Despite the importance of P. palmivora diseases, genetic resources to allow for disease resistance breeding and identification of microbial effectors are scarce.

Results. We employed the model plant Nicotiana benthamiana to study the P. palmivora root infections at the cellular and molecular level. Time-resolved dual transcriptomics revealed different pathogen and host transcriptome dynamics. De novo assembly of P. palmivora transcriptome and semi-automated prediction and annotation of the secretome enabled robust identification of conserved infection-promoting effectors. We show that one of them, REX3, suppresses plant secretion processes. In a survey for early transcriptionally activated plant genes we identified a N. benthamiana gene specifically induced at infected root tips that encodes a peptide with danger-associated molecular features.

Conclusions. These results constitute a major advance in our understanding of P. palmivora diseases and establish extensive resources for P. palmivora pathogenomics, effector-aided resistance breeding and the generation of induced resistance to Phytophthora root infections. Furthermore, our approach to find infection relevant secreted genes is transferable to other pathogen-host interactions and not restricted to plants.

The oomycete Phytophthora infestans was the causal agent of the Irish Great Famine and is a recurring threat to global food security. The pathogen can reproduce both sexually and asexually, with high potential to adapt to various environments and great risk to break disease resistance genes in potato. As other oomycetes, P. infestans is regarded to be diploid during the vegetative phase of its life cycle, although some studies reported trisomy, and polyploidy. Using microsatellite fingerprinting, genome-wide assessment of SNP polymorphism, nuclear DNA quantification, and microscopic counting of chromosome numbers we assessed the ploidy level of isolates. All progeny from sexual populations of P. infestans in nature were found to be diploid, in contrast nearly all dominant asexual lineages, including the most important pandemic clonal lineages US-1 and 13_A2 were triploid. Such triploids possess significantly more allelic variation than diploids. We observed that triploid genotype can change to a diploid genome constitution when exposed to artificial stress conditions. This study reveals that fluctuations in the ploidy level maybe a key factor in the adaptation process of this notorious plant destroyer and imposes an extra challenge to control this disease.

Since the late 1990s, UK farmers growing barley have seen the yields and quality of their harvests hurt by an emerging disease called Ramularia leaf spot. The disease is caused by the pathogenic fungus Ramularia collo-cygni. Now a team of scientists studying this fungus have sequenced and explored its genome. Their work helps to illuminate how the fungus causes disease in barley, and enables future studies to investigate why Ramularia leaf spot has become a threat to barley production and a serious economic problem. The scientists performed the work at Scotland’s Rural College (SRUC), The University of Edinburgh's Institute of Evolutionary Biology and Edinburgh Genomics facility, and Rothamsted Research, which receives strategic funding from BBSRC. The study is published this week in the journal BMC Genomics.

Pathogenicity islands are sets of successive genes in a genome that determine the virulence of a bacterium. In a growing number of studies, bacterial virulence appears to be determined by multiple islands scattered along the genome. This is the case in a family of seven plant pathogens and a human pathogen that, under KdgR regulation, massively secrete enzymes such as pectinases that degrade plant cell wall. Here we show that their multiple pathogenicity islands form together a coherently organized, single “archipelago” at the genome scale. Furthermore, in half of the species, most genes encoding secreted pectinases are expressed from the same DNA strand (transcriptional co-orientation). This genome architecture favors DNA conformations that are conducive to genes spatial co-localization, sometimes complemented by co-orientation. As proteins tend to be synthetized close to their encoding genes in bacteria, we propose that this architecture would favor the efficient funneling of pectinases at convergent points within the cell. The underlying functional hypothesis is that this convergent funneling of the full blend of pectinases constitutes a crucial strategy for successful degradation of the plant cell wall. Altogether, our work provides a new approach to describe and predict, at the genome scale, the full virulence complement.

Sequencing unsampled fungal diversity

Phytophthora megakarya (Pmeg) and P. palmivora (Ppal) are closely related species causing cacao black pod rot. While Ppal is a cosmopolitan pathogen, cacao is the only known host of economic importance for Pmeg. Pmeg is more virulent on cacao than Ppal. We sequenced and compared the Pmeg and Ppal genomes and identified virulence-related putative gene models (PGeneM) that may be responsible for their differences in host specificities and virulence. Pmeg and Ppal have estimated genome sizes of 126.88 and 151.23 Mb and PGeneM numbers of 42,036 and 44,327, respectively. The evolutionary histories of Pmeg and Ppal appear quite different. Post-speciation, Ppal underwent whole-genome duplication whereas Pmeg has undergone selective increases in PGeneM numbers, likely through accelerated transposable element driven duplications. Many PGeneMs in both species failed to match transcripts and may represent pseudogenes or cryptic genetic reservoirs. Pmeg appears to have amplified specific gene families, some of which are virulence-related. Analysis of mycelium, zoospore and in planta transcriptome expression profiles using neural network self-organizing map analysis generated 24 multi-variate and non-linear self-organizing map classes. Many members of the RxLR, NPP, and pectinase genes families were specifically induced in planta. Pmeg displays a diverse virulence-related gene complement similar in size to and potentially of greater diversity than Ppal but it remains likely that the specific functions of the genes determine each species’ unique characteristics as pathogens.

Background - Oomycetes are a group of fungus-like eukaryotes with diverse microorganisms living in marine, freshwater and terrestrial environments. Many of them are important pathogens of plants and animals, causing severe economic losses. Based on previous study, gene expression in eukaryotic cells is regulated by epigenetic mechanisms such as DNA methylation and histone modification. However, little is known about epigenetic mechanisms of oomycetes.

Results - In this study, we investigated the candidate genes in regulating histone acetylation in oomycetes genomes through bioinformatics approaches and identified a group of diverse histone acetyltransferases (HATs) and histone deacetylases (HDACs), along with three putative novel HATs. Phylogenetic analyses suggested that most of these oomycetes HATs and HDACs derived from distinct evolutionary ancestors. Phylogenetic based analysis revealed the complex and distinct patterns of duplications and losses of HATs and HDACs in oomycetes. Moreover, gene expression analysis unveiled the specific expression patterns of the 33 HATs and 11 HDACs of Phytophthora infestans during the stages of development, infection and stress response.

Conclusions - In this study, we reveal the structure, diversity and the phylogeny of HATs and HDACs of oomycetes. By analyzing the expression data, we provide an overview of the specific biological stages of these genes involved. Our datasets provide useful inputs to help explore the epigenetic mechanisms and the relationship between genomes and phenotypes of oomycetes.

Smut fungi are plant pathogens mostly parasitizing wild species of grasses as well as domesticated cereal crops. Genome analysis of several smut fungi including Ustilago maydis revealed a singular clustered organization of genes encoding secreted effectors. In U. maydis, many of these clusters have a role in virulence. Reconstructing the evolutionary history of clusters of effector genes is difficult because of their intrinsically fast evolution, which erodes the phylogenetic signal and homology relationships. Here, we describe the use of comparative evolutionary analyses of quality draft assemblies of genomes to study the mechanisms of this evolution. We report the genome sequence of a South African isolate of Sporisorium scitamineum, a smut fungus parasitizing sugar cane with a phylogenetic position intermediate to the two previously sequenced species U. maydisand Sporisorium reilianum. We show that the genome of S. scitamineumcontains more and larger gene clusters encoding secreted effectors than any previously described species in this group. We trace back the origin of the clusters and find that their evolution is mainly driven by tandem gene duplication. In addition, transposable elements play a major role in the evolution of the clustered genes. Transposable elements are significantly associated with clusters of genes encoding fast evolving secreted effectors. This suggests that such clusters represent a case of genome compartmentalization that restrains the activity of transposable elements on genes under diversifying selection for which this activity is potentially beneficial, while protecting the rest of the genome from its deleterious effect.

Ralstonia solanacearum, a species complex of bacterial plant pathogens divided into four monophyletic phylotypes, causes plant diseases in tropical climates around the world. Some strains exhibit a broad host range on solanaceous hosts, while others are highly host-specific as for example some banana-pathogenic strains. Previous studies showed that transcription activator-like (TAL) effectors from Ralstonia, termed RipTALs, are capable of activating reporter genes in planta, if these are preceded by a matching effector binding element (EBE). RipTALs target DNA via their central repeat domain (CRD), where one repeat pairs with one DNA-base of the given EBE. The repeat variable diresidue dictates base repeat specificity in a predictable fashion, known as the TALE code. In this work, we analyze RipTALs across all phylotypes of the Ralstonia solanacearum species complex. We find that RipTALs are prevalent in phylotypes I and IV but absent from most phylotype III and II strains (10/12, 8/14, 1/24, and 1/5 strains contained a RipTAL, respectively). RipTALs originating from strains of the same phylotype show high levels of sequence similarity (>98%) in the N-terminal and C-terminal regions, while RipTALs isolated from different phylotypes show 47–91% sequence similarity in those regions, giving rise to four RipTAL classes. We show that, despite sequence divergence, the base preference for guanine, mediated by the N-terminal region, is conserved across RipTALs of all classes. Using the number and order of repeats found in the CRD, we functionally sub-classify RipTALs, introduce a new simple nomenclature, and predict matching EBEs for all seven distinct RipTALs identified. We experimentally study RipTAL EBEs and uncover that some RipTALs are able to target the EBEs of other RipTALs, referred to as cross-reactivity. In particular, RipTALs from strains with a broad host range on solanaceous hosts cross-react on each other’s EBEs. Investigation of sequence divergence between RipTAL repeats allows for a reconstruction of repeat array biogenesis, for example through slipped strand mispairing or gene conversion. Using these studies we show how RipTALs of broad host range strains evolved convergently toward a shared target sequence. Finally, we discuss the differences between TALE-likes of plant pathogens in the context of disease ecology.

Three members of the Puccini genus, P. triticina (Pt), P. striiformis f.sp. tritici (Pst), and P. graminis f.sp. tritici (Pgt), cause the most common and often most significant foliar diseases of wheat. While similar in biology and life cycle, each species is uniquely adapted and specialized. The genomes of Pt and Pst were sequenced and compared to that of Pgt to identify common and distinguishing gene content, to determine gene variation among wheat rust pathogens, other rust fungi and basidiomycetes, and to identify genes of significance for infection. Pt had the largest genome of the three, estimated at 135 Mb with expansion due to mobile elements and repeats encompassing 50.9% of contig bases; by comparison repeats occupy 31.5% for Pst and 36.5% for Pgt. We find all three genomes are highly heterozygous, with Pst (5.97 SNPs/kb) nearly twice the level detected in Pt (2.57 SNPs/kb) and that previously reported for Pgt. Of 1,358 predicted effectors in Pt, 784 were found expressed across diverse life cycle stages including the sexual stage. Comparison to related fungi highlighted the expansion of gene families involved in transcriptional regulation and nucleotide binding, protein modification, and carbohydrate enzyme degradation. Two allelic homeodomain, HD1 and HD2, pairs and three pheromone receptor (STE3) mating-type genes were identified in each dikaryotic Pucciniaspecies. The HD proteins were active in a heterologous Ustilago maydis mating assay and host induced gene silencing of the HD and STE3 alleles reduced wheat host infection.

Background. Understanding how plants and pathogens modulate gene expression during the host-pathogen interaction is key to uncovering the molecular mechanisms that regulate disease progression. Recent advances in sequencing technologies have provided new opportunities to decode the complexity of such interactions. In this study, we used an RNA-based sequencing approach (RNA-seq) to assess the global expression profiles of the wheat yellow rust pathogen Puccinia striiformis f. sp. tritici (PST) and its host during infection.

Results. We performed a detailed RNA-seq time-course for a susceptible and a resistant wheat host infected with PST. This study (i) defined the global gene expression profiles for PST and its wheat host, (ii) substantially improved the gene models for PST, (iii) evaluated the utility of several programmes for quantification of global gene expression for PST and wheat, and (iv) identified clusters of differentially expressed genes in the host and pathogen. By focusing on components of the defence response in susceptible and resistant hosts, we were able to visualise the effect of PST infection on the expression of various defence components and host immune receptors.

Conclusions. Our data showed sequential, temporally coordinated activation and suppression of expression of a suite of immune-response regulators that varied between compatible and incompatible interactions. These findings provide the framework for a better understanding of how PST causes disease and support the idea that PST can suppress the expression of defence components in wheat to successfully colonize a susceptible host.

Background - The protist Plasmodiophora brassicae is a soil-borne pathogen of cruciferous species and the causal agent of clubroot disease of Brassicas including agriculturally important crops such as canola/rapeseed (Brassica napus). P. brassicae has remained an enigmatic plant pathogen and is a rare example of an obligate biotroph that resides entirely inside the host plant cell. The pathogen is the cause of severe yield losses and can render infested fields unsuitable for Brassica crop growth due to the persistence of resting spores in the soil for up to 20 years.

Results - To provide insight into the biology of the pathogen and its interaction with its primary host B. napus, we produced a draft genome of P. brassicae pathotypes 3 and 6 (Pb3 and Pb6) that differ in their host range. Pb3 is highly virulent on B. napus (but also infects other Brassica species) while Pb6 infects only vegetable Brassica crops. Both the Pb3 and Pb6 genomes are highly compact, each with a total size of 24.2 Mb, and contain less than 2 % repetitive DNA. Clustering of genome-wide single nucleotide polymorphisms (SNP) of Pb3, Pb6 and three additional re-sequenced pathotypes (Pb2, Pb5 and Pb8) shows a high degree of correlation of cluster grouping with host range. The Pb3 genome features significant reduction of intergenic space with multiple examples of overlapping untranslated regions (UTRs). Dependency on the host for essential nutrients is evident from the loss of genes for the biosynthesis of thiamine and some amino acids and the presence of a wide range of transport proteins, including some unique to P. brassicae. The annotated genes of Pb3 include those with a potential role in the regulation of the plant growth hormones cytokinin and auxin. The expression profile of Pb3 genes, including putative effectors, during infection and their potential role in manipulation of host defence is discussed.

Conclusion - The P. brassicae genome sequence reveals a compact genome, a dependency of the pathogen on its host for some essential nutrients and a potential role in the regulation of host plant cytokinin and auxin. Genome annotation supported by RNA sequencing reveals significant reduction in intergenic space which, in addition to low repeat content, has likely contributed to the P. brassicae compact genome.