DNA methylation changes in Down syndrome derived neural iPSCs uncover co-dysregulation of ZNF and HOX3 families of transcription factors

Integrative multi-omic analysis reveals conserved cell-projection deficits in human Down syndrome brains

Down syndrome (DS) is the most common genetic cause of cognitive disability. However, it is largely unclear how triplication of a small gene subset may impinge on diverse aspects of DS brain physiopathology. Here, we took a multi-omic approach and simultaneously analyzed by RNA-seq and proteomics the expression signatures of two diverse regions of human postmortem DS brains. We found that the overexpression of triplicated genes triggered global expression dysregulation, differentially affecting transcripts, miRNAs, and proteins involved in both known and novel biological candidate pathways. Among the latter, we observed an alteration in RNA splicing, specifically modulating the expression of genes involved in cytoskeleton and axonal dynamics in DS brains. Accordingly, we found an alteration in axonal polarization in neurons from DS human iPSCs and mice. Thus, our study provides an integrated multilayer expression database capable of identifying new potential targets to aid in designing future clinical interventions for DS.

Genome-wide epigenetic signatures facilitated the variant classification of the PURA gene and uncovered the pathomechanism of PURA-related neurodevelopmental disorders

PURPOSE

Rare genetic variants in the PURA gene cause the PURA-related neurodevelopmental disorder (PURA-NDD), characterized by neonatal abnormalities and developmental delay. Using genome-wide DNA methylation analysis on patients with PURA variants, we aim to establish a PURA-NDD-specific methylation profile and provide further insights on the molecular basis of the PURA-NDD.

METHODS

Twenty three individuals (including 12 unpublished) carrying PURA variants were enrolled. We conducted the Illumina Infinium EPIC microarray analysis in 17 PURA-NDD individuals. In vitro experiments were performed to examine how PURA variants affect Pur-a expression.

RESULTS

Additional phenotypes in 12 newly identified patients were described in this study. Genome-wide DNA methylation analysis unveiled distinctive methylation profiles to PURA-NDD, and the established classifier can reclassify PURA variants of uncertain significance. Patients bearing PURA hapoloinsufficient and missense variants have comparable DNA methylation profiles, and cells expressing these PURA variants showed consistent Pur-a downregulation, suggesting a haploinsufficiency mechanism.

CONCLUSION

Patients with PURA-NDD exhibit a specific episignature, which has potential to aid identification and diagnosis of PURA-NDD patients and offer implications for further functional investigations.

Identification of lncRNA-miRNA-mRNA Networks in the Lenticular Nucleus Region of the Brain Contributes to Hepatolenticular Degeneration Pathogenesis and Therapy

Long non-coding RNAs (lncRNAs) are a recently discovered group of non-coding RNAs that play a crucial role in the regulation of various human diseases, especially in the study of nervous system diseases which has garnered significant attention. However, there is limited knowledge on the identification and function of lncRNAs in hepatolenticular degeneration (HLD). The objective of this study was to identify novel lncRNAs and determine their involvement in the networks associated with HLD. We conducted a comprehensive analysis of RNA sequencing (RNA-seq) data, reverse transcription-quantitative polymerase chain reaction (RT-qPCR), and computational biology to identify novel lncRNAs and explore their potential mechanisms in HLD. We identified 212 differently expressed lncRNAs, with 98 upregulated and 114 downregulated. Additionally, 32 differently expressed mRNAs were found, with 15 upregulated and 17 downregulated. We obtained a total of 1131 pairs of co-expressed lncRNAs and mRNAs by Pearson correlation test and prediction and annotation of the lncRNA-targeted miRNA-mRNA network. The differential lncRNAs identified in this study were found to be involved in various biological functions and signaling pathways. These include translational initiation, motor learning, locomotors behavior, dioxygenase activity, integral component of postsynaptic membrane, neuroactive ligand-receptor interaction, nuclear factor-kappa B (NF-κB) signaling pathway, cholinergic synapse, sphingolipid signaling pathway, and Parkinson's disease signaling pathway, as revealed by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Six lncRNAs, including XR_001782921.1 (P < 0.01), XR_ 001780581.1 (P < 0.01), ENSMUST_00000207119 (P < 0.01), XR_865512.2 (P < 0.01), TCONS_00005916 (P < 0.01), and TCONS_00020683 (P < 0.01), showed significant differences in expression levels between the model group and normal group by RT-qPCR. Among these, four lncRNAs (TCONS_00020683, XR_865512.2, XR_001780581.1, and ENSMUST00000207119) displayed a high degree of conservation. This study provides a unique perspective for the pathogenesis and therapy of HLD by constructing the lncRNA-miRNA-mRNA network. This insight provides a foundation for future exploration in this field.

