Arabidopsis Topless-related 1 mitigates physiological damage and growth penalties of induced immunity

Transcriptional corepressors of the Topless (TPL) family regulate plant hormone and immunity signaling. The lack of a genome-wide profile of their chromatin associations limits understanding of the TPL family roles in transcriptional regulation. Chromatin immunoprecipitation with sequencing (ChIP-Seq) was performed on Arabidopsis thaliana lines expressing GFP-tagged Topless-related 1 (TPR1-GFP) with and without constitutive immunity via Enhanced Disease Susceptibility 1 (EDS1). RNA-Seq profiling of the TPR1-GFP lines and pathogen-infected tpl/tpr mutants, combined with measuring immunity, growth, and physiological parameters was employed to investigate TPL/TPR roles in immunity and defense homeostasis. TPR1 was enriched at promoter regions of c. 1400 genes and c. 10% of the detected binding required EDS1 immunity signaling. In a tpr1 tpl tpr4 (t3) mutant, resistance to bacteria was slightly compromised, and defense-related transcriptional reprogramming was weakly reduced or enhanced, respectively, at early (< 1 h) and late 24 h stages of bacterial infection. The t3 plants challenged with bacteria or pathogen-associated molecular pattern nlp24 displayed photosystem II dysfunctions. Also, t3 plants were hypersensitive to phytocytokine pep1 at the level of root growth inhibition. Transgenic expression of TPR1 rescued these t3 physiological defects. We propose that TPR1 and TPL family proteins function in Arabidopsis to reduce detrimental effects associated with activated transcriptional immunity.

A rodent model for Dirofilaria immitis, canine heartworm: parasite growth, development, and drug sensitivity in NSG mice

Heartworm disease, caused by Dirofilaria immitis, remains a significant threat to canines and felines. The development of parasites resistant to macrocyclic lactones (ML) has created a significant challenge to the control of the infection. The goal of this study was to determine if mice lacking a functional immune response would be susceptible to D. immitis. Immunodeficient NSG mice were susceptible to the infection, sustaining parasites for at least 15 weeks, with infective third-stage larvae molting and developing into the late fourth-stage larvae. Proteomic analysis of host responses to the infection revealed a complex pattern of changes after infection, with at least some of the responses directed at reducing immune control mechanisms that remain in NSG mice. NSG mice were infected with isolates of D. immitis that were either susceptible or resistant to MLs, as a population. The susceptible isolate was killed by ivermectin whereas the resistant isolate had improved survivability, while both isolates were affected by moxidectin. It was concluded that D. immitis survives in NSG mice for at least 15 weeks. NSG mice provide an ideal model for monitoring host responses to the infection and for testing parasites in vivo for susceptibility to direct chemotherapeutic activity of new agents.

The leucine-rich repeats in allelic barley MLA immune receptors define specificity towards sequence-unrelated powdery mildew avirulence effectors with a predicted common RNase-like fold

Nucleotide-binding domain leucine-rich repeat-containing receptors (NLRs) in plants can detect avirulence (AVR) effectors of pathogenic microbes. The Mildew locus a (Mla) NLR gene has been shown to confer resistance against diverse fungal pathogens in cereal crops. In barley, Mla has undergone allelic diversification in the host population and confers isolate-specific immunity against the powdery mildew-causing fungal pathogen Blumeria graminis forma specialis hordei (Bgh). We previously isolated the Bgh effectors AVRA1, AVRA7, AVRA9, AVRA13, and allelic AVRA10/AVRA22, which are recognized by matching MLA1, MLA7, MLA9, MLA13, MLA10 and MLA22, respectively. Here, we extend our knowledge of the Bgh effector repertoire by isolating the AVRA6 effector, which belongs to the family of catalytically inactive RNase-Like Proteins expressed in Haustoria (RALPHs). Using structural prediction, we also identified RNase-like folds in AVRA1, AVRA7, AVRA10/AVRA22, and AVRA13, suggesting that allelic MLA recognition specificities could detect structurally related avirulence effectors. To better understand the mechanism underlying the recognition of effectors by MLAs, we deployed chimeric MLA1 and MLA6, as well as chimeric MLA10 and MLA22 receptors in plant co-expression assays, which showed that the recognition specificity for AVRA1 and AVRA6 as well as allelic AVRA10 and AVRA22 is largely determined by the receptors' C-terminal leucine-rich repeats (LRRs). The design of avirulence effector hybrids allowed us to identify four specific AVRA10 and five specific AVRA22 aa residues that are necessary to confer MLA10- and MLA22-specific recognition, respectively. This suggests that the MLA LRR mediates isolate-specific recognition of structurally related AVRA effectors. Thus, functional diversification of multi-allelic MLA receptors may be driven by a common structural effector scaffold, which could be facilitated by proliferation of the RALPH effector family in the pathogen genome.

