Alu elements as regulators of gene expression

A transposable element prevents severe hemophilia B and provides insights into the evolution of new- and old world primates

Alu-elements comprise a large part of the human genome and some insertions have been shown to cause diseases. Here, we illuminate the protective role of an Alu-element in the 3'UTR of the human Factor 9 gene and its ability to ameliorate a poly(A) site mutation in a hemophilia B patient, preventing him from developing a severe disease. Using a minigene, we examined the disease-causing mutation and the modifying effect of the transposon in cellulo. Further, we simulated evolutionary scenarios regarding alternative polyadenylation before and after Alu insertion. A sequence analysis revealed that Old World monkeys displayed a highly conserved polyadenylation sites in this Alu-element, whereas New World monkeys lacked this motif, indicating a selective pressure. We conclude that this transposon has inserted shortly before the separation of Old and New World monkeys and thus also serves as a molecular landmark in primate evolution.

Detecting Alu Element Insertion Variant in RP1 Gene Using Whole Genome Sequencing in Patients with Retinitis Pigmentosa

Background/Objectives:Alu element insertion in the exon 4 of the RP1 gene was newly identified through whole genome sequencing (WGS). This was not detected in previous next-generation sequencing (NGS) analysis. We report three cases of Korean retinitis pigmentosa (RP) patients with compound heterozygous variants including Alu element insertion in the RP1 gene, indicating that Alu element insertion could be a cause of RP; Methods: Among patients diagnosed with RP having variants in the RP1 gene in the Asan Medical Center, WGS was additionally performed for genetically unsolved cases in previous NGS analysis to detect any presence of Alu element insertion. For cases detected to have Alu element insertion in the exon 4 of the RP1 gene, genetic and clinical characteristics were analyzed; Results: Among 16 patients with RP, 3 patients were detected to have Alu element insertion in the RP1 gene. Alu element insertion in the RP1 gene was also detected using WGS. It was revealed to be a pathogenic variant. Therefore, RP1 gene mutation was the confirmed genetic cause of RP for these three cases and genetic counseling was enabled for them; Conclusions: Alu element insertion in the RP1 gene could be a genetic cause of autosomal recessive RP patients with compound heterozygous variants. Through WGS, the identification of this pathogenic variant was possible. Confirmation is needed to check the presence of Alu element insertion in patients with compound heterozygous variants in the RP1 gene.

Enhancing in situ cancer vaccines using delivery technologies

In situ cancer vaccination refers to any approach that exploits tumour antigens available at a tumour site to induce tumour-specific adaptive immune responses. These approaches hold great promise for the treatment of many solid tumours, with numerous candidate drugs under preclinical or clinical evaluation and several products already approved. However, there are challenges in the development of effective in situ cancer vaccines. For example, inadequate release of tumour antigens from tumour cells limits antigen uptake by immune cells; insufficient antigen processing by antigen-presenting cells restricts the generation of antigen-specific T cell responses; and the suppressive immune microenvironment of the tumour leads to exhaustion and death of effector cells. Rationally designed delivery technologies such as lipid nanoparticles, hydrogels, scaffolds and polymeric nanoparticles are uniquely suited to overcome these challenges through the targeted delivery of therapeutics to tumour cells, immune cells or the extracellular matrix. Here, we discuss delivery technologies that have the potential to reduce various clinical barriers for in situ cancer vaccines. We also provide our perspective on this emerging field that lies at the interface of cancer vaccine biology and delivery technologies.

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).

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.

Massively parallel jumping assay decodes Alu retrotransposition activity

The human genome contains millions of retrotransposons, several of which could become active due to somatic mutations having phenotypic consequences, including disease. However, it is not thoroughly understood how nucleotide changes in retrotransposons affect their jumping activity. Here, we developed a novel massively parallel jumping assay (MPJA) that can test the jumping potential of thousands of transposons en masse. We generated nucleotide variant library of selected four Alu retrotransposons containing 165,087 different haplotypes and tested them for their jumping ability using MPJA. We found 66,821 unique jumping haplotypes, allowing us to pinpoint domains and variants vital for transposition. Mapping these variants to the Alu-RNA secondary structure revealed stem-loop features that contribute to jumping potential. Combined, our work provides a novel high-throughput assay that assesses the ability of retrotransposons to jump and identifies nucleotide changes that have the potential to reactivate them in the human genome.

A novel AluYb8 insertion-associated non-coding RNA, lncMUTYH, impairs mitochondrial function and dampens the M2-like polarization of macrophages

An inverted AluYb8 insertion in the MUTYH intron 15 (AluYb8MUTYH variant) has been reported to be associated with reduced MUTYH1 expression and mitochondrial dysfunction with age. However, the underlying mechanism remains unknown. In this study, we identified a novel transcript associated with the AluYb8MUTYH variant, which revealed that this transcript is about 780 nucleotides in length with a poly-A tail, lacks protein-coding potential, referred to as lncMUTYH. The results from the reporter gene system confirmed that the lncMUTYH down-regulates MUTYH1 expression at the translational level. Site-directed mutagenesis on the 5'-terminal exon sequences of α-MUTYH and lncMUTYH constructs revealed that lncMUTYH can act as a trans-regulator that depends on the partial base pairing between its exonized AluYb8 sequence and the 5'UTR of α-MUTYH to impede MUTYH 1 expression. Furthermore, we have demonstrated a correlation between decreased mitochondrion-localized MUTYH1 caused by lncMUTYH and lowered levels of mitochondrial biological function indicators, such as mtDNA content, mitochondrial regulatory gene expression, oxygen consumption rate, ATP product, and mitochondrial respiratory capacity. Notably, we found that lncMUTYH inhibited the M2-like polarization of macrophages, and CD68/CD206-positive cell fractions were significantly lower in lncMUTYH ectopically expressing cells. The results confirmed that the AluYb8MUTYH-associated lncMUTYH, derived from an AluYb8 insertion mutation, acts as a trans-regulatory factor that inhibits the MUTYH1 protein expression, leading to a progressive mitochondrial dysfunction that may disrupt macrophage differentiation. In summary, lncMUTYH can contribute to AluYb8MUTYH-associated mitochondrial dysfunction with age and hamper the macrophage polarization process, potentially increasing the risk of developing age-related diseases.

Molecular Mechanisms of Alu and LINE-1 Interspersed Repetitive Sequences Reveal Diseases of Visual System Dysfunction

BACKGROUND

Short interspersed nuclear elements (SINEs) and long interspersed nuclear elements (LINE-1s) are the abundant and well-characterized repetitive elements in the human genome.

METHODS

For this review, all relevant original research studies were assessed by searching electronic databases, including PubMed, Google Scholar, and Web of Science, by using relevant keywords. Accumulating evidence indicates that the disorder of gene expression regulated by these repetitive sequences is one of the causes of the diseases of visual system dysfunction, including retinal degenerations, glaucoma, retinitis punctata albescens, retinitis pigmentosa, geographic atrophy, and age-related macular degeneration, suggesting that SINEs and LINE-1s may have great potential implications in ophthalmology.

RESULTS

Alu elements belonging to the SINEs are present in more than one million copies, comprising 10% of the human genome.

CONCLUSION

This study offers recent advances in Alu and LINE-1 mechanisms in the development of eye diseases. The current study could advance our knowledge of the roles of SINEs and LINE-1s in the developing process of eye diseases, suggesting new diagnostic biomarkers, therapeutic strategies, and significant points for future studies.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

Exposure to smoky coal combustion emissions and leukocyte Alu retroelement copy number

