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

Long noncoding RNA MALAT1 mediates fibrous topography-driven pathologic calcification through trans-differentiation of myoblasts

Prosthesis-induced pathological calcification is a significant challenge in biomaterial applications and is often associated with various reconstructive medical procedures. It is uncertain whether the fibrous extracellular matrix (ECM) adjacent to biomaterials directly triggers osteogenic trans-differentiation in nearby cells. To investigate this possibility, we engineered a heterogeneous polystyrene fibrous matrix (PSF) designed to mimic the ECM. Our findings revealed that the myoblasts grown on this PSF acquired osteogenic properties, resulting in mineralization both in vitro and in vivo. Transcriptomic analyses indicated a notable upregulation in the expression of the long noncoding RNA metastsis-associated lung adenocarcinoma transcript 1 (Malat1) in the C2C12 myoblasts cultured on PSF. Intriguingly, silencing Malat1 curtailed the PSF-induced mineralization and downregulated the expression of bone morphogenetic proteins (Bmps) and osteogenic markers. Further, we found that PSF prompted the activation of Yap1 signaling and epigenetic modifications in the Malat1 promoter, crucial for the expression of Malat1. These results indicate that the fibrous matrix adjacent to biomaterials can instigate Malat1 upregulation, subsequently driving osteogenic trans-differentiation in myoblasts and ectopic calcification through its transcriptional regulation of osteogenic genes, including Bmps. Our findings point to a novel therapeutic avenue for mitigating prosthesis-induced pathological calcification, heralding new possibilities in the field of biomaterial-based therapies.

Single generation epigenetic change in captivity and reinforcement in subsequent generations in a delta smelt (Hypomesus transpacificus) conservation hatchery

A refugial population of the endangered delta smelt (Hypomesus transpacificus) has been maintained at the Fish Conservation and Culture Laboratory (FCCL) at UC Davis since 2008. Despite intense genetic management, fitness differences between wild and cultured fish have been observed at the FCCL. To investigate the molecular underpinnings of hatchery domestication, we used whole-genome bisulfite sequencing to quantify epigenetic differences between wild and hatchery-origin delta smelt. Differentially methylated regions (DMRs) were identified from 104 individuals by comparing the methylation patterns in different generations of hatchery fish (G1, G2, G3) with their wild parents (G0). We discovered a total of 132 significant DMRs (p < .05) between G0 and G1, 132 significant DMRs between G0 and G2, and 201 significant DMRs between G0 and G3. Our results demonstrate substantial differences in methylation patterns emerged between the wild and hatchery-reared fish in the early generations in the hatchery, with a higher proportion of hypermethylated DMRs in hatchery-reared fish. The rearing environment was found to be a stronger predictor of individual clustering based on methylation patterns than family, sex or generation. Our study indicates a reinforcement of the epigenetic status with successive generations in the hatchery environment, as evidenced by an increase in methylation in hypermethylated DMRs and a decrease in methylation in hypomethylated DMRs over time. Lastly, our results demonstrated heterogeneity in inherited methylation pattern in families across generations. These insights highlight the long-term consequences of hatchery practices on the epigenetic landscape, potentially impacting wild fish populations.

The genomic mosaic of mitochondrial dysfunction: Decoding nuclear and mitochondrial epigenetic contributions to maternally inherited diabetes and deafness pathogenesis

Aims

Maternally inherited diabetes and deafness (MIDD) is a complex disorder characterized by multiorgan clinical manifestations, including diabetes, hearing loss, and ophthalmic complications. This pilot study aimed to elucidate the intricate interplay between nuclear and mitochondrial genetics, epigenetic modifications, and their potential implications in the pathogenesis of MIDD.

Main methods

A comprehensive genomic approach was employed to analyze a Sicilian family affected by clinically characterized MIDD, negative to the only known causative m.3243 A > G variant, integrating whole-exome sequencing and whole-genome bisulfite sequencing of both nuclear and mitochondrial analyses.

Key findings

