An RNA structure-mediated, posttranscriptional model of human α-1-antitrypsin expression

A new mRNA structure prediction based approach to identifying improved signal peptides for bone morphogenetic protein 2

BACKGROUND

Signal peptide (SP) engineering has proven able to improve production of many proteins yet is a laborious process that still relies on trial and error. mRNA structure around the translational start site is important in translation initiation and has rarely been considered in this context, with recent improvements in in silico mRNA structure potentially rendering it a useful predictive tool for SP selection. Here we attempt to create a method to systematically screen candidate signal peptide sequences in silico based on both their nucleotide and amino acid sequences. Several recently released computational tools were used to predict signal peptide activity (SignalP), localization target (DeepLoc) and predicted mRNA structure (MXFold2). The method was tested with Bone Morphogenetic Protein 2 (BMP2), an osteogenic growth factor used clinically for bone regeneration. It was hoped more effective BMP2 SPs could improve BMP2-based gene therapies and reduce the cost of recombinant BMP2 production.

RESULTS

Amino acid sequence analysis indicated 2,611 SPs from the TGF-β superfamily were predicted to function when attached to BMP2. mRNA structure prediction indicated structures at the translational start site were likely highly variable. The five sequences with the most accessible translational start sites, a codon optimized BMP2 SP variant and the well-established hIL2 SP sequence were taken forward to in vitro testing. The top five candidates showed non-significant improvements in BMP2 secretion in HEK293T cells. All showed reductions in secretion versus the native sequence in C2C12 cells, with several showing large and significant decreases. None of the tested sequences were able to increase alkaline phosphatase activity above background in C2C12s. The codon optimized control sequence and hIL2 SP showed reasonable activity in HEK293T but very poor activity in C2C12.

CONCLUSIONS

These results support the use of peptide sequence based in silico tools for basic predictions around signal peptide activity in a synthetic biology context. However, mRNA structure prediction requires improvement before it can produce reliable predictions for this application. The poor activity of the codon optimized BMP2 SP variant in C2C12 emphasizes the importance of codon choice, mRNA structure, and cellular context for SP activity.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

PCR Analysis to Identify AAT Gene Promoters and Splice Variants

This chapter explores the methods used for the analysis of alpha1-antitrypsin gene expression. This includes the use of the polymerase chain reaction (PCR), reverse transcriptase-PCR (RT-PCR), and whole transcriptome analysis combined with parallel DNA sequencing to understand the processes involved in AAT expression.

Testing Alpha-1 Antitrypsin Deficiency in Black Populations

Alpha-1 antitrypsin (AAT) deficiency (AATD) is an under-recognized hereditary disorder and a significant cause of chronic obstructive pulmonary disease (COPD), a disease that contributes to global mortality. AAT is encoded by the SERPINA1 gene, and severe mutation variants of this gene increase the risk of developing COPD. AATD is more frequently screened for in non-Hispanic White populations. However, AATD is also observed in other ethnic groups and very few studies have documented the mutation frequency in these other ethnic populations. Here, we review the current literature on AATD and allele frequency primarily in Black populations and discuss the possible clinical outcomes of low screening rates in a population that experiences poor health outcomes and whether the low frequency of AATD is related to a lack of screening in this population or a truly low frequency of mutations causing AATD. This review also outlines the harmful SERPINA1 variants, the current epidemiology knowledge of AATD, health inequity in Black populations, AATD prevalence in Black populations, the clinical implications of low screening of AATD in this population, and the possible dangers of not diagnosing or treating AATD.

Unraveling the influences of sequence and position on yeast uORF activity using massively parallel reporter systems and machine learning

Upstream open-reading frames (uORFs) are potent cis-acting regulators of mRNA translation and nonsense-mediated decay (NMD). While both AUG- and non-AUG initiated uORFs are ubiquitous in ribosome profiling studies, few uORFs have been experimentally tested. Consequently, the relative influences of sequence, structural, and positional features on uORF activity have not been determined. We quantified thousands of yeast uORFs using massively parallel reporter assays in wildtype and ∆upf1 yeast. While nearly all AUG uORFs were robust repressors, most non-AUG uORFs had relatively weak impacts on expression. Machine learning regression modeling revealed that both uORF sequences and locations within transcript leaders predict their effect on gene expression. Indeed, alternative transcription start sites highly influenced uORF activity. These results define the scope of natural uORF activity, identify features associated with translational repression and NMD, and suggest that the locations of uORFs in transcript leaders are nearly as predictive as uORF sequences.

Structural analysis of MALAT1 long noncoding RNA in cells and in evolution

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

SERPINA1 long transcripts produce non-secretory alpha1-antitrypsin isoform: In vitro translation in living cells

SERPINA1 is a well-studied serpin gene due to its dramatic impact on human health. Translation initiation at the main SERPINA1 start codon produces the only known alpha1-antitrypsin (AAT) isoform intended for secretion. AAT performs essential functions by inhibiting proteases and modulating immunity. However, SERPINA1 expression at the level of translation is not sufficiently studied. Here we hypothesize that the main SERPINA1 ORF can be alternatively translated, producing a non-secretory AAT isoform by either masking or excluding a signal peptide. We defined SERPINA1 long mRNA isoforms specific for prostate (DU145) and liver (HepG2) cell lines and studied their individual expression by in vitro assay. We found that all long transcripts produce both glycosylated secretory AAT-eGFP fusion protein and non-glycosylated intracellular AAT-eGFP (initiated from an alternative AUG-2 start codon), with the proportion regulated by the SERPINA1 5'-UTR. Both fusion proteins localize to distinct cellular compartments: in contrast to a fusion with the secretory AAT accumulating in the ER, the intracellular one exhibits nuclear-cytoplasmic shuttling. We detected putative endogenous AAT isoform enriching the nuclear speckles. CONCLUSION: Alternative translation initiation might be a mechanism through which SERPINA1 expands the biological diversity of its protein products. Our findings open up new prospects for the study of SERPINA1 gene expression.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

