Genome-wide prediction and functional validation of promoter motifs regulating gene expression in spore and infection stages of Phytophthora infestans

Promoter Analysis and Dissection Using Reporter Genes, Comparative Genomics, and Gel Shift Assays in Phytophthora

Transcriptional regulation allows cells to execute developmental programs, maintain homeostasis, and respond to intra- and extracellular signals. Central to these processes are promoters, which in eukaryotes are sequences upstream of genes that bind transcription factors (TFs) and which recruit RNA polymerase to initiate mRNA synthesis. Valuable tools for studying promoters include reporter genes, which can be used to indicate when and where genes are activated. Moreover, functional regions within promoters (typically TF binding sites) can be identified by integrating reporter assays with promoter mutagenesis. These sites may also be revealed through comparative genomics, or by the DNA-protein binding procedure known as a gel shift or electrophoretic mobility shift assay (EMSA). The latter can also be used to test if a specific TF binds a DNA target or assess the binding kinetics or affinity of the complex. In this chapter, we describe procedures for expressing reporter genes in Phytophthora, assaying reporter activity, identifying functional sites within promoters, and testing purified TFs or proteins within nuclear extracts for DNA binding.

Specific Transcriptional Regulation Controls Plant Organ-Specific Infection by the Oomycete Pathogen Phytophthora sojae

The organs of a plant species vary in cell structure, metabolism and defence responses. However, the mechanisms that enable a single pathogen to colonise different plant organs remain unclear. Here we compared the transcriptome of the oomycete pathogen Phytophthora sojae during infection of roots versus leaves of soybeans. We found differences in the transcript levels of hundreds of pathogenicity-related genes, particularly genes encoding carbohydrate-active enzymes, secreted (effector) proteins, oxidoreductase-related proteins and transporters. To identify the key regulator for root-specific infection, we knocked out root-specific transcription factors (TFs) and found the mutants of PsBZPc29, which encodes a member of an oomycete-specific class of basic leucine zipper (bZIP) TFs, displayed reduced virulence on soybean roots but not on leaves. More than 60% of the root-specific genes showed reduced expression in the mutants during root infection. The results suggest that transcriptional regulation underlies the organ-specific infection by P. sojae, and that a bZIP TF plays a key role in root-specific transcriptional regulation.

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

A MYB-related transcription factor regulates effector gene expression in an oomycete pathogen

Phytophthora pathogens possess hundreds of effector genes that exhibit diverse expression patterns during infection, yet how the expression of effector genes is precisely regulated remains largely elusive. Previous studies have identified a few potential conserved transcription factor binding sites (TFBSs) in the promoters of Phytophthora effector genes. Here, we report a MYB-related protein, PsMyb37, in Phytophthora sojae, the major causal agent of root and stem rot in soybean. Yeast one-hybrid and electrophoretic mobility shift assays showed that PsMyb37 binds to the TACATGTA motif, the most prevalent TFBS in effector gene promoters. The knockout mutant of PsMyb37 exhibited significantly reduced virulence on soybean and was more sensitive to oxidative stress. Consistently, transcriptome analysis showed that numerous effector genes associated with suppressing plant immunity or scavenging reactive oxygen species were down-regulated in the PsMyb37 knockout mutant during infection compared to the wild-type P. sojae. Several promoters of effector genes were confirmed to drive the expression of luciferase in a reporter assay. These results demonstrate that a MYB-related transcription factor contributes to the expression of effector genes in P. sojae.

High-throughput time series expression profiling of Plasmopara halstedii infecting Helianthus annuus reveals conserved sequence motifs upstream of co-expressed genes

Downy mildew disease of sunflower, caused by the obligate biotrophic oomycete Plasmopara halstedii, can have significant economic impact on sunflower cultivation. Using high-throughput whole transcriptome sequencing, four developmental phases in 16 time-points of Pl. halstedii infecting Helianthus annuus were investigated. With the aim of identifying potential functional and regulatory motifs upstream of co-expressed genes, time-series derived gene expression profiles were clustered based on their time-course similarity, and their upstream regulatory gene sequences were analyzed here. Several conserved motifs were found upstream of co-expressed genes, which might be involved in binding specific transcription factors. Such motifs were also found associated with virulence related genes, and could be studied on a genetically tractable model to clarify, if these are involved in regulating different stages of pathogenesis.

