A massively parallel 3' UTR reporter assay reveals relationships between nucleotide content, sequence conservation, and mRNA destabilization

A new perspective on microRNA-guided gene regulation specificity, and its potential generalization to transcription factors and RNA-binding proteins

Our conception of gene regulation specificity has undergone profound changes over the last 20 years. Previously, regulators were considered to control few genes, recognized with exquisite specificity by a 'lock and key' mechanism. However, recently genome-wide exploration of regulator binding site occupancy (whether on DNA or RNA targets) revealed extensive lists of molecular targets for every studied regulator. Such poor biochemical specificity suggested that each regulator controls many genes, collectively contributing to biological phenotypes. Here, I propose a third model, whereby regulators' biological specificity is only partially due to 'lock and key' biochemistry. Rather, regulators affect many genes at the microscopic scale, but biological consequences for most interactions are attenuated at the mesoscopic scale: only a few regulatory events propagate from microscopic to macroscopic scale; others are made inconsequential by homeostatic mechanisms. This model is well supported by the microRNA literature, and data suggest that it extends to other regulators. It reconciles contradicting observations from biochemistry and comparative genomics on one hand and in vivo genetics on the other hand, but this conceptual unification is obscured by common misconceptions and counter-intuitive modes of graphical display. Profound understanding of gene regulation requires conceptual clarification, and better suited statistical analyses and graphical representation.

Identification, Design, and Application of Noncoding Cis-Regulatory Elements

Cis-regulatory elements (CREs) play a pivotal role in orchestrating interactions with trans-regulatory factors such as transcription factors, RNA-binding proteins, and noncoding RNAs. These interactions are fundamental to the molecular architecture underpinning complex and diverse biological functions in living organisms, facilitating a myriad of sophisticated and dynamic processes. The rapid advancement in the identification and characterization of these regulatory elements has been marked by initiatives such as the Encyclopedia of DNA Elements (ENCODE) project, which represents a significant milestone in the field. Concurrently, the development of CRE detection technologies, exemplified by massively parallel reporter assays, has progressed at an impressive pace, providing powerful tools for CRE discovery. The exponential growth of multimodal functional genomic data has necessitated the application of advanced analytical methods. Deep learning algorithms, particularly large language models, have emerged as invaluable tools for deconstructing the intricate nucleotide sequences governing CRE function. These advancements facilitate precise predictions of CRE activity and enable the de novo design of CREs. A deeper understanding of CRE operational dynamics is crucial for harnessing their versatile regulatory properties. Such insights are instrumental in refining gene therapy techniques, enhancing the efficacy of selective breeding programs, pushing the boundaries of genetic innovation, and opening new possibilities in microbial synthetic biology.

Predicting functional UTR variants by integrating region-specific features

The untranslated region (UTR) of messenger ribonucleic acid (mRNA), including the 5'UTR and 3'UTR, plays a critical role in regulating gene expression and translation. Variants within the UTR can lead to changes associated with human traits and diseases; however, computational prediction of UTR variant effect is challenging. Current noncoding variant prediction mainly focuses on the promoters and enhancers, neglecting the unique sequence of the UTR and thereby limiting their predictive accuracy. In this study, using consolidated datasets of UTR variants from disease databases and large-scale experimental data, we systematically analyzed more than 50 region-specific features of UTR, including functional elements, secondary structure, sequence composition and site conservation. Our analysis reveals that certain features, such as C/G-related sequence composition in 5'UTR and A/T-related sequence composition in 3'UTR, effectively differentiate between nonfunctional and functional variant sets, unveiling potential sequence determinants of functional UTR variants. Leveraging these insights, we developed two classification models to predict functional UTR variants using machine learning, achieving an area under the curve (AUC) value of 0.94 for 5'UTR and 0.85 for 3'UTR, outperforming all existing methods. Our models will be valuable for enhancing clinical interpretation of genetic variants, facilitating the prediction and management of disease risk.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

Computational design of mRNA vaccines