SERPINA1 and More? A Putative Genetic Contributor to Pulmonary Dysfunction in Alpha-1 Antitrypsin Deficiency

Alpha-1 antitrypsin deficiency (AATD) is a common inherited disorder associated with an increased risk of pulmonary disease. Its clinical presentation, including the nature and severity of organ involvement, is highly variable and unpredictable and is not as strongly linked to genotype and environmental exposure (e.g., smoking history) as might be expected. Significant differences were observed within matched populations of severe AATD patients regarding risk of complications, age at onset, and disease course, including the dynamics of lung function decline. Genetic factors are among the putative modifiers contributing to the clinical variability in AATD, yet their role remains elusive. Here, we review and summarise our current understanding of epigenetic and genetic modifiers of pulmonary dysfunction in subjects with AATD.

Biological and Genetic Mechanisms of COPD, Its Diagnosis, Treatment, and Relationship with Lung Cancer

Chronic obstructive pulmonary disease (COPD) is one of the most prevalent chronic adult diseases, with significant worldwide morbidity and mortality. Although long-term tobacco smoking is a critical risk factor for this global health problem, its molecular mechanisms remain unclear. Several phenomena are thought to be involved in the evolution of emphysema, including airway inflammation, proteinase/anti-proteinase imbalance, oxidative stress, and genetic/epigenetic modifications. Furthermore, COPD is one main risk for lung cancer (LC), the deadliest form of human tumor; formation and chronic inflammation accompanying COPD can be a potential driver of malignancy maturation (0.8-1.7% of COPD cases develop cancer/per year). Recently, the development of more research based on COPD and lung cancer molecular analysis has provided new light for understanding their pathogenesis, improving the diagnosis and treatments, and elucidating many connections between these diseases. Our review emphasizes the biological factors involved in COPD and lung cancer, the advances in their molecular mechanisms' research, and the state of the art of diagnosis and treatments. This work combines many biological and genetic elements into a single whole and strongly links COPD with lung tumor features.

Global 5'-UTR RNA structure regulates translation of a SERPINA1 mRNA

SERPINA1 mRNAs encode the protease inhibitor α-1-antitrypsin and are regulated through post-transcriptional mechanisms. α-1-antitrypsin deficiency leads to chronic obstructive pulmonary disease (COPD) and liver cirrhosis, and specific variants in the 5'-untranslated region (5'-UTR) are associated with COPD. The NM_000295.4 transcript is well expressed and translated in lung and blood and features an extended 5'-UTR that does not contain a competing upstream open reading frame (uORF). We show that the 5'-UTR of NM_000295.4 folds into a well-defined multi-helix structural domain. We systematically destabilized mRNA structure across the NM_000295.4 5'-UTR, and measured changes in (SHAPE quantified) RNA structure and cap-dependent translation relative to a native-sequence reporter. Surprisingly, despite destabilizing local RNA structure, most mutations either had no effect on or decreased translation. Most structure-destabilizing mutations retained native, global 5'-UTR structure. However, those mutations that disrupted the helix that anchors the 5'-UTR domain yielded three groups of non-native structures. Two of these non-native structure groups refolded to create a stable helix near the translation initiation site that decreases translation. Thus, in contrast to the conventional model that RNA structure in 5'-UTRs primarily inhibits translation, complex folding of the NM_000295.4 5'-UTR creates a translation-optimized message by promoting accessibility at the translation initiation site.

Quantitative prediction of variant effects on alternative splicing in MAPT using endogenous pre-messenger RNA structure probing

Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor mRNA (pre-mRNA) structure and downstream function are particularly challenging. Here, we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on pre-mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the B^act complex with pre-mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a β-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and 6 newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.

The Potential Role of RNA Structure in Crop Molecular Breeding

The continually growing human population creates a concomitantly increasing demand for nutritious crops with high yields. Advances in high throughput sequencing technologies have revealed the genetic architecture of major crops. This includes extensive information enabling comprehensive genetic markers for breeding selection, new gene discoveries, and novel gene regulatory strategies for crop editing. RNA structure is an important type of genetic feature, essential for post-transcriptional regulation of gene expression. Here, we summarize recent advances in genome-wide RNA structure studies in crops and review the associated RNA structure-mediated regulation of gene expression. We also discuss the functional importance of those single nucleotide variations that induce large RNA structure disparities. Lastly, we discuss the potential role of RNA structure in crop molecular breeding.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.

A Review of Alpha-1 Antitrypsin Binding Partners for Immune Regulation and Potential Therapeutic Application

