Next-generation yeast-two-hybrid analysis with Y2H-SCORES identifies novel interactors of the MLA immune receptor

Protein-protein interaction networks are one of the most effective representations of cellular behavior. In order to build these models, high-throughput techniques are required. Next-generation interaction screening (NGIS) protocols that combine yeast two-hybrid (Y2H) with deep sequencing are promising approaches to generate interactome networks in any organism. However, challenges remain to mining reliable information from these screens and thus, limit its broader implementation. Here, we present a computational framework, designated Y2H-SCORES, for analyzing high-throughput Y2H screens. Y2H-SCORES considers key aspects of NGIS experimental design and important characteristics of the resulting data that distinguish it from RNA-seq expression datasets. Three quantitative ranking scores were implemented to identify interacting partners, comprising: 1) significant enrichment under selection for positive interactions, 2) degree of interaction specificity among multi-bait comparisons, and 3) selection of in-frame interactors. Using simulation and an empirical dataset, we provide a quantitative assessment to predict interacting partners under a wide range of experimental scenarios, facilitating independent confirmation by one-to-one bait-prey tests. Simulation of Y2H-NGIS enabled us to identify conditions that maximize detection of true interactors, which can be achieved with protocols such as prey library normalization, maintenance of larger culture volumes and replication of experimental treatments. Y2H-SCORES can be implemented in different yeast-based interaction screenings, with an equivalent or superior performance than existing methods. Proof-of-concept was demonstrated by discovery and validation of novel interactions between the barley nucleotide-binding leucine-rich repeat (NLR) immune receptor MLA6, and fourteen proteins, including those that function in signaling, transcriptional regulation, and intracellular trafficking.

[1]  P. Schulze-Lefert,et al.  Structure-Function Analysis of Barley NLR Immune Receptor MLA10 Reveals Its Cell Compartment Specific Activity in Cell Death and Disease Resistance , 2012, PLoS pathogens.

[2]  W. Huber,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets , 2011 .

[3]  P. Schulze-Lefert,et al.  The MLA6 coiled-coil, NBS-LRR protein confers AvrMla6-dependent resistance specificity to Blumeria graminis f. sp. hordei in barley and wheat. , 2001, The Plant journal : for cell and molecular biology.

[4]  John K. McCooke,et al.  A chromosome conformation capture ordered sequence of the barley genome , 2017, Nature.

[5]  Mike Tyers,et al.  BioGRID: a general repository for interaction datasets , 2005, Nucleic Acids Res..

[6]  Ulrich Schlecht,et al.  A scalable double-barcode sequencing platform for characterization of dynamic protein-protein interactions , 2017, Nature Communications.

[7]  Y. Qi,et al.  Reconstitution and structure of a plant NLR resistosome conferring immunity , 2019, Science.

[8]  R. Wing,et al.  Genome Dynamics and Evolution of the Mla (Powdery Mildew) Resistance Locus in Barley Online version contains Web-only data. Article, publication date, and citation information can be found at www.plantcell.org/cgi/doi/10.1105/tpc.002238. , 2002, The Plant Cell Online.

[9]  Xueping Zhou,et al.  The Rice Dynamin-Related Protein OsDRP1E Negatively Regulates Programmed Cell Death by Controlling the Release of Cytochrome c from Mitochondria , 2017, PLoS pathogens.

[10]  K. Rahman,et al.  Development and application of a recombination-based library versus library high- throughput yeast two-hybrid (RLL-Y2H) screening system , 2017, Nucleic acids research.

[11]  Hsin-Hung Lin,et al.  Transcriptome dynamics of developing maize leaves and genomewide prediction of cis elements and their cognate transcription factors , 2015, Proceedings of the National Academy of Sciences.

[12]  Kai Li,et al.  A DnaJ protein that interacts with soybean mosaic virus coat protein serves as a key susceptibility factor for viral infection. , 2020, Virus research.

[13]  A. Nebenführ,et al.  Kinesins and Myosins: Molecular Motors that Coordinate Cellular Functions in Plants. , 2018, Annual review of plant biology.

[14]  Xing-Ming Zhao,et al.  PPIM : A protein-protein interaction database for Maize 11 12 , 2015 .

[15]  Roger P. Wise,et al.  NGPINT: A Next-generation protein-protein interaction software , 2020, bioRxiv.

[16]  Roger P Wise,et al.  Interchromosomal Transfer of Immune Regulation During Infection of Barley with the Powdery Mildew Pathogen , 2017, G3: Genes, Genomes, Genetics.

[17]  Anthony Skjellum,et al.  Expression-based network biology identifies immune-related functional modules involved in plant defense , 2014, BMC Genomics.

[18]  David E Hill,et al.  Pooled‐matrix protein interaction screens using Barcode Fusion Genetics , 2016, Molecular systems biology.

