Recombination of repeat elements generates somatic complexity in human genomes.

Non-allelic recombination between homologous repetitive elements contributes to evolution and human genetic disorders. Here, we combine short- and long-DNA read sequencing of repeat elements with a new bioinformatics pipeline to show that somatic recombination of Alu and L1 elements is widespread in the human genome. Our analysis uncovers tissue-specific non-allelic homologous recombination hallmarks; moreover, we find that centromeres and cancer-associated genes are enriched for retroelements that may act as recombination hotspots. We compare recombination profiles in human-induced pluripotent stem cells and differentiated neurons and find that the neuron-specific recombination of repeat elements accompanies chromatin changes during cell-fate determination. Finally, we report that somatic recombination profiles are altered in Parkinson's and Alzheimer's disease, suggesting a link between retroelement recombination and genomic instability in neurodegeneration. This work highlights a significant contribution of the somatic recombination of repeat elements to genomic diversity in health and disease.

Breaching Self-Tolerance to Alu Duplex RNA Underlies MDA5-Mediated Inflammation.

Aberrant activation of innate immune receptors can cause a spectrum of immune disorders, such as Aicardi-Goutières syndrome (AGS). One such receptor is MDA5, a viral dsRNA sensor that induces antiviral immune response. Using a newly developed RNase-protection/RNA-seq approach, we demonstrate here that constitutive activation of MDA5 in AGS results from the loss of tolerance to cellular dsRNAs formed by Alu retroelements. While wild-type MDA5 cannot efficiently recognize Alu-dsRNAs because of its limited filament formation on imperfect duplexes, AGS variants of MDA5 display reduced sensitivity to duplex structural irregularities, assembling signaling-competent filaments on Alu-dsRNAs. Moreover, we identified an unexpected role of an RNA-rich cellular environment in suppressing aberrant MDA5 oligomerization, highlighting context dependence of self versus non-self discrimination. Overall, our work demonstrates that the increased efficiency of MDA5 in recognizing dsRNA comes at a cost of self-recognition and implicates a unique role of Alu-dsRNAs as virus-like elements that shape the primate immune system.

Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

Human ADAR1 Prevents Endogenous RNA from Triggering Translational Shutdown.

Type I interferon (IFN) is produced when host sensors detect foreign nucleic acids, but how sensors differentiate self from nonself nucleic acids, such as double-stranded RNA (dsRNA), is incompletely understood. Mutations in ADAR1, an adenosine-to-inosine editing enzyme of dsRNA, cause Aicardi-Goutières syndrome, an autoinflammatory disorder associated with spontaneous interferon production and neurologic sequelae. We generated ADAR1 knockout human cells to explore ADAR1 substrates and function. ADAR1 primarily edited Alu elements in RNA polymerase II (pol II)-transcribed mRNAs, but not putative pol III-transcribed Alus. During the IFN response, ADAR1 blocked translational shutdown by inhibiting hyperactivation of PKR, a dsRNA sensor. ADAR1 dsRNA binding and catalytic activities were required to fully prevent endogenous RNA from activating PKR. Remarkably, ADAR1 knockout neuronal progenitor cells exhibited MDA5 (dsRNA sensor)-dependent spontaneous interferon production, PKR activation, and cell death. Thus, human ADAR1 regulates sensing of self versus nonself RNA, allowing pathogen detection while avoiding autoinflammation.

hnRNPM protects against the dsRNA-mediated interferon response by repressing LINE-associated cryptic splicing.

RNA splicing is pivotal in post-transcriptional gene regulation, yet the exponential expansion of intron length in humans poses a challenge for accurate splicing. Here, we identify hnRNPM as an essential RNA-binding protein that suppresses cryptic splicing through binding to deep introns, maintaining human transcriptome integrity. Long interspersed nuclear elements (LINEs) in introns harbor numerous pseudo splice sites. hnRNPM preferentially binds at intronic LINEs to repress pseudo splice site usage for cryptic splicing. Remarkably, cryptic exons can generate long dsRNAs through base-pairing of inverted ALU transposable elements interspersed among LINEs and consequently trigger an interferon response, a well-known antiviral defense mechanism. Significantly, hnRNPM-deficient tumors show upregulated interferon-associated pathways and elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity by repressing cryptic splicing and suggest that targeting hnRNPM in tumors may be used to trigger an inflammatory immune response, thereby boosting cancer surveillance.