Rare and deleterious variants were identified across multiple nuclear genes involved in retinal homeostasis, mitochondrial function, and epigenetic regulation, while complementary mitochondrial DNA analysis revealed a rich tapestry of genetic diversity across genes encoding components of the electron transport chain and ATP synthesis machinery. Epigenetic analyses uncovered significant differentially methylated regions across the genome and within the mitochondrial genome, suggesting a nuanced landscape of epigenetic modulation.

Significance

The integration of genetic and epigenetic data highlighted the potential crosstalk between nuclear and mitochondrial regulation, with specific mtDNA variants influencing methylation patterns and potentially impacting the expression and regulation of mitochondrial genes. This pilot study provides valuable insights into the complex molecular mechanisms underlying MIDD, emphasizing the interplay between nucleus and mitochondrion, tracing the way for future research into targeted therapeutic interventions and personalized approaches for disease management.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.

Mitochondrial DNA copy number reduction via in vitro TFAM knockout remodels the nuclear epigenome and transcriptome

Mitochondrial DNA copy number (mtDNA-CN) is associated with several age-related chronic diseases and is a predictor of all-cause mortality. Here, we examine site-specific differential nuclear DNA (nDNA) methylation and differential gene expression resulting from in vitro reduction of mtDNA-CN to uncover shared genes and biological pathways mediating the effect of mtDNA-CN on disease. Epigenome and transcriptome profiles were generated for three independent human embryonic kidney (HEK293T) cell lines harbouring a mitochondrial transcription factor A (TFAM) heterozygous knockout generated via CRISPR-Cas9, and matched control lines. We identified 4,242 differentially methylated sites, 228 differentially methylated regions, and 179 differentially expressed genes associated with mtDNA-CN. Integrated analysis uncovered 381 Gene-CpG pairs. GABAA receptor genes and related pathways, the neuroactive ligand receptor interaction pathway, ABCD1/2 gene activity, and cell signalling processes were overrepresented, providing insight into the underlying biological mechanisms facilitating these associations. We also report evidence implicating chromatin state regulatory mechanisms as modulators of mtDNA-CN effect on gene expression. We demonstrate that mitochondrial DNA variation signals to the nuclear DNA epigenome and transcriptome and may lead to nuclear remodelling relevant to development, aging, and complex disease.

Artificial intelligence for prenatal chromosome analysis

This review article delves into the rapidly advancing domain of prenatal diagnostics, with a primary focus on the detection and management of chromosomal abnormalities such as trisomy 13 ("Patau syndrome)", "trisomy 18 (Edwards syndrome)", and "trisomy 21 (Down syndrome)". The objective of the study is to examine the utilization and effectiveness of novel computational methodologies, such as "machine learning (ML)", "deep learning (DL)", and data analysis, in enhancing the detection rates and accuracy of these prenatal conditions. The contribution of the article lies in its comprehensive examination of advancements in "Non-Invasive Prenatal Testing (NIPT)", prenatal screening, genomics, and medical imaging. It highlights the potential of these techniques for prenatal diagnosis and the contributions of ML and DL to these advancements. It highlights the application of ensemble models and transfer learning to improving model performance, especially with limited datasets. This also delves into optimal feature selection and fusion of high-dimensional features, underscoring the need for future research in these areas. The review finds that ML and DL have substantially improved the detection and management of prenatal conditions, despite limitations such as small sample sizes and issues related to model generalizability. It recognizes the promising results achieved through the use of ensemble models and transfer learning in prenatal diagnostics. The review also notes the increased importance of feature selection and high-dimensional feature fusion in the development and training of predictive models. The findings underline the crucial role of AI and machine learning techniques in early detection and improved therapeutic strategies in prenatal diagnostics, highlighting a pressing need for further research in this area.

TET1 regulates responses to house dust mite by altering chromatin accessibility, DNA methylation, and gene expression in airway epithelial cells

Background

Previous studies have identified TET1 as a potential key regulator of genes linked to asthma. TET1 has been shown to transcriptionally respond to house dust mite extract, an allergen known to directly cause allergic asthma development, and regulate the expression of genes involved in asthma. How TET1 regulates expression of these genes, however, is unknown. TET1 is a DNA demethylase; therefore, most prior research on TET1-based gene regulation has focused on how TET1 affects methylation. However, TET1 can also interact directly with transcription factors and histone modifiers to regulate gene expression. Understanding how TET1 regulates expression to contribute to allergic responses and asthma development thus requires a comprehensive approach. To this end, we measured mRNA expression, DNA methylation, chromatin accessibility and histone modifications in control and TET1 knockdown human bronchial epithelial cells treated or untreated with house dust mite extract.

Results

Throughout our analyses, we detected strong similarities between the effects of TET1 knockdown alone and the effects of HDM treatment alone. One especially striking pattern was that both TET1 knockdown and HDM treatment generally led to decreased chromatin accessibility at largely the same genomic loci. Transcription factor enrichment analyses indicated that altered chromatin accessibility following the loss of TET1 may affect, or be affected by, CTCF and CEBP binding. TET1 loss also led to changes in DNA methylation, but these changes were generally in regions where accessibility was not changing.

Conclusions

TET1 regulates gene expression through different mechanisms (DNA methylation and chromatin accessibility) in different parts of the genome in the airway epithelial cells, which mediates inflammatory responses to allergen. Collectively, our data suggest novel molecular mechanisms through which TET1 regulates critical pathways following allergen challenges and contributes to the development of asthma.

Epigenomic signature of major congenital heart defects in newborns with Down syndrome

BACKGROUND

Congenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS), but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.

