Massively parallel functional dissection of schizophrenia-associated noncoding genetic variants.

Schizophrenia (SCZ) is a highly heritable mental disorder with thousands of associated genetic variants located mostly in the noncoding space of the genome. Translating these associations into insights regarding the underlying pathomechanisms has been challenging because the causal variants, their mechanisms of action, and their target genes remain largely unknown. We implemented a massively parallel variant annotation pipeline (MVAP) to perform SCZ variant-to-function mapping at scale in disease-relevant neural cell types. This approach identified 620 functional variants (1.7%) that operate in a highly developmental context and neuronal-activity-dependent manner. Multimodal integration of epigenomic and CRISPRi screening data enabled us to link these functional variants to target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a multistage prioritization strategy to map functional single-nucleotide polymorphism (SNP)-to-gene-to-endophenotype relations and offer biological insights into the context-dependent molecular processes modulated by SCZ-associated genetic variation.

Long Noncoding RNAs: Molecular Modalities to Organismal Functions.

We have known for decades that long noncoding RNAs (lncRNAs) can play essential functions across most forms of life. The maintenance of chromosome length requires an lncRNA (e.g., hTERC) and two lncRNAs in the ribosome that are required for protein synthesis. Thus, lncRNAs can represent powerful RNA machines. More recently, it has become clear that mammalian genomes encode thousands more lncRNAs. Thus, we raise the question: Which, if any, of these lncRNAs could also represent RNA-based machines? Here we synthesize studies that are beginning to address this question by investigating fundamental properties of lncRNA genes, revealing new insights into the RNA structure-function relationship, determining cis- and trans-acting lncRNAs in vivo, and generating new developments in high-throughput screening used to identify functional lncRNAs. Overall, these findings provide a context toward understanding the molecular grammar underlying lncRNA biology.

Persistent immune imprinting occurs after vaccination with the COVID-19 XBB.1.5 mRNA booster in humans.

Immune imprinting describes how the first exposure to a virus shapes immunological outcomes of subsequent exposures to antigenically related strains. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Omicron breakthrough infections and bivalent COVID-19 vaccination primarily recall cross-reactive memory B cells induced by prior Wuhan-Hu-1 spike mRNA vaccination rather than priming Omicron-specific naive B cells. These findings indicate that immune imprinting occurs after repeated Wuhan-Hu-1 spike exposures, but whether it can be overcome remains unclear. To understand the persistence of immune imprinting, we investigated memory and plasma antibody responses after administration of the updated XBB.1.5 COVID-19 mRNA vaccine booster. We showed that the XBB.1.5 booster elicited neutralizing antibody responses against current variants that were dominated by recall of pre-existing memory B cells previously induced by the Wuhan-Hu-1 spike. Therefore, immune imprinting persists after multiple exposures to Omicron spikes through vaccination and infection, including post XBB.1.5 booster vaccination, which will need to be considered to guide future vaccination.

Scalable Functional Assays for the Interpretation of Human Genetic Variation.

Scalable sequence-function studies have enabled the systematic analysis and cataloging of hundreds of thousands of coding and noncoding genetic variants in the human genome. This has improved clinical variant interpretation and provided insights into the molecular, biophysical, and cellular effects of genetic variants at an astonishing scale and resolution across the spectrum of allele frequencies. In this review, we explore current applications and prospects for the field and outline the principles underlying scalable functional assay design, with a focus on the study of single-nucleotide coding and noncoding variants.

High-throughput characterization of functional variants highlights heterogeneity and polygenicity underlying lung cancer susceptibility.

Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility.

Viral T-cell epitopes - Identification, characterization and clinical application.