Comparative genomics among cyst nematodes reveals distinct evolutionary histories among effector families and an irregular distribution of effector-associated promoter motifs

Potato cyst nematodes (PCNs), an umbrella term used for two species, Globodera pallida and G. rostochiensis, belong worldwide to the most harmful pathogens of potato. Pathotype-specific host plant resistances are essential for PCN control. However, the poor delineation of G. pallida pathotypes has hampered the efficient use of available host plant resistances. Long-read sequencing technology allowed us to generate a new reference genome of G. pallida population D383 and, as compared to the current reference, the new genome assembly is 42 times less fragmented. For comparison of diversification patterns of six effector families between G. pallida and G. rostochiensis, an additional reference genome was generated for an outgroup, the beet cyst nematode Heterodera schachtii (IRS population). Large evolutionary contrasts in effector family topologies were observed. While VAPs (venom allergen-like proteins) diversified before the split between the three cyst nematode species, the families GLAND5 and GLAND13 only expanded in PCNs after their separation from the genus Heterodera. Although DNA motifs in the promoter regions thought to be involved in the orchestration of effector expression ("DOG boxes") were present in all three cyst nematode species, their presence is not a necessity for dorsal gland-produced effectors. Notably, DOG box dosage was only loosely correlated with the expression level of individual effector variants. Comparison of the G. pallida genome with those of two other cyst nematodes underlined the fundamental differences in evolutionary history between effector families. Resequencing of PCN populations with different virulence characteristics will allow for the linking of these characteristics to the composition of the effector repertoire as well as for the mapping of PCN diversification patterns resulting from extreme anthropogenic range expansion.

RNA silencing proteins and small RNAs in oomycete plant pathogens and biocontrol agents

Introduction

Oomycetes cause several damaging diseases of plants and animals, and some species also act as biocontrol agents on insects, fungi, and other oomycetes. RNA silencing is increasingly being shown to play a role in the pathogenicity of Phytophthora species, either through trans-boundary movement of small RNAs (sRNAs) or through expression regulation of infection promoting effectors.

Methods

To gain a wider understanding of RNA silencing in oomycete species with more diverse hosts, we mined genome assemblies for Dicer-like (DCL), Argonaute (AGO), and RNA dependent RNA polymerase (RDRP) proteins from Phytophthora plurivora, Ph. cactorum, Ph. colocasiae, Pythium oligandrum, Py. periplocum, and Lagenidium giganteum. Moreover, we sequenced small RNAs from the mycelium stage in each of these species.

Results and discussion

Each of the species possessed a single DCL protein, but they differed in the number and sequence of AGOs and RDRPs. SRNAs of 21nt, 25nt, and 26nt were prevalent in all oomycetes analyzed, but the relative abundance and 5' base preference of these classes differed markedly between genera. Most sRNAs mapped to transposons and other repeats, signifying that the major role for RNA silencing in oomycetes is to limit the expansion of these elements. We also found that sRNAs may act to regulate the expression of duplicated genes. Other sRNAs mapped to several gene families, and this number was higher in Pythium spp., suggesting a role of RNA silencing in regulating gene expression. Genes for most effector classes were the source of sRNAs of variable size, but some gene families showed a preference for specific classes of sRNAs, such as 25/26 nt sRNAs targeting RxLR effector genes in Phytophthora species. Novel miRNA-like RNAs (milRNAs) were discovered in all species, and two were predicted to target transcripts for RxLR effectors in Ph. plurivora and Ph. cactorum, indicating a putative role in regulating infection. Moreover, milRNAs from the biocontrol Pythium species had matches in the predicted transcriptome of Phytophthora infestans and Botrytis cinerea, and L. giganteum milRNAs matched candidate genes in the mosquito Aedes aegypti. This suggests that trans-boundary RNA silencing may have a role in the biocontrol action of these oomycetes.

ATAC-Seq Reveals the Landscape of Open Chromatin and cis-Regulatory Elements in the Phytophthora sojae Genome