METHODS

We used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: (1) 45 DS-CHD (27 female, 18 male) and (2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD versus DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell-type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS versus typical development (TD) WGBS NDBS samples.

RESULTS

We found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS versus TD samples.

CONCLUSIONS

A sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

RUN(X) out of blood: emerging RUNX1 functions beyond hematopoiesis and links to Down syndrome

BACKGROUND

RUNX1 is a transcription factor and a master regulator for the specification of the hematopoietic lineage during embryogenesis and postnatal megakaryopoiesis. Mutations and rearrangements on RUNX1 are key drivers of hematological malignancies. In humans, this gene is localized to the 'Down syndrome critical region' of chromosome 21, triplication of which is necessary and sufficient for most phenotypes that characterize Trisomy 21.

MAIN BODY

Individuals with Down syndrome show a higher predisposition to leukemias. Hence, RUNX1 overexpression was initially proposed as a critical player on Down syndrome-associated leukemogenesis. Less is known about the functions of RUNX1 in other tissues and organs, although growing reports show important implications in development or homeostasis of neural tissues, muscle, heart, bone, ovary, or the endothelium, among others. Even less is understood about the consequences on these tissues of RUNX1 gene dosage alterations in the context of Down syndrome. In this review, we summarize the current knowledge on RUNX1 activities outside blood/leukemia, while suggesting for the first time their potential relation to specific Trisomy 21 co-occurring conditions.

CONCLUSION

Our concise review on the emerging RUNX1 roles in different tissues outside the hematopoietic context provides a number of well-funded hypotheses that will open new research avenues toward a better understanding of RUNX1-mediated transcription in health and disease, contributing to novel potential diagnostic and therapeutic strategies for Down syndrome-associated conditions.

Epigenomic signatures reveal mechanistic clues and predictive markers for autism spectrum disorder

Autism spectrum disorder (ASD) comprises a heterogeneous group of neurodevelopmental outcomes in children with a commonality in deficits in social communication and language combined with repetitive behaviors and interests. The etiology of ASD is heterogeneous, as several hundred genes have been implicated as well as multiple in utero environmental exposures. Over the past two decades, epigenetic investigations, including DNA methylation, have emerged as a novel way to capture the complex interface of multivariate ASD etiologies. More recently, epigenome-wide association studies using human brain and surrogate accessible tissues have revealed some convergent genes that are epigenetically altered in ASD, many of which overlap with known genetic risk factors. Unlike transcriptomes, epigenomic signatures defined by DNA methylation from surrogate tissues such as placenta and cord blood can reflect past differences in fetal brain gene transcription, transcription factor binding, and chromatin. For example, the discovery of NHIP (neuronal hypoxia inducible, placenta associated) through an epigenome-wide association in placenta, identified a common genetic risk for ASD that was modified by prenatal vitamin use. While epigenomic signatures are distinct between different genetic syndromic causes of ASD, bivalent chromatin and some convergent gene pathways are consistently epigenetically altered in both syndromic and idiopathic ASD, as well as some environmental exposures. Together, these epigenomic signatures hold promising clues towards improved early prediction and prevention of ASD as well genes and gene pathways to target for pharmacological interventions. Future advancements in single cell and multi-omic technologies, machine learning, as well as non-invasive screening of epigenomic signatures during pregnancy or newborn periods are expected to continue to impact the translatability of the recent discoveries in epigenomics to precision public health.

RUNX1 Upregulation Causes Mitochondrial Dysfunction via Regulating the PI3K-Akt Pathway in iPSC from Patients with Down Syndrome

Down syndrome (DS) is the most common autosomal aneuploidy caused by trisomy of chromosome 21. Previous studies demonstrated that DS affected mitochondrial functions, which may be associated with the abnormal development of the nervous system in patients with DS. Runt-related transcription factor 1 (RUNX1) is an encoding gene located on chromosome 21. It has been reported that RUNX1 may affect cell apoptosis via the mitochondrial pathway. The present study investigated whether RUNX1 plays a critical role in mitochondrial dysfunction in DS and explored the mechanism by which RUNX1 affects mitochondrial functions. Expression of RUNX1 was detected in induced pluripotent stem cells of patients with DS (DS-iPSCs) and normal iPSCs (N-iPSCs), and the mitochondrial functions were investigated in the current study. Subsequently, RUNX1 was overexpressed in N-iPSCs and inhibited in DS-iPSCs. The mitochondrial functions were investigated thoroughly, including reactive oxygen species levels, mitochondrial membrane potential, ATP content and lysosomal activity. Finally, RNA-sequencing was used to explore the global expression pattern. It was observed that the expression levels of RUNX1 in DS-iPSCs were significantly higher than those in normal controls. Impaired mitochondrial functions were observed in DS-iPSCs. Of note, overexpression of RUNX1 in N-iPSCs resulted in mitochondrial dysfunction, while inhibition of RUNX1 expression could improve the mitochondrial function in DS-iPSCs. Global gene expression analysis indicated that overexpression of RUNX1 may promote the induction of apoptosis in DS-iPSCs by activating the PI3K/Akt signaling pathway. The present findings indicate that abnormal expression of RUNX1 may play a critical role in mitochondrial dysfunction in DS-iPSCs.

Integrative epigenome profiling of 47XXY provides insights into whole genomic DNA hypermethylation and active chromatin accessibility

Klinefelter syndrome (KS, 47XXY) is a disorder characterized by sex chromosomal aneuploidy, which may lead to changes in epigenetic regulations of gene expression. To define epigenetic architectures in 47XXY, we annotated DNA methylation in euploid males (46XY) and females (46XX), and 47XXY individuals using whole genome bisulfite sequencing (WGBS) and integrated chromatin accessbilty, and detected abnormal hypermethylation in 47XXY. Furthermore, we detected altered chromatin accessibility in 47XXY, in particular in chromosome X, using Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) in cultured amniotic cells. Our results construct the whole genome-wide DNA methylation map in 47XXY, and provide new insights into the early epigenomic dysregulation resulting from an extra chromosome X in 47XXY.