Household air pollution (HAP) from indoor combustion of solid fuel is a global health burden that has been linked to multiple diseases including lung cancer. In Xuanwei, China, lung cancer rate for non-smoking women is among the highest in the world and largely attributed to high levels of polycyclic aromatic hydrocarbons (PAHs) that are produced from combustion of smoky (bituminous) coal. Alu retroelements, repetitive mobile DNA sequences that can somatically multiply and promote genomic instability have been associated with risk of lung cancer and diesel engine exhaust exposure. We conducted analyses for 160 non-smoking women in an exposure assessment study in Xuanwei, China with a repeat sample from 49 subjects. Quantitative PCR was used to measure Alu repeat copy number relative to albumin gene copy number (Alu/ALB ratio). Associations between clusters derived from predicted levels of 43 HAP constituents, 5-methylchrysene (5-MC), a PAH previously associated with lung cancer in Xuanwei and was selected a priori for analysis, and Alu repeats were analyzed using generalized estimating equations. A cluster of 31 PAHs reflecting current exposure was associated with increased Alu copy number (β:0.03 per standard deviation change; 95% confidence interval (CI):0.01,0.04; P-value = 2E-04). One compound within this cluster, 5-MC, was also associated with increased Alu copy number (P-value = 0.02). Our findings suggest that exposure to PAHs due to indoor smoky coal combustion may contribute to genomic instability. Additionally, our study provides further support for 5-MC as a prominent carcinogenic component of smoky coal emissions. Further studies are needed to replicate our findings.

Classification of Promoter Sequences from Human Genome

We have developed a new method for promoter sequence classification based on a genetic algorithm and the MAHDS sequence alignment method. We have created four classes of human promoters, combining 17,310 sequences out of the 29,598 present in the EPD database. We searched the human genome for potential promoter sequences (PPSs) using dynamic programming and position weight matrices representing each of the promoter sequence classes. A total of 3,065,317 potential promoter sequences were found. Only 1,241,206 of them were located in unannotated parts of the human genome. Every other PPS found intersected with either true promoters, transposable elements, or interspersed repeats. We found a strong intersection between PPSs and Alu elements as well as transcript start sites. The number of false positive PPSs is estimated to be 3 × 10-8 per nucleotide, which is several orders of magnitude lower than for any other promoter prediction method. The developed method can be used to search for PPSs in various eukaryotic genomes.

Identification of a psychiatric risk gene NISCH at 3p21.1 GWAS locus mediating dendritic spine morphogenesis and cognitive function

BACKGROUND

Schizophrenia and bipolar disorder (BD) are believed to share clinical symptoms, genetic risk, etiological factors, and pathogenic mechanisms. We previously reported that single nucleotide polymorphisms spanning chromosome 3p21.1 showed significant associations with both schizophrenia and BD, and a risk SNP rs2251219 was in linkage disequilibrium with a human specific Alu polymorphism rs71052682, which showed enhancer effects on transcriptional activities using luciferase reporter assays in U251 and U87MG cells.

METHODS

CRISPR/Cas9-directed genome editing, real-time quantitative PCR, and public Hi-C data were utilized to investigate the correlation between the Alu polymorphism rs71052682 and NISCH. Primary neuronal culture, immunofluorescence staining, co-immunoprecipitation, lentiviral vector production, intracranial stereotaxic injection, behavioral assessment, and drug treatment were used to examine the physiological impacts of Nischarin (encoded by NISCH).

RESULTS

Deleting the Alu sequence in U251 and U87MG cells reduced mRNA expression of NISCH, the gene locates 180 kb from rs71052682, and Hi-C data in brain tissues confirmed the extensive chromatin contacts. These data suggested that the genetic risk of schizophrenia and BD predicted elevated NISCH expression, which was also consistent with the observed higher NISCH mRNA levels in the brain tissues from psychiatric patients compared with controls. We then found that overexpression of NISCH resulted in a significantly decreased density of mushroom dendritic spines with a simultaneously increased density of thin dendritic spines in primary cultured neurons. Intriguingly, elevated expression of this gene in mice also led to impaired spatial working memory in the Y-maze. Given that Nischarin is the target of anti-hypertensive agents clonidine and tizanidine, which have shown therapeutic effects in patients with schizophrenia and patients with BD in preliminary clinical trials, we demonstrated that treatment with those antihypertensive drugs could reduce NISCH mRNA expression and rescue the impaired working memory in mice.

CONCLUSIONS

We identify a psychiatric risk gene NISCH at 3p21.1 GWAS locus influencing dendritic spine morphogenesis and cognitive function, and Nischarin may have potentials for future therapeutic development.

The RNA polymerase III-RIG-I axis in antiviral immunity and inflammation

Nucleic acid sensors survey subcellular compartments for atypical or mislocalized RNA or DNA, ultimately triggering innate immune responses. Retinoic acid-inducible gene-I (RIG-I) is part of the family of cytoplasmic RNA receptors that can detect viruses. A growing literature demonstrates that mammalian RNA polymerase III (Pol III) transcribes certain viral or cellular DNA sequences into immunostimulatory RIG-I ligands, which elicits antiviral or inflammatory responses. Dysregulation of the Pol III-RIG-I sensing axis can lead to human diseases including severe viral infection outcomes, autoimmunity, and tumor progression. Here, we summarize the newly emerging role of viral and host-derived Pol III transcripts in immunity and also highlight recent advances in understanding how mammalian cells prevent unwanted immune activation by these RNAs to maintain homeostasis.

Exonized Alu repeats in the 3'UTR of a CYP20A1_Alu-LT transcript act as a miRNA sponge

OBJECTIVE

Alu repeats have gained huge importance in the creation and modification of regulatory networks. We previously reported a unique isoform of human CYP20A1 i.e. CYP20A1_Alu-LT with 23 Alu repeats exonized in its 9 kb long 3'UTR with 4742 potential binding sites for 994 miRNAs. The role of this transcript was hypothesized as a potential miRNA sponge in primary neurons as its expression correlated with that of 380 genes having shared miRNA sites and enriched in neuro-coagulopathy. This study provides experimental evidence for the miRNA sponge activity of CYP20A1_Alu-LT in neuronal cell lines.

RESULTS

We studied the Alu-rich fragment of the CYP20A1_Alu-LT extended 3'UTR with > 10 binding sites for miR-619-5p and miR-3677-3p. Enrichment of the Alu-rich fragment with Ago2 confirmed miRNA association of this transcript. Cloning the fragment downstream of a reporter gene led to a 90% decrease in luciferase activity. Overexpression and knockdown studies revealed a positive correlation between the expression of CYP20A1_Alu-LT and miR-619-5p / miR-3677-3p target genes. GAP43, one of the key modulators of nerve regeneration, was significantly altered by the expression of CYP20A1_Alu-LT. This study, for the first time, provides evidence for a unique regulatory function of exonized Alu repeats as miRNA sponges.

RNA Polymerases I and III in development and disease

Ribosomes are macromolecular machines that are globally required for the translation of all proteins in all cells. Ribosome biogenesis, which is essential for cell growth, proliferation and survival, commences with transcription of a variety of RNAs by RNA Polymerases I and III. RNA Polymerase I (Pol I) transcribes ribosomal RNA (rRNA), while RNA Polymerase III (Pol III) transcribes 5S ribosomal RNA and transfer RNAs (tRNA) in addition to a wide variety of small non-coding RNAs. Interestingly, despite their global importance, disruptions in Pol I and Pol III function result in tissue-specific developmental disorders, with craniofacial anomalies and leukodystrophy/neurodegenerative disease being among the most prevalent. Furthermore, pathogenic variants in genes encoding subunits shared between Pol I and Pol III give rise to distinct syndromes depending on whether Pol I or Pol III function is disrupted. In this review, we discuss the global roles of Pol I and III transcription, the consequences of disruptions in Pol I and III transcription, disorders arising from pathogenic variants in Pol I and Pol III subunits, and mechanisms underpinning their tissue-specific phenotypes.

