High-throughput Functional Genomics Research Articles

A benchmarked, high-efficiency prime editing platform for multiplexed dropout screening.

Prime editing installs precise edits into the genome with minimal unwanted byproducts, but low and variable editing efficiencies have complicated application of the approach to high-throughput functional genomics. Here we assembled a prime editing platform capable of high-efficiency substitution editing suitable for functional interrogation of small genetic variants. We benchmarked this platform for pooled, loss-of-function screening using a library of ~240,000 engineered prime editing guide RNAs (epegRNAs) targeting ~17,000 codons with 1-3 bp substitutions. Comparing the abundance of these epegRNAs across screen samples identified negative selection phenotypes for 7,996 nonsense mutations targeted to 1,149 essential genes and for synonymous mutations that disrupted splice site motifs at 3' exon boundaries. Rigorous evaluation of codon-matched controls demonstrated that these phenotypes were highly specific to the intended edit. Altogether, we established a prime editing approach for multiplexed, functional characterization of genetic variants with simple readouts.

Fine-tuning of a CRISPRi screen in the seventh pandemic Vibrio cholerae

BackgroundVibrio cholerae O1 El Tor, the etiological agent responsible for the last cholera pandemic, has become a well-established model organism for which some genetic tools are available. While CRISPRi technology has been applied to V. cholerae, improvements were necessary to upscale it and enable pooled screening by high-throughput sequencing in this bacterium.ResultsIn this study, we present a genome-wide CRISPR-dCas9 screen specifically optimized for the N16961 El Tor model strain of V. cholerae. This approach is characterized by a tight control of dCas9 expression and activity, as well as a streamlined experimental setup. Our library allows the depletion of 3,674 (98.9%) annotated genes from the V. cholerae genome. To confirm its effectiveness, we screened for genes that are essential during exponential growth in rich medium and identified 369 genes for which guides were significantly depleted from the library (log2FC < -2). Remarkably, 82% of these genes had previously been described as hypothetical essential genes in V. cholerae or in a closely related bacterium, V. natriegens.ConclusionWe thus validated the robustness and accuracy of our CRISPRi-based approach for assessing gene fitness in a given condition. Our findings highlight the efficacy of the developed CRISPRi platform as a powerful tool for high-throughput functional genomics studies of V. cholerae.

High-throughput fitness experiments reveal specific vulnerabilities of human-adapted Salmonella during stress and infection

Salmonella enterica is comprised of genetically distinct ‘serovars’ that together provide an intriguing model for exploring the genetic basis of pathogen evolution. Although the genomes of numerous Salmonella isolates with broad variations in host range and human disease manifestations have been sequenced, the functional links between genetic and phenotypic differences among these serovars remain poorly understood. Here, we conduct high-throughput functional genomics on both generalist (Typhimurium) and human-restricted (Typhi and Paratyphi A) Salmonella at unprecedented scale in the study of this enteric pathogen. Using a comprehensive systems biology approach, we identify gene networks with serovar-specific fitness effects across 25 host-associated stresses encountered at key stages of human infection. By experimentally perturbing these networks, we characterize previously undescribed pseudogenes in human-adapted Salmonella. Overall, this work highlights specific vulnerabilities encoded within human-restricted Salmonella that are linked to the degradation of their genomes, shedding light into the evolution of this enteric pathogen.

Genetics and Glaucoma: the state of the art.

Glaucoma is the second leading cause of irreversible blindness worldwide. Although genetic background contributes differently to rare early-onset glaucoma (before age 40) or common adult-onset glaucoma, it is now considered an important factor in all major forms of the disease. Genetic and genomic studies, including GWAS, are contributing to identifying novel loci associated with glaucoma or to endophenotypes across ancestries to enrich the knowledge about glaucoma genetic susceptibility. Moreover, new high-throughput functional genomics contributes to defining the relevance of genetic results in the biological pathways and processes involved in glaucoma pathogenesis. Such studies are expected to advance significantly our understanding of glaucoma's genetic basis and provide new druggable targets to treat glaucoma. This review gives an overview of the role of genetics in the pathogenesis or risk of glaucoma.

GARP and EARP are required for efficient BoHV-1 replication as identified by a genome wide CRISPR knockout screen.

