Advances in molecular pathogenesis of hidradenitis suppurativa: Dysregulated keratins and ECM signaling.

Hidradenitis suppurativa (HS) is characterized by deep-seated, highly inflamed, and painful lumps/abscesses, fistulae, and sinus tracts that grow extensively deep in the dermis and are highly immunogenic in nature. In about one-third of the HS patients there is strong evidence for the role of γ-secretase mutations along with dysregulated Notch signaling. However, the contribution of dysregulated Notch signaling in HS pathogenesis in relation to hair follicle alterations and hyper-activation of the immune system remains undefined. A genome-wide association study (GWAS), proteomic data and functional investigations of identified sequence variants in HS pathology are not fully revealing. The disease initiation or progression may involve bacterial infection besides intrinsic functional defects in keratinocytes, which may be key to further exacerbate immune cell infiltration and cytokine production in and around the lesional tissue. The absence of a suitable animal model that could fully recapitulate the pathogenesis of HS is a major impediment for proper understanding the underlying mechanisms and development of effective treatments. The presence of extracellular matrix (ECM) degradation products along with dysregulation in keratinocytes and, dermal fibroblasts ultimately affect immune regulation and are various components of HS pathogenesis. Bacterial infection further exacerbates the complexity of the disease progression. While anti-TNFα therapy shows partial efficacy, treatment to cure HS is absent. Multiple clinical trials targeting various cytokines, complement C5a and ECM products are in progress. This review provides state-of-the-art information on these aspects with a focus on dysregulated keratinocyte and immune cells; and role of ECM, and Keratin functions in this regard.

NPR1 in Plant Defense: It's Not over 'til It's Turned over.

NPR1 is a key transcriptional coregulator in plant defense responses. In this issue, Spoel et al. (2009) demonstrate that proteasome-mediated degradation of NPR1 in the nucleus promotes efficient expression of defense response genes following infection and prevents spurious activation of defensive responses in the absence of infection.

Publisher Correction: An extracellular network of Arabidopsis leucine-rich repeat receptor kinases.

In this Letter, an incorrect version of the Supplementary Information file was inadvertently used, which contained several errors. The details of references 59-65 were missing from the end of the Supplementary Discussion section on page 4. In addition, the section 'Text 3. Y2H on ICD interactions' incorrectly referred to 'Extended Data Fig. 4d' instead of 'Extended Data Fig. 3d' on page 3. Finally, the section 'Text 4. Interaction network analysis' incorrectly referred to 'Fig. 1b and Extended Data Fig. 6' instead of 'Fig. 2b and Extended Data Fig. 7' on page 3. These errors have all been corrected in the Supplementary Information.

An extracellular network of Arabidopsis leucine-rich repeat receptor kinases.

The cells of multicellular organisms receive extracellular signals using surface receptors. The extracellular domains (ECDs) of cell surface receptors function as interaction platforms, and as regulatory modules of receptor activation. Understanding how interactions between ECDs produce signal-competent receptor complexes is challenging because of their low biochemical tractability. In plants, the discovery of ECD interactions is complicated by the massive expansion of receptor families, which creates tremendous potential for changeover in receptor interactions. The largest of these families in Arabidopsis thaliana consists of 225 evolutionarily related leucine-rich repeat receptor kinases (LRR-RKs), which function in the sensing of microorganisms, cell expansion, stomata development and stem-cell maintenance. Although the principles that govern LRR-RK signalling activation are emerging, the systems-level organization of this family of proteins is unknown. Here, to address this, we investigated 40,000 potential ECD interactions using a sensitized high-throughput interaction assay, and produced an LRR-based cell surface interaction network (CSILRR) that consists of 567 interactions. To demonstrate the power of CSILRR for detecting biologically relevant interactions, we predicted and validated the functions of uncharacterized LRR-RKs in plant growth and immunity. In addition, we show that CSILRR operates as a unified regulatory network in which the LRR-RKs most crucial for its overall structure are required to prevent the aberrant signalling of receptors that are several network-steps away. Thus, plants have evolved LRR-RK networks to process extracellular signals into carefully balanced responses.

