We evaluated interactions of SNP-by-ACE-I/ARB and SNP-by-TD on serum potassium (K+) among users of antihypertensive treatments (anti-HTN). Our study included seven European-ancestry (EA) (N = 4835) and four African-ancestry (AA) cohorts (N = 2016). We performed race-stratified, fixed-effect, inverse-variance-weighted meta-analyses of 2.5 million SNP-by-drug interaction estimates; race-combined meta-analysis; and trans-ethnic fine-mapping. Among EAs, we identified 11 significant SNPs (P < 5 × 10-8) for SNP-ACE-I/ARB interactions on serum K+ that were located between NR2F1-AS1 and ARRDC3-AS1 on chromosome 5 (top SNP rs6878413 P = 1.7 × 10-8; ratio of serum K+ in ACE-I/ARB exposed compared to unexposed is 1.0476, 1.0280, 1.0088 for the TT, AT, and AA genotypes, respectively). Trans-ethnic fine mapping identified the same group of SNPs on chromosome 5 as genome-wide significant for the ACE-I/ARB analysis. In conclusion, SNP-by-ACE-I /ARB interaction analyses uncovered loci that, if replicated, could have future implications for the prevention of arrhythmias due to anti-HTN treatment-related hyperkalemia. Before these loci can be identified as clinically relevant, future validation studies of equal or greater size in comparison to our discovery effort are needed.

Data from the Encyclopedia of DNA Elements (ENCODE) project show over 9640 human genome loci classified as long noncoding RNAs (lncRNAs), yet only ~100 have been deeply characterized to determine their role in the cell. To measure the protein-coding output from these RNAs, we jointly analyzed two recent data sets produced in the ENCODE project: tandem mass spectrometry (MS/MS) data mapping expressed peptides to their encoding genomic loci, and RNA-seq data generated by ENCODE in long polyA+ and polyA- fractions in the cell lines K562 and GM12878. We used the machine-learning algorithm RuleFit3 to regress the peptide data against RNA expression data. The most important covariate for predicting translation was, surprisingly, the Cytosol polyA- fraction in both cell lines. LncRNAs are ~13-fold less likely to produce detectable peptides than similar mRNAs, indicating that ~92% of GENCODE v7 lncRNAs are not translated in these two ENCODE cell lines. Intersecting 9640 lncRNA loci with 79,333 peptides yielded 85 unique peptides matching 69 lncRNAs. Most cases were due to a coding transcript misannotated as lncRNA. Two exceptions were an unprocessed pseudogene and a bona fide lncRNA gene, both with open reading frames (ORFs) compromised by upstream stop codons. All potentially translatable lncRNA ORFs had only a single peptide match, indicating low protein abundance and/or false-positive peptide matches. We conclude that with very few exceptions, ribosomes are able to distinguish coding from noncoding transcripts and, hence, that ectopic translation and cryptic mRNAs are rare in the human lncRNAome.

Previous research suggested that a novel compound PNT2258 inhibits B-cell lymphoma 2 (BCL-2) transcription by DNA interference (DNAi) and demonstrated its activity in preclinical xenograft models and in a pilot Phase II clinical trial in non-Hodgkin's lymphoma (NHL). While the drug downregulates BCL-2 at the promoter, mRNA, and protein levels, there is a significant homology (13-16 bases) between PNT100 and a number of promoters of genes involved in cell cycle regulation and survival. In this study, we identify cyclin-dependent kinase-4 (CDK4) as an unintended target gene of PNT2258 and examine its relevance to NHL.