Detection and Characterization of a De Novo Alu Retrotransposition Event Causing NKX2-1-Related Disorder

BACKGROUND

Heterozygous NKX2-1 loss-of-function variants cause combinations of hyperkinetic movement disorders (MDs, particularly childhood-onset chorea), pulmonary dysfunction, and hypothyroidism. Mobile element insertions (MEIs) are potential disease-causing structural variants whose detection in routine diagnostics remains challenging.

OBJECTIVE

To establish the molecular diagnosis of two first-degree relatives with clinically suspected NKX2-1-related disorder who had negative NKX2-1 Sanger (SS), whole-exome (WES), and whole-genome (WGS) sequencing.

METHODS

The proband's WES was analyzed for MEIs. A candidate MEI in NKX2-1 underwent optimized SS after plasmid cloning. Functional studies exploring NKX2-1 haploinsufficiency at RNA and protein levels were performed.

RESULTS

A 347-bp AluYa5 insertion with a 65-bp poly-A tail followed by a 16-bp duplication of the pre-insertion wild-type sequence in exon 3 of NKX2-1 (ENST00000354822.7:c.556_557insAlu541_556dup) segregated with the disease phenotype.

CONCLUSIONS

We identified a de novo exonic AluYa5 insertion causing NKX2-1-related disorder in SS/WES/WGS-negative cases, suggesting that MEI analysis of short-read sequencing data or targeted long-read sequencing could unmask the molecular diagnosis of unsolved MD cases. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Z-DNA and Z-RNA: Methods-Past and Future

A quote attributed to Yogi Berra makes the observation that "It's tough to make predictions, especially about the future," highlighting the difficulties posed to an author writing a manuscript like the present. The history of Z-DNA shows that earlier postulates about its biology have failed the test of time, both those from proponents who were wildly enthusiastic in enunciating roles that till this day still remain elusive to experimental validation and those from skeptics within the larger community who considered the field a folly, presumably because of the limitations in the methods available at that time. If anything, the biological roles we now know for Z-DNA and Z-RNA were not anticipated by anyone, even when those early predictions are interpreted in the most favorable way possible. The breakthroughs in the field were made using a combination of methods, especially those based on human and mouse genetic approaches informed by the biochemical and biophysical characterization of the Zα family of proteins. The first success was with the p150 Zα isoform of ADAR1 (adenosine deaminase RNA specific), with insights into the functions of ZBP1 (Z-DNA-binding protein 1) following soon after from the cell death community. Just as the replacement of mechanical clocks by more accurate designs changed expectations about navigation, the discovery of the roles assigned by nature to alternative conformations like Z-DNA has forever altered our view of how the genome operates. These recent advances have been driven by better methodology and by better analytical approaches. This article will briefly describe the methods that were key to these discoveries and highlight areas where new method development is likely to further advance our knowledge.

cGAS inhibition alleviates Alu RNA-induced immune responses and cytotoxicity in retinal pigmented epithelium

Background

The degeneration of retinal pigmented epithelium (RPE) cells results in severe diseases, such as age-related macular degeneration (AMD) that causes blindness in millions of individuals.

Results

We report that targeting GMP-AMP (cGAMP) synthase (cGAS) alleviates Alu RNA-induced immune responses and cytotoxicity in RPE. We find that the deletion of cGAS in RPE inhibits the Alu RNA-stimulated interferon production. cGAS deficiency also protects RPE from cell death triggered by Alu RNA. Importantly, two natural chemicals, epigallocatechin gallate (EGCG) and resveratrol (RSVL), are effective in suppressing the immunogenic and cytotoxic effect of Alu RNA in RPE.

Conclusions

Our findings further demonstrate the crucial role of cGAS in the Alu RNA-induced RPE damage and present EGCG and RSVL as potential therapies for AMD and other RPE degeneration-related conditions.

Reduction of A-to-I RNA editing in the failing human heart regulates formation of circular RNAs

Alterations of RNA editing that affect the secondary structure of RNAs can cause human diseases. We therefore studied RNA editing in failing human hearts. Transcriptome sequencing showed that adenosine-to-inosine (A-to-I) RNA editing was responsible for 80% of the editing events in the myocardium. Failing human hearts were characterized by reduced RNA editing. This was primarily attributable to Alu elements in introns of protein-coding genes. In the failing left ventricle, 166 circRNAs were upregulated and 7 circRNAs were downregulated compared to non-failing controls. Most of the upregulated circRNAs were associated with reduced RNA editing in the host gene. ADAR2, which binds to RNA regions that are edited from A-to-I, was decreased in failing human hearts. In vitro, reduction of ADAR2 increased circRNA levels suggesting a causal effect of reduced ADAR2 levels on increased circRNAs in the failing human heart. To gain mechanistic insight, one of the identified upregulated circRNAs with a high reduction of editing in heart failure, AKAP13, was further characterized. ADAR2 reduced the formation of double-stranded structures in AKAP13 pre-mRNA, thereby reducing the stability of Alu elements and the circularization of the resulting circRNA. Overexpression of circAKAP13 impaired the sarcomere regularity of human induced pluripotent stem cell-derived cardiomyocytes. These data show that ADAR2 mediates A-to-I RNA editing in the human heart. A-to-I RNA editing represses the formation of dsRNA structures of Alu elements favoring canonical linear mRNA splicing and inhibiting the formation of circRNAs. The findings are relevant to diseases with reduced RNA editing and increased circRNA levels and provide insights into the human-specific regulation of circRNA formation.

RNA transcription and degradation of Alu retrotransposons depends on sequence features and evolutionary history

Alu elements are one of the most successful groups of RNA retrotransposons and make up 11% of the human genome with over 1 million individual loci. They are linked to genetic defects, increases in sequence diversity, and influence transcriptional activity. Still, their RNA metabolism is poorly understood yet. It is even unclear whether Alu elements are mostly transcribed by RNA Polymerase II or III. We have conducted a transcription shutoff experiment by α-amanitin and metabolic RNA labeling by 4-thiouridine combined with RNA fragmentation (TT-seq) and RNA-seq to shed further light on the origin and life cycle of Alu transcripts. We find that Alu RNAs are more stable than previously thought and seem to originate in part from RNA Polymerase II activity, as previous reports suggest. Their expression however seems to be independent of the transcriptional activity of adjacent genes. Furthermore, we have developed a novel statistical test for detecting the expression of quantitative trait loci in Alu elements that relies on the de Bruijn graph representation of all Alu sequences. It controls for both statistical significance and biological relevance using a tuned k-mer representation, discovering influential sequence features missed by regular motif search. In addition, we discover several point mutations using a generalized linear model, and motifs of interest, which also match transcription factor-binding motifs.

Genetic variants affecting chemical mediated skin immunotoxicity

The skin is an immune-competent organ and this function may be impaired by exposure to chemicals, which may ultimately result in immune-mediated dermal disorders. Interindividual variability to chemical-induced skin immune reactions is associated with intrinsic individual characteristics and their genomes. In the last 30-40 years, several genes influencing susceptibility to skin immune reactions were identified. The aim of this review is to provide information regarding common genetic variations affecting skin immunotoxicity. The polymorphisms selected for this review are related to xenobiotic-metabolizing enzymes (CYPA1 and CYPB1 genes), antioxidant defense (GSTM1, GSTT1, and GSTP1 genes), aryl hydrocarbon receptor signaling pathway (AHR and ARNT genes), skin barrier function transepidermal water loss (FLG, CASP14, and SPINK5 genes), inflammation (TNF, IL10, IL6, IL18, IL31, and TSLP genes), major histocompatibility complex (MHC) and neuroendocrine system peptides (CALCA, TRPV1, ACE genes). These genes present variants associated with skin immune responses and diseases, as well as variants associated with protecting skin immune homeostasis following chemical exposure. The molecular and association studies focusing on these genetic variants may elucidate their functional consequences and contribution in the susceptibility to skin immunotoxicity. Providing information on how genetic variations affect the skin immune system may reduce uncertainties in estimating chemical hazards/risks for human health in the future.

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 < 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.