Epigenetic regulation in the pathogenesis of non-melanoma skin cancer.

Understanding of cancer with the help of ever-expanding cutting edge technological tools and bioinformatics is revolutionizing modern cancer research by broadening the space of discovery window of various genomic and epigenomic processes. Genomics data integrated with multi-omics layering have advanced cancer research. Uncovering such layers of genetic mutations/modifications, epigenetic regulation and their role in the complex pathophysiology of cancer progression could lead to novel therapeutic interventions. Although a plethora of literature is available in public domain defining the role of various tumor driver gene mutations, understanding of epigenetic regulation of cancer is still emerging. This review focuses on epigenetic regulation association with the pathogenesis of non-melanoma skin cancer (NMSC). NMSC has higher prevalence in Caucasian populations compared to other races. Due to lack of proper reporting to cancer registries, the incidence rates for NMSC worldwide cannot be accurately estimated. However, this is the most common neoplasm in humans, and millions of new cases per year are reported in the United States alone. In organ transplant recipients, the incidence of NMSC particularly of squamous cell carcinoma (SCC) is very high and these SCCs frequently become metastatic and lethal. Understanding of solar ultraviolet (UV) light-induced damage and impaired DNA repair process leading to DNA mutations and nuclear instability provide an insight into the pathogenesis of metastatic neoplasm. This review discusses the recent advances in the field of epigenetics of NMSCs. Particularly, the role of DNA methylation, histone hyperacetylation and non-coding RNA such as long-chain noncoding (lnc) RNAs, circular RNAs and miRNA in the disease progression are summarized.

Phosphoprotein-based biomarkers as predictors for cancer therapy.

Disparities in cancer patient responses have prompted widespread searches to identify differences in sensitive vs. nonsensitive populations and form the basis of personalized medicine. This customized approach is dependent upon the development of pathway-specific therapeutics in conjunction with biomarkers that predict patient responses. Here, we show that Cdk5 drives growth in subgroups of patients with multiple types of neuroendocrine neoplasms. Phosphoproteomics and high throughput screening identified phosphorylation sites downstream of Cdk5. These phosphorylation events serve as biomarkers and effectively pinpoint Cdk5-driven tumors. Toward achieving targeted therapy, we demonstrate that mouse models of neuroendocrine cancer are responsive to selective Cdk5 inhibitors and biomimetic nanoparticles are effective vehicles for enhanced tumor targeting and reduction of drug toxicity. Finally, we show that biomarkers of Cdk5-dependent tumors effectively predict response to anti-Cdk5 therapy in patient-derived xenografts. Thus, a phosphoprotein-based diagnostic assay combined with Cdk5-targeted therapy is a rational treatment approach for neuroendocrine malignancies.

Network Biology Analyses and Dynamic Modeling of Gene Regulatory Networks under Drought Stress Reveal Major Transcriptional Regulators in Arabidopsis.

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs.

Ranking Plant Network Nodes Based on Their Centrality Measures.

Biological networks are often large and complex, making it difficult to accurately identify the most important nodes. Node prioritization algorithms are used to identify the most influential nodes in a biological network by considering their relationships with other nodes. These algorithms can help us understand the functioning of the network and the role of individual nodes. We developed CentralityCosDist, an algorithm that ranks nodes based on a combination of centrality measures and seed nodes. We applied this and four other algorithms to protein-protein interactions and co-expression patterns in Arabidopsis thaliana using pathogen effector targets as seed nodes. The accuracy of the algorithms was evaluated through functional enrichment analysis of the top 10 nodes identified by each algorithm. Most enriched terms were similar across algorithms, except for DIAMOnD. CentralityCosDist identified more plant-pathogen interactions and related functions and pathways compared to the other algorithms.

Combined inhibition of BET bromodomain and mTORC1/2 provides therapeutic advantage for rhabdomyosarcoma by switching cell death mechanism.