Nucleosome-free open chromatin often harbors transcription factor (TF)-binding sites that are associated with active cis-regulatory elements. However, analysis of open chromatin regions has rarely been applied to oomycete or fungal plant pathogens. In this study, we performed the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to identify open chromatin and cis-regulatory elements in Phytophthora sojae at the mycelial stage. We identified 10,389 peaks representing nucleosome-free regions (NFRs). The peaks were enriched in gene-promoter regions and associated with 40% of P. sojae genes; transcription levels were higher for genes with multiple peaks than genes with a single peak and were higher for genes with a single peak than genes without peak. Chromatin accessibility was positively correlated with gene transcription level. Through motif discovery based on NFR peaks in core promoter regions, 25 candidate cis-regulatory motifs with evidence of TF-binding footprints were identified. These motifs exhibited various preferences for location in the promoter region and associations with the transcription level of their target genes, which included some putative pathogenicity-related genes. As the first study revealing the landscape of open chromatin and the correlation between chromatin accessibility and gene transcription level in oomycetes, the results provide a technical reference and data resources for future studies on the regulatory mechanisms of gene transcription.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Identification of Cis-Regulatory Sequences Controlling Pollen-Specific Expression of Hydroxyproline-Rich Glycoprotein Genes in Arabidopsis thaliana

Hydroxyproline-rich glycoproteins (HRGPs) are a superfamily of plant cell wall structural proteins that function in various aspects of plant growth and development, including pollen tube growth. We have previously characterized protein sequence signatures for three family members in the HRGP superfamily: the hyperglycosylated arabinogalactan-proteins (AGPs), the moderately glycosylated extensins (EXTs), and the lightly glycosylated proline-rich proteins (PRPs). However, the mechanism of pollen-specific HRGP gene expression remains unexplored. To this end, we developed an integrative analysis pipeline combining RNA-seq gene expression and promoter sequences to identify cis-regulatory motifs responsible for pollen-specific expression of HRGP genes in Arabidopsis thaliana. Specifically, we mined the public RNA-seq datasets and identified 13 pollen-specific HRGP genes. Ensemble motif discovery identified 15 conserved promoter elements between A.thaliana and A. lyrata. Motif scanning revealed two pollen related transcription factors: GATA12 and brassinosteroid (BR) signaling pathway regulator BZR1. Finally, we performed a regression analysis and demonstrated that the 15 motifs provided a good model of HRGP gene expression in pollen (R = 0.61). In conclusion, we performed the first integrative analysis of cis-regulatory motifs in pollen-specific HRGP genes, revealing important insights into transcriptional regulation in pollen tissue.

The Genome of Peronospora belbahrii Reveals High Heterozygosity, a Low Number of Canonical Effectors, and TC-Rich Promoters

Along with Plasmopara destructor, Peronosopora belbahrii has arguably been the economically most important newly emerging downy mildew pathogen of the past two decades. Originating from Africa, it has started devastating basil production throughout the world, most likely due to the distribution of infested seed material. Here, we present the genome of this pathogen and results from comparisons of its genomic features to other oomycetes. The assembly of the nuclear genome was around 35.4 Mbp in length, with an N₅₀ scaffold length of around 248 kbp and an L₅₀ scaffold count of 46. The circular mitochondrial genome consisted of around 40.1 kbp. From the repeat-masked genome, 9,049 protein-coding genes were predicted, out of which 335 were predicted to have extracellular functions, representing the smallest secretome so far found in peronosporalean oomycetes. About 16% of the genome consists of repetitive sequences, and, based on simple sequence repeat regions, we provide a set of microsatellites that could be used for population genetic studies of P. belbahrii. P. belbahrii has undergone a high degree of convergent evolution with other obligate parasitic pathogen groups, reflecting its obligate biotrophic lifestyle. Features of its secretome, signaling networks, and promoters are presented, and some patterns are hypothesized to reflect the high degree of host specificity in Peronospora species. In addition, we suggest the presence of additional virulence factors apart from classical effector classes that are promising candidates for future functional studies.

Sex chromosome evolution in parasitic nematodes of humans

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

Promoter tools for further development of Aspergillus oryzae as a platform for fungal secondary metabolite production

Background

The filamentous fungus Aspergillus oryzae is widely used for secondary metabolite production by heterologous expression; thus, a wide variety of promoter tools is necessary to broaden the application of this species. Here we built a procedure to survey A. flavus genes constitutively highly expressed in 83 transcriptome datasets obtained under various conditions affecting secondary metabolite production, to find promoters useful for heterologous expression of genes in A. oryzae.