Long non-coding RNAs (lncRNAs) are transcripts of a recently discovered class of genes which do not code for proteins. LncRNA genes are approximately as numerous as protein-coding genes in the human genome. However, comparatively little remains known about lncRNA functions. We globally interrogated changes in the lncRNA transcriptome of oestrogen receptor positive human breast cancer cells following treatment with oestrogen, and identified 127 oestrogen-responsive lncRNAs. Consistent with the emerging evidence that most human lncRNA genes lack homologues outside of primates, our evolutionary analysis revealed primate-specific lncRNAs downstream of oestrogen signalling. We demonstrate, using multiple functional assays to probe gain- and loss-of-function phenotypes in two oestrogen receptor positive human breast cancer cell lines, that two primate-specific oestrogen-responsive lncRNAs identified in this study (the oestrogen-repressed lncRNA BC041455, which reduces cell viability, and the oestrogen-induced lncRNA CR593775, which increases cell viability) exert previously unrecognized functions in cell proliferation and growth factor signalling pathways. The results suggest that oestrogen-responsive lncRNAs are capable of altering the proliferation and viability of human breast cancer cells. No effects on cellular phenotypes were associated with control transfections. As heretofore unappreciated components of key signalling pathways in cancers, including the MAP kinase pathway, lncRNAs hence represent a novel mechanism of action for oestrogen effects on cellular proliferation and viability phenotypes. This finding warrants further investigation in basic and translational studies of breast and potentially other types of cancers, has broad relevance to lncRNAs in other nuclear hormone receptor pathways, and should facilitate exploiting and targeting these cell viability modulating lncRNAs in post-genomic therapeutics.

In the post-genomic era, our understanding of the molecular regulators of physiologic and pathologic processes in pregnancy is expanding at the whole-genome level. Longitudinal changes in the known protein-coding transcriptome during normal pregnancy, which we recently reported (Gomez-Lopez et al., 2019), have improved our definition of the major operant networks, yet pregnancy-related functions of the non-coding RNA transcriptome remain poorly understood. A key finding of the ENCODE (Encyclopedia of DNA Elements) Consortium, the successor of the Human Genome Project, was that the human genome contains approximately 60,000 genes, the majority of which do not encode proteins. The total transcriptional output of non-protein-coding RNA genes, collectively referred to as the non-coding transcriptome, is comprised mainly of long non-coding RNA (lncRNA) transcripts (Derrien et al., 2012). Although the ncRNA transcriptome eclipses its protein-coding counterpart in abundance, it has until recently lacked a comprehensive, unbiased, genome-scale characterization over the timecourse of normal human pregnancy. Here, we annotated, characterized, and selectively validated the longitudinal changes in the non-coding transcriptome of maternal whole blood during normal pregnancy to term. We identified nine long non-coding RNAs (lncRNAs), including long intergenic non-coding RNAs (lincRNAs) as well as lncRNAs antisense to or otherwise in the immediate vicinity of protein-coding genes, that were differentially expressed with advancing gestation in normal pregnancy: AL355711, BC039551 (expressed mainly in the placenta), JHDM1D-AS1, A2M-AS1, MANEA-AS1, NR_034004, LINC00649, LINC00861, and LINC01094. By cross-referencing our dataset against major public pseudogene catalogs, we also identified six transcribed pseudogenes that were differentially expressed over time during normal pregnancy in maternal blood: UBBP4, FOXO3B, two Makorin (MKRN) pseudogenes (MKRN9P and LOC441455), PSME2P2, and YBX3P1. We also identified three non-coding RNAs belonging to other classes that were modulated during gestation: the microRNA MIR4439, the small nucleolar RNA (snoRNA) SNORD41, and the small Cajal-body specific ncRNA SCARNA2. The expression profiles of most hits were broadly suggestive of functions in pregnancy. These time-dependent changes of the non-coding transcriptome during normal pregnancy, which may confer specific regulatory impacts on their protein-coding gene targets, will facilitate a deeper molecular understanding of pregnancy and lncRNA-mediated molecular pathways at the maternal-fetal interface and of how these pathways impact maternal and fetal health.

