Sequence variation between 462 human individuals fine-tunes functional sites of RNA processing

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We find frequent splicing order differences between alleles and uncover significant single-nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This includes SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also have a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, show the most allele-specific splicing orders, which frequently co-occur with allele-specific AS, APA, or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Systematically developing a registry of splice-site creating variants utilizing massive publicly available transcriptome sequence data

Genomic variants causing abnormal splicing play important roles in genetic disorders and cancer development. Among them, variants that cause the formation of novel splice-sites (splice-site creating variants, SSCVs) are particularly difficult to identify and often overlooked in genomic studies. Additionally, these SSCVs are frequently considered promising candidates for treatment with splice-switching antisense oligonucleotides (ASOs). To leverage massive transcriptome sequence data such as those available from the Sequence Read Archive, we develop a novel framework to screen for SSCVs solely using transcriptome data. We apply it to 322,072 publicly available transcriptomes and identify 30,130 SSCVs. Among them, 5121 SSCVs affect disease-causing variants. By utilizing this extensive collection of SSCVs, we reveal the characteristics of Alu exonization via SSCVs, especially the hotspots of SSCVs within Alu sequences and their evolutionary relationships. We discover novel gain-of-function SSCVs in the deep intronic region of the NOTCH1 gene and demonstrate that their activation can be suppressed using splice-switching ASOs. Collectively, we provide a systematic approach for automatically acquiring a registry of SSCVs, which facilitates the elucidation of novel biological mechanisms underlying splicing and serves as a valuable resource for drug discovery. The catalogs of SSCVs identified in this study are accessible on the SSCV DB ( https://sscvdb.io ).

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Gene Self-Expressive Networks as a Generalization-Aware Tool to Model Gene Regulatory Networks

Self-expressiveness is a mathematical property that aims at characterizing the relationship between instances in a dataset. This property has been applied widely and successfully in computer-vision tasks, time-series analysis, and to infer underlying network structures in domains including protein signaling interactions and social-networks activity. Nevertheless, despite its potential, self-expressiveness has not been explicitly used to infer gene networks. In this article, we present Generalizable Gene Self-Expressive Networks, a new, interpretable, and generalization-aware formalism to model gene networks, and we propose two methods: GXN•EN and GXN•OMP, based respectively on ElasticNet and OMP (Orthogonal Matching Pursuit), to infer and assess Generalizable Gene Self-Expressive Networks. We evaluate these methods on four Microarray datasets from the DREAM5 benchmark, using both internal and external metrics. The results obtained by both methods are comparable to those obtained by state-of-the-art tools, but are fast to train and exhibit high levels of sparsity, which make them easier to interpret. Moreover we applied these methods to three complex datasets containing RNA-seq informations from different mammalian tissues/cell-types. Lastly, we applied our methodology to compare a normal vs. a disease condition (Alzheimer), which allowed us to detect differential expression of genes’ sub-networks between these two biological conditions. Globally, the gene networks obtained exhibit a sparse and modular structure, with inner communities of genes presenting statistically significant over/under-expression on specific cell types, as well as significant enrichment for some anatomical GO terms, suggesting that such communities may also drive important functional roles.

Use of potassium ion channel and spliceosome proteins as diagnostic biomarkers for sudden unexplained death in schizophrenia

Sudden unexplained death in schizophrenia (SUD-SCZ) is not uncommon and its incidence is approximately three times higher than that in the general population. However, diagnosis of SUD-SCZ remains a great challenge in forensic pathology. This study designed a two-phase study to investigate whether three proteins, namely two potassium ion channel proteins (KCNJ3 and KCNAB1) and one spliceosome protein (SF3B3) that were identified in our previous work, could be applied in the postmortem diagnosis of SUD-SCZ. Immunohistochemical staining of the three biomarkers, followed by a rigorous quantitative analysis, was performed on heart specimens from both SUD-SCZ and control groups. A diagnostic software based on the logistic regression formula derived from the test phase data was then constructed. In the test phase, we found that the staining intensities of KCNJ3, KCNAB1, and SF3B3 were all significantly lower in the SUD-SCZ group (n = 20) as compared with the control group that died from non-natural causes (n = 25), with fold-changes being 14.85 (p < 0.001), 4.13 (p = 0.028) and 2.12 (p = 0.048), respectively. Receiver operating characteristic analysis further illustrated that combination of the three biomarkers achieved the optimal diagnostic specificity (92%) and area under the curve (0.886). In the validation phase, the diagnostic software was confirmed to be a promising tool for predicting the risk of SUD-SCZ in authentic cases. Our study provided a valid strategy towards the practical diagnosis of SUD-SCZ by using KCNJ3, KCNAB1, and SF3B3 proteins as diagnostic biomarkers.

EpisomiR, a New Family of miRNAs, and Its Possible Roles in Human Diseases

MicroRNAs (miRNAs) are synthesized through a canonical pathway and play a role in human diseases, such as cancers and cardiovascular, neurodegenerative, psychiatric, and chronic inflammatory diseases. The development of sequencing technologies has enabled the identification of variations in noncoding miRNAs. These miRNA variants, called isomiRs, are generated through a non-canonical pathway, by several enzymes that alter the length and sequence of miRNAs. The isomiR family is, now, expanding further to include episomiRs, which are miRNAs with different modifications. Since recent findings have shown that isomiRs reflect the cell-specific biological function of miRNAs, knowledge about episomiRs and isomiRs can, possibly, contribute to the optimization of diagnosis and therapeutic technology for precision medicine.

Novel Indel Variation of NPC1 Gene Associates With Risk of Sudden Cardiac Death

Background and Aims: Sudden cardiac death (SCD) was defined as an unexpected death from cardiac causes during a very short duration. It has been reported that Niemann-Pick type C1 (NPC1) gene mutations might be related to cardiovascular diseases. The purpose of the study is to investigate whether common genetic variants of NPC1 is involved in SCD susceptibility.

Methods: Based on a candidate-gene-based approach and systematic screening strategy, this study analyzed an 8-bp insertion/deletion polymorphism (rs150703258) within downstream of NPC1 for the association with SCD risk in Chinese populations using 158 SCD cases and 524 controls. The association of rs150703258 and SCD susceptibility was analyzed using logistic regression. Genotype-phenotype correlation analysis was performed using public database including 1000G, expression quantitative trait loci (eQTL), and further validated by human heart tissues using PCR. Dual-luciferase assay was used to explore the potential regulatory role of rs150703258. Gene expression profiling interactive analysis and transcription factors prediction were performed.

Results: Logistic regression analysis exhibited that the deletion allele of rs150703258 significantly increased the risk of SCD [odds ratio (OR) = 1.329; 95% confidence interval (95%CI):1.03–1.72; p = 0.0289]. Genotype-phenotype correlation analysis showed that the risk allele was significantly associated with higher expression of NPC1 at mRNA and protein expressions level in human heart tissues. eQTL analysis showed NPC1 and C18orf8 (an adjacent gene to NPC1) are both related to rs150703258 and have higher expression level in the samples with deletion allele. Dual-luciferase activity assays indicate a significant regulatory role for rs150703258. Gene expression profiling interactive analysis revealed that NPC1 and C18orf8 seemed to be co-regulated in human blood, arteries and heart tissues. In silico analysis showed that the rs150703258 deletion variant may create transcription factor binding sites. In addition, a rare 12-bp allele (4-bp longer than the insertion allele) of rs150703258 was discovered in the current cohort.

Conclusion: In summary, our study revealed that rs150703258 might contribute to SCD susceptibility by regulating NPC1 and C18orf8 expression. This indel may be a potential marker for risk stratification and molecular diagnosis of SCD. Validations in different ethnic groups with larger sample size and mechanism explorations are warranted to confirm our findings.

A comparison of transcriptome analysis methods with reference genome

Background

The application of RNA-seq technology has become more extensive and the number of analysis procedures available has increased over the past years. Selecting an appropriate workflow has become an important issue for researchers in the field.

Methods

In our study, six popular analytical procedures/pipeline were compared using four RNA-seq datasets from mouse, human, rat, and macaque, respectively. The gene expression value, fold change of gene expression, and statistical significance were evaluated to compare the similarities and differences among the six procedures. qRT-PCR was performed to validate the differentially expressed genes (DEGs) from all six procedures.

Results

Cufflinks-Cuffdiff demands the highest computing resources and Kallisto-Sleuth demands the least. Gene expression values, fold change, p and q values of differential expression (DE) analysis are highly correlated among procedures using HTseq for quantification. For genes with medium expression abundance, the expression values determined using the different procedures were similar. Major differences in expression values come from genes with particularly high or low expression levels. HISAT2-StringTie-Ballgown is more sensitive to genes with low expression levels, while Kallisto-Sleuth may only be useful to evaluate genes with medium to high abundance. When the same thresholds for fold change and p value are chosen in DE analysis, StringTie-Ballgown produce the least number of DEGs, while HTseq-DESeq2, -edgeR or -limma generally produces more DEGs. The performance of Cufflinks-Cuffdiff and Kallisto-Sleuth varies in different datasets. For DEGs with medium expression levels, the biological verification rates were similar among all procedures.

Conclusion

Results are highly correlated among RNA-seq analysis procedures using HTseq for quantification. Difference in gene expression values mainly come from genes with particularly high or low expression levels. Moreover, biological validation rates of DEGs from all six procedures were similar for genes with medium expression levels. Investigators can choose analytical procedures according to their available computer resources, or whether genes of high or low expression levels are of interest. If computer resources are abundant, one can utilize multiple procedures to obtain the intersection of results to get the most reliable DEGs, or to obtain a combination of results to get a more comprehensive DE profile for transcriptomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08465-0.

De novo Splice Site Mutation of the CHD7 Gene in a Chinese Patient with Typical CHARGE Syndrome

INTRODUCTION

CHARGE syndrome (CS, OMIM 214800) is a rare genetic disease characterized by multiple congenital abnormalities, including coloboma, heart defect, atresia of the choanae, retardation of development, genital anomalies, and ear anomalies/deafness. The syndrome is mainly caused by a heterozygous variant in the chromodomain helicase DNA-binding protein 7 (CHD7) gene that encodes the CHD7 protein, involved in the ATP-dependent remodeling of chromatin.

METHODS

In this study, the next-generation sequencing targeted panel was used to detect a de novo variant c.3523-2A>G in the CHD7 gene in a patient with severe CS, congenital heart disease, left coloboma of the choroid, cryptorchidism, and congenital deafness. The Sanger sequencing confirmed the variant and clarified it as de novo variant by short tandem repeat analysis in the patient family. We analyzed the effect of a variant by Minigene assay to evaluate the pathogenicity of the variant.

RESULTS

In summary, cDNA analysis confirmed that c.3523-2A>G variant activates a cryptic splice site, resulting in 172 base pair missing in exon 15, leading to the premature truncation of the CHD7 protein (p.V1175Afs*11).

CONCLUSION

The present study functionally characterized the novel c.3523-2A>G variant in CHD7, providing further confirmatory evidence that it is associated with CS.

Modulation of STIM1 by a risk insertion/deletion polymorphism underlying genetics susceptibility to sudden cardiac death originated from coronary artery disease

Stromal interaction molecule 1 (STIM1), as a dynamic calcium signal transducer and key regulator of cardiomyocyte Ca²⁺ homeostasis, has been implicated in various pathological processes related to sudden cardiac death originated from coronary artery disease (SCD-CAD). In this study, we performed a systematic variant screening on promoter region of STIM1 to filter potential functional genetic variations. Based on the screening results, a 5-bp insertion/deletion (indel) polymorphism (rs3061890) in promoter region of STIM1 was selected as the candidate variant. We investigated the association of rs3061890 with SCD-CAD susceptibility in Chinese Han populations. The homozygote del/del genotype significantly increased risk for SCD-CAD as compared with the ins/ins genotype (odds ratio, 2.86 [95% confidence interval, 1.69-4.29]; P = 2.3 × 10-5). Compared with the common allele, the 5-bp deletion risk allele exhibited lower transcriptional capacity in luciferase assays. Intriguingly, genotype-phenotype correlation studies using human myocardium tissue samples revealed that the expression of STIM1 was associated with the genotype of rs3061890. Computational prediction combined with electrophoretic mobility shift (EMSA) and chromatin immunoprecipitation (ChIP) assays provided convincing evidence for stronger binding affinity of ELF1 (E74 like ETS transcription factor 1) with the deletion allele promoter. Taken together, our findings implied an allele-specific mechanism of regulating the transcription of STIM1 via ELF1, which contribute to SCD-CAD susceptibility. rs3061890 may thus considered as a candidate genetic marker for SCD-CAD prediction and prevention.

A Functional Indel Polymorphism Within MIR155HG Is Associated With Sudden Cardiac Death Risk in a Chinese Population

Sudden cardiac death (SCD) is a devastating complication of multiple disease processes and has gradually became a major public health issue. miR-155 is one of the best characterized miRNAs and plays a critical role in several physiological and pathological process, including cardiovascular diseases. In this study, we systematically screened the whole region of miR-155 host gene (MIR155HG) and identified a 4-bp insertion/deletion variant (rs72014506) residing in the intron region of MIR155HG as the candidate polymorphism. The association of rs72014506 with SCD susceptibility was evaluated using 166 SCD cases and 830 healthy controls in a Chinese population. Logistic regression analysis suggested that the homozygote del/del genotype significantly decreased the risk of SCD [odds ratio (OR) = 0.29; 95% confidence interval (CI) = 0.12–0.74; P_trend = 0.0004]. Further genotype–expression association study using human myocardium tissue samples suggested that the deletion allele was intimately linked to lower the expression of both MIR155HG and mature miR155. Luciferase activity assay also revealed that the deletion allele of rs72014506 inhibited gene transcriptional activity. Finally, we performed electrophoretic mobility shift assay and verified the preferential binding affinity of the deletion allele with POU2F1 (POU domain class 2 transcription factor 1). Collectively, we have successfully identified a SCD risk conferring polymorphism in the MIR155HG gene and a likely biological mechanism for the decreased risk of SCD associated with the deletion allele. This novel variant may thus serve as a potential genetic marker for SCD diagnosis and prevention in natural populations, if validated by further studies with a larger sample size.

Population-scale study of eRNA transcription reveals bipartite functional enhancer architecture

Enhancer RNAs (eRNA) are unstable non-coding RNAs, transcribed bidirectionally from active regulatory sequences, whose expression levels correlate with enhancer activity. We use capped-nascent-RNA sequencing to efficiently capture bidirectional transcription initiation across several human lymphoblastoid cell lines (Yoruba population) and detect ~75,000 eRNA transcription sites with high sensitivity and specificity. The use of nascent-RNA sequencing sidesteps the confounding effect of eRNA instability. We identify quantitative trait loci (QTLs) associated with the level and directionality of eRNA expression. High-resolution analyses of these two types of QTLs reveal distinct positions of enrichment at the central transcription factor (TF) binding regions and at the flanking eRNA initiation regions, both of which are associated with mRNA expression QTLs. These two regions-the central TF-binding footprint and the eRNA initiation cores-define a bipartite architecture of enhancers, inform enhancer function, and can be used as an indicator of the significance of non-coding regulatory variants.

A molecular-level perspective on the frequency, distribution, and consequences of messenger RNA modifications

Cells use chemical modifications to alter the sterics, charge, and conformations of large biomolecules, modulating their biogenesis, function, and stability. Until recently post-transcriptional RNA modifications were thought to be largely limited to nonprotein coding RNA species. However, this dogma has rapidly transformed with the discovery of a host of modifications in protein coding messenger RNAs (mRNAs). Recent advancements in genome-wide sequencing technologies have enabled the identification of mRNA modifications as a potential new frontier in gene regulation-leading to the development of the epitranscriptome field. As a result, there has been a flurry of multiple groundbreaking discoveries, including new modifications, nucleoside modifying enzymes ("writers" and "erasers"), and RNA binding proteins that recognize chemical modifications ("readers"). These discoveries opened the door to understanding how post-transcriptional mRNA modifications can modulate the mRNA lifecycle, and established a link between the epitranscriptome and human health and disease. Despite a rapidly growing recognition of their importance, fundamental questions regarding the identity, prevalence, and functional consequences of mRNA modifications remain to be answered. Here, we highlight quantitative studies that characterize mRNA modification abundance, frequency, and interactions with cellular machinery. As the field progresses, we see a need for the further integration of quantitative and reductionist approaches to complement transcriptome wide studies in order to establish a molecular-level framework for understanding the consequences of mRNA chemical modifications on biological processes. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Processing > RNA Editing and Modification.

A novel splice site indel alteration in the EIF2AK3 gene is responsible for the first cases of Wolcott-Rallison syndrome in Hungary

BACKGROUND

Wolcott-Rallison Syndrome (WRS) is a rare autosomal recessive disease that is the most common cause of neonatal diabetes in consanguineous families. WRS is caused by various genetic alterations of the Eukaryotic Translation Initiation Factor 2-Alpha Kinase 3 (EIF2AK3) gene.

METHODS

Genetic analysis of a consanguineous family where two children were diagnosed with WRS was performed by Sanger sequencing. The altered protein was investigated by in vitro cloning, expression and immunohistochemistry.

RESULTS

The first cases in Hungary, - two patients in one family, where the parents were fourth-degree cousins - showed the typical clinical features of WRS: early onset diabetes mellitus with hyperglycemia, growth retardation, infection-induced multiple organ failure. The genetic background of the disease was a novel alteration in the EIF2AK3 gene involving the splice site of exon 11- intron 11-12 boundary: g.53051_53062delinsTG. According to cDNA sequencing this created a new splice site and resulted in a frameshift and the development of an early termination codon at amino acid position 633 (p.Pro627AspfsTer7). Based on in vitro cloning and expression studies, the truncated protein was functionally inactive. Immunohistochemistry revealed that the intact protein was absent in the islets of pancreas, furthermore insulin expressing cells were also dramatically diminished. Elevated GRP78 and reduced CHOP protein expression were observed in the liver.

CONCLUSIONS

The novel genetic alteration causing the absence of the EIF2AK3 protein resulted in insufficient handling of severe endoplasmic reticulum stress, leading to liver failure and demise of the patients.

The Length of the Expressed 3' UTR Is an Intermediate Molecular Phenotype Linking Genetic Variants to Complex Diseases

In the last decades, genome-wide association studies (GWAS) have uncovered tens of thousands of associations between common genetic variants and complex diseases. However, these statistical associations can rarely be interpreted functionally and mechanistically. As the majority of the disease-associated variants are located far from coding sequences, even the relevant gene is often unclear. A way to gain insight into the relevant mechanisms is to study the genetic determinants of intermediate molecular phenotypes, such as gene expression and transcript structure. We propose a computational strategy to discover genetic variants affecting the relative expression of alternative 3′ untranslated region (UTR) isoforms, generated through alternative polyadenylation, a widespread posttranscriptional regulatory mechanism known to have relevant functional consequences. When applied to a large dataset in which whole genome and RNA sequencing data are available for 373 European individuals, 2,530 genes with alternative polyadenylation quantitative trait loci (apaQTL) were identified. We analyze and discuss possible mechanisms of action of these variants, and we show that they are significantly enriched in GWAS hits, in particular those concerning immune-related and neurological disorders. Our results point to an important role for genetically determined alternative polyadenylation in affecting predisposition to complex diseases, and suggest new ways to extract functional information from GWAS data.

Pathogenicity and selective constraint on variation near splice sites

Mutations that perturb normal pre-mRNA splicing are significant contributors to human disease. We used exome sequencing data from 7833 probands with developmental disorders (DDs) and their unaffected parents, as well as more than 60,000 aggregated exomes from the Exome Aggregation Consortium, to investigate selection around the splice sites and quantify the contribution of splicing mutations to DDs. Patterns of purifying selection, a deficit of variants in highly constrained genes in healthy subjects, and excess de novo mutations in patients highlighted particular positions within and around the consensus splice site of greater functional relevance. By using mutational burden analyses in this large cohort of proband–parent trios, we could estimate in an unbiased manner the relative contributions of mutations at canonical dinucleotides (73%) and flanking noncanonical positions (27%), and calculate the positive predictive value of pathogenicity for different classes of mutations. We identified 18 patients with likely diagnostic de novo mutations in dominant DD-associated genes at noncanonical positions in splice sites. We estimate 35%–40% of pathogenic variants in noncanonical splice site positions are missing from public databases.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

Targeting the Polyadenylation Signal of Pre-mRNA: A New Gene Silencing Approach for Facioscapulohumeral Dystrophy

Facioscapulohumeral dystrophy (FSHD) is characterized by the contraction of the D4Z4 array located in the sub-telomeric region of the chromosome 4, leading to the aberrant expression of the DUX4 transcription factor and the mis-regulation of hundreds of genes. Several therapeutic strategies have been proposed among which the possibility to target the polyadenylation signal to silence the causative gene of the disease. Indeed, defects in mRNA polyadenylation leads to an alteration of the transcription termination, a disruption of mRNA transport from the nucleus to the cytoplasm decreasing the mRNA stability and translation efficiency. This review discusses the polyadenylation mechanisms, why alternative polyadenylation impacts gene expression, and how targeting polyadenylation signal may be a potential therapeutic approach for FSHD.

Structure-mediated modulation of mRNA abundance by A-to-I editing

RNA editing introduces single nucleotide changes to RNA, thus potentially diversifying gene expression. Recent studies have reported significant changes in RNA editing profiles in disease and development. The functional consequences of these widespread alterations remain elusive because of the unknown function of most RNA editing sites. Here, we carry out a comprehensive analysis of A-to-I editomes in human populations. Surprisingly, we observe highly similar editing profiles across populations despite striking differences in the expression levels of ADAR genes. Striving to explain this discrepancy, we uncover a functional mechanism of A-to-I editing in regulating mRNA abundance. We show that A-to-I editing stabilizes RNA secondary structures and reduces the accessibility of AGO2-miRNA to target sites in mRNAs. The editing-dependent stabilization of mRNAs in turn alters the observed editing levels in the stable RNA repertoire. Our study provides valuable insights into the functional impact of RNA editing in human cells.

A to I editing in disease is not fake news

Adenosine deaminases acting on RNA (ADARs) are zinc-containing enzymes that deaminate adenosine bases to inosines within dsRNA regions in transcripts. In short, structured dsRNA hairpins individual adenosine bases may be targeted specifically and edited with up to one hundred percent efficiency, leading to the production of alternative protein variants. However, the majority of editing events occur within longer stretches of dsRNA formed by pairing of repetitive sequences. Here, many different adenosine bases are potential targets but editing efficiency is usually much lower. Recent work shows that ADAR-mediated RNA editing is also required to prevent aberrant activation of antiviral innate immune sensors that detect viral dsRNA in the cytoplasm. Missense mutations in the ADAR1 RNA editing enzyme cause a fatal auto-inflammatory disease, Aicardi-Goutières syndrome (AGS) in affected children. In addition RNA editing by ADARs has been observed to increase in many cancers and also can contribute to vascular disease. Thus the role of RNA editing in the progression of various diseases can no longer be ignored. The ability of ADARs to alter the sequence of RNAs has also been used to artificially target model RNAs in vitro and in cells for RNA editing. Potentially this approach may be used to repair genetic defects and to alter genetic information at the RNA level. In this review we focus on the role of ADARs in disease development and progression and on their potential use to artificially modify RNAs in a targeted manner.

Cleavage and polyadenylation: Ending the message expands gene regulation

Cleavage and polyadenylation (pA) is a fundamental step that is required for the maturation of primary protein encoding transcripts into functional mRNAs that can be exported from the nucleus and translated in the cytoplasm. 3'end processing is dependent on the assembly of a multiprotein processing complex on the pA signals that reside in the pre-mRNAs. Most eukaryotic genes have multiple pA signals, resulting in alternative cleavage and polyadenylation (APA), a widespread phenomenon that is important to establish cell state and cell type specific transcriptomes. Here, we review how pA sites are recognized and comprehensively summarize how APA is regulated and creates mRNA isoform profiles that are characteristic for cell types, tissues, cellular states and disease.