Effect of Genetic and Dietary Perturbation of Glycine Metabolism on Atherosclerosis in Humans and Mice

Objective

Epidemiological and genetic studies have reported inverse associations between circulating glycine levels and risk of coronary artery disease (CAD). However, these findings have not been consistently observed in all studies. We sought to evaluate the causal relationship between circulating glycine levels and atherosclerosis using large-scale genetic analyses in humans and dietary supplementation experiments in mice.

Methods

Serum glycine levels were evaluated for association with prevalent and incident CAD in the UK Biobank. A multi-ancestry genome-wide association study (GWAS) meta-analysis was carried out to identify genetic determinants for circulating glycine levels, which were then used to evaluate the causal relationship between glycine and risk of CAD by Mendelian randomization (MR). A glycine feeding study was carried out with atherosclerosis-prone apolipoprotein E deficient (ApoE-/-) mice to determine the effects of increased circulating glycine levels on amino acid metabolism, metabolic traits, and aortic lesion formation.

Results

Among 105,718 subjects from the UK Biobank, elevated serum glycine levels were associated with significantly reduced risk of prevalent CAD (Quintile 5 vs. Quintile 1 OR=0.76, 95% CI 0.67-0.87; P<0.0001) and incident CAD (Quintile 5 vs. Quintile 1 HR=0.70, 95% CI 0.65-0.77; P<0.0001) in models adjusted for age, sex, ethnicity, anti-hypertensive and lipid-lowering medications, blood pressure, kidney function, and diabetes. A meta-analysis of 13 GWAS datasets (total n=230,947) identified 61 loci for circulating glycine levels, of which 26 were novel. MR analyses provided modest evidence that genetically elevated glycine levels were causally associated with reduced systolic blood pressure and risk of type 2 diabetes, but did provide evidence for an association with risk of CAD. Furthermore, glycine-supplementation in ApoE-/- mice did not alter cardiometabolic traits, inflammatory biomarkers, or development of atherosclerotic lesions.

Conclusions

Circulating glycine levels were inversely associated with risk of prevalent and incident CAD in a large population-based cohort. While substantially expanding the genetic architecture of circulating glycine levels, MR analyses and in vivo feeding studies in humans and mice, respectively, did not provide evidence that the clinical association of this amino acid with CAD represents a causal relationship, despite being associated with two correlated risk factors.

Trisomies Reorganize Human 3D Genome

Trisomy is the presence of one extra copy of an entire chromosome or its part in a cell nucleus. In humans, autosomal trisomies are associated with severe developmental abnormalities leading to embryonic lethality, miscarriage or pronounced deviations of various organs and systems at birth. Trisomies are characterized by alterations in gene expression level, not exclusively on the trisomic chromosome, but throughout the genome. Here, we applied the high-throughput chromosome conformation capture technique (Hi-C) to study chromatin 3D structure in human chorion cells carrying either additional chromosome 13 (Patau syndrome) or chromosome 16 and in cultured fibroblasts with extra chromosome 18 (Edwards syndrome). The presence of extra chromosomes results in systematic changes of contact frequencies between small and large chromosomes. Analyzing the behavior of individual chromosomes, we found that a limited number of chromosomes change their contact patterns stochastically in trisomic cells and that it could be associated with lamina-associated domains (LAD) and gene content. For trisomy 13 and 18, but not for trisomy 16, the proportion of compacted loci on a chromosome is correlated with LAD content. We also found that regions of the genome that become more compact in trisomic cells are enriched in housekeeping genes, indicating a possible decrease in chromatin accessibility and transcription level of these genes. These results provide a framework for understanding the mechanisms of pan-genome transcription dysregulation in trisomies in the context of chromatin spatial organization.