Alpha-1 antitrypsin (AAT) is the canonical serine protease inhibitor of neutrophil-derived proteases and can modulate innate immune mechanisms through its anti-inflammatory activities mediated by a broad spectrum of protein, cytokine, and cell surface interactions. AAT contains a reactive methionine residue that is critical for its protease-specific binding capacity, whereby AAT entraps the protease on cleavage of its reactive centre loop, neutralises its activity by key changes in its tertiary structure, and permits removal of the AAT-protease complex from the circulation. Recently, however, the immunomodulatory role of AAT has come increasingly to the fore with several prominent studies focused on lipid or protein-protein interactions that are predominantly mediated through electrostatic, glycan, or hydrophobic potential binding sites. The aim of this review was to investigate the spectrum of AAT molecular interactions, with newer studies supporting a potential therapeutic paradigm for AAT augmentation therapy in disorders in which a chronic immune response is strongly linked.

Translation and emerging functions of non-coding RNAs in inflammation and immunity

Regulatory non‐coding RNAs (ncRNAs) including small non‐coding RNAs (sRNAs), long non‐coding RNAs (lncRNAs), and circular RNAs (circRNAs) have gained considerable attention in the last few years. This is mainly due to their condition‐ and tissue‐specific expression and their various modes of action, which suggests them as promising biomarkers and therapeutic targets. One important mechanism of ncRNAs to regulate gene expression is through translation of short open reading frames (sORFs). These sORFs can be located in lncRNAs, in non‐translated regions of mRNAs where upstream ORFs (uORFs) represent the majority, or in circRNAs. Regulation of their translation can function as a quick way to adapt protein production to changing cellular or environmental cues, and can either depend solely on the initiation and elongation of translation, or on the roles of the produced functional peptides. Due to the experimental challenges to pinpoint translation events and to detect the produced peptides, translational regulation through regulatory RNAs is not well studied yet. In the case of circRNAs, they have only recently started to be recognized as regulatory molecules instead of mere artifacts of RNA biosynthesis. Of the many roles described for regulatory ncRNAs, we will focus here on their regulation during inflammation and in immunity.

A novel algorithm for ranking RNA structure candidates

RNA research is advancing at an ever increasing pace. The newest and most state-of-the-art instruments and techniques have made possible the discoveries of new RNAs, and they have carried the field to new frontiers of disease research, vaccine development, therapeutics, and architectonics. Like proteins, RNAs show a marked relationship between structure and function. A deeper grasp of RNAs requires a finer understanding of their elaborate structures. In pursuit of this, cutting-edge experimental and computational structure-probing techniques output several candidate geometries for a given RNA, each of which is perfectly aligned with experimentally determined parameters. Identifying which structure is the most accurate, however, remains a major obstacle. In recent years, several algorithms have been developed for ranking candidate RNA structures in order from most to least probable, though their levels of accuracy and transparency leave room for improvement. Most recently, advances in both areas are demonstrated by rsRNASP, a novel algorithm proposed by Tan et al. rsRNASP is a residue-separation-based statistical potential for three-dimensional structure evaluation, and it outperforms the leading algorithms in the field.

Alternative poly-adenylation modulates α1-antitrypsin expression in chronic obstructive pulmonary disease

α1-anti-trypsin (A1AT), encoded by SERPINA1, is a neutrophil elastase inhibitor that controls the inflammatory response in the lung. Severe A1AT deficiency increases risk for Chronic Obstructive Pulmonary Disease (COPD), however, the role of A1AT in COPD in non-deficient individuals is not well known. We identify a 2.1-fold increase (p = 2.5x10-6) in the use of a distal poly-adenylation site in primary lung tissue RNA-seq in 82 COPD cases when compared to 64 controls and replicate this in an independent study of 376 COPD and 267 controls. This alternative polyadenylation event involves two sites, a proximal and distal site, 61 and 1683 nucleotides downstream of the A1AT stop codon. To characterize this event, we measured the distal ratio in human primary tissue short read RNA-seq data and corroborated our results with long read RNA-seq data. Integrating these results with 3' end RNA-seq and nanoluciferase reporter assay experiments we show that use of the distal site yields mRNA transcripts with over 50-fold decreased translation efficiency and A1AT expression. We identified seven RNA binding proteins using enhanced CrossLinking and ImmunoPrecipitation precipitation (eCLIP) with one or more binding sites in the SERPINA1 3' UTR. We combined these data with measurements of the distal ratio in shRNA knockdown experiments, nuclear and cytoplasmic fractionation, and chemical RNA structure probing. We identify Quaking Homolog (QKI) as a modulator of SERPINA1 mRNA translation and confirm the role of QKI in SERPINA1 translation with luciferase reporter assays. Analysis of single-cell RNA-seq showed differences in the distribution of the SERPINA1 distal ratio among hepatocytes, macrophages, αβ-Tcells and plasma cells in the liver. Alveolar Type 1,2, dendritic cells and macrophages also vary in their distal ratio in the lung. Our work reveals a complex post-transcriptional mechanism that regulates alternative polyadenylation and A1AT expression in COPD.

The Ends Dictate the Means: Promoter Switching in Herpesvirus Gene Expression

Herpesvirus gene expression is dynamic and complex, with distinct complements of viral genes expressed at specific times in different infection contexts. These complex patterns of viral gene expression arise in part from the integration of multiple cellular and viral signals that affect the transcription of viral genes. The use of alternative promoters provides an increased level of control, allowing different promoters to direct the transcription of the same gene in response to distinct temporal and contextual cues. While once considered rare, herpesvirus alternative promoter usage was recently found to be far more pervasive and impactful than previously thought. Here we review several examples of promoter switching in herpesviruses and discuss the functional consequences on the transcriptional and post-transcriptional regulation of viral gene expression.