Alu RNA Structural Features Modulate Immune Cell Activation and A-to-I Editing of Alu RNAs Is Diminished in Human Inflammatory Bowel Disease

Alu retrotransposons belong to the class of short interspersed nuclear elements (SINEs). Alu RNA is abundant in cells and its repetitive structure forms double-stranded RNAs (dsRNA) that activate dsRNA sensors and trigger innate immune responses with significant pathological consequences. Mechanisms to prevent innate immune activation include deamination of adenosines to inosines in dsRNAs, referred to as A-to-I editing, degradation of Alu RNAs by endoribonucleases, and sequestration of Alu RNAs by RNA binding proteins. We have previously demonstrated that widespread loss of Alu RNA A-to-I editing is associated with diverse human diseases including viral (COVID-19, influenza) and autoimmune diseases (multiple sclerosis). Here we demonstrate loss of A-to-I editing in leukocytes is also associated with inflammatory bowel diseases. Our structure-function analysis demonstrates that ability to activate innate immune responses resides in the left arm of Alu RNA, requires a 5’-PPP, RIG-I is the major Alu dsRNA sensor, and A-to-I editing disrupts both structure and function. Further, edited Alu RNAs inhibit activity of unedited Alu RNAs. Altering Alu RNA nucleotide sequence increases biological activity. Two classes of Alu RNAs exist, one class stimulates both IRF and NF-kB transcriptional activity and a second class only stimulates IRF transcriptional activity. Thus, Alu RNAs play important roles in human disease but may also have therapeutic potential.

Anti-aging and anti-oxidant activities of murine short interspersed nuclear element antisense RNA

Short interspersed nuclear elements (SINEs) play a key role in regulating gene expression, and SINE RNAs are involved in age-related diseases. We investigated the anti-aging effects of a genetically engineered murine SINE B1 antisense RNA (B1as RNA) and explored its mechanism of action in naturally senescent BALB/c (≥14 months) and moderately senscent C57BL/6N (≥9 months) mice. After tail vein injection, B1as RNA was available in the blood of mice for approximately 30 min, persisted for approximately 2-4 h in most detected tissues and persisted approximately 48 h in lungs. We found that treatment with B1as RNA improved stamina and promoted hair re-growth in aged mice. Treatment with B1as RNA also partially rescued the increase in mitochondrial DNA copy number in liver and spleen tissues observed in aged and moderately senescent mice. Finally, treatment with B1as RNA increased the activities of superoxide dismutase and glutathione peroxidase in aged and moderately senescent mice, reduced these animals' malondialdehyde and reactive oxygen species levels, and modulated the expression of several aging-associated genes, including Sirtuin 1, p21, p16Ink4a, p15Ink4b and p19Arf, and anti-oxidant genes (Sesn1 and Sesn 2). These data suggest that B1as RNA inhibits the aging process by enhancing antioxidant activity, promoting the scavenging of free radicals, and modulating the expression of aging-associated genes. This is the first report describing the anti-aging activity of SINE antisense RNA, which may serve as an effective nucleic acid drug for the treatment of age-related diseases.

Alu RNA and their roles in human disease states

Alu RNA are implicated in the poor prognosis of several human disease states. These RNA are transcription products of primate specific transposable elements called Alu elements. These elements are extremely abundant, comprising over 10% of the human genome, and 100 to 1000 cytoplasmic copies of Alu RNA per cell. Alu RNA do not have a single universal functional role aside from selfish self-propagation. Despite this, Alu RNA have been found to operate in a diverse set of translational and transcriptional mechanisms. This review will focus on the current knowledge of Alu RNA involved in human disease states and known mechanisms of action. Examples of Alu RNA that are transcribed in a variety of contexts such as introns, mature mRNA, and non-coding transcripts will be discussed. Past and present challenges in studying Alu RNA, and the future directions of Alu RNA in basic and clinical research will also be examined.

Complementary Square-Wave Voltammetry and LC-MS/MS Analysis to Elucidate Induced Damaged and Mutated Mitochondrial and Nuclear DNA from in Vivo Knockdown of the BRCA1 Gene in the Mouse Skeletal Muscle

Breast cancer 1 gene (BRCA1) DNA mutations impact skeletal muscle functions. Inducible skeletal muscle specific Brca1 homozygote knockout (Brca1KOsmi, KO) mice accumulate mitochondrial DNA (mtDNA) mutations resulting in loss of muscle quality.1 Complementary electrochemical andmass spectrometry analyses were utilized to rapidly assess mtDNA or nuclear DNA (nDNA) extracted directly from mouse skeletal muscles. Oxidative peak currents (Ip) from DNA immobilized layer by layer (LbL) were monitored using square-wave voltammetry (SWV) via Ru(bpy)32+ electrocatalysis. Ip significantly decreased (p < 0.05) for KO mtDNA compared to heterozygous KO (Het) or wild type (WT), indicative of decreases in the guanine content. nDNA Ip significantly increased in KO compared to WT (p < 0.05), suggesting an accumulation of damaged nDNA. Guanine or oxidatively damaged guanine content was monitored via appropriate m/z mass transitions using liquid chromatography-tandem mass spectroscopy (LC-MS/MS). Guanine in both KO mtDNA and nDNA was significantly lower, while oxidatively damaged guanine in KO nDNA was significantly elevated versus WT. These data demonstrate a loss of guanine content consistent with mtDNA mutation accumulation. Oxidative damage in KO nDNA suggests that repair processes associated with Brca1 are impacted. Overall, electrochemical and LC-MS/MS analysis can provide chemical-level answers to biological model phenotypic responses as a rapid and cost-effective analysis alternative to established assays.

DeepGRP: engineering a software tool for predicting genomic repetitive elements using Recurrent Neural Networks with attention

BACKGROUND

Repetitive elements contribute a large part of eukaryotic genomes. For example, about 40 to 50% of human, mouse and rat genomes are repetitive. So identifying and classifying repeats is an important step in genome annotation. This annotation step is traditionally performed using alignment based methods, either in a de novo approach or by aligning the genome sequence to a species specific set of repetitive sequences. Recently, Li (Bioinformatics 35:4408-4410, 2019) developed a novel software tool dna-brnn to annotate repetitive sequences using a recurrent neural network trained on sample annotations of repetitive elements.

RESULTS

We have developed the methods of dna-brnn further and engineered a new software tool DeepGRP. This combines the basic concepts of Li (Bioinformatics 35:4408-4410, 2019) with current techniques developed for neural machine translation, the attention mechanism, for the task of nucleotide-level annotation of repetitive elements. An evaluation on the human genome shows a 20% improvement of the Matthews correlation coefficient for the predictions delivered by DeepGRP, when compared to dna-brnn. DeepGRP predicts two additional classes of repeats (compared to dna-brnn) and is able to transfer repeat annotations, using RepeatMasker-based training data to a different species (mouse). Additionally, we could show that DeepGRP predicts repeats annotated in the Dfam database, but not annotated by RepeatMasker. DeepGRP is highly scalable due to its implementation in the TensorFlow framework. For example, the GPU-accelerated version of DeepGRP is approx. 1.8 times faster than dna-brnn, approx. 8.6 times faster than RepeatMasker and over 100 times faster than HMMER searching for models of the Dfam database.

CONCLUSIONS