Alu transposable elements rewire enhancer-promoter network through RNA pairing.

A recent study by Liang et al.1 reveals that interacting enhancer RNAs (eRNAs) and promoter-transcribed upstream antisense RNAs (uaRNAs) can identify enhancer-promoter interactions. Complementary sequences within the interacting eRNAs and uaRNAs, predominantly Alu sequences, confer the specificity for eRNA-uaRNA pairing and hence enhancer-promoter recognition.

HIV-1 integration landscape during latent and active infection.

The barrier to curing HIV-1 is thought to reside primarily in CD4(+) T cells containing silent proviruses. To characterize these latently infected cells, we studied the integration profile of HIV-1 in viremic progressors, individuals receiving antiretroviral therapy, and viremic controllers. Clonally expanded T cells represented the majority of all integrations and increased during therapy. However, none of the 75 expanded T cell clones assayed contained intact virus. In contrast, the cells bearing single integration events decreased in frequency over time on therapy, and the surviving cells were enriched for HIV-1 integration in silent regions of the genome. Finally, there was a strong preference for integration into, or in close proximity to, Alu repeats, which were also enriched in local hotspots for integration. The data indicate that dividing clonally expanded T cells contain defective proviruses and that the replication-competent reservoir is primarily found in CD4(+) T cells that remain relatively quiescent.

On the genetic basis of tail-loss evolution in humans and apes.

The loss of the tail is among the most notable anatomical changes to have occurred along the evolutionary lineage leading to humans and to the 'anthropomorphous apes'1-3, with a proposed role in contributing to human bipedalism4-6. Yet, the genetic mechanism that facilitated tail-loss evolution in hominoids remains unknown. Here we present evidence that an individual insertion of an Alu element in the genome of the hominoid ancestor may have contributed to tail-loss evolution. We demonstrate that this Alu element-inserted into an intron of the TBXT gene7-9-pairs with a neighbouring ancestral Alu element encoded in the reverse genomic orientation and leads to a hominoid-specific alternative splicing event. To study the effect of this splicing event, we generated multiple mouse models that express both full-length and exon-skipped isoforms of Tbxt, mimicking the expression pattern of its hominoid orthologue TBXT. Mice expressing both Tbxt isoforms exhibit a complete absence of the tail or a shortened tail depending on the relative abundance of Tbxt isoforms expressed at the embryonic tail bud. These results support the notion that the exon-skipped transcript is sufficient to induce a tail-loss phenotype. Moreover, mice expressing the exon-skipped Tbxt isoform develop neural tube defects, a condition that affects approximately 1 in 1,000 neonates in humans10. Thus, tail-loss evolution may have been associated with an adaptive cost of the potential for neural tube defects, which continue to affect human health today.

Circular DNA in the human germline and its association with recombination.

Extrachromosomal circular DNA (eccDNA) is common in somatic tissue, but its existence and effects in the human germline are unexplored. We used microscopy, long-read DNA sequencing, and new analytic methods to document thousands of eccDNAs from human sperm. EccDNAs derived from all genomic regions and mostly contained a single DNA fragment, although some consisted of multiple fragments. The generation of eccDNA inversely correlates with the meiotic recombination rate, and chromosomes with high coding-gene density and Alu element abundance form the least eccDNA. Analysis of insertions in human genomes further indicates that eccDNA can persist in the human germline when the circular molecules reinsert themselves into the chromosomes. Our results suggest that eccDNA has transient and permanent effects on the germline. They explain how differences in the physical and genetic map might arise and offer an explanation of how Alu elements coevolved with genes to protect genome integrity against deleterious mutations producing eccDNA.

Complementary sequence-mediated exon circularization.