The advances in gene editing bring unprecedented opportunities in high throughput functional genomics to animal research. Here we describe a genome wide CRISPR knockout library, btCRISPRko.v1, targeting all protein coding genes in the cattle genome. Using it, we conducted genome wide screens during Bovine Herpes Virus type 1 (BoHV-1) replication and compiled a list of pro-viral and anti-viral candidates. These candidates might influence multiple aspects of BoHV-1 biology such as viral entry, genome replication and transcription, viral protein trafficking and virion maturation in the cytoplasm. Some of the most intriguing examples are VPS51, VPS52 and VPS53 that code for subunits of two membrane tethering complexes, the endosome-associated recycling protein (EARP) complex and the Golgi-associated retrograde protein (GARP) complex. These complexes mediate endosomal recycling and retrograde trafficking to the trans Golgi Network (TGN). Simultaneous loss of both complexes in MDBKs resulted in greatly reduced production of infectious BoHV-1 virions. We also found that viruses released by these deficient cells severely lack VP8, the most abundant tegument protein of BoHV-1 that are crucial for its virulence. In combination with previous reports, our data suggest vital roles GARP and EARP play during viral protein packaging and capsid re-envelopment in the cytoplasm. It also contributes to evidence that both the TGN and the recycling endosomes are recruited in this process, mediated by these complexes. The btCRISPRko.v1 library generated here has been controlled for quality and shown to be effective in host gene discovery. We hope it will facilitate efforts in the study of other pathogens and various aspects of cell biology in cattle.

Disentangling the complexity of psoriasis in the post-genome-wide association era.

In recent years, genome-wide association studies (GWAS) have been instrumental in unraveling the genetic architecture of complex diseases, including psoriasis. The application of large-scale GWA studies in psoriasis has illustrated several associated loci that participate in the cutaneous inflammation, however explaining a fraction of the disease heritability. With the advent of high-throughput sequencing technologies and functional genomics approaches, the post-GWAS era aims to unravel the functional mechanisms underlying the inter-individual variability in psoriasis patients. In this review, we present the key advances of psoriasis GWAS in under-represented populations, rare, non-coding and structural variants and epistatic phenomena that orchestrate the interplay between different cell types. We further review the gene-gene and gene-environment interactions contributing to the disease predisposition and development of comorbidities through Mendelian randomization studies and pleiotropic effects of psoriasis-associated loci. We finally examine the holistic approaches conducted in psoriasis through system genetics and state-of-the-art transcriptomic analyses, discussing their potential implication in the expanding field of precision medicine and characterization of comorbidities.

High-throughput functional genomics: A (myco)bacterial perspective.

Advances in sequencing technologies have enabled unprecedented insights into bacterial genome composition and dynamics. However, the disconnect between the rapid acquisition of genomic data and the (much slower) confirmation of inferred genetic function threatens to widen unless techniques for fast, high-throughput functional validation can be applied at scale. This applies equally to Mycobacterium tuberculosis, the leading infectious cause of death globally and a pathogen whose genome, despite being among the first to be sequenced two decades ago, still contains many genes of unknown function. Here, we summarize the evolution of bacterial high-throughput functional genomics, focusing primarily on transposon (Tn)-based mutagenesis and the construction of arrayed mutant libraries in diverse bacterial systems. We also consider the contributions of CRISPR interference as a transformative technique for probing bacterial gene function at scale. Throughout, we situate our analysis within the context of functional genomics of mycobacteria, focusing specifically on the potential to yield insights into M. tuberculosis pathogenicity and vulnerabilities for new drug and regimen development. Finally, we offer suggestions for future approaches that might be usefully applied in elucidating the complex cellular biology of this major human pathogen.

Genome-scale CRISPR screening in a single mouse liver

A complete understanding of the genetic determinants underlying mammalian physiology and disease is limited by the capacity for high-throughput genetic dissection in the living organism. Genome-wide CRISPR screening is a powerful method for uncovering the genetic regulation of cellular processes, but the need to stably deliver single guide RNAs to millions of cells has largely restricted its implementation to ex vivo systems. There thus remains a need for accessible high-throughput functional genomics in vivo. Here, we establish genome-wide screening in the liver of a single mouse and use this approach to uncover regulation of hepatocyte fitness. We uncover pathways not identified in cell culture screens, underscoring the power of genetic dissection in the organism. The approach we developed is accessible, scalable, and adaptable to diverse phenotypes and applications. We have hereby established a foundation for high-throughput functional genomics in a living mammal, enabling comprehensive investigation of physiology and disease.

