Noncoding RNA in the transcriptional landscape of human neural progenitor cell differentiation

Transition of allele-specific DNA hydroxymethylation at regulatory loci is associated with phenotypic variation in monozygotic twins discordant for psychiatric disorders

BACKGROUND

Major psychiatric disorders such as schizophrenia (SCZ) and bipolar disorder (BPD) are complex genetic mental illnesses. Their non-Mendelian features, such as those observed in monozygotic twins discordant for SCZ or BPD, are likely complicated by environmental modifiers of genetic effects. 5-Hydroxymethylcytosine (5hmC) is an important epigenetic mark in gene regulation, and whether it is linked to genetic variants that contribute to non-Mendelian features remains largely unexplored.

METHODS

We combined the 5hmC-selective chemical labeling method (5hmC-seq) and whole-genome sequencing (WGS) analysis of peripheral blood DNA obtained from monozygotic (MZ) twins discordant for SCZ or BPD to identify allelic imbalances in hydroxymethylome maps, and examined association of allele-specific hydroxymethylation (AShM) transition with disease susceptibility based on Bayes factors (BF) derived from the Bayesian generalized additive linear mixed model. We then performed multi-omics integrative analysis to determine the molecular pathogenic basis of those AShM sites. We finally employed luciferase reporter, CRISPR/Cas9 technology, electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP), PCR, FM4-64 imaging analysis, and RNA sequencing to validate the function of interested AShM sites in the human neuroblastoma SK-N-SH cells and human embryonic kidney 293T (HEK293T) cells.

RESULTS

We identified thousands of genetic variants associated with AShM imbalances that exhibited phenotypic variation-associated AShM changes at regulatory loci. These AShM marks showed plausible associations with SCZ or BPD based on their effects on interactions among transcription factors (TFs), DNA methylation levels, or other epigenomic marks and thus contributed to dysregulated gene expression, which ultimately increased disease susceptibility. We then validated that competitive binding of POU3F2 on the alternative allele at the AShM site rs4558409 (G/T) in PLLP-enhanced PLLP expression, while the hydroxymethylated alternative allele, which alleviated the POU3F2 binding activity at the rs4558409 site, might be associated with the downregulated PLLP expression observed in BPD or SCZ. Moreover, disruption of rs4558409 promoted neural development and vesicle trafficking.

CONCLUSION

Our study provides a powerful strategy for prioritizing regulatory risk variants and contributes to our understanding of the interplay between genetic and epigenetic factors in mediating SCZ or BPD susceptibility.

Integrated RNA Sequencing Reveals Epigenetic Impacts of Diesel Particulate Matter Exposure in Human Cerebral Organoids

Autism spectrum disorder (ASD) manifests early in childhood. While genetic variants increase risk for ASD, a growing body of literature has established that in utero chemical exposures also contribute to ASD risk. These chemicals include air-based pollutants like diesel particulate matter (DPM). A combination of single-cell and direct transcriptomics of DPM-exposed human-induced pluripotent stem cell-derived cerebral organoids revealed toxicogenomic effects of DPM exposure during fetal brain development. Direct transcriptomics, sequencing RNA bases via Nanopore, revealed that cerebral organoids contain extensive RNA modifications, with DPM-altering cytosine methylation in oxidative mitochondrial transcripts expressed in outer radial glia cells. Single-cell transcriptomics further confirmed an oxidative phosphorylation change in cell groups such as outer radial glia upon DPM exposure. This approach highlights how DPM exposure perturbs normal mitochondrial function and cellular respiration during early brain development, which may contribute to developmental disorders like ASD by altering neurodevelopment.

The Long Noncoding RNA RPS10P2-AS1 Is Implicated in Autism Spectrum Disorder Risk and Modulates Gene Expression in Human Neuronal Progenitor Cells

Most of the genetic risk for autism spectrum disorder (ASD) is inherited as common genetic variants, although some rare mutations have been identified in individuals with ASD. Common genetic variants are most parsimoniously identified by genome wide association studies. Genome wide association studies have identified several genetic loci with genome wide association with ASD. However, genome wide association studies only identify regions of the genome associated with phenotypic traits. Identification of the functional elements requires additional experimental evidence. Here, we demonstrate that a genome wide association study locus for ASD on chromosome 20p12.1, rs4141463, implicates a noncoding RNA as a functional element. Although rs4141463 lies within an intron of the protein-coding MACROD2 (MACRO domain containing 2) gene, expression of MACROD2 is neither altered in postmortem temporal cortex of individuals with ASD nor correlated with rs4141463 genotype. Our bioinformatics approaches revealed a noncoding RNA transcript near the autism susceptibility signal, RPS10P2-AS1 (ribosomal protein S10 pseudogene 2 anti-sense 1). In a panel of 15 human tissues, RPS10P2-AS1 was expressed at higher levels than the protein-coding MACROD2 in both fetal temporal cortex and adult peripheral blood. In postmortem temporal cortex, expression of RPS10P2-AS1 was increased 7-fold in individuals with ASD (P = 0.02) and increased 8-fold in individuals with the ASD-associated rs4141463 genotype (P = 0.01). Further, RPS10P2-AS1 expression was increased in human neural progenitor cells exposed to model air pollutants, indicating that both genetic and environmental factors that contribute to ASD increased RPS10P2-AS1 expression. Overexpression of RPS10P2-AS1 in human neural progenitor cells indicated substantial changes in neuronal gene expression. These data indicate that genome-wide significant associations with ASD implicate long noncoding RNAs. Because long noncoding RNAs are more abundant in human brain than protein-coding RNAs, this class of molecules is likely to contribute to ASD risk.