[19]  S. Dudoit,et al.  Normalization of RNA-seq data using factor analysis of control genes or samples , 2014, Nature Biotechnology.

[20]  Nicolas Servant,et al.  A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis , 2013, Briefings Bioinform..

[21]  Santosh B. Satbhai,et al.  An extracellular network of Arabidopsis leucine-rich repeat receptor kinases , 2017, Nature.

[22]  Li Yang,et al.  Convergent targeting of a common host protein-network by pathogen effectors from three kingdoms of life. , 2014, Cell host & microbe.

[23]  Zhaojun Ding,et al.  A P-Loop NTPase Regulates Quiescent Center Cell Division and Distal Stem Cell Identity through the Regulation of ROS Homeostasis in Arabidopsis Root , 2016, PLoS genetics.

[24]  D. Nettleton,et al.  Correction to: small RNA discovery in the interaction between barley and the powdery mildew pathogen , 2019, BMC Genomics.

[25]  M. S. Mukhtar,et al.  Independently Evolved Virulence Effectors Converge onto Hubs in a Plant Immune System Network , 2011, Science.

[26]  Johanna Hardin,et al.  Selecting between‐sample RNA‐Seq normalization methods from the perspective of their assumptions , 2016, Briefings Bioinform..

[27]  Takaya Saito,et al.  The Precision-Recall Plot Is More Informative than the ROC Plot When Evaluating Binary Classifiers on Imbalanced Datasets , 2015, PloS one.

[28]  P. Bayrak-Toydemir,et al.  Hereditary hemorrhagic telangiectasia: genetics and molecular diagnostics in a new era , 2015, Front. Genet..

[29]  Joseph R. Ecker,et al.  CrY2H-seq: a massively-multiplexed assay for deep coverage interactome mapping , 2017, Nature Methods.

[30]  D. Tang,et al.  A mutation in the GTP hydrolysis site of Arabidopsis dynamin-related protein 1E confers enhanced cell death in response to powdery mildew infection. , 2006, The Plant journal : for cell and molecular biology.

[31]  S. Fields,et al.  The yeast two-hybrid assay: still finding connections after 25 years , 2014, Nature Methods.

[32]  Günter P. Wagner,et al.  Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples , 2012, Theory in Biosciences.

[33]  P. Schulze-Lefert,et al.  Barley MLA Immune Receptors Directly Interfere with Antagonistically Acting Transcription Factors to Initiate Disease Resistance Signaling[C][W] , 2013, Plant Cell.

[34]  Ulrich Stelzl,et al.  A Y2H-seq approach defines the human protein methyltransferase interactome , 2013, Nature Methods.

[35]  Tatsuya Akutsu,et al.  Determining the minimum number of protein-protein interactions required to support known protein complexes , 2018, PloS one.

[36]  M. Luo,et al.  The wheat Sr50 gene reveals rich diversity at a cereal disease resistance locus , 2015, Nature Plants.

[37]  D. Golenbock,et al.  The TLR4 adaptor TRAM controls the phagocytosis of Gram-negative bacteria by interacting with the Rab11-family interacting protein 2 , 2019, PLoS pathogens.

[38]  S. Lukyanov,et al.  Normalizing cDNA Libraries , 2010, Current protocols in molecular biology.

[39]  Uwe Schlattner,et al.  Yeast Two-Hybrid, a Powerful Tool for Systems Biology , 2009, International journal of molecular sciences.

[40]  E. Pohl,et al.  A unique dynamin-related protein is essential for mitochondrial fission in Toxoplasma gondii , 2019, PLoS pathogens.

[41]  Darrell Desveaux,et al.  Quantitative Interactor Screening with next-generation Sequencing (QIS-Seq) identifies Arabidopsis thaliana MLO2 as a target of the Pseudomonas syringae type III effector HopZ2 , 2012, BMC Genomics.

[42]  E. Lagudah,et al.  Cytosolic activation of cell death and stem rust resistance by cereal MLA-family CC–NLR proteins , 2016, Proceedings of the National Academy of Sciences.

[43]  Nischal Karki,et al.  The effect of light conditions on interpreting oil composition engineering in Arabidopsis seeds , 2018, Plant direct.

[44]  Bernhard Suter,et al.  Next-Generation Sequencing for Binary Protein–Protein Interactions , 2015, Front. Genet..

[45]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[46]  F. Martinon,et al.  Intracellular Trafficking of Interleukin-1 Receptor I Requires Tollip , 2006, Current Biology.

[47]  Jae-Seong Yang,et al.  rec-YnH enables simultaneous many-by-many detection of direct protein–protein and protein–RNA interactions , 2018, Nature Communications.

[48]  P. Boevink,et al.  Localization and domain characterization of Arabidopsis golgin candidates. , 2007, Journal of experimental botany.