Genetics and Molecular Basis of Congenital Heart Defects in Down Syndrome: Role of Extracellular Matrix Regulation

Down syndrome (DS), a complex disorder that is caused by the trisomy of chromosome 21 (Hsa21), is a major cause of congenital heart defects (CHD). Interestingly, only about 50% of individuals with Hsa21 trisomy manifest CHD. Here we review the genetic basis of CHD in DS, focusing on genes that regulate extracellular matrix (ECM) organization. The overexpression of Hsa21 genes likely underlies the molecular mechanisms that contribute to CHD, even though the genes responsible for CHD could only be located in a critical region of Hsa21. A role in causing CHD has been attributed not only to protein-coding Hsa21 genes, but also to genes on other chromosomes, as well as miRNAs and lncRNAs. It is likely that the contribution of more than one gene is required, and that the overexpression of Hsa21 genes acts in combination with other genetic events, such as specific mutations or polymorphisms, amplifying their effect. Moreover, a key function in determining alterations in cardiac morphogenesis might be played by ECM. A large number of genes encoding ECM proteins are overexpressed in trisomic human fetal hearts, and many of them appear to be under the control of a Hsa21 gene, the RUNX1 transcription factor.

The Genetics of Neurodevelopment in Congenital Heart Disease

Congenital heart disease (CHD) is the most common birth anomaly, affecting almost 1% of infants. Neurodevelopmental delay is the most common extracardiac feature in people with CHD. Many factors may contribute to neurodevelopmental risk, including genetic factors, CHD physiology, and the prenatal/postnatal environment. Damaging variants are most highly enriched among individuals with extracardiac anomalies or neurodevelopmental delay in addition to CHD, indicating that genetic factors have an impact beyond cardiac tissues in people with CHD. Potential sources of genetic risk include large deletions or duplications that affect multiple genes, such as 22q11 deletion syndrome, single genes that alter both heart and brain development, such as CHD7, and common variants that affect neurodevelopmental resiliency, such as APOE. Increased use of genome-sequencing technologies in studies of neurodevelopmental outcomes in people with CHD will improve our ability to detect relevant genes and variants. Ultimately, such knowledge can lead to improved and more timely intervention of learning support for affected children.

