Alu insertion variants alter mRNA splicing

Canine hyper-sociability structural variants associated with altered three-dimensional chromatin state

Strong selection on complex traits can lead to skewed trait means and reduced trait variability in populations. An example of this phenomenon can be evidenced in allele frequency changes and skewed trait distributions driven by persistent human-directed selective pressures in domesticated species. Dog domestication is linked to several genomic variants; however, the functional impacts of these variants may not always be straightforward when found in non-coding regions of the genome. Four polymorphic transposable elements (TE) found within non-coding sites along a 5 Mb region on canine CFA6 have evolved due to directional selection associated with heightened human-directed hyper-sociability in domesticated dogs. We found that the polymorphic TE in intron 17 of the canine GTF2I gene, which was previously reported to be negatively correlated with canid human-directed hyper-sociability, is associated with altered chromatin looping and hence distinct cis-regulatory landscapes. We reported supporting evidence of an E2F1-DNA binding peak concordant with the altered loop and higher expression of GTF2I exon 18, indicative of alternative splicing. Globally, we discovered differences in pathways regulating the extra-cellular matrix with respect to TE copy number. Overall, we reported evidence suggesting an intriguing molecular convergence between the emergence of hypersocial behaviors in dogs and the same genes that, when hemizygous, produce human Williams Beuren Syndrome characterized by cranio-facial defects and heightened social behaviors. Our results additionally emphasize the often-overlooked potential role of chromatin architecture in social evolution.

Uncovering hidden genetic variations: long-read sequencing reveals new insights into tuberous sclerosis complex

Background

Tuberous sclerosis is a multi-system disorder caused by mutations in either TSC1 or TSC2. The majority of affected patients (85%-90%) have heterozygous variants, and a smaller number (around 5%) have mosaic variants. Despite using various techniques, some patients still have "no mutation identified" (NMI).

Methods

We hypothesized that the causal variants of patients with NMI may be structural variants or deep intronic variants. To investigate this, we sequenced the DNA of 26 tuberous sclerosis patients with NMI using targeted long-read sequencing.

Results

We identified likely pathogenic/pathogenic variants in 13 of the cases, of which 6 were large deletions, four were InDels, two were deep intronic variants, one had retrotransposon insertion in either TSC1 or TSC2, and one was complex rearrangement. Furthermore, there was a de novo Alu element insertion with a high suspicion of pathogenicity that was classified as a variant of unknown significance.

Conclusion

Our findings expand the current knowledge of known pathogenic variants related to tuberous sclerosis, particularly uncovering mosaic complex structural variations and retrotransposon insertions that have not been previously reported in tuberous sclerosis. Our findings suggest a higher prevalence of mosaicism among tuberous sclerosis patients than previously recognized. Our results indicate that long-read sequencing is a valuable approach for tuberous sclerosis cases with no mutation identified (NMI).

KPNA2 promotes the progression of gastric cancer by regulating the alternative splicing of related genes

RNA-binding proteins (RBPs) play critical roles in genome regulation. In this study, we explored the latent function of KPNA2, which is an essential member of the RBP family, in the regulation of alternative splicing (AS) in gastric cancer (GC). We analyzed the role of KPNA2 in regulating differential expression and AS via RNA sequencing (RNA-seq) and improved RNA immunoprecipitation sequencing (iRIP-seq). Clinical specimens were used to analyze the associations between KPNA2 expression and clinicopathological characteristics. CCK8 assays, transwell assays and wound healing assays were performed to explore the effect of KPNA2/WDR62 on GC cell progression. KPNA2 was shown to be highly expressed in GC cells and tissues and associated with lymph node metastases. KPNA2 promoted the proliferation, migration and invasion of GC cells and primarily regulated exon skipping, alternative 3's splice sites (A3SSs), alternative 5' splice sites (A5SSs), and cassette exons. We further revealed that KPNA2 participated in biological processes related to cell proliferation, and the immune response in GC via the regulation of transcription. In addition, KPNA2 preferentially bound to intron regions. Notably, KPNA2 regulated the A3SS AS mode of WDR62, and upregulation of WDR62 reversed the KPNA2 downregulation-induced inhibition of GC cell proliferation, migration and invasion. Finally, we discovered that the AS of immune-related molecules could be regulated by KPNA2. Overall, our results demonstrated for the first time that KPNA2 functions as an oncogenic splicing factor in GC that regulated the AS and differential expression of GC-related genes, and KPNA2 may be a potential target for GC treatment.

Clinically relevant immune subtypes based on alternative splicing landscape of immune-related genes for lung cancer advanced PPPM approach

Background

Alternative splicing (AS) occurs in the process of gene post-transcriptional process, which is very important for the correct synthesis and function of protein. The change of AS pattern may lead to the change of expression level or function of lung cancer-related genes, and then affect the occurrence and development of lung cancers. The specific AS pattern might be used as a biomarker for early warning and prognostic assessment of a cancer in the framework of predictive, preventive, and personalized medicine (PPPM; 3PM). AS events of immune-related genes (IRGs) were closely associated with tumor progression and immunotherapy. We hypothesize that IRG-AS events are significantly different in lung adenocarcinomas (LUADs) vs. controls or in lung squamous cell carcinomas (LUSCs) vs. controls. IRG-AS alteration profiling was identified to construct IRG-differentially expressed AS (IRG-DEAS) signature models. Study on the selective AS events of specific IRGs in lung cancer patients might be of great significance for further exploring the pathogenesis of lung cancer, realizing early detection and effective monitoring of lung cancer, finding new therapeutic targets, overcoming drug resistance, and developing more effective therapeutic strategies, and better used for the prediction, diagnosis, prevention, and personalized medicine of lung cancer.

Methods

The transcriptomic, clinical, and AS data of LUADs and LUSCs were downloaded from TCGA and its SpliceSeq databases. IRG-DEAS events were identified in LUAD and LUSC, followed by their functional characteristics, and overall survival (OS) analyses. OS-related IRG-DEAS prognostic models were constructed for LUAD and LUSC with Lasso regression, which were used to classify LUADs and LUSCs into low- and high-risk score groups. Furthermore, the immune cell distribution, immune-related scores, drug sensitivity, mutation status, and GSEA/GSVA status were analyzed between low- and high-risk score groups. Also, low- and high-immunity clusters and AS factor (SF)-OS-related-AS co-expression network and verification of cell function of CELF6 were analyzed in LUAD and LUSC.

Results

Comprehensive analysis of transcriptomic, clinical, and AS data of LUADs and LUSCs identified IRG-AS events in LUAD (n = 1607) and LUSC (n = 1656), including OS-related IRG-AS events in LUAD (n = 127) and LUSC (n = 105). A total of 66 IRG-DEAS events in LUAD and 89 IRG-DEAS events in LUSC were identified compared to controls. The overlapping analysis between IRG-DEASs and OS-related IRG-AS events revealed 14 OS-related IRG-DEAS events for LUAD and 16 OS-related IRG-DEAS events for LUSC, which were used to identify and optimize a 12-OS-related-IRG-DEAS signature prognostic model for LUAD and an 11-OS-related-IRG-DEAS signature prognostic model for LUSC. These two prognostic models effectively divided LUAD or LUSC samples into low- and high-risk score groups that were closely associated with OS, clinical characteristics, and tumor immune microenvironment, with significant gene sets and pathways enriched in the two groups. Moreover, weighted gene co-expression network (WGCNA) and nonnegative matrix factorization method (NMF) analyses identified four OS-relevant subtypes of LUAD and six OS-relevant subtypes of LUSC, and ssGSEA identified five immunity-relevant subtypes of LUAD and five immunity-relevant subtypes of LUSC. Interestingly, splicing factors-OS-related-AS network revealed hub molecule CELF6 was significantly related to the malignant phenotype in lung cancer cells.