Transition from Animal-Based to Human Induced Pluripotent Stem Cells (iPSCs)-Based Models of Neurodevelopmental Disorders: Opportunities and Challenges

Neurodevelopmental disorders (NDDs) arise from the disruption of highly coordinated mechanisms underlying brain development, which results in impaired sensory, motor and/or cognitive functions. Although rodent models have offered very relevant insights to the field, the translation of findings to clinics, particularly regarding therapeutic approaches for these diseases, remains challenging. Part of the explanation for this failure may be the genetic differences-some targets not being conserved between species-and, most importantly, the differences in regulation of gene expression. This prompts the use of human-derived models to study NDDS. The generation of human induced pluripotent stem cells (hIPSCs) added a new suitable alternative to overcome species limitations, allowing for the study of human neuronal development while maintaining the genetic background of the donor patient. Several hIPSC models of NDDs already proved their worth by mimicking several pathological phenotypes found in humans. In this review, we highlight the utility of hIPSCs to pave new paths for NDD research and development of new therapeutic tools, summarize the challenges and advances of hIPSC-culture and neuronal differentiation protocols and discuss the best way to take advantage of these models, illustrating this with examples of success for some NDDs.

Integrative identification of hub genes in development of atrial fibrillation related stroke

BACKGROUND

As the most common arrhythmia, atrial fibrillation (AF) is associated with a significantly increased risk of stroke, which causes high disability and mortality. To date, the underlying mechanism of stroke occurring after AF remains unclear. Herein, we studied hub genes and regulatory pathways involved in AF and secondary stroke and aimed to reveal biomarkers and therapeutic targets of AF-related stroke.

METHODS

The GSE79768 and GSE58294 datasets were used to analyze AF- and stroke-related differentially expressed genes (DEGs) to obtain a DEG1 dataset. Weighted correlation network analysis (WGCNA) was used to identify modules associated with AF-related stroke in GSE66724 (DEG2). DEG1 and DEG2 were merged, and hub genes were identified based on protein-protein interaction networks. Gene Ontology terms were used to analyze the enriched pathways. The GSE129409 and GSE70887 were applied to construct a circRNA-miRNA-mRNA network in AF-related stroke. Hub genes were verified in patients using quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

We identified 3,132 DEGs in blood samples and 253 DEGs in left atrial specimens. Co-expressed hub genes of EIF4E3, ZNF595, ZNF700, MATR3, ACKR4, ANXA3, SEPSECS-AS1, and RNF166 were significantly associated with AF-related stroke. The hsa_circ_0018657/hsa-miR-198/EIF4E3 pathway was explored as the regulating axis in AF-related stroke. The qRT-PCR results were consistent with the bioinformatic analysis.

CONCLUSIONS

Hub genes EIF4E3, ZNF595, ZNF700, MATR3, ACKR4, ANXA3, SEPSECS-AS1, and RNF166 have potential as novel biomarkers and therapeutic targets in AF-related stroke. The hsa_circ_0018657/hsa-miR-198/EIF4E3 axis could play an important role regulating the development of AF-related stroke.

Future Perspectives in Oxidative Stress in Trisomy 13 and 18 Evaluation