Multi-omic brain and behavioral correlates of cell-free fetal DNA methylation in macaque maternal obesity models

Maternal obesity during pregnancy is associated with neurodevelopmental disorder (NDD) risk. We utilized integrative multi-omics to examine maternal obesity effects on offspring neurodevelopment in rhesus macaques by comparison to lean controls and two interventions. Differentially methylated regions (DMRs) from longitudinal maternal blood-derived cell-free fetal DNA (cffDNA) significantly overlapped with DMRs from infant brain. The DMRs were enriched for neurodevelopmental functions, methylation-sensitive developmental transcription factor motifs, and human NDD DMRs identified from brain and placenta. Brain and cffDNA methylation levels from a large region overlapping mir-663 correlated with maternal obesity, metabolic and immune markers, and infant behavior. A DUX4 hippocampal co-methylation network correlated with maternal obesity, infant behavior, infant hippocampal lipidomic and metabolomic profiles, and maternal blood measurements of DUX4 cffDNA methylation, cytokines, and metabolites. We conclude that in this model, maternal obesity was associated with changes in the infant brain and behavior, and these differences were detectable in pregnancy through integrative analyses of cffDNA methylation with immune and metabolic factors.

Prenatal vitamin intake in first month of pregnancy and DNA methylation in cord blood and placenta in two prospective cohorts

BACKGROUND

Prenatal vitamin use is recommended before and during pregnancies for normal fetal development. Prenatal vitamins do not have a standard formulation, but many contain calcium, folic acid, iodine, iron, omega-3 fatty acids, zinc, and vitamins A, B6, B12, and D, and usually they contain higher concentrations of folic acid and iron than regular multivitamins in the US Nutrient levels can impact epigenetic factors such as DNA methylation, but relationships between maternal prenatal vitamin use and DNA methylation have been relatively understudied. We examined use of prenatal vitamins in the first month of pregnancy in relation to cord blood and placenta DNA methylation in two prospective pregnancy cohorts: the Early Autism Risk Longitudinal Investigation (EARLI) and Markers of Autism Risk Learning Early Signs (MARBLES) studies.