By incorporating methods from neural machine translation, DeepGRP achieves a consistent improvement of the quality of the predictions compared to dna-brnn. Improved running times are obtained by employing TensorFlow as implementation framework and the use of GPUs. By incorporating two additional classes of repeats, DeepGRP provides more complete annotations, which were evaluated against three state-of-the-art tools for repeat annotation.

Haplotypic Associations and Differentiation of MHC Class II Polymorphic Alu Insertions at Five Loci With HLA-DRB1 Alleles in 12 Minority Ethnic Populations in China

The analysis of polymorphic variations in the human major histocompatibility complex (MHC) class II genomic region on the short-arm of chromosome 6 is a scientific enquiry to better understand the diversity in population structure and the effects of evolutionary processes such as recombination, mutation, genetic drift, demographic history, and natural selection. In order to investigate associations between the polymorphisms of HLA-DRB1 gene and recent Alu insertions (POALINs) in the HLA class II region, we genotyped HLA-DRB1 and five Alu loci (AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10), and determined their allele frequencies and haplotypic associations in 12 minority ethnic populations in China. There were 42 different HLA-DRB1 alleles for ethnic Chinese ranging from 12 alleles in the Jinuo to 28 in the Yugur with only DRB1^∗08:03, DRB1^∗09:01, DRB1^∗12:02, DRB1^∗14:01, DRB1^∗15:01, and DRB1^∗15:02 present in all ethnic groups. The POALINs varied in frequency between 0.279 and 0.514 for AluDPB2, 0 and 0.127 for AluDQA2, 0.777 and 0.995 for AluDQA1, 0.1 and 0.455 for AluDRB1 and 0.084 and 0.368 for AluORF10. By comparing the data of the five-loci POALIN in 13 Chinese ethnic populations (including Han-Yunnan published data) against Japanese and Caucasian published data, marked differences were observed between the populations at the allelic or haplotypic levels. Five POALIN loci were in significant linkage disequilibrium with HLA-DRB1 in different populations and AluDQA1 had the highest percentage association with most of the HLA-DRB1 alleles, whereas the nearby AluDRB1 indel was strongly haplotypic for only DRB1^∗01, DRB1^∗10, DRB1^∗15 and DRB1^∗16. There were 30 five-locus POALIN haplotypes inferred in all populations with H5 (no Alu insertions except for AluDQA1) and H21 (only AluDPB2 and AluDQA1 insertions) as the two predominant haplotypes. Neighbor joining trees and principal component analyses of the Alu and HLA-DRB1 polymorphisms showed that genetic diversity of these genomic markers is associated strongly with the population characteristics of language family, migration and sociality. This comparative study of HLA-DRB1 alleles and multilocus, lineage POALIN frequencies of Chinese ethnic populations confirmed that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic and evolutionary markers for the identification of allele and haplotype lineages and genetic variations within the same and/or different populations.

Elevated Alu retroelement copy number among workers exposed to diesel engine exhaust

BACKGROUND

Millions of workers worldwide are exposed to diesel engine exhaust (DEE), a known genotoxic carcinogen. Alu retroelements are repetitive DNA sequences that can multiply and compromise genomic stability. There is some evidence linking altered Alu repeats to cancer and elevated mortality risks. However, whether Alu repeats are influenced by environmental pollutants is unexplored. In an occupational setting with high DEE exposure levels, we investigated associations with Alu repeat copy number.

METHODS

A cross-sectional study of 54 male DEE-exposed workers from an engine testing facility and a comparison group of 55 male unexposed controls was conducted in China. Personal air samples were assessed for elemental carbon, a DEE surrogate, using NIOSH Method 5040. Quantitative PCR (qPCR) was used to measure Alu repeat copy number relative to albumin (Alb) single-gene copy number in leucocyte DNA. The unitless Alu/Alb ratio reflects the average quantity of Alu repeats per cell. Linear regression models adjusted for age and smoking status were used to estimate relations between DEE-exposed workers versus unexposed controls, DEE tertiles (6.1-39.0, 39.1-54.5 and 54.6-107.7 µg/m³) and Alu/Alb ratio.

RESULTS

DEE-exposed workers had a higher average Alu/Alb ratio than the unexposed controls (p=0.03). Further, we found a positive exposure-response relationship (p=0.02). The Alu/Alb ratio was highest among workers exposed to the top tertile of DEE versus the unexposed controls (1.12±0.08 SD vs 1.06±0.07 SD, p=0.01).

CONCLUSION

Our findings suggest that DEE exposure may contribute to genomic instability. Further investigations of environmental pollutants, Alu copy number and carcinogenesis are warranted.

Early Treatment and IL1RN Single-Nucleotide Polymorphisms Affect Response to Anakinra in Systemic Juvenile Idiopathic Arthritis

OBJECTIVE

To evaluate the impact of early treatment and IL1RN genetic variants on the response to anakinra in systemic juvenile idiopathic arthritis (JIA).

METHODS

Response to anakinra was defined as achievement of clinically inactive disease (CID) at 6 months without glucocorticoid treatment. Demographic, clinical, and laboratory characteristics of 56 patients were evaluated in univariate and multivariate analyses as predictors of response to treatment. Six single-nucleotide polymorphisms (SNPs) in the IL1RN gene, previously demonstrated to be associated with a poor response to anakinra, were genotyped by quantitative polymerase chain reaction (qPCR) or Sanger sequencing. Haplotype mapping was performed with Haploview software. IL1RN messenger RNA (mRNA) expression in whole blood from patients, prior to anakinra treatment initiation, was assessed by qPCR.

RESULTS

After 6 months of anakinra treatment, 73.2% of patients met the criteria for CID without receiving glucocorticoids. In the univariate analysis, the variable most strongly related to the response was disease duration from onset to initiation of anakinra treatment, with an optimal cutoff at 3 months (area under the curve 84.1%). Patients who started anakinra treatment ≥3 months after disease onset had an 8-fold higher risk of nonresponse at 6 months of treatment. We confirmed that the 6 IL1RN SNPs were inherited as a common haplotype. We found that homozygosity for ≥1 high-expression SNP correlated with higher IL1RN mRNA levels and was associated with a 6-fold higher risk of nonresponse, independent of disease duration.

CONCLUSION

Our findings on patients with systemic JIA confirm the important role of early interleukin-1 inhibition and suggest that genetic IL1RN variants predict nonresponse to therapy with anakinra.

Immunomodulation by epigenome alterations in Mycobacterium tuberculosis infection

Mycobacterium tuberculosis (MTB) has co-evolved with humans for decades and developed several mechanisms to evade host immunity. It can efficiently alter the host epigenome, thus playing a major role in immunomodulation by either activating or suppressing genes responsible for mounting an immune response against the pathogen. Epigenetic modifications such as DNA methylation and chromatin remodelling regulate gene expression and influence several cellular processes. The involvement of epigenetic factors in disease onset and development had been overlooked upon in comparison to genetic mutations. It is now believed that assessment of epigenetic changes hold great potential in diagnosis, prevention and treatment strategies for a wide range of diseases. In this review, we unravel the principles of epigenetics and the numerous ways by which MTB re-shapes the host epigenetic landscape as a strategy to overpower the host immune system for its survival and persistence.

The transcriptional trajectories of pluripotency and differentiation comprise genes with antithetical architecture and repetitive-element content

Background

Extensive molecular differences exist between proliferative and differentiated cells. Here, we conduct a meta-analysis of publicly available transcriptomic datasets from preimplantation and differentiation stages examining the architectural properties and content of genes whose abundance changes significantly across developmental time points.

Results