Exon circularization has been identified from many loci in mammals, but the detailed mechanism of its biogenesis has remained elusive. By using genome-wide approaches and circular RNA recapitulation, we demonstrate that exon circularization is dependent on flanking intronic complementary sequences. Such sequences and their distribution exhibit rapid evolutionary changes, showing that exon circularization is evolutionarily dynamic. Strikingly, exon circularization efficiency can be regulated by competition between RNA pairing across flanking introns or within individual introns. Importantly, alternative formation of inverted repeated Alu pairs and the competition between them can lead to alternative circularization, resulting in multiple circular RNA transcripts produced from a single gene. Collectively, exon circularization mediated by complementary sequences in human introns and the potential to generate alternative circularization products extend the complexity of mammalian posttranscriptional regulation.

Complementary Alu sequences mediate enhancer-promoter selectivity.

Enhancers determine spatiotemporal gene expression programs by engaging with long-range promoters1-4. However, it remains unknown how enhancers find their cognate promoters. We recently developed a RNA in situ conformation sequencing technology to identify enhancer-promoter connectivity using pairwise interacting enhancer RNAs and promoter-derived noncoding RNAs5,6. Here we apply this technology to generate high-confidence enhancer-promoter RNA interaction maps in six additional cell lines. Using these maps, we discover that 37.9% of the enhancer-promoter RNA interaction sites are overlapped with Alu sequences. These pairwise interacting Alu and non-Alu RNA sequences tend to be complementary and potentially form duplexes. Knockout of Alu elements compromises enhancer-promoter looping, whereas Alu insertion or CRISPR-dCasRx-mediated Alu tethering to unregulated promoter RNAs can create new loops to homologous enhancers. Mapping 535,404 noncoding risk variants back to the enhancer-promoter RNA interaction maps enabled us to construct variant-to-function maps for interpreting their molecular functions, including 15,318 deletions or insertions in 11,677 Alu elements that affect 6,497 protein-coding genes. We further demonstrate that polymorphic Alu insertion at the PTK2 enhancer can promote tumorigenesis. Our study uncovers a principle for determining enhancer-promoter pairing specificity and provides a framework to link noncoding risk variants to their molecular functions.

A-to-I editing of coding and non-coding RNAs by ADARs.

Adenosine deaminases acting on RNA (ADARs) convert adenosine to inosine in double-stranded RNA. This A-to-I editing occurs not only in protein-coding regions of mRNAs, but also frequently in non-coding regions that contain inverted Alu repeats. Editing of coding sequences can result in the expression of functionally altered proteins that are not encoded in the genome, whereas the significance of Alu editing remains largely unknown. Certain microRNA (miRNA) precursors are also edited, leading to reduced expression or altered function of mature miRNAs. Conversely, recent studies indicate that ADAR1 forms a complex with Dicer to promote miRNA processing, revealing a new function of ADAR1 in the regulation of RNA interference.

Direct competition between hnRNP C and U2AF65 protects the transcriptome from the exonization of Alu elements.

There are ~650,000 Alu elements in transcribed regions of the human genome. These elements contain cryptic splice sites, so they are in constant danger of aberrant incorporation into mature transcripts. Despite posing a major threat to transcriptome integrity, little is known about the molecular mechanisms preventing their inclusion. Here, we present a mechanism for protecting the human transcriptome from the aberrant exonization of transposable elements. Quantitative iCLIP data show that the RNA-binding protein hnRNP C competes with the splicing factor U2AF65 at many genuine and cryptic splice sites. Loss of hnRNP C leads to formation of previously suppressed Alu exons, which severely disrupt transcript function. Minigene experiments explain disease-associated mutations in Alu elements that hamper hnRNP C binding. Thus, by preventing U2AF65 binding to Alu elements, hnRNP C plays a critical role as a genome-wide sentinel protecting the transcriptome. The findings have important implications for human evolution and disease.

RNA editing underlies genetic risk of common inflammatory diseases.