Conclusions

This study established two reliable IRG-DEAS signature prognostic models and constructed interesting splicing factor-splicing event networks in LUAD and LUSC, which can be used to construct clinically relevant immune subtypes, patient stratification, prognostic prediction, and personalized medical services in the PPPM practice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-024-00366-4.

Inverted Alu repeats: friends or foes in the human transcriptome

Alu elements are highly abundant primate-specific short interspersed nuclear elements that account for ~10% of the human genome. Due to their preferential location in gene-rich regions, especially in introns and 3' UTRs, Alu elements can exert regulatory effects on the expression of both host and neighboring genes. When two Alu elements with inverse orientations are positioned in close proximity, their transcription results in the generation of distinct double-stranded RNAs (dsRNAs), known as inverted Alu repeats (IRAlus). IRAlus are key immunogenic self-dsRNAs and post-transcriptional cis-regulatory elements that play a role in circular RNA biogenesis, as well as RNA transport and stability. Recently, IRAlus dsRNAs have emerged as regulators of transcription and activators of Z-DNA-binding proteins. The formation and activity of IRAlus can be modulated through RNA editing and interactions with RNA-binding proteins, and misregulation of IRAlus has been implicated in several immune-associated disorders. In this review, we summarize the emerging functions of IRAlus dsRNAs, the regulatory mechanisms governing IRAlus activity, and their relevance in the pathogenesis of human diseases.

FAME4-associating YEATS2 knockdown impairs dopaminergic synaptic integrity and leads to seizure-like behaviours in Drosophila melanogaster

Familial adult myoclonus epilepsy (FAME) is a neurological disorder caused by a TTTTA/TTTCA intronic repeat expansion. FAME4 is one of the six types of FAME that results from the repeat expansion in the first intron of the gene YEATS2. Although the RNA toxicity is believed to be the primary mechanism underlying FAME, the role of genes where repeat expansions reside is still unclear, particularly in the case of YEATS2 in neurons. This study used Drosophila to explore the effects of reducing YEATS2 expression. Two pan-neuronally driven dsDNA were used for knockdown of Drosophila YEATS2 (dYEATS2), and the resulting molecular and behavioural outcomes were evaluated. Drosophila with reduced dYEATS2 expression exhibited decreased tolerance to acute stress, disturbed locomotion, abnormal social behaviour, and decreased motivated activity. Additionally, reducing dYEATS2 expression negatively affected tyrosine hydroxylase (TH) gene expression, resulting in decreased dopamine biosynthesis. Remarkably, seizure-like behaviours induced by knocking down dYEATS2 were rescued by the administration of L-DOPA. This study reveals a novel role of YEATS2 in neurons in regulating acute stress responses, locomotion, and complex behaviours, and suggests that haploinsufficiency of YEATS2 may play a role in FAME4.

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.

AAGGG repeat expansions trigger RFC1-independent synaptic dysregulation in human CANVAS Neurons

Cerebellar ataxia with neuropathy and vestibular areflexia syndrome (CANVAS) is a late onset, recessively inherited neurodegenerative disorder caused by biallelic, non-reference pentameric AAGGG(CCCTT) repeat expansions within the second intron of replication factor complex subunit 1 (RFC1). To investigate how these repeats cause disease, we generated CANVAS patient induced pluripotent stem cell (iPSC) derived neurons (iNeurons) and utilized calcium imaging and transcriptomic analysis to define repeat-elicited gain-of-function and loss-of-function contributions to neuronal toxicity. AAGGG repeat expansions do not alter neuronal RFC1 splicing, expression, or DNA repair pathway functions. In reporter assays, AAGGG repeats are translated into pentapeptide repeat proteins that selectively accumulate in CANVAS patient brains. However, neither these proteins nor repeat RNA foci were detected in iNeurons, and overexpression of these repeats in isolation did not induce neuronal toxicity. CANVAS iNeurons exhibit defects in neuronal development and diminished synaptic connectivity that is rescued by CRISPR deletion of a single expanded allele. These phenotypic deficits were not replicated by knockdown of RFC1 in control neurons and were not rescued by ectopic expression of RFC1. These findings support a repeat-dependent but RFC1-independent cause of neuronal dysfunction in CANVAS, with important implications for therapeutic development in this currently untreatable condition.

Proceedings of the inaugural Dark Genome Symposium: November 2022

In November 2022 the first Dark Genome Symposium was held in Boston, USA. The meeting was hosted by Rome Therapeutics and Enara Bio, two biotechnology companies working on translating our growing understanding of this vast genetic landscape into therapies for human disease. The spirit and ambition of the meeting was one of shared knowledge, looking to strengthen the network of researchers engaged in the field. The meeting opened with a welcome from Rosana Kapeller and Kevin Pojasek followed by a first session of field defining talks from key academics in the space. A series of panels, bringing together academia and industry views, were then convened covering a wide range of pertinent topics. Finally, Richard Young and David Ting gave their views on the future direction and promise for patient impact inherent in the growing understanding of the Dark Genome.

SARS-CoV-2 infection severity and mortality is modulated by repeat-mediated regulation of alternative splicing

Like single-stranded RNA viruses, SARS-CoV-2 hijacks the host transcriptional machinery for its own replication. Numerous traditional differential gene expression-based investigations have examined the diverse clinical symptoms caused by SARS-CoV-2 infection. The virus, on the other hand, also affects the host splicing machinery, causing host transcriptional dysregulation, which can lead to diverse clinical outcomes. Hence, in this study, we performed host transcriptome sequencing of 125 hospital-admitted COVID-19 patients to understand the transcriptomic differences between the severity sub-phenotypes of mild, moderate, severe, and mortality. We performed transcript-level differential expression analysis, investigated differential isoform usage, looked at the splicing patterns within the differentially expressed transcripts (DET), and elucidated the possible genome regulatory features. Our DTE analysis showed evidence of diminished transcript length and diversity as well as altered promoter site usage in the differentially expressed protein-coding transcripts in the COVID-19 mortality patients. We also investigated the potential mechanisms driving the alternate splicing and discovered a compelling differential enrichment of repeats in the promoter region and a specific enrichment of SINE (Alu) near the splicing sites of differentially expressed transcripts. These findings suggested a repeat-mediated plausible regulation of alternative splicing as a potential modulator of COVID-19 disease severity. In this work, we emphasize the role of scarcely elucidated functional role of alternative splicing in influencing COVID-19 disease severity sub-phenotypes, clinical outcomes, and its putative mechanism. IMPORTANCE The wide range of clinical symptoms reported during the COVID-19 pandemic inherently highlights the numerous factors that influence the progression and prognosis of SARS-CoV-2 infection. While several studies have investigated the host response and discovered immunological dysregulation during severe infection, most of them have the common theme of focusing only up to the gene level. Viruses, especially RNA viruses, are renowned for hijacking the host splicing machinery for their own proliferation, which inadvertently puts pressure on the host transcriptome, exposing another side of the host response to the pathogen challenge. Therefore, in this study, we examine host response at the transcript-level to discover a transcriptional difference that culminates in differential gene-level expression. Importantly, this study highlights diminished transcript diversity and possible regulation of transcription by differentially abundant repeat elements near the promoter region and splicing sites in COVID-19 mortality patients, which together with differentially expressed isoforms hold the potential to elaborate disease severity and outcome.