Alpha-1 Antitrypsin and Hepatocellular Carcinoma in Liver Cirrhosis: SERPINA1 MZ or MS Genotype Carriage Decreases the Risk

Heterozygotes for Z or S alleles of alpha-1-antrypsin (AAT) have low serum AAT levels. Our aim was to compare the risk of hepatocellular carcinoma (HCC) in patients with liver cirrhosis carrying the SERPINA1 MM, MZ and MS genotypes. The study groups consisted of 1119 patients with liver cirrhosis of various aetiologies, and 3240 healthy individuals served as population controls. The MZ genotype was significantly more frequent in the study group (55/1119 vs. 87/3240, p < 0.0001). The MS genotype frequency was comparable in controls (32/119 vs. 101/3240, p = 0.84). MZ and MS heterozygotes had lower serum AAT level than MM homozygotes (medians: 0.90 g/L; 1.40 g/L and 1.67 g/L; p < 0.001 for both). There were significantly fewer patients with HCC in the cirrhosis group among MZ and MS heterozygotes than in MM homozygotes (5/55 and 1/32 respectively, vs. 243/1022, p < 0.01 for both). The risk of HCC was lower in MZ and MS heterozygotes than in MM homozygotes (OR 0.3202; 95% CI 0.1361–0.7719 and OR 0.1522; 95% CI 0.02941–0.7882, respectively). Multivariate analysis of HCC risk factors identified MZ or MS genotype carriage as a protective factor, whereas age, male sex, BMI and viral aetiology of cirrhosis increased HCC risk.

Targeting RNA structures in diseases with small molecules

RNA is crucial for gene expression and regulation. Recent advances in understanding of RNA biochemistry, structure and molecular biology have revealed the importance of RNA structure in cellular processes and diseases. Various approaches to discovering drug-like small molecules that target RNA structure have been developed. This review provides a brief introduction to RNA structural biology and how RNA structures function as disease regulators. We summarize approaches to targeting RNA with small molecules and highlight their advantages, shortcomings and therapeutic potential.

SHAPE Directed Discovery of New Functions in Large RNAs

RNA lies upstream of nearly all biology and functions as the central conduit of information exchange in all cells. RNA molecules encode information both in their primary sequences and in complex structures that form when an RNA folds back on itself. From the time of discovery of mRNA in the late 1950s until quite recently, we had only a rudimentary understanding of RNA structure across vast regions of most messenger and noncoding RNAs. This deficit is now rapidly being addressed, especially by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, mutational profiling (MaP), and closely related platform technologies that, collectively, create chemical microscopes for RNA. These technologies make it possible to interrogate RNA structure, quantitatively, at nucleotide resolution, and at large scales, for entire mRNAs, noncoding RNAs, and viral RNA genomes. By applying comprehensive structure probing to diverse problems, we and others are showing that control of biological function mediated by RNA structure is ubiquitous across prokaryotic and eukaryotic organisms.Work over the past decade using SHAPE-based analyses has clarified key principles. First, the method of RNA structure probing matters. SHAPE-MaP, with its direct and one-step readout that probes nearly every nucleotide by reaction at the 2'-hydroxyl, gives a more detailed and accurate readout than alternatives. Second, comprehensive chemical probing is essential. Focusing on fragments of large RNAs or using meta-gene or statistical analyses to compensate for sparse data sets misses critical features and often yields structure models with poor predictive power. Finally, every RNA has its own internal structural personality. There are myriad ways in which RNA structure modulates sequence accessibility, protein binding, translation, splice-site choice, phase separation, and other fundamental biological processes. In essentially every instance where we have applied rigorous and quantitative SHAPE technologies to study RNA structure-function interrelationships, new insights regarding biological regulatory mechanisms have emerged. RNA elements with more complex higher-order structures appear more likely to contain high-information-content clefts and pockets that bind small molecules, broadly informing a vigorous field of RNA-targeted drug discovery.The broad implications of this collective work are twofold. First, it is long past time to abandon depiction of large RNAs as simple noodle-like or gently flowing molecules. Instead, we need to emphasize that nearly all RNAs are punctuated with distinctive internal structures, a subset of which modulate function in profound ways. Second, structure probing should be an integral component of any effort that seeks to understand the functional nexuses and biological roles of large RNAs.

An analysis of the Kozak consensus in retinal genes and its relevance to gene therapy

PURPOSE

The classic Kozak consensus is a critical genetic element included in gene therapy transgenes to encourage the translation of the therapeutic coding sequence. Despite optimizations of other transgene elements, the Kozak consensus has not yet been considered for potential tissue-specific sequence refinement. We screened the -9 to -1 region relative to the AUG start codon of retina-specific genes to identify whether a Kozak consensus that is different from the classic sequence may be more appropriate for inclusion in gene therapy transgenes that treat inherited retinal disease.

METHODS

Sequences for 135 genes known to cause nonsyndromic inherited retinal disease were extracted from the NCBI database, and the -9 to -1 nucleotides were compared. This panel was then refined to 75 genes with specific retinal functions, for which the -9 to -1 nucleotides were placed in front of a GFP transcript sequence and RNAfold predictions performed. These were compared with a GFP sequence with the classic Kozak consensus (GCCGCCACC), and sequences from retinal genes with minimum free energy (MFE) predictions greater than the reference sequence were selected to generate an optimized Kozak consensus sequence. The original Kozak consensus and the refined retina Kozak consensus were placed upstream of the Renilla luciferase coding sequence, which were used to transfect retinoblastoma cell lines Y-79 and WERI-RB-1 and HEK 293T/17 cells.