A major challenge in human genetics is to identify the molecular mechanisms of trait-associated and disease-associated variants. To achieve this, quantitative trait locus (QTL) mapping of genetic variants with intermediate molecular phenotypes such as gene expression and splicing have been widely adopted1,2. However, despite successes, the molecular basis for a considerable fraction of trait-associated and disease-associated variants remains unclear3,4. Here we show that ADAR-mediated adenosine-to-inosine RNA editing, a post-transcriptional event vital for suppressing cellular double-stranded RNA (dsRNA)-mediated innate immune interferon responses5-11, is an important potential mechanism underlying genetic variants associated with common inflammatory diseases. We identified and characterized 30,319 cis-RNA editing QTLs (edQTLs) across 49 human tissues. These edQTLs were significantly enriched in genome-wide association study signals for autoimmune and immune-mediated diseases. Colocalization analysis of edQTLs with disease risk loci further pinpointed key, putatively immunogenic dsRNAs formed by expected inverted repeat Alu elements as well as unexpected, highly over-represented cis-natural antisense transcripts. Furthermore, inflammatory disease risk variants, in aggregate, were associated with reduced editing of nearby dsRNAs and induced interferon responses in inflammatory diseases. This unique directional effect agrees with the established mechanism that lack of RNA editing by ADAR1 leads to the specific activation of the dsRNA sensor MDA5 and subsequent interferon responses and inflammation7-9. Our findings implicate cellular dsRNA editing and sensing as a previously underappreciated mechanism of common inflammatory diseases.

TFIIIC Binding to Alu Elements Controls Gene Expression via Chromatin Looping and Histone Acetylation.

How repetitive elements, epigenetic modifications, and architectural proteins interact ensuring proper genome expression remains poorly understood. Here, we report regulatory mechanisms unveiling a central role of Alu elements (AEs) and RNA polymerase III transcription factor C (TFIIIC) in structurally and functionally modulating the genome via chromatin looping and histone acetylation. Upon serum deprivation, a subset of AEs pre-marked by the activity-dependent neuroprotector homeobox Protein (ADNP) and located near cell-cycle genes recruits TFIIIC, which alters their chromatin accessibility by direct acetylation of histone H3 lysine-18 (H3K18). This facilitates the contacts of AEs with distant CTCF sites near promoter of other cell-cycle genes, which also become hyperacetylated at H3K18. These changes ensure basal transcription of cell-cycle genes and are critical for their re-activation upon serum re-exposure. Our study reveals how direct manipulation of the epigenetic state of AEs by a general transcription factor regulates 3D genome folding and expression.

Alu insertion-mediated dsRNA structure formation with pre-existing Alu elements as a disease-causing mechanism.

We previously identified a homozygous Alu insertion variant (Alu_Ins) in the 3'-untranslated region (3'-UTR) of SPINK1 as the cause of severe infantile isolated exocrine pancreatic insufficiency. Although we established that Alu_Ins leads to the complete loss of SPINK1 mRNA expression, the precise mechanisms remained elusive. Here, we aimed to elucidate these mechanisms through a hypothesis-driven approach. Initially, we speculated that, owing to its particular location, Alu_Ins could independently disrupt mRNA 3' end formation and/or affect other post-transcriptional processes such as nuclear export and translation. However, employing a 3'-UTR luciferase reporter assay, Alu_Ins was found to result in only an ∼50% reduction in luciferase activity compared to wild type, which is insufficient to account for the severe pancreatic deficiency in the Alu_Ins homozygote. We then postulated that double-stranded RNA (dsRNA) structures formed between Alu elements, an upstream mechanism regulating gene expression, might be responsible. Using RepeatMasker, we identified two Alu elements within SPINK1's third intron, both oriented oppositely to Alu_Ins. Through RNAfold predictions and full-length gene expression assays, we investigated orientation-dependent interactions between these Alu repeats. We provide compelling evidence to link the detrimental effect of Alu_Ins to extensive dsRNA structures formed between Alu_Ins and pre-existing intronic Alu sequences, including the restoration of SPINK1 mRNA expression by aligning all three Alu elements in the same orientation. Given the widespread presence of Alu elements in the human genome and the potential for new Alu insertions at almost any locus, our findings have important implications for detecting and interpreting Alu insertions in disease genes.

Human transposon tectonics.

