Functional viromic screens uncover regulatory RNA elements.

The number of sequenced viral genomes has surged recently, presenting an opportunity to understand viral diversity and uncover unknown regulatory mechanisms. Here, we conducted a screening of 30,367 viral segments from 143 species representing 96 genera and 37 families. Using a library of viral segments in 3' UTR, we identified hundreds of elements impacting RNA abundance, translation, and nucleocytoplasmic distribution. To illustrate the power of this approach, we investigated K5, an element conserved in kobuviruses, and found its potent ability to enhance mRNA stability and translation in various contexts, including adeno-associated viral vectors and synthetic mRNAs. Moreover, we identified a previously uncharacterized protein, ZCCHC2, as a critical host factor for K5. ZCCHC2 recruits the terminal nucleotidyl transferase TENT4 to elongate poly(A) tails with mixed sequences, delaying deadenylation. This study provides a unique resource for virus and RNA research and highlights the potential of the virosphere for biological discoveries.

GPC3-Unc5 receptor complex structure and role in cell migration.

Neural migration is a critical step during brain development that requires the interactions of cell-surface guidance receptors. Cancer cells often hijack these mechanisms to disseminate. Here, we reveal crystal structures of Uncoordinated-5 receptor D (Unc5D) in complex with morphogen receptor glypican-3 (GPC3), forming an octameric glycoprotein complex. In the complex, four Unc5D molecules pack into an antiparallel bundle, flanked by four GPC3 molecules. Central glycan-glycan interactions are formed by N-linked glycans emanating from GPC3 (N241 in human) and C-mannosylated tryptophans of the Unc5D thrombospondin-like domains. MD simulations, mass spectrometry and structure-based mutants validate the crystallographic data. Anti-GPC3 nanobodies enhance or weaken Unc5-GPC3 binding and, together with mutant proteins, show that Unc5/GPC3 guide migrating pyramidal neurons in the mouse cortex, and cancer cells in an embryonic xenograft neuroblastoma model. The results demonstrate a conserved structural mechanism of cell guidance, where finely balanced Unc5-GPC3 interactions regulate cell migration.

Neutralizing immunity in vaccine breakthrough infections from the SARS-CoV-2 Omicron and Delta variants.

Virus-like particle (VLP) and live virus assays were used to investigate neutralizing immunity against Delta and Omicron SARS-CoV-2 variants in 259 samples from 128 vaccinated individuals. Following Delta breakthrough infection, titers against WT rose 57-fold and 3.1-fold compared with uninfected boosted and unboosted individuals, respectively, versus only a 5.8-fold increase and 3.1-fold decrease for Omicron breakthrough infection. Among immunocompetent, unboosted patients, Delta breakthrough infections induced 10.8-fold higher titers against WT compared with Omicron (p = 0.037). Decreased antibody responses in Omicron breakthrough infections relative to Delta were potentially related to a higher proportion of asymptomatic or mild breakthrough infections (55.0% versus 28.6%, respectively), which exhibited 12.3-fold lower titers against WT compared with moderate to severe infections (p = 0.020). Following either Delta or Omicron breakthrough infection, limited variant-specific cross-neutralizing immunity was observed. These results suggest that Omicron breakthrough infections are less immunogenic than Delta, thus providing reduced protection against reinfection or infection from future variants.

A single-embryo, single-cell time-resolved model for mouse gastrulation.

Mouse embryonic development is a canonical model system for studying mammalian cell fate acquisition. Recently, single-cell atlases comprehensively charted embryonic transcriptional landscapes, yet inference of the coordinated dynamics of cells over such atlases remains challenging. Here, we introduce a temporal model for mouse gastrulation, consisting of data from 153 individually sampled embryos spanning 36 h of molecular diversification. Using algorithms and precise timing, we infer differentiation flows and lineage specification dynamics over the embryonic transcriptional manifold. Rapid transcriptional bifurcations characterize the commitment of early specialized node and blood cells. However, for most lineages, we observe combinatorial multi-furcation dynamics rather than hierarchical transcriptional transitions. In the mesoderm, dozens of transcription factors combinatorially regulate multifurcations, as we exemplify using time-matched chimeric embryos of Foxc1/Foxc2 mutants. Our study rejects the notion of differentiation being governed by a series of binary choices, providing an alternative quantitative model for cell fate acquisition.

Precise and Programmable Detection of Mutations Using Ultraspecific Riboregulators.

