Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DNA sequencing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

Empirical evaluation of all unique Cas9 protospacers in E. coli reveal widespread functionality and rules for gRNA design

The Cas9 nuclease has become central to modern methods and technologies in synthetic biology, largely due to the ease with which it can be targeted to specific DNA loci via guide RNAs (gRNAs). Reports vary widely on the actual specificity of this targeting, with some studies observing 60% of gRNAs possessing no activity against the genome, yet an assumption persists within the E. coli community that inactive gRNAs are rare. To resolve these contradictions, we evaluated the activity of 463 000 unique gRNAs in the E. coli K12 MG1655 genome. We show that the overwhelming majority (at least 93%) of unique gRNAs are functional while only 0.3% are nonfunctional. These nonfunctional gRNAs exhibit strong spacer self-interaction, which can either be excluded using a simple design rule or “repaired” during library design. Finally, this work provides the greater microbial synthetic biology community both a set of nearly half a million empirically evaluated E. coli gRNAs as well as a thoroughly evaluated experimental procedure, complete with appropriate controls for Cas9 activity, for conducting Cas9 assays in E. coli specifically and bacteria more generally. Lastly, we have produced a webapp to allow users to easily browse and extract gRNA sequences from the E. coli genome, which can be accessed at https://grna.ornl.gov.