mRNA technology has emerged as a successful vaccine platform that offered a swift response to the COVID-19 pandemic. Accumulating evidence shows that vaccine efficacy, thermostability, and other important properties, are largely impacted by intrinsic properties of the mRNA molecule, such as RNA sequence and structure, both of which can be optimized. Designing mRNA sequence for vaccines presents a combinatorial problem due to an extremely large selection space. For instance, due to the degeneracy of the genetic code, there are over 10632 possible mRNA sequences that could encode the spike protein, the COVID-19 vaccines' target. Moreover, designing different elements of the mRNA sequence simultaneously against multiple objectives such as translational efficiency, reduced reactogenicity, and improved stability requires an efficient and sophisticated optimization strategy. Recently, there has been a growing interest in utilizing computational tools to redesign mRNA sequences to improve vaccine characteristics and expedite discovery timelines. In this review, we explore important biophysical features of mRNA to be considered for vaccine design and discuss how computational approaches can be applied to rapidly design mRNA sequences with desirable characteristics.

Critical view on oligo(dT)-based RNA-seq: bias arising, modeling, and mitigating

The precise biological interpretation of oligo(dT)-based RNA sequencing (RNA-seq) datasets, particularly in single-cell RNA-seq (scRNA-seq), is invaluable for understanding complex biological systems. However, the presence of biases can lead to misleading results in downstream analysis. This study has now identified two additional biases that are not accounted for in established bias models: poly(A)-tail length bias and fixed-position GC-content bias. These biases have a significant negative impact on the overall quality of oligo(dT)-based RNA-seq data. To address these biases, we have developed a universal bias-mitigating method based on the lower-affinity binding of short and nonanchored oligo(dT) primers to poly(A) tails. This method significantly reduces poly(A) length bias and completely eliminates fixed-position GC bias. Furthermore, the use of short oligo(dT) with impartial binding behavior toward the diverse poly(A) tails renders RNA-seq with more reliable measurements. The findings of this study are particularly beneficial for scRNA-seq datasets, where accurate benchmarking is critical.

A Cre-dependent massively parallel reporter assay allows for cell-type specific assessment of the functional effects of non-coding elements in vivo

The function of regulatory elements is highly dependent on the cellular context, and thus for understanding the function of elements associated with psychiatric diseases these would ideally be studied in neurons in a living brain. Massively Parallel Reporter Assays (MPRAs) are molecular genetic tools that enable functional screening of hundreds of predefined sequences in a single experiment. These assays have not yet been adapted to query specific cell types in vivo in a complex tissue like the mouse brain. Here, using a test-case 3'UTR MPRA library with genomic elements containing variants from autism patients, we developed a method to achieve reproducible measurements of element effects in vivo in a cell type-specific manner, using excitatory cortical neurons and striatal medium spiny neurons as test cases. This targeted technique should enable robust, functional annotation of genetic elements in the cellular contexts most relevant to psychiatric disease.

CLT-seq as a universal homopolymer-sequencing concept reveals poly(A)-tail-tuned ncRNA regulation

Dynamic tuning of the poly(A) tail is a crucial mechanism for controlling translation and stability of eukaryotic mRNA. Achieving a comprehensive understanding of how this regulation occurs requires unbiased abundance quantification of poly(A)-tail transcripts and simple poly(A)-length measurement using high-throughput sequencing platforms. Current methods have limitations due to complicated setups and elaborate library preparation plans. To address this, we introduce central limit theorem (CLT)-managed RNA-seq (CLT-seq), a simple and straightforward homopolymer-sequencing method. In CLT-seq, an anchor-free oligo(dT) primer rapidly binds to and unbinds from anywhere along the poly(A) tail string, leading to position-directed reverse transcription with equal probability. The CLT mechanism enables the synthesized poly(T) lengths, which correspond to the templated segment of the poly(A) tail, to distribute normally. Based on a well-fitted pseudogaussian-derived poly(A)-poly(T) conversion model, the actual poly(A)-tail profile is reconstructed from the acquired poly(T)-length profile through matrix operations. CLT-seq follows a simple procedure without requiring RNA-related pre-treatment, enrichment or selection, and the CLT-shortened poly(T) stretches are more compatible with existing sequencing platforms. This proof-of-concept approach facilitates direct homopolymer base-calling and features unbiased RNA-seq. Therefore, CLT-seq provides unbiased, robust and cost-efficient transcriptome-wide poly(A)-tail profiling. We demonstrate that CLT-seq on the most common Illumina platform delivers reliable poly(A)-tail profiling at a transcriptome-wide scale in human cellular contexts. We find that the poly(A)-tail-tuned ncRNA regulation undergoes a dynamic, complex process similar to mRNA regulation. Overall, CLT-seq offers a simplified, effective and economical approach to investigate poly(A)-tail regulation, with potential implications for understanding gene expression and identifying therapeutic targets.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Multi-level functional genomics reveals molecular and cellular oncogenicity of patient-based 3' untranslated region mutations