Mobile DNAs have had a central role in shaping our genome. More than half of our DNA is comprised of interspersed repeats resulting from replicative copy and paste events of retrotransposons. Although most are fixed, incapable of templating new copies, there are important exceptions to retrotransposon quiescence. De novo insertions cause genetic diseases and cancers, though reliably detecting these occurrences has been difficult. New technologies aimed at uncovering polymorphic insertions reveal that mobile DNAs provide a substantial and dynamic source of structural variation. Key questions going forward include how and how much new transposition events affect human health and disease.

DICER1 loss and Alu RNA induce age-related macular degeneration via the NLRP3 inflammasome and MyD88.

Alu RNA accumulation due to DICER1 deficiency in the retinal pigmented epithelium (RPE) is implicated in geographic atrophy (GA), an advanced form of age-related macular degeneration that causes blindness in millions of individuals. The mechanism of Alu RNA-induced cytotoxicity is unknown. Here we show that DICER1 deficit or Alu RNA exposure activates the NLRP3 inflammasome and triggers TLR-independent MyD88 signaling via IL18 in the RPE. Genetic or pharmacological inhibition of inflammasome components (NLRP3, Pycard, Caspase-1), MyD88, or IL18 prevents RPE degeneration induced by DICER1 loss or Alu RNA exposure. These findings, coupled with our observation that human GA RPE contains elevated amounts of NLRP3, PYCARD, and IL18 and evidence of increased Caspase-1 and MyD88 activation, provide a rationale for targeting this pathway in GA. Our findings also reveal a function of the inflammasome outside the immune system and an immunomodulatory action of mobile elements.

lncRNA BC200 is processed into a stable Alu monomer.

The noncoding RNA BC200 is elevated in human cancers and is implicated in translation regulation as well as cell survival and proliferation. Upon BC200 overexpression, we observed correlated expression of a second, smaller RNA species. This RNA is expressed endogenously and exhibits cell-type-dependent variability relative to BC200. Aptamer-tagged expression constructs confirmed that the RNA is a truncated form of BC200, and sequencing revealed a modal length of 120 nt; thus, we refer to the RNA fragment as BC120. We present a methodology for accurate and specific detection of BC120 and establish that BC120 is expressed in several normal human tissues and is also elevated in ovarian cancer. BC120 exhibits remarkable stability relative to BC200 and is resistant to knockdown strategies that target the 3' unique sequence of BC200. Combined knockdown of BC200 and BC120 exhibits greater phenotypic impacts than knockdown of BC200 alone, and overexpression of BC120 negatively impacts translation of a GFP reporter, providing insight into a potential translational regulatory role for this RNA. The presence of a novel, truncated, and stable form of BC200 adds complexity to the investigation of this noncoding RNA that must be considered in future studies of BC200 and other related Alu RNAs.

Epigenetic therapy induces transcription of inverted SINEs and ADAR1 dependency.

Cancer therapies that target epigenetic repressors can mediate their effects by activating retroelements within the human genome. Retroelement transcripts can form double-stranded RNA (dsRNA) that activates the MDA5 pattern recognition receptor1-6. This state of viral mimicry leads to loss of cancer cell fitness and stimulates innate and adaptive immune responses7,8. However, the clinical efficacy of epigenetic therapies has been limited. To find targets that would synergize with the viral mimicry response, we sought to identify the immunogenic retroelements that are activated by epigenetic therapies. Here we show that intronic and intergenic SINE elements, specifically inverted-repeat Alus, are the major source of drug-induced immunogenic dsRNA. These inverted-repeat Alus are frequently located downstream of 'orphan' CpG islands9. In mammals, the ADAR1 enzyme targets and destabilizes inverted-repeat Alu dsRNA10, which prevents activation of the MDA5 receptor11. We found that ADAR1 establishes a negative-feedback loop, restricting the viral mimicry response to epigenetic therapy. Depletion of ADAR1 in patient-derived cancer cells potentiates the efficacy of epigenetic therapy, restraining tumour growth and reducing cancer initiation. Therefore, epigenetic therapies trigger viral mimicry by inducing a subset of inverted-repeats Alus, leading to an ADAR1 dependency. Our findings suggest that combining epigenetic therapies with ADAR1 inhibitors represents a promising strategy for cancer treatment.