LGG-59. Identifying hidden drivers of low-grade glioma tumor growth

Genomic drivers of pediatric low-grade gliomas (pLGGs) converge on alterations that activate the MAPK pathway. However, expression of individual driver oncogenes fails to induce tumor formation with high penetrance and, paradoxically, expression of these oncogenes suppresses growth in vitro. This, combined with the non-monotonic tumor growth rate in patients, suggests that there are “hidden drivers” beyond a single driver oncogene that are necessary to support tumor growth. The goal of this project is to leverage high-throughput functional genomics strategies to identify these hidden drivers of pLGG. Our preliminary data indicates that genes which modulate differentiation are required for the survival of LGG cells, suggesting that these genes may be hidden drivers of LGG tumor growth. Additionally, we hypothesize that secreted factors in the tumor microenvironment regulate pLGG tumor growth, potentially by modulating differentiation. In total, genes which cooperate with pLGG driver oncogenes to promote tumor growth may represent a new class of therapeutic targets and may explain the complex patterns of tumor growth that are observed in patients.

Pediatric Sarcomas: The Next Generation of Molecular Studies.

Simple SummaryThere has been an incredible amount of discovery in pediatric sarcomas, but much remains to be accomplished. Clinical challenges include diagnostic heterogeneity and the poor outcome of patients with high risk, metastatic, and relapsed disease. The emergence of single cell sequencing has allowed the ability to document tumor cell heterogeneity in amazing detail, but it does not allow the ability to visualize spatial orientation. This problem has been solved by spatial multi-omics, which can be used to map tumors and visualize the distribution of critical transcripts, mutations, and proteins. However, these tools only offer observational data. High-throughput functional genomics provides a powerful way to highlight oncogenic drivers and potential therapy opportunities. Research has been hamstrung by a need for annotated specimens, particularly in post-therapy, relapsed, and metastatic disease, and initial biopsies offer only limited data opportunities. Data complexity, variability, and inconsistency present problems best approached with AI/machine learning. We stand on the threshold of a revolution in cancer cell biology that has the potential for translation into more effective and more directed therapies, particularly for previously recalcitrant diseases.Pediatric sarcomas constitute one of the largest groups of childhood cancers, following hematopoietic, neural, and renal lesions. Partly because of their diversity, they continue to offer challenges in diagnosis and treatment. In spite of the diagnostic, nosologic, and therapeutic gains made with genetic technology, newer means for investigation are needed. This article reviews emerging technology being used to study human neoplasia and how these methods might be applicable to pediatric sarcomas. Methods reviewed include single cell RNA sequencing (scRNAseq), spatial multi-omics, high-throughput functional genomics, and clustered regularly interspersed short palindromic sequence-Cas9 (CRISPR-Cas9) technology. In spite of these advances, the field continues to be challenged by a dearth of properly annotated materials, particularly from recurrences and metastases and pre- and post-treatment samples.

Virus-Induced Gene Silencing in Wheat and Related Monocot Species.

Advances made in genome sequencing projects and structural genomics are generating large repertoire of candidate genes in plants associated with specific agronomic traits. Rapid and high-throughput functional genomics approaches are therefore needed to validate the biological function of these genes especially for agronomically important crops beyond the few model plant species. This can be achieved by utilizing available gene knockout or transgenic methodologies, but these can take considerable time and effort particularly in crops with large and complex genomes such as wheat. Therefore, any tool that expedites the validation of gene function is of particular benefit especially in cereal crop plants that are genetically difficult to transform. One such reverse genetics tool is virus-induced gene silencing (VIGS) which relies on the plants' natural antiviral RNA silencing defence mechanism. VIGS is used to downregulate target gene expression in a transient manner which persists long enough to determine its effect on a specific trait. VIGS based on Barley stripe mosaic virus (BSMV) is rapid, powerful, efficient, and relatively inexpensive tool for the analysis of gene function in cereal species. Here we present detailed protocols for BSMV-mediated VIGS for robust gene silencing in bread wheat and related species.

Capillary Electrophoresis-Based Functional Genomics Screening to Discover Novel Archaeal DNA Modifying Enzymes.