Aberrant activation of multiple complex signaling pathways underlies the pathogenesis of rhabdomyosarcoma (RMS), which remains a cause of mortality in approximately 30% of children with RMS. Bromodomain and extraterminal (BET) domain chromatin remodeling regulates several of these pathways. Here, we targeted bromodomain 4 (BRD4) in combination with another molecular metabolic tumor driver, the Akt/mTOR signaling pathway, to provide a highly effective treatment for this neoplasm. We demonstrated that a nexus of these two molecular pathways underlies RMS pathogenesis. Our data show that the combined inhibition of the BET bromodomain and mTORC1/2 signaling abrogates aggressive RMS growth. Thus, the bromodomain inhibitor RVX-208 significantly augmented the therapeutic effects of the dual mTORC1/2 inhibitors, OSI-027 and PP242, both in vitro and in a human xenograft murine model. Drug-treated residual tumors showed a decrease in the activation of underlying signaling mechanisms characterized by a reduction in the expression of p-AKT, p-mTOR, p-p70S6K, cyclin D1, and proliferation. Our ChIP-seq data demonstrated that RVX-208 effectively blocked BRD4 occupancy on its target promoters. ChIP-qPCR assays further confirmed that RVX-208 treatment resulted in a significant decrease in H3K27ac and H4K8ac signals at their target loci. While single RVX-208 treatment induces apoptosis and a single mTORC1/2 inhibitor induces macropinocytosis, their combined treatment led to necroptosis-mediated cell death. These data suggest that combined treatment with drugs targeting BRD4 and mTORC1/2 may be an effective therapeutic intervention for drug-resistant RMS.

Molecular insight into cotton leaf curl geminivirus disease resistance in cultivated cotton (Gossypium hirsutum).

Cultivated cotton (Gossypium hirsutum) is the most important fibre crop in the world. Cotton leaf curl disease (CLCuD) is the major limiting factor and a threat to textile industry in India and Pakistan. All the local cotton cultivars exhibit moderate to no resistance against CLCuD. In this study, we evaluated an exotic cotton accession Mac7 as a resistance source to CLCuD by challenging it with viruliferous whiteflies and performing qPCR to evaluate the presence/absence and relative titre of CLCuD-associated geminiviruses/betasatellites. The results indicated that replication of pathogenicity determinant betasatellite is significantly attenuated in Mac7 and probably responsible for resistance phenotype. Afterwards, to decipher the genetic basis of CLCuD resistance in Mac7, we performed RNA sequencing on CLCuD-infested Mac7 and validated RNA-Seq data with qPCR on 24 independent genes. We performed co-expression network and pathway analysis for regulation of geminivirus/betasatellite-interacting genes. We identified nine novel modules with 52 hubs of highly connected genes in network topology within the co-expression network. Analysis of these hubs indicated the differential regulation of auxin stimulus and cellular localization pathways in response to CLCuD. We also analysed the differential regulation of geminivirus/betasatellite-interacting genes in Mac7. We further performed the functional validation of selected candidate genes via virus-induced gene silencing (VIGS). Finally, we evaluated the genomic context of resistance responsive genes and found that these genes are not specific to A or D sub-genomes of G. hirsutum. These results have important implications in understanding CLCuD resistance mechanism and developing a durable resistance in cultivated cotton.

Systems Biology and Machine Learning in Plant-Pathogen Interactions.

Systems biology is an inclusive approach to study the static and dynamic emergent properties on a global scale by integrating multiomics datasets to establish qualitative and quantitative associations among multiple biological components. With an abundance of improved high throughput -omics datasets, network-based analyses and machine learning technologies are playing a pivotal role in comprehensive understanding of biological systems. Network topological features reveal most important nodes within a network as well as prioritize significant molecular components for diverse biological networks, including coexpression, protein-protein interaction, and gene regulatory networks. Machine learning techniques provide enormous predictive power through specific feature extraction from biological data. Deep learning, a subtype of machine learning, has plausible future applications because a domain expert for feature extraction is not needed in this algorithm. Inspired by diverse domains of biology, we here review classic systems biology techniques applied in plant immunity thus far. We also discuss additional advanced approaches in both graph theory and machine learning, which may provide new insights for understanding plant-microbe interactions. Finally, we propose a hybrid approach in plant immune systems that harnesses the power of both network biology and machine learning, with a potential to be applicable to both model systems and agronomically important crop plants.