3' untranslated region (3' UTR) somatic mutations represent a largely unexplored avenue of alternative oncogenic gene dysregulation. To determine the significance of 3' UTR mutations in disease, we identify 3' UTR somatic variants across 185 advanced prostate tumors, discovering 14,497 single-nucleotide mutations enriched in oncogenic pathways and 3' UTR regulatory elements. By developing two complementary massively parallel reporter assays, we measure how thousands of patient-based mutations affect mRNA translation and stability and identify hundreds of functional variants that allow us to define determinants of mutation significance. We demonstrate the clinical relevance of these mutations, observing that CRISPR-Cas9 endogenous editing of distinct variants increases cellular stress resistance and that patients harboring oncogenic 3' UTR mutations have a particularly poor prognosis. This work represents an expansive view of the extent to which disease-relevant 3' UTR mutations affect mRNA stability, translation, and cancer progression, uncovering principles of regulatory functionality and potential therapeutic targets in previously unexplored regulatory regions.

Genomic alterations involved in fluoroquinolone resistance development in Staphylococcus aureus

AIM

Fluoroquinolone (FQ) is a potent antibiotic class. However, resistance to this class emerges quickly which hinders its application. In this study, mechanisms leading to the emergence of multidrug-resistant (MDR) Staphylococcus aureus (S. aureus) strains under FQ exposure were investigated.

METHODOLOGY

S. aureus ATCC 29213 was serially exposed to ciprofloxacin (CIP), ofloxacin (OFL), or levofloxacin (LEV) at sub-minimum inhibitory concentrations (sub-MICs) for 12 days to obtain S. aureus -1 strains and antibiotic-free cultured for another 10 days to obtain S. aureus-2 strains. The whole genome (WGS) and target sequencing were applied to analyze genomic alterations; and RT-qPCR was used to access the expressions of efflux-related genes, alternative sigma factors, and genes involved in FQ resistance.

RESULTS

A strong and irreversible increase of MICs was observed in all applied FQs (32 to 128 times) in all S. aureus-1 and remained 16 to 32 times in all S. aureus-2. WGS indicated 10 noticeable mutations occurring in all FQ-exposed S. aureus including 2 insdel mutations in SACOL0573 and rimI; a synonymous mutation in hslO; and 7 missense mutations located in an untranslated region. GrlA, was found mutated (R570H) in all S. aureus-1 and -2. Genes encoding for efflux pumps and their regulator (norA, norB, norC, and mgrA); alternative sigma factors (sigB and sigS); acetyltransferase (rimI); methicillin resistance (fmtB); and hypothetical protein BJI72_0645 were overexpressed in FQ-exposed strains.

CONCLUSION

The emergence of MDR S. aureus was associated with the mutations in the FQ-target sequences and the overexpression of efflux pump systems and their regulators.

Expression of targets of the RNA-binding protein AUF-1 in human airway epithelium indicates its role in cellular senescence and inflammation

Introduction

The RNA-binding protein AU-rich-element factor-1 (AUF-1) participates to posttranscriptional regulation of genes involved in inflammation and cellular senescence, two pathogenic mechanisms of chronic obstructive pulmonary disease (COPD). Decreased AUF-1 expression was described in bronchiolar epithelium of COPD patients versus controls and in vitro cytokine- and cigarette smoke-challenged human airway epithelial cells, prompting the identification of epithelial AUF-1-targeted transcripts and function, and investigation on the mechanism of its loss.

Results