Autosomal aneuploidies are the most frequently occurring congenital abnormalities and are related to many metabolic disorders, hormonal dysfunctions, neurotransmitter abnormalities, and intellectual disabilities. Trisomies are generated by an error of chromosomal segregation during cell division. Accumulating evidence has shown that deregulated gene expression resulting from the triplication of chromosomes 13 and 18 is associated with many disturbed cellular processes. Moreover, a disturbed oxidative stress status may be implicated in the occurrence of fetal malformations. Therefore, a literature review was undertaken to provide novel insights into the evaluation of trisomy 13 (T13) and 18 (T18) pathogeneses, with a particular concern on the oxidative stress. Corresponding to the limited literature data focused on factors leading to T13 and T18 phenotype occurrence, the importance of oxidative stress evaluation in T13 and T18 could enable the determination of subsequent disturbed metabolic pathways, highlighting the related role of mitochondrial dysfunction or epigenetics. This review illustrates up-to-date T13 and T18 research and discusses the strengths, limitations, and possible directions for future studies. The progressive unification of trisomy-related research protocols might provide potential medical targets in the future along with the implementation of the foundation of modern prenatal medicine.

Cell models for Down syndrome-Alzheimer’s disease research

Down syndrome (DS) is the most common chromosomal abnormality and leads to intellectual disability, increased risk of cardiac defects, and an altered immune response. Individuals with DS have an extra full or partial copy of chromosome 21 (trisomy 21) and are more likely to develop early-onset Alzheimer’s disease (AD) than the general population. Changes in expression of human chromosome 21 (Hsa21)-encoded genes, such as amyloid precursor protein (APP), play an important role in the pathogenesis of AD in DS (DS-AD). However, the mechanisms of DS-AD remain poorly understood. To date, several mouse models with an extra copy of genes syntenic to Hsa21 have been developed to characterise DS-AD-related phenotypes. Nonetheless, due to genetic and physiological differences between mouse and human, mouse models cannot faithfully recapitulate all features of DS-AD. Cells differentiated from human-induced pluripotent stem cells (iPSCs), isolated from individuals with genetic diseases, can be used to model disease-related cellular and molecular pathologies, including DS. In this review, we will discuss the limitations of mouse models of DS and how these can be addressed using recent advancements in modelling DS using human iPSCs and iPSC-mouse chimeras, and potential applications of iPSCs in preclinical studies for DS-AD.

The Role of Janus Kinase/STAT3 Pathway in Hematologic Malignancies With an Emphasis on Epigenetics

The Janus kinase (JAK)/signal transducer and activator of transcription (STAT) signaling pathway has been known to be involved in cell growth, cellular differentiation processes development, immune cell survival, and hematopoietic system development. As an important member of the STAT family, STAT3 participates as a major regulator of cellular development and differentiation-associated genes. Prolonged and persistent STAT3 activation has been reported to be associated with tumor cell survival, proliferation, and invasion. Therefore, the JAK-STAT pathway can be a potential target for drug development to treat human cancers, e.g., hematological malignancies. Although STAT3 upregulation has been reported in hematopoietic cancers, protein-level STAT3 mutations have also been reported in invasive leukemias/lymphomas. The principal role of STAT3 in tumor cell growth clarifies the importance of approaches that downregulate this molecule. Epigenetic modifications are a major regulatory mechanism controlling the activity and function of STAT3. So far, several compounds have been developed to target epigenetic regulatory enzymes in blood malignancies. Here, we discuss the current knowledge about STAT3 abnormalities and carcinogenic functions in hematopoietic cancers, novel STAT3 inhibitors, the role of epigenetic mechanisms in STAT3 regulation, and targeted therapies, by focusing on STAT3-related epigenetic modifications.

Promoter methylation-mediated repression of UNC5 receptors and the associated clinical significance in human colorectal cancer

Background

Deregulated methylation of tumor suppressor genes is a hallmark event in colorectal cancer (CRC) carcinogenesis. UNC5 receptors, down-regulated in various human malignancies due to epigenetic alterations, have been proposed as putative tumor suppressor genes. In this study, we focused on the methylation-mediated inhibition of UNC5 receptors and the associated clinical significance in CRC.

Methods