Results

To test the ability of the promoters of the top 6 genes to induce production of a fungal secondary metabolite, ustiloxin B, we inserted the promoters before the start codon of ustR, which encodes the transcription factor of the gene cluster responsible for ustiloxin B biosynthesis, in A. oryzae. Four of the 6 promoters induced ustiloxin B production in all tested media (solid maize, liquid V8 and PDB media), and also ustR expression. Two of the 4 promoters were those of tef1 and gpdA, which are well characterized in A. oryzae and A. nidulans, respectively, whereas the other two, those of AFLA_030930 and AFLA_113120, are newly reported here and show activities comparable to that of the gpdA promoter with respect to induction of gene expression and ustiloxin B production.

Conclusion

We newly reported two sequences as promoter tools for secondary metabolite production in A. oryzae. Our results demonstrate that our simple strategy of surveying for constitutively highly expressed genes in large-scale transcriptome datasets is useful for finding promoter sequences that can be used as heterologous expression tools in A. oryzae.

Defining Transgene Insertion Sites and Off-Target Effects of Homology-Based Gene Silencing Informs the Application of Functional Genomics Tools in Phytophthora infestans

DNA transformation and homology-based transcriptional silencing are frequently used to assess gene function in Phytophthora spp. Since unplanned side-effects of these tools are not well-characterized, we used P. infestans to study plasmid integration sites and whether knockdowns caused by homology-dependent silencing spread to other genes. Insertions occurred both in gene-dense and gene-sparse regions but disproportionately near the 5' ends of genes, which disrupted native coding sequences. Microhomology at the recombination site between plasmid and chromosome was common. Studies of transformants silenced for 12 different gene targets indicated that neighbors within 500 nt were often cosilenced, regardless of whether hairpin or sense constructs were employed and the direction of transcription of the target. However, this cis spreading of silencing did not occur in all transformants obtained with the same plasmid. Genome-wide studies indicated that unlinked genes with partial complementarity with the silencing-inducing transgene were not usually down-regulated. We learned that hairpin or sense transgenes were not cosilenced with the target in all transformants, which informs how screens for silencing should be performed. We conclude that transformation and gene silencing can be reliable tools for functional genomics in Phytophthora spp. but must be used carefully, especially by testing for the spread of silencing to genes flanking the target.

A High-Quality Grapevine Downy Mildew Genome Assembly Reveals Rapidly Evolving and Lineage-Specific Putative Host Adaptation Genes

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high contiguity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant–pathogen interactions. Plasmopara viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.

Ancestral Function and Diversification of a Horizontally Acquired Oomycete Carboxylic Acid Transporter

Horizontal gene transfer (HGT) can equip organisms with novel genes, expanding the repertoire of genetic material available for evolutionary innovation and allowing recipient lineages to colonize new environments. However, few studies have characterized the functions of HGT genes experimentally or examined postacquisition functional divergence. Here, we report the use of ancestral sequence reconstruction and heterologous expression in Saccharomyces cerevisiae to examine the evolutionary history of an oomycete transporter gene family that was horizontally acquired from fungi. We demonstrate that the inferred ancestral oomycete HGT transporter proteins and their extant descendants transport dicarboxylic acids which are intermediates of the tricarboxylic acid cycle. The substrate specificity profile of the most ancestral protein has largely been retained throughout the radiation of oomycetes, including in both plant and animal pathogens and in a free-living saprotroph, indicating that the ancestral HGT transporter function has been maintained by selection across a range of different lifestyles. No evidence of neofunctionalization in terms of substrate specificity was detected for different HGT transporter paralogues which have different patterns of temporal expression. However, a striking expansion of substrate range was observed for one plant pathogenic oomycete, with a HGT derived paralogue from Pythium aphanidermatum encoding a protein that enables tricarboxylic acid uptake in addition to dicarboxylic acid uptake. This demonstrates that HGT acquisitions can provide functional additions to the recipient proteome as well as the foundation material for the evolution of expanded protein functions.

A High-Quality Grapevine Downy Mildew Genome Assembly Reveals Rapidly Evolving and Lineage-Specific Putative Host Adaptation Genes

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high contiguity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant-pathogen interactions. Plasmopara viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.

