Search

The selective pressure of pathogen-host symbiosis drives adaptations. How these interactions shape the metabolism of pathogens is largely unknown. Here, we use comparative genomics to systematically analyze the metabolic networks of oomycetes, a diverse group of eukaryotes that includes saprotrophs as well as animal and plant pathogens, with the latter causing devastating diseases with significant economic and/or ecological impacts. In our analyses of 44 oomycete species, we uncover considerable variation in metabolism that can be linked to lifestyle differences. Comparisons of metabolic gene content reveal that plant pathogenic oomycetes have a bipartite metabolism consisting of a conserved core and an accessory set. The accessory set can be associated with the degradation of defense compounds produced by plants when challenged by pathogens. Obligate biotrophic oomycetes have smaller metabolic networks, and taxonomically distantly related biotrophic lineages display convergent evolution by repeated gene losses in both the conserved as well as the accessory set of metabolisms. When investigating to what extent the metabolic networks in obligate biotrophs differ from those in hemibiotrophic plant pathogens, we observe that the losses of metabolic enzymes in obligate biotrophs are not random and that gene losses predominantly influence the terminal branches of the metabolic networks. Our analyses represent the first metabolism-focused comparison of oomycetes at this scale and will contribute to a better understanding of the evolution of oomycete metabolism in relation to lifestyle adaptation. Numerous oomycete species are devastating plant pathogens that cause major damage in crops and natural ecosystems. Their interactions with hosts are shaped by strong selection, but how selection affects adaptation of the primary metabolism to a pathogenic lifestyle is not yet well established. By pan-genome and metabolic network analyses of distantly related oomycete pathogens and their nonpathogenic relatives, we reveal considerable lifestyle- and lineage-specific adaptations. This study contributes to a better understanding of metabolic adaptations in pathogenic oomycetes in relation to lifestyle, host, and environment, and the findings will help in pinpointing potential targets for disease control. [Graphic: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Abstract Cell-surface receptors form the front line of plant immunity. The leucine-rich repeat (LRR)-receptor-like kinases SOBIR1 and BAK1 are required for the functionality of the tomato LRR-receptor-like protein Cf-4, which detects the secreted effector Avr4 of the pathogenic fungus Fulvia fulva. Here, we show that the kinase domains of SOBIR1 and BAK1 directly phosphorylate each other and that residues Thr522 and Tyr469 of the kinase domain of Nicotiana benthamiana SOBIR1 are required for its kinase activity and for interacting with signalling partners, respectively. By knocking out multiple genes belonging to different receptor-like cytoplasmic kinase (RLCK)-VII subfamilies in N. benthamiana:Cf-4, we show that members of RLCK-VII-6, −7, and −8 differentially regulate the Avr4/Cf-4-triggered biphasic burst of reactive oxygen species. In addition, members of RLCK-VII-7 play an essential role in resistance against the oomycete pathogen Phytophthora palmivora. Our study provides molecular evidence for the specific roles of RLCKs downstream of SOBIR1/BAK1-containing immune complexes.

The recent incursion of Fusarium wilt of banana (FWB) caused by TR4 threatens banana production in the Chira Valley of Peru. To develop a management strategy, we mapped the entire production area from the air. During 12 flights at an altitude of 610 m (2,000 ft), we gathered 133,700 images in a timeframe of 2 weeks and constructed an orthomosaic map of 73,000 ha. This unveiled the complex logistic network across the banana-producing region, comprising 150 km of primary roads, as well as numerous secondary and tertiary unpaved paths through and bordering banana plantations. Moreover, the Poechos reservoir—with a total length of 450 km—feeds the entire Chira Valley irrigation system, which could further exacerbate TR4 expansion. Using georeferenced landscape details (digital terrain model), we determined areas prone to flooding, which is also relevant for disseminating TR4. Analyses of the maps resulted in the identification of many suspect areas for direct sampling and subsequent analyses to study the expansion of TR4. At four locations, we confirmed two cases, and by June 8, 2023, a total of 207 cases were reported. We conclude that TR4 is a serious threat to the entire region, which exports approximately 25% of the global organic bananas. [Figure: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Abstract The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Abstract Polycomb group proteins, as part of the Polycomb repressive complexes, are essential in gene repression through chromatin compaction by canonical PRC1, mono-ubiquitylation of histone H2A by non-canonical PRC1 and tri-methylation of histone H3K27 by PRC2. Despite prevalent models emphasizing tight functional coupling between PRC1 and PRC2, it remains unclear whether this paradigm indeed reflects the evolution and functioning of these complexes. Here, we conduct a comprehensive analysis of the presence or absence of cPRC1, nPRC1 and PRC2 across the entire eukaryotic tree of life, and find that both complexes were present in the Last Eukaryotic Common Ancestor (LECA). Strikingly, ~42% of organisms contain only PRC1 or PRC2, showing that their evolution since LECA is largely uncoupled. The identification of ncPRC1-defining subunits in unicellular relatives of animals and fungi suggests ncPRC1 originated before cPRC1, and we propose a scenario for the evolution of cPRC1 from ncPRC1. Together, our results suggest that crosstalk between these complexes is a secondary development in evolution.

Abstract The Arabidopsis thaliana Receptor-Like Protein RLP30 contributes to immunity against the fungal pathogen Sclerotinia sclerotiorum. Here we identify the RLP30-ligand as a small cysteine-rich protein (SCP) that occurs in many fungi and oomycetes and is also recognized by the Nicotiana benthamiana RLP RE02. However, RLP30 and RE02 share little sequence similarity and respond to different parts of the native/folded protein. Moreover, some Brassicaceae other than Arabidopsis also respond to a linear SCP peptide instead of the folded protein, suggesting that SCP is an eminent immune target that led to the convergent evolution of distinct immune receptors in plants. Surprisingly, RLP30 shows a second ligand specificity for a SCP-nonhomologous protein secreted by bacterial Pseudomonads. RLP30 expression in N. tabacum results in quantitatively lower susceptibility to bacterial, fungal and oomycete pathogens, thus demonstrating that detection of immunogenic patterns by Arabidopsis RLP30 is involved in defense against pathogens from three microbial kingdoms.

Zymoseptoria tritici is an important fungal pathogen of wheat which has spread globally. Here, the authors perform genomic analyses on a collection of ~1100 Z. tritici samples from 42 countries to describe its global spread and elucidate mechanisms of adaptation to different environmental conditions.

Through the association of protein complexes to DNA, the eukaryotic nuclear genome is broadly organized into open euchromatin that is accessible for enzymes acting on DNA and condensed heterochromatin that is inaccessible. Chemical and physical alterations to chromatin may impact its organization and functionality and are therefore important regulators of nuclear processes. Studies in various fungal plant pathogens have uncovered an association between chromatin organization and expression of in planta-induced genes that are important for pathogenicity. This review discusses chromatin-based regulation mechanisms as determined in the fungal plant pathogen Verticillium dahliae and relates the importance of epigenetic transcriptional regulation and other nuclear processes more broadly in fungal plant pathogens.

Summary: Polycomb Repressive Complex 2 (PRC2) catalyzes the mono-, di-, and trimethylation of histone protein H3 on lysine 27 (H3K27), which is strongly associated with transcriptionally silent chromatin. The functional core of PRC2 is highly conserved in animals and consists of four subunits. One of these, SUZ12, has not been identified in the genetic model Caenorhabditis elegans, whereas C. elegans PRC2 contains the clade-specific MES-3 protein. Through unbiased sensitive sequence similarity searches complemented by high-quality structure predictions of monomers and multimers, we here demonstrate that MES-3 is a highly divergent ortholog of SUZ12. MES-3 shares protein folds and conserved residues of key domains with SUZ12 and is predicted to interact with core PRC2 members similar to SUZ12 in human PRC2. Thus, in agreement with previous genetic and biochemical studies, we provide evidence that C. elegans contains a diverged yet evolutionary conserved core PRC2, like other animals.

Abstract Background DNA methylation is an important epigenetic control mechanism that in many fungi is restricted to genomic regions containing transposable elements (TEs). Two DNA methyltransferases, Dim2 and Dnmt5, are known to perform methylation at cytosines in fungi. While most ascomycete fungi encode both Dim2 and Dnmt5, only few functional studies have been performed in species containing both. Methods In this study, we report functional analysis of both Dim2 and Dnmt5 in the plant pathogenic fungus Verticillium dahliae. Results Our results show that Dim2, but not Dnmt5 or the putative sexual-cycle-related DNA methyltransferase Rid, is responsible for the majority of DNA methylation under the tested conditions. Single or double DNA methyltransferase mutants did not show altered development, virulence, or transcription of genes or TEs. In contrast, Hp1 and Dim5 mutants that are impacted in chromatin-associated processes upstream of DNA methylation are severely affected in development and virulence and display transcriptional reprogramming in specific hypervariable genomic regions (so-called adaptive genomic regions) that contain genes associated with host colonization. As these adaptive genomic regions are largely devoid of DNA methylation and of Hp1- and Dim5-associated heterochromatin, the differential transcription is likely caused by pleiotropic effects rather than by differential DNA methylation. Conclusion Overall, our study suggests that Dim2 is the main DNA methyltransferase in V. dahliae and, in conjunction with work on other fungi, is likely the main active DNMT in ascomycetes, irrespective of Dnmt5 presence. We speculate that Dnmt5 and Rid act under specific, presently enigmatic, conditions or, alternatively, act in DNA-associated processes other than DNA methylation.

AbstractMembers of the fungal genus PhyllostictaCitrusCitrus sinensisCitrus limonCitrus maximaPhyllostictaPhyllostictaPhyllosticta citrichinaensisPhyllostictaCitrusPhyllosticta citrichinaensisPhyllosticta citrichinaensisPhyllosticta

ABSTRACT Differential growth conditions typically trigger global transcriptional responses in filamentous fungi. Such fungal responses to environmental cues involve epigenetic regulation, including chemical histone modifications. It has been proposed that conditionally expressed genes, such as those that encode secondary metabolites but also effectors in pathogenic species, are often associated with a specific histone modification, lysine27 methylation of H3 (H3K27me3). However, thus far, no analyses on the global H3K27me3 profiles have been reported under differential growth conditions in order to assess if H3K27me3 dynamics govern differential transcription. Using chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing data from the plant-pathogenic fungus Verticillium dahliae grown in three in vitro cultivation media, we now show that a substantial number of the identified H3K27me3 domains globally display stable profiles among these growth conditions. However, we observe local quantitative differences in H3K27me3 ChIP-seq signals that are associated with a subset of differentially transcribed genes between media. Comparing the in vitro results to expression during plant infection suggests that in planta-induced genes may require chromatin remodeling to achieve expression. Overall, our results demonstrate that some loci display H3K27me3 dynamics associated with concomitant transcriptional variation, but many differentially expressed genes are associated with stable H3K27me3 domains. Thus, we conclude that while H3K27me3 is required for transcriptional repression, it does not appear that transcriptional activation requires the global erasure of H3K27me3. We propose that the H3K27me3 domains that do not undergo dynamic methylation may contribute to transcription through other mechanisms or may serve additional genomic regulatory functions. IMPORTANCE In many organisms, including filamentous fungi, epigenetic mechanisms that involve chemical and physical modifications of DNA without changing the genetic sequence have been implicated in transcriptional responses upon developmental or environmental cues. In fungi, facultative heterochromatin that can decondense to allow transcription in response to developmental changes or environmental stimuli is characterized by the trimethylation of lysine 27 on histone H3 (H3K27me3), and H3K27me3 has been implicated in transcriptional regulation, although the precise mechanisms and functions remain enigmatic. Based on ChIP and RNA sequencing data, we show for the soilborne broad-host-range vascular wilt plant-pathogenic fungus Verticillium dahliae that although some loci display H3K27me3 dynamics that can contribute to transcriptional variation, other loci do not show such a dependence. Thus, although we recognize that H3K27me3 is required for transcriptional repression, we also conclude that this mark is not a conditionally responsive global regulator of differential transcription upon responses to environmental cues.

In the last century, Fusarium wilt of banana (FWB) destroyed the banana cultivar Gros Michel. The Cavendish cultivars saved the global banana industry, and currently they dominate global production (~50%) and the export trade (~95%). However, a new strain called Tropical Race 4 (TR4) surfaced in the late 1960’s, spread globally and greatly damages Cavendish plantations as well as manifold local varieties that are primarily grown by small holders. Presently, there is no commercially available replacement for Cavendish and hence control strategies must be developed and implemented to manage FWB. Here, we studied whether it is possible to induce resistance to TR4 by pre-inoculations with different Fusarium spp. Only pre-treatments with an avirulent Race 1 strain significantly reduced disease development of TR4 in a Cavendish genotype and this effect was stable at various nutritional and pH conditions. We then used transcriptome analysis to study the molecular basis of this response. Several genes involved in plant defence responses were up-regulated during the initial stages of individual infections with TR4 and Race 1, as well as in combined treatments. In addition, a number of genes in the ethylene and jasmonate response pathways as well as several gibberellin synthesis associated genes were induced. We observed upregulation of RGA2 like genes in all treatments. Hence, RGA2 could be a key factor involved in both R1 and TR4 resistance. The data support the hypothesis that activating resistance to Race 1 in Cavendish bananas affects TR4 development and provide a first insight of gene expression during the interaction between various Fusarium spp. and banana.

Metabolism is the set of biochemical reactions of an organism that enables it to assimilate nutrients from its environment and to generate building blocks for growth and proliferation. It forms a complex network that is intertwined with the many molecular and cellular processes that take place within cells. Systems biology aims to capture the complexity of cells, organisms, or communities by reconstructing models based on information gathered by high-throughput analyses (omics data) and prior knowledge. One type of model is a genome-scale metabolic model (GEM) that allows studying the distributions of metabolic fluxes, i.e., the “mass-flow” through the network of biochemical reactions. GEMs are nowadays widely applied and have been reconstructed for various microbial pathogens, either in a free-living state or in interaction with their hosts, with the aim to gain insight into mechanisms of pathogenicity. In this review, we first introduce the principles of systems biology and GEMs. We then describe how metabolic modeling can contribute to unraveling microbial pathogenesis and host–pathogen interactions, with a specific focus on oomycete plant pathogens and in particular Phytophthora infestans. Subsequently, we review achievements obtained so far and identify and discuss potential pitfalls of current models. Finally, we propose a workflow for reconstructing high-quality GEMs and elaborate on the resources needed to advance a system biology approach aimed at untangling the intimate interactions between plants and pathogens.

ABSTRACT Hybridization is an important evolutionary mechanism that can enable organisms to adapt to environmental challenges. It has previously been shown that the fungal allodiploid species Verticillium longisporum, the causal agent of verticillium stem striping in rapeseed, originated from at least three independent hybridization events between two haploid Verticillium species. To reveal the impact of genome duplication as a consequence of hybridization, we studied the genome and transcriptome dynamics upon two independent V. longisporum hybridization events, represented by the hybrid lineages “A1/D1” and “A1/D3.” We show that V. longisporum genomes are characterized by extensive chromosomal rearrangements, including between parental chromosomal sets. V. longisporum hybrids display signs of evolutionary dynamics that are typically associated with the aftermath of allodiploidization, such as haploidization and more relaxed gene evolution. The expression patterns of the two subgenomes within the two hybrid lineages are more similar than those of the shared A1 parent between the two lineages, showing that the expression patterns of the parental genomes homogenized within a lineage. However, as genes that display differential parental expression in planta do not typically display the same pattern in vitro, we conclude that subgenome-specific responses occur in both lineages. Overall, our study uncovers genomic and transcriptomic plasticity during the evolution of the filamentous fungal hybrid V. longisporum and illustrates its adaptive potential. IMPORTANCE Verticillium is a genus of plant-associated fungi that includes a few plant pathogens that collectively affect a wide range of hosts. On several occasions, haploid Verticillium species hybridized into the stable allodiploid species Verticillium longisporum, which is, in contrast to haploid Verticillium species, a Brassicaceae specialist. Here, we studied the evolutionary genome and transcriptome dynamics of V. longisporum and the impact of the hybridization. V. longisporum genomes display a mosaic structure due to genomic rearrangements between the parental chromosome sets. Similar to other allopolyploid hybrids, V. longisporum displays an ongoing loss of heterozygosity and more relaxed gene evolution. Also, differential parental gene expression is observed, with enrichment for genes that encode secreted proteins. Intriguingly, the majority of these genes display subgenome-specific responses under differential growth conditions. In conclusion, hybridization has incited the genomic and transcriptomic plasticity that enables adaptation to environmental changes in a parental allele-specific fashion.

Genomes store information at scales beyond the linear nucleotide sequence, which impacts genome function at the level of an individual, while influences on populations and long-term genome function remains unclear. Here, we addressed how physical and chemical DNA characteristics influence genome evolution in the plant pathogenic fungus Verticillium dahliae. We identified incomplete DNA methylation of repetitive elements, associated with specific genomic compartments originally defined as Lineage-Specific (LS) regions that contain genes involved in host adaptation. Further chromatin characterization revealed associations with features such as H3 Lys-27 methylated histones (H3K27me3) and accessible DNA. Machine learning trained on chromatin data identified twice as much LS DNA as previously recognized, which was validated through orthogonal analysis, and we propose to refer to this DNA as adaptive genomic regions. Our results provide evidence that specific chromatin profiles define adaptive genomic regions, and highlight how different epigenetic factors contribute to the organization of these regions.

ABSTRACT Centromeres are chromosomal regions that are crucial for chromosome segregation during mitosis and meiosis, and failed centromere formation can contribute to chromosomal anomalies. Despite this conserved function, centromeres differ significantly between and even within species. Thus far, systematic studies into the organization and evolution of fungal centromeres remain scarce. In this study, we identified the centromeres in each of the 10 species of the fungal genus Verticillium and characterized their organization and evolution. Chromatin immunoprecipitation of the centromere-specific histone CenH3 (ChIP-seq) and chromatin conformation capture (Hi-C) followed by high-throughput sequencing identified eight conserved, large (∼150-kb), AT-, and repeat-rich regional centromeres that are embedded in heterochromatin in the plant pathogen Verticillium dahliae. Using Hi-C, we similarly identified repeat-rich centromeres in the other Verticillium species. Strikingly, a single degenerated long terminal repeat (LTR) retrotransposon is strongly associated with centromeric regions in some but not all Verticillium species. Extensive chromosomal rearrangements occurred during Verticillium evolution, of which some could be linked to centromeres, suggesting that centromeres contributed to chromosomal evolution. The size and organization of centromeres differ considerably between species, and centromere size was found to correlate with the genome-wide repeat content. Overall, our study highlights the contribution of repetitive elements to the diversity and rapid evolution of centromeres within the fungal genus Verticillium. IMPORTANCE The genus Verticillium contains 10 species of plant-associated fungi, some of which are notorious pathogens. Verticillium species evolved by frequent chromosomal rearrangements that contribute to genome plasticity. Centromeres are instrumental for separation of chromosomes during mitosis and meiosis, and failed centromere functionality can lead to chromosomal anomalies. Here, we used a combination of experimental techniques to identify and characterize centromeres in each of the Verticillium species. Intriguingly, we could strongly associate a single repetitive element to the centromeres of some of the Verticillium species. The presence of this element in the centromeres coincides with increased centromere sizes and genome-wide repeat expansions. Collectively, our findings signify a role of repetitive elements in the function, organization, and rapid evolution of centromeres in a set of closely related fungal species.

The Alternaria genus consists of saprophytic fungi as well as plant-pathogenic species that have significant economic impact. To date, the genomes of multiple Alternaria species have been sequenced. These studies have yielded valuable data for molecular studies on Alternaria fungi. However, most of the current Alternaria genome assemblies are highly fragmented, thereby hampering the identification of genes that are involved in causing disease. Here, we report a gapless genome assembly of A. solani, the causal agent of early blight in tomato and potato. The genome assembly is a significant step toward a better understanding of pathogenicity of A. solani.

ABSTRACT The oomycete pathogen Phytophthora infestans causes potato and tomato late blight, a disease that is a serious threat to agriculture. P. infestans is a hemibiotrophic pathogen, and during infection, it scavenges nutrients from living host cells for its own proliferation. To date, the nutrient flux from host to pathogen during infection has hardly been studied, and the interlinked metabolisms of the pathogen and host remain poorly understood. Here, we reconstructed an integrated metabolic model of P. infestans and tomato (Solanum lycopersicum) by integrating two previously published models for both species. We used this integrated model to simulate metabolic fluxes from host to pathogen and explored the topology of the model to study the dependencies of the metabolism of P. infestans on that of tomato. This showed, for example, that P. infestans, a thiamine auxotroph, depends on certain metabolic reactions of the tomato thiamine biosynthesis. We also exploited dual-transcriptome data of a time course of a full late blight infection cycle on tomato leaves and integrated the expression of metabolic enzymes in the model. This revealed profound changes in pathogen-host metabolism during infection. As infection progresses, P. infestans performs less de novo synthesis of metabolites and scavenges more metabolites from tomato. This integrated metabolic model for the P. infestans-tomato interaction provides a framework to integrate data and generate hypotheses about in planta nutrition of P. infestans throughout its infection cycle. IMPORTANCE Late blight disease caused by the oomycete pathogen Phytophthora infestans leads to extensive yield losses in tomato and potato cultivation worldwide. To effectively control this pathogen, a thorough understanding of the mechanisms shaping the interaction with its hosts is paramount. While considerable work has focused on exploring host defense mechanisms and identifying P. infestans proteins contributing to virulence and pathogenicity, the nutritional strategies of the pathogen are mostly unresolved. Genome-scale metabolic models (GEMs) can be used to simulate metabolic fluxes and help in unravelling the complex nature of metabolism. We integrated a GEM of tomato with a GEM of P. infestans to simulate the metabolic fluxes that occur during infection. This yields insights into the nutrients that P. infestans obtains during different phases of the infection cycle and helps in generating hypotheses about nutrition in planta.

Banana is the most popular and most exported fruit and also a major food crop for millions of people around the world. Despite its importance and the presence of serious disease threats, research into this crop is limited. One of those is Panama disease or Fusarium wilt. In the previous century Fusarium wilt wiped out the “Gros Michel” based banana industry in Central America. The epidemic was eventually quenched by planting “Cavendish” bananas. However, 50 years ago the disease recurred, but now on “Cavendish” bananas. Since then the disease has spread across South-East Asia, to the Middle-East and the Indian subcontinent and leaped into Africa. Here, we report the presence of Fusarium oxysporum f.sp. cubense Tropical Race 4 (Foc TR4) in “Cavendish” plantations in Laos, Myanmar, and Vietnam. A combination of classical morphology, DNA sequencing, and phenotyping assays revealed a very close relationship between the Foc TR4 strains in the entire Greater Mekong Subregion (GMS), which is increasingly prone to intensive banana production. Analyses of single-nucleotide polymorphisms enabled us to initiate a phylogeography of Foc TR4 across three geographical areas—GMS, Indian subcontinent, and the Middle East revealing three distinct Foc TR4 sub-lineages. Collectively, our data place these new incursions in a broader agroecological context and underscore the need for awareness campaigns and the implementation of validated quarantine measures to prevent further international dissemination of Foc TR4.

ABSTRACT Sensing external signals and transducing these into intracellular responses requires a molecular signaling system that is crucial for every living organism. Two important eukaryotic signal transduction pathways that are often interlinked are G-protein signaling and phospholipid signaling. Heterotrimeric G-protein subunits activated by G-protein-coupled receptors (GPCRs) are typical stimulators of phospholipid signaling enzymes such as phosphatidylinositol phosphate kinases (PIPKs) or phospholipase C (PLC). However, a direct connection between the two pathways likely exists in oomycetes and slime molds, as they possess a unique class of GPCRs that have a PIPK as an accessory domain. In principle, these so-called GPCR-PIPKs have the capacity of perceiving an external signal (via the GPCR domain) that, via PIPK, directly activates downstream phospholipid signaling. Here we reveal the sporadic occurrence of GPCR-PIPKs in all eukaryotic supergroups, except for plants. Notably, all species having GPCR-PIPKs are unicellular microorganisms that favor aquatic environments. Phylogenetic analysis revealed that GPCR-PIPKs are likely ancestral to eukaryotes and significantly expanded in the last common ancestor of oomycetes. In addition to GPCR-PIPKs, we identified five hitherto-unknown classes of GPCRs with accessory domains, four of which are universal players in signal transduction. Similarly to GPCR-PIPKs, this enables a direct coupling between extracellular sensing and downstream signaling. Overall, our findings point to an ancestral signaling system in eukaryotes where GPCR-mediated sensing is directly linked to downstream responses. IMPORTANCE G-protein-coupled receptors (GPCRs) are central sensors that activate eukaryotic signaling and are the primary targets of human drugs. In this report, we provide evidence for the widespread though limited presence of a novel class of GPCRs in a variety of unicellular eukaryotes. These include free-living organisms and organisms that are pathogenic for plants, animals, and humans. The novel GPCRs have a C-terminal phospholipid kinase domain, pointing to a direct link between sensing external signals via GPCRs and downstream intracellular phospholipid signaling. Genes encoding these receptors were likely present in the last common eukaryotic ancestor and were lost during the evolution of higher eukaryotes. We further describe five other types of GPCRs with a catalytic accessory domain, the so-called GPCR-bigrams, four of which may potentially have a role in signaling. These findings shed new light onto signal transduction in microorganisms and provide evidence for alternative eukaryotic signaling pathways.

Vascular wilts caused by Verticillium spp. are destructive plant diseases affecting hundreds of hosts. Only a few Verticillium spp. are causal agents of vascular wilt diseases, of which V. dahliae is the most notorious pathogen, and several V. dahliae genomes are available. In contrast, V. tricorpus is mainly known as a saprophyte and causal agent of opportunistic infections. Based on a hybrid approach that combines second and third generation sequencing, a near-gapless V. tricorpus genome assembly was obtained. With comparative genomics, we sought to identify genomic features in V. dahliae that confer the ability to cause vascular wilt disease. Unexpectedly, both species encode similar effector repertoires and share a genomic structure with genes encoding secreted proteins clustered in genomic islands. Intriguingly, V. tricorpus contains significantly fewer repetitive elements and an extended spectrum of secreted carbohydrate- active enzymes when compared with V. dahliae. In conclusion, we highlight the technical advances of a hybrid sequencing and assembly approach and show that the saprophyte V. tricorpus shares many hallmark features with the pathogen V. dahliae.

Bacterial flagellin molecules are strong inducers of innate immune responses in both mammals and plants. The opportunistic pathogen Pseudomonas aeruginosa secretes an alkaline protease called AprA that degrades flagellin monomers. Here, we show that AprA is widespread among a wide variety of bacterial species. In addition, we investigated the role of AprA in virulence of the bacterial plant pathogen P. syringae pv. tomato DC3000. The AprA-deficient DC3000 ΔaprA knockout mutant was significantly less virulent on both tomato and Arabidopsis thaliana. Moreover, infiltration of A. thaliana Col-0 leaves with DC3000 ΔaprA evoked a significantly higher level of expression of the defense-related genes FRK1 and PR-1 than did wild-type DC3000. In the flagellin receptor mutant fls2, pathogen virulence and defense-related gene activation did not differ between DC3000 and DC3000 ΔaprA. Together, these results suggest that AprA of DC3000 is important for evasion of recognition by the FLS2 receptor, allowing wild-type DC3000 to be more virulent on its host plant than AprA-deficient DC3000 ΔaprA. To provide further evidence for the role of DC3000 AprA in host immune evasion, we overexpressed the AprA inhibitory peptide AprI of DC3000 in A. thaliana to counteract the immune evasive capacity of DC3000 AprA. Ectopic expression of aprI in A. thaliana resulted in an enhanced level of resistance against wild-type DC3000, while the already elevated level of resistance against DC3000 ΔaprA remained unchanged. Together, these results indicate that evasion of host immunity by the alkaline protease AprA is important for full virulence of strain DC3000 and likely acts by preventing flagellin monomers from being recognized by its cognate immune receptor.

ABSTRACT Next-generation sequencing (NGS) technologies have increased the scalability, speed, and resolution of genomic sequencing and, thus, have revolutionized genomic studies. However, eukaryotic genome sequencing initiatives typically yield considerably fragmented genome assemblies. Here, we assessed various state-of-the-art sequencing and assembly strategies in order to produce a contiguous and complete eukaryotic genome assembly, focusing on the filamentous fungus Verticillium dahliae. Compared with Illumina-based assemblies of the V. dahliae genome, hybrid assemblies that also include PacBio-generated long reads establish superior contiguity. Intriguingly, provided that sufficient sequence depth is reached, assemblies solely based on PacBio reads outperform hybrid assemblies and even result in fully assembled chromosomes. Furthermore, the addition of optical map data allowed us to produce a gapless and complete V. dahliae genome assembly of the expected eight chromosomes from telomere to telomere. Consequently, we can now study genomic regions that were previously not assembled or poorly assembled, including regions that are populated by repetitive sequences, such as transposons, allowing us to fully appreciate an organism's biological complexity. Our data show that a combination of PacBio-generated long reads and optical mapping can be used to generate complete and gapless assemblies of fungal genomes. IMPORTANCE Studying whole-genome sequences has become an important aspect of biological research. The advent of next-generation sequencing (NGS) technologies has nowadays brought genomic science within reach of most research laboratories, including those that study nonmodel organisms. However, most genome sequencing initiatives typically yield (highly) fragmented genome assemblies. Nevertheless, considerable relevant information related to genome structure and evolution is likely hidden in those nonassembled regions. Here, we investigated a diverse set of strategies to obtain gapless genome assemblies, using the genome of a typical ascomycete fungus as the template. Eventually, we were able to show that a combination of PacBio-generated long reads and optical mapping yields a gapless telomere-to-telomere genome assembly, allowing in-depth genome analyses to facilitate functional studies into an organism's biology.

Cytokinins, a class of phytohormones, are adenine derivatives common to many different organisms. In plants, these play a crucial role as regulators of plant development and the reaction to abiotic and biotic stress. Key enzymes in the cytokinin synthesis and degradation in modern land plants are the isopentyl transferases and the cytokinin dehydrogenases, respectively. Their encoding genes have been probably introduced into the plant lineage during the primary endosymbiosis. To shed light on the evolution of these proteins, the genes homologous to plant adenylate isopentenyl transferase and cytokinin dehydrogenase were amplified from the genomic DNA of cyanobacterium Nostoc sp. PCC 7120 and expressed in Escherichia coli. The putative isopentenyl transferase was shown to be functional in a biochemical assay. In contrast, no enzymatic activity was detected for the putative cytokinin dehydrogenase, even though the principal domains necessary for its function are present. Several mutant variants, in which conserved amino acids in land plant cytokinin dehydrogenases had been restored, were inactive. A combination of experimental data with phylogenetic analysis indicates that adenylate-type isopentenyl transferases might have evolved several times independently. While the Nostoc genome contains a gene coding for protein with characteristics of cytokinin dehydrogenase, the organism is not able to break down cytokinins in the way shown for land plants.

The biotrophic plant pathogen Hyaloperonospora arabidopsidis produces a set of putative effector proteins that contain the conserved RXLR motif. For most of these RXLR proteins the role during infection is unknown. Thirteen RXLR proteins from H. arabidopsidis strain Waco9 were analyzed for sequence similarities and tested for a role in virulence. The thirteen RXLR proteins displayed conserved N-termini and this N-terminal conservation was also found in the 134 predicted RXLR genes from the genome of H. arabidopsidis strain Emoy2. To investigate the effects of single RXLR effector proteins on plant defense responses, thirteen H. arabidopsidis Waco9 RXLR genes were expressed in Arabidopsis thaliana. Subsequently, these plants were screened for altered susceptibility to the oomycetes H. arabidopsidis and Phytophthora capsici, and the bacterial pathogen Pseudomonas syringae. Additionally, the effect of the RXLR proteins on flg22-triggered basal immune responses was assessed. Multifactorial analysis of results collated from all experiments revealed that, except for RXLR20, all RXLR effector proteins tested affected plant immunity. For RXLR9 this was confirmed using a P. syringae ΔCEL-mediated effector delivery system. Together, the results show that many H. arabidopsidis RXLR effectors have small effects on the plant immune response, suggesting that suppression of host immunity by this biotrophic pathogen is likely to be caused by the combined actions of effectors.

Oomycetes in the class Saprolegniomycetidae of the Eukaryotic kingdom Stramenopila have evolved as severe pathogens of amphibians, crustaceans, fish and insects, resulting in major losses in aquaculture and damage to aquatic ecosystems. We have sequenced the 63 Mb genome of the fresh water fish pathogen, Saprolegnia parasitica. Approximately 1/3 of the assembled genome exhibits loss of heterozygosity, indicating an efficient mechanism for revealing new variation. Comparison of S. parasitica with plant pathogenic oomycetes suggests that during evolution the host cellular environment has driven distinct patterns of gene expansion and loss in the genomes of plant and animal pathogens. S. parasitica possesses one of the largest repertoires of proteases (270) among eukaryotes that are deployed in waves at different points during infection as determined from RNA-Seq data. In contrast, despite being capable of living saprotrophically, parasitism has led to loss of inorganic nitrogen and sulfur assimilation pathways, strikingly similar to losses in obligate plant pathogenic oomycetes and fungi. The large gene families that are hallmarks of plant pathogenic oomycetes such as Phytophthora appear to be lacking in S. parasitica, including those encoding RXLR effectors, Crinkler's, and Necrosis Inducing-Like Proteins (NLP). S. parasitica also has a very large kinome of 543 kinases, 10% of which is induced upon infection. Moreover, S. parasitica encodes several genes typical of animals or animal-pathogens and lacking from other oomycetes, including disintegrins and galactose-binding lectins, whose expression and evolutionary origins implicate horizontal gene transfer in the evolution of animal pathogenesis in S. parasitica.

The gene composition of present-day genomes has been shaped by a complicated evolutionary history, resulting in diverse distributions of genes across genomes. The pattern of presence and absence of a gene in different genomes is called its phylogenetic profile. It has been shown that proteins whose encoding genes have highly similar profiles tend to be functionally related: As these genes were gained and lost together, their encoded proteins can probably only perform their full function if both are present. However, a large proportion of genes encoding interacting proteins do not have matching profiles. In this study, we analysed one possible reason for this, namely that phylogenetic profiles can be affected by multi-functional proteins such as shared subunits of two or more protein complexes. We found that by considering triplets of proteins, of which one protein is multi-functional, a large fraction of disturbed co-occurrence patterns can be explained.