Genetic variation and microRNA targeting of A-to-I RNA editing fine tune human tissue transcriptomes.

BACKGROUND: A-to-I RNA editing diversifies the transcriptome and has multiple downstream functional effects. Genetic variation contributes to RNA editing variability between individuals and has the potential to impact phenotypic variability. RESULTS: We analyze matched genetic and transcriptomic data in 49 tissues across 437 individuals to identify RNA editing events that are associated with genetic variation. Using an RNA editing quantitative trait loci (edQTL) mapping approach, we identify 3117 unique RNA editing events associated with a cis genetic polymorphism. Fourteen percent of these edQTL events are also associated with genetic variation in their gene expression. A subset of these events are associated with genome-wide association study signals of complex traits or diseases. We determine that tissue-specific levels of ADAR and ADARB1 are able to explain a subset of tissue-specific edQTL events. We find that certain microRNAs are able to differentiate between the edited and unedited isoforms of their targets. Furthermore, microRNAs can generate an expression quantitative trait loci (eQTL) signal from an edQTL locus by microRNA-mediated transcript degradation in an editing-specific manner. By integrative analyses of edQTL, eQTL, and microRNA expression profiles, we computationally discover and experimentally validate edQTL-microRNA pairs for which the microRNA may generate an eQTL signal from an edQTL locus in a tissue-specific manner. CONCLUSIONS: Our work suggests a mechanism in which RNA editing variability can influence the phenotypes of complex traits and diseases by altering the stability and steady-state level of critical RNA molecules.

The E3 ubiquitin ligase Cul4b promotes CD4+ T cell expansion by aiding the repair of damaged DNA.

The capacity for T cells to become activated and clonally expand during pathogen invasion is pivotal for protective immunity. Our understanding of how T cell receptor (TCR) signaling prepares cells for this rapid expansion remains limited. Here we provide evidence that the E3 ubiquitin ligase Cullin-4b (Cul4b) regulates this process. The abundance of total and neddylated Cul4b increased following TCR stimulation. Disruption of Cul4b resulted in impaired proliferation and survival of activated T cells. Additionally, Cul4b-deficient CD4+ T cells accumulated DNA damage. In T cells, Cul4b preferentially associated with the substrate receptor DCAF1, and Cul4b and DCAF1 were found to interact with proteins that promote the sensing or repair of damaged DNA. While Cul4b-deficient CD4+ T cells showed evidence of DNA damage sensing, downstream phosphorylation of SMC1A did not occur. These findings reveal an essential role for Cul4b in promoting the repair of damaged DNA to allow survival and expansion of activated T cells.

Detecting Allele-Specific Alternative Splicing from Population-Scale RNA-Seq Data.

RNA sequencing (RNA-seq) is a powerful technology for studying human transcriptome variation. We introduce PAIRADISE (Paired Replicate Analysis of Allelic Differential Splicing Events), a method for detecting allele-specific alternative splicing (ASAS) from RNA-seq data. Unlike conventional approaches that detect ASAS events one sample at a time, PAIRADISE aggregates ASAS signals across multiple individuals in a population. By treating the two alleles of an individual as paired, and multiple individuals sharing a heterozygous SNP as replicates, we formulate ASAS detection using PAIRADISE as a statistical problem for identifying differential alternative splicing from RNA-seq data with paired replicates. PAIRADISE outperforms alternative statistical models in simulation studies. Applying PAIRADISE to replicate RNA-seq data of a single individual and to population-scale RNA-seq data across many individuals, we detect ASAS events associated with genome-wide association study (GWAS) signals of complex traits or diseases. Additionally, PAIRADISE ASAS analysis detects the effects of rare variants on alternative splicing. PAIRADISE provides a useful computational tool for elucidating the genetic variation and phenotypic association of alternative splicing in populations.

miR-128 Restriction of LINE-1 (L1) Retrotransposition Is Dependent on Targeting hnRNPA1 mRNA.