Kammerdiener, Elise K. [Oak Ridge National Laborat↗

The TLK-ASF1 histone chaperone pathway plays a critical role in IL-1β–mediated AML progression

Identifying and targeting microenvironment-driven pathways that are active across acute myeloid leukemia (AML) genetic subtypes should allow the development of more broadly effective therapies. The proinflammatory cytokine interleukin-1β (IL-1β) is abundant in the AML microenvironment and promotes leukemic growth. Through RNA-sequencing analysis, we identify that IL-1β–upregulated ASF1B (antisilencing function-1B), a histone chaperone, in AML progenitors compared with healthy progenitors. ASF1B, along with its paralogous protein ASF1A, recruits H3-H4 histones onto the replication fork during S-phase, a process regulated by Tousled-like kinase 1 and 2 (TLKs). Although ASF1s and TLKs are known to be overexpressed in multiple solid tumors and associated with poor prognosis, their functional roles in hematopoiesis and inflammation-driven leukemia remain unexplored. In this study, we identify that ASF1s and TLKs are overexpressed in multiple genetic subtypes of AML. We demonstrate that depletion of ASF1s significantly reduces leukemic cell growth in both in vitro and in vivo models using human cells. Using a murine model, we show that overexpression of ASF1B accelerates leukemia progression. Moreover, Asf1b or Tlk2 deletion delayed leukemia progression, whereas these proteins are dispensable for normal hematopoiesis. Through proteomics and phosphoproteomics analyses, we uncover that the TLK-ASF1 pathway promotes leukemogenesis by affecting the cell cycle and DNA damage pathways. Collectively, our findings identify the TLK1-ASF1 pathway as a novel mediator of inflammatory signaling and a promising therapeutic target for AML treatment across diverse genetic subtypes. Finally, selective inhibition of this pathway offers potential opportunities to intervene effectively, address intratumoral heterogeneity, and ultimately improve clinical outcomes in AML.

60 APPLIED LIFE SCIENCES↗

Phylogenomic Insights into the Evolution and Origin of Nematoda

Abstract The phylum Nematoda represents one of the most cosmopolitan and abundant metazoan groups on Earth. In this study, we reconstructed the phylogenomic tree for phylum Nematoda. A total of 60 genomes, belonging to 8 nematode orders, were newly sequenced, providing the first low-coverage genomes for the orders Dorylaimida, Mononchida, Monhysterida, Chromadorida, Triplonchida, and Enoplida. The resulting phylogeny is well-resolved across most clades, with topologies remaining consistent across various reconstruction parameters. The subclass Enoplia is placed as a sister group to the rest of Nematoda, agreeing with previously published phylogenies. While the order Triplonchida is monophyletic, it is not well-supported, and the order Enoplida is paraphyletic. Taxa possessing a stomatostylet form a monophyletic group; however, the superfamily Aphelenchoidea does not constitute a monophyletic clade. The genera Trichinella and Trichuris are inferred to have shared a common ancestor approximately 202 millions of years ago (Ma), a considerably later period than previously suggested. All stomatostylet-bearing nematodes are proposed to have originated ~305 Ma, corresponding to the transition from the Devonian to the Permian period. The genus Thornia is placed outside of Dorylaimina and Nygolaimina, disagreeing with its position in previous studies. In addition, we tested the whole genome amplification method and demonstrated that it is a promising strategy for obtaining sufficient DNA for phylogenomic studies of microscopic eukaryotes. This study significantly expanded the current nematode genome dataset, and the well-resolved phylogeny enhances our understanding of the evolution of Nematoda.

Qing, Xue (ORCID:0000000203559956)↗

gRNA-SeqRET: a universal tool for targeted and genome-scale gRNA design and sequence extraction for prokaryotes and eukaryotes

High-throughput genetic screening is frequently employed to rapidly associate gene with phenotype and establish sequence-function relationships. With the advent of CRISPR technology, and the ability to functionally interrogate previously genetically recalcitrant organisms, non-model organisms can be investigated using pooled guide RNA (gRNA) libraries and sequencing-based assays to quantitatively assess fitness of every targeted locus in parallel. To aid the construction of pooled gRNA assemblies, we have developed an in silico design workflow for gRNA selection using the gRNA Sequence Region Extraction Tool (gRNA-SeqRET). Built upon the previously developed CCTop, gRNA-SeqRET enables automated, scalable design of gRNA libraries that target user-specified regions or whole genomes of any prokaryote or eukaryote. Additionally, gRNA-SeqRET automates the bulk extraction of any regions of sequence relative to genes or other features, aiding in the design of homology arms for insertion or deletion constructs. We also assess in silico the application of a designed gRNA library to other closely related genomes and demonstrate that for very closely related organisms Average Nucleotide Identity (ANI) > 95% a large fraction of the library may be of relevance. The gRNA-SeqRET web application pipeline can be accessed at https://grna.jgi.doe.gov. The source code is comprised of freely available software tools and customized Python scripts, and is available at https://bitbucket.org/berkeleylab/grnadesigner/src/master/ under a modified BSD open-source license (https://bitbucket.org/berkeleylab/grnadesigner).

59 BASIC BIOLOGICAL SCIENCES↗

Isolation of Histone from Sorghum Leaf Tissue for Top Down Mass Spectrometry Profiling of Potential Epigenetic Markers

Histones belong to a family of highly conserved proteins in eukaryotes. They pack DNA into nucleosomes as functional units of chromatin. Post-translational modifications (PTMs) of histones, which are highly dynamic and can be added or removed by enzymes, play critical roles in regulating gene expression. In plants, epigenetic factors including histone PTMs are related to their adaptive responses to the environment. Understanding the molecular mechanisms of epigenetic control can bring unprecedented opportunities for engineering solutions to increase the resilience of crops to climate change. Herein, we describe a protocol to isolate the nuclei and purify histones from sorghum leaf tissue. The extracted histones can be analyzed as their intact forms by top-down mass spectrometry (MS) coupled to online reversed-phase (RP) liquid chromatography (LC). Combinations and stoichiometry of multiple PTMs on the same histone proteoform can be readily identified. In addition, histone tail clipping can be detected using the top-down LC-MS workflow thus yielding the global PTM profile of core histones (H4, H2A, H2B, H3). By comparing the PTM profile among samples corresponding to different conditions (e.g. drought vs. control), potential epigenetic marks can be discovered as targets for further characterization using approaches such as chromatin immunoprecipitation – sequencing (ChIP-seq).

59 BASIC BIOLOGICAL SCIENCES↗

Effects of warming on bacterial growth rates in a peat soil under ambient and elevated CO 2

Boreal peatlands are important global carbon reservoirs that are particularly vulnerable to predicted climate changes such as increasing CO 2 and temperature. Since microbial activities regulate the balance of carbon sequestered into soil organic matter or remineralized to CO 2 , characterizing their response to these environmental factors is critical to predicting how peatland ecosystems will affect climate-carbon cycle feedbacks. Here we examined in-situ taxon-specific variation in microbial growth under long-term elevated CO 2 and across a gradient of warming treatments in a northern Minnesota peat bog using quantitative stable isotope probing with 18 O-water. Across temperatures, bacterial taxa were grouped according to the excess atom fraction 18 O (EAF) of their genomes, a proxy for DNA replication and hence, growth. Taxon-specific growth across CO 2 and temperature treatments clustered into relatively few response patterns. While a large portion of taxa showed little to no growth under ambient CO 2 , many of the same taxa grew rapidly under elevated CO 2 . We found support for phylogenetic conservation of response patterns among Acidobacteria and Proteobacteria, the two most abundant phyla in our data. Our results suggest certain taxa may be primed for new climate conditions and have a greater influence on carbon cycling with implications for future climate mitigation strategies.

16S amplicon sequencing, Carbon Dioxide (CO2), pea↗

Molecular Insight into the β-Sheet Twist and Related Morphology of Self-Assembled Peptide Amphiphile Ribbons

Self-assembly of high-aspect-ratio filaments containing β-sheets has attracted much attention due to potential use in bioengineering and biomedicine. However, precisely predicting the assembled morphologies remains a grand challenge because of insufficient understanding of the self-assembly process. We employed an atomistic model to study the self-assembly of peptide amphiphiles (PAs) containing valine–glutamic acid (VE) dimeric repeats. By changing of the sequence length, the assembly morphology changes from flat ribbon to left-handed twisted ribbon, implying a relationship between β-sheet twist and strength of interstrand hydrogen bonds. The calculations are used to quantify this relationship including both magnitude and sign of the ribbon twist angle. Interestingly, a change in chirality is observed when we introduce the RGD epitope into the C-terminal of VE repeats, suggesting arginine and glycine’s role in suppressing right-handed β-sheet formation. Furthermore, this study provides insight into the relationship between β-sheet twist and self-assembled nanostructures including a possible design rule for PA self-assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Respiration data, microbial community assembly data, and FTICR-MS data associated with: “Disturbance Triggers Non-Linear Microbe-Environment Feedbacks. Sengupta et al., 2021, Biogeosciences”

This data package is associated with the manuscript “Disturbance Triggers Non-Linear Microbe-Environment Feedbacks, in revision in Biogeosciences (Sengupta et al. and 2021;https://bg.copernicus.org/preprints/bg-2021-51/). The study used hyporheic zone sediments as a model system to provide an integrated view of how disturbance modulates linkages among microbial ecology, biogeochemistry, and organic matter thermodynamics. Laboratory experiments exposed hyporheic sediment to varying wetting/drying dynamics. Data types include dissolved oxygen rates used to derive respiration rates, Fourier-transform ion cyclotron resonance mass spectrometry (FTICR-MS) data used to derive thermodynamic properties of organic matter, and microbial community assembly metrics derived from amplicon-sequence data of putatively active (cDNA) and whole community (gDNA). The outcomes of the study are condensed into a broadly applicable conceptual model linking external forcing, internal dynamics, and history. This data package is comprised of a file-level metadata (FLMD) csv, metadata csv, and seven folders that contain csv files, R scripts, xml files, and associated documentation: (1) Rates, (2) bNTI, and (3) FTICR, (4) Statistics_Analyses, (5) Raw OTU Beta dispersion Analysis, (6) bMNTD Randomizations, and (7) Data Dictionaries. The FLMD file has a description of each file included in the data package.

54 ENVIRONMENTAL SCIENCES↗

Data for An Orphan Gene BOOSTER Enhances Photosynthetic Efficiency and Plant Productivity

Seeds of Col-0 wild type, sig6 T-DNA mutants (CS877785, ABRC), PRL-1-OE, and sig6 T-DNA mutants transfected with PRL-1 (sig6::PRL-1) were planted in 1/2 MS media. Seedlings growth including chlorophyll development defects were investigated across the genotypes. Four-days-old-post-light exposure seedlings were harvested and performed RNAseq analysis with four biological replicates.

Biomass Analytics↗

Chromatin structure predicts survival in glioma patients

The pathological changes in epigenetics and gene regulation that accompany the progression of low-grade to high-grade gliomas are under-studied. The authors use a large set of paired atac-seq and RNA-seq data from surgically resected glioma specimens to infer gene regulatory relationships in glioma. Thirty-eight glioma patient samples underwent atac-seq sequencing and 16 samples underwent additional RNA-seq analysis. Using an atac-seq/RNA-seq correlation matrix, atac-seq peaks were paired with genes based on high correlation values (|r 2 | > 0.6). Samples clustered by IDH1 status but not by grade. Surprisingly there was a trend for IDH1 mutant samples to have more peaks. The majority of peaks are positively correlated with survival and positively correlated with gene expression. Constructing a model of the top six atac-seq peaks created a highly accurate survival prediction model (r 2 = 0.68). Four of these peaks were still significant after controlling for age, grade, pathology, IDH1 status and gender. Grade II, III, and IV (primary) samples have similar transcription factors and gene modules. However, grade IV (recurrent) samples have strikingly few peaks. Patient-derived glioma cultures showed decreased peak counts following radiation indicating that this may be radiation-induced. This study supports the notion that IDH1 mutant and IDH1 wildtype gliomas have different epigenetic landscapes and that accessible chromatin sites mapped by atac-seq peaks tend to be positively correlated with expression. The data in this study leads to a new model of treatment response wherein glioma cells respond to radiation therapy by closing open regions of DNA.

59 BASIC BIOLOGICAL SCIENCES↗

Pooled PPIseq: Screening the SARS-CoV-2 and human interface with a scalable multiplexed protein-protein interaction assay platform

Protein-Protein Interactions (PPIs) are a key interface between virus and host, and these interactions are important to both viral reprogramming of the host and to host restriction of viral infection. In particular, viral-host PPI networks can be used to further our understanding of the molecular mechanisms of tissue specificity, host range, and virulence. At higher scales, viral-host PPI screening could also be used to screen for small-molecule antivirals that interfere with essential viral-host interactions, or to explore how the PPI networks between interacting viral and host genomes co-evolve. Current high-throughput PPI assays have screened entire viral-host PPI networks. However, these studies are time consuming, often require specialized equipment, and are difficult to further scale. Here, we develop methods that make larger-scale viral-host PPI screening more accessible. This approach combines the mDHFR split-tag reporter with the iSeq2 interaction-barcoding system to permit massively-multiplexed PPI quantification by simple pooled engineering of barcoded constructs, integration of these constructs into budding yeast, and fitness measurements by pooled cell competitions and barcode-sequencing. We applied this method to screen for PPIs between SARS-CoV-2 proteins and human proteins, screening in triplicate >180,000 ORF-ORF combinations represented by >1,000,000 barcoded lineages. Our results complement previous screens by identifying 74 putative PPIs, including interactions between ORF7A with the taste receptors TAS2R41 and TAS2R7, and between NSP4 with the transmembrane KDELR2 and KDELR3. We show that this PPI screening method is highly scalable, enabling larger studies aimed at generating a broad understanding of how viral effector proteins converge on cellular targets to effect replication.

60 APPLIED LIFE SCIENCES↗

Intestinal microbiome profile of the brown rock sea cucumber ( Holothuria glaberrima ) using ITS and 16S rDNA amplicons from direct mechanical, enzymatic, and chemical metagenomic extraction

Using direct mechanical, enzymatic, and chemical extraction methods, the intestinal microbiome of the marine invertebrate Holothuria glaberrima was obtained. ITS and 16S rDNA regions were sequenced to enrich and investigate the prokaryotic and fungal diversity profiles from different anatomical regions within the sea cucumber’s intestinal biology.

16S↗

Essential amino acids in the Plant-Conserved and Class-Specific Regions of cellulose synthases

The Plant-Conserved Region (P-CR) and the Class-Specific Region (CSR) are two plant-unique sequences in the catalytic core of cellulose synthases (CESAs) for which specific functions have not been established. Here, we used site-directed mutagenesis to replace amino acids and motifs within these sequences predicted to be essential for assembly and function of CESAs. We developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis (Arabidopsis thaliana) temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the CSR, which blocks dimerization in vitro, rendered the AtCesA1 transgene unable to complement the rsw1 mutation. Examination of the CSR sequences from 33 diverse angiosperm species showed domains of high-sequence conservation in a class-specific manner but with variation in the degrees of disorder, indicating a nonredundant role of the CSR structures in different CESA isoform classes. The Cys residue essential for dimerization was not always located in domains of intrinsic disorder. Expression of AtCesA1 transgene constructs, in which Pro417 and Arg453 were substituted for Ala or Lys in the coiled-coil of the P-CR, were also unable to complement the rsw1 mutation. Despite an expected role for Arg457 in trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Our data support that Cys662 within the CSR and Pro417 and Arg453 within the P-CR of Arabidopsis CESA1 are essential residues for functional synthase complex formation, but our data do not support a specific role for Arg457 in trimerization in native CESA complexes.

59 BASIC BIOLOGICAL SCIENCES↗

Hyporheic zone, river, and groundwater metagenome resolved genomes and rpS3 genes in East River Watershed, Colorado USA Summer 2020, 2021

Here we present metagenome assembled genomes (MAGs) for the bacterial and archaeal communities from water filter collected across 8 locations along the East River Watershed, CO, and 1 nearby groundwater well. The purpose was to look for connectivity and similarities across the network and to see the impact of the groundwater. As a part of Lawrence Berkeley National Laboratory (LBNL) Watershed Science Focus Area (SFA), we assessed community composition and strain similarities between the sites and we also compared it to previous metagenomic studies within the watershed looking at floodplain (Matheus Carnevali et al. 2021) and hillslope (Lavy et al. 2019) microbiomes. Here we present metagenome assembled genomes (MAGs) for the bacterial and archaeal communities from filters across 8 locations during August 2020 and July 2021. This resulted in 32 samples. The groundwater sample was sequenced at UC Berkley's QB3. The other 31 samples were sequenced at University of Maryland. Metagenomes were assembled using four autobinners and the best bins were selected using dasTool. The genomes were dereplicated at 95% with dRep and the subset of winning genomes were manually curated based on visual inspection of taxonomic profile, GC content, coverage, and a set of 51 bacterial single copy genes (BSCG), and 38 archaeal signal copy genes (ASCG). The dataset includes a zip file of 311 genomes (HZ_River_SW_MAGS_Dereplicated_95.zip). The dataset additionally includes a zipped file of ribosomal protein small subunit 3 (rpS3) proteins from the hyporheic zone and river data (rpS3_Proteins_HZ_River.zip), a metadata file used to register associated samples with IGSNs (International Generic Sample Numbers) (samples.csv), a location metadata file (locations.csv). This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

DNA↗

Integration of ultra-low coverage whole-genome sequences for reconstructing the evolutionary history of Galapagos giant tortoises

Genomic data from contemporary and historical samples often need to be coupled for evolutionary reconstructions of multitaxon complexes. However, the genetic data recovered from historical samples may result only in ultra-low coverage whole-genome sequences (ulcWGS; <0.15× depth), leading to inaccurate evolutionary inferences given a preponderance of missing data. Using the Galapagos giant tortoise radiation as a study system (Chelonoidis spp., composed of 13 extant and four extinct lineages), we assembled a novel methodological pipeline that removes potential noise introduced by the missing data and enhances the evolutionary signal from ulcWGS samples. We leveraged existing tools for phylogenomic placement (EPA-ng), population genomic structure (smartsnp) and admixture (Admixfrog, NGSadmix) to demonstrate that the evolutionary history of samples can be uncovered with sequencing depths as low as 0.008–0.139×. Importantly, these approaches do not use genotype imputation of the ulcWGS samples, which would require extensive reference datasets. Our application to two cases of extinct lineages of Galapagos giant tortoises, with and without references from the same lineage, demonstrates the general value of the approach. We confirm where the extinct lineages from San Cristóbal and Santa Fe islands fit into the Galapagos giant tortoise radiation, and that these lineages were evolutionarily distinct entities.

ancient DNA↗

Montane Conifer, Aspen, Meadow, and Sagebrush Metagenome Resolved Genomes and Traits in East River Watershed, Colorado, USA

Climate change is driving vegetation shifts in mountain watersheds, with unknown impacts on biogeochemical cycles. We hypothesize that these shifts will reshape soil microbiomes and associated biogeochemical processes. As a part of Lawrence Berkeley National Laboratory (LBNL) Watershed Science Focus Area (SFA), we assessed microbiome and microbial functional trait differences between soils under conifer, aspen, forby meadows, and sagebrush across the East River Watershed, CO, controlling for elevation and aspect.Here we present metagenome assembled genomes (MAGs) for the bacterial and archaeal communities from soils 0-20cm in depth across three locations in the watershed—Headwaters, Upper Reaches, and Lower Reaches from August 3-11th 2016. Each location was further subdivided into two blocks, with one block on a west facing aspect, and two on the east aspect of the valley. Within blocks, two samples per vegetation type were taken (one at each depth). This resulted in 66 samples, which were sequenced at JGI and can be found under the Joint Genome Institute (JGI) Genomes Online Database (GOLD) sequencing project Gs0118068. Metagenomes were assembled through an inhouse pipeline (see methods), binned using four autobinners (concoct, maxbin2, metabat2, and vamb) and consolidated using dastool. The consolidated bins from all metagenomes were pooled, filtered by completeness (>75%) and contamination (<25%), and dereplicated at 95% ANI using drep. The dataset includes a zip file of 687 genomes (Vegtype_MAGS.zip), the accession numbers for the underlying metagenomes, a csv file with MAG quality metrics and taxonomy from Genome Taxonomy Database (GTDB) and National Center for Biotechnology Information (NCBI) taxonomic representative genome proteins (EastRiver_Vegtype_drep_genome_info.csv), and a file containing MAG quality metrics and taxonomy (gtdb_drep_bin_taxonomy.csv). The dataset additionally includes a sample metadata file (EastRiver_Vegtype_sample_metadata.csv), a metadata file used to register associated samples with IGSNs (International Generic Sample Numbers) (samples.csv), a Google KML file for the sampled locations (sample_collection_sites.kml), a location metadata file (locations.csv), a file-level metadata file (flmd.csv), and a data dictionary (dd.csv) file.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Development, optimization, and application of an episomal plasmid system for Rhodotorula toruloides

Rhodotorula toruloides is an emerging oleaginous yeast with strong potential as a microbial cell factory for the production of acetyl-CoA-derived bioproducts. However, engineering of this organism has been limited by the absence of a functional episomal plasmid system, a foundational genetic tool for rapid gene expression, pathway testing, and CRISPR-based genome engineering. Here, we report the first episomal plasmid system for R. toruloides . Through systematic screening of candidate autonomously replicating sequences (ARSs) from diverse sources, we identified multiple functional ARS elements and selected C63F4, a fragment derived from Contig 63 of R. toruloides CBS14, because of its stable performance. The resulting pC63F4 plasmid was maintained episomally, supported GFP reporter expression, exhibited a copy number of 2.39 ± 0.13, and showed good stability during long term cultivation. To overcome poor transformation efficiency, we developed a Cre- loxP -mediated in vivo re-circularization strategy that enabled reliable delivery of the episomal plasmid. Using this improved system, we demonstrated functional episomal expression of metabolic engineering genes and multi-gene pathways for the production of triacetic acid lactone, fatty alcohols, and limonene. Finally, we leveraged this platform to establish a redesigned CRISPR system that enables seamless genome editing in R. toruloides for the first time, while also simplifying marker recycling. Together, this work establishes a long-needed episomal plasmid platform and associated CRISPR toolkit that will accelerate metabolic engineering, synthetic biology, and fundamental studies in R. toruloides .

CRISPR-Cas9↗

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.

59 BASIC BIOLOGICAL SCIENCES↗