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.
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.
Abstract
When female host feeding parasitoids encounter a potential host, they face a complicated trade-off between either laying an egg for investing in current reproduction or feeding on or killing ...the host for future reproduction. Few studies have measured these behavioral shift patterns in a given host-parasitoid association thus far. We systematically assessed the behavioral shifts and life history traits of a host feeding parasitoid,
Necremnus tutae,
on different instars of its host
Tuta absoluta. N. tutae
females, as idiobiont host feeding parasitoids, can act on the 1st–4th instar larvae of
T. absoluta
by either host feeding, parasitizing or host killing. Moreover, a significant behavioral shift was observed on different instar hosts.
N. tutae
preferred to feed on the young hosts (1st and 2nd instars), lay eggs on middle-aged hosts (3rd instars) and kill old hosts (4th instars) by ovipositor-mediated stinging. The offspring of
N. tutae
showed a significant female-biased sex ratio, with the number of instars of
T. absoluta
larvae that were parasitized increasing. Specifically, nonreproductive host mortality induced by host feeding and host killing accounted for high percentages of the total mortality (ranging from 70% on 3rd instar hosts to 88% on 1st instar and 4th instar hosts). We hypothesize that
N. tutae
could be not merely a parasitoid but also a predator. Our results shed light on the nonreproductive abilities of a host feeding parasitoid that should be given further attention, especially when evaluating the efficacy of parasitoids.
To date, little is known about the pleiotropic genetic variants among depression, cognition, and memory. The current research aimed to identify the potential pleiotropic single nucleotide ...polymorphisms (SNPs), genes, and pathways of the three phenotypes by conducting a multivariate genome-wide association study and an additional pleiotropy analysis among Chinese individuals and further validate the top variants in the UK Biobank (UKB). In the discovery phase, the participants were 139 pairs of dizygotic twins from the Qingdao Twins Registry. The genome-wide efficient mixed-model analysis identified 164 SNPs reaching suggestive significance (P < 1 × 10
). Among them, rs3967317 (P = 1.21 × 10
) exceeded the genome-wide significance level (P < 5 × 10
) and was also demonstrated to be associated with depression and memory in pleiotropy analysis, followed by rs9863698, rs3967316, and rs9261381 (P = 7.80 × 10
-5.68 × 10
), which were associated with all three phenotypes. After imputation, a total of 457 SNPs reached suggestive significance. The top SNP chr6:24597173 was located in the KIAA0319 gene, which had biased expression in brain tissues. Genes and pathways related to metabolism, immunity, and neuronal systems demonstrated nominal significance (P < 0.05) in gene-based and pathway enrichment analyses. In the validation phase, 12 of the abovementioned SNPs reached the nominal significance level (P < 0.05) in the UKB. Among them, three SNPs were located in the KIAA0319 gene, and four SNPs were identified as significant expression quantitative trait loci in brain tissues. These findings may provide evidence for pleiotropic variants among depression, cognition, and memory and clues for further exploring the shared genetic pathogenesis of depression with Alzheimer's disease.
Elevated fasting plasma glucose (FPG) levels can increase morbidity and mortality even when it is below the diagnostic threshold of type 2 diabetes mellitus (T2DM). We conducted a genome-wide DNA ...methylation analysis to detect DNA methylation (DNAm) variants potentially related to FPG in Chinese monozygotic twins.
Genome-wide DNA methylation profiling in whole blood of twins was performed using Reduced Representation Bisulfite Sequencing (RRBS), yielding 551,447 raw CpGs. Association between DNAm of single CpG and FPG was tested using a generalized estimation equation. Differentially methylated regions (DMRs) were identified using comb-P approach. ICE FALCON method was utilized to perform the causal inference. Candidate CpGs were quantified and validated using Sequenom MassARRAY platform in a community population. Weighted gene co-expression network analysis (WGCNA) was conducted using gene expression data from twins.
The mean age of 52 twin pairs was 52 years (SD: 7). The relationship between DNAm of 142 CpGs and FPG reached the genome-wide significance level. Thirty-two DMRs within 24 genes were identified, including TLCD1, MRPS31P5, CASZ1, and CXADRP3. The causal relationship of top CpGs mapped to TLCD1, MZF1, PTPRN2, SLC6A18, ASTN2, IQCA1, GRIN1, and PDE2A genes with FPG were further identified using ICE FALCON method. Pathways potentially related to FPG were also identified, such as phospholipid-hydroperoxide glutathione peroxidase activity and mitogen-activated protein kinase p38 binding. Three CpGs mapped to SLC6A18 gene were validated in a community population, with a hypermethylated direction in diabetic patients. The expression levels of 18 genes (including SLC6A18 and TLCD1) were positively correlated with FPG levels.
We detect many DNAm variants that may be associated with FPG in whole blood, particularly the loci within SLC6A18 gene. Our findings provide important reference for the epigenetic regulation of elevated FPG levels and diabetes.
Hypertension is a crucial risk factor for developing cardiovascular disease and reducing life expectancy. We aimed to detect DNA methylation (DNAm) variants potentially related to systolic blood ...pressure (SBP) and diastolic blood pressure (DBP) by conducting epigenome-wide association studies in 60 and 59 Chinese monozygotic twin pairs, respectively.
Genome-wide DNA methylation profiling in whole blood of twins was performed using Reduced Representation Bisulfite Sequencing, yielding 551,447 raw CpGs. Association between DNAm of single CpG and blood pressure was tested by applying generalized estimation equation. Differentially methylated regions (DMRs) were identified by comb-P approach. Inference about Causation through Examination of Familial Confounding was utilized to perform the causal inference. Ontology enrichment analysis was performed using Genomic Regions Enrichment of Annotations Tool. Candidate CpGs were quantified using Sequenom MassARRAY platform in a community population. Weighted gene co-expression network analysis (WGCNA) was conducted using gene expression data.
The median age of twins was 52 years (95% range 40, 66). For SBP, 31 top CpGs (p < 1 × 10
) and 8 DMRs were identified, with several DMRs within NFATC1, CADM2, IRX1, COL5A1, and LRAT. For DBP, 43 top CpGs (p < 1 × 10
) and 12 DMRs were identified, with several DMRs within WNT3A, CNOT10, and DAB2IP. Important pathways, such as Notch signaling pathway, p53 pathway by glucose deprivation, and Wnt signaling pathway, were significantly enriched for SBP and DBP. Causal inference analysis suggested that DNAm at top CpGs within NDE1, MYH11, SRRM1P2, and SMPD4 influenced SBP, while SBP influenced DNAm at CpGs within TNK2. DNAm at top CpGs within WNT3A influenced DBP, while DBP influenced DNAm at CpGs within GNA14. Three CpGs mapped to WNT3A and one CpG mapped to COL5A1 were validated in a community population, with a hypermethylated and hypomethylated direction in hypertension cases, respectively. Gene expression analysis by WGCNA further identified some common genes and enrichment terms.
We detect many DNAm variants that may be associated with blood pressure in whole blood, particularly the loci within WNT3A and COL5A1. Our findings provide new clues to the epigenetic modification underlying hypertension pathogenesis.
Dyslipidemia represents a strong and independent risk factor for cardiovascular disease. Plasma cholesterol, such as total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), and high ...density lipoprotein cholesterol (HDL-C), is the common indicator of diagnosing dyslipidemia. Here based on 382 Chinese twin pairs, we explored the magnitude of genetic impact on TC, HDL-C, LDL-C variation and further searched for genetic susceptibility loci for them using genome-wide association study (GWAS). The ACE model was the best fit model with additive genetic parameter (A) accounting for 26.6%, common or shared environmental parameter (C) accounting for 47.8%, unique/non-shared environmental parameter (E) accounting for 25.6% for the variance in HDL-C. The AE model was the best fit model for TC (A: 61.4%; E: 38.6%) and LDL-C (A: 65.5%; E: 34.5%). While no SNPs reached the genome-wide significance level (
< 5 × 10
), 8, 14, 9 SNPs exceeded the suggestive significance level (
< 1 × 10
) for TC, HDL-C, LDL-C, respectively. The promising genetic regions for TC, HDL-C, LDL-C were on chromosome 11 around rs7107698, chromosome 5 around rs12518218, chromosome 2 around rs10490120, respectively. Gene-based analysis found 1038, 1033 and 1090 genes nominally associated with TC, HDL-C, LDL-C (
< 0.05), especially
for TC,
for HDL-C, and
for LDL-C, respectively. The number of common related genes among TC, HDL-C and LDL-C was 71, including
, etc. Pathway enrichment analysis discovered known related pathways-zinc transporters, metal ion SLC transporters for TC, cell adhesion molecules CAMs, IL-6 signaling for HDL, FC epsilon RI signaling pathway, NFAT pathway for LDL, respectively. In conclusion, the TC and LDL-C level is moderately heritable and the HDL-C level is lowly heritable in Chinese population. The genomic loci, functional genes and pathways are identified to account for the heritability of plasma cholesterol level. Our findings provide important insights into plasma cholesterol molecular physiology and expect future research to replicate and validate our results.