T-cell immunity, mediated by CD4+ and CD8+ T cells, represents a cornerstone in the control of viral infections. Virus-derived T-cell epitopes are represented by human leukocyte antigen (HLA)-presented viral peptides on the surface of virus-infected cells. They are the prerequisite for the recognition of infected cells by T cells. Knowledge of viral T-cell epitopes provides on the one hand a diagnostic tool to decipher protective T-cell immune responses in the human population and on the other hand various prophylactic and therapeutic options including vaccination approaches and the transfer of virus-specific T cells. Such approaches have already been proven to be effective against various viral infections, particularly in immunocompromised patients lacking sufficient humoral, antibody-based immune response. This review provides an overview on the state of the art as well as current studies regarding the identification and characterization of viral T-cell epitopes and approaches of clinical application. In the first chapter in silico prediction tools and direct, mass spectrometry-based identification of viral T-cell epitopes is compared. The second chapter provides an overview of commonly used assays for further characterization of T-cell responses and phenotypes. The final chapter presents an overview of clinical application of viral T-cell epitopes with a focus on human immunodeficiency virus (HIV), human cytomegalovirus (HCMV) and severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), being representatives of relevant viruses.

The functional impact of 1,570 individual amino acid substitutions in human OTC.

Deleterious mutations in the X-linked gene encoding ornithine transcarbamylase (OTC) cause the most common urea cycle disorder, OTC deficiency. This rare but highly actionable disease can present with severe neonatal onset in males or with later onset in either sex. Individuals with neonatal onset appear normal at birth but rapidly develop hyperammonemia, which can progress to cerebral edema, coma, and death, outcomes ameliorated by rapid diagnosis and treatment. Here, we develop a high-throughput functional assay for human OTC and individually measure the impact of 1,570 variants, 84% of all SNV-accessible missense mutations. Comparison to existing clinical significance calls, demonstrated that our assay distinguishes known benign from pathogenic variants and variants with neonatal onset from late-onset disease presentation. This functional stratification allowed us to identify score ranges corresponding to clinically relevant levels of impairment of OTC activity. Examining the results of our assay in the context of protein structure further allowed us to identify a 13 amino acid domain, the SMG loop, whose function appears to be required in human cells but not in yeast. Finally, inclusion of our data as PS3 evidence under the current ACMG guidelines, in a pilot reclassification of 34 variants with complete loss of activity, would change the classification of 22 from variants of unknown significance to clinically actionable likely pathogenic variants. These results illustrate how large-scale functional assays are especially powerful when applied to rare genetic diseases.

Direct and indirect control of Rho-dependent transcription termination by the Escherichia coli lysC riboswitch.

Bacterial riboswitches are molecular structures that play a crucial role in controlling gene expression to maintain cellular balance. The Escherichia coli lysC riboswitch has been previously shown to regulate gene expression through translation initiation and mRNA decay. Recent research suggests that lysC gene expression is also influenced by Rho-dependent transcription termination. Through a series of in silico, in vitro, and in vivo experiments, we provide experimental evidence that the lysC riboswitch directly and indirectly modulates Rho transcription termination. Our study demonstrates that Rho-dependent transcription termination plays a significant role in the cotranscriptional regulation of lysC expression. Together with previous studies, our work suggests that lysC expression is governed by a lysine-sensing riboswitch that regulates translation initiation, transcription termination, and mRNA degradation. Notably, both Rho and RNase E target the same region of the RNA molecule, implying that RNase E may degrade Rho-terminated transcripts, providing a means to selectively eliminate these incomplete messenger RNAs. Overall, this study sheds light on the complex regulatory mechanisms used by bacterial riboswitches, emphasizing the role of transcription termination in the control of gene expression and mRNA stability.

Mapping Degradation Signals and Pathways in a Eukaryotic N-terminome.

Most eukaryotic proteins are N-terminally acetylated. This modification can be recognized as a signal for selective protein degradation (degron) by the N-end rule pathways. However, the prevalence and specificity of such degrons in the proteome are unclear. Here, by systematically examining how protein turnover is affected by N-terminal sequences, we perform a comprehensive survey of degrons in the yeast N-terminome. We find that approximately 26% of nascent protein N termini encode cryptic degrons. These degrons exhibit high hydrophobicity and are frequently recognized by the E3 ubiquitin ligase Doa10, suggesting a role in protein quality control. In contrast, N-terminal acetylation rarely functions as a degron. Surprisingly, we identify two pathways where N-terminal acetylation has the opposite function and blocks protein degradation through the E3 ubiquitin ligase Ubr1. Our analysis highlights the complexity of N-terminal degrons and argues that hydrophobicity, not N-terminal acetylation, is the predominant feature of N-terminal degrons in nascent proteins.

