Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression.

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Integrated multiplexed assays of variant effect reveal cis-regulatory determinants of catechol-O-methyltransferase gene expression.

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase and decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Single-cell quantification of ribosome occupancy in early mouse development.

Translation regulation is critical for early mammalian embryonic development1. However, previous studies had been restricted to bulk measurements2, precluding precise determination of translation regulation including allele-specific analyses. Here, to address this challenge, we developed a novel microfluidic isotachophoresis (ITP) approach, named RIBOsome profiling via ITP (Ribo-ITP), and characterized translation in single oocytes and embryos during early mouse development. We identified differential translation efficiency as a key mechanism regulating genes involved in centrosome organization and N6-methyladenosine modification of RNAs. Our high-coverage measurements enabled, to our knowledge, the first analysis of allele-specific ribosome engagement in early development. These led to the discovery of stage-specific differential engagement of zygotic RNAs with ribosomes and reduced translation efficiency of transcripts exhibiting allele-biased expression. By integrating our measurements with proteomics data, we discovered that ribosome occupancy in germinal vesicle-stage oocytes is the predominant determinant of protein abundance in the zygote. The Ribo-ITP approach will enable numerous applications by providing high-coverage and high-resolution ribosome occupancy measurements from ultra-low input samples including single cells.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect.

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

Genes with 5' terminal oligopyrimidine tracts preferentially escape global suppression of translation by the SARS-CoV-2 Nsp1 protein.

Viruses rely on the host translation machinery to synthesize their own proteins. Consequently, they have evolved varied mechanisms to co-opt host translation for their survival. SARS-CoV-2 relies on a nonstructural protein, Nsp1, for shutting down host translation. However, it is currently unknown how viral proteins and host factors critical for viral replication can escape a global shutdown of host translation. Here, using a novel FACS-based assay called MeTAFlow, we report a dose-dependent reduction in both nascent protein synthesis and mRNA abundance in cells expressing Nsp1. We perform RNA-seq and matched ribosome profiling experiments to identify gene-specific changes both at the mRNA expression and translation levels. We discover that a functionally coherent subset of human genes is preferentially translated in the context of Nsp1 expression. These genes include the translation machinery components, RNA binding proteins, and others important for viral pathogenicity. Importantly, we uncovered a remarkable enrichment of 5' terminal oligo-pyrimidine (TOP) tracts among preferentially translated genes. Using reporter assays, we validated that 5' UTRs from TOP transcripts can drive preferential expression in the presence of Nsp1. Finally, we found that LARP1, a key effector protein in the mTOR pathway, may contribute to preferential translation of TOP transcripts in response to Nsp1 expression. Collectively, our study suggests fine-tuning of host gene expression and translation by Nsp1 despite its global repressive effect on host protein synthesis.

Genes with 5' terminal oligopyrimidine tracts preferentially escape global suppression of translation by the SARS-CoV-2 Nsp1 protein.

Viruses rely on the host translation machinery to synthesize their own proteins. Consequently, they have evolved varied mechanisms to co-opt host translation for their survival. SARS-CoV-2 relies on a non-structural protein, Nsp1, for shutting down host translation. However, it is currently unknown how viral proteins and host factors critical for viral replication can escape a global shutdown of host translation. Here, using a novel FACS-based assay called MeTAFlow, we report a dose-dependent reduction in both nascent protein synthesis and mRNA abundance in cells expressing Nsp1. We perform RNA-Seq and matched ribosome profiling experiments to identify gene-specific changes both at the mRNA expression and translation level. We discover a functionally-coherent subset of human genes are preferentially translated in the context of Nsp1 expression. These genes include the translation machinery components, RNA binding proteins, and others important for viral pathogenicity. Importantly, we uncovered a remarkable enrichment of 5' terminal oligo-pyrimidine (TOP) tracts among preferentially translated genes. Using reporter assays, we validated that 5' UTRs from TOP transcripts can drive preferential expression in the presence of NSP1. Finally, we found that LARP1, a key effector protein in the mTOR pathway may contribute to preferential translation of TOP transcripts in response to Nsp1 expression. Collectively, our study suggests fine tuning of host gene expression and translation by Nsp1 despite its global repressive effect on host protein synthesis.

Use of Spiked Normalizers to More Precisely Quantify Tumor Markers and Viral Genomes by Massive Parallel Sequencing of Plasma DNA.

A problematic aspect of massive parallel sequencing is that somatic mutations and viral loads are typically quantified as a fraction relative to wild-type human DNA, yet wild-type levels vary with diverse biologic and preanalytic interferences. A novel strategy was devised to quantify target analytes in copies per mL of plasma after normalizing for read counts of spiked DNAs. Five synthetic DNAs (called EndoGenus spikes) were added to plasma before library preparation (modified ArcherDX LiquidPlex 28). By normalizing to the fractional recovery of EndoGenus spike reads, numerical values for each disease marker were reportable in units of copies per mL. To show how well this system operates, replicate assays were performed on 40 mock plasmas having 23 engineered mutations and on 21 natural plasmas. Reads for all five EndoGenus spikes were recovered (means, 313 and 376 copies/mL in mock and natural plasmas, respectively). Normalizing read counts for the proportional recovery of spikes helped control for variables in the multistep protocol, reducing the CV in replicate tests from 34% to 22% for mutations and from 25% to 7% for viral loads. In conclusion, the EndoGenus system is useful for evaluating efficiency of the total test system and for precisely quantifying target molecules. This system may benefit patients being monitored for disease burden while also tracking emerging subclones.

Detecting Gene Rearrangements in Patient Populations Through a 2-Step Diagnostic Test Comprised of Rapid IHC Enrichment Followed by Sensitive Next-Generation Sequencing.

Targeted therapy combined with companion diagnostics has led to the advancement of next-generation sequencing (NGS) for detection of molecular alterations. However, using a diagnostic test to identify patient populations with low prevalence molecular alterations, such as gene rearrangements, poses efficiency, and cost challenges. To address this, we have developed a 2-step diagnostic test to identify NTRK1, NTRK2, NTRK3, ROS1, and ALK rearrangements in formalin-fixed paraffin-embedded clinical specimens. This test is comprised of immunohistochemistry screening using a pan-receptor tyrosine kinase cocktail of antibodies to identify samples expressing TrkA (encoded by NTRK1), TrkB (encoded by NTRK2), TrkC (encoded by NTRK3), ROS1, and ALK followed by an RNA-based anchored multiplex polymerase chain reaction NGS assay. We demonstrate that the NGS assay is accurate and reproducible in identification of gene rearrangements. Furthermore, implementation of an RNA quality control metric to assess the presence of amplifiable nucleic acid input material enables a measure of confidence when an NGS result is negative for gene rearrangements. Finally, we demonstrate that performing a pan-receptor tyrosine kinase immunohistochemistry staining enriches detection of the patient population for gene rearrangements from 4% to 9% and has a 100% negative predictive value. Together, this 2-step assay is an efficient method for detection of gene rearrangements in both clinical testing and studies of archival formalin-fixed paraffin-embedded specimens.

Alphaherpesvirus DNA replication in dissociated human trigeminal ganglia.

Analysis of the frequency and PCR-quantifiable abundance of herpes simplex virus type 1 (HSV-1) and varicella zoster virus (VZV) DNA in multiple biological replicates of cells from dissociated randomly distributed human trigeminal ganglia (TG) of four subjects revealed an increase in both parameters and in both viruses during 5 days of culture, with no further change by 10 days. Dissociated TG provides a platform to analyze initiation of latent virus DNA replication within 5 days of culture.