The ability to identify single-nucleotide mutations is critical for probing cell biology and for precise detection of disease. However, the small differences in hybridization energy provided by single-base changes makes identification of these mutations challenging in living cells and complex reaction environments. Here, we report a class of de novo-designed prokaryotic riboregulators that provide ultraspecific RNA detection capabilities in vivo and in cell-free transcription-translation reactions. These single-nucleotide-specific programmable riboregulators (SNIPRs) provide over 100-fold differences in gene expression in response to target RNAs differing by a single nucleotide in E. coli and resolve single epitranscriptomic marks in vitro. By exploiting the programmable SNIPR design, we implement an automated design algorithm to develop riboregulators for a range of mutations associated with cancer, drug resistance, and genetic disorders. Integrating SNIPRs with portable paper-based cell-free reactions enables convenient isothermal detection of cancer-associated mutations from clinical samples and identification of Zika strains through unambiguous colorimetric reactions.

SARS-CoV-2 Reverse Genetics Reveals a Variable Infection Gradient in the Respiratory Tract.

The mode of acquisition and causes for the variable clinical spectrum of coronavirus disease 2019 (COVID-19) remain unknown. We utilized a reverse genetics system to generate a GFP reporter virus to explore severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogenesis and a luciferase reporter virus to demonstrate sera collected from SARS and COVID-19 patients exhibited limited cross-CoV neutralization. High-sensitivity RNA in situ mapping revealed the highest angiotensin-converting enzyme 2 (ACE2) expression in the nose with decreasing expression throughout the lower respiratory tract, paralleled by a striking gradient of SARS-CoV-2 infection in proximal (high) versus distal (low) pulmonary epithelial cultures. COVID-19 autopsied lung studies identified focal disease and, congruent with culture data, SARS-CoV-2-infected ciliated and type 2 pneumocyte cells in airway and alveolar regions, respectively. These findings highlight the nasal susceptibility to SARS-CoV-2 with likely subsequent aspiration-mediated virus seeding to the lung in SARS-CoV-2 pathogenesis. These reagents provide a foundation for investigations into virus-host interactions in protective immunity, host susceptibility, and virus pathogenesis.

Horizontal Cell Biology: Monitoring Global Changes of Protein Interaction States with the Proteome-Wide Cellular Thermal Shift Assay (CETSA).

The cellular thermal shift assay (CETSA) is a biophysical technique allowing direct studies of ligand binding to proteins in cells and tissues. The proteome-wide implementation of CETSA with mass spectrometry detection (MS-CETSA) has now been successfully applied to discover targets for orphan clinical drugs and hits from phenotypic screens, to identify off-targets, and to explain poly-pharmacology and drug toxicity. Highly sensitive multidimensional MS-CETSA implementations can now also access binding of physiological ligands to proteins, such as metabolites, nucleic acids, and other proteins. MS-CETSA can thereby provide comprehensive information on modulations of protein interaction states in cellular processes, including downstream effects of drugs and transitions between different physiological cell states. Such horizontal information on ligandmodulation in cells is largely orthogonal to vertical information on the levels of different proteins and therefore opens novel opportunities to understand operational aspects of cellular proteomes.

Promoter-Intrinsic and Local Chromatin Features Determine Gene Repression in LADs.

It is largely unclear whether genes that are naturally embedded in lamina-associated domains (LADs) are inactive due to their chromatin environment or whether LADs are merely secondary to the lack of transcription. We show that hundreds of human promoters become active when moved from their native LAD position to a neutral context in the same cells, indicating that LADs form a repressive environment. Another set of promoters inside LADs is able to "escape" repression, although their transcription elongation is attenuated. By inserting reporters into thousands of genomic locations, we demonstrate that escaper promoters are intrinsically less sensitive to LAD repression. This is not simply explained by promoter strength but by the interplay between promoter sequence and local chromatin features that vary strongly across LADs. Enhancers also differ in their sensitivity to LAD chromatin. This work provides a general framework for the systematic understanding of gene regulation by repressive chromatin.

Illuminating G-Protein-Coupling Selectivity of GPCRs.