Multiple pairs of allelic MLA immune receptor-powdery mildew AVRA effectors argue for a direct recognition mechanism

Nucleotide-binding domain and leucine-rich repeat (NLR)-containing proteins in plants and animals mediate intracellular pathogen sensing. Plant NLRs typically detect strain-specific pathogen effectors and trigger immune responses often linked to localized host cell death. The barley Mla disease resistance locus has undergone extensive functional diversification in the host population and encodes numerous allelic NLRs each detecting a matching isolate-specific avirulence effector (AVRA) of the fungal pathogen Blumeria graminis f. sp. hordei (Bgh). We report here the isolation of Bgh AVRa7, AVRa9, AVRa10, and AVRa22, which encode small secreted proteins recognized by allelic MLA7, MLA9, MLA10, and MLA22 receptors, respectively. These effectors are sequence-unrelated, except for allelic AVRa10 and AVRa22 that are co-maintained in pathogen populations in the form of a balanced polymorphism. Contrary to numerous examples of indirect recognition of bacterial effectors by plant NLRs, co-expression experiments with matching Mla-AVRa pairs indicate direct detection of the sequence-unrelated fungal effectors by MLA receptors.

Representing dynamic biological networks with multi-scale probabilistic models

Dynamic models analyzing gene regulation and metabolism face challenges when adapted to modeling signal transduction networks. During signal transduction, molecular reactions and mechanisms occur in different spatial and temporal frames and involve feedbacks. This impedes the straight-forward use of methods based on Boolean networks, Bayesian approaches, and differential equations. We propose a new approach, ProbRules, that combines probabilities and logical rules to represent the dynamics of a system across multiple scales. We demonstrate that ProbRules models can represent various network motifs of biological systems. As an example of a comprehensive model of signal transduction, we provide a Wnt network that shows remarkable robustness under a range of phenotypical and pathological conditions. Its simulation allows the clarification of controversially discussed molecular mechanisms of Wnt signaling by predicting wet-lab measurements. ProbRules provides an avenue in current computational modeling by enabling systems biologists to integrate vast amounts of available data on different scales.

Subfamily-Specific Specialization of RGH1/MLA Immune Receptors in Wild Barley

The barley disease resistance (R) gene locus mildew locus A (Mla) provides isolate-specific resistance against the powdery mildew fungus Blumeria graminis hordei and has been introgressed into modern cultivars from diverse germplasms, including the wild relative Hordeum spontaneum. Known Mla disease resistance specificities to B. graminis hordei appear to encode allelic variants of the R gene homolog 1 (RGH1) family of nucleotide-binding domain and leucine-rich repeat (NLR) proteins. Here, we sequenced and assembled the transcriptomes of 50 H. spontaneum accessions representing nine populations distributed throughout the Fertile Crescent. The assembled Mla transcripts exhibited rich sequence diversity, linked neither to geographic origin nor population structure, and could be grouped into two similar-sized subfamilies based on two major N-terminal coiled-coil (CC) signaling domains that are both capable of eliciting cell death. The presence of positively selected sites located mainly in the C-terminal leucine-rich repeats of both MLA subfamilies, together with the fact that both CC signaling domains mediate cell death, implies that the two subfamilies are actively maintained in the population. Unexpectedly, known MLA receptor variants that confer B. graminis hordei resistance belong exclusively to one subfamily. Thus, signaling domain divergence, potentially as adaptation to distinct pathogen populations, is an evolutionary signature of functional diversification of an immune receptor. Copyright © 2018 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Principles and characteristics of the Arabidopsis WRKY regulatory network during early MAMP-triggered immunity