RESULTS

The nucleotide frequencies of the original panel of genes were determined to be comparable to the classic Kozak consensus. RNAfold analysis of a GFP transcript with the classic Kozak sequence in the 5' untranslated region (UTR) generated an MFE prediction of -503.3 kcal/mol. RNAfold analysis was then performed with a GFP transcript containing each -9 to -1 Kozak sequence of 75 retinal genes. Thirty-eight of the 75 genes provided a greater MFE value than -503.3 kcal/mol and exhibited an absence of stable secondary structures before the AUG codon. The -9 to -1 nucleotide frequencies of these genes identified a Kozak consensus of ACCGAGACC, differing from the classic Kozak consensus at positions -9, -5, and -4. Applying this sequence to the GFP transcript increased the MFE prediction to -500.1 kcal/mol. The newly identified retina Kozak sequence was also applied to Renilla luciferase plus the REP1 and RPGR transcripts used in current clinical trials. In all examples, the predicted transcript MFE score increased when compared with the current transcript sequences containing classic Kozak consensus sequences. In vitro transfections identified a 7%-9% increase in Renilla activity when incorporating the optimized Kozak sequence.

CONCLUSIONS

The Kozak consensus is a critical element of eukaryotic genes; therefore, it is a required feature of gene therapy transgenes. To date, the classic sequence of GCCRCC (-6 to -1) has typically been incorporated in gene therapy transgenes, but the analysis described here suggests that, for vectors targeting the retina, using a Kozak consensus derived from retinal genes can provide increased expression of the target product.

Therapeutic Potential of Alpha-1 Antitrypsin in Type 1 and Type 2 Diabetes Mellitus

Alpha-1 antitrypsin (AAT) has established anti-inflammatory and immunomodulatory effects in chronic obstructive pulmonary disease but there is increasing evidence of its role in other inflammatory and immune-mediated conditions, like diabetes mellitus (DM). AAT activity is altered in both developing and established type 1 diabetes mellitus (T1DM) as well in established type 2 DM (T2DM). Augmentation therapy with AAT appears to favorably impact T1DM development in mice models and to affect β-cell function and inflammation in humans with T1DM. The role of AAT in T2DM is less clear, but AAT activity appears to be reduced in T2DM. This article reviews these associations and emerging therapeutic strategies using AAT to treat DM.

Visualization of lncRNA and mRNA Structure Models Within the Integrative Genomics Viewer

Every class of RNA forms base-paired structures that impact biological functions. Chemical probing of RNA structure, especially with the advent of strategies such as SHAPE-MaP, vastly expands the scale and quantitative accuracy over which RNA structure can be examined. These methods have enabled large-scale structural studies of mRNAs and lncRNAs, but the length and complexity of these RNAs makes interpretation of the data challenging. We have created modules available through the open-source Integrative Genomics Viewer (IGV) for straightforward visualization of RNA structures along with complementary experimental data. Here we present detailed and stepwise strategies for exploring and visualizing complex RNA structures in IGV. Individuals can use these instructions and supplied sample data to become adept at using IGV to visualize RNA structure models in conjunction with useful allied information.

Plant Immune Mechanisms: From Reductionistic to Holistic Points of View

After three decades of the amazing progress made on molecular studies of plant-microbe interactions (MPMI), we have begun to ask ourselves "what are the major questions still remaining?" as if the puzzle has only a few pieces missing. Such an exercise has ultimately led to the realization that we still have many more questions than answers. Therefore, it would be an impossible task for us to project a coherent "big picture" of the MPMI field in a single review. Instead, we provide our opinions on where we would like to go in our research as an invitation to the community to join us in this exploration of new MPMI frontiers.

LINC00675 Suppresses Cell Proliferation and Migration via Downregulating the H3K4me2 Level at the SPRY4 Promoter in Gastric Cancer

Accumulating evidence indicates that long noncoding RNAs (lncRNAs) are dysregulated in diverse tumors and take a pivotal role in modulating biological processes. In our study, a decreased expression level of LINC00675 in gastric cancer (GC) was first determined by data from The Cancer Genome Atlas (TCGA) and was identified using specimens from GC patients. Then, in vitro and in vivo functional experiments elaborated that LINC00675 could suppress cell proliferation and migration in GC. Multiple differentially expressed genes (DEGs) in LINC00675-overexpressing cells were identified through RNA sequencing analysis. An RNA-binding protein immunoprecipitation (RIP) assay was conducted to reveal that LINC00675 competitively bound with lysine-specific demethylase 1 (LSD1). A coimmunoprecipitation (coIP) assay indicated that LINC00675 overexpression may strengthen the binding of LSD1 and H3K4me2, whereas the chromatin immunoprecipitation (ChIP) assay results verified lower expression of H3K4me2 at the sprouty homolog 4 (SPRY4) promoter region. Together, our research identified that LINC00675 was remarkably downregulated in GC tissues and cells relative to nontumor tissues and cells. LINC00675 could repress GC tumorigenesis and metastasis via competitively binding with LSD1 and intensifying the binding of LSD1 and its target H3K4me2. Importantly, this contributed to attenuated binding of H3K4me2 at the promoter region of oncogene SPRY4 and suppressed SPRY4 transcription, thus suppressing GC cell proliferation and migration.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Control of translation by eukaryotic mRNA transcript leaders-Insights from high-throughput assays and computational modeling