The human brain has evolved to have extraordinary capabilities, enabling complex behaviors. The uniqueness of the human brain is increasingly posited to be due in part to the functions of primate-specific, including human-specific, long non-coding RNA (lncRNA) genes, systemically less conserved than protein-coding genes in evolution. Patients who have surgery for drug-resistant epilepsy are subjected to extensive electrical recordings of the brain tissue that is subsequently removed in order to treat their epilepsy. Precise localization of brain tissues with distinct electrical properties offers a rare opportunity to explore the effects of brain activity on gene expression. Here, we identified 231 co-regulated, activity-dependent lncRNAs within the human MAPK signaling cascade. Six lncRNAs, four of which were antisense to known protein-coding genes, were further examined because of their high expression and potential impact on the disease phenotype. Using a model of repeated depolarizations in human neuronal-like cells (Sh-SY5Y), we show that five out of six lncRNAs were electrical activity-dependent, with three of four antisense lncRNAs having reciprocal expression patterns relative to their protein-coding gene partners. Some were directly regulated by MAPK signaling, while others effectively downregulated the expression of the protein-coding genes encoded on the opposite strands of their genomic loci. These lncRNAs, therefore, likely contribute to highly evolved and primate-specific human brain regulatory functions that could be therapeutically modulated to treat epilepsy.

Yin Yang 1 (YY1) regulates gene transcription in a variety of biological processes. In this study, we aim to determine the role of YY1 in vascular smooth muscle cell (VSMC) phenotypic modulation both in vivo and in vitro. Here we show that vascular injury in rodent carotid arteries induces YY1 expression along with reduced expression of smooth muscle differentiation markers in the carotids. Consistent with this finding, YY1 expression is induced in differentiated VSMCs in response to serum stimulation. To determine the underlying molecular mechanisms, we found that YY1 suppresses the transcription of CArG box-dependent SMC-specific genes including SM22α, SMα-actin and SMMHC. Interestingly, YY1 suppresses the transcriptional activity of the SM22α promoter by hindering the binding of serum response factor (SRF) to the proximal CArG box. YY1 also suppresses the transcription and the transactivation of myocardin (MYOCD), a master regulator for SMC-specific gene transcription by binding to SRF to form the MYOCD/SRF/CArG box triad (known as the ternary complex). Mechanistically, YY1 directly interacts with MYOCD to competitively displace MYOCD from SRF. This is the first evidence showing that YY1 inhibits SMC differentiation by directly targeting MYOCD. These findings provide new mechanistic insights into the regulatory mechanisms that govern SMC phenotypic modulation in the pathogenesis of vascular diseases.

Long non-coding RNAs (lncRNAs) play important roles in the maintenance of metabolic homeostasis. Recently, many studies have suggested that lncRNAs, such as Metastasis Associated Lung Adenocarcinoma Transcript 1 (MALAT1) and Imprinted Maternally Expressed Transcript (H19), might participate in the pathogenesis of metabolic disorders such as obesity. We conducted a case-control study with 150 Russian children and adolescents aged between 5 and 17 years old in order to assess the statistical association between the single nucleotide polymorphisms (SNPs) rs3200401 in MALAT1 and rs217727 in H19, and the risk of developing obesity in this population. We further explored the possible association of rs3200401 and rs217727 with BMI Z-score and insulin resistance. The MALAT1 rs3200401 and H19 rs217727 SNPs were genotyped using Taqman SNP genotyping assay. The MALAT1 rs3200401 SNP was identified as a risk factor for childhood obesity (p < 0.05) under the dominant and allelic models, and the CT heterozygous genotype was associated with the risk of increased BMI and with insulin resistance. The H19 rs217727 SNP had no significant association with obesity risk (all p > 0.05). Our findings thus suggest that MALAT1 SNP rs3200401 is a potential indicator of obesity susceptibility and pathogenesis in children and adolescents.

