Genome-wide DNA methylation and gene expression analyses in monozygotic twins identify potential biomarkers of depression

Depression is currently the leading cause of disability around the world. We conducted an epigenome-wide association study (EWAS) in a sample of 58 depression score-discordant monozygotic twin pairs, aiming to detect specific epigenetic variants potentially related to depression and further integrate with gene expression profile data. Association between the methylation level of each CpG site and depression score was tested by applying a linear mixed effect model. Weighted gene co-expression network analysis (WGCNA) was performed for gene expression data. The association of DNA methylation levels of 66 CpG sites with depression score reached the level of P < 1 x 10(-4). These top CpG sites were located at 34 genes, especially PTPRN2, HES5, GATA2, PRDM7, and KCNIP1. Many ontology enrichments were highlighted, including Notch signaling pathway, Huntington disease, p53 pathway by glucose deprivation, hedgehog signaling pathway, DNA binding, and nucleic acid metabolic process. We detected 19 differentially methylated regions (DMRs), some of which were located at GRIK2, DGKA, and NIPA2. While integrating with gene expression data, HELZ2, PTPRN2, GATA2, and ZNF624 were differentially expressed. In WGCNA, one specific module was positively correlated with depression score (r = 0.62, P = 0.002). Some common genes (including BMP2, PRDM7, KCNIP1, and GRIK2) and enrichment terms (including complement and coagulation cascades pathway, DNA binding, neuron fate specification, glial cell differentiation, and thyroid gland development) were both identified in methylation analysis and WGCNA. Our study identifies specific epigenetic variations which are significantly involved in regions, functional genes, biological function, and pathways that mediate depression disorder.Peer reviewe

Differential regulation of the DNA methylome in adults born during the Great Chinese Famine in 1959-1961

Background: Extensive epidemiological studies have established the association between exposure to early-life adversity and health status and diseases in adults. Epigenetic regulation is considered as a key mediator for this phenomenon but analysis on humans is sparse. The Great Chinese Famine lasting from 1958 to 1961 is a natural string of disasters offering a precious opportunity for elucidating the underlying epigenetic mechanism of the long-term effect of early adversity. Methods: Using a high-throughput array platform for DNA methylome profiling, we conducted a case-control epigenome-wide association study on early-life exposure to Chinese famine in 79 adults born during 1959-1961 and compared to 105 unexposed subjects born 1963-1964. Results: The single CpG site analysis of whole epigenome revealed a predominant pattern of decreased DNA methylation levels associated with fetal exposure to famine. Four CpG sites were detected with p < 1e-06 (linked to EHMT1, CNR1, UBXN7 and ESM1 genes), 16 CpGs detected with 1e-06 < p < 1e-05 and 157 CpGs with 1e-05 < p < 1e-04, with a predominant pattern of hypomethylation. Functional annotation to genes and their enriched biological pathways mainly involved neurodevelopment, neuropsychological disorders and metabolism. Multiple sites analysis detected two top-rank differentially methylated regions harboring RNF39 on chromosome 6 and PTPRN2 on chromosome 7, both showing epigenetic association with stress-related conditions. Conclusion: Early-life exposure to famine could mediate DNA methylation regulations that persist into adulthood with broad impacts in the activities of genes and biological pathways. Results from this study provide new clues to the epigenetic embedding of early-life adversity and its impacts on adult health.Peer reviewe

Global Gene Expression Profiling of Body-Mass Index in Middle-Aged Danish Twins

Objective: The body mass index (BMI) measured as weight in relation to height is an important monitor for obesity and diabetes, with individual variation under control by genetic and environmental factors. In transcriptome-wide association studies on BMI, the genetic contribution calls for controlling of genetic confounding that affects both BMI and gene expression. We performed a global gene expression profiling of BMI on peripheral blood cells using monozygotic twins for efficient handling of genetic make-ups. Methods: We applied a generalized association method to genome-wide gene expression data on 229 pairs of monozygotic twins (age 56-80 years) for detecting diverse patterns of correlation between BMI and gene expression. Results: We detected seven probes associated with BMI with p<1e-04, among them two probes with p<1e-05 (p=2.83e-06 AAK1; p=7.83e-06 LILRA3). In total, the analysis found 1579 probes with nominal p<0.05. Biological pathway analysis of enriched pathways found 50 KEGG and 45 Reactome pathways (FDR<0.05). The identified top functional pathways included immune function, JAK-STAT signalling, insulin signalling and regulation of energy metabolism. Conclusion: This transcriptome-wide association study using monozygotic twins and generalized correlation identified differentially expressed genes and a broad spectrum of enriched biological pathways that may implicate metabolic health