The MADS-box Transcription Factor PsMAD1 Is Involved in Zoosporogenesis and Pathogenesis of Phytophthora sojae

Transcriptional regulation is critical for plant pathogen development and virulence. MADS-box transcription factors belong to a highly conserved transcriptional regulator family in eukaryotic organisms that are involved in various important biological processes. Only one predicted MADS-box gene, PsMAD1, was identified in Phytophthora sojae, which was highly expressed during the sporangia and infection stages. To investigate its function, we generated PsMAD1 knockout mutants using the CRISPR/Cas9 system. Compared with the wild-type strain, the mutants showed no changes in vegetative growth, oospore production, or no differences in sensitivity to various abiotic stresses. Although sporangia production was normal, no zoospore release was detected in PsMAD1 mutants. Microscopy analyses revealed failure of cleavage of the cytoplasm into uninucleate zoospores in the mutants. In addition, the mutants showed reduced virulence in soybean. RNA-seq data indicated that PsMAD1 may regulate many zoospore development and infection associated genes. Thus, PsMAD1 may be a major regulator of P. sojae involved in zoosporogenesis and pathogenesis.

Metabolic Diversity and Novelties in the Oomycetes

The eukaryotic microbes called oomycetes include many important saprophytes and pathogens, with the latter exhibiting necrotrophy, biotrophy, or obligate biotrophy. Understanding oomycete metabolism is fundamental to understanding these lifestyles. Genome mining and biochemical studies have shown that oomycetes, which belong to the kingdom Stramenopila, secrete suites of carbohydrate- and protein-degrading enzymes adapted to their environmental niches and produce unusual lipids and energy storage compounds. Despite having limited secondary metabolism, many oomycetes make chemicals for communicating within their species or with their hosts. Horizontal and endosymbiotic gene transfer events have diversified oomycete metabolism, resulting in biochemical pathways that often depart from standard textbook descriptions by amalgamating enzymes from multiple sources. Gene fusions and duplications have further shaped the composition and expression of the enzymes. Current research is helping us learn how oomycetes interact with host and environment, understand eukaryotic diversity and evolution, and identify targets for drugs and crop protection chemicals.

Gene Expression and Silencing Studies in Phytophthora infestans Reveal Infection-Specific Nutrient Transporters and a Role for the Nitrate Reductase Pathway in Plant Pathogenesis

To help learn how phytopathogens feed from their hosts, genes for nutrient transporters from the hemibiotrophic potato and tomato pest Phytophthora infestans were annotated. This identified 453 genes from 19 families. Comparisons with a necrotrophic oomycete, Pythium ultimum var. ultimum, and a hemibiotrophic fungus, Magnaporthe oryzae, revealed diversity in the size of some families although a similar fraction of genes encoded transporters. RNA-seq of infected potato tubers, tomato leaves, and several artificial media revealed that 56 and 207 transporters from P. infestans were significantly up- or down-regulated, respectively, during early infection timepoints of leaves or tubers versus media. About 17 were up-regulated >4-fold in both leaves and tubers compared to media and expressed primarily in the biotrophic stage. The transcription pattern of many genes was host-organ specific. For example, the mRNA level of a nitrate transporter (NRT) was about 100-fold higher during mid-infection in leaves, which are nitrate-rich, than in tubers and three types of artificial media, which are nitrate-poor. The NRT gene is physically linked with genes encoding nitrate reductase (NR) and nitrite reductase (NiR), which mobilize nitrate into ammonium and amino acids. All three genes were coregulated. For example, the three genes were expressed primarily at mid-stage infection timepoints in both potato and tomato leaves, but showed little expression in potato tubers. Transformants down-regulated for all three genes were generated by DNA-directed RNAi, with silencing spreading from the NR target to the flanking NRT and NiR genes. The silenced strains were nonpathogenic on leaves but colonized tubers. We propose that the nitrate assimilation genes play roles both in obtaining nitrogen for amino acid biosynthesis and protecting P. infestans from natural or fertilization-induced nitrate and nitrite toxicity.