Autoantibodies against mitochondrial-derived antigens play a key role in chronic tissue inflammation in autoimmune disorders and cancers. Here, we identify autoreactive nuclear genomic DNA (nDNA)-encoded mitochondrial gene products (GAPDH, PKM2, GSTP1, SPATA5, MFF, TSPOAP1, PHB2, COA4, and HAGH) recognized by breast cancer (BC) patients' sera as nonself, supporting a direct relationship of mitochondrial autoimmunity to breast carcinogenesis. Autoreactivity of multiple nDNA-encoded mitochondrial gene products was mapped to protein-coding regions, 3' untranslated regions (UTRs), as well as introns. In addition, autoantibodies in BC sera targeted intergenic sequences that may be parts of long non-coding RNA (lncRNA) genes, including LINC02381 and other putative lncRNA neighbors of the protein-coding genes ERCC4, CXCL13, SOX3, PCDH1, EDDM3B, and GRB2. Increasing evidence indicates that lncRNAs play a key role in carcinogenesis. Consistent with this, our findings suggest that lncRNAs, as well as mRNAs of nDNA-encoded mitochondrial genes, mechanistically contribute to BC progression. This work supports a new paradigm of breast carcinogenesis based on a globally dysfunctional genome with altered function of multiple mitochondrial and non-mitochondrial oncogenic pathways caused by the effects of autoreactivity-induced dysregulation of multiple genes and their products. This autoimmunity-based model of carcinogenesis will open novel avenues for BC treatment.

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol l(-1), P=3.4 × 10(-12)), T2D risk (OR[95%CI]=0.86[0.76-0.96], P=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose(-1), P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l(-1), P=4.3 × 10(-4)). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8 × 10(-6)) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l(-1), P=1.3 × 10(-8)). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

The maintenance of pluripotency and specification of cellular lineages during embryonic development are controlled by transcriptional regulatory networks, which coordinate specific sets of genes through both activation and repression. The transcriptional repressor RE1-silencing transcription factor (REST) plays important but distinct regulatory roles in embryonic (ESC) and neural (NSC) stem cells. We investigated how these distinct biological roles are effected at a genomic level. We present integrated, comparative genome- and transcriptome-wide analyses of transcriptional networks governed by REST in mouse ESC and NSC. The REST recruitment profile has dual components: a developmentally independent core that is common to ESC, NSC, and differentiated cells; and a large, ESC-specific set of target genes. In ESC, the REST regulatory network is highly integrated into that of pluripotency factors Oct4-Sox2-Nanog. We propose that an extensive, pluripotency-specific recruitment profile lends REST a key role in the maintenance of the ESC phenotype.

The human genome contains many thousands of long noncoding RNAs (lncRNAs). While several studies have demonstrated compelling biological and disease roles for individual examples, analytical and experimental approaches to investigate these genes have been hampered by the lack of comprehensive lncRNA annotation. Here, we present and analyze the most complete human lncRNA annotation to date, produced by the GENCODE consortium within the framework of the ENCODE project and comprising 9277 manually annotated genes producing 14,880 transcripts. Our analyses indicate that lncRNAs are generated through pathways similar to that of protein-coding genes, with similar histone-modification profiles, splicing signals, and exon/intron lengths. In contrast to protein-coding genes, however, lncRNAs display a striking bias toward two-exon transcripts, they are predominantly localized in the chromatin and nucleus, and a fraction appear to be preferentially processed into small RNAs. They are under stronger selective pressure than neutrally evolving sequences-particularly in their promoter regions, which display levels of selection comparable to protein-coding genes. Importantly, about one-third seem to have arisen within the primate lineage. Comprehensive analysis of their expression in multiple human organs and brain regions shows that lncRNAs are generally lower expressed than protein-coding genes, and display more tissue-specific expression patterns, with a large fraction of tissue-specific lncRNAs expressed in the brain. Expression correlation analysis indicates that lncRNAs show particularly striking positive correlation with the expression of antisense coding genes. This GENCODE annotation represents a valuable resource for future studies of lncRNAs.