RNA immunoprecipitation-sequencing (RIP-Seq) identified, in the human airway epithelial cell line BEAS-2B, 494 AUF-1-bound mRNAs enriched in their 3'-untranslated regions for a Guanine-Cytosine (GC)-rich binding motif. AUF-1 association with selected transcripts and with a synthetic GC-rich motif were validated by biotin pulldown. AUF-1-targets' steady-state levels were equally affected by partial or near-total AUF-1 loss induced by cytomix (TNFα/IL1β/IFNγ/10 nM each) and siRNA, respectively, with differential transcript decay rates. Cytomix-mediated decrease in AUF-1 levels in BEAS-2B and primary human small-airways epithelium (HSAEC) was replicated by treatment with the senescence- inducer compound etoposide and associated with readouts of cell-cycle arrest, increase in lysosomal damage and senescence-associated secretory phenotype (SASP) factors, and with AUF-1 transfer in extracellular vesicles, detected by transmission electron microscopy and immunoblotting. Extensive in-silico and genome ontology analysis found, consistent with AUF-1 functions, enriched RIP-Seq-derived AUF-1-targets in COPD-related pathways involved in inflammation, senescence, gene regulation and also in the public SASP proteome atlas; AUF-1 target signature was also significantly represented in multiple transcriptomic COPD databases generated from primary HSAEC, from lung tissue and from single-cell RNA-sequencing, displaying a predominant downregulation of expression.

Discussion

Loss of intracellular AUF-1 may alter posttranscriptional regulation of targets particularly relevant for protection of genomic integrity and gene regulation, thus concurring to airway epithelial inflammatory responses related to oxidative stress and accelerated aging. Exosomal-associated AUF-1 may in turn preserve bound RNA targets and sustain their function, participating to spreading of inflammation and senescence to neighbouring cells.

The long and short of EJC-independent nonsense-mediated RNA decay

Nonsense-mediated RNA decay (NMD) plays a dual role as an RNA surveillance mechanism against aberrant transcripts containing premature termination codons and as a gene regulatory mechanism for normal physiological transcripts. This dual function is possible because NMD recognizes its substrates based on the functional definition of a premature translation termination event. An efficient mode of NMD target recognition involves the presence of exon-junction complexes (EJCs) downstream of the terminating ribosome. A less efficient, but highly conserved, mode of NMD is triggered by long 3' untranslated regions (UTRs) that lack EJCs (termed EJC-independent NMD). While EJC-independent NMD plays an important regulatory role across organisms, our understanding of its mechanism, especially in mammalian cells, is incomplete. This review focuses on EJC-independent NMD and discusses the current state of knowledge and factors that contribute to the variability in the efficiency of this mechanism.

Cumulative effects of weakly repressive regulatory regions in the 3' UTR maintain PD-1 expression homeostasis in mammals

PD-1 has become a common target for cancer treatment. However, the molecular regulation of PD-1 expression homeostasis remains unclear. Here we report the PD-1 3' UTR can dramatically repress gene expression via promoting mRNA decay. Deletion of the PD-1 3' UTR inhibits T cell activity and promotes T-ALL cell proliferation. Interestingly, the robust repression is attributable to cumulative effects of many weak regulatory regions, which we show together are better able to maintain PD-1 expression homeostasis. We further identify several RNA binding proteins (RBPs) that modulate PD-1 expression via the 3' UTR, including IGF2BP2, RBM38, SRSF7, and SRSF4. Moreover, despite rapid evolution, PD-1 3' UTRs are functionally conserved and strongly repress gene expression through many common RBP binding sites. These findings reveal a previously unrecognized mechanism of maintaining PD-1 expression homeostasis and might represent a general model for how small regulatory effects play big roles in regulation of gene expression and biology.

Advances in Vertebrate (Cyto)Genomics Shed New Light on Fish Compositional Genome Evolution

Cytogenetic and compositional studies considered fish genomes rather poor in guanine-cytosine content (GC%) because of a putative "sharp increase in genic GC% during the evolution of higher vertebrates". However, the available genomic data have not been exploited to confirm this viewpoint. In contrast, further misunderstandings in GC%, mostly of fish genomes, originated from a misapprehension of the current flood of data. Utilizing public databases, we calculated the GC% in animal genomes of three different, technically well-established fractions: DNA (entire genome), cDNA (complementary DNA), and cds (exons). Our results across chordates help set borders of GC% values that are still incorrect in literature and show: (i) fish in their immense diversity possess comparably GC-rich (or even GC-richer) genomes as higher vertebrates, and fish exons are GC-enriched among vertebrates; (ii) animal genomes generally show a GC-enrichment from the DNA, over cDNA, to the cds level (i.e., not only the higher vertebrates); (iii) fish and invertebrates show a broad(er) inter-quartile range in GC%, while avian and mammalian genomes are more constrained in their GC%. These results indicate no sharp increase in the GC% of genes during the transition to higher vertebrates, as stated and numerously repeated before. We present our results in 2D and 3D space to explore the compositional genome landscape and prepared an online platform to explore the AT/GC compositional genome evolution.