A Genome-Wide Screen for the Exonisation of Reference SINE-VNTR-Alus and Their Expression in CNS Tissues of Individuals with Amyotrophic Lateral Sclerosis

The hominid-specific retrotransposon SINE-VNTR-Alu (SVA) is a composite element that has contributed to the genetic variation between individuals and influenced genomic structure and function. SVAs are involved in modulating gene expression and splicing patterns, altering mRNA levels and sequences, and have been associated with the development of disease. We evaluated the genome-wide effects of SVAs present in the reference genome on transcript sequence and expression in the CNS of individuals with and without the neurodegenerative disorder Amyotrophic Lateral Sclerosis (ALS). This study identified SVAs in the exons of 179 known transcripts, several of which were expressed in a tissue-specific manner, as well as 92 novel exonisation events occurring in the motor cortex. An analysis of 65 reference genome SVAs polymorphic for their presence/absence in the ALS consortium cohort did not identify any elements that were significantly associated with disease status, age at onset, and survival. However, there were transcripts, such as transferrin and HLA-A, that were differentially expressed between those with or without disease, and expression levels were associated with the genotype of proximal SVAs. This study demonstrates the functional consequences of several SVA elements altering mRNA splicing patterns and expression levels in tissues of the CNS.

The Immunological Conundrum of Endogenous Retroelements

Our defenses against infection rely on the ability of the immune system to distinguish invading pathogens from self. This task is exceptionally challenging, if not seemingly impossible, in the case of retroviruses that have integrated almost seamlessly into the host. This review examines the limits of innate and adaptive immune responses elicited by endogenous retroviruses and other retroelements, the targets of immune recognition, and the consequences for host health and disease. Contrary to theoretical expectation, endogenous retroelements retain substantial immunogenicity, which manifests most profoundly when their epigenetic repression is compromised, contributing to autoinflammatory and autoimmune disease and age-related inflammation. Nevertheless, recent evidence suggests that regulated immune reactivity to endogenous retroelements is integral to immune system development and function, underpinning cancer immunosurveillance, resistance to infection, and responses to the microbiota. Elucidation of the interaction points with endogenous retroelements will therefore deepen our understanding of immune system function and contribution to disease.

Calling and Phasing of Single-Nucleotide and Structural Variants of the LDLR Gene Using Oxford Nanopore MinION