Analysis of preimplantation embryos from human and mouse showed that short genes whose introns are enriched in Alu (human) and B (mouse) elements, respectively, have higher abundance in the blastocyst compared to the zygote. These highly expressed genes encode ribosomal proteins or metabolic enzymes. On the other hand, long genes whose introns are depleted in repetitive elements have lower abundance in the blastocyst and include genes from signaling pathways. Additionally, the sequences of the genes that are differentially expressed between the blastocyst and the zygote contain distinct collections of pyknon motifs that differ between up- and down-regulated genes. Further examination of the genes that participate in the stem cell-specific protein interaction network shows that their introns are short and enriched in Alu (human) and B (mouse) elements. As organogenesis progresses, in both human and mouse, we find that the primarily short and repeat-rich expressed genes make way for primarily longer, repeat-poor genes. With that in mind, we used a machine learning-based approach to identify gene signatures able to classify human adult tissues: we find that the most discriminatory genes comprising these signatures have long introns that are repeat-poor and include transcription factors and signaling-cascade genes. The introns of widely expressed genes across human tissues, on the other hand, are short and repeat-rich, and coincide with those with the highest expression at the blastocyst stage.

Conclusions

Protein-coding genes that are characteristic of each trajectory, i.e., proliferation/pluripotency or differentiation, exhibit antithetical biases in their intronic and exonic lengths and in their repetitive-element content. While the respective human and mouse gene signatures are functionally and evolutionarily conserved, their introns and exons are enriched or depleted in organism-specific repetitive elements. We posit that these organism-specific repetitive sequences found in exons and introns are used to effect the corresponding genes’ regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-020-00928-8.

A Mechanism Leading to Changes in Copy Number Variations Affected by Transcriptional Level Might Be Involved in Evolution, Embryonic Development, Senescence, and Oncogenesis Mediated by Retrotransposons

In recent years, more and more evidence has emerged showing that changes in copy number variations (CNVs) correlated with the transcriptional level can be found during evolution, embryonic development, and oncogenesis. However, the underlying mechanisms remain largely unknown. The success of the induced pluripotent stem cell suggests that genome changes could bring about transformations in protein expression and cell status; conversely, genome alterations generated during embryonic development and senescence might also be the result of genome changes. With rapid developments in science and technology, evidence of changes in the genome affected by transcriptional level has gradually been revealed, and a rational and concrete explanation is needed. Given the preference of the HIV-1 genome to insert into transposons of genes with high transcriptional levels, we propose a mechanism based on retrotransposons facilitated by specific pre-mRNA splicing style and homologous recombination (HR) to explain changes in CNVs in the genome. This mechanism is similar to that of the group II intron that originated much earlier. Under this proposed mechanism, CNVs on genome are dynamically and spontaneously extended in a manner that is positively correlated with transcriptional level or contract as the cell divides during evolution, embryonic development, senescence, and oncogenesis, propelling alterations in them. Besides, this mechanism explains several critical puzzles in these processes. From evidence collected to date, it can be deduced that the message contained in genome is not just three-dimensional but will become four-dimensional, carrying more genetic information.

Neovascular Macular Degeneration: A Review of Etiology, Risk Factors, and Recent Advances in Research and Therapy

Neovascular age-related macular degeneration (exudative or wet AMD) is a prevalent, progressive retinal degenerative macular disease that is characterized by neovascularization of the choroid, mainly affecting the elderly population causing gradual vision impairment. Risk factors such as age, race, genetics, iris color, smoking, drinking, BMI, and diet all play a part in nvAMD's progression, with anti-vascular endothelial growth factor (anti-VEGF) therapy being the mainstay of treatment. Current therapeutic advancements slow the progression of the disease but do not cure or reverse its course. Newer therapies such as gene therapies, Rho-kinase inhibitors, and levodopa offer potential new targets for treatment.

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

Dicer1 promotes Aβ clearance via blocking B2 RNA-mediated repression of apolipoprotein E

Metabolism of β-amyloid is critical for healthy brain. Decreased clearance of β-amyloid is associated with ensued accumulation of amyloid peptide, culminating in formation of senile plaques, a neuropathological hallmark of Alzheimer's disease(AD). Apolipoprotein E (APOE), a lipoprotein for phospholipid and cholesterol metabolism, is predominantly synthesized by glia in the central nervous system, controlling Aβ aggregation and metabolism. By use of stereotactic injection and a Morris water maze, we found that delivery of Dicer1-expressing adenovirus into the hippocampus of an animal model of AD mice APPswe/PSEN1deltaE9 significantly improved spatial memory. The effect was associated with reduced amyloid peptides in the hippocampus which were analyzed with immunofluorescence and enzyme-linked immunosorbent assay. With western blot, quantitative real-time PCR, fluorescence in situ hybridization, and northern blot,Dicer1 overexpression increased apolipoprotein E (APOE) and concomitantly decreased B2 RNA in the hippocampus of the AD mice and in astrocyte cultures whereas transfection of B2 Mm2 RNA decreased APOE mRNA and protein levels in astrocyte cultures. Further, human or mouse APOE mRNA was found containing Alu RNA or its equivalent, B2 Mm2 RNA, locating downstream of its 3'-untranslated region (UTR), respectively. The 3'-UTR or 3'-UTR in conjunction with the downstream Alu/B2 RNA were cloned into a luciferase reporter; with dual-luciferase assay, we found that simultaneous transfection of Dicer1 siRNA or Alu/B2 RNA decreased the corresponding luciferase activities which suggest that Alu RNA mediated APOE mRNA degradation. Altogether, Dicer1 expression mediated amyloid peptide clearance by increasing APOE via blocking B2 RNA-mediated APOE mRNA degradation.

Multiple Alu Exonization in 3'UTR of a Primate-Specific Isoform of CYP20A1 Creates a Potential miRNA Sponge