Harnessing introns to monitor gene expression

A number of new technologies have utilized synthetic RNAs which leverage the cell’s RNA splicing machinery to drive expression of gene products. Now a new study reports a technique to dynamically and non-invasively monitor gene expression by embedding reporters within introns contained in the parent gene.

Multiplexed functional genomic assays to decipher the noncoding genome

Linkage disequilibrium and the incomplete regulatory annotation of the noncoding genome complicates the identification of functional noncoding genetic variants and their causal association with disease. Current computational methods for variant prioritization have limited predictive value, necessitating the application of highly parallelized experimental assays to efficiently identify functional noncoding variation. Here, we summarize two distinct approaches, massively parallel reporter assays and CRISPR-based pooled screens and describe their flexible implementation to characterize human noncoding genetic variation at unprecedented scale. Each approach provides unique advantages and limitations, highlighting the importance of multimodal methodological integration. These multiplexed assays of variant effects are undoubtedly poised to play a key role in the experimental characterization of noncoding genetic risk, informing our understanding of the underlying mechanisms of disease-associated loci and the development of more robust predictive classification algorithms.

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

The most stable structure of DNA is the canonical right-handed double helix termed B DNA. However, certain environments and sequence motifs favor alternative conformations, termed non-canonical secondary structures. The roles of DNA and RNA secondary structures in transcriptional regulation remain incompletely understood. However, advances in high-throughput assays have enabled genome wide characterization of some secondary structures. Here, we describe their regulatory functions in promoters and 3’UTRs, providing insights into key mechanisms through which they regulate gene expression. We discuss their implication in human disease, and how advances in molecular technologies and emerging high-throughput experimental methods could provide additional insights.

Context-specific effects of sequence elements on subcellular localization of linear and circular RNAs

Long RNAs vary extensively in their post-transcriptional fates, and this variation is attributed in part to short sequence elements. We used massively parallel RNA assays to study how sequences derived from noncoding RNAs influence the subcellular localization and stability of circular and linear RNAs, including spliced and unspliced forms. We find that the effects of sequence elements strongly depend on the host RNA context, with limited overlap between sequences that drive nuclear enrichment of linear and circular RNAs. Binding of specific RNA binding proteins underpins some of these differences-SRSF1 binding leads to nuclear enrichment of circular RNAs; SAFB binding is associated with nuclear enrichment of predominantly unspliced linear RNAs; and IGF2BP1 promotes export of linear spliced RNA molecules. The post-transcriptional fate of long RNAs is thus dictated by combinatorial contributions of specific sequence elements, of splicing, and of the presence of the terminal features unique to linear RNAs.

Massively parallel analysis of human 3' UTRs reveals that AU-rich element length and registration predict mRNA destabilization

AU-rich elements (AREs) are 3' UTR cis-regulatory elements that regulate the stability of mRNAs. Consensus ARE motifs have been determined, but little is known about how differences in 3' UTR sequences that conform to these motifs affect their function. Here, we use functional annotation of sequences from 3' UTRs (fast-UTR), a massively parallel reporter assay (MPRA), to investigate the effects of 41,288 3' UTR sequence fragments from 4653 transcripts on gene expression and mRNA stability in Jurkat and Beas2B cells. Our analyses demonstrate that the length of an ARE and its registration (the first and last nucleotides of the repeating ARE motif) have significant effects on gene expression and stability. Based on this finding, we propose improved ARE classification and concomitant methods to categorize and predict the effect of AREs on gene expression and stability. Finally, to investigate the advantages of our general experimental design we examine other motifs including constitutive decay elements (CDEs), where we show that the length of the CDE stem-loop has a significant impact on steady-state expression and mRNA stability. We conclude that fast-UTR, in conjunction with our analytical approach, can produce improved yet simple sequence-based rules for predicting the activity of human 3' UTRs.

Sequence determinants as key regulators in gene expression of T cells