Spatiotemporal heterogeneity and impact factors of hepatitis B and C in China from 2010 to 2018: Bayesian space–time hierarchy model

IntroductionViral hepatitis is a global public health problem, and China still faces great challenges to achieve the WHO goal of eliminating hepatitis.MethodsThis study focused on hepatitis B and C, aiming to explore the long-term spatiotemporal heterogeneity of hepatitis B and C incidence in China from 2010 to 2018 and quantify the impact of socioeconomic factors on their risk through Bayesian spatiotemporal hierarchical model.ResultsThe results showed that the risk of hepatitis B and C had significant spatial and temporal heterogeneity. The risk of hepatitis B showed a slow downward trend, and the high-risk provinces were mainly distributed in the southeast and northwest regions, while the risk of hepatitis C had a clear growth trend, and the high-risk provinces were mainly distributed in the northern region. In addition, for hepatitis B, illiteracy and hepatitis C prevalence were the main contributing factors, while GDP per capita, illiteracy rate and hepatitis B prevalence were the main contributing factors to hepatitis C.DisussionThis study analyzed the spatial and temporal heterogeneity of hepatitis B and C and their contributing factors, which can serve as a basis for monitoring efforts. Meanwhile, the data provided by this study will contribute to the effective allocation of resources to eliminate viral hepatitis and the design of interventions at the provincial level

Novel DNA methylation marker discovery by assumption-free genome-wide association analysis of cognitive function in twins

Privileged by rapid increase in available epigenomic data, epigenome-wide association studies (EWAS) are to make a profound contribution to understand the molecular mechanism of DNA methylation in cognitive aging. Current statistical methods used in EWAS are dominated by models based on multiple assumptions, for example, linear relationship between molecular profiles and phenotype, normal distribution for the methylation data and phenotype. In this study, we applied an assumption-free method, the generalized correlation coefficient (GCC), and compare it to linear models, namely the linear mixed model and kinship model. We use DNA methylation associated with a cognitive score in 400 and 206 twins as discovery and replication samples respectively. DNA methylation associated with cognitive function using GCC, linear mixed model, and kinship model, identified 65 CpGs (p < 1e-04) from discovery sample displaying both nonlinear and linear correlations. Replication analysis successfully replicated 9 of these top CpGs. When combining results of GCC and linear models to cover diverse patterns of relationships, we identified genes like KLHDC4, PAPSS2, and MRPS18B as well as pathways including focal adhesion, axon guidance, and some neurological signaling. Genomic region-based analysis found 15 methylated regions harboring 11 genes, with three verified in gene expression analysis, also the 11 genes were related to top functional clusters including neurohypophyseal hormone and maternal aggressive behaviors. The GCC approach detects valuable methylation sites missed by traditional linear models. A combination of methylation markers from GCC and linear models enriched biological pathways sensible in neurological function that could implicate cognitive performance and cognitive aging.Peer reviewe

Global Gene Expression Profiling and Transcription Factor Network Analysis of Cognitive Aging in Monozygotic Twins

Cognitive aging is one of the major problems worldwide, especially as people get older. This study aimed to perform global gene expression profiling of cognitive function to identify associated genes and pathways and a novel transcriptional regulatory network analysis to identify important regulons. We performed single transcript analysis on 400 monozygotic twins using an assumption-free generalized correlation coefficient (GCC), linear mixed-effect model (LME) and kinship model and identified six probes (one significant at the standard FDR < 0.05 while the other results were suggestive with 0.18 ≤ FDR ≤ 0.28). We combined the GCC and linear model results to cover diverse patterns of relationships, and meaningful and novel genes like APOBEC3G, H6PD, SLC45A1, GRIN3B, and PDE4D were detected. Our exploratory study showed the downregulation of all these genes with increasing cognitive function or vice versa except the SLC45A1 gene, which was upregulated with increasing cognitive function. Linear models found only H6PD and SLC45A1, the other genes were captured by GCC. Significant functional pathways (FDR < 3.95e-10) such as focal adhesion, ribosome, cysteine and methionine metabolism, Huntington's disease, eukaryotic translation elongation, nervous system development, influenza infection, metabolism of RNA, and cell cycle were identified. A total of five regulons (FDR< 1.3e-4) were enriched in a transcriptional regulatory analysis in which CTCF and REST were activated and SP3, SRF, and XBP1 were repressed regulons. The genome-wide transcription analysis using both assumption-free GCC and linear models identified important genes and biological pathways implicated in cognitive performance, cognitive aging, and neurological diseases. Also, the regulatory network analysis revealed significant activated and repressed regulons on cognitive function.Peer reviewe