The majority of the human genome is made of transposable elements, giving rise to interspaced repeats, including Long INterspersed Element-1s (LINE-1s or L1s). L1s are active human transposable elements involved in genomic diversity and evolution; however, they can also contribute to genomic instability and diseases. L1s require host factors to complete their life cycles, whereas the host has evolved numerous mechanisms to restrict L1-induced mutagenesis. Restriction mechanisms in somatic cells include methylation of the L1 promoter, anti-viral factors and RNA-mediated processes such as small RNAs. microRNAs (miRNAs or miRs) are small non-coding RNAs that post-transcriptionally repress multiple target genes often found in the same cellular pathways. We have recently established that miR-128 functions as a novel restriction factor inhibiting L1 mobilization in somatic cells. We have further demonstrated that miR-128 functions through a dual mechanism; by directly targeting L1 RNA for degradation and indirectly by inhibiting a cellular co-factor which L1 is dependent on to transpose to new genomic locations (TNPO1). Here, we add another piece to the puzzle of the enigmatic L1 lifecycle. We show that miR-128 also inhibits another key cellular factor, hnRNPA1 (heterogeneous nuclear ribonucleoprotein A1), by significantly reducing mRNA and protein levels through direct interaction with the coding sequence (CDS) of hnRNPA1 mRNA. In addition, we demonstrate that repression of hnRNPA1 using hnRNPA1-shRNA significantly decreases de novo L1 retro-transposition and that induced hnRNPA1 expression enhances L1 mobilization. Furthermore, we establish that hnRNPA1 is a functional target of miR-128. Finally, we determine that induced hnRNPA1 expression in miR-128-overexpressing cells can partly rescue the miR-128-induced repression of L1's ability to transpose to different genomic locations. Thus, we have identified an additional mechanism by which miR-128 represses L1 retro-transposition and mediates genomic stability.

A panoramic view of RNA modifications: exploring new frontiers.

Meeting report on the Cold Spring Harbor Asia conference on RNA Modifications and Epitranscriptomics, held in Suzhou, China, 13-17 November, 2017.

The Expanding Landscape of Alternative Splicing Variation in Human Populations.

Alternative splicing is a tightly regulated biological process by which the number of gene products for any given gene can be greatly expanded. Genomic variants in splicing regulatory sequences can disrupt splicing and cause disease. Recent developments in sequencing technologies and computational biology have allowed researchers to investigate alternative splicing at an unprecedented scale and resolution. Population-scale transcriptome studies have revealed many naturally occurring genetic variants that modulate alternative splicing and consequently influence phenotypic variability and disease susceptibility in human populations. Innovations in experimental and computational tools such as massively parallel reporter assays and deep learning have enabled the rapid screening of genomic variants for their causal impacts on splicing. In this review, we describe technological advances that have greatly increased the speed and scale at which discoveries are made about the genetic variation of alternative splicing. We summarize major findings from population transcriptomic studies of alternative splicing and discuss the implications of these findings for human genetics and medicine.

Population and allelic variation of A-to-I RNA editing in human transcriptomes.

BACKGROUND: A-to-I RNA editing is an important step in RNA processing in which specific adenosines in some RNA molecules are post-transcriptionally modified to inosines. RNA editing has emerged as a widespread mechanism for generating transcriptome diversity. However, there remain significant knowledge gaps about the variation and function of RNA editing. RESULTS: In order to determine the influence of genetic variation on A-to-I RNA editing, we integrate genomic and transcriptomic data from 445 human lymphoblastoid cell lines by combining an RNA editing QTL (edQTL) analysis with an allele-specific RNA editing (ASED) analysis. We identify 1054 RNA editing events associated with cis genetic polymorphisms. Additionally, we find that a subset of these polymorphisms is linked to genome-wide association study signals of complex traits or diseases. Finally, compared to random cis polymorphisms, polymorphisms associated with RNA editing variation are located closer spatially to their respective editing sites and have a more pronounced impact on RNA secondary structure. CONCLUSIONS: Our study reveals widespread cis variation in RNA editing among genetically distinct individuals and sheds light on possible phenotypic consequences of such variation on complex traits and diseases.

Mutations of human NARS2, encoding the mitochondrial asparaginyl-tRNA synthetase, cause nonsyndromic deafness and Leigh syndrome.

Here we demonstrate association of variants in the mitochondrial asparaginyl-tRNA synthetase NARS2 with human hearing loss and Leigh syndrome. A homozygous missense mutation ([c.637G>T; p.Val213Phe]) is the underlying cause of nonsyndromic hearing loss (DFNB94) and compound heterozygous mutations ([c.969T>A; p.Tyr323*] + [c.1142A>G; p.Asn381Ser]) result in mitochondrial respiratory chain deficiency and Leigh syndrome, which is a neurodegenerative disease characterized by symmetric, bilateral lesions in the basal ganglia, thalamus, and brain stem. The severity of the genetic lesions and their effects on NARS2 protein structure cosegregate with the phenotype. A hypothetical truncated NARS2 protein, secondary to the Leigh syndrome mutation p.Tyr323* is not detectable and p.Asn381Ser further decreases NARS2 protein levels in patient fibroblasts. p.Asn381Ser also disrupts dimerization of NARS2, while the hearing loss p.Val213Phe variant has no effect on NARS2 oligomerization. Additionally we demonstrate decreased steady-state levels of mt-tRNAAsn in fibroblasts from the Leigh syndrome patients. In these cells we show that a decrease in oxygen consumption rates (OCR) and electron transport chain (ETC) activity can be rescued by overexpression of wild type NARS2. However, overexpression of the hearing loss associated p.Val213Phe mutant protein in these fibroblasts cannot complement the OCR and ETC defects. Our findings establish lesions in NARS2 as a new cause for nonsyndromic hearing loss and Leigh syndrome.