RESULTS

In placenta, prenatal vitamin intake was marginally associated with -0.52% (95% CI -1.04, 0.01) lower mean array-wide DNA methylation in EARLI, and associated with -0.60% (-1.08, -0.13) lower mean array-wide DNA methylation in MARBLES. There was little consistency in the associations between prenatal vitamin intake and single DNA methylation site effect estimates across cohorts and tissues, with only a few overlapping sites with correlated effect estimates. However, the single DNA methylation sites with p-value < 0.01 (EARLI cord nCpGs = 4068, EARLI placenta nCpGs = 3647, MARBLES cord nCpGs = 4068, MARBLES placenta nCpGs = 9563) were consistently enriched in neuronal developmental pathways.

CONCLUSIONS

Together, our findings suggest that prenatal vitamin intake in the first month of pregnancy may be related to lower placental global DNA methylation and related to DNA methylation in brain-related pathways in both placenta and cord blood.

Innovative in Silico Approaches for Characterization of Genes and Proteins

Bioinformatics is an amalgamation of biology, mathematics and computer science. It is a science which gathers the information from biology in terms of molecules and applies the informatic techniques to the gathered information for understanding and organizing the data in a useful manner. With the help of bioinformatics, the experimental data generated is stored in several databases available online like nucleotide database, protein databases, GENBANK and others. The data stored in these databases is used as reference for experimental evaluation and validation. Till now several online tools have been developed to analyze the genomic, transcriptomic, proteomics, epigenomics and metabolomics data. Some of them include Human Splicing Finder (HSF), Exonic Splicing Enhancer Mutation taster, and others. A number of SNPs are observed in the non-coding, intronic regions and play a role in the regulation of genes, which may or may not directly impose an effect on the protein expression. Many mutations are thought to influence the splicing mechanism by affecting the existing splice sites or creating a new sites. To predict the effect of mutation (SNP) on splicing mechanism/signal, HSF was developed. Thus, the tool is helpful in predicting the effect of mutations on splicing signals and can provide data even for better understanding of the intronic mutations that can be further validated experimentally. Additionally, rapid advancement in proteomics have steered researchers to organize the study of protein structure, function, relationships, and dynamics in space and time. Thus the effective integration of all of these technological interventions will eventually lead to steering up of next-generation systems biology, which will provide valuable biological insights in the field of research, diagnostic, therapeutic and development of personalized medicine.

Comethyl: a network-based methylome approach to investigate the multivariate nature of health and disease

Health outcomes are frequently shaped by difficult to dissect inter-relationships between biological, behavioral, social and environmental factors. DNA methylation patterns reflect such multivariate intersections, providing a rich source of novel biomarkers and insight into disease etiologies. Recent advances in whole-genome bisulfite sequencing enable investigation of DNA methylation over all genomic CpGs, but existing bioinformatic approaches lack accessible system-level tools. Here, we develop the R package Comethyl, for weighted gene correlation network analysis of user-defined genomic regions that generates modules of comethylated regions, which are then tested for correlations with multivariate sample traits. First, regions are defined by CpG genomic location or regulatory annotation and filtered based on CpG count, sequencing depth and variability. Next, correlation networks are used to find modules of interconnected nodes using methylation values within the selected regions. Each module containing multiple comethylated regions is reduced in complexity to a single eigennode value, which is then tested for correlations with experimental metadata. Comethyl has the ability to cover the noncoding regulatory regions of the genome with high relevance to interpretation of genome-wide association studies and integration with other types of epigenomic data. We demonstrate the utility of Comethyl on a dataset of male cord blood samples from newborns later diagnosed with autism spectrum disorder (ASD) versus typical development. Comethyl successfully identified an ASD-associated module containing regions mapped to genes enriched for brain glial functions. Comethyl is expected to be useful in uncovering the multivariate nature of health disparities for a variety of common disorders. Comethyl is available at github.com/cemordaunt/comethyl with complete documentation and example analyses.