During microbe-associated molecular pattern-triggered immunity more than 5000 Arabidopsis genes are significantly altered in their expression, and the question arises, how such an enormous reprogramming of the transcriptome can be regulated in a safe and robust manner? For the WRKY transcription factors (TFs), which are important regulators of numerous defense responses, it appears that they act in a complex regulatory sub-network rather than in a linear fashion, which would be much more vulnerable to gene function loss either by pathogen-derived effectors or by mutations. In this study we employed RNA-seq, mass spectrometry and chromatin immunoprecipitation-seq to find evidence for and uncover principles and characteristics of this network. Upon flg22-treatment, one can distinguish between two sets of WRKY genes: constitutively expressed and induced WRKY genes. Prior to elicitation the induced WRKY genes appear to be maintained in a repressed state mainly by the constitutively expressed WRKY factors, which themselves appear to be regulated by non-WRKY TFs. Upon elicitation, induced WRKYs rapidly bind to induced WRKY gene promoters and by auto- and cross-regulation build up the regulatory network. Maintenance of this flg22-induced network appears highly robust as removal of three key WRKY factors can be physically and functionally compensated for by other WRKY family members.

The Defense Phytohormone Signaling Network Enables Rapid, High-Amplitude Transcriptional Reprogramming during Effector-Triggered Immunity

The phytohormone network consisting of jasmonate, ethylene, PHYTOALEXIN-DEFICIENT4, and salicylic acid signaling is required for the two modes of plant immunity, pattern-triggered immunity (PTI), and effector-triggered immunity (ETI). A previous study showed that during PTI, the transcriptional responses of over 5000 genes qualitatively depend on complex interactions between the network components. However, the role of the network in transcriptional reprogramming during ETI and whether it differs between PTI and ETI remain elusive. Here, we generated time-series RNA-sequencing data of Arabidopsis thaliana wild-type and combinatorial mutant plants deficient in components of the network upon challenge with virulent or ETI-triggering avirulent strains of the foliar bacterial pathogen Pseudomonas syringae Resistant plants such as the wild type achieved high-amplitude transcriptional reprogramming 4 h after challenge with avirulent strains and sustained this transcriptome response. Strikingly, susceptible plants including the quadruple network mutant showed almost identical transcriptome responses to resistant plants but with several hours delay. Furthermore, gene coexpression network structure was highly conserved between the wild type and quadruple mutant. Thus, in contrast to PTI, the phytohormone network is required only for achieving high-amplitude transcriptional reprogramming within the early time window of ETI against this bacterial pathogen.

Signatures of host specialization and a recent transposable element burst in the dynamic one-speed genome of the fungal barley powdery mildew pathogen

Background

Powdery mildews are biotrophic pathogenic fungi infecting a number of economically important plants. The grass powdery mildew, Blumeria graminis, has become a model organism to study host specialization of obligate biotrophic fungal pathogens. We resolved the large-scale genomic architecture of B. graminis forma specialis hordei (Bgh) to explore the potential influence of its genome organization on the co-evolutionary process with its host plant, barley (Hordeum vulgare).

Results

The near-chromosome level assemblies of the Bgh reference isolate DH14 and one of the most diversified isolates, RACE1, enabled a comparative analysis of these haploid genomes, which are highly enriched with transposable elements (TEs). We found largely retained genome synteny and gene repertoires, yet detected copy number variation (CNV) of secretion signal peptide-containing protein-coding genes (SPs) and locally disrupted synteny blocks. Genes coding for sequence-related SPs are often locally clustered, but neither the SPs nor the TEs reside preferentially in genomic regions with unique features. Extended comparative analysis with different host-specific B. graminis formae speciales revealed the existence of a core suite of SPs, but also isolate-specific SP sets as well as congruence of SP CNV and phylogenetic relationship. We further detected evidence for a recent, lineage-specific expansion of TEs in the Bgh genome.

Conclusions

The characteristics of the Bgh genome (largely retained synteny, CNV of SP genes, recently proliferated TEs and a lack of significant compartmentalization) are consistent with a “one-speed” genome that differs in its architecture and (co-)evolutionary pattern from the “two-speed” genomes reported for several other filamentous phytopathogens.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4750-6) contains supplementary material, which is available to authorized users.