Advances in computational and experimental approaches for deciphering transcriptional regulatory networks: Understanding the roles of cis-regulatory elements is essential, and recent research utilizing MPRAs, STARR-seq, CRISPR-Cas9, and machine learning has yielded valuable insights.

Understanding the influence of cis-regulatory elements on gene regulation poses numerous challenges given complexities stemming from variations in transcription factor (TF) binding, chromatin accessibility, structural constraints, and cell-type differences. This review discusses the role of gene regulatory networks in enhancing understanding of transcriptional regulation and covers construction methods ranging from expression-based approaches to supervised machine learning. Additionally, key experimental methods, including MPRAs and CRISPR-Cas9-based screening, which have significantly contributed to understanding TF binding preferences and cis-regulatory element functions, are explored. Lastly, the potential of machine learning and artificial intelligence to unravel cis-regulatory logic is analyzed. These computational advances have far-reaching implications for precision medicine, therapeutic target discovery, and the study of genetic variations in health and disease.

Dual effector population modification gene-drive strains of the African malaria mosquitoes, Anopheles gambiae and Anopheles coluzzii.

Proposed genetic approaches for reducing human malaria include population modification, which introduces genes into vector mosquitoes to reduce or prevent parasite transmission. We demonstrate the potential of Cas9/guide RNA (gRNA)-based gene-drive systems linked to dual antiparasite effector genes to spread rapidly through mosquito populations. Two strains have an autonomous gene-drive system coupled to dual anti-Plasmodium falciparum effector genes comprising single-chain variable fragment monoclonal antibodies targeting parasite ookinetes and sporozoites in the African malaria mosquitoes Anopheles gambiae (AgTP13) and Anopheles coluzzii (AcTP13). The gene-drive systems achieved full introduction within 3 to 6 mo after release in small cage trials. Life-table analyses revealed no fitness loads affecting AcTP13 gene-drive dynamics but AgTP13 males were less competitive than wild types. The effector molecules reduced significantly both parasite prevalence and infection intensities. These data supported transmission modeling of conceptual field releases in an island setting that shows meaningful epidemiological impacts at different sporozoite threshold levels (2.5 to 10 k) for human infection by reducing malaria incidence in optimal simulations by 50 to 90% within as few as 1 to 2 mo after a series of releases, and by ≥90% within 3 mo. Modeling outcomes for low sporozoite thresholds are sensitive to gene-drive system fitness loads, gametocytemia infection intensities during parasite challenges, and the formation of potentially drive-resistant genome target sites, extending the predicted times to achieve reduced incidence. TP13-based strains could be effective for malaria control strategies following validation of sporozoite transmission threshold numbers and testing field-derived parasite strains. These or similar strains are viable candidates for future field trials in a malaria-endemic region.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation.

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Compartmentalization of galactan biosynthesis in mycobacteria.

Galactan polymer is a prominent component of the mycobacterial cell wall core. Its biogenesis starts at the cytoplasmic side of the plasma membrane by a build-up of the linker disaccharide [rhamnosyl (Rha) - N-acetyl-glucosaminyl (GlcNAc) phosphate] on the decaprenyl-phosphate carrier. This decaprenyl-P-P-GlcNAc-Rha intermediate is extended by two bifunctional galactosyl transferases, GlfT1 and GlfT2, and then it is translocated to the periplasmic space by an ABC transporter Wzm-Wzt. The cell wall core synthesis is finalized by the action of an array of arabinosyl transferases, mycolyl transferases, and ligases that catalyze an attachment of the arabinogalactan polymer to peptidoglycan through the linker region. Based on visualization of the GlfT2 enzyme fused with fluorescent tags it was proposed that galactan polymerization takes place in a specific compartment of the mycobacterial cell envelope, the intracellular membrane domain, representing pure plasma membrane free of cell wall components (previously denoted as the "PMf" domain), which localizes to the polar region of mycobacteria. In this work, we examined the activity of the galactan-producing cellular machine in the cell-wall containing cell envelope fraction and in the cell wall-free plasma membrane fraction prepared from Mycobacterium smegmatis by the enzyme assays using radioactively labeled substrate UDP-[14C]-galactose as a tracer. We found that despite a high abundance of GlfT2 in both of these fractions as confirmed by their thorough proteomic analyses, galactan is produced only in the reaction mixtures containing the cell wall components. Our findings open the discussion about the distribution of GlfT2 and the regulation of its activity in mycobacteria.