Genetic and environmental determinants of O6-methylguanine DNA-methyltransferase (MGMT) gene methylation: a 10-year longitudinal study of Danish twins

Background: Epigenetic inactivation of O6-methylguanine DNA-methyltransferase (MGMT) is associated with increased sensitivity to alkylating chemotherapeutic agents in glioblastoma patients. The genetic background underlying MGMT gene methylation may explain individual differences in treatment response and provide a clue to a personalized treatment strategy. Making use of the longitudinal twin design, we aimed, for the first time, to estimate the genetic contributions to MGMT methylation in a Danish twin cohort. Methods: DNA-methylation from whole blood (18 monozygotic (MZ) and 25 dizygotic (DZ) twin pairs) repeated 10 years apart from the Longitudinal Study of Aging Danish Twins (LSADT) were used to search for genetic and environmental contributions to DNA-methylation at 170 CpG sites of across the MGMT gene. Both univariate and bivariate twin models were applied. The intraclass correlations, performed on cross-sectional data (246 MZ twin pairs) from an independent study population, the Middle-Aged Danish Twins (MADT), were used to assess the genetic influence at each CpG site of MGMT for replication. Results: Univariate twin model revealed twelve CpG sites showing significantly high heritability at intake (wave 1, h2 > 0.43), and seven CpG sites with significant heritability estimates at end of follow-up (wave 2, h2 > 0.5). There were six significant CpG sites, located at the gene body region, that overlapped among the two waves (h2 > 0.5), of which five remained significant in the bivariate twin model, which was applied to both waves. Within MZ pair correlation in these six CpGs from MADT demarks top level of genetic influence. There were 11 CpGs constantly have substantial common environmental component over the 10 years. Conclusions: We have identified 6 CpG sites linked to the MGMT gene with strong and persistent genetic control based on their DNA methylation levels. The genetic basis of MGMT gene methylation could help to explain individual differences in glioblastoma treatment response and most importantly, provide references for mapping the methylation Quantitative Trait Loci (meQTL) underlying the genetic regulation.Peer reviewe

Huobahuagen tablet improves renal function in diabetic kidney disease: a real-world retrospective cohort study

ObjectiveWe aimed to explore the value of Huobahuagen tablet (HBT) in improving decreased renal function for patients with diabetic kidney disease (DKD) over time.MethodsThis was a single-center, retrospective, real-world study on eligible 122 DKD patients who continued to use HBT + Huangkui capsule (HKC) therapy or HKC therapy without interruption or alteration in Jiangsu Province Hospital of Chinese Medicine from July 2016 to March 2022. The primary observation outcomes included estimated glomerular filtration rate (eGFR) at baseline and 1-, 3-, 6-, 9-, and 12-month follow-up visits and changes in eGFR from baseline (ΔeGFR). Propensity score (PS) and inverse probability treatment weighting (IPTW) were used to control for confounders.ResultseGFR was significantly higher in the HBT + HKC group than in the HKC alone group at the 6-, 9-, and 12-month follow-up visits (p = 0.0448, 0.0002, and 0.0037, respectively), indicating the superiority of HBT + HKC over HBT alone. Furthermore, the ΔeGFR of the HBT + HKC group was significantly higher than that of the HKC alone group at the 6- and 12-month follow-up visits (p = 0.0369 and 0.0267, respectively). In the DKD G4 patients, eGFR was higher in the HBT + HKC group at the 1-, 3-, 6-, 9-, and 12-month follow-up visits compared with baseline, with statistically significant differences at the 1-, 3-, and 6- month follow-up visits (p = 0.0256, 0.0069, and 0.0252, respectively). The fluctuations in ΔeGFR ranged from 2.54 ± 4.34 to 5.01 ± 5.55 ml/min/1.73 m2. Change in the urinary albumin/creatinine ratio from baseline did not exhibit a significant difference between the two groups at any of the follow-up visits (p &gt; 0.05 for all). Adverse event incidence was low in both groups.ConclusionThe findings of this study based on real-world clinical practice indicate that HBT + HKC therapy exhibited better efficacy in improving and protecting renal function with a favorable safety profile than HKC therapy alone. However, further large-scale prospective randomized controlled trials are warranted to confirm these results