Placenta and fetal brain share a neurodevelopmental disorder DNA methylation profile in a mouse model of prenatal PCB exposure

Polychlorinated biphenyls (PCBs) are developmental neurotoxicants implicated as environmental risk factors for neurodevelopmental disorders (NDDs). Here, we report the effects of prenatal exposure to a human-relevant mixture of PCBs on the DNA methylation profiles of mouse placenta and fetal brain. Thousands of differentially methylated regions (DMRs) distinguish placenta and fetal brain from PCB-exposed mice from sex-matched vehicle controls. In both placenta and fetal brain, PCB-associated DMRs are enriched for functions related to neurodevelopment and cellular signaling and enriched within regions of bivalent chromatin. The placenta and brain PCB DMRs overlap significantly and map to a shared subset of genes enriched for Wnt signaling, Slit/Robo signaling, and genes differentially expressed in NDD models. The consensus PCB DMRs also significantly overlap with DMRs from human NDD brain and placenta. These results demonstrate that PCB-exposed placenta contains a subset of DMRs that overlap fetal brain DMRs relevant to an NDD.

Placental methylome reveals a 22q13.33 brain regulatory gene locus associated with autism

BACKGROUND

Autism spectrum disorder (ASD) involves complex genetics interacting with the perinatal environment, complicating the discovery of common genetic risk. The epigenetic layer of DNA methylation shows dynamic developmental changes and molecular memory of in utero experiences, particularly in placenta, a fetal tissue discarded at birth. However, current array-based methods to identify novel ASD risk genes lack coverage of the most structurally and epigenetically variable regions of the human genome.

RESULTS

We use whole genome bisulfite sequencing in placenta samples from prospective ASD studies to discover a previously uncharacterized ASD risk gene, LOC105373085, renamed NHIP. Out of 134 differentially methylated regions associated with ASD in placental samples, a cluster at 22q13.33 corresponds to a 118-kb hypomethylated block that replicates in two additional cohorts. Within this locus, NHIP is functionally characterized as a nuclear peptide-encoding transcript with high expression in brain, and increased expression following neuronal differentiation or hypoxia, but decreased expression in ASD placenta and brain. NHIP overexpression increases cellular proliferation and alters expression of genes regulating synapses and neurogenesis, overlapping significantly with known ASD risk genes and NHIP-associated genes in ASD brain. A common structural variant disrupting the proximity of NHIP to a fetal brain enhancer is associated with NHIP expression and methylation levels and ASD risk, demonstrating a common genetic influence.

CONCLUSIONS

Together, these results identify and initially characterize a novel environmentally responsive ASD risk gene relevant to brain development in a hitherto under-characterized region of the human genome.

Update on Etiology and Pathogenesis of Biliary Atresia

Biliary atresia is a rare inflammatory sclerosing obstructive cholangiopathy that initiates in infancy as complete choledochal blockage and progresses to the involvement of intrahepatic biliary epithelium. Growing evidence shows that biliary atresia is not a single entity with a single etiology but a phenotype resulting from multifactorial events whose common path is obliterative cholangiopathy. The etiology of biliary atresia has been explained as resulting from genetic variants, toxins, viral infection, chronic inflammation or bile duct lesions mediated by autoimmunity, abnormalities in the development of the bile ducts, and defects in embryogenesis, abnormal fetal or prenatal circulation and susceptibility factors. It is increasingly evident that the genetic and epigenetic predisposition combined with the environmental factors to which the mother is exposed are potential triggers for biliary atresia. There is also an indication that a progressive thickening of the arterial middle layer occurs in this disease, suggestive of vascular remodeling and disappearance of the interlobular bile ducts. It is suggested that the hypoxia/ischemia process can affect portal structures in biliary atresia and is associated with both the extent of biliary proliferation and the thickening of the medial layer.

Long-term effects of wildfire smoke exposure during early life on the nasal epigenome in rhesus macaques

BACKGROUND