Little is known of how plant pathogens adapt to different growth conditions and host tissues. To understand the interaction between the filamentous eukaryotic microbe Phytophthora infestans and its potato and tomato hosts, we mined the genome for genes encoding proteins involved in nutrient uptake and measured their expression in leaves, tubers, and three artificial media. We observed dynamic changes between the growth conditions, and identified transporters expressed mainly in the biotrophic stage, leaves, tubers, or artificial media. When we blocked the expression of a nitrate transporter and two other genes involved in assimilating nitrate, we observed that those genes were required for successful colonization of nitrate-rich leaves but not nitrate-poor tissues, and that nitrate had become toxic to the silenced strains. We therefore hypothesize that the nitrate assimilation pathway may help the pathogen use inorganic nitrogen for nutrition and/or detoxify nitrate when its levels may become damaging.

Comparisons of Ribosomal Protein Gene Promoters Indicate Superiority of Heterologous Regulatory Sequences for Expressing Transgenes in Phytophthora infestans

Molecular genetics approaches in Phytophthora research can be hampered by the limited number of known constitutive promoters for expressing transgenes and the instability of transgene activity. We have therefore characterized genes encoding the cytoplasmic ribosomal proteins of Phytophthora and studied their suitability for expressing transgenes in P. infestans. Phytophthora spp. encode a standard complement of 79 cytoplasmic ribosomal proteins. Several genes are duplicated, and two appear to be pseudogenes. Half of the genes are expressed at similar levels during all stages of asexual development, and we discovered that the majority share a novel promoter motif named the PhRiboBox. This sequence is enriched in genes associated with transcription, translation, and DNA replication, including tRNA and rRNA biogenesis. Promoters from the three P. infestans genes encoding ribosomal proteins S9, L10, and L23 and their orthologs from P. capsici were tested for their ability to drive transgenes in stable transformants of P. infestans. Five of the six promoters yielded strong expression of a GUS reporter, but the stability of expression was higher using the P. capsici promoters. With the RPS9 and RPL10 promoters of P. infestans, about half of transformants stopped making GUS over two years of culture, while their P. capsici orthologs conferred stable expression. Since cross-talk between native and transgene loci may trigger gene silencing, we encourage the use of heterologous promoters in transformation studies.

Genome analyses of the sunflower pathogen Plasmopara halstedii provide insights into effector evolution in downy mildews and Phytophthora

BACKGROUND

Downy mildews are the most speciose group of oomycetes and affect crops of great economic importance. So far, there is only a single deeply-sequenced downy mildew genome available, from Hyaloperonospora arabidopsidis. Further genomic resources for downy mildews are required to study their evolution, including pathogenicity effector proteins, such as RxLR effectors. Plasmopara halstedii is a devastating pathogen of sunflower and a potential pathosystem model to study downy mildews, as several Avr-genes and R-genes have been predicted and unlike Arabidopsis downy mildew, large quantities of almost contamination-free material can be obtained easily.

RESULTS

Here a high-quality draft genome of Plasmopara halstedii is reported and analysed with respect to various aspects, including genome organisation, secondary metabolism, effector proteins and comparative genomics with other sequenced oomycetes. Interestingly, the present analyses revealed further variation of the RxLR motif, suggesting an important role of the conservation of the dEER-motif. Orthology analyses revealed the conservation of 28 RxLR-like core effectors among Phytophthora species. Only six putative RxLR-like effectors were shared by the two sequenced downy mildews, highlighting the fast and largely independent evolution of two of the three major downy mildew lineages. This is seemingly supported by phylogenomic results, in which downy mildews did not appear to be monophyletic.

CONCLUSIONS

The genome resource will be useful for developing markers for monitoring the pathogen population and might provide the basis for new approaches to fight Phytophthora and downy mildew pathogens by targeting core pathogenicity effectors.

Identification of putative TAL effector targets of the citrus canker pathogens shows functional convergence underlying disease development and defense response

Background

Transcriptional activator-like (TAL) effectors, formerly known as the AvrBs3/PthA protein family, are DNA-binding effectors broadly found in Xanthomonas spp. that transactivate host genes upon injection via the bacterial type three-secretion system. Biologically relevant targets of TAL effectors, i.e. host genes whose induction is vital to establish a compatible interaction, have been reported for xanthomonads that colonize rice and pepper; however, citrus genes modulated by the TAL effectors PthA“s” and PthC“s” of the citrus canker bacteria Xanthomonas citri (Xc) and Xanthomonas aurantifolii pathotype C (XaC), respectively, are poorly characterized. Of particular interest, XaC causes canker disease in its host lemon (Citrus aurantifolia), but triggers a defense response in sweet orange.