Improving the Accuracy of Bulk Fitness Assays by Correcting Barcode Processing Biases.

Measuring the fitnesses of genetic variants is a fundamental objective in evolutionary biology. A standard approach for measuring microbial fitnesses in bulk involves labeling a library of genetic variants with unique sequence barcodes, competing the labeled strains in batch culture, and using deep sequencing to track changes in the barcode abundances over time. However, idiosyncratic properties of barcodes can induce nonuniform amplification or uneven sequencing coverage that causes some barcodes to be over- or under-represented in samples. This systematic bias can result in erroneous read count trajectories and misestimates of fitness. Here, we develop a computational method, named REBAR (Removing the Effects of Bias through Analysis of Residuals), for inferring the effects of barcode processing bias by leveraging the structure of systematic deviations in the data. We illustrate this approach by applying it to two independent data sets, and demonstrate that this method estimates and corrects for bias more accurately than standard proxies, such as GC-based corrections. REBAR mitigates bias and improves fitness estimates in high-throughput assays without introducing additional complexity to the experimental protocols, with potential applications in a range of experimental evolution and mutation screening contexts.

Intricacies of single-cell multi-omics data integration.

A wealth of single-cell protocols makes it possible to characterize different molecular layers at unprecedented resolution. Integrating the resulting multimodal single-cell data to find cell-to-cell correspondences remains a challenge. We argue that data integration needs to happen at a meaningful biological level of abstraction and that it is necessary to consider the inherent discrepancies between modalities to strike a balance between biological discovery and noise removal. A survey of current methods reveals that a distinction between technical and biological origins of presumed unwanted variation between datasets is not yet commonly considered. The increasing availability of paired multimodal data will aid the development of improved methods by providing a ground truth on cell-to-cell matches.

High-resolution mapping of function and protein binding in an RNA nuclear enrichment sequence.

The functions of long RNAs, including mRNAs and long noncoding RNAs (lncRNAs), critically depend on their subcellular localization. The identity of the sequences that dictate subcellular localization and their high-resolution anatomy remain largely unknown. We used a suite of massively parallel RNA assays and libraries containing thousands of sequence variants to pinpoint the functional features within the SIRLOIN element, which dictates nuclear enrichment through hnRNPK recruitment. In addition, we profiled the endogenous SIRLOIN RNA-nucleoprotein complex and identified the nuclear RNA-binding proteins SLTM and SNRNP70 as novel SIRLOIN binders. Taken together, using massively parallel assays, we identified the features that dictate binding of hnRNPK, SLTM, and SNRNP70 to SIRLOIN and found that these factors are jointly required for SIRLOIN activity. Our study thus provides a roadmap for high-throughput dissection of functional sequence elements in long RNAs.

The functional impact of BRCA1 BRCT domain variants using multiplexed DNA double-strand break repair assays.