Wildfire smoke is responsible for around 20% of all particulate emissions in the U.S. and affects millions of people worldwide. Children are especially vulnerable, as ambient air pollution exposure during early childhood is associated with reduced lung function. Most studies, however, have focused on the short-term impacts of wildfire smoke exposures. We aimed to identify long-term baseline epigenetic changes associated with early-life exposure to wildfire smoke. We collected nasal epithelium samples for whole genome bisulfite sequencing (WGBS) from two groups of adult female rhesus macaques: one group born just before the 2008 California wildfire season and exposed to wildfire smoke during early-life (n = 8), and the other group born in 2009 with no wildfire smoke exposure during early-life (n = 14). RNA-sequencing was also performed on a subset of these samples.

RESULTS

We identified 3370 differentially methylated regions (DMRs) (difference in methylation ≥ 5%, empirical p < 0.05) and 1 differentially expressed gene (FLOT2) (FDR < 0.05, fold of change ≥ 1.2). The DMRs were annotated to genes significantly enriched for synaptogenesis signaling, protein kinase A signaling, and a variety of immune processes, and some DMRs significantly correlated with gene expression differences. DMRs were also significantly enriched within regions of bivalent chromatin (top odds ratio = 1.46, q-value < 3 × 10^-6) that often silence key developmental genes while keeping them poised for activation in pluripotent cells.

CONCLUSIONS

These data suggest that early-life exposure to wildfire smoke leads to long-term changes in the methylome over genes impacting the nervous and immune systems. Follow-up studies will be required to test whether these changes influence transcription following an immune/respiratory challenge.

Ten Years of EWAS

Epigenome-wide association study (EWAS) has been applied to analyze DNA methylation variation in complex diseases for a decade, and epigenome as a research target has gradually become a hot topic of current studies. The DNA methylation microarrays, next-generation, and third-generation sequencing technologies have prepared a high-quality platform for EWAS. Here, the progress of EWAS research is reviewed, its contributions to clinical applications, and mainly describe the achievements of four typical diseases. Finally, the challenges encountered by EWAS and make bold predictions for its future development are presented.

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.

Increased risk of leukaemia in children with Down syndrome: a somatic evolutionary view

Children show a higher incidence of leukaemia compared with young adolescents, yet their cells are less damaged because of their young age. Children with Down syndrome (DS) have an even higher risk of developing leukaemia during the first years of life. The presence of a constitutive trisomy of chromosome 21 (T21) in DS acts as a genetic driver for leukaemia development, however, additional oncogenic mutations are required. Therefore, T21 provides the opportunity to better understand leukaemogenesis in children. Here, we describe the increased risk of leukaemia in DS during childhood from a somatic evolutionary view. According to this idea, cancer is caused by a variation in inheritable phenotypes within cell populations that are subjected to selective forces within the tissue context. We propose a model in which the increased risk of leukaemia in DS children derives from higher rates of mutation accumulation, already present during fetal development, which is further enhanced by changes in selection dynamics within the fetal liver niche. This model could possibly be used to understand the rate-limiting steps of leukaemogenesis early in life.

The role of DNA methylation in syndromic and non-syndromic congenital heart disease

Congenital heart disease (CHD) is a common structural birth defect worldwide, and defects typically occur in the walls and valves of the heart or enlarged blood vessels. Chromosomal abnormalities and genetic mutations only account for a small portion of the pathogenic mechanisms of CHD, and the etiology of most cases remains unknown. The role of epigenetics in various diseases, including CHD, has attracted increased attention. The contributions of DNA methylation, one of the most important epigenetic modifications, to CHD have not been illuminated. Increasing evidence suggests that aberrant DNA methylation is related to CHD. Here, we briefly introduce DNA methylation and CHD and then review the DNA methylation profiles during cardiac development and in CHD, abnormalities in maternal genome-wide DNA methylation patterns are also described. Whole genome methylation profile and important differentially methylated genes identified in recent years are summarized and clustered according to the sample type and methodologies. Finally, we discuss the novel technology for and prospects of CHD-related DNA methylation.

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.