Using a chromatin immunoprecipitation-paired end diTag cloning and sequencing strategy, we mapped estrogen receptor alpha (ERalpha) binding sites in MCF-7 breast cancer cells. We identified 1,234 high confidence binding clusters of which 94% are projected to be bona fide ERalpha binding regions. Only 5% of the mapped estrogen receptor binding sites are located within 5 kb upstream of the transcriptional start sites of adjacent genes, regions containing the proximal promoters, whereas vast majority of the sites are mapped to intronic or distal locations (>5 kb from 5' and 3' ends of adjacent transcript), suggesting transcriptional regulatory mechanisms over significant physical distances. Of all the identified sites, 71% harbored putative full estrogen response elements (EREs), 25% bore ERE half sites, and only 4% had no recognizable ERE sequences. Genes in the vicinity of ERalpha binding sites were enriched for regulation by estradiol in MCF-7 cells, and their expression profiles in patient samples segregate ERalpha-positive from ERalpha-negative breast tumors. The expression dynamics of the genes adjacent to ERalpha binding sites suggest a direct induction of gene expression through binding to ERE-like sequences, whereas transcriptional repression by ERalpha appears to be through indirect mechanisms. Our analysis also indicates a number of candidate transcription factor binding sites adjacent to occupied EREs at frequencies much greater than by chance, including the previously reported FOXA1 sites, and demonstrate the potential involvement of one such putative adjacent factor, Sp1, in the global regulation of ERalpha target genes. Unexpectedly, we found that only 22%-24% of the bona fide human ERalpha binding sites were overlapping conserved regions in whole genome vertebrate alignments, which suggest limited conservation of functional binding sites. Taken together, this genome-scale analysis suggests complex but definable rules governing ERalpha binding and gene regulation.

Non-coding RNAs (ncRNAs) are essential players in many cellular processes, from normal development to oncogenic transformation. Initially, ncRNAs were defined as transcripts that lacked an open reading frame (ORF). However, multiple lines of evidence suggest that certain ncRNAs encode small peptides of less than 100 amino acids. The sequences encoding these peptides are known as small open reading frames (smORFs), many initiating with the traditional AUG start codon but terminating with atypical stop codons, suggesting a different biogenesis. The ncRNA-encoded peptides (ncPEPs) are gradually becoming appreciated as a new class of functional molecules that contribute to diverse cellular processes, and are deregulated in different diseases contributing to pathogenesis. As multiple publications have identified unique ncPEPs, we appreciated the need for assembling a new web resource that could gather information about these functional ncPEPs. We developed FuncPEP, a new database of functional ncRNA encoded peptides, containing all experimentally validated and functionally characterized ncPEPs. Currently, FuncPEP includes a comprehensive annotation of 112 functional ncPEPs and specific details regarding the ncRNA transcripts that encode these peptides. We believe that FuncPEP will serve as a platform for further deciphering the biologic significance and medical use of ncPEPs.

Transcriptional regulation of protein-coding genes is increasingly well-understood on a global scale, yet no comparable information exists for long non-coding RNA (lncRNA) genes, which were recently recognized to be as numerous as protein-coding genes in mammalian genomes. We performed a genome-wide comparative analysis of the promoters of human lncRNA and protein-coding genes, finding global differences in specific genetic and epigenetic features relevant to transcriptional regulation. These two groups of genes are hence subject to separate transcriptional regulatory programs, including distinct transcription factor (TF) proteins that significantly favor lncRNA, rather than coding-gene, promoters. We report a specific signature of promoter-proximal transcriptional regulation of lncRNA genes, including several distinct transcription factor binding sites (TFBS). Experimental DNase I hypersensitive site profiles are consistent with active configurations of these lncRNA TFBS sets in diverse human cell types. TFBS ChIP-seq datasets confirm the binding events that we predicted using computational approaches for a subset of factors. For several TFs known to be directly regulated by lncRNAs, we find that their putative TFBSs are enriched at lncRNA promoters, suggesting that the TFs and the lncRNAs may participate in a bidirectional feedback loop regulatory network. Accordingly, cells may be able to modulate lncRNA expression levels independently of mRNA levels via distinct regulatory pathways. Our results also raise the possibility that, given the historical reliance on protein-coding gene catalogs to define the chromatin states of active promoters, a revision of these chromatin signature profiles to incorporate expressed lncRNA genes is warranted in the future.