Genome-wide identification and functional analysis of Apobec-1-mediated C-to-U RNA editing in mouse small intestine and liver.

BACKGROUND: RNA editing encompasses a post-transcriptional process in which the genomically templated sequence is enzymatically altered and introduces a modified base into the edited transcript. Mammalian C-to-U RNA editing represents a distinct subtype of base modification, whose prototype is intestinal apolipoprotein B mRNA, mediated by the catalytic deaminase Apobec-1. However, the genome-wide identification, tissue-specificity and functional implications of Apobec-1-mediated C-to-U RNA editing remain incompletely explored. RESULTS: Deep sequencing, data filtering and Sanger-sequence validation of intestinal and hepatic RNA from wild-type and Apobec-1-deficient mice revealed 56 novel editing sites in 54 intestinal mRNAs and 22 novel sites in 17 liver mRNAs, all within 3' untranslated regions. Eleven of 17 liver RNAs shared editing sites with intestinal RNAs, while 6 sites are unique to liver. Changes in RNA editing lead to corresponding changes in intestinal mRNA and protein levels for 11 genes. Analysis of RNA editing in vivo following tissue-specific Apobec-1 adenoviral or transgenic Apobec-1 overexpression reveals that a subset of targets identified in wild-type mice are restored in Apobec-1-deficient mouse intestine and liver following Apobec-1 rescue. We find distinctive polysome profiles for several RNA editing targets and demonstrate novel exonic editing sites in nuclear preparations from intestine but not hepatic apolipoprotein B RNA. RNA editing is validated using cell-free extracts from wild-type but not Apobec-1-deficient mice, demonstrating that Apobec-1 is required. CONCLUSIONS: These studies define selective, tissue-specific targets of Apobec-1-dependent RNA editing and show the functional consequences of editing are both transcript- and tissue-specific.

RNA editing in the human ENCODE RNA-seq data.

RNA-seq data can be mined for sequence differences relative to the reference genome to identify both genomic SNPs and RNA editing events. We analyzed the long, polyA-selected, unstranded, deeply sequenced RNA-seq data from the ENCODE Project across 14 human cell lines for candidate RNA editing events. On average, 43% of the RNA sequencing variants that are not in dbSNP and are within gene boundaries are A-to-G(I) RNA editing candidates. The vast majority of A-to-G(I) edits are located in introns and 3' UTRs, with only 123 located in protein-coding sequence. In contrast, the majority of non-A-to-G variants (60%-80%) map near exon boundaries and have the characteristics of splice-mapping artifacts. After filtering out all candidates with evidence of private genomic variation using genome resequencing or ChIP-seq data, we find that up to 85% of the high-confidence RNA variants are A-to-G(I) editing candidates. Genes with A-to-G(I) edits are enriched in Gene Ontology terms involving cell division, viral defense, and translation. The distribution and character of the remaining non-A-to-G variants closely resemble known SNPs. We find no reproducible A-to-G(I) edits that result in nonsynonymous substitutions in all three lymphoblastoid cell lines in our study, unlike RNA editing in the brain. Given that only a fraction of sites are reproducibly edited in multiple cell lines and that we find a stronger association of editing and specific genes suggests that the editing of the transcript is more important than the editing of any individual site.

The spreader flap in primary rhinoplasty.

BACKGROUND: In a primary rhinoplasty that requires a humpectomy, the dorsal aspect of the upper lateral cartilages is commonly discarded. Many of these patients need spreader grafts to reconstruct the middle third of the nose. However, it is possible to reconstruct the upper lateral cartilages into "spreader flaps" that act much like spreader grafts. METHODS: A tunnel is created on the underside of the upper lateral cartilage, which is released from the cartilaginous septum and also from its attachment to the nasal bone (medially). It is then rolled on itself to make a spreader flap, which is secured with sutures. Scoring along the dorsal edge of the upper lateral cartilage may be necessary. The flap is then secured to the dorsal edge of the reduced dorsal septum. RESULTS: In 21 patients who underwent an open approach (and four patients who underwent the closed approach), the spreader flap almost always reconstructed the middle third of the nose. It was easy to execute in the open approach but difficult in the closed approach. At surgery, two patients undergoing the open approach and one patient undergoing the closed approach needed spreader grafts because the flaps were too narrow. Postoperatively, only one patient (operated on by the open approach) exhibited inadequate nasal width. CONCLUSIONS: Spreader grafts are the standard for reconstructing the middle third of the nose. However, the spreader flap avoids harvesting and carving cartilage for those grafts. In the open approach, the technique is easy to execute. Conclusions could not be drawn regarding the long-term success with the closed approach.