A dominant-interfering camta3 mutation compromises primary transcriptional outputs mediated by both cell surface and intracellular immune receptors in Arabidopsis thaliana

Pattern recognition receptors (PRRs) and nucleotide-binding domain and leucine-rich repeat (LRR)-containing proteins (NLRs) initiate pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), respectively, each associated with the activation of an overlapping set of defence genes. The regulatory mechanism behind this convergence of PTI- and ETI-mediated defence gene induction remains elusive. We generated transgenic Arabidopsis plants that enable conditional NLR activation without pathogen infection to dissect NLR- and PRR-mediated transcriptional signals. A comparative analysis of over 40 transcriptome datasets linked calmodulin-binding transcription activators (CAMTAs) to the activation of overlapping defence genes in PTI and ETI. We used a dominant camta3 mutant (camta3-D) to assess CAMTA functions in the corresponding transcriptional regulation. Transcriptional regulation by NLRs, although highly similar to PTI responses, can be established independently of pathogen-associated molecular pattern (PAMP) perception, defence phytohormones and host cell death. Conditional expression of the N-terminal coiled-coil domain of the barley MLA (Mildew resistance locus A) NLR is sufficient to trigger similar transcriptional reprogramming as full-length NLRs. CAMTA-binding motifs are overrepresented in the 5' regulatory regions of the identified primary immune response genes, consistent with their altered expression and disease resistance responses in camta3-D plants. We propose that CAMTA-mediated transcriptional regulation defines an early convergence point in NLR- and PRR-mediated signalling.

The Root Growth-Regulating Brevicompanine Natural Products Modulate the Plant Circadian Clock

Plant growth regulating properties of brevicompanines (Brvs), natural products of the fungus Penicillium brevicompactum, have been known for several years, but further investigations into the molecular mechanism of their bioactivity have not been performed. Following chemical synthesis of brevicompanine derivatives, we studied their activity in the model plant Arabidopsis by a combination of plant growth assays, transcriptional profiling, and numerous additional bioassays. These studies demonstrated that brevicompanines cause transcriptional misregulation of core components of the circadian clock, whereas other biological read-outs were not affected. Brevicompanines thus represent promising chemical tools for investigating the regulation of the plant circadian clock. In addition, our study also illustrates the potential of an unbiased -omics-based characterization of bioactive compounds for identifying the often cryptic modes of action of small molecules.

A core function of EDS1 with PAD4 is to protect the salicylic acid defense sector in Arabidopsis immunity

Plant defenses induced by salicylic acid (SA) are vital for resistance against biotrophic pathogens. In basal and receptor-triggered immunity, SA accumulation is promoted by Enhanced Disease Susceptibility1 with its co-regulator Phytoalexin Deficient4 (EDS1/PAD4). Current models position EDS1/PAD4 upstream of SA but their functional relationship remains unclear. In a genetic and transcriptomic analysis of Arabidopsis autoimmunity caused by constitutive or conditional EDS1/PAD4 overexpression, intrinsic EDS1/PAD4 signaling properties and their relation to SA were uncovered. A core EDS1/PAD4 pathway works in parallel with SA in basal and effector-triggered bacterial immunity. It protects against disabled SA-regulated gene expression and pathogen resistance, and is distinct from a known SA-compensatory route involving MAPK signaling. Results help to explain previously identified EDS1/PAD4 regulated SA-dependent and SA-independent gene expression sectors. Plants have evolved an alternative route for preserving SA-regulated defenses against pathogen or genetic perturbations. In a proposed signaling framework, EDS1 with PAD4, besides promoting SA biosynthesis, maintains important SA-related resistance programs, thereby increasing robustness of the innate immune system.

Induced Genome-Wide Binding of Three Arabidopsis WRKY Transcription Factors during Early MAMP-Triggered Immunity