Eukaryotic gene expression is tightly regulated during translation of mRNA to protein. Mis-regulation of translation can lead to aberrant proteins which accumulate in cancers and cause neurodegenerative diseases. Foundational studies on model genes established fundamental roles for mRNA 5' transcript leader (TL) sequences in controlling ribosome recruitment, scanning, and initiation. TL cis-regulatory elements and their corresponding trans-acting factors control cap-dependent initiation under unstressed conditions. Under stress, cap-dependent initiation is suppressed, and specific mRNA structures and sequences promote translation of stress-responsive transcripts to remodel the proteome. In this review, we summarize current knowledge of TL functions in translation initiation. We focus on insights from high-throughput analyses of ribosome occupancy, mRNA structure, RNA Binding Protein occupancy, and massively parallel reporter assays. These data-driven approaches, coupled with computational analyses and modeling, have paved the way for a comprehensive understanding of TL functions. Finally, we will discuss areas of future research on the roles of mRNA sequences and structures in translation. This article is categorized under: Translation > Translation Mechanisms RNA Evolution and Genomics > Computational Analyses of RNA RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Two sides of the same medal: Noncoding mutations reveal new pathological mechanisms and insights into the regulation of gene expression

Noncoding sequences constitute the major part of the human genome and also of pre-mRNAs. Single nucleotide variants in these regions are often overlooked, but may be responsible for much of the variation of phenotypes observed. Mutations in the noncoding part of pre-mRNAs often reveal new and meaningful insights into the regulation of cellular gene expression. Thus, the mechanistic analysis of the pathological mechanism of such mutations will both foster a deeper understanding of the disease and the underlying cellular pathways. Even synonymous mutations can cause diseases, since the primary mRNA sequence not only encodes amino acids, but also encrypts information on RNA-binding proteins and secondary structure. In fact, the RNA sequence directs assembly of a specific mRNP complex, which in turn dictates the fate of the mRNA or regulates its biogenesis. The accumulation of genomic sequence information is increasing at a rapid pace. However, much of the diversity uncovered may not explain the phenotype of a certain syndrome or disease. For this reason, we also emphasize the value of mechanistic studies on pathological mechanisms being complementary to genome-wide studies and bioinformatic approaches. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

SERPINA1 gene identified in RNA-Seq showed strong association with milk protein concentration in Chinese Holstein cows

The detection of candidate genes and mutations associated with phenotypic traits is important for livestock animals. A previous RNA-Seq study revealed that SERPINA1 gene was a functional candidate that may affect milk protein concentration in dairy cows. To further confirm the genetic effect of SERPINA1 on milk protein traits, genetic polymorphisms were identified and genotype-phenotype associations were performed in a large Chinese Holstein cattle population. The entire coding region and the 5′-regulatory region (5′-UTR) of SERPINA1 was sequenced using pooled DNA of 17 unrelated sires. Association studies for five milk production traits were performed using a mixed model with a population encompassing 1,027 Chinese Holstein cows. A total of four SNPs were identified in SERPINA1, among which rs210222822 and rs41257068 presented in exons, rs207601878 presented in an intron, and rs208607693 was in the 5′-UTR. Analyses of pairwise D′ measures of linkage disequilibrium (LD) showed strong linkage among these four SNPs (D′ = 0.99–1.00), and a 9 Kb haplotype block involving three main haplotypes with GTGT, CCCC and CCGT was inferred. An association study revealed that all four single SNPs and their haplotypes had significant genetic effects on milk protein percentage, milk protein yield and milk yield (P = 0.0458 −  < 0.0001). The phenotypic variance ratio for all 11 significant SNP-trait pairs ranged from 1.01% to 7.54%. The candidate gene of SERPINA1 revealed by our previous RNA-Seq study was confirmed to have pronounced effect on milk protein traits on a genome level. Two SNPs (rs208607693 and rs210222822) presented phenotypic variances of approximately 7% and may be used as key or potential markers to assist selection for new lines of cows with high protein concentration.

The Effects of Rare SERPINA1 Variants on Lung Function and Emphysema in SPIROMICS

Rationale: The role of PI (protease inhibitor) type Z heterozygotes and additional rare variant genotypes in the gene encoding alpha-1 antitrypsin, SERPINA1 (serpin peptidase inhibitor, clade A, member 1), in determining chronic obstructive pulmonary disease risk and severity is controversial.Objectives: To comprehensively evaluate the effects of rare SERPINA1 variants on lung function and emphysema phenotypes in subjects with significant tobacco smoke exposure using deep gene resequencing and alpha-1 antitrypsin concentrations.Methods: DNA samples from 1,693 non-Hispanic white individuals, 385 African Americans, and 90 Hispanics with ≥20 pack-years smoking were resequenced for the identification of rare variants (allele frequency < 0.05) in 16.9 kB of SERPINA1.Measurements and Main Results: White PI Z heterozygotes confirmed by sequencing (MZ; n = 74) had lower post-bronchodilator FEV1 (P = 0.007), FEV1/FVC (P = 0.003), and greater computed tomography-based emphysema (P = 0.02) compared with 1,411 white individuals without PI Z, S, or additional rare variants denoted as VR. PI Z-containing compound heterozygotes (ZS/ZVR; n = 7) had lower FEV1/FVC (P = 0.02) and forced expiratory flow, midexpiratory phase (P = 0.009). Nineteen white heterozygotes for five non-S/Z coding variants associated with lower alpha-1 antitrypsin had greater computed tomography-based emphysema compared with those without rare variants. In African Americans, a 5' untranslated region insertion (rs568223361) was associated with lower alpha-1 antitrypsin and functional small airway disease (P = 0.007).Conclusions: In this integrative deep sequencing study of SERPINA1 with alpha-1 antitrypsin concentrations in a heavy smoker and chronic obstructive pulmonary disease cohort, we confirmed the effects of PI Z heterozygote and compound heterozygote genotypes. We demonstrate the cumulative effects of multiple SERPINA1 variants on alpha-1 antitrypsin deficiency, lung function, and emphysema, thus significantly increasing the frequency of SERPINA1 variation associated with respiratory disease in at-risk smokers.