T cell homeostasis, T cell differentiation, and T cell effector function rely on the constant fine-tuning of gene expression. To alter the T cell state, substantial remodeling of the proteome is required. This remodeling depends on the intricate interplay of regulatory mechanisms, including post-transcriptional gene regulation. In this review, we discuss how the sequence of a transcript influences these post-transcriptional events. In particular, we review how sequence determinants such as sequence conservation, GC content, and chemical modifications define the levels of the mRNA and the protein in a T cell. We describe the effect of different forms of alternative splicing on mRNA expression and protein production, and their effect on subcellular localization. In addition, we discuss the role of sequences and structures as binding hubs for miRNAs and RNA-binding proteins in T cells. The review thus highlights how the intimate interplay of post-transcriptional mechanisms dictate cellular fate decisions in T cells.

Effects of sequence motifs in the yeast 3' untranslated region determined from massively parallel assays of random sequences

BACKGROUND

The 3' untranslated region (UTR) plays critical roles in determining the level of gene expression through effects on activities such as mRNA stability and translation. Functional elements within this region have largely been identified through analyses of native genes, which contain multiple co-evolved sequence features.

RESULTS

To explore the effects of 3' UTR sequence elements outside of native sequence contexts, we analyze hundreds of thousands of random 50-mers inserted into the 3' UTR of a reporter gene in the yeast Saccharomyces cerevisiae. We determine relative protein expression levels from the fitness of transformants in a growth selection. We find that the consensus 3' UTR efficiency element significantly boosts expression, independent of sequence context; on the other hand, the consensus positioning element has only a small effect on expression. Some sequence motifs that are binding sites for Puf proteins substantially increase expression in the library, despite these proteins generally being associated with post-transcriptional downregulation of native mRNAs. Our measurements also allow a systematic examination of the effects of point mutations within efficiency element motifs across diverse sequence backgrounds. These mutational scans reveal the relative in vivo importance of individual bases in the efficiency element, which likely reflects their roles in binding the Hrp1 protein involved in cleavage and polyadenylation.

CONCLUSIONS

The regulatory effects of some 3' UTR sequence features, like the efficiency element, are consistent regardless of sequence context. In contrast, the consequences of other 3' UTR features appear to be strongly dependent on their evolved context within native genes.

Genome-wide functional screen of 3'UTR variants uncovers causal variants for human disease and evolution

3' untranslated region (3'UTR) variants are strongly associated with human traits and diseases, yet few have been causally identified. We developed the massively parallel reporter assay for 3'UTRs (MPRAu) to sensitively assay 12,173 3'UTR variants. We applied MPRAu to six human cell lines, focusing on genetic variants associated with genome-wide association studies (GWAS) and human evolutionary adaptation. MPRAu expands our understanding of 3'UTR function, suggesting that simple sequences predominately explain 3'UTR regulatory activity. We adapt MPRAu to uncover diverse molecular mechanisms at base pair resolution, including an adenylate-uridylate (AU)-rich element of LEPR linked to potential metabolic evolutionary adaptations in East Asians. We nominate hundreds of 3'UTR causal variants with genetically fine-mapped phenotype associations. Using endogenous allelic replacements, we characterize one variant that disrupts a miRNA site regulating the viral defense gene TRIM14 and one that alters PILRB abundance, nominating a causal variant underlying transcriptional changes in age-related macular degeneration.

The three as: Alternative splicing, alternative polyadenylation and their impact on apoptosis in immune function

The latest advances in next-generation sequencing studies and transcriptomic profiling over the past decade have highlighted a surprising frequency of genes regulated by RNA processing mechanisms in the immune system. In particular, two control steps in mRNA maturation, namely alternative splicing and alternative polyadenylation, are now recognized to occur in the vast majority of human genes. Both have the potential to alter the identity of the encoded protein, as well as control protein abundance or even protein localization or association with other factors. In this review, we will provide a summary of the general mechanisms by which alternative splicing (AS) and alternative polyadenylation (APA) occur, their regulation within cells of the immune system, and their impact on immunobiology. In particular, we will focus on how control of apoptosis by AS and APA is used to tune cell fate during an immune response.

Rare and common variant discovery by whole-genome sequencing of 101 Thoroughbred racehorses