Results

Based on, 1) the TAL effector-DNA binding code, 2) gene expression data of Xc and XaC-infiltrated sweet orange leaves, and 3) citrus hypocotyls transformed with PthA2, PthA4 or PthC1, we have identified a collection of Citrus sinensis genes potentially targeted by Xc and XaC TAL effectors. Our results suggest that similar with other strains of Xanthomonas TAL effectors, PthA2 and PthA4, and PthC1 to some extent, functionally converge. In particular, towards induction of genes involved in the auxin and gibberellin synthesis and response, cell division, and defense response. We also present evidence indicating that the TAL effectors act as transcriptional repressors and that the best scoring predicted DNA targets of PthA“s” and PthC“s” in citrus promoters predominantly overlap with or localize near to TATA boxes of core promoters, supporting the idea that TAL effectors interact with the host basal transcriptional machinery to recruit the RNA pol II and start transcription.

Conclusions

The identification of PthA“s” and PthC“s” targets, such as the LOB (LATERAL ORGAN BOUNDARY) and CCNBS genes that we report here, is key for the understanding of the canker symptoms development during host susceptibility, or the defenses of sweet orange against the canker bacteria. We have narrowed down candidate targets to a few, which pointed out the host metabolic pathways explored by the pathogens.

Patterns of genomic variation in the poplar rust fungus Melampsora larici-populina identify pathogenesis-related factors

Melampsora larici-populina is a fungal pathogen responsible for foliar rust disease on poplar trees, which causes damage to forest plantations worldwide, particularly in Northern Europe. The reference genome of the isolate 98AG31 was previously sequenced using a whole genome shotgun strategy, revealing a large genome of 101 megabases containing 16,399 predicted genes, which included secreted protein genes representing poplar rust candidate effectors. In the present study, the genomes of 15 isolates collected over the past 20 years throughout the French territory, representing distinct virulence profiles, were characterized by massively parallel sequencing to assess genetic variation in the poplar rust fungus. Comparison to the reference genome revealed striking structural variations. Analysis of coverage and sequencing depth identified large missing regions between isolates related to the mating type loci. More than 611,824 single-nucleotide polymorphism (SNP) positions were uncovered overall, indicating a remarkable level of polymorphism. Based on the accumulation of non-synonymous substitutions in coding sequences and the relative frequencies of synonymous and non-synonymous polymorphisms (i.e., PN/PS ), we identify candidate genes that may be involved in fungal pathogenesis. Correlation between non-synonymous SNPs in genes encoding secreted proteins (SPs) and pathotypes of the studied isolates revealed candidate genes potentially related to virulences 1, 6, and 8 of the poplar rust fungus.

Phylogenetic and transcriptional analysis of an expanded bZIP transcription factor family in Phytophthora sojae

BACKGROUND

Basic leucine zipper (bZIP) transcription factors are present exclusively in eukaryotes and constitute one of the largest and most diverse transcription factor families. The proteins are responsible for central developmental and physiological processes in plants, animals, and fungi, including the pathogenicity of fungal plant pathogens. However, there is limited understanding of bZIPs in oomycetes, which are fungus-like organisms in the kingdom Stramenopila. Oomycetes include many destructive plant pathogens, including the well-studied species Phytophthora sojae, which causes soybean stem and root rot.

RESULTS

Candidate bZIPs encoded in the genomes of P. sojae and four other oomycetes, two diatoms, and two fungal species were predicted using bioinformatic methods. Comparative analysis revealed expanded numbers of bZIP candidates in oomycetes, especially the Phytophthora species, due to the expansion of several novel bZIP classes whose highly conserved asparagines in basic DNA-binding regions were substituted by other residues such as cysteine. The majority of these novel bZIP classes were mostly restricted to oomycetes. The large number of novel bZIPs appears to be the result of widespread gene duplications during oomycete evolution. The majority of P. sojae bZIP candidates, including both conventional and novel bZIP classes, were predicted to contain canonical protein secondary structures. Detection of gene transcripts using digital gene expression profiling and qRT-PCR suggested that most of the candidates were not pseudogenes. The major transcriptional shifts of bZIPs occurred during the zoosporangia/zoospore/cyst and host infection stages. Several infection-associated bZIP genes were identified that were positively regulated by H2O2 exposure.