Methylation and expression analysis was performed in TCGA datasets. And the results were confirmed in vitro in CRC cell lines treated with 5-aza-deoxycytidine. Then, the expression and epigenetic alterations of UNC5 receptors were evaluated in clinical specimens. Moreover, the diagnostic and prognostic values of the methylation alterations were also analyzed.

Results

Methylation-mediated repression was observed in UNC5C and UNC5D, but not in UNC5A and UNC5B, which was confirmed in CRC cell lines. Except for UNC5B, significantly elevated methylation was observed in UNC5A, UNC5C, and UNC5D in CRC. The discrimination efficiency of the three receptors was comparable with that of SEPT9. Kaplan–Meier curve survival analysis showed that hypermethylation of UNC5A, UNC5C and UNC5D was associated with poor progression-free and overall survival. Moreover, methylation levels of UNC5C and UNC5D were independent predictors of CRC progression-free (P = 0.001, P = 0.003, respectively) and overall survival (P = 0.008, P = 0.004, respectively).

Conclusions

Hypermethylation of UNC5C and UNC5D mediates the repression and has promising diagnostic and prognostic values in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01211-5.

DNA Methylation Alterations in Blood Cells of Toddlers with Down Syndrome

Recent research has provided evidence on genome-wide alterations in DNA methylation patterns due to trisomy 21, which have been detected in various tissues of individuals with Down syndrome (DS) across different developmental stages. Here, we report new data on the systematic genome-wide DNA methylation perturbations in blood cells of individuals with DS from a previously understudied age group—young children. We show that the study findings are highly consistent with those from the prior literature. In addition, utilizing relevant published data from two other developmental stages, neonatal and adult, we track a quasi-longitudinal trend in the DS-associated DNA methylation patterns as a systematic epigenomic destabilization with age.

Impact of the Olig Family on Neurodevelopmental Disorders

The Olig genes encode members of the basic helix-loop-helix (bHLH) family of transcription factors. Olig1, Olig2, and Olig3 are expressed in both the developing and mature central nervous system (CNS) and strictly regulate cellular specification and differentiation. Extensive studies have established functional roles of Olig1 and Olig2 in directing neuronal and glial formation during different stages in development. Recently, Olig2 overexpression was implicated in neurodevelopmental disorders down syndrome (DS) and autism spectrum disorder (ASD) but its influence on cognitive and intellectual defects remains unknown. In this review, we summarize the biological functions of the Olig family and how it uniquely promotes cellular diversity in the CNS. This is followed up with a discussion on how abnormal Olig2 expression impacts brain development and function in DS and ASD. Collectively, the studies described here emphasize vital features of the Olig members and their distinctive potential roles in neurodevelopmental disease states.

Stable DNMT3L overexpression in SH-SY5Y neurons recreates a facet of the genome-wide Down syndrome DNA methylation signature

Background

Down syndrome (DS) is characterized by a genome-wide profile of differential DNA methylation that is skewed towards hypermethylation in most tissues, including brain, and includes pan-tissue differential methylation. The molecular mechanisms involve the overexpression of genes related to DNA methylation on chromosome 21. Here, we stably overexpressed the chromosome 21 gene DNA methyltransferase 3L (DNMT3L) in the human SH-SY5Y neuroblastoma cell line and assayed DNA methylation at over 26 million CpGs by whole genome bisulfite sequencing (WGBS) at three different developmental phases (undifferentiated, differentiating, and differentiated).

Results

DNMT3L overexpression resulted in global CpG and CpG island hypermethylation as well as thousands of differentially methylated regions (DMRs). The DNMT3L DMRs were skewed towards hypermethylation and mapped to genes involved in neurodevelopment, cellular signaling, and gene regulation. Consensus DNMT3L DMRs showed that cell lines clustered by genotype and then differentiation phase, demonstrating sets of common genes affected across neuronal differentiation. The hypermethylated DNMT3L DMRs from all pairwise comparisons were enriched for regions of bivalent chromatin marked by H3K4me3 as well as differentially methylated sites from previous DS studies of diverse tissues. In contrast, the hypomethylated DNMT3L DMRs from all pairwise comparisons displayed a tissue-specific profile enriched for regions of heterochromatin marked by H3K9me3 during embryonic development.