Pathogenic variants in BRCA1 are associated with a greatly increased risk of hereditary breast and ovarian cancer (HBOC). With the increased availability and affordability of genetic testing, many individuals have been identified with BRCA1 variants of uncertain significance (VUSs), which are individually detected in the population too infrequently to ascertain a clinical risk. Functional assays can be used to experimentally assess the effects of these variants. In this study, we used multiplexed DNA repair assays of variants in the BRCA1 carboxyl terminus to functionally characterize 2,271 variants for homology-directed repair function (HDR) and 1,427 variants for cisplatin resistance (CR). We found a high level of consistent results (Pearson's r = 0.74) in the two multiplexed functional assays with non-functional variants located within regions of the BRCA1 protein necessary for its tumor suppression activity. In addition, functional categorizations of variants tested in the multiplex HDR and CR assays correlated with known clinical significance and with other functional assays for BRCA1 (Pearson's r = 0.53 to 0.71). The results of the multiplex HDR and CR assays are useful resources for characterizing large numbers of BRCA1 VUSs.

Microtubules under mechanical pressure can breach dense actin networks.

The crosstalk between the actin network and microtubules is essential for cell polarity. It orchestrates microtubule organization within the cell, driven by the asymmetry of actin architecture along the cell periphery. The physical intertwining of these networks regulates spatial organization and force distribution in the microtubule network. Although their biochemical interactions are becoming clearer, the mechanical aspects remain less understood. To explore this mechanical interplay, we developed an in vitro reconstitution assay to investigate how dynamic microtubules interact with various actin filament structures. Our findings revealed that microtubules can align and move along linear actin filament bundles through polymerization force. However, they are unable to pass through when encountering dense branched actin meshworks, similar to those present in the lamellipodium along the periphery of the cell. Interestingly, immobilizing microtubules through crosslinking with actin or other means allow the buildup of pressure, enabling them to breach these dense actin barriers. This mechanism offers insights into microtubule progression towards the cell periphery, with them overcoming obstacles within the denser parts of the actin network and ultimately contributing to cell polarity establishment.

Kinesin-1 transports morphologically distinct intracellular virions during vaccinia infection.

Intracellular mature viruses (IMVs) are the first and most abundant infectious form of vaccinia virus to assemble during its replication cycle. IMVs can undergo microtubule-based motility, but their directionality and the motor involved in their transport remain unknown. Here, we demonstrate that IMVs, like intracellular enveloped viruses (IEVs), the second form of vaccinia that are wrapped in Golgi-derived membranes, recruit kinesin-1 and undergo anterograde transport. In vitro reconstitution of virion transport in infected cell extracts revealed that IMVs and IEVs move toward microtubule plus ends with respective velocities of 0.66 and 0.56 µm/s. Quantitative imaging established that IMVs and IEVs recruit an average of 139 and 320 kinesin-1 motor complexes, respectively. In the absence of kinesin-1, there was a near-complete loss of in vitro motility and reduction in the intracellular spread of both types of virions. Our observations demonstrate that kinesin-1 transports two morphologically distinct forms of vaccinia. Reconstitution of vaccinia-based microtubule motility in vitro provides a new model to elucidate how motor number and regulation impacts transport of a bona fide kinesin-1 cargo.

Hammerhead ribozyme-based U-insertion and deletion RNA editing assays for multiplexing in HTS applications.

Untranslatable mitochondrial transcripts in kinetoplastids are decrypted post-transcriptionally through an RNA editing process that entails uridine insertion/deletion. This unique stepwise process is mediated by the editosome, a multiprotein complex that is a validated drug target of considerable interest in addressing the unmet medical needs for kinetoplastid diseases. With that objective, several in vitro RNA editing assays have been developed, albeit with limited success in discovering potent inhibitors. This manuscript describes the development of three hammerhead ribozyme (HHR) FRET reporter-based RNA editing assays for precleaved deletion, insertion, and ligation assays that bypass the rate-limiting endonucleolytic cleavage step, providing information on U-deletion, U-insertion, and ligation activities. These assays exhibit higher editing efficiencies in shorter incubation times while requiring significantly less purified editosome and 10,000-fold less ATP than the previously published full round of in vitro RNA editing assay. Moreover, modifications in the reporter ribozyme sequence enable the feasibility of multiplexing a ribozyme-based insertion/deletion editing (RIDE) assay that simultaneously surveils U-insertion and deletion editing suitable for HTS. These assays can be used to find novel chemical compounds with chemotherapeutic applications or as probes for studying the editosome machinery.