Up-regulated LINC01234 promotes non-small-cell lung cancer cell metastasis by activating VAV3 and repressing BTG2 expression

Background

Long noncoding RNAs (lncRNAs) are known to regulate tumorigenesis and cancer progression, but their contributions to non-small-cell lung cancer (NSCLC) metastasis remain poorly understood. Our previous and other studies have revealed the involvement of upregulated LINC01234 in regulating gastric cancer and colon cancer cells proliferation, and we aimed to investigate whether LINC01234 overexpression also contribute to cancer cells metastasis in this study.

Methods

We collect the NSCLC tissues and adjacent non-tumor tissues and analyzed expression levels of LINC01234 by quantitative reverse-transcription PCR. LINC01234 were knocked down by using siRNAs or shRNAs, and overexpressed by transfection with overexpression vector; RNA levels of miRNA were downregulated or upregulated with inhibitors or mimics. Transwell assays were used to evaluate cell migration and invasive ability; in vivo metastasis experiments were performed to investigate the effect of LINC01234 on NSCLC cells metastasis. Luciferase reporter, RIP, and ChIP assays were used to determine the regulation of LINC01234 on its targets.

Results

LINC01234 expression is increased in NSCLC tissues, and its upregulation is associated with metastasis and shorter survival in NSCLC. Downregulation of LINC01234 impairs cell migration and invasion in vitro, and inhibits cells metastasis in vivo by acting as a competing endogenous RNA for the miR-340-5p and miR-27b-3p. LINC01234 also interacts with the RNA-binding proteins LSD1 and EZH2, leading to histone modification and transcriptional repression of the anti-proliferative genes BTG2.

Conclusions

Taken together, our findings identify two oncogenic regulatory axes in NSCLC centering on LINC01234: one involving miR-340-5p/miR-27b-3p in the cytoplasm and the second involving EZH2, LSD1, and BTG2 in the nucleus. Our study indicates that these genes may be targeted to reduce or prevent NSCLC metastasis.

Impacts of uORF codon identity and position on translation regulation

Translation regulation plays an important role in eukaryotic gene expression. Upstream open reading frames (uORFs) are potent regulatory elements located in 5′ mRNA transcript leaders. Translation of uORFs usually inhibit the translation of downstream main open reading frames, but some enhance expression. While a minority of uORFs encode conserved functional peptides, the coding regions of most uORFs are not conserved. Thus, the importance of uORF coding sequences on their regulatory functions remains largely unknown. We investigated the impact of an uORF coding region on gene regulation by assaying the functions of thousands of variants in the yeast YAP1 uORF. Varying uORF codons resulted in a wide range of functions, including repressing and enhancing expression of the downstream ORF. The presence of rare codons resulted in the most inhibitory YAP1 uORF variants. Inhibitory functions of such uORFs were abrogated by overexpression of complementary tRNA. Finally, regression analysis of our results indicated that both codon identity and position impact uORF function. Our results support a model in which a uORF coding sequence impacts its regulatory functions by altering the speed of uORF translation.

G-quadruplex located in the 5'UTR of the BAG-1 mRNA affects both its cap-dependent and cap-independent translation through global secondary structure maintenance

The anti-apoptotic BAG-1 protein isoforms are known to be overexpressed in colorectal tumors and are considered to be potential therapeutic targets. The isoforms are derived from alternative translation initiations occuring at four in-frame start codons of a single mRNA transcript. Its 5′UTR also contains an internal ribosome entry site (IRES) regulating the cap-independent translation of the transcript. An RNA G-quadruplex (rG4) is located at the 5′end of the BAG-1 5′UTR, upstream of the known cis-regulatory elements. Herein, we observed that the expression of BAG-1 isoforms is post-transcriptionally regulated in colorectal cancer cells and tumors, and that stabilisation of the rG4 by small molecules ligands reduces the expression of endogenous BAG-1 isoforms. We demonstrated a critical role for the rG4 in the control of both cap-dependent and independent translation of the BAG-1 mRNA in colorectal cancer cells. Additionally, we found an upstream ORF that also represses BAG-1 mRNA translation. The structural probing of the complete 5′UTR showed that the rG4 acts as a steric block which controls the initiation of translation at each start codon of the transcript and also maintains the global 5′UTR secondary structure required for IRES-dependent translation.

Guidelines for SHAPE Reagent Choice and Detection Strategy for RNA Structure Probing Studies