Heterotrimetic G proteins consist of four subfamilies (Gs, Gi/o, Gq/11, and G12/13) that mediate signaling via G-protein-coupled receptors (GPCRs), principally by receptors binding Gα C termini. G-protein-coupling profiles govern GPCR-induced cellular responses, yet receptor sequence selectivity determinants remain elusive. Here, we systematically quantified ligand-induced interactions between 148 GPCRs and all 11 unique Gα subunit C termini. For each receptor, we probed chimeric Gα subunit activation via a transforming growth factor-α (TGF-α) shedding response in HEK293 cells lacking endogenous Gq/11 and G12/13 proteins, and complemented G-protein-coupling profiles through a NanoBiT-G-protein dissociation assay. Interrogation of the dataset identified sequence-based coupling specificity features, inside and outside the transmembrane domain, which we used to develop a coupling predictor that outperforms previous methods. We used the predictor to engineer designer GPCRs selectively coupled to G12. This dataset of fine-tuned signaling mechanisms for diverse GPCRs is a valuable resource for research in GPCR signaling.

A Deep Neural Network for Predicting and Engineering Alternative Polyadenylation.

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3'UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3' end processing, and integrates these features into a comprehensive, interpretable, cis-regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease.

Cell-Penetrating Peptide Mediates Intracellular Membrane Passage of Human Papillomavirus L2 Protein to Trigger Retrograde Trafficking.

Cell-penetrating peptides (CPPs) are short protein segments that can transport cargos into cells. Although CPPs are widely studied as potential drug delivery tools, their role in normal cell physiology is poorly understood. Early during infection, the L2 capsid protein of human papillomaviruses binds retromer, a cytoplasmic trafficking factor required for delivery of the incoming non-enveloped virus into the retrograde transport pathway. Here, we show that the C terminus of HPV L2 proteins contains a conserved cationic CPP that drives passage of a segment of the L2 protein through the endosomal membrane into the cytoplasm, where it binds retromer, thereby sorting the virus into the retrograde pathway for transport to the trans-Golgi network. These experiments define the cell-autonomous biological role of a CPP in its natural context and reveal how a luminal viral protein engages an essential cytoplasmic entry factor.

Modulation of Protein-Interaction States through the Cell Cycle.

Global profiling of protein expression through the cell cycle has revealed subsets of periodically expressed proteins. However, expression levels alone only give a partial view of the biochemical processes determining cellular events. Using a proteome-wide implementation of the cellular thermal shift assay (CETSA) to study specific cell-cycle phases, we uncover changes of interaction states for more than 750 proteins during the cell cycle. Notably, many protein complexes are modulated in specific cell-cycle phases, reflecting their roles in processes such as DNA replication, chromatin remodeling, transcription, translation, and disintegration of the nuclear envelope. Surprisingly, only small differences in the interaction states were seen between the G1 and the G2 phase, suggesting similar hardwiring of biochemical processes in these two phases. The present work reveals novel molecular details of the cell cycle and establishes proteome-wide CETSA as a new strategy to study modulation of protein-interaction states in intact cells.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants.

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Enhancer Function and Evolutionary Roles of Human Accelerated Regions.

Human accelerated regions (HARs) are the fastest-evolving sequences in the human genome. When HARs were discovered in 2006, their function was mysterious due to scant annotation of the noncoding genome. Diverse technologies, from transgenic animals to machine learning, have consistently shown that HARs function as gene regulatory enhancers with significant enrichment in neurodevelopment. It is now possible to quantitatively measure the enhancer activity of thousands of HARs in parallel and model how each nucleotide contributes to gene expression. These strategies have revealed that many human HAR sequences function differently than their chimpanzee orthologs, though individual nucleotide changes in the same HAR may have opposite effects, consistent with compensatory substitutions. To fully evaluate the role of HARs in human evolution, it will be necessary to experimentally and computationally dissect them across more cell types and developmental stages.

Impact of chromatin context on Cas9-induced DNA double-strand break repair pathway balance.

DNA double-strand break (DSB) repair is mediated by multiple pathways. It is thought that the local chromatin context affects the pathway choice, but the underlying principles are poorly understood. Using a multiplexed reporter assay in combination with Cas9 cutting, we systematically measure the relative activities of three DSB repair pathways as a function of chromatin context in >1,000 genomic locations. This reveals that non-homologous end-joining (NHEJ) is broadly biased toward euchromatin, while the contribution of microhomology-mediated end-joining (MMEJ) is higher in specific heterochromatin contexts. In H3K27me3-marked heterochromatin, inhibition of the H3K27 methyltransferase EZH2 reverts the balance toward NHEJ. Single-stranded template repair (SSTR), often used for precise CRISPR editing, competes with MMEJ and is moderately linked to chromatin context. These results provide insight into the impact of chromatin on DSB repair pathway balance and guidance for the design of Cas9-mediated genome editing experiments.

Revealing protein-protein interactions at the transcriptome scale by sequencing.

