Abstract
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 × 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.
Age-related hearing loss is a complex disease caused by a combination of genetic and environmental factors, and a study have conducted animal experiments to explore the association between BCL11B ...heterozygosity and age-related hearing loss. The present study used established genetic models to examine the association between BCL11B gene polymorphisms and age-related hearing loss. A total of 410 older adults from two communities in Qingdao, China, participated in this study. The case group comprised individuals aged greater than or equal to 60 years with age-related hearing loss, and the control group comprised individuals without age-related hearing loss from the same communities. The groups were matched 1:1 for age and sex. The individual characteristics of the participants were analyzed descriptively using the Mann-Whitney U test and the chi-square test. To explore the association between BCL11B gene polymorphisms and age-related hearing loss, conditional logistic regression was performed to construct genetic models for two single-nucleotide-polymorphisms (SNPs) of BCL11B, and haplotype analysis was conducted to construct their haplotype domains. Two SNP sites of the BCL11B gene, four genetic models of rs1152781 (additive, dominant, recessive, and codominant), and five genetic models of rs1152783 (additive, dominant, recessive, codominant, and over dominant) were significantly associated with age-related hearing loss in the models both unadjusted and adjusted for all covariates (P < 0.05). Additionally, a linkage disequilibrium between rs1152781 and rs1152783 was revealed through haplotype analysis. Our study revealed that BCL11B gene polymorphisms were significantly associated with age-related hearing loss.
Smoke-free government(SFG), as a key tobacco control measure, has been added in Healthy China 2030 blueprint and Qingdao started the establishment of the demonstrative SFG in 2020.This study examined ...the effects of SFG policy on smoking and smoke-free(SF) environment after establishing the demonstrative SFG. This cross-sectional survey selected participants by simple random sampling from party and government agencies in Qingdao (N = 3625) and the participants filled in questionnaires online from November 31 to December 15, 2020. We utilized AMOS to set up models to analyze the direct and indirect effects of SFG policy. The findings showed that knowledge of SFG policy was positively associated with SF environment(β = 0.29, P<0.001) and negatively associated with smoking(β = -0.14,P<0.001). Knowledge of SFG policy had indirect effects on SF environment through four channels: independent mediation of discouraging smoking and attitude towards SFG policy, as well as chain mediation of these two factors, and perception of tobacco hazards and discouraging smoking, with indirect effects accounting for 33.5% of the total effect. Knowledge of SFG policy had indirect effects on smoking reduction via SF environment and two chain mediation: SF environment and attitude towards SFG policy, perception of tobacco hazards and intention to quit smoking, with indirect effects accounting for 50.2% of the total effect. The results provided the evidence that SFG policy not only had positive effects on creating SF environment but also on reducing smoking. The efficient policy infiltration to individuals played a vital role in the establishment of SFG. Attitude towards SFG policy, discouraging smoking and SF environment were the potential mediators for SFG policy. Findings in this study added more evidence related to effect mechanism of SFG policy and had a positive influence on promoting the implementation SFG policies for China and other countries.
Grip strength is an important biomarker reflecting muscle strength, and depression is a psychiatric disorder all over the world. Several studies found a significant inverse association between grip ...strength and depression, and there is also evidence for common physiological mechanisms between them. We used twin data from Qingdao, China to calculate genetic correlations, and we performed a bivariate GWAS to explore potential SNPs, genes, and pathways in common between grip strength and depression. 139 pairs of Dizygotic twins were used for bivariate GWAS. VEAGSE2 and PASCAL software were used for gene-based analysis and pathway enrichment analysis, respectively. And the resulting SNPs were subjected to eQTL analysis and pleiotropy analysis. The genetic correlation coefficient between grip strength and depression was -0.41 (-0.96, -0.15). In SNP-based analysis, 7 SNPs exceeded the genome-wide significance level (P<5×10-8) and a total of 336 SNPs reached the level of suggestive significance (P<1×10-5). Gene-based analysis and pathway-based analysis identified genes and pathways related to muscle strength and the nervous system. The results of eQTL analysis were mainly enriched in tissues such as the brain, thyroid, and skeletal muscle. Pleiotropy analysis shows that 9 of the 15 top SNPs were associated with both grip strength and depression. In conclusion, this bivariate GWAS identified potentially common pleiotropic SNPs, genes, and pathways in grip strength and depression.
Recently, new loci related to body mass index (BMI) or blood pressure (BP) have been identified respectively in genome-wide association studies (GWAS). However, limited studies focused on jointly ...associated genetic variance between systolic pressure (SBP), diastolic pressure (DBP) and BMI. Therefore, a bivariate twin study was performed to explore the genetic variants associated with BMI-SBP, BMI-DBP and SBP-DBP. A total of 380 twin pairs (137 dizygotic pairs and 243 monozygotic pairs) recruited from Qingdao Twin Registry system were used to access the genetic correlations (0.2108 for BMI-SBP, 0.2345 for BMI-DBP, and 0.6942 for SBP-DBP, respectively) by bivariate Cholesky decomposition model. Bivariate GWAS in 137 dizygotic pairs nominated 27 single identified 27 quantitative trait nucleotides (QTNs) for BMI and SBP, 27 QTNs for BMI and DBP, and 25 QTNs for SBP and DBP with the suggestive P-value threshold of 1×10-5. After imputation, we found eight SNPs, one for both BMI-SBP and SBP-DBP, and eight for SBP-DBP, exceed significant statistic level. Expression quantitative trait loci analysis identified rs4794029 as new significant eQTL in tissues related to BMI and SBP. Also, we found 6 new significant eQTLs (rs4400367, rs10113750, rs11776003, rs3739327, rs55978930, and rs4794029) in tissues were related to SBP and DBP. Gene-based analysis identified nominally associated genes (P < 0.05) with BMI-SBP, BMI-DBP, and SBP-DBP, respectively, such as PHOSPHO1, GNGT2, KEAP1, and S1PR5. In the pathway analysis, we found some pathways associated with BMI-SBP, BMI-DBP and SBP-DBP, such as prion diseases, IL5 pathway, cyclin E associated events during G1/S transition, TGF beta signaling pathway, G βγ signaling through PI3Kγ, prolactin receptor signaling etc. These findings may enrich the results of genetic variants related to BMI and BP traits, and provide some evidences to future study the pathogenesis of hypertension and obesity in the northern Chinese population.
Currently, diabetes has become one of the leading causes of death worldwide. Fasting plasma glucose (FPG) levels that are higher than optimal, even if below the diagnostic threshold of diabetes, can ...also lead to increased morbidity and mortality. Here we intend to study the magnitude of the genetic influence on FPG variation by conducting structural equation modelling analysis and to further identify specific genetic variants potentially related to FPG levels by performing a genome-wide association study (GWAS) in Chinese twins.
The final sample included 382 twin pairs: 139 dizygotic (DZ) pairs and 243 monozygotic (MZ) pairs. The DZ twin correlation for the FPG level (r
= 0.20, 95% CI: 0.04-0.36) was much lower than half that of the MZ twin correlation (r
= 0.68, 95% CI: 0.62-0.74). For the variation in FPG level, the AE model was the better fitting model, with additive genetic parameters (A) accounting for 67.66% (95% CI: 60.50-73.62%) and unique environmental or residual parameters (E) accounting for 32.34% (95% CI: 26.38-39.55%), respectively. In the GWAS, although no genetic variants reached the genome-wide significance level (P < 5 × 10
), 28 SNPs exceeded the level of a suggestive association (P < 1 × 10
). One promising genetic region (2q33.1) around rs10931893 (P = 1.53 × 10
) was found. After imputing untyped SNPs, we found that rs60106404 (P = 2.38 × 10
) located at SPATS2L reached the genome-wide significance level, and 216 SNPs exceeded the level of a suggestive association. We found 1007 genes nominally associated with the FPG level (P < 0.05), including SPATS2L, KCNK5, ADCY5, PCSK1, PTPRA, and SLC26A11. Moreover, C1orf74 (P = 0.014) and SLC26A11 (P = 0.021) were differentially expressed between patients with impaired fasting glucose and healthy controls. Some important enriched biological pathways, such as β-alanine metabolism, regulation of insulin secretion, glucagon signaling in metabolic regulation, IL-1 receptor pathway, signaling by platelet derived growth factor, cysteine and methionine metabolism pathway, were identified.
The FPG level is highly heritable in the Chinese population, and genetic variants are significantly involved in regulatory domains, functional genes and biological pathways that mediate FPG levels. This study provides important clues for further elucidating the molecular mechanism of glucose homeostasis and discovering new diagnostic biomarkers and therapeutic targets for diabetes.
Previous studies have determined the epigenetic association between DNA methylation and pulmonary function among various ethnics, whereas this association is largely unknown in Chinese adults. Thus, ...we aimed to explore epigenetic relationships between genome-wide DNA methylation levels and pulmonary function among middle-aged Chinese monozygotic twins.
The monozygotic twin sample was drawn from the Qingdao Twin Registry. Pulmonary function was measured by three parameters including forced expiratory volume the first second (FEV1), forced vital capacity (FVC), and FEV1/FVC ratio. Linear mixed effect model was used to regress the methylation level of CpG sites on pulmonary function. After that, we applied Genomic Regions Enrichment of Annotations Tool (GREAT) to predict the genomic regions enrichment, and used comb-p python library to detect differentially methylated regions (DMRs). Gene expression analysis was conducted to validate the results of differentially methylated analyses.
We identified 112 CpG sites with the level of P < 1 × 10
which were annotated to 40 genes. We identified 12 common enriched pathways of three pulmonary function parameters. We detected 39 DMRs located at 23 genes, of which PRDM1 was related to decreased pulmonary function, and MPL, LTB4R2, and EPHB3 were related to increased pulmonary function. The gene expression analyses validated DIP2C, ASB2, SLC6A5, and GAS6 related to decreased pulmonary function.
Our DNA methylation sequencing analysis on identical twins provides new references for the epigenetic regulation on pulmonary function. Several CpG sites, genes, biological pathways and DMRs are considered as possible crucial to pulmonary function.
The effects of fine particulate matter (PM2.5) on acute myocardial infarction (AMI) have been widely recognized. However, no studies have comprehensively evaluated future PM2.5-attributed AMI burdens ...under different climate mitigation and population change scenarios. We aimed to quantify the PM2.5-AMI association and estimate the future change in PM2.5-attributed AMI incident cases under six integrated scenarios in 2030 and 2060 in Shandong Province, China.
Daily AMI incident cases and air pollutant data were collected from 136 districts/counties in Shandong Province from 2017 − 2019. A two-stage analysis with a distributed lag nonlinear model was conducted to quantify the baseline PM2.5-AMI association. The future change in PM2.5-attributed AMI incident cases was estimated by combining the fitted PM2.5-AMI association with the projected daily PM2.5 concentrations under six integrated scenarios. We further analyzed the factors driving changes in PM2.5-related AMI incidence using a decomposition method.
Each 10 μg/m3 increase in PM2.5 exposure at lag05 was related to an excess risk of 1.3 % (95 % confidence intervals: 0.9 %, 1.7 %) for AMI incidence from 2017 − 2019 in Shandong Province. The estimated total PM2.5-attributed AMI incident cases would increase by 10.9−125.9 % and 6.4–244.6 % under Scenarios 1 − 3 in 2030 and 2060, whereas they would decrease by 0.9–5.2 % and 33.0–46.2 % under Scenarios 5 – 6 in 2030 and 2060, respectively. Furthermore, the percentage increases in PM2.5-attributed female cases (2030: −0.3 % to 135.1 %; 2060: −33.2 % to 321.5 %) and aging cases (2030: 15.2–171.8 %; 2060: −21.5 % to 394.2 %) would wholly exceed those in male cases (2030: −1.8 % to 133.2 %; 2060: −41.1 % to 264.3 %) and non-aging cases (2030: −41.0 % to 45.7 %; 2060: −89.5 % to −17.0 %) under six scenarios in 2030 and 2060. Population aging is the main driver of increased PM2.5-related AMI incidence under Scenarios 1 − 3 in 2030 and 2060, while improved air quality can offset these negative effects of population aging under the implementation of the carbon neutrality and 1.5 °C targets.
The combination of ambitious climate policies (i.e., 1.5 °C warming limits and carbon neutrality targets) with stringent clean air policies is necessary to reduce the health impacts of air pollution in Shandong Province, China, regardless of population aging.
•Short-term exposure to PM2.5 might be related to increased risks of AMI incidence.•PM2.5-attributed AMI incident cases would be varied under six integrated scenarios.•Women and the elderly would suffer from increased burden of PM2.5-related AMI incident.•More coordinated policy efforts are needed for the effective disease burden reduction.
The therapeutic management of obesity is challenging, hence further elucidating the underlying mechanisms of obesity development and identifying new diagnostic biomarkers and therapeutic targets are ...urgent and necessary. Here, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) to identify significant genes and specific modules related to BMI based on gene expression profile data of 7 discordant monozygotic twins.
In the differential gene expression analysis, it appeared that 32 differentially expressed genes (DEGs) were with a trend of up-regulation in twins with higher BMI when compared to their siblings. Categories of positive regulation of nitric-oxide synthase biosynthetic process, positive regulation of NF-kappa B import into nucleus, and peroxidase activity were significantly enriched within GO database and NF-kappa B signaling pathway within KEGG database. DEGs of NAMPT, TLR9, PTGS2, HBD, and PCSK1N might be associated with obesity. In the WGCNA, among the total 20 distinct co-expression modules identified, coral1 module (68 genes) had the strongest positive correlation with BMI (r = 0.56, P = 0.04) and disease status (r = 0.56, P = 0.04). Categories of positive regulation of phospholipase activity, high-density lipoprotein particle clearance, chylomicron remnant clearance, reverse cholesterol transport, intermediate-density lipoprotein particle, chylomicron, low-density lipoprotein particle, very-low-density lipoprotein particle, voltage-gated potassium channel complex, cholesterol transporter activity, and neuropeptide hormone activity were significantly enriched within GO database for this module. And alcoholism and cell adhesion molecules pathways were significantly enriched within KEGG database. Several hub genes, such as GAL, ASB9, NPPB, TBX2, IL17C, APOE, ABCG4, and APOC2 were also identified. The module eigengene of saddlebrown module (212 genes) was also significantly correlated with BMI (r = 0.56, P = 0.04), and hub genes of KCNN1 and AQP10 were differentially expressed.
We identified significant genes and specific modules potentially related to BMI based on the gene expression profile data of monozygotic twins. The findings may help further elucidate the underlying mechanisms of obesity development and provide novel insights to research potential gene biomarkers and signaling pathways for obesity treatment. Further analysis and validation of the findings reported here are important and necessary when more sample size is acquired.
The association of calcium signaling pathway gene variants, bone mineral density (BMD) and mild cognitive impairment (MCI) is poorly understood so far. A total of 878 participants from Qingdao city ...were recruited in this study. According to the candidate gene selection method, 58 single nucleotide polymorphisms (SNPs) in eight calcium signaling genes were selected. The association between gene polymorphisms and MCI was revealed by using multiple genetic models. Polygenic risk scores (PRS) were used to summarize the effects of the whole gene. Logistic regression was used to analyze the association between each PRS and MCI. The multiplicative interaction term in the regression models was used to estimate the interaction effects between the PRS and BMD. We observed significant associations of rs6877893 (NR3C1), rs6448456 (CCKAR), and rs723672 (CACNA1C) polymorphisms with MCI. The PRSs of NR3C1 (OR = 4.012, 95% CI = 1.722-9.347,
< 0.001), PRKCA (OR = 1.414, 95% CI = 1.083-1.845,
= 0.011) and TRPM1 (OR = 3.253, 95% CI = 1.116-9.484,
= 0.031) were associated with an increased risk of developing MCI, and the PRS of total genes (OR = 0.330, 95% CI = 0.224-0.485,
< 0.001) was associated with a decreased risk of developing MCI. In interaction effect analysis, the interaction effect of PRKCA and BMD was significant. Genetic variations of the calcium signaling pathway were associated with MCI in older people. There was an interaction effect between PRKCA gene variants and BMD on MCI.