Long noncoding RNA TCONS_00024652 regulates vascular endothelial cell proliferation and angiogenesis via microRNA-21

Acute coronary syndrome caused by the rupture of atherosclerotic plaques is one of the primary causes of major cardiovascular events, and neovascularization within the plaque is closely associated with its stability. Long noncoding RNA (lncRNAs) is a type of noncoding RNA that serves a crucial role in regulating vascular endothelial cells (VECs). The aim of the present study was to investigate the effect of lncRNA TCONS_00024652 on the proliferation and angiogenesis of VECs following stimulation with TNF-α. The expression of lncRNA and miRNA was measured in human umbilical vein endothelial cells (HUVECs) by reverse transcription-quantitative polymerase chain reaction. Cell proliferation was measured using a Cell Counting Kit-8 assay. Wound healing and tube formation assays were performed to determine cell migration and angiogenesis. Interactions between TCONS_00024652 and miR-21 were determined using bioinformatics and a dual-luciferase reporter assay. The results demonstrated that TCONS_00024652 is highly expressed in TNF-α-induced HUVECs. Functional assays demonstrated that the dysregulated expression of TCONS_00024652 promotes endothelial cell proliferation and angiogenesis, whereas TCONS_00024652 knockdown induces the opposite effects. Bioinformatics analysis using starBase predicted putative binding at the 3′-untranslated region of TCONS_00024652 and miR-21 and luciferase reporter assays further verified this interaction. The results of the present study suggest that the targeting of TCONS_00024652 by miR-21 may be a potential method of improving vascular endothelial dysfunction, neovascularization maturation and plaque stabilization.

Expression profiling and functional annotation of noncoding genes across 11 distinct organs in rat development

Accumulating evidence suggests that noncoding RNAs (ncRNAs) have important regulatory functions. However, lacking of functional annotations for ncRNAs hampered us from carrying out the subsequent functional or predictive research. Here we dissected the expression profiles of 3,458 rat noncoding genes using rat bodymap RNA-sequencing data consisting of 11 solid organs over four developmental stages (juvenile, adolescent, adult and aged) from both sexes, and conducted a comprehensive analysis of differentially expressed noncoding genes (DEnGs) between various conditions. We then constructed a co-expression network between protein-coding and noncoding genes to infer biological functions of noncoding genes. Modules of interest were linked to online databases including DAVID for functional annotation and pathway analysis. Our results indicated that noncoding genes are functionally enriched through pathways similar to those of protein-coding genes. Terms about development of the immune system were enriched with genes from age-related modules, whereas terms about sexual reproduction were enriched with genes in sex-related modules. We also built connection networks on some significant modules to visualize the interactions and regulatory relationship between protein-coding and noncoding genes. Our study could improve our understanding and facilitate a deeper investigation on organ/age/sex-related regulatory events of noncoding genes, which may lead to a superior preclinical model for drug development and translational medicine.

Structure Prediction: New Insights into Decrypting Long Noncoding RNAs

Long noncoding RNAs (lncRNAs), which form a diverse class of RNAs, remain the least understood type of noncoding RNAs in terms of their nature and identification. Emerging evidence has revealed that a small number of newly discovered lncRNAs perform important and complex biological functions such as dosage compensation, chromatin regulation, genomic imprinting, and nuclear organization. However, understanding the wide range of functions of lncRNAs related to various processes of cellular networks remains a great experimental challenge. Structural versatility is critical for RNAs to perform various functions and provides new insights into probing the functions of lncRNAs. In recent years, the computational method of RNA structure prediction has been developed to analyze the structure of lncRNAs. This novel methodology has provided basic but indispensable information for the rapid, large-scale and in-depth research of lncRNAs. This review focuses on mainstream RNA structure prediction methods at the secondary and tertiary levels to offer an additional approach to investigating the functions of lncRNAs.