[49]  M. Kabbage,et al.  An inhibitor of apoptosis (SfIAP) interacts with SQUAMOSA promoter‐binding protein (SBP) transcription factors that exhibit pro‐cell death characteristics , 2018, bioRxiv.

[50]  S. Cesari,et al.  Multiple strategies for pathogen perception by plant immune receptors. , 2018, The New phytologist.

[51]  M. Somerville.  A postmodern moral tale: the ethics of research relationships , 2002, Nature Reviews Drug Discovery.

[52]  M. Kemp,et al.  Oral, Nasal and Pharyngeal Exposure to Lipopolysaccharide Causes a Fetal Inflammatory Response in Sheep , 2015, PloS one.

[53]  Ming Chen,et al.  PRIN: a predicted rice interactome network , 2011, BMC Bioinformatics.

[54]  Jure Leskovec,et al.  Evolution of resilience in protein interactomes across the tree of life , 2018, Proceedings of the National Academy of Sciences.

[55]  M. Vidal,et al.  Identification of potential interaction networks using sequence-based searches for conserved protein-protein interactions or "interologs". , 2001, Genome research.

[56]  Javier De Las Rivas,et al.  Protein–Protein Interactions Essentials: Key Concepts to Building and Analyzing Interactome Networks , 2010, PLoS Comput. Biol..

[57]  P. Schulze-Lefert,et al.  Conservation of NLR-triggered immunity across plant lineages , 2012, Proceedings of the National Academy of Sciences.

[58]  P. Schulze-Lefert,et al.  Diversity at the Mla powdery mildew resistance locus from cultivated barley reveals sites of positive selection. , 2010, Molecular plant-microbe interactions : MPMI.

[59]  P. D. Nagy,et al.  The role of co-opted ESCRT proteins and lipid factors in protection of tombusviral double-stranded RNA replication intermediate against reconstituted RNAi in yeast , 2017, PLoS pathogens.

[60]  I. Somssich,et al.  Nuclear Activity of MLA Immune Receptors Links Isolate-Specific and Basal Disease-Resistance Responses , 2007, Science.

[61]  James E. Allen,et al.  Ensembl Genomes 2020—enabling non-vertebrate genomic research , 2019, Nucleic Acids Res..

[62]  D. Eeckhout,et al.  POLAR-guided signalling complex assembly and localization drive asymmetric cell division , 2018, Nature.

[63]  Patrick Breheny,et al.  DEEPN as an Approach for Batch Processing of Yeast 2-Hybrid Interactions. , 2016, Cell reports.

[64]  A. Hastie,et al.  Yeast two-hybrid interaction partner screening through in vivo Cre-mediated Binary Interaction Tag generation , 2007, Nucleic acids research.

[65]  I. Somssich,et al.  The wheat Mla homologue TmMla1 exhibits an evolutionarily conserved function against powdery mildew in both wheat and barley. , 2011, The Plant journal : for cell and molecular biology.

[66]  Qianhua Shen,et al.  An E3 Ligase Affects the NLR Receptor Stability and Immunity to Powdery Mildew1 , 2016, Plant Physiology.

[67]  S. Munro,et al.  Golgins , 2018, Current Biology.

[68]  Alain Goossens,et al.  A user-friendly platform for yeast two-hybrid library screening using next generation sequencing , 2018, PloS one.

[69]  Shili Lin,et al.  Rank aggregation methods , 2010 .

[70]  F. Ghadessy,et al.  Compartmentalized linkage of genes encoding interacting protein pairs , 2011, Proteomics.

[71]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[72]  S. Frank,et al.  The role of dynamin-related protein 1, a mediator of mitochondrial fission, in apoptosis. , 2001, Developmental cell.

[73]  D. Nettleton,et al.  Small RNA discovery in the interaction between barley and the powdery mildew pathogen , 2019, BMC Genomics.

[74]  Sankar Ghosh,et al.  Negative Regulation of Toll-like Receptor-mediated Signaling by Tollip* , 2002, The Journal of Biological Chemistry.

[75]  Tina T. Hu,et al.  A Genomic Map of the Effects of Linked Selection in Drosophila , 2014, PLoS genetics.

[76]  Mark Goadrich,et al.  The relationship between Precision-Recall and ROC curves , 2006, ICML.

[77]  Venkatramanan Krishnamani,et al.  MALTA: A calculator for estimating the coverage with shRNA, CRISPR, and cDNA libraries , 2019, SoftwareX.

[78]  A. Nakano,et al.  Arabidopsis dynamin-related proteins DRP3A and DRP3B are functionally redundant in mitochondrial fission, but have distinct roles in peroxisomal fission. , 2009, The Plant journal : for cell and molecular biology.

[79]  David Baker,et al.  High-throughput characterization of protein–protein interactions by reprogramming yeast mating , 2017, Proceedings of the National Academy of Sciences.