Direct Regulation of the EFR-Dependent Immune Response by Arabidopsis TCP Transcription Factors.

One layer of the innate immune system allows plants to recognize pathogen-associated molecular patterns (PAMPS), activating a defense response known as PAMP-triggered immunity (PTI). Maintaining an active immune response, however, comes at the cost of plant growth and development; accordingly, optimization of the balance between defense and development is critical to plant fitness. The TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) transcription factor family consists of well-characterized transcriptional regulators of plant development and morphogenesis. The three closely related class I TCP transcription factors TCP8, TCP14, and TCP15 have also been implicated in the regulation of effector-triggered immunity, but there has been no previous characterization of PTI-related phenotypes. To identify TCP targets involved in PTI, we screened a PAMP-induced gene promoter library in a yeast one-hybrid assay and identified interactions of these three TCPs with the EF-Tu RECEPTOR (EFR) promoter. The direct interactions between TCP8 and EFR were confirmed to require an intact TCP binding site in planta. A tcp8 tcp14 tcp15 triple mutant was impaired in EFR-dependent PTI and exhibited reduced levels of PATHOGENESIS-RELATED PROTEIN 2 and induction of EFR expression after elicitation with elf18 but also increased production of reactive oxygen species relative to Col-0. Our data support an increasingly complex role for TCPs at the nexus of plant development and defense.

Pseudomonas syringae type III effector HopAF1 suppresses plant immunity by targeting methionine recycling to block ethylene induction.

HopAF1 is a type III effector protein of unknown function encoded in the genomes of several strains of Pseudomonas syringae and other plant pathogens. Structural modeling predicted that HopAF1 is closely related to deamidase proteins. Deamidation is the irreversible substitution of an amide group with a carboxylate group. Several bacterial virulence factors are deamidases that manipulate the activity of specific host protein substrates. We identified Arabidopsis methylthioadenosine nucleosidase proteins MTN1 and MTN2 as putative targets of HopAF1 deamidation. MTNs are enzymes in the Yang cycle, which is essential for the high levels of ethylene biosynthesis in Arabidopsis We hypothesized that HopAF1 inhibits the host defense response by manipulating MTN activity and consequently ethylene levels. We determined that bacterially delivered HopAF1 inhibits ethylene biosynthesis induced by pathogen-associated molecular patterns and that Arabidopsis mtn1 mtn2 mutant plants phenocopy the effect of HopAF1. Furthermore, we identified two conserved asparagines in MTN1 and MTN2 from Arabidopsis that confer loss of function phenotypes when deamidated via site-specific mutation. These residues are potential targets of HopAF1 deamidation. HopAF1-mediated manipulation of Yang cycle MTN proteins is likely an evolutionarily conserved mechanism whereby HopAF1 orthologs from multiple plant pathogens contribute to disease in a large variety of plant hosts.

Cytokinin Response Factor 6 Represses Cytokinin-Associated Genes during Oxidative Stress.

Cytokinin is a phytohormone that is well known for its roles in numerous plant growth and developmental processes, yet it has also been linked to abiotic stress response in a less defined manner. Arabidopsis (Arabidopsis thaliana) Cytokinin Response Factor 6 (CRF6) is a cytokinin-responsive AP2/ERF-family transcription factor that, through the cytokinin signaling pathway, plays a key role in the inhibition of dark-induced senescence. CRF6 expression is also induced by oxidative stress, and here we show a novel function for CRF6 in relation to oxidative stress and identify downstream transcriptional targets of CRF6 that are repressed in response to oxidative stress. Analysis of transcriptomic changes in wild-type and crf6 mutant plants treated with H2O2 identified CRF6-dependent differentially expressed transcripts, many of which were repressed rather than induced. Moreover, many repressed genes also show decreased expression in 35S:CRF6 overexpressing plants. Together, these findings suggest that CRF6 functions largely as a transcriptional repressor. Interestingly, among the H2O2 repressed CRF6-dependent transcripts was a set of five genes associated with cytokinin processes: (signaling) ARR6, ARR9, ARR11, (biosynthesis) LOG7, and (transport) ABCG14. We have examined mutants of these cytokinin-associated target genes to reveal novel connections to oxidative stress. Further examination of CRF6-DNA interactions indicated that CRF6 may regulate its targets both directly and indirectly. Together, this shows that CRF6 functions during oxidative stress as a negative regulator to control this cytokinin-associated module of CRF6-dependent genes and establishes a novel connection between cytokinin and oxidative stress response.