ABSTRACTIt has been predicted that 30 to 80% of archaeal genomes remain annotated as hypothetical proteins with no assigned gene function. Further, many archaeal organisms are difficult to grow or are unculturable. To overcome these technical and experimental hurdles, we developed a high-throughput functional genomics screen that utilizes capillary electrophoresis (CE) to identify nucleic acid modifying enzymes based on activity rather than sequence homology. Here, we describe a functional genomics screening workflow to find DNA modifying enzyme activities encoded by the hyperthermophile Thermococcus kodakarensis (T. kodakarensis). Large DNA insert fosmid libraries representing an ∼5-fold average coverage of the T. kodakarensis genome were prepared in Escherichia coli. RNA-seq showed a high fraction (84%) of T. kodakarensis genes were transcribed in E. coli despite differences in promoter structure and translational machinery. Our high-throughput screening workflow used fluorescently labeled DNA substrates directly in heat-treated lysates of fosmid clones with capillary electrophoresis detection of reaction products. Using this method, we identified both a new DNA endonuclease activity for a previously described RNA endonuclease (Nob1) and a novel AP lyase DNA repair enzyme family (termed 'TK0353') that is found only in a small subset of Thermococcales. The screening methodology described provides a fast and efficient way to explore the T. kodakarensis genome for a variety of nucleic acid modifying activities and may have implications for similar exploration of enzymes and pathways that underlie core cellular processes in other Archaea.IMPORTANCE This study provides a rapid, simple, high-throughput method to discover novel archaeal nucleic acid modifying enzymes by utilizing a fosmid genomic library, next-generation sequencing, and capillary electrophoresis. The method described here provides the details necessary to create 384-well fosmid library plates from Thermococcus kodakarensis genomic DNA, sequence 384-well fosmids plates using Illumina next-generation sequencing, and perform high-throughput functional read-out assays using capillary electrophoresis to identify a variety of nucleic acid modifying activities, including DNA cleavage and ligation. We used this approach to identify a new DNA endonuclease activity for a previously described RNA endonuclease (Nob1) and identify a novel AP lyase enzyme (TK0353) that lacks sequence homology to known nucleic acid modifying enzymes.

SorTn-seq: a high-throughput functional genomics approach to discovering regulators of bacterial gene expression

We recently developed a high-throughput functional genomics approach, named 'SorTn-seq', to identify factors affecting expression of any gene of interest in bacteria. Our approach facilitates high-throughput screening of complex mutant pools, a task previously hindered by a lack of suitable techniques. SorTn-seq combines high-density, Tn5-like transposon mutagenesis with fluorescence-activated cell sorting of a strain harboring a promoter-fluorescent reporter fusion, to isolate mutants with altered gene expression. The transposon mutant pool is sorted into different bins on the basis of fluorescence, and mutants are deep-sequenced to identify transposon insertions. DNA is prepared for sequencing by using commercial kits augmented with custom primers, enhancing ease of use and reproducibility. Putative regulators are identified by comparing the number of insertions per genomic feature in the different sort bins, by using existing bioinformatic pipelines and software packages. SorTn-seq can be completed in 1-2 weeks and requires general microbiology skills and basic flow cytometry experience.

Multiplexed functional genomic analysis of 5\u2019 untranslated region mutations across the spectrum of prostate cancer

The functional consequences of genetic variants within 5’ untranslated regions (UTRs) on a genome-wide scale are poorly understood in disease. Here we develop a high-throughput multi-layer functional genomics method called PLUMAGE (Pooled full-length UTR Multiplex Assay on Gene Expression) to quantify the molecular consequences of somatic 5’ UTR mutations in human prostate cancer. We show that 5’ UTR mutations can control transcript levels and mRNA translation rates through the creation of DNA binding elements or RNA-based cis-regulatory motifs. We discover that point mutations can simultaneously impact transcript and translation levels of the same gene. We provide evidence that functional 5’ UTR mutations in the MAP kinase signaling pathway can upregulate pathway-specific gene expression and are associated with clinical outcomes. Our study reveals the diverse mechanisms by which the mutational landscape of 5’ UTRs can co-opt gene expression and demonstrates that single nucleotide alterations within 5’ UTRs are functional in cancer.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course.

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at −1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the −1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Three-in-One Simultaneous Extraction of Proteins, Metabolites and Lipids for Multi-Omics.