The LDLR locus has clinical significance for lipid metabolism, Mendelian familial hypercholesterolemia (FH), and common lipid metabolism-related diseases (coronary artery disease and Alzheimer's disease), but its intronic and structural variants are underinvestigated. The aim of this study was to design and validate a method for nearly complete sequencing of the LDLR gene using long-read Oxford Nanopore sequencing technology (ONT). Five PCR amplicons from LDLR of three patients with compound heterozygous FH were analyzed. We used standard workflows of EPI2ME Labs for variant calling. All rare missense and small deletion variants detected previously by massively parallel sequencing and Sanger sequencing were identified using ONT. One patient had a 6976 bp deletion (exons 15 and 16) that was detected by ONT with precisely located breakpoints between AluY and AluSx1. Trans-heterozygous associations between mutation c.530C>T and c.1054T>C, c.2141-966_2390-330del, and c.1327T>C, and between mutations c.1246C>T and c.940+3_940+6del of LDLR, were confirmed. We demonstrated the ability of ONT to phase variants, thereby enabling haplotype assignment for LDLR with personalized resolution. The ONT-based method was able to detect exonic variants with the additional benefit of intronic analysis in one run. This method can serve as an efficient and cost-effective tool for diagnosing FH and conducting research on extended LDLR haplotype reconstruction.

Retrotransposon insertions associated with risk of neurologic and psychiatric diseases

Transposable elements (TEs) are active in neuronal cells raising the question whether TE insertions contribute to risk of neuropsychiatric disease. While genome-wide association studies (GWAS) serve as a tool to discover genetic loci associated with neuropsychiatric diseases, unfortunately GWAS do not directly detect structural variants such as TEs. To examine the role of TEs in psychiatric and neurologic disease, we evaluated 17,000 polymorphic TEs and find 76 are in linkage disequilibrium with disease haplotypes (P < 10-6 ) defined by GWAS. From these 76 polymorphic TEs, we identify potentially causal candidates based on having insertions in genomic regions of regulatory chromatin and on having associations with altered gene expression in brain tissues. We show that lead candidate insertions have regulatory effects on gene expression in human neural stem cells altering the activity of a minimal promoter. Taken together, we identify 10 polymorphic TE insertions that are potential candidates on par with other variants for having a causal role in neurologic and psychiatric disorders.

Alternative splicing and genetic variation of mhc-e: implications for rhesus cytomegalovirus-based vaccines

Rhesus cytomegalovirus (RhCMV)-based vaccination against Simian Immunodeficiency virus (SIV) elicits MHC-E-restricted CD8+ T cells that stringently control SIV infection in ~55% of vaccinated rhesus macaques (RM). However, it is unclear how accurately the RM model reflects HLA-E immunobiology in humans. Using long-read sequencing, we identified 16 Mamu-E isoforms and all Mamu-E splicing junctions were detected among HLA-E isoforms in humans. We also obtained the complete Mamu-E genomic sequences covering the full coding regions of 59 RM from a RhCMV/SIV vaccine study. The Mamu-E gene was duplicated in 32 (54%) of 59 RM. Among four groups of Mamu-E alleles: three ~5% divergent full-length allele groups (G1, G2, G2_LTR) and a fourth monomorphic group (G3) with a deletion encompassing the canonical Mamu-E exon 6, the presence of G2_LTR alleles was significantly (p = 0.02) associated with the lack of RhCMV/SIV vaccine protection. These genomic resources will facilitate additional MHC-E targeted translational research.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

Transcriptome innovations in primates revealed by single-molecule long-read sequencing

Transcriptomic diversity greatly contributes to the fundamentals of disease, lineage-specific biology, and environmental adaptation. However, much of the actual isoform repertoire contributing to shaping primate evolution remains unknown. Here, we combined deep long- and short-read sequencing complemented with mass spectrometry proteomics in a panel of lymphoblastoid cell lines (LCLs) from human, three other great apes, and rhesus macaque, producing the largest full-length isoform catalog in primates to date. Around half of the captured isoforms are not annotated in their reference genomes, significantly expanding the gene models in primates. Furthermore, our comparative analyses unveil hundreds of transcriptomic innovations and isoform usage changes related to immune function and immunological disorders. The confluence of these evolutionary innovations with signals of positive selection and their limited impact in the proteome points to changes in alternative splicing in genes involved in immune response as an important target of recent regulatory divergence in primates.

3D genome, on repeat: Higher-order folding principles of the heterochromatinized repetitive genome

Nearly half of the human genome is comprised of diverse repetitive sequences ranging from satellite repeats to retrotransposable elements. Such sequences are susceptible to stepwise expansions, duplications, inversions, and recombination events which can compromise genome function. In this review, we discuss the higher-order folding mechanisms of compartmentalization and loop extrusion and how they shape, and are shaped by, heterochromatin. Using primarily mammalian model systems, we contrast mechanisms governing H3K9me3-mediated heterochromatinization of the repetitive genome and highlight emerging links between repetitive elements and chromatin folding.

Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells

Domestication of transposable elements (TEs) into functional cis-regulatory elements is a widespread phenomenon. However, the mechanisms behind why some TEs are co-opted as functional enhancers while others are not are underappreciated. SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where ~3,700 copies are annotated, nearly half of which are human-specific. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we leveraged CRISPR-interference (CRISPRi), computational and functional genomics to elucidate the genomic features that underlie SVA domestication into human stem-cell gene regulation. We found that ~750 SVAs are co-opted as functional cis-regulatory elements in human induced pluripotent stem cells. These SVAs are significantly closer to genes and harbor more transcription factor binding sites than non-co-opted SVAs. We show that a long DNA motif composed of flanking YY1/2 and OCT4 binding sites is enriched in the co-opted SVAs and that these two transcription factors bind consecutively on the TE sequence. We used CRISPRi to epigenetically repress active SVAs in stem cell-like NCCIT cells. Epigenetic perturbation of active SVAs strongly attenuated YY1/OCT4 binding and influenced neighboring gene expression. Ultimately, SVA repression resulted in ~3,000 differentially expressed genes, 131 of which were the nearest gene to an annotated SVA. In summary, we demonstrated that SVAs modulate human gene expression, and uncovered that location and sequence composition contribute to SVA domestication into gene regulatory networks.

SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where ~3,700 copies are annotated. Nearly half of the SVAs annotated in the human genome are exclusive to our species. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we filled this knowledge-gap by focusing on human induced pluripotent stem cells (iPSCs) and on a pluripotent-like cell line (NCCITs). Through the analysis of histone marks, gene expression profiles, and by means of genome editing (CRISPR-interference), we identified ~750 SVAs that work as enhancers and promoters in human pluripotent cells, and characterized a mechanism for SVA co-option involving the transcription factors OCT4 and YY1. With our CRISPR approach, we demonstrated that repressing the 750 active SVAs leads to alteration in the expression of ~3,000 genes.

An antisense Alu transposon insertion/deletion polymorphism of ALDH1A1 may functionally associate with Parkinson's disease

BACKGROUND

Aldehyde dehydrogenase 1 (encoded by ALDH1A1) has been shown to protect against Parkinson's disease (PD) by reducing toxic metabolites of dopamine. We herein revealed an antisense Alu element insertion/deletion polymorphism in intron 4 of ALDH1A1, and hypothesized that it might play a role in PD.  METHODS: A Han Chinese cohort comprising 488 PD patients and 515 controls was recruited to validate the Alu insertion/deletion polymorphism following a previous study of tag-single nucleotide polymorphisms, where rs7043217 was shown to be significantly associated with PD. Functional analyses of the Alu element insertion were performed.

RESULTS

The Alu element of ALDH1A1 was identified to be a variant of Yb8 subfamily and termed as Yb8c4. The antisense Yb8c4 insertion/deletion polymorphism (named asYb8c4^ins and asYb8c4^del, respectively) appeared to be in a complete linkage disequilibrium with rs7043217 and was validated to be significantly associated with PD susceptibility with asYb8c4^ins serving as a risk allele (P = 0.030, OR = 1.224, 95% CI = 1.020-1.470). Multiple functional analyses including ALDH1A1 mRNA expression in blood cells of carriers, and reporters of EGFP and luciferase showed that the asYb8c4^ins had a suppressive activity on gene transcription. Mechanistic explorations suggested that the asYb8c4^ins induced no changes in CpG methylation and mRNA splicing of ALDH1A1 and appeared no binding of transcription factors.

CONCLUSIONS

Our results consolidate an involvement of ALDH1 in PD pathogenesis. The asYb8c4 polymorphism may be a functional output of its linkage disequilibrium-linked single nucleotide polymorphisms.

Global chromosome rearrangement induced by CRISPR-Cas9 reshapes the genome and transcriptome of human cells

Chromosome rearrangement plays important roles in development, carcinogenesis and evolution. However, its mechanism and subsequent effects are not fully understood. Large-scale chromosome rearrangement has been performed in the simple eukaryote, wine yeast, but the relative research in mammalian cells remains at the level of individual chromosome rearrangement due to technical limitations. In this study, we used CRISPR-Cas9 to target the highly repetitive human endogenous retrotransposons, LINE-1 and Alu, resulting in a large number of DNA double-strand breaks in the chromosomes. While this operation killed the majority of the cells, we eventually obtained live cell groups. Karyotype analysis and genome re-sequencing proved that we have achieved global chromosome rearrangement (GCR) in human cells. The copy number variations of the GCR genomes showed typical patterns observed in tumor genomes. The ATAC-seq and RNA-seq further revealed that the epigenetic and transcriptomic landscapes were deeply reshaped by GCR. Gene expressions related to p53 pathway, DNA repair, cell cycle and apoptosis were greatly altered to facilitate the cell survival. Our study provided a new application of CRISPR-Cas9 and a practical approach for GCR in complex mammalian genomes.

Protein-coding tRNA sequences?

Transfer RNAs (tRNAs) are ancient molecules likely predating the translation machinery. These extremely conserved RNA molecules transfer amino acids to the ribosome for the synthesis of proteins encoded by mRNAs, but canonical tRNAs are not protein-coding RNAs. Surprisely, when virtually translated, I observed that peptides derived from tRNA sequences match thousands of protein entries in databases. The analysis of these sequences indicates that the vast majority of these tRNA-derived proteins are annotated as small hypothetical peptides, likely arising from sequencing, prediction and/or annotation errors. But life often surpasses fiction. Importantly, tRNA-encoded amino acid domains were also found embedded in large functional proteins. Phylogenetic analysis of representative tRNA-derived protein domains may provide new insights into the origin, plasticity, and evolution of protein-coding genes.

Transposable Elements and Human Diseases: Mechanisms and Implication in the Response to Environmental Pollutants

Transposable elements (TEs) are recognized as major players in genome plasticity and evolution. The high abundance of TEs in the human genome, especially the Alu and Long Interspersed Nuclear Element-1 (LINE-1) repeats, makes them responsible for the molecular origin of several diseases. This involves several molecular mechanisms that are presented in this review: insertional mutation, DNA recombination and chromosomal rearrangements, modification of gene expression, as well as alteration of epigenetic regulations. This literature review also presents some of the more recent and/or more classical examples of human diseases in which TEs are involved. Whether through insertion of LINE-1 or Alu elements that cause chromosomal rearrangements, or through epigenetic modifications, TEs are widely implicated in the origin of human cancers. Many other human diseases can have a molecular origin in TE-mediated chromosomal recombination or alteration of gene structure and/or expression. These diseases are very diverse and include hemoglobinopathies, metabolic and neurological diseases, and common diseases. Moreover, TEs can also have an impact on aging. Finally, the exposure of individuals to stresses and environmental contaminants seems to have a non-negligible impact on the epigenetic derepression and mobility of TEs, which can lead to the development of diseases. Thus, improving our knowledge of TEs may lead to new potential diagnostic markers of diseases.

An Alu insertion map of the Indian population: identification and analysis in 1021 genomes of the IndiGen project

Actively retrotransposing primate-specific Alu repeats display insertion-deletion (InDel) polymorphism through their insertion at new loci. In the global datasets, Indian populations remain under-represented and so do their Alu InDels. Here, we report the genomic landscape of Alu InDels from the recently released 1021 Indian Genomes (IndiGen) (available at https://clingen.igib.res.in/indigen). We identified 9239 polymorphic Alu insertions that include private (3831), rare (3974) and common (1434) insertions with an average of 770 insertions per individual. We achieved an 89% PCR validation of the predicted genotypes in 94 samples tested. About 60% of identified InDels are unique to IndiGen when compared to other global datasets; 23% of sites were shared with both SGDP and HGSVC; among these, 58% (1289 sites) were common polymorphisms in IndiGen. The insertions not only show a bias for genic regions, with a preference for introns but also for the associated genes showing enrichment for processes like cell morphogenesis and neurogenesis (P-value &lt; 0.05). Approximately, 60% of InDels mapped to genes present in the OMIM database. Finally, we show that 558 InDels can serve as ancestry informative markers to segregate global populations. This study provides a valuable resource for baseline Alu InDels that would be useful in population genomics.

Emerging roles for endogenous retroviruses in immune epigenetic regulation

In recent years, there has been significant progress toward understanding the transcriptional networks underlying mammalian immune responses, fueled by advances in regulatory genomic technologies. Epigenomic studies profiling immune cells have generated detailed genome-wide maps of regulatory elements that will be key to deciphering the regulatory networks underlying cellular immune responses and autoimmune disorders. Unbiased analyses of these genomic maps have uncovered endogenous retroviruses as an unexpected ally in the regulation of human immune systems. Despite their parasitic origins, studies are finding an increasing number of examples of retroviral sequences having been co-opted for beneficial immune function and regulation by the host cell. Here, we review how endogenous retroviruses have given rise to numerous regulatory elements that shape the epigenetic landscape of host immune responses. We will discuss the implications of these elements on the function, dysfunction, and evolution of innate immunity.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Transposable elements that have recently been mobile in the human genome

Background

Transposable elements (TE) comprise nearly half of the human genome and their insertions have profound effects to human genetic diversification and as well as disease. Despite their abovementioned significance, there is no consensus on the TE subfamilies that remain active in the human genome. In this study, we therefore developed a novel statistical test for recently mobile subfamilies (RMSs), based on patterns of overlap with > 100,000 polymorphic indels.

Results

Our analysis produced a catalogue of 20 high-confidence RMSs, which excludes many false positives in public databases. Intriguingly though, it includes HERV-K, an LTR subfamily previously thought to be extinct. The RMS catalogue is strongly enriched for contributions to germline genetic disorders (P = 1.1e-10), and thus constitutes a valuable resource for diagnosing disorders of unknown aetiology using targeted TE-insertion screens. Remarkably, RMSs are also highly enriched for somatic insertions in diverse cancers (P = 2.8e-17), thus indicating strong correlations between germline and somatic TE mobility. Using CRISPR/Cas9 deletion, we show that an RMS-derived polymorphic TE insertion increased the expression of RPL17, a gene associated with lower survival in liver cancer. More broadly, polymorphic TE insertions from RMSs were enriched near genes with allele-specific expression, suggesting widespread effects on gene regulation.

Conclusions

By using a novel statistical test we have defined a catalogue of 20 recently mobile transposable element subfamilies. We illustrate the gene regulatory potential of RMS-derived polymorphic TE insertions, using CRISPR/Cas9 deletion in vitro on a specific candidate, as well as by genome wide analysis of allele-specific expression. Our study presents novel insights into TE mobility and regulatory potential and provides a key resource for human disease genetics and population history studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08085-0.

Dynamic Variations of 3'UTR Length Reprogram the mRNA Regulatory Landscape

This paper concerns 3′-untranslated regions (3′UTRs) of mRNAs, which are non-coding regulatory platforms that control stability, fate and the correct spatiotemporal translation of mRNAs. Many mRNAs have polymorphic 3′UTR regions. Controlling 3′UTR length and sequence facilitates the regulation of the accessibility of functional effectors (RNA binding proteins, miRNAs or other ncRNAs) to 3′UTR functional boxes and motifs and the establishment of different regulatory landscapes for mRNA function. In this context, shortening of 3′UTRs would loosen miRNA or protein-based mechanisms of mRNA degradation, while 3′UTR lengthening would strengthen accessibility to these effectors. Alterations in the mechanisms regulating 3′UTR length would result in widespread deregulation of gene expression that could eventually lead to diseases likely linked to the loss (or acquisition) of specific miRNA binding sites. Here, we will review the mechanisms that control 3′UTR length dynamics and their alterations in human disorders. We will discuss, from a mechanistic point of view centered on the molecular machineries involved, the generation of 3′UTR variability by the use of alternative polyadenylation and cleavage sites, of mutually exclusive terminal alternative exons (exon skipping) as well as by the process of exonization of Alu cassettes to generate new 3′UTRs with differential functional features.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

Angiotensin-Converting Enzyme (ACE) 1 Gene Polymorphism and Phenotypic Expression of COVID-19 Symptoms

The renin–angiotensin–aldosterone system (RAAS) appears to play an important role in SARS-CoV-2 infection. Polymorphisms within the genes that control this enzymatic system are candidates for elucidating the pathogenesis of COVID-19, since COVID-19 is not only a pulmonary disease but also affects many organs and systems throughout the body in multiple ways. Most striking is the fact that ACE2, one of the major components of the RAAS, is a prerequisite for SARS-COV-2 infection. Recently, we and other groups reported an association between a polymorphism of the ACE1 gene (a homolog of ACE2) and the phenotypic expression of COVID-19, particularly in its severity. The ethnic difference in ACE1 insertion (I)/deletion (D) polymorphism seems to explain the apparent difference in mortality between the West and East Asia. The purpose of this review was to further evaluate the evidence linking ACE1 polymorphisms to COVID-19. We searched the Medline database (2019–2021) for reference citations of relevant articles and selected studies on the clinical outcome of COVID-19 related to ACE1 I/D polymorphism. Although the numbers of patients are not large enough yet, most available evidence supports the notion that the DD genotype adversely influences COVID-19 symptoms. Surprisingly, small studies conducted in several countries yielded opposite results, suggesting that the ACE1 II genotype is a risk factor. This contradictory result may be the case in certain geographic areas, especially in subgroups of patients. It may also be due to interactions with other genes or to yet unexplained biochemical mechanisms. According to our hypothesis, such candidates are genes that are functionally involved in the pathophysiology of COVID-19, can act in concert with the ACE1 DD genotype, and that show differences in their frequency between the West and East Asia. For this, we conducted research focusing on Alu-related genes. The current study on the ACE1 genotype will provide potentially new clues to the pathogenesis, treatment, and diagnosis of SARS-CoV-2 infections.

Primate-specific retrotransposons and the evolution of circadian networks in the human brain

The circadian rhythm of the human brain is attuned to sleep-wake cycles that entail global alterations in neuronal excitability. This periodicity involves a highly coordinated regulation of gene expression. A growing number of studies are documenting a fascinating connection between primate-specific retrotransposons (Alu elements) and key epigenetic regulatory processes in the primate brain. Collectively, these studies indicate that Alu elements embedded in the human neuronal genome mediate post-transcriptional processes that unite human-specific neuroepigenetic landscapes and circadian rhythm. Here, we review evidence linking Alu retrotransposon-mediated posttranscriptional pathways to circadian gene expression. We hypothesize that Alu retrotransposons participate in the organization of circadian brain function through multidimensional neuroepigenetic pathways. We anticipate that these pathways are closely tied to the evolution of human cognition and their perturbation contributes to the manifestation of human-specific neurological diseases. Finally, we address current challenges and accompanying opportunities in studying primate- and human-specific transposable elements.

A 192 bp ERV fragment insertion in the first intron of porcine TLR6 may act as an enhancer associated with the increased expressions of TLR6 and TLR1

Background

Toll-like receptors (TLRs) play important roles in building innate immune and inducing adaptive immune responses. Associations of the TLR genes polymorphisms with disease susceptibility, which are the basis of molecular breeding for disease resistant animals, have been reported extensively. Retrotransposon insertion polymorphisms (RIPs), as a new type of molecular markers developed recently, have great potential in population genetics and quantitative trait locus mapping. In this study, bioinformatic prediction combined with PCR-based amplification was employed to screen for RIPs in porcine TLR genes. Their population distribution was examined, and for one RIP the impact on gene activity and phenotype was further evaluated.

Results

Five RIPs, located at the 3' flank of TLR3, 5' flank of TLR5, intron 1 of TLR6, intron 1 of TLR7, and 3' flank of TLR8 respectively, were identified. These RIPs were detected in different breeds with an uneven distribution among them. By using the dual luciferase activity assay a 192 bp endogenous retrovirus (ERV) in the intron 1 of TLR6 was shown to act as an enhancer increasing the activities of TLR6 putative promoter and two mini-promoters. Furthermore, real-time quantitative polymerase chain reaction (qPCR) analysis revealed significant association (p < 0.05) of the ERV insertion with increased mRNA expression of TLR6, the neighboring gene TLR1, and genes downstream in the TLR signaling pathway such as MyD88 (Myeloid differentiation factor 88), Rac1 (Rac family small GTPase 1), TIRAP (TIR domain containing adaptor protein), Tollip (Toll interacting protein) as well as the inflammatory factors IL6 (Interleukin 6), IL8 (Interleukin 8), and TNFα (Tumor necrosis factor alpha) in tissues of 30 day-old piglet. In addition, serum IL6 and TNFα concentrations were also significantly upregulated by the ERV insertion (p < 0.05).

Conclusions

A total of five RIPs were identified in five different TLR loci. The 192 bp ERV insertion in the first intron of TLR6 was associated with higher expression of TLR6, TLR1, and several genes downstream in the signaling cascade. Thus, the ERV insertion may act as an enhancer affecting regulation of the TLR signaling pathways, and can be potentially applied in breeding of disease resistant animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00248-w.

The Simons Genome Diversity Project: A Global Analysis of Mobile Element Diversity

Ongoing retrotransposition of Alu, LINE-1, and SINE–VNTR–Alu elements generates diversity and variation among human populations. Previous analyses investigating the population genetics of mobile element insertions (MEIs) have been limited by population ascertainment bias or by relatively small numbers of populations and low sequencing coverage. Here, we use 296 individuals representing 142 global populations from the Simons Genome Diversity Project (SGDP) to discover and characterize MEI diversity from deeply sequenced whole-genome data. We report 5,742 MEIs not originally reported by the 1000 Genomes Project and show that high sampling diversity leads to a 4- to 7-fold increase in MEI discovery rates over the original 1000 Genomes Project data. As a result of negative selection, nonreference polymorphic MEIs are underrepresented within genes, and MEIs within genes are often found in the transcriptional orientation opposite that of the gene. Globally, 80% of Alu subfamilies predate the expansion of modern humans from Africa. Polymorphic MEIs show heterozygosity gradients that decrease from Africa to Eurasia to the Americas, and the number of MEIs found uniquely in a single individual are also distributed in this general pattern. The maximum fraction of MEI diversity partitioned among the seven major SGDP population groups (F_ST) is 7.4%, similar to, but slightly lower than, previous estimates and likely attributable to the diverse sampling strategy of the SGDP. Finally, we utilize these MEIs to extrapolate the primary Native American shared ancestry component to back to Asia and provide new evidence from genome-wide identical-by-descent genetic markers that add additional support for a southeastern Siberian origin for most Native Americans.

CD58 Immunobiology at a Glance

The glycoprotein CD58, also known as lymphocyte-function antigen 3 (LFA-3), is a costimulatory receptor distributed on a broad range of human tissue cells. Its natural ligand CD2 is primarily expressed on the surface of T/NK cells. The CD2-CD58 interaction is an important component of the immunological synapse (IS) that induces activation and proliferation of T/NK cells and triggers a series of intracellular signaling in T/NK cells and target cells, respectively, in addition to promoting cell adhesion and recognition. Furthermore, a soluble form of CD58 (sCD58) is also present in cellular supernatant in vitro and in local tissues in vivo. The sCD58 is involved in T/NK cell-mediated immune responses as an immunosuppressive factor by affecting CD2-CD58 interaction. Altered accumulation of sCD58 may lead to immunosuppression of T/NK cells in the tumor microenvironment, allowing sCD58 as a novel immunotherapeutic target. Recently, the crucial roles of costimulatory molecule CD58 in immunomodulation seem to be reattracting the interests of investigators. In particular, the CD2-CD58 interaction is involved in the regulation of antiviral responses, inflammatory responses in autoimmune diseases, immune rejection of transplantation, and immune evasion of tumor cells. In this review, we provide a comprehensive summary of CD58 immunobiology.

Survival-associated alternative splicing events interact with the immune microenvironment in stomach adenocarcinoma

BACKGROUND

Alternative splicing (AS) increases the diversity of mRNA during transcription; it might play a role in alteration of the immune microenvironment, which could influence the development of immunotherapeutic strategies against cancer.

AIM

To obtain the transcriptomic and clinical features and AS events in stomach adenocarcinoma (STAD) from the database. The overall survival data associated with AS events were used to construct a signature prognostic model for STAD.

METHODS

Differentially expressed immune-related genes were identified between subtypes on the basis of the prognostic model. In STAD, 2042 overall-survival-related AS events were significantly enriched in various pathways and influenced several cellular functions. Furthermore, the network of splicing factors and overall-survival-associated AS events indicated potential regulatory mechanisms underlying the AS events in STAD.

RESULTS

An eleven-AS-signature prognostic model (CD44|14986|ES, PPHLN1|21214|AT, RASSF4|11351|ES, KIAA1147|82046|AP, PPP2R5D|76200|ES, LOH12CR1|20507|ES, CDKN3|27569|AP, UBA52|48486|AD, CADPS|65499|AT, SRSF7| 53276|RI, and WEE1|14328|AP) was constructed and significantly related to STAD overall survival, immune cells, and cancer-related pathways. The differentially expressed immune-related genes between the high- and low-risk score groups were significantly enriched in cancer-related pathways.

CONCLUSION

This study provided an AS-related prognostic model, potential mechanisms for AS, and alterations in the immune microenvironment (immune cells, genes, and pathways) for future research in STAD.

Principles and Practical Considerations for the Analysis of Disease-Associated Alternative Splicing Events Using the Gateway Cloning-Based Minigene Vectors pDESTsplice and pSpliceExpress

Splicing is an important RNA processing step. Genetic variations can alter the splicing process and thereby contribute to the development of various diseases. Alterations of the splicing pattern can be examined by gene expression analyses, by computational tools for predicting the effects of genetic variants on splicing, and by splicing reporter minigene assays for studying alternative splicing events under defined conditions. The minigene assay is based on transient transfection of cells with a vector containing a genomic region of interest cloned between two constitutive exons. Cloning can be accomplished by the use of restriction enzymes or by site-specific recombination using Gateway cloning. The vectors pDESTsplice and pSpliceExpress represent two minigene systems based on Gateway cloning, which are available through the Addgene plasmid repository. In this review, we describe the features of these two splicing reporter minigene systems. Moreover, we provide an overview of studies in which determinants of alternative splicing were investigated by using pDESTsplice or pSpliceExpress. The studies were reviewed with regard to the investigated splicing regulatory events and the experimental strategy to construct and perform a splicing reporter minigene assay. We further elaborate on how analyses on the regulation of RNA splicing offer promising prospects for gaining important insights into disease mechanisms.

Biology of the mRNA Splicing Machinery and Its Dysregulation in Cancer Providing Therapeutic Opportunities

Dysregulation of messenger RNA (mRNA) processing—in particular mRNA splicing—is a hallmark of cancer. Compared to normal cells, cancer cells frequently present aberrant mRNA splicing, which promotes cancer progression and treatment resistance. This hallmark provides opportunities for developing new targeted cancer treatments. Splicing of precursor mRNA into mature mRNA is executed by a dynamic complex of proteins and small RNAs called the spliceosome. Spliceosomes are part of the supraspliceosome, a macromolecular structure where all co-transcriptional mRNA processing activities in the cell nucleus are coordinated. Here we review the biology of the mRNA splicing machinery in the context of other mRNA processing activities in the supraspliceosome and present current knowledge of its dysregulation in lung cancer. In addition, we review investigations to discover therapeutic targets in the spliceosome and give an overview of inhibitors and modulators of the mRNA splicing process identified so far. Together, this provides insight into the value of targeting the spliceosome as a possible new treatment for lung cancer.

Unannotated splicing regulatory elements in deep intron space

Deep intron space harbors a diverse array of splicing regulatory elements that cooperate with better-known exon-proximal elements to enforce proper tissue-specific and development-specific pre-mRNA processing. Many deep intron elements have been highly conserved through vertebrate evolution, yet remain poorly annotated in the human genome. Recursive splicing exons (RS-exons) and intraexons promote noncanonical, multistep resplicing pathways in long introns, involving transient intermediate structures that are greatly underrepresented in RNA-seq datasets. Decoy splice sites and decoy exons act at a distance to inhibit splicing catalysis at annotated splice sites, with functional consequences such as exon skipping and intron retention. RNA:RNA bridges can juxtapose distant sequences within or across introns to activate deep intron splicing enhancers and silencers, to loop out exons to be skipped, or to select one member of a mutually exclusive set of exons. Similarly, protein bridges mediated by interactions among transcript-bound RNA binding proteins (RBPs) can modulate splicing outcomes. Experimental disruption of deep intron elements serving any of these functions can abrogate normal splicing, strongly suggesting that natural mutations of deep intron elements can do likewise to cause human disease. Understanding noncanonical splicing pathways and discovering deep intron regulatory signals, many of which map hundreds to many thousands of nucleotides from annotated splice junctions, is of great academic interest for basic scientists studying alternative splicing mechanisms. Hopefully, this knowledge coupled with increased analysis of deep intron sequences will also have important medical applications, as better interpretation of deep intron mutations may reveal new disease mechanisms and suggest new therapies. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing.

The "missing heritability"-Problem in psychiatry: Is the interaction of genetics, epigenetics and transposable elements a potential solution?

Psychiatric disorders exhibit an enormous burden on the health care systems worldwide accounting for around one-third of years lost due to disability among adults. Their etiology is largely unknown and diagnostic classification is based on symptomatology and course of illness and not on objective biomarkers. Most psychiatric disorders are moderately to highly heritable. However, it is still unknown what mechanisms may explain the discrepancy between heritability estimates and the present data from genetic analysis. In addition to genetic differences also epigenetic modifications are considered as potentially relevant in the transfer of susceptibility to psychiatric diseases. Though, whether or not epigenetic alterations can be inherited for many generations is highly controversial. In the present article, we will critically summarize both the genetic findings and the results from epigenetic analyses, including also those of noncoding RNAs. We will argue that one possible solution to the "missing heritability" problem in psychiatry is a potential role of retrotransposons, the exploration of which is presently only in its beginnings.

Insertion of an Alu-like element in MLH1 intron 7 as a novel cause of Lynch syndrome

Background

Lynch Syndrome (LS) is caused by germline mutations in the DNA mismatch repair (MMR) genes with mutations in MLH1 accounting for ~40% of LS‐related alterations.

Methods

MSK‐IMPACT analysis was performed on peripheral blood from a patient with early‐ onset colorectal cancer. Subsequently PCR and sequencing was performed to characterize the insertion. Immunohistochemistry for MMR genes and MLH1 promoter methylation were analyzed on patient's tumor.

Results

MSK‐IMPACT germline testing revealed an insertion into c.588+8_588+9 of MLH1 intron 7. The insertion was further characterized as an AluSx‐like element with ~115 bp in length. Functional studies demonstrated that the AluSx‐like element led to complete disruption of mRNA splicing and probably resulted in transcriptional termination at the poly (A) region of the AluSx‐like insertion.

Conclusions

The insertion of a truncated AluSx like element into MLH1 intron 7 results in aberrant splicing and transcription, thereby causing Lynch syndrome. This study confirms that retrotransposon insertions may be an important mechanism for cancer predisposition.

Polymorphic mobile element insertions contribute to gene expression and alternative splicing in human tissues

BACKGROUND

Mobile elements are a major source of structural variants in the human genome, and some mobile elements can regulate gene expression and transcript splicing. However, the impact of polymorphic mobile element insertions (pMEIs) on gene expression and splicing in diverse human tissues has not been thoroughly studied. The multi-tissue gene expression and whole genome sequencing data generated by the Genotype-Tissue Expression (GTEx) project provide a great opportunity to systematically evaluate the role of pMEIs in regulating gene expression in human tissues.

RESULTS

Using the GTEx whole genome sequencing data, we identify 20,545 high-quality pMEIs from 639 individuals. Coupling pMEI genotypes with gene expression profiles, we identify pMEI-associated expression quantitative trait loci (eQTLs) and splicing quantitative trait loci (sQTLs) in 48 tissues. Using joint analyses of pMEIs and other genomic variants, pMEIs are predicted to be the potential causal variant for 3522 eQTLs and 3717 sQTLs. The pMEI-associated eQTLs and sQTLs show a high level of tissue specificity, and these pMEIs are enriched in the proximity of affected genes and in regulatory elements. Using reporter assays, we confirm that several pMEIs associated with eQTLs and sQTLs can alter gene expression levels and isoform proportions, respectively.

CONCLUSION

Overall, our study shows that pMEIs are associated with thousands of gene expression and splicing variations, indicating that pMEIs could have a significant role in regulating tissue-specific gene expression and transcript splicing. Detailed mechanisms for the role of pMEIs in gene regulation in different tissues will be an important direction for future studies.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Our Conflict with Transposable Elements and Its Implications for Human Disease

Our genome is a historic record of successive invasions of mobile genetic elements. Like other eukaryotes, we have evolved mechanisms to limit their propagation and minimize the functional impact of new insertions. Although these mechanisms are vitally important, they are imperfect, and a handful of retroelement families remain active in modern humans. This review introduces the intrinsic functions of transposons, the tactics employed in their restraint, and the relevance of this conflict to human pathology. The most straightforward examples of disease-causing transposable elements are germline insertions that disrupt a gene and result in a monogenic disease allele. More enigmatic are the abnormal patterns of transposable element expression in disease states. Changes in transposon regulation and cellular responses to their expression have implicated these sequences in diseases as diverse as cancer, autoimmunity, and neurodegeneration. Distinguishing their epiphenomenal from their pathogenic effects may provide wholly new perspectives on our understanding of disease.

TypeTE: a tool to genotype mobile element insertions from whole genome resequencing data

Alu retrotransposons account for more than 10% of the human genome, and insertions of these elements create structural variants segregating in human populations. Such polymorphic Alus are powerful markers to understand population structure, and they represent variants that can greatly impact genome function, including gene expression. Accurate genotyping of Alus and other mobile elements has been challenging. Indeed, we found that Alu genotypes previously called for the 1000 Genomes Project are sometimes erroneous, which poses significant problems for phasing these insertions with other variants that comprise the haplotype. To ameliorate this issue, we introduce a new pipeline - TypeTE - which genotypes Alu insertions from whole-genome sequencing data. Starting from a list of polymorphic Alus, TypeTE identifies the hallmarks (poly-A tail and target site duplication) and orientation of Alu insertions using local re-assembly to reconstruct presence and absence alleles. Genotype likelihoods are then computed after re-mapping sequencing reads to the reconstructed alleles. Using a high-quality set of PCR-based genotyping of >200 loci, we show that TypeTE improves genotype accuracy from 83% to 92% in the 1000 Genomes dataset. TypeTE can be readily adapted to other retrotransposon families and brings a valuable toolbox addition for population genomics.

Contribution of unfixed transposable element insertions to human regulatory variation

Thousands of unfixed transposable element (TE) insertions segregate in the human population, but little is known about their impact on genome function. Recently, a few studies associated unfixed TE insertions to mRNA levels of adjacent genes, but the biological significance of these associations, their replicability across cell types and the mechanisms by which they may regulate genes remain largely unknown. Here, we performed a TE-expression QTL analysis of 444 lymphoblastoid cell lines (LCL) and 289 induced pluripotent stem cells using a newly developed set of genotypes for 2743 polymorphic TE insertions. We identified 211 and 176 TE-eQTL acting in cis in each respective cell type. Approximately 18% were shared across cell types with strongly correlated effects. Furthermore, analysis of chromatin accessibility QTL in a subset of the LCL suggests that unfixed TEs often modulate the activity of enhancers and other distal regulatory DNA elements, which tend to lose accessibility when a TE inserts within them. We also document a case of an unfixed TE likely influencing gene expression at the post-transcriptional level. Our study points to broad and diverse cis-regulatory effects of unfixed TEs in the human population and underscores their plausible contribution to phenotypic variation. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Fine mapping and identification of serum urate loci in American Indians: The Strong Heart Family Study

While studies have reported genetic loci affecting serum urate (SU) concentrations, few studies have been conducted in minority populations. Our objective for this study was to identify genetic loci regulating SU in a multigenerational family-based cohort of American Indians, the Strong Heart Family Study (SHFS). We genotyped 162,718 single nucleotide polymorphisms (SNPs) in 2000 SHFS participants using an Illumina MetaboChip array. A genome-wide association analysis of SU was conducted using measured genotype analysis approach accounting for kinships in SOLAR, and meta-analysis in METAL. Our results showed strong association of SU with rs4481233, rs9998811, rs7696092 and rs13145758 (minor allele frequency (MAF) = 25–44%; P < 3 × 10⁻¹⁴) of solute carrier family 2, member 9 (SLC2A9) and rs41481455, rs2231142 and rs1481012 (MAF = 29%; p < 3 × 10⁻⁹) of ATP-binding cassette protein, subfamily G, member 2 (ABCG2). Carriers of G alleles of rs9998811, rs4148155 and rs1481012 and A alleles of rs4481233, rs7696092 and rs13145758 and rs2231142 had lower SU concentrations as compared to non-carriers. Genetic analysis of SU conditional on significant SLC2A9 and ABCG2 SNPs revealed new loci, nucleobindin 1 (NUCB1) and neuronal PAS domain protein 4 (NPAS4) (p <6× 10⁻⁶). To identify American Indian-specific SNPs, we conducted targeted sequencing of key regions of SLC2A9. A total of 233 SNPs were identified of which 89 were strongly associated with SU (p < 7.1 × 10⁻¹⁰) and 117 were American Indian specific. Analysis of key SNPs in cohorts of Mexican-mestizos, European, Indian and East Asian ancestries showed replication of common SNPs, including our lead SNPs. Our results demonstrate the association of SU with uric acid transporters in a minority population of American Indians and potential novel associations of SU with neuronal-related genes which warrant further investigation.

Disease modeling of core pre-mRNA splicing factor haploinsufficiency

The craniofacial disorder mandibulofacial dysostosis Guion-Almeida type is caused by haploinsufficiency of the U5 snRNP gene EFTUD2/SNU114. However, it is unclear how reduced expression of this core pre-mRNA splicing factor leads to craniofacial defects. Here we use a CRISPR-Cas9 nickase strategy to generate a human EFTUD2-knockdown cell line and show that reduced expression of EFTUD2 leads to diminished proliferative ability of these cells, increased sensitivity to endoplasmic reticulum (ER) stress and the mis-expression of several genes involved in the ER stress response. RNA-Seq analysis of the EFTUD2-knockdown cell line revealed transcriptome-wide changes in gene expression, with an enrichment for genes associated with processes involved in craniofacial development. Additionally, our RNA-Seq data identified widespread mis-splicing in EFTUD2-knockdown cells. Analysis of the functional and physical characteristics of mis-spliced pre-mRNAs highlighted conserved properties, including length and splice site strengths, of retained introns and skipped exons in our disease model. We also identified enriched processes associated with the affected genes, including cell death, cell and organ morphology and embryonic development. Together, these data support a model in which EFTUD2 haploinsufficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on gene expression, including altering the expression of ER stress response genes and genes involved in the development of the craniofacial region. The increased burden of unfolded proteins in the ER resulting from mis-splicing would exceed the capacity of the defective ER stress response, inducing apoptosis in cranial neural crest cells that would result in craniofacial abnormalities during development.