Alu repeats contribute to phylogenetic novelties in conserved regulatory networks in primates. Our study highlights how exonized Alus could nucleate large-scale mRNA-miRNA interactions. Using a functional genomics approach, we characterize a transcript isoform of an orphan gene, CYP20A1 (CYP20A1_Alu-LT) that has exonization of 23 Alus in its 3'UTR. CYP20A1_Alu-LT, confirmed by 3'RACE, is an outlier in length (9 kb 3'UTR) and widely expressed. Using publically available data sets, we demonstrate its expression in higher primates and presence in single nucleus RNA-seq of 15,928 human cortical neurons. miRanda predicts ∼4,700 miRNA recognition elements (MREs) for ∼1,000 miRNAs, primarily originated within these 3'UTR-Alus. CYP20A1_Alu-LT could be a potential multi-miRNA sponge as it harbors ≥10 MREs for 140 miRNAs and has cytosolic localization. We further tested whether expression of CYP20A1_Alu-LT correlates with mRNAs harboring similar MRE targets. RNA-seq with conjoint miRNA-seq analysis was done in primary human neurons where we observed CYP20A1_Alu-LT to be downregulated during heat shock response and upregulated in HIV1-Tat treatment. In total, 380 genes were positively correlated with its expression (significantly downregulated in heat shock and upregulated in Tat) and they harbored MREs for nine expressed miRNAs which were also enriched in CYP20A1_Alu-LT. MREs were significantly enriched in these 380 genes compared with random sets of differentially expressed genes (P = 8.134e-12). Gene ontology suggested involvement of these genes in neuronal development and hemostasis pathways thus proposing a novel component of Alu-miRNA-mediated transcriptional modulation that could govern specific physiological outcomes in higher primates.

3'READS + RIP defines differential Staufen1 binding to alternative 3'UTR isoforms and reveals structures and sequence motifs influencing binding and polysome association

Staufen1 (STAU1) is an RNA-binding protein (RBP) that interacts with double-stranded RNA structures and has been implicated in regulating different aspects of mRNA metabolism. Previous studies have indicated that STAU1 interacts extensively with RNA structures in coding regions (CDSs) and 3'-untranslated regions (3'UTRs). In particular, duplex structures formed within 3'UTRs by inverted-repeat Alu elements (IRAlus) interact with STAU1 through its double-stranded RNA-binding domains (dsRBDs). Using 3' region extraction and deep sequencing coupled to ribonucleoprotein immunoprecipitation (3'READS + RIP), together with reanalyzing previous STAU1 binding and RNA structure data, we delineate STAU1 interactions transcriptome-wide, including binding differences between alternative polyadenylation (APA) isoforms. Consistent with previous reports, RNA structures are dominant features for STAU1 binding to CDSs and 3'UTRs. Overall, relative to short 3'UTR counterparts, longer 3'UTR isoforms of genes have stronger STAU1 binding, most likely due to a higher frequency of RNA structures, including specific IRAlus sequences. Nevertheless, a sizable fraction of genes express transcripts showing the opposite trend, attributable to AU-rich sequences in their alternative 3'UTRs that may recruit antagonistic RBPs and/or destabilize RNA structures. Using STAU1-knockout cells, we show that strong STAU1 binding to mRNA 3'UTRs generally enhances polysome association. However, IRAlus generally have little impact on STAU1-mediated polysome association despite having strong interactions with the protein. Taken together, our work reveals complex interactions of STAU1 with its cognate RNA substrates. Our data also shed light on distinct post-transcriptional fates for the widespread APA isoforms in mammalian cells.

DICER1 in the Pathogenesis of Age-related Macular Degeneration (AMD) - Alu RNA Accumulation versus miRNA Dysregulation

DICER1 deficiency in the retinal pigment epithelium (RPE) was associated with the accumulation of Alu transcripts and implicated in geographic atrophy (GA), a form of age-related macular degeneration (AMD), an eye disease leading to blindness in millions of people. Although the exact mechanism of this association is not fully known, the activation of the NLRP3 inflammasome, maturation of caspase-1 and disruption in mitochondrial homeostasis in RPE cells were shown as critical for it. DICER1 deficiency results in dysregulation of miRNAs and changes in the expression of many genes important for RPE homeostasis, which may also contribute to AMD. DICER1 deficiency can change the functions of the miR-183/96/182 cluster that regulates photoreceptors and their synaptic transmission. Aging, the main AMD risk factor, is associated with decreased expression of DICER1 and changes in its diurnal pattern that are not synchronized with circadian regulation in the retina. The initial insult inducing DICER1 deficiency in AMD may be oxidative stress, another major risk factor of AMD, but further studies on the role of deficient DICER1 in AMD pathogenesis and its therapeutic potential are needed.

The preparation of endotoxin-free genetically engineered murine B1 antisense RNA

One of the major limitations in the production of genetically engineered RNA from Escherichia coli (E. coli) is contamination by endotoxin. Here we report the first method that is capable of removing endotoxin from genetically engineered RNA. As a proof of concept, we transformed E. coli with a plasmid containing a tandem short interspersed nuclear elements from the mouse genome (SINE B1 elements). We then evaluated several extraction methods (SDS-NaCl centrifugation, SDS-NaCl filtration, TRIzol and SDS hot-phenol) and refinements thereof, and measured the resulting RNA yield, RNA purity, RNA integrity and endotoxin content. SDS-NaCl filtration with 2 mol/L NaCl, incorporating DEPC as an RNA protective agent, effectively removed endotoxin and resulted in a good RNA yield. Triton X-114 phase separation further reduced the endotoxin content of SDS-NaCl filtration-extracted RNA. RNA extracted by SDS-NaCl filtration with Triton X-114 phase separation did not cause adverse reactions in BALB/c mice and did not induce fever in rabbits when injected into these animals. The RNA met the requirements of nucleic acid reagents for in vivo experiments on animals.

Evolution of a splice variant that acts as an endogenous antagonist of the original INSL3 in primates

Alu sequences are the most abundant repetitive elements in the human genome, and have proliferated to more than one million copies in the human genome. Primate-specific Alu sequences account for ~10% of the human genome, and their spread within the genome has the potential to generate new exons. The new exons produced by Alu elements appear in various primate genes, and their functions have been elucidated. Here, we identified a new exon in the insulin-like 3 gene (INSL3), which evolved ~50 million years ago, and led to a splicing variant with 31 extra amino acid residues in addition to the original 95 nucleotides (NTs) of INSL3. The Alu-INSL3 isoform underwent diverse changes during primate evolution; we identified that human Alu-INSL3 might be on its way to functionality and has potential to antagonize LGR8-INSL3 function. Therefore, the present study is designed to provide an example of the evolutionary trajectory of a variant peptide hormone antagonist that caused by the insertion of an Alu element in primates.

A single mutation in the ACTR8 gene associated with lineage-specific expression in primates

Background

Alternative splicing (AS) generates various transcripts from a single gene and thus plays a significant role in transcriptomic diversity and proteomic complexity. Alu elements are primate-specific transposable elements (TEs) and can provide a donor or acceptor site for AS. In a study on TE-mediated AS, we recently identified a novel AluSz6-exonized ACTR8 transcript of the crab-eating monkey (Macaca fascicularis). In the present study, we sought to determine the molecular mechanism of AluSz6 exonization of the ACTR8 gene and investigate its evolutionary and functional consequences in the crab-eating monkey.

Results

We performed RT-PCR and genomic PCR to analyze AluSz6 exonization in the ACTR8 gene and the expression of the AluSz6-exonized transcript in nine primate samples, including prosimians, New world monkeys, Old world monkeys, and hominoids. AluSz6 integration was estimated to have occurred before the divergence of simians and prosimians. The Alu-exonized transcript obtained by AS was lineage-specific and expressed only in Old world monkeys and apes, and humans. This lineage-specific expression was caused by a single G duplication in AluSz6, which provides a new canonical 5′ splicing site. We further identified other alternative transcripts that were unaffected by the AluSz6 insertion. Finally, we observed that the alternative transcripts were transcribed into new isoforms with C-terminus deletion, and in silico analysis showed that these isoforms do not have a destructive function.

Conclusions

The single G duplication in the TE sequence is the source of TE exonization and AS, and this mutation may suffer a different fate of ACTR8 gene expression during primate evolution.

Identification of a functional human-unique 351-bp Alu insertion polymorphism associated with major depressive disorder in the 1p31.1 GWAS risk loci

Genome-wide association studies (GWAS) have reported substantial single-nucleotide polymorphisms (SNPs) associated with major depressive disorder (MDD), but the underlying functional variations in the GWAS risk loci are unclear. Here we show that the European MDD genome-wide risk-associated allele of rs12129573 at 1p31.1 is associated with MDD in Han Chinese, and this SNP is in strong linkage disequilibrium (LD) with a human-unique Alu insertion polymorphism (rs70959274) in the 5' flanking region of a long non-coding RNA (lncRNA) LINC01360 (Long Intergenic Non-Protein Coding RNA 1360), which is preferably expressed in human testis in the currently available expression datasets. The risk allele at rs12129573 is almost completely linked with the absence of this Alu insertion. The Alu insertion polymorphism (rs70959274) is significantly associated with a lower RNA level of LINC01360 and acts as a transcription silencer likely through modulating the methylation of its internal CpG sites. Luciferase assays confirm that the presence of Alu insertion at rs70959274 suppresses transcriptional activities in human cells, and deletion of the Alu insertion through CRISPR/Cas9-directed genome editing increases RNA expression of LINC01360. Deletion of the Alu insertion in human cells also leads to dysregulation of gene expression, biological processes and pathways relevant to MDD, such as the alterations of mRNA levels of DRD2 and FLOT1, transcription of genes involved in synaptic transmission, neurogenesis, learning or memory, and the PI3K-Akt signaling pathway. In summary, we identify a human-unique DNA repetitive polymorphism in robust LD with the MDD risk-associated SNP at the prominent 1p31.1 GWAS loci, and offer insights into the molecular basis of the illness.

The Developmental Gene Hypothesis for Punctuated Equilibrium: Combined Roles of Developmental Regulatory Genes and Transposable Elements

Theories of the genetics underlying punctuated equilibrium (PE) have been vague to date. Here the developmental gene hypothesis is proposed, which states that: 1) developmental regulatory (DevReg) genes are responsible for the orchestration of metazoan morphogenesis and their extreme conservation and mutation intolerance generates the equilibrium or stasis present throughout much of the fossil record and 2) the accumulation of regulatory elements and recombination within these same genes-often derived from transposable elements-drives punctuated bursts of morphological divergence and speciation across metazoa. This two-part hypothesis helps to explain the features that characterize PE, providing a theoretical genetic basis for the once-controversial theory. Also see the video abstract here https://youtu.be/C-fu-ks5yDs.