Elucidation of complex molecular networks requires integrative analysis of molecular features and changes at different levels of information flow and regulation. Accordingly, high throughput functional genomics tools such as transcriptomics, proteomics, metabolomics, and lipidomics have emerged to provide system-wide investigations. Unfortunately, analysis of different types of biomolecules requires specific sample extraction procedures in combination with specific analytical instrumentation. The most efficient extraction protocols often only cover a restricted type of biomolecules due to their different physicochemical properties. Therefore, several sets/aliquots of samples are needed for extracting different molecules. Here we adapted a biphasic fractionation method to extract proteins, metabolites, and lipids from the same sample (3-in-1) for liquid chromatography-tandem mass spectrometry (LC-MS/MS) multi-omics. To demonstrate utility of the improved method, we used bacteria-primed Arabidopsis leaves to generate multi-omics datasets from the same sample. In total, we were able to analyze 1849 proteins, 1967 metabolites, and 424 lipid species in single samples. The molecules cover a wide range of biological and molecular processes, and allow quantitative analyses of different molecules and pathways. Our results have shown the clear advantages of the multi-omics method, including sample conservation, high reproducibility, and tight correlation between different types of biomolecules.

Finding function with base editing screens.

Two new studies in Cell report the use of CRISPR–Cas cytosine base editors (CBEs) in mammalian cells for high-throughput functional genomics screens.

GeneWalk identifies relevant gene functions for a biological context using network representation learning

A bottleneck in high-throughput functional genomics experiments is identifying the most important genes and their relevant functions from a list of gene hits. Gene Ontology (GO) enrichment methods provide insight at the gene set level. Here, we introduce GeneWalk (github.com/churchmanlab/genewalk) that identifies individual genes and their relevant functions critical for the experimental setting under examination. After the automatic assembly of an experiment-specific gene regulatory network, GeneWalk uses representation learning to quantify the similarity between vector representations of each gene and its GO annotations, yielding annotation significance scores that reflect the experimental context. By performing gene- and condition-specific functional analysis, GeneWalk converts a list of genes into data-driven hypotheses.

Heterozygous KCNH2 variant phenotyping using Flp-In HEK293 and high-throughput automated patch clamp electrophysiology.

KCNH2 is one of the 59 medically actionable genes recommended by the American College of Medical Genetics for reporting of incidental findings from clinical genomic sequencing. However, half of the reported KCNH2 variants in the ClinVar database are classified as variants of uncertain significance. In the absence of strong clinical phenotypes, there is a need for functional phenotyping to help decipher the significance of variants identified incidentally. Here, we report detailed methods for assessing the molecular phenotype of any KCNH2 missense variant. The key components of the assay include quick and cost-effective generation of a bi-cistronic vector to co-express Wild-type (WT) and any KCNH2 variant allele, generation of stable Flp-In HEK293 cell lines and high-throughput automated patch clamp electrophysiology analysis of channel function. Stable cell lines take 3–4 weeks to produce and can be generated in bulk, which will then allow up to 30 variants to be phenotyped per week after 48 h of channel expression. This high-throughput functional genomics assay will enable a much more rapid assessment of the extent of loss of function of any KCNH2 variant.

Transcriptomic Analysis of Respiratory Tissue and Cell Line Models to Examine Glycosylation Machinery during SARS-CoV-2 Infection.

Glycosylation, being the most abundant post-translational modification, plays a profound role affecting expression, localization and function of proteins and macromolecules in immune response to infection. Presented are the findings of a transcriptomic analysis performed using high-throughput functional genomics data from public repository to examine the altered transcription of the human glycosylation machinery in response to SARS-CoV-2 stimulus and infection. In addition to the conventional in silico functional enrichment analysis methods we also present results from the manual analysis of biomedical literature databases to bring about the biological significance of glycans and glycan-binding proteins in modulating the host immune response during SARS-CoV-2 infection. Our analysis revealed key immunomodulatory lectins, proteoglycans and glycan epitopes implicated in exerting both negative and positive downstream inflammatory signaling pathways, in addition to its vital role as adhesion receptors for SARS-CoV-2 pathogen. A hypothetical correlation of the differentially expressed human glycogenes with the altered host inflammatory response and the cytokine storm-generated in response to SARS-CoV-2 pathogen is proposed. These markers can provide novel insights into the diverse roles and functioning of glycosylation pathways modulated by SARS-CoV-2, provide avenues of stratification, treatment, and targeted approaches for COVID-19 immunity and other viral infectious agents.