The Thoroughbred breed was formed by crossing Oriental horse breeds and British native horses and is currently used in horseracing worldwide. In this study, we constructed a single-nucleotide variant (SNV) database using data from 101 Thoroughbred racehorses. Whole genome sequencing (WGS) revealed 11,570,312 and 602,756 SNVs in autosomal (1–31) and X chromosomes, respectively, yielding a total of 12,173,068 SNVs. About 6.9% of identified SNVs were rare variants observed only in one allele in 101 horses. The number of SNVs detected in individual horses ranged from 4.8 to 5.3 million. Individual horses had a maximum of 25,554 rare variants; several of these were functional variants, such as non-synonymous substitutions, start-gained, start-lost, stop-gained, and stop-lost variants. Therefore, these rare variants may affect differences in traits and phenotypes among individuals. When observing the distribution of rare variants among horses, one breeding stallion had a smaller number of rare variants compared to other horses, suggesting that the frequency of rare variants in the Japanese Thoroughbred population increases through breeding. In addition, our variant database may provide useful basic information for industrial applications, such as the detection of genetically modified racehorses in gene-doping control and pedigree-registration of racehorses using SNVs as markers.

A prometastatic splicing program regulated by SNRPA1 interactions with structured RNA elements

Aberrant alternative splicing is a hallmark of cancer, yet the underlying regulatory programs that control this process remain largely unknown. Here, we report a systematic effort to decipher the RNA structural code that shapes pathological splicing during breast cancer metastasis. We discovered a previously unknown structural splicing enhancer that is enriched near cassette exons with increased inclusion in highly metastatic cells. We show that the spliceosomal protein small nuclear ribonucleoprotein polypeptide A' (SNRPA1) interacts with these enhancers to promote cassette exon inclusion. This interaction enhances metastatic lung colonization and cancer cell invasion, in part through SNRPA1-mediated regulation of PLEC alternative splicing, which can be counteracted by splicing modulating morpholinos. Our findings establish a noncanonical regulatory role for SNRPA1 as a prometastatic splicing enhancer in breast cancer.

Factors affecting the rapid changes of protein under short-term heat stress

BACKGROUND

Protein content determines the state of cells. The variation in protein abundance is crucial when organisms are in the early stages of heat stress, but the reasons affecting their changes are largely unknown.

RESULTS

We quantified 47,535 mRNAs and 3742 proteins in the filling grains of wheat in two different thermal environments. The impact of mRNA abundance and sequence features involved in protein translation and degradation on protein expression was evaluated by regression analysis. Transcription, codon usage and amino acid frequency were the main drivers of changes in protein expression under heat stress, and their combined contribution explains 58.2 and 66.4% of the protein variation at 30 and 40 °C (20 °C as control), respectively. Transcription contributes more to alterations in protein content at 40 °C (31%) than at 30 °C (6%). Furthermore, the usage of codon AAG may be closely related to the rapid alteration of proteins under heat stress. The contributions of AAG were 24 and 13% at 30 and 40 °C, respectively.

CONCLUSION

In this study, we analyzed the factors affecting the changes in protein expression in the early stage of heat stress and evaluated their influence.

Massively Parallel Reporter Assays: Defining Functional Psychiatric Genetic Variants Across Biological Contexts

Neuropsychiatric phenotypes have long been known to be influenced by heritable risk factors, directly confirmed by the past decade of genetic studies that have revealed specific genetic variants enriched in disease cohorts. However, the initial hope that a small set of genes would be responsible for a given disorder proved false. The more complex reality is that a given disorder may be influenced by myriad small-effect noncoding variants and/or by rare but severe coding variants, many de novo. Noncoding genomic sequences-for which molecular functions cannot usually be inferred-harbor a large portion of these variants, creating a substantial barrier to understanding higher-order molecular and biological systems of disease. Fortunately, novel genetic technologies-scalable oligonucleotide synthesis, RNA sequencing, and CRISPR (clustered regularly interspaced short palindromic repeats)-have opened novel avenues to experimentally identify biologically significant variants en masse. Massively parallel reporter assays (MPRAs) are an especially versatile technique resulting from such innovations. MPRAs are powerful molecular genetics tools that can be used to screen thousands of untranscribed or untranslated sequences and their variants for functional effects in a single experiment. This approach, though underutilized in psychiatric genetics, has several useful features for the field. We review methods for assaying putatively functional genetic variants and regions, emphasizing MPRAs and the opportunities they hold for dissection of psychiatric polygenicity. We discuss literature applying functional assays in neurogenetics, highlighting strengths, caveats, and design considerations-especially regarding disease-relevant variables (cell type, neurodevelopment, and sex), and we ultimately propose applications of MPRA to both computational and experimental neurogenetics of polygenic disease risk.