Conclusions

Taken together, these results support a mechanism whereby regions of bivalent chromatin that lose H3K4me3 during neuronal differentiation are targeted by excess DNMT3L and become hypermethylated. Overall, these findings demonstrate that DNMT3L overexpression during neurodevelopment recreates a facet of the genome-wide DS DNA methylation signature by targeting known genes and gene clusters that display pan-tissue differential methylation in DS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00387-7.

DNMT3L promotes neural differentiation by enhancing STAT1 and STAT3 phosphorylation independent of DNA methylation

Previously, we reported global hypermethylation in DS might be attributed to the overexpression of HSA21 gene DNMT3L, which can enhance DNMT3A and DNMT3B activities in DNA methylation. To test this hypothesis, we compared the DNA methylation and RNA expression profiles of early-differentiated human neuroprogenitors with and without DNMT3L overexpression. We found DNMT3L overexpression only moderately increased the DNA methylation of limited genes, yet significantly altered global RNA expression of genes involved in neural differentiation. We further found that DNMT3L bound STAT1 or STAT3, and increased its phosphorylation and nuclear translocation, which in turn activated the expression of transcription factors including HES3, ASCL1, NEUROD2 and NEUROG2 and CDK inhibitor CDKN1A, which promoted cell cycle exit and neural differentiation. This phenomenon was also confirmed in Dnmt3l conditional knockin mice, which could be rescued by STAT1 and STAT3 phosphorylation inhibitors (Fludarabine and SH-4-54) but not DNA methylation inhibitor (Decitabine). These results suggest that DNMT3L play an important role during neurodevelopment independent of DNA methylation, which may contribute to the abnormal phenotypes observed in Down syndrome cortex.

Generation of two human iPSC lines (UUIGPi013-A and UUIPGi014-A) from cases with Down syndrome and full trisomy for chromosome 21 (T21)

Down syndrome (DS) is caused by trisomy for chromosome 21 (T21). We generated two induced pluripotent stem cell (iPSC) lines from skin fibroblasts of two males with DS using Sendai virus delivery of OCT4, SOX2, KLF4, and c-MYC. Characterization of the two iPSC lines, UUIGPi013-A and UUIPGi014-A, showed that they are genetically stable with a 47,XY,+21 karyotype. Both lines displayed expression of pluripotency markers and trilineage differentiation capacity. These two iPSC lines provide a useful resource for DS modeling and pharmacological interventions.

Low-pass whole genome bisulfite sequencing of neonatal dried blood spots identifies a role for RUNX1 in Down syndrome DNA methylation profiles

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables retrospective investigation into early life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q < 0.05) differentially methylated regions (DMRs) that distinguished DS from typical development and idiopathic developmental delay were identified. Machine learning feature selection refined these DMRs to 22 loci. The DS DMRs mapped to genes involved in neurodevelopment, metabolism, and transcriptional regulation. Based on comparisons with previous DS methylation studies and reference epigenomes, the hypermethylated DS DMRs were significantly (q < 0.05) enriched across tissues while the hypomethylated DS DMRs were significantly (q < 0.05) enriched for blood-specific chromatin states. A ~28 kb block of hypermethylation was observed on chromosome 21 in the RUNX1 locus, which encodes a hematopoietic transcription factor whose binding motif was the most significantly enriched (q < 0.05) overall and specifically within the hypomethylated DMRs. Finally, we also identified DMRs that distinguished DS NDBS based on the presence or absence of congenital heart disease (CHD). Together, these results not only demonstrate the utility of low-pass WGBS on NDBS samples for epigenome-wide association studies, but also provide new insights into the early life mechanisms of epigenomic dysregulation resulting from trisomy 21.