CONCLUSIONS

The identification of large classes of bZIP proteins in oomycetes with novel bZIP motif variants, that are conserved and developmentally regulated and thus presumably functional, extends our knowledge of this important family of eukaryotic transcription factors. It also lays the foundation for detailed studies of the roles of these proteins in development and infection in P. sojae and other oomycetes.

Phenotypic diversification by gene silencing in Phytophthora plant pathogens

Advances in genome sequencing technologies have enabled generation of unprecedented information on genome content and organization. Eukaryote genomes in particular may contain large populations of transposable elements (TEs) and other repeated sequences. Active TEs can result in insertional mutations, altered transcription levels and ectopic recombination of DNA. The genome of the oomycete plant pathogen, Phytophthora infestans, contains vast numbers of TE sequences. There are also hundreds of predicted disease-promoting effector proteins, predominantly located in TE-rich genomic regions. Expansion of effector gene families is also a genomic signature of related oomycetes such as P. sojae. Deep sequencing of small RNAs (sRNAs) from P. infestans has identified sRNAs derived from all families of transposons, highlighting the importance of RNA silencing for maintaining these genomic invaders in an inactive form. Small RNAs were also identified from specific effector encoding genes, possibly leading to RNA silencing of these genes and variation in pathogenicity and virulence toward plant resistance genes. Similar findings have also recently been made for the distantly related species, P. sojae. Small RNA “hotspots” originating from arrays of amplified gene sequences, or from genes displaying overlapping antisense transcription, were also identified in P. infestans. These findings suggest a major role for RNA silencing processes in the adaptability and diversification of these economically important plant pathogens. Here we review the latest progress and understanding of gene silencing in oomycetes with emphasis on transposable elements and sRNA-associated events.

From pathogen genomes to host plant processes: the power of plant parasitic oomycetes

Recent pathogenomic research on plant parasitic oomycete effector function and plant host responses has resulted in major conceptual advances in plant pathology, which has been possible thanks to the availability of genome sequences.

Novel core promoter elements in the oomycete pathogen Phytophthora infestans and their influence on expression detected by genome-wide analysis

Background

The core promoter is the region flanking the transcription start site (TSS) that directs formation of the pre-initiation complex. Core promoters have been studied intensively in mammals and yeast, but not in more diverse eukaryotes. Here we investigate core promoters in oomycetes, a group within the Stramenopile kingdom that includes important plant and animal pathogens. Prior studies of a small collection of genes proposed that oomycete core promoters contain a 16 to 19 nt motif bearing an Initiator-like sequence (INR) flanked by a novel sequence named FPR, but this has not been extended to whole-genome analysis.

Results

We used expectation maximization to find over-represented motifs near TSSs of Phytophthora infestans, the potato blight pathogen. The motifs corresponded to INR, FPR, and a new element found about 25 nt downstream of the TSS called DPEP. TATA boxes were not detected. Assays of DPEP function by mutagenesis were consistent with its role as a core motif. Genome-wide searches found a well-conserved combined INR+FPR in only about 13% of genes after correcting for false discovery, which contradicted prior reports that INR and FPR are found together in most genes. INR or FPR were found alone near TSSs in 18% and 7% of genes, respectively. Promoters lacking the motifs had pyrimidine-rich regions near the TSS. The combined INR+FPR motif was linked to higher than average mRNA levels, developmentally-regulated transcription, and functions related to plant infection, while DPEP and FPR were over-represented in constitutively-expressed genes. The INR, FPR, and combined INR+FPR motifs were detected in other oomycetes including Hyaloperonospora arabidopsidis, Phytophthora sojae, Pythium ultimum, and Saprolegnia parasitica, while DPEP was found in all but S. parasitica. Only INR seemed present in a non-oomycete stramenopile.

Conclusions

The absence of a TATA box and presence of novel motifs show that the oomycete core promoter is diverged from that of model systems, and likely explains the lack of activity of non-oomycete promoters in Phytophthora transformants. The association of the INR+FPR motif with developmentally-regulated genes shows that oomycete core elements influence stage-specific transcription in addition to regulating formation of the pre-initiation complex.