An Alu Element Insertion in Intron 1 Results in Aberrant Alternative Splicing of APOBEC3G Pre-mRNA in Northern Pig-Tailed Macaques (Macaca leonina) That May Reduce APOBEC3G-Mediated Hypermutation Pressure on HIV-1

APOBEC3 family members, particularly APOBEC3F and APOBEC3G, inhibit the replication and spread of various retroviruses by inducing hypermutation in newly synthesized viral DNA. Viral hypermutation by APOBEC3 is associated with viral evolution, viral transmission, and disease progression. In recent years, increasing attention has been paid to targeting APOBEC3G for AIDS therapy. Thus, a controllable model system using species such as macaques, which provide a relatively ideal in vivo system, is needed for the study of APOBEC3-related issues. To appropriately utilize this animal model for biomedical research, important differences between human and macaque APOBEC3s must be considered. In this study, we found that the ratio of APOBEC3G-mediated/APOBEC3-mediated HIV-1 hypermutation footprints was much lower in peripheral blood mononuclear cells (PBMCs) from northern pig-tailed macaques than in PBMCs from humans. Next, we identified a novel and conserved APOBEC3G pre-mRNA alternative splicing pattern in macaques, which differed from that in humans and resulted from an Alu element insertion into macaque APOBEC3G gene intron 1. This alternative splicing pattern generating an aberrant APOBEC3G mRNA isoform may significantly dilute full-length APOBEC3G and reduce APOBEC3G-mediated hypermutation pressure on HIV-1 in northern pig-tailed macaques, which was supported by the elimination of other possibilities accounting for this hypermutation difference between the two hosts.IMPORTANCE APOBEC3 family members, particularly APOBEC3F and APOBEC3G, are important cellular antiviral factors. Recently, more attention has been paid to targeting APOBEC3G for AIDS therapy. To appropriately utilize macaque animal models for the study of APOBEC3-related issues, it is important that the differences between human and macaque APOBEC3s are clarified. In this study, we identified a novel and conserved APOBEC3G pre-mRNA alternative splicing pattern in macaques, which differed from that in humans and which may reduce the APOBEC3G-mediated hypermutation pressure on HIV-1 in northern pig-tailed macaques (NPMs). Our work provides important information for the proper application of macaque animal models for APOBEC3-related issues in AIDS research and a better understanding of the biological functions of APOBEC3 proteins.

The genome-wide risk alleles for psychiatric disorders at 3p21.1 show convergent effects on mRNA expression, cognitive function, and mushroom dendritic spine

Schizophrenia and bipolar disorder (BPD) are believed to share clinical features, etiological factors, and disease pathologies (such as impaired cognitive functions and dendritic spine pathology). Meanwhile, there is growing evidence of shared genetic risk between schizophrenia and BPD, despite that our knowledge of the functional risk variations and biological mechanisms is still limited. Here, we conduct summary data-based Mendelian randomization (SMR) analyses through combining the statistical data from genome-wide association studies (GWAS) of both schizophrenia and BPD and multiple expression quantitative trait loci (eQTL) datasets of the human brain dorsolateral prefrontal cortex (DLPFC) tissues. These integrative investigations identify a lead risk locus at the chromosome 3p21.1 region, which contains numerous single-nucleotide polymorphisms (SNPs) in varied linkage disequilibrium (LD) and encompasses more than 20 genes. Further analyses suggest that many SNPs at 3p21.1 are significantly associated with both schizophrenia and BPD, and even depression, and the psychiatric risk alleles at 3p21.1 are correlated with mRNA expression of multiple genes such as NEK4, GNL3, and PBRM1. We also identify a 335-bp functional Alu polymorphism rs71052682 in significant LD with the psychiatric GWAS risk SNP rs2251219, and confirm the regulatory effects of this Alu polymorphism on transcription activities. We then explore the involvement of the 3p21.1 locus in the common clinical features and etiology of these illnesses. We reveal that psychiatric risk alleles at 3p21.1 in low-to-high LD consistently predict worse cognitive functions in humans, and manipulating the gene expression (NEK4, GNL3, and PBRM1) linked with higher genetic risk could reduce the density of mushroom dendritic spines in rat primary cortical neurons, mirroring the spine pathology in the prefrontal cortex of psychiatric patients. Our results find that, although the risk alleles at 3p21.1 are in low-to-moderate LD spanning a large genomic area, their underlying biological mechanisms in psychiatric disorders likely converge. These results provide essential insights into the neural mechanisms underlying the chromosome 3p21.1 risk locus in the shared pathological and etiological features of both schizophrenia and BPD.

Chromatin organization modulates the origin of heritable structural variations in human genome

Structural variations (SVs) in the human genome originate from different mechanisms related to DNA repair, replication errors, and retrotransposition. Our analyses of 26 927 SVs from the 1000 Genomes Project revealed differential distributions and consequences of SVs of different origin, e.g. deletions from non-allelic homologous recombination (NAHR) are more prone to disrupt chromatin organization while processed pseudogenes can create accessible chromatin. Spontaneous double stranded breaks (DSBs) are the best predictor of enrichment of NAHR deletions in open chromatin. This evidence, along with strong physical interaction of NAHR breakpoints belonging to the same deletion suggests that majority of NAHR deletions are non-meiotic i.e. originate from errors during homology directed repair (HDR) of spontaneous DSBs. In turn, the origin of the spontaneous DSBs is associated with transcription factor binding in accessible chromatin revealing the vulnerability of functional, open chromatin. The chromatin itself is enriched with repeats, particularly fixed Alu elements that provide the homology required to maintain stability via HDR. Through co-localization of fixed Alus and NAHR deletions in open chromatin we hypothesize that old Alu expansion had a stabilizing role on the human genome.

The impact of short tandem repeat variation on gene expression

Short tandem repeats (STRs) have been implicated in a variety of complex traits in humans. However, genome-wide studies of the effects of STRs on gene expression thus far have had limited power to detect associations and provide insights into putative mechanisms. Here, we leverage whole-genome sequencing and expression data for 17 tissues from the Genotype-Tissue Expression Project to identify more than 28,000 STRs for which repeat number is associated with expression of nearby genes (eSTRs). We use fine-mapping to quantify the probability that each eSTR is causal and characterize the top 1,400 fine-mapped eSTRs. We identify hundreds of eSTRs linked with published genome-wide association study signals and implicate specific eSTRs in complex traits, including height, schizophrenia, inflammatory bowel disease and intelligence. Overall, our results support the hypothesis that eSTRs contribute to a range of human phenotypes, and our data should serve as a valuable resource for future studies of complex traits.