We describe PROPER-seq (protein-protein interaction sequencing) to map protein-protein interactions (PPIs) en masse. PROPER-seq first converts transcriptomes of input cells into RNA-barcoded protein libraries, in which all interacting protein pairs are captured through nucleotide barcode ligation, recorded as chimeric DNA sequences, and decoded at once by sequencing and mapping. We applied PROPER-seq to human embryonic kidney cells, T lymphocytes, and endothelial cells and identified 210,518 human PPIs (collected in the PROPER v.1.0 database). Among these, 1,365 and 2,480 PPIs are supported by published co-immunoprecipitation (coIP) and affinity purification-mass spectrometry (AP-MS) data, 17,638 PPIs are predicted by the prePPI algorithm without previous experimental validation, and 100 PPIs overlap human synthetic lethal gene pairs. In addition, four previously uncharacterized interaction partners with poly(ADP-ribose) polymerase 1 (PARP1) (a critical protein in DNA repair) known as XPO1, MATR3, IPO5, and LEO1 are validated in vivo. PROPER-seq presents a time-effective technology to map PPIs at the transcriptome scale, and PROPER v.1.0 provides a rich resource for studying PPIs.

Cryo-EM structure of the ß3-adrenergic receptor reveals the molecular basis of subtype selectivity.

The ß3-adrenergic receptor (ß3AR) is predominantly expressed in adipose tissue and urinary bladder and has emerged as an attractive drug target for the treatment of type 2 diabetes, obesity, and overactive bladder (OAB). Here, we report the cryogenic electron microscopy structure of the ß3AR-Gs signaling complex with the selective agonist mirabegron, a first-in-class drug for OAB. Comparison of this structure with the previously reported ß1AR and ß2AR structures reveals a receptor activation mechanism upon mirabegron binding to the orthosteric site. Notably, the narrower exosite in ß3AR creates a perpendicular pocket for mirabegron. Mutational analyses suggest that a combination of both the exosite shape and the amino-acid-residue substitutions defines the drug selectivity of the ßAR agonists. Our findings provide a molecular basis for ßAR subtype selectivity, allowing the design of more-selective agents with fewer adverse effects.

Stress Can Induce Transcription of Toxin-Antitoxin Systems without Activating Toxin.

Toxin-antitoxin (TA) systems are ubiquitous genetic elements in bacterial genomes, but their functions are controversial. Although they are frequently postulated to regulate cell growth following stress, few null phenotypes for TA systems have been reported. Here, we show that TA transcript levels can increase substantially in response to stress, but toxin is not liberated. We find that the growth of an Escherichia coli strain lacking ten TA systems encoding endoribonuclease toxins is not affected following exposure to six stresses that each trigger TA transcription. Additionally, using RNA sequencing, we find no evidence of mRNA cleavage following stress. Stress-induced transcription arises from antitoxin degradation and relief of transcriptional autoregulation. Importantly, although free antitoxin is readily degraded in vivo, antitoxin bound to toxin is protected from proteolysis, preventing release of active toxin. Thus, transcription is not a reliable marker of TA activity, and TA systems do not strongly promote survival following individual stresses.

XACT-Seq Comprehensively Defines the Promoter-Position and Promoter-Sequence Determinants for Initial-Transcription Pausing.

Pausing by RNA polymerase (RNAP) during transcription elongation, in which a translocating RNAP uses a "stepping" mechanism, has been studied extensively, but pausing by RNAP during initial transcription, in which a promoter-anchored RNAP uses a "scrunching" mechanism, has not. We report a method that directly defines the RNAP-active-center position relative to DNA with single-nucleotide resolution (XACT-seq; "crosslink-between-active-center-and-template sequencing"). We apply this method to detect and quantify pausing in initial transcription at 411 (∼4,000,000) promoter sequences in vivo in Escherichia coli. The results show initial-transcription pausing can occur in each nucleotide addition during initial transcription, particularly the first 4 to 5 nucleotide additions. The results further show initial-transcription pausing occurs at sequences that resemble the consensus sequence element for transcription-elongation pausing. Our findings define the positional and sequence determinants for initial-transcription pausing and establish initial-transcription pausing is hard coded by sequence elements similar to those for transcription-elongation pausing.

Transcription and replication: walking a genomic tightrope hand-in-hand.

A recent study by Fenstermaker et al. in Nature describes how transcriptionally active RNA polymerase II (Pol II) clings to the genomic tightrope during the passage of the replication fork and rapidly resumes transcription of immature RNA from both strands of nascent DNA, facilitated by protein-protein interactions between the replication and transcription machineries.