During microbial-associated molecular pattern-triggered immunity (MTI), molecules derived from microbes are perceived by cell surface receptors and upon signaling to the nucleus initiate a massive transcriptional reprogramming critical to mount an appropriate host defense response. WRKY transcription factors play an important role in regulating these transcriptional processes. Here, we determined on a genome-wide scale the flg22-induced in vivo DNA binding dynamics of three of the most prominent WRKY factors, WRKY18, WRKY40, and WRKY33. The three WRKY factors each bound to more than 1000 gene loci predominantly at W-box elements, the known WRKY binding motif. Binding occurred mainly in the 500-bp promoter regions of these genes. Many of the targeted genes are involved in signal perception and transduction not only during MTI but also upon damage-associated molecular pattern-triggered immunity, providing a mechanistic link between these functionally interconnected basal defense pathways. Among the additional targets were genes involved in the production of indolic secondary metabolites and in modulating distinct plant hormone pathways. Importantly, among the targeted genes were numerous transcription factors, encoding predominantly ethylene response factors, active during early MTI, and WRKY factors, supporting the previously hypothesized existence of a WRKY subregulatory network. Transcriptional analysis revealed that WRKY18 and WRKY40 function redundantly as negative regulators of flg22-induced genes often to prevent exaggerated defense responses.

Survival trade-offs in plant roots during colonization by closely related beneficial and pathogenic fungi

The sessile nature of plants forced them to evolve mechanisms to prioritize their responses to simultaneous stresses, including colonization by microbes or nutrient starvation. Here, we compare the genomes of a beneficial root endophyte, Colletotrichum tofieldiae and its pathogenic relative C. incanum, and examine the transcriptomes of both fungi and their plant host Arabidopsis during phosphate starvation. Although the two species diverged only 8.8 million years ago and have similar gene arsenals, we identify genomic signatures indicative of an evolutionary transition from pathogenic to beneficial lifestyles, including a narrowed repertoire of secreted effector proteins, expanded families of chitin-binding and secondary metabolism-related proteins, and limited activation of pathogenicity-related genes in planta. We show that beneficial responses are prioritized in C. tofieldiae-colonized roots under phosphate-deficient conditions, whereas defense responses are activated under phosphate-sufficient conditions. These immune responses are retained in phosphate-starved roots colonized by pathogenic C. incanum, illustrating the ability of plants to maximize survival in response to conflicting stresses.

Root Endophyte Colletotrichum tofieldiae Confers Plant Fitness Benefits that Are Phosphate Status Dependent

A staggering diversity of endophytic fungi associate with healthy plants in nature, but it is usually unclear whether these represent stochastic encounters or provide host fitness benefits. Although most characterized species of the fungal genus Colletotrichum are destructive pathogens, we show here that C. tofieldiae (Ct) is an endemic endophyte in natural Arabidopsis thaliana populations in central Spain. Colonization by Ct initiates in roots but can also spread systemically into shoots. Ct transfers the macronutrient phosphorus to shoots, promotes plant growth, and increases fertility only under phosphorus-deficient conditions, a nutrient status that might have facilitated the transition from pathogenic to beneficial lifestyles. The host's phosphate starvation response (PSR) system controls Ct root colonization and is needed for plant growth promotion (PGP). PGP also requires PEN2-dependent indole glucosinolate metabolism, a component of innate immune responses, indicating a functional link between innate immunity and the PSR system during beneficial interactions with Ct.

Root Endophyte Colletotrichum tofieldiae Confers Plant Fitness Benefits that Are Phosphate Status Dependent

A staggering diversity of endophytic fungi associate with healthy plants in nature, but it is usually unclear whether these represent stochastic encounters or provide host fitness benefits. Although most characterized species of the fungal genus Colletotrichum are destructive pathogens, we show here that C. tofieldiae (Ct) is an endemic endophyte in natural Arabidopsis thaliana populations in central Spain. Colonization by Ct initiates in roots but can also spread systemically into shoots. Ct transfers the macronutrient phosphorus to shoots, promotes plant growth, and increases fertility only under phosphorus-deficient conditions, a nutrient status that might have facilitated the transition from pathogenic to beneficial lifestyles. The host’s phosphate starvation response (PSR) system controls Ct root colonization and is needed for plant growth promotion (PGP). PGP also requires PEN2-dependent indole glucosinolate metabolism, a component of innate immune responses, indicating a functional link between innate immunity and the PSR system during beneficial interactions with Ct.