Chemical probing is an important tool for characterizing the complex folded structures of RNA molecules, many of which play key cellular roles. Electrophilic SHAPE reagents create adducts at the 2'-hydroxyl position on the RNA backbone of flexible ribonucleotides with relatively little dependence on nucleotide identity. Strategies for adduct detection such as mutational profiling (MaP) allow accurate, automated calculation of relative adduct frequencies for each nucleotide in a given RNA or group of RNAs. A number of alternative reagents and adduct detection strategies have been proposed, especially for use in living cells. Here we evaluate five SHAPE reagents: three previously well-validated reagents 1M7 (1-methyl-7-nitroisatoic anhydride), 1M6 (1-methyl-6-nitroisatoic anhydride), and NMIA ( N-methylisatoic anhydride), one more recently proposed NAI (2-methylnicotinic acid imidazolide), and one novel reagent 5NIA (5-nitroisatoic anhydride). We clarify the importance of carefully designed software in reading out SHAPE experiments using massively parallel sequencing approaches. We examine SHAPE modification in living cells in diverse cell lines, compare MaP and reverse transcription-truncation as SHAPE adduct detection strategies, make recommendations for SHAPE reagent choice, and outline areas for future development.

The Expression of Matryoshka Gene Encoding a Homologue of Kunitz Peptidase Inhibitor Is Regulated Both at the Level of Transcription and Translation

The gene for Kunitz peptidase inhibitor-like protein (KPILP) contains nested alternative open reading frame (aORF) that controls expression of the maternal mRNA. The content of NbKPILP mRNA in intact leaves of Nicotiana benthamiana plant is low but increases significantly upon extended dark exposure or when foreign nucleic acid is overexpressed in the cells. The NbKPILP gene promoter along with the expressed nested aORF are likely to play an important role in maintaining the levels of NbKPILP mRNA. To elucidate the role of NbKPILP promoter, we isolated a fragment of N. benthamiana chromosomal DNA upstream of the NbKPILP transcription start, sequenced it, and created constructs in which reporter E. coli uidA gene coding for β-D-glucuronidase (GUS) was placed under control of the NbKPILP promoter. By assessing the efficacy of uidA mRNA synthesis directed by the NbKPILP promoter and 35S promoter of the cauliflower mosaic virus in a transient expression system, we showed that the levels of GUS accumulation were comparable for both promoters. Prolonged incubation of the agroinjected plants in the darkness stimulated accumulation of the uidA mRNA directed by the NbKPILP promoter. Our experiments indicate that along with regulation at the transcriptional level, expression of NbKPILP mRNA can be affected by expression of the nested aORF controlled by the polypurine block (PPB) located upstream of its start codon, since introduction of mutations in the PPB resulted in significant accumulation of the NbKPILP mRNA. Nucleotide replacement in the aORF start codon led to the drastic increase in the amounts of NbKPILP mRNA and its protein product.

SERPINA1 methylation and lung function in tobacco-smoke exposed European children and adults: a meta-analysis of ALEC population-based cohorts

BACKGROUND

The pathophysiological role of SERPINA1 in respiratory health may be more strongly determined by the regulation of its expression than by common genetic variants. A family based study of predominantly smoking adults found methylation at two Cytosine-phosphate-Guanine sites (CpGs) in SERPINA1 gene to be associated with chronic obstructive pulmonary disease risk. The objective of this study was to confirm the association of lung function with SERPINA1 methylation in general population samples by testing a comprehensive set of CpGs in the SERPINA gene cluster. We considered lung function level and decline in adult smokers from three European population-based cohorts and lung function level and growth in tobacco-smoke exposed children from a birth cohort.

METHODS

DNA methylation using Illumina Infinium Human Methylation 450 k and EPIC beadchips and lung function were measured at two time points in 1076 SAPALDIA, ECRHS and NFBC adult cohort participants and 259 ALSPAC children. Associations of methylation at 119 CpG sites in the SERPINA gene cluster (PP4R4-SERPINA13P) with lung functions and circulating alpha-1-antitripsin (AAT) were assessed using multivariable cross-sectional and longitudinal regression models.

RESULTS

Methylation at cg08257009 in the SERPINA gene cluster, located 32 kb downstream of SERPINA1, not annotated to a gene, was associated with FEV1/FVC at the Bonferroni corrected level in adults, but not in children. None of the methylation signals in the SERPINA1 gene showed associations with lung function after correcting for multiple testing.

CONCLUSIONS

The results do not support a role of SERPINA1 gene methylation as determinant of lung function across the life course in the tobacco smoke exposed general population exposed.

Principles for targeting RNA with drug-like small molecules

Recent studies have indicated the potential to develop small-molecule drugs that act on RNA targets, leading to burgeoning interest in the field. This article discusses general principles for discovering small-molecule drugs that target RNA and argues that the overarching challenge is to identify appropriate target structures in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets.

RNA molecules are essential for cellular information transfer and gene regulation, and RNAs have been implicated in many human diseases. Messenger and non-coding RNAs contain highly structured elements, and evidence suggests that many of these structures are important for function. Targeting these RNAs with small molecules offers opportunities to therapeutically modulate numerous cellular processes, including those linked to 'undruggable' protein targets. Despite this promise, there is currently only a single class of human-designed small molecules that target RNA used clinically — the linezolid antibiotics. However, a growing number of small-molecule RNA ligands are being identified, leading to burgeoning interest in the field. Here, we discuss principles for discovering small-molecule drugs that target RNA and argue that the overarching challenge is to identify appropriate target structures — namely, in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets. If focus is placed on such druggable binding sites in RNA, extensive knowledge of the typical physicochemical properties of drug-like small molecules could then enable small-molecule drug discovery for RNA targets to become (only) roughly as difficult as for protein targets.