Despite their abundance, the molecular functions of long non-coding RNAs in mammalian nervous systems remain poorly understood. Here we show that the long non-coding RNA, NEAT1, directly modulates neuronal excitability and is associated with pathological seizure states. Specifically, NEAT1 is dynamically regulated by neuronal activity in vitro and in vivo, binds epilepsy-associated potassium channel-interacting proteins including KCNAB2 and KCNIP1, and induces a neuronal hyper-potentiation phenotype in iPSC-derived human cortical neurons following antisense oligonucleotide knockdown. Next generation sequencing reveals a strong association of NEAT1 with increased ion channel gene expression upon activation of iPSC-derived neurons following NEAT1 knockdown. Furthermore, we show that while NEAT1 is acutely down-regulated in response to neuronal activity, repeated stimulation results in NEAT1 becoming chronically unresponsive in independent in vivo rat model systems relevant to temporal lobe epilepsy. We extended previous studies showing increased NEAT1 expression in resected cortical tissue from high spiking regions of patients suffering from intractable seizures. Our results indicate a role for NEAT1 in modulating human neuronal activity and suggest a novel mechanistic link between an activity-dependent long non-coding RNA and epilepsy.

The human brain and human cognitive abilities are strikingly different from those of other great apes despite relatively modest genome sequence divergence. However, little is presently known about the interspecies divergence in gene structure and transcription that might contribute to these phenotypic differences. To date, most comparative studies of gene structure in the brain have examined humans, chimpanzees, and macaque monkeys. To add to this body of knowledge, we analyze here the brain transcriptome of the western lowland gorilla (Gorilla gorilla gorilla), an African great ape species that is phylogenetically closely related to humans, but with a brain that is approximately one-third the size. Manual transcriptome curation from a sample of the planum temporale region of the neocortex revealed 12 protein-coding genes and one noncoding-RNA gene with exons in the gorilla unmatched by public transcriptome data from the orthologous human loci. These interspecies gene structure differences accounted for a total of 134 amino acids in proteins found in the gorilla that were absent from protein products of the orthologous human genes. Proteins varying in structure between human and gorilla were involved in immunity and energy metabolism, suggesting their relevance to phenotypic differences. This gorilla neocortical transcriptome comprises an empirical, not homology- or prediction-driven, resource for orthologous gene comparisons between human and gorilla. These findings provide a unique repository of the sequences and structures of thousands of genes transcribed in the gorilla brain, pointing to candidate genes that may contribute to the traits distinguishing humans from other closely related great apes.

Pangolins are scale-covered mammals, containing eight endangered species. Maintaining pangolins in captivity is a significant challenge, in part because little is known about their genetics. Here we provide the first large-scale sequencing of the critically endangered Manis javanica transcriptomes from eight different organs using Illumina HiSeq technology, yielding ~75 Giga bases and 89,754 unigenes. We found some unigenes involved in the insect hormone biosynthesis pathway and also 747 lipids metabolism-related unigenes that may be insightful to understand the lipid metabolism system in pangolins. Comparative analysis between M. javanica and other mammals revealed many pangolin-specific genes significantly over-represented in stress-related processes, cell proliferation and external stimulus, probably reflecting the traits and adaptations of the analyzed pregnant female M. javanica. Our study provides an invaluable resource for future functional works that may be highly relevant for the conservation of pangolins.

Abdominal aortic aneurysm (AAA) is a complex disease with both genetic and environmental risk factors. Together, 6 previously identified risk loci only explain a small proportion of the heritability of AAA.

The link between DNA methylation, obesity, and adiposity-related diseases in the general population remains uncertain.