•

Colletotrichum tofieldiae (Ct) is a fungal root endophyte of Arabidopsis

•

Ct transfers the macronutrient phosphorus to Arabidopsis shoots

•

Ct-mediated plant growth promotion needs an intact phosphate starvation response

•

A branch of the plant innate immune system is essential for beneficial Ct activities

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified &gt;1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified &gt;1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified &gt;1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Mosaic genome structure of the barley powdery mildew pathogen and conservation of transcriptional programs in divergent hosts

Barley powdery mildew, Blumeria graminis f. sp. hordei (Bgh), is an obligate biotrophic ascomycete fungal pathogen that can grow and reproduce only on living cells of wild or domesticated barley (Hordeum sp.). Domestication and deployment of resistant barley cultivars by humans selected for amplification of Bgh isolates with different virulence combinations. We sequenced the genomes of two European Bgh isolates, A6 and K1, for comparative analysis with the reference genome of isolate DH14. This revealed a mosaic genome structure consisting of large isolate-specific DNA blocks with either high or low SNP densities. Some of the highly polymorphic blocks likely accumulated SNPs for over 10,000 years, well before the domestication of barley. These isolate-specific blocks of alternating monomorphic and polymorphic regions imply an exceptionally large standing genetic variation in the Bgh population and might be generated and maintained by rare outbreeding and frequent clonal reproduction. RNA-sequencing experiments with isolates A6 and K1 during four early stages of compatible and incompatible interactions on leaves of partially immunocompromised Arabidopsis mutants revealed a conserved Bgh transcriptional program during pathogenesis compared with the natural host barley despite ~200 million years of reproductive isolation of these hosts. Transcripts encoding candidate-secreted effector proteins are massively induced in successive waves. A specific decrease in candidate-secreted effector protein transcript abundance in the incompatible interaction follows extensive transcriptional reprogramming of the host transcriptome and coincides with the onset of localized host cell death, suggesting a host-inducible defense mechanism that targets fungal effector secretion or production.

Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses

Colletotrichum species are fungal pathogens that devastate crop plants worldwide. Host infection involves the differentiation of specialized cell types that are associated with penetration, growth inside living host cells (biotrophy) and tissue destruction (necrotrophy). We report here genome and transcriptome analyses of Colletotrichum higginsianum infecting Arabidopsis thaliana and Colletotrichum graminicola infecting maize. Comparative genomics showed that both fungi have large sets of pathogenicity-related genes, but families of genes encoding secreted effectors, pectin-degrading enzymes, secondary metabolism enzymes, transporters and peptidases are expanded in C. higginsianum. Genome-wide expression profiling revealed that these genes are transcribed in successive waves that are linked to pathogenic transitions: effectors and secondary metabolism enzymes are induced before penetration and during biotrophy, whereas most hydrolases and transporters are upregulated later, at the switch to necrotrophy. Our findings show that preinvasion perception of plant-derived signals substantially reprograms fungal gene expression and indicate previously unknown functions for particular fungal cell types.

Inferring Boolean network structure via correlation

MOTIVATION

Accurate, context-specific regulation of gene expression is essential for all organisms. Accordingly, it is very important to understand the complex relations within cellular gene regulatory networks. A tool to describe and analyze the behavior of such networks are Boolean models. The reconstruction of a Boolean network from biological data requires identification of dependencies within the network. This task becomes increasingly computationally demanding with large amounts of data created by recent high-throughput technologies. Thus, we developed a method that is especially suited for network structure reconstruction from large-scale data. In our approach, we took advantage of the fact that a specific transcription factor often will consistently either activate or inhibit a specific target gene, and this kind of regulatory behavior can be modeled using monotone functions.

RESULTS

To detect regulatory dependencies in a network, we examined how the expression of different genes correlates to successive network states. For this purpose, we used Pearson correlation as an elementary correlation measure. Given a Boolean network containing only monotone Boolean functions, we prove that the correlation of successive states can identify the dependencies in the network. This method not only finds dependencies in randomly created artificial networks to very high percentage, but also reconstructed large fractions of both a published Escherichia coli regulatory network from simulated data and a yeast cell cycle network from real microarray data.