The human PTGR1 gene expression is controlled by TE-derived Z-DNA forming sequence cooperating with miR-6867-5p.

Z-DNA, a well-known non-canonical form of DNA involved in gene regulation, is often found in gene promoters. Transposable elements (TEs), which make up 45% of the human genome, can move from one location to another within the genome. TEs play various biological roles in host organisms, and like Z-DNA, can influence transcriptional regulation near promoter regions. MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that play a critical role in the regulation of gene expression. Although TEs can generate Z-DNA and miRNAs can bind to Z-DNA, how these factors affect gene transcription has yet to be elucidated. Here, we identified potential Z-DNA forming sequence (ZFS), including TE-derived ZFS, in the promoter of prostaglandin reductase 1 (PTGR1) by data analysis. The transcriptional activity of these ZFS in PTGR1 was confirmed using dual-luciferase reporter assays. In addition, we discovered a novel ZFS-binding miRNA (miR-6867-5p) that suppressed PTGR1 expression by targeting to ZFS. In conclusion, these findings suggest that ZFS, including TE-derived ZFS, can regulate PTGR1 gene expression and that miR-6867-5p can suppress PTGR1 by interacting with ZFS.

Human Endogenous Retrovirus-H-Derived miR-4454 Inhibits the Expression of DNAJB4 and SASH1 in Non-Muscle-Invasive Bladder Cancer.

Although most human endogenous retroviruses (HERVs) have been silenced and lost their ability to translocate because of accumulated mutations during evolution, they still play important roles in human biology. Several studies have demonstrated that HERVs play pathological roles in numerous human diseases, especially cancer. A few studies have revealed that long non-coding RNAs that are transcribed from HERV sequences affect cancer progression. However, there is no study on microRNAs derived from HERVs related to cancer. In this study, we identified 29 microRNAs (miRNAs) derived from HERV sequences in the human genome. In particular, we discovered that miR-4454, which is HERV-H-derived miRNA, was upregulated in non-muscle-invasive bladder cancer (NMIBC) cells. To figure out the effects of upregulated miR-4454 in NMIBC, genes whose expression was downregulated in NMIBC, as well as tumor suppressor genes, were selected as putative target genes of miR-4454. The dual-luciferase assay was used to determine the negative relationship between miR-4454 and its target genes, DNAJB4 and SASH1, and they were confirmed to be promising target genes of miR-4454. Taken together, this study suggests that the upregulation of miR-4454 derived from HERV-H in NMIBC reduces the expression of the tumor suppressor genes, DNAJB4 and SASH1, to promote NMIBC progression.

m1A and m6A modifications function cooperatively to facilitate rapid mRNA degradation.

N6-Methyladenosine (m6A), the most abundant internal mRNA modification, affects multiple steps in gene expression. Mechanistically, the binding of YTHDF2 to m6A on mRNAs elicits rapid mRNA degradation by recruiting several RNA degrading enzymes. Here, we show that N1-methyladenosine (m1A), another type of RNA modification, accelerates rapid m6A RNA degradation. We identify HRSP12 as an RNA-binding protein that recognizes m1A. The binding of HRSP12 to m1A promotes efficient interaction of YTHDF2 with m6A, consequently facilitating endoribonucleolytic cleavage via the RNase P/MRP complex. Transcriptome-wide analyses also reveal that mRNAs harboring both m1A and m6A are downregulated in an HRSP12-dependent manner compared with mRNAs harboring m6A only. Accordingly, a subset of endogenous circular RNAs that harbor m6A and associate with YTHDF2 in an HRSP12-dependent manner is also subjected to m1A-facilitated rapid degradation. Together, our observations provide compelling evidence for crosstalk between different RNA modifications.

Genomic Analyses of Non-Coding RNAs Overlapping Transposable Elements and Its Implication to Human Diseases.

It is estimated that up to 80% of the human genome is transcribed into RNA molecules but less than 2% of the genome encodes the proteins, and the rest of the RNA transcripts that are not translated into protein are called non-coding RNAs (ncRNAs). Many studies have revealed that ncRNAs have biochemical activities as epigenetic regulators at the post-transcriptional level. Growing evidence has demonstrated that transposable elements (TEs) contribute to a large percentage of ncRNAs' transcription. The TEs inserted into certain parts of the genome can act as alternative promoters, enhancers, and insulators, and the accumulation of TEs increases genetic diversity in the human genome. The TEs can also generate microRNAs, so-called miRNA-derived from transposable elements (MDTEs), and are also implicated in disease progression, such as infectious diseases and cancer. Here, we analyzed the origin of ncRNAs and reviewed the published literature on MDTEs related to disease progression.

Z-DNA-Containing Long Terminal Repeats of Human Endogenous Retrovirus Families Provide Alternative Promoters for Human Functional Genes.

Transposable elements (TEs) account for approximately 45% of the human genome. TEs have proliferated randomly and integrated into functional genes during hominoid radiation. They appear as right-handed B-DNA double helices and slightly elongated left-handed Z-DNAs. Human endogenous retrovirus (HERV) families are widely distributed in human chromosomes at a ratio of 8%. They contain a 5'-long terminal repeat (LTR)-gag-pol-env-3'-LTR structure. LTRs contain the U3 enhancer and promoter region, transcribed R region, and U5 region. LTRs can influence host gene expression by acting as regulatory elements. In this review, we describe the alternative promoters derived from LTR elements that overlap Z-DNA by comparing Z-hunt and DeepZ data for human functional genes. We also present evidence showing the regulatory activity of LTR elements containing Z-DNA in GSDML. Taken together, the regulatory activity of LTR elements with Z-DNA allows us to understand gene function in relation to various human diseases.

UPF1 promotes rapid degradation of m6A-containing RNAs.

N6-methyladenosine (m6A) is the most prevalent internal modification in eukaryotic mRNAs and affects RNA processing and metabolism. When YTHDF2, an m6A-recognizing protein, binds to m6A, it facilitates the destabilization of m6A-containing RNAs (m6A RNAs). Here, we demonstrate that upstream frameshift 1 (UPF1), a key factor for nonsense-mediated mRNA decay, interacts with YTHDF2, thereby triggering rapid degradation of m6A RNAs. The UPF1-mediated m6A RNA degradation depends on a specific interaction between UPF1 and N-terminal residues 101-168 of YTHDF2, UPF1 ATPase/helicase activities, and UPF1 interaction with proline-rich nuclear receptor coactivator 2 (PNRC2), a decapping-promoting factor preferentially involved in nonsense-mediated mRNA decay. Furthermore, transcriptome-wide analyses show that YTHDF2-bound mRNAs that are not substrates for HRSP12-RNase P/MRP-mediated endoribonucleolytic cleavage are destabilized with a higher dependency on UPF1. Collectively, our data indicate dynamic and multilayered regulation of the stability of m6A RNAs and highlight the multifaceted role of UPF1 in mRNA decay.

LC3B is an RNA-binding protein to trigger rapid mRNA degradation during autophagy.

LC3/ATG8 has long been appreciated to play a central role in autophagy, by which a variety of cytoplasmic materials are delivered to lysosomes and eventually degraded. However, information on the molecular functions of LC3 in RNA biology is very limited. Here, we show that LC3B is an RNA-binding protein that directly binds to mRNAs with a preference for a consensus AAUAAA motif corresponding to a polyadenylation sequence. Autophagic activation promotes an association between LC3B and target mRNAs and triggers rapid degradation of target mRNAs in a CCR4-NOT-dependent manner before autolysosome formation. Furthermore, our transcriptome-wide analysis reveals that PRMT1 mRNA, which encodes a negative regulator of autophagy, is one of the major substrates. Rapid degradation of PRMT1 mRNA by LC3B facilitates autophagy. Collectively, we demonstrate that LC3B acts as an RNA-binding protein and an mRNA decay factor necessary for efficient autophagy.

The pioneer round of translation ensures proper targeting of ER and mitochondrial proteins.

The pioneer (or first) round of translation of newly synthesized mRNAs is largely mediated by a nuclear cap-binding complex (CBC). In a transcriptome-wide analysis of polysome-associated and CBC-bound transcripts, we identify RN7SL1, a noncoding RNA component of a signal recognition particle (SRP), as an interaction partner of the CBC. The direct CBC-SRP interaction safeguards against abnormal expression of polypeptides from a ribosome-nascent chain complex (RNC)-SRP complex until the latter is properly delivered to the endoplasmic reticulum. Failure of this surveillance causes abnormal expression of misfolded proteins at inappropriate intracellular locations, leading to a cytosolic stress response. This surveillance pathway also blocks protein synthesis through RNC-SRP misassembled on an mRNA encoding a mitochondrial protein. Thus, our results reveal a surveillance pathway in which pioneer translation ensures proper targeting of endoplasmic reticulum and mitochondrial proteins.

Integrated Mobile Element Scanning (ME-Scan) method for identifying multiple types of polymorphic mobile element insertions.

BACKGROUND: Mobile elements are ubiquitous components of mammalian genomes and constitute more than half of the human genome. Polymorphic mobile element insertions (pMEIs) are a major source of human genomic variation and are gaining research interest because of their involvement in gene expression regulation, genome integrity, and disease. RESULTS: Building on our previous Mobile Element Scanning (ME-Scan) protocols, we developed an integrated ME-Scan protocol to identify three major active families of human mobile elements, AluYb, L1HS, and SVA. This approach selectively amplifies insertion sites of currently active retrotransposons for Illumina sequencing. By pooling the libraries together, we can identify pMEIs from all three mobile element families in one sequencing run. To demonstrate the utility of the new ME-Scan protocol, we sequenced 12 human parent-offspring trios. Our results showed high sensitivity (> 90%) and accuracy (> 95%) of the protocol for identifying pMEIs in the human genome. In addition, we also tested the feasibility of identifying somatic insertions using the protocol. CONCLUSIONS: The integrated ME-Scan protocol is a cost-effective way to identify novel pMEIs in the human genome. In addition, by developing the protocol to detect three mobile element families, we demonstrate the flexibility of the ME-Scan protocol. We present instructions for the library design, a sequencing protocol, and a computational pipeline for downstream analyses as a complete framework that will allow researchers to easily adapt the ME-Scan protocol to their own projects in other genomes.

Pedigree-based estimation of human mobile element retrotransposition rates.

Germline mutation rates in humans have been estimated for a variety of mutation types, including single-nucleotide and large structural variants. Here, we directly measure the germline retrotransposition rate for the three active retrotransposon elements: L1, Alu, and SVA. We used three tools for calling mobile element insertions (MEIs) (MELT, RUFUS, and TranSurVeyor) on blood-derived whole-genome sequence (WGS) data from 599 CEPH individuals, comprising 33 three-generation pedigrees. We identified 26 de novo MEIs in 437 births. The retrotransposition rate estimates for Alu elements, one in 40 births, is roughly half the rate estimated using phylogenetic analyses, a difference in magnitude similar to that observed for single-nucleotide variants. The L1 retrotransposition rate is one in 63 births and is within range of previous estimates (1:20-1:200 births). The SVA retrotransposition rate, one in 63 births, is much higher than the previous estimate of one in 900 births. Our large, three-generation pedigrees allowed us to assess parent-of-origin effects and the timing of insertion events in either gametogenesis or early embryonic development. We find a statistically significant paternal bias in Alu retrotransposition. Our study represents the first in-depth analysis of the rate and dynamics of human retrotransposition from WGS data in three-generation human pedigrees.

Staufen1 and UPF1 exert opposite actions on the replacement of the nuclear cap-binding complex by eIF4E at the 5' end of mRNAs.

Newly synthesized mRNAs are exported from the nucleus to cytoplasm with a 5'-cap structure bound by the nuclear cap-binding complex (CBC). During or after export, the CBC should be properly replaced by cytoplasmic cap-binding protein eIF4E for efficient protein synthesis. Nonetheless, little is known about how the replacement takes place. Here, we show that double-stranded RNA-binding protein staufen1 (STAU1) promotes efficient replacement by facilitating an association between the CBC-importin α complex and importin ß. Our transcriptome-wide analyses and artificial tethering experiments also reveal that the replacement occurs more efficiently when an mRNA associates with STAU1. This event is inhibited by a key nonsense-mediated mRNA decay factor, UPF1, which directly interacts with STAU1. Furthermore, we find that cellular apoptosis that is induced by ionizing radiation is accompanied by inhibition of the replacement via increased association between STAU1 and hyperphosphorylated UPF1. Altogether, our data highlight the functional importance of STAU1 and UPF1 in the course of the replacement of the CBC by eIF4E, adding a previously unappreciated layer of post-transcriptional gene regulation.

Endoribonucleolytic Cleavage of m6A-Containing RNAs by RNase P/MRP Complex.

N6-methyladenosine (m6A) is the most abundant internal modification in RNAs and plays regulatory roles in a variety of biological and physiological processes. Despite its important roles, the molecular mechanism underlying m6A-mediated gene regulation is poorly understood. Here, we show that m6A-containing RNAs are subject to endoribonucleolytic cleavage via YTHDF2 (m6A reader protein), HRSP12 (adaptor protein), and RNase P/MRP (endoribonucleases). We demonstrate that HRSP12 functions as an adaptor to bridge YTHDF2 and RNase P/MRP, eliciting rapid degradation of YTHDF2-bound RNAs. Transcriptome-wide analyses show that m6A RNAs that are preferentially targeted for endoribonucleolytic cleavage have an HRSP12-binding site and a RNase P/MRP-directed cleavage site upstream and downstream of the YTHDF2-binding site, respectively. We also find that a subset of m6A-containing circular RNAs associates with YTHDF2 in an HRSP12-dependent manner and is selectively downregulated by RNase P/MRP. Thus, our data expand the known functions of RNase P/MRP to endoribonucleolytic cleavage of m6A RNAs.

eIF4A3 Phosphorylation by CDKs Affects NMD during the Cell Cycle.

Exon junction complexes (EJCs) loaded onto spliced mRNAs during splicing serve as molecular markers for various post-transcriptional gene-regulatory processes, including nonsense-mediated mRNA decay (NMD). Although the composition and structure of EJCs are well characterized, the mechanism regulating EJC deposition remains unknown. Here we find that threonine 163 (T163) within the RNA-binding motif of eIF4A3 (a core EJC component) is phosphorylated by cyclin-dependent protein kinases 1 and 2 in a cell cycle-dependent manner. T163 phosphorylation hinders binding of eIF4A3 to spliced mRNAs and other EJC components. Instead, it promotes association of eIF4A3 with CWC22, which guides eIF4A3 to an active spliceosome. These molecular events ensure the fidelity of specific deposition of the EJC ∼20-24 nt upstream of an exon-exon junction. Accordingly, NMD is affected by T163 phosphorylation. Collectively, our data provide evidence that T163 phosphorylation affects EJC formation and, consequently, NMD efficiency in a cell cycle-dependent manner.

DLK regulates a distinctive transcriptional regeneration program after peripheral nerve injury.

Following damage to a peripheral nerve, injury signaling pathways converge in the cell body to generate transcriptional changes that support axon regeneration. Here, we demonstrate that dual leucine zipper kinase (DLK), a central regulator of injury responses including axon regeneration and neuronal apoptosis, is required for the induction of the pro-regenerative transcriptional program in response to peripheral nerve injury. Using a sensory neuron-conditional DLK knockout mouse model, we show a time course for the dependency of gene expression changes on the DLK pathway after sciatic nerve injury. Gene ontology analysis reveals that DLK-dependent gene sets are enriched for specific functional annotations such as ion transport and immune response. A series of comparative analyses shows that the DLK-dependent transcriptional program is distinct from that promoted by the importin-dependent retrograde signaling pathway, while it is partially shared between PNS and CNS injury responses. We suggest that DLK-dependency might provide a selective filter for regeneration-associated genes among the injury-responsive transcriptome.

HuR stabilizes a polyadenylated form of replication-dependent histone mRNAs under stress conditions.

All metazoan mRNAs have a poly(A) tail at the 3' end with the exception of replication-dependent histone (RDH) mRNAs, which end in a highly conserved stem-loop (SL) structure. However, a subset of RDH mRNAs are reported to be polyadenylated under physiologic conditions. The molecular details of the biogenesis of polyadenylated RDH [poly(A)+ RDH] mRNAs remain unknown. In this study, our genome-wide analyses reveal that puromycin treatment or UVC irradiation stabilizes poly(A)+ RDH mRNAs, relative to canonical RDH mRNAs, which end in an SL structure. We demonstrate that the stabilization of poly(A)+ RDH mRNAs occurs in a translation-independent manner and is regulated via human antigen R (HuR) binding to the extended 3' UTR under stress conditions. Our data suggest that HuR regulates the expression of poly(A)+ RDH mRNAs.-Ryu, I., Park, Y., Seo, J.-W., Park, O. H., Ha, H., Nam, J.-W., Kim, Y. K. HuR stabilizes a polyadenylated form of replication-dependent histone mRNAs under stress conditions.

Comparative analysis of the transcriptome of injured nerve segments reveals spatiotemporal responses to neural damage in mice.

Nerves are particularly vulnerable to damage due to their unique structure with meter-long axons. In the peripheral nervous system neurons and Schwann cells can activate the injury-response program that directs axons to either regenerate or degenerate after traumatic nerve injury. However, the differences between the genetic programs driving nerve regeneration and degeneration have not yet been described extensively. To understand these differences, in this study we have compared the injury-induced transcriptomic changes between the regenerating proximal segment and the degenerating distal segment of a transected nerve, at different post-injury time points. We analyzed the spatiotemporal dynamics of the mouse transcriptome using a sciatic nerve-injury model by means of RNA sequencing. The results of the differentially regulated genes (DEGs) analysis show that some DEG groups are similarly regulated in both proximal and distal segments, and primarily display a positive correlation. However, some DEG groups are exclusively regulated in either the proximal or the distal segment, suggesting that these DEG groups constitute a genetic network for distinguishing the regenerative and degenerative responses. In addition, our gene ontology analysis revealed an enrichment of particular biological processes in different phases and locations. Thus, our data provide a spatiotemporal profile of the transcriptomes that are differentially regulated in either regenerating or degenerating nerves, in vivo. The specific biological processes enriched in the DEG groups might delineate the injury-responsive program that induces contrasting regenerative and degenerative responses in different nerve segments.

Long non-coding RNAs transcribed by ERV-9 LTR retrotransposon act in cis to modulate long-range LTR enhancer function.

LTR retrotransposons are repetitive DNA elements comprising ∼10% of the human genome. However, LTR sequences are disproportionately present in human long, non-coding RNAs (lncRNAs). Whether and how the LTR lncRNAs serve biological functions are largely unknown. Here we show that in primary human erythroblasts, lncRNAs transcribed from the LTR retrotransposons of ERV-9 human endogenous retrovirus activated transcription of key erythroid genes and modulated ex vivo erythropoiesis. To dissect the functional mechanism of ERV-9 lncRNAs, we performed genome-wide RNA and ChIRP analyses before and after global knockdown or locus-specific deletion of ERV-9 lncRNAs in human erythroblasts carrying ∼4000 copies of the ERV-9 LTRs and in transgenic mouse erythroblasts carrying a single copy of the primate-specific ERV-9 LTR in the 100 kb human ß-globin gene locus. We found that ERV-9 lncRNAs acted in cis to stabilize assembly of the ERV-9 LTR enhancer complex and facilitate long-range LTR enhancer function in activating transcription of downstream, cis-linked globin genes. Our findings suggested that LTR lncRNAs transcribed from many of the 4000 copies of ERV-9 LTR retrotransposons acted by a similar cis mechanism to modulate LTR enhancer function in activating transcription of downstream genes critical to cellular processes including erythropoiesis.

Library Construction for High-Throughput Mobile Element Identification and Genotyping.

Mobile genetic elements are discrete DNA elements that can move around and copy themselves in a genome. As a ubiquitous component of the genome, mobile elements contribute to both genetic and epigenetic variation. Therefore, it is important to determine the genome-wide distribution of mobile elements. Here we present a targeted high-throughput sequencing protocol called Mobile Element Scanning (ME-Scan) for genome-wide mobile element detection. We will describe oligonucleotides design, sequencing library construction, and computational analysis for the ME-Scan protocol.

Identification of polymorphic SVA retrotransposons using a mobile element scanning method for SVA (ME-Scan-SVA).

BACKGROUND: Mobile element insertions are a major source of human genomic variation. SVA (SINE-R/VNTR/Alu) is the youngest retrotransposon family in the human genome and a number of diseases are known to be caused by SVA insertions. However, inter-individual genomic variations generated by SVA insertions and their impacts have not been studied extensively due to the difficulty in identifying polymorphic SVA insertions. RESULTS: To systematically identify SVA insertions at the population level and assess their genomic impact, we developed a mobile element scanning (ME-Scan) protocol we called ME-Scan-SVA. Using a nested SVA-specific PCR enrichment method, ME-Scan-SVA selectively amplify the 5' end of SVA elements and their flanking genomic regions. To demonstrate the utility of the protocol, we constructed and sequenced a ME-Scan-SVA library of 21 individuals and analyzed the data using a new analysis pipeline designed for the protocol. Overall, the method achieved high SVA-specificity and over >90 % of the sequenced reads are from SVA insertions. The method also had high sensitivity (>90 %) for fixed SVA insertions that contain the SVA-specific primer-binding sites in the reference genome. Using candidate locus selection criteria that are expected to have a 90 % sensitivity, we identified 151 and 29 novel polymorphic SVA candidates under relaxed and stringent cutoffs, respectively (average 12 and 2 per individual). For six polymorphic SVAs that we were able to validate by PCR, the average individual genotype accuracy is 92 %, demonstrating a high accuracy of the computational genotype calling pipeline. CONCLUSIONS: The new approach allows identifying novel SVA insertions using high-throughput sequencing. It is cost-effective and can be applied in large-scale population study. It also can be applied for detecting potential active SVA elements, and somatic SVA retrotransposition events in different tissues or developmental stages.