Heterotrimeric G-proteins in Picea abies and their regulation in response to Heterobasidion annosum s.l. infection.

BACKGROUND: Heterotrimeric G-proteins are important signalling switches, present in all eukaryotic kingdoms. In plants they regulate several developmental functions and play an important role in plant-microbe interactions. The current knowledge on plant G-proteins is mostly based on model angiosperms and little is known about the G-protein repertoire and function in other lineages. In this study we investigate the heterotrimeric G-protein subunit repertoire in Pinaceae, including phylogenetic relationships, radiation and sequence diversity levels in relation to other plant linages. We also investigate functional diversification of the G-protein complex in Picea abies by analysing transcriptional regulation of the G-protein subunits in different tissues and in response to pathogen infection. RESULTS: A full repertoire of G-protein subunits in several conifer species were identified in silico. The full-length P. abies coding regions of one Gα-, one Gß- and four Gγ-subunits were cloned and sequenced. The phylogenetic analysis of the Gγ-subunits showed that PaGG1 clustered with A-type-like subunits, PaGG3 and PaGG4 clustered with C-type-like subunits, while PaGG2 and its orthologs represented a novel conifer-specific putative Gγ-subunit type. Gene expression analyses by quantitative PCR of P. abies G-protein subunits showed specific up-regulation of the Gα-subunit gene PaGPA1 and the Gγ-subunit gene PaGG1 in response to Heterobasidion annosum sensu lato infection. CONCLUSIONS: Conifers possess a full repertoire of G-protein subunits. The differential regulation of PaGPA1 and PaGG1 indicates that the heterotrimeric G-protein complex represents a critical linchpin in Heterobasidion annosum s.l. perception and downstream signaling in P. abies.

Expression-based network biology identifies immune-related functional modules involved in plant defense.

BACKGROUND: Plants respond to diverse environmental cues including microbial perturbations by coordinated regulation of thousands of genes. These intricate transcriptional regulatory interactions depend on the recognition of specific promoter sequences by regulatory transcription factors. The combinatorial and cooperative action of multiple transcription factors defines a regulatory network that enables plant cells to respond to distinct biological signals. The identification of immune-related modules in large-scale transcriptional regulatory networks can reveal the mechanisms by which exposure to a pathogen elicits a precise phenotypic immune response. RESULTS: We have generated a large-scale immune co-expression network using a comprehensive set of Arabidopsis thaliana (hereafter Arabidopsis) transcriptomic data, which consists of a wide spectrum of immune responses to pathogens or pathogen-mimicking stimuli treatments. We employed both linear and non-linear models to generate Arabidopsis immune co-expression regulatory (AICR) network. We computed network topological properties and ascertained that this newly constructed immune network is densely connected, possesses hubs, exhibits high modularity, and displays hallmarks of a "real" biological network. We partitioned the network and identified 156 novel modules related to immune functions. Gene Ontology (GO) enrichment analyses provided insight into the key biological processes involved in determining finely tuned immune responses. We also developed novel software called OCCEAN (One Click Cis-regulatory Elements ANalysis) to discover statistically enriched promoter elements in the upstream regulatory regions of Arabidopsis at a whole genome level. We demonstrated that OCCEAN exhibits higher precision than the existing promoter element discovery tools. In light of known and newly discovered cis-regulatory elements, we evaluated biological significance of two key immune-related functional modules and proposed mechanism(s) to explain how large sets of diverse GO genes coherently function to mount effective immune responses. CONCLUSIONS: We used a network-based, top-down approach to discover immune-related modules from transcriptomic data in Arabidopsis. Detailed analyses of these functional modules reveal new insight into the topological properties of immune co-expression networks and a comprehensive understanding of multifaceted plant defense responses. We present evidence that our newly developed software, OCCEAN, could become a popular tool for the Arabidopsis research community as well as potentially expand to analyze other eukaryotic genomes.

Duckweed: Beyond an Efficient Plant Model System.

Duckweed (Lemnaceae) rises as a crucial model system due to its unique characteristics and wide-ranging utility. The significance of physiological research and phytoremediation highlights the intricate potential of duckweed in the current era of plant biology. Special attention to duckweed has been brought due to its distinctive features of nutrient uptake, ion transport dynamics, detoxification, intricate signaling, and stress tolerance. In addition, duckweed can alleviate environmental pollutants and enhance sustainability by participating in bioremediation processes and wastewater treatment. Furthermore, insights into the genomic complexity of Lemnaceae species and the flourishing field of transgenic development highlight the opportunities for genetic manipulation and biotechnological innovations. Novel methods for the germplasm conservation of duckweed can be adopted to preserve genetic diversity for future research endeavors and breeding programs. This review centers around prospects in duckweed research promoting interdisciplinary collaborations and technological advancements to drive its full potential as a model organism.

Molecular heterogeneity of quiescent melanocyte stem cells revealed by single-cell RNA-sequencing.

Melanocyte stem cells (McSCs) of the hair follicle are a rare cell population within the skin and are notably underrepresented in whole-skin, single-cell RNA sequencing (scRNA-seq) datasets. Using a cell enrichment strategy to isolate KIT+/CD45- cells from the telogen skin of adult female C57BL/6J mice, we evaluated the transcriptional landscape of quiescent McSCs (qMcSCs) at high resolution. Through this evaluation, we confirmed existing molecular signatures for qMcCS subpopulations (e.g., Kit+, Cd34+/-, Plp1+, Cd274+/-, Thy1+, Cdh3+/-) and identified novel qMcSC subpopulations, including two that differentially regulate their immune privilege status. Within qMcSC subpopulations, we also predicted melanocyte differentiation potential, neural crest potential, and quiescence depth. Taken together, the results demonstrate that the qMcSC population is heterogeneous and future studies focused on investigating changes in qMcSCs should consider changes in subpopulation composition.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes.

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Dynamic evolution of pathogenicity revealed by sequencing and comparative genomics of 19 Pseudomonas syringae isolates.

Closely related pathogens may differ dramatically in host range, but the molecular, genetic, and evolutionary basis for these differences remains unclear. In many Gram- negative bacteria, including the phytopathogen Pseudomonas syringae, type III effectors (TTEs) are essential for pathogenicity, instrumental in structuring host range, and exhibit wide diversity between strains. To capture the dynamic nature of virulence gene repertoires across P. syringae, we screened 11 diverse strains for novel TTE families and coupled this nearly saturating screen with the sequencing and assembly of 14 phylogenetically diverse isolates from a broad collection of diseased host plants. TTE repertoires vary dramatically in size and content across all P. syringae clades; surprisingly few TTEs are conserved and present in all strains. Those that are likely provide basal requirements for pathogenicity. We demonstrate that functional divergence within one conserved locus, hopM1, leads to dramatic differences in pathogenicity, and we demonstrate that phylogenetics-informed mutagenesis can be used to identify functionally critical residues of TTEs. The dynamism of the TTE repertoire is mirrored by diversity in pathways affecting the synthesis of secreted phytotoxins, highlighting the likely role of both types of virulence factors in determination of host range. We used these 14 draft genome sequences, plus five additional genome sequences previously reported, to identify the core genome for P. syringae and we compared this core to that of two closely related non-pathogenic pseudomonad species. These data revealed the recent acquisition of a 1 Mb megaplasmid by a sub-clade of cucumber pathogens. This megaplasmid encodes a type IV secretion system and a diverse set of unknown proteins, which dramatically increases both the genomic content of these strains and the pan-genome of the species.