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.

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.

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.

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.

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.

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.

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.

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.

piRNAs derived from ancient viral processed pseudogenes as transgenerational sequence-specific immune memory in mammals.

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences within vertebrate genomes derived from reverse transcription and integration of ancient bornaviral nucleoprotein mRNA via the host retrotransposon machinery. While species with EBLNs appear relatively resistant to bornaviral disease, the nature of this association is unclear. We hypothesized that EBLNs could give rise to antiviral interfering RNA in the form of PIWI-interacting RNAs (piRNAs), a class of small RNA known to silence transposons but not exogenous viruses. We found that in both rodents and primates, which acquired their EBLNs independently some 25-40 million years ago, EBLNs are present within piRNA-generating regions of the genome far more often than expected by chance alone (ℙ = 8 × 10(-3)-6 × 10(-8)). Three of the seven human EBLNs fall within annotated piRNA clusters and two marmoset EBLNs give rise to bona fide piRNAs. In both rats and mice, at least two of the five EBLNs give rise to abundant piRNAs in the male gonad. While no EBLNs are syntenic between rodent and primate, some of the piRNA clusters containing EBLNs are; thus we deduce that EBLNs were integrated into existing piRNA clusters. All true piRNAs derived from EBLNs are antisense relative to the proposed ancient bornaviral nucleoprotein mRNA. These observations are consistent with a role for EBLN-derived piRNA-like RNAs in interfering with ancient bornaviral infection. They raise the hypothesis that retrotransposon-dependent virus-to-host gene flow could engender RNA-mediated, sequence-specific antiviral immune memory in metazoans analogous to the CRISPR/Cas system in prokaryotes.

Whole exome sequencing identifies novel candidate genes that modify chronic obstructive pulmonary disease susceptibility.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is characterized by an irreversible airflow limitation in response to inhalation of noxious stimuli, such as cigarette smoke. However, only 15-20 % smokers manifest COPD, suggesting a role for genetic predisposition. Although genome-wide association studies have identified common genetic variants that are associated with susceptibility to COPD, effect sizes of the identified variants are modest, as is the total heritability accounted for by these variants. In this study, an extreme phenotype exome sequencing study was combined with in vitro modeling to identify COPD candidate genes. RESULTS: We performed whole exome sequencing of 62 highly susceptible smokers and 30 exceptionally resistant smokers to identify rare variants that may contribute to disease risk or resistance to COPD. This was a cross-sectional case-control study without therapeutic intervention or longitudinal follow-up information. We identified candidate genes based on rare variant analyses and evaluated exonic variants to pinpoint individual genes whose function was computationally established to be significantly different between susceptible and resistant smokers. Top scoring candidate genes from these analyses were further filtered by requiring that each gene be expressed in human bronchial epithelial cells (HBECs). A total of 81 candidate genes were thus selected for in vitro functional testing in cigarette smoke extract (CSE)-exposed HBECs. Using small interfering RNA (siRNA)-mediated gene silencing experiments, we showed that silencing of several candidate genes augmented CSE-induced cytotoxicity in vitro. CONCLUSIONS: Our integrative analysis through both genetic and functional approaches identified two candidate genes (TACC2 and MYO1E) that augment cigarette smoke (CS)-induced cytotoxicity and, potentially, COPD susceptibility.

Mobile element scanning (ME-Scan) identifies thousands of novel Alu insertions in diverse human populations.

Alu retrotransposons are the most numerous and active mobile elements in humans, causing genetic disease and creating genomic diversity. Mobile element scanning (ME-Scan) enables comprehensive and affordable identification of mobile element insertions (MEI) using targeted high-throughput sequencing of multiplexed MEI junction libraries. In a single experiment, ME-Scan identifies nearly all AluYb8 and AluYb9 elements, with high sensitivity for both rare and common insertions, in 169 individuals of diverse ancestry. ME-Scan detects heterozygous insertions in single individuals with 91% sensitivity. Insertion presence or absence states determined by ME-Scan are 95% concordant with those determined by locus-specific PCR assays. By sampling diverse populations from Africa, South Asia, and Europe, we are able to identify 5799 Alu insertions, including 2524 novel ones, some of which occur in exons. Sub-Saharan populations and a Pygmy group in particular carry numerous intermediate-frequency Alu insertions that are absent in non-African groups. There is a significant dearth of exon-interrupting insertions among common Alu polymorphisms, but the density of singleton Alu insertions is constant across exonic and nonexonic regions. In one case, a validated novel singleton Alu interrupts a protein-coding exon of FAM187B. This implies that exonic Alu insertions are generally deleterious and thus eliminated by natural selection, but not so quickly that they cannot be observed as extremely rare variants.

A comprehensive analysis of piRNAs from adult human testis and their relationship with genes and mobile elements.

BACKGROUND: Piwi-interacting RNAs (piRNAs) are a recently discovered class of small non-coding RNAs whose best-understood function is to repress mobile element (ME) activity in animal germline. To date, nearly all piRNA studies have been conducted in model organisms and little is known about piRNA diversity, target specificity and biological function in human. RESULTS: Here we performed high-throughput sequencing of piRNAs from three human adult testis samples. We found that more than 81% of the ~17 million putative piRNAs mapped to ~6,000 piRNA-producing genomic clusters using a relaxed definition of clusters. A set of human protein-coding genes produces a relatively large amount of putative piRNAs from their 3'UTRs, and are significantly enriched for certain biological processes, suggestive of non-random sampling by the piRNA biogenesis machinery. Up to 16% of putative piRNAs mapped to a few hundred annotated long non-coding RNA (lncRNA) genes, suggesting that some lncRNA genes can act as piRNA precursors. Among major ME families, young families of LTR and endogenous retroviruses have a greater association with putative piRNAs than other MEs. In addition, piRNAs preferentially mapped to specific regions in the consensus sequences of several ME (sub)families and some piRNA mapping peaks showed patterns consistent with the "ping-pong" cycle of piRNA targeting and amplification. CONCLUSIONS: Overall our data provide a comprehensive analysis and improved annotation of human piRNAs in adult human testes and shed new light into the relationship of piRNAs with protein-coding genes, lncRNAs, and mobile genetic elements in human.

Genome-wide analysis of DNA methylation patterns in horse.

BACKGROUND: DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known. RESULTS: We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density. CONCLUSIONS: We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.

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.

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.

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.