Genetic association of GWAS-supported MAD1L1 gene polymorphism rs12666575 with schizophrenia susceptibility in a Chinese population

Maternal genetic polymorphisms in the major mitotic checkpoint genes MAD1L1 and MAD2L1 associated with the risk of survival in abnormal chromosomal fetuses

Background: The genetic etiology of fetal chromosome abnormalities remains unknown, which brings about an enormous burden for patients, families, and society. The spindle assembly checkpoint (SAC) controls the normal procedure of chromosome disjunction and may take part in the process.

Objective: The aim of this study was to explore the association between polymorphisms of MAD1L1 rs1801368 and MAD2L1 rs1283639804, involved in SAC and fetal chromosome abnormalities.

Methods: The case–control study collected 563 cases and 813 health controls to test the genotypes of MAD1L1 rs1801368 and MAD2L1 rs1283639804 polymorphisms by polymerase chain reaction–restrictive fragment length polymorphism methods (PCR-RFLP).

Results:MAD1L1 rs1801368 polymorphism was associated with fetal chromosome abnormalities alone or combined to lower homocysteine (HCY) levels (alone: dominant: OR: 1.75, 95%CI: 1.19–2.57, and p = 0.005; CT vs. CC: OR = 0.73, 95%CI: 0.57–0.94, and p = 0.016; lower HCY: C vs. T: OR = 0.74, 95%CI: 0.57–0.95, and p = 0.02; dominant: OR = 1.75, 95%CI: 0.79–1.92, and p = 0.005). No significant differences were found in other genetic models or subgroups (p > 0.05, respectively). MAD2L1 rs1283639804 polymorphism revealed a sole genotype in the studied population. HCY is significantly associated with fetal chromosome abnormalities in younger groups (OR: 1.78, 95%CI: 1.28–2.47, and p = 0.001).

Conclusion: The results implied that the polymorphism of MAD1L1 rs1801368 may become the susceptibility factor to fetal chromosome abnormalities alone or combined to lower HCY levels but not to MAD2L1 rs1283639804 polymorphism. In addition, HCY significantly affects fetal chromosomal abnormalities in younger women.

Methylation in MAD1L1 is associated with the severity of suicide attempt and phenotypes of depression

Depression is a multifactorial disorder representing a significant public health burden. Previous studies have linked multiple single nucleotide polymorphisms with depressive phenotypes and suicidal behavior. MAD1L1 is a mitosis metaphase checkpoint protein that has been linked to depression in GWAS. Using a longitudinal EWAS approach in an adolescent cohort at two time points (n = 216 and n = 154), we identified differentially methylated sites that were associated with depression-related genetic variants in MAD1L1. Three methylation loci (cg02825527, cg18302629, and cg19624444) were consistently hypomethylated in the minor allele carriers, being cross-dependent on several SNPs. We further investigated whether DNA methylation at these CpGs is associated with depressive psychiatric phenotypes in independent cohorts. The first site (cg02825527) was hypomethylated in blood (exp(β) = 84.521, p value ~ 0.003) in participants with severe suicide attempts (n = 88). The same locus showed increased methylation in glial cells (exp(β) = 0.041, p value ~ 0.004) in the validation cohort, involving 29 depressed patients and 29 controls, and showed a trend for association with suicide (n = 40, p value ~ 0.089) and trend for association with depression treatment (n = 377, p value ~ 0.075). The second CpG (cg18302629) was significantly hypomethylated in depressed participants (exp(β) = 56.374, p value ~ 0.023) in glial cells, but did not show associations in the discovery cohorts. The last methylation site (cg19624444) was hypomethylated in the whole blood of severe suicide attempters; however, this association was at the borderline for statistical significance (p value ~ 0.061). This locus, however, showed a strong association with depression treatment in the validation cohort (exp(β) = 2.237, p value ~ 0.003) with 377 participants. The direction of associations between psychiatric phenotypes appeared to be different in the whole blood in comparison with brain samples for cg02825527 and cg19624444. The association analysis between methylation at cg18302629 and cg19624444 and MAD1L1 transcript levels in CD¹⁴⁺ cells shows a potential link between methylation at these CpGs and MAD1L1 expression. This study suggests evidence that methylation at MAD1L1 is important for psychiatric health as supported by several independent cohorts.

DNA methylation as a potential mediator of the association between prenatal tobacco and alcohol exposure and child neurodevelopment in a South African birth cohort

Prenatal tobacco exposure (PTE) and prenatal alcohol exposure (PAE) have been associated with an increased risk of delayed neurodevelopment in children as well as differential newborn DNA methylation (DNAm). However, the biological mechanisms connecting PTE and PAE, DNAm, and neurodevelopment are largely unknown. Here we aim to determine whether differential DNAm mediates the association between PTE and PAE and neurodevelopment at 6 (N = 112) and 24 months (N = 184) in children from the South African Drakenstein Child Health Study. PTE and PAE were assessed antenatally using urine cotinine measurements and the ASSIST questionnaire, respectively. Cord blood DNAm was measured using the EPIC and 450 K BeadChips. Neurodevelopment (cognitive, language, motor, adaptive behavior, socioemotional) was measured using the Bayley Scales of Infant and Toddler Development, Third Edition. We constructed methylation risk scores (MRS) for PTE and PAE and conducted causal mediation analysis (CMA) with these MRS as mediators. Next, we conducted a high-dimensional mediation analysis to identify individual CpG sites as potential mediators, followed by a CMA to estimate the average causal mediation effects (ACME) and total effect (TE). PTE and PAE were associated with neurodevelopment at 6 but not at 24 months. PTE MRS reached a prediction accuracy (R2) of 0.23 but did not significantly mediate the association between PTE and neurodevelopment. PAE MRS was not predictive of PAE (R2 = 0.006). For PTE, 31 CpG sites and eight CpG sites were identified as significant mediators (ACME and TE P < 0.05) for the cognitive and motor domains at 6 months, respectively. For PAE, 16 CpG sites and 1 CpG site were significant mediators for the motor and adaptive behavior domains at 6 months, respectively. Several of the associated genes, including MAD1L1, CAMTA1, and ALDH1A2 have been implicated in neurodevelopmental delay, suggesting that differential DNAm may partly explain the biological mechanisms underlying the relationship between PTE and PAE and child neurodevelopment.

Prenatal metal exposure, cord blood DNA methylation and persistence in childhood: an epigenome-wide association study of 12 metals

BACKGROUND

Prenatal exposure to essential and non-essential metals impacts birth and child health, including fetal growth and neurodevelopment. DNA methylation (DNAm) may be involved in pathways linking prenatal metal exposure and health. In the Project Viva cohort, we analyzed the extent to which metals (As, Ba, Cd, Cr, Cs, Cu, Hg, Mg, Mn, Pb, Se, and Zn) measured in maternal erythrocytes were associated with differentially methylated positions (DMPs) and regions (DMRs) in cord blood and tested if associations persisted in blood collected in mid-childhood. We measured metal concentrations in first-trimester maternal erythrocytes, and DNAm in cord blood (N = 361) and mid-childhood blood (N = 333, 6-10 years) with the Illumina HumanMethylation450 BeadChip. For each metal individually, we tested for DMPs using linear models (considered significant at FDR < 0.05), and for DMRs using comb-p (Sidak p < 0.05). Covariates included biologically relevant variables and estimated cell-type composition. We also performed sex-stratified analyses.

RESULTS

Pb was associated with decreased methylation of cg20608990 (CASP8) (FDR = 0.04), and Mn was associated with increased methylation of cg02042823 (A2BP1) in cord blood (FDR = 9.73 × 10-6). Both associations remained significant but attenuated in blood DNAm collected at mid-childhood (p < 0.01). Two and nine Mn-associated DMPs were identified in male and female infants, respectively (FDR < 0.05), with two and six persisting in mid-childhood (p < 0.05). All metals except Ba and Pb were associated with ≥ 1 DMR among all infants (Sidak p < 0.05). Overlapping DMRs annotated to genes in the human leukocyte antigen (HLA) region were identified for Cr, Cs, Cu, Hg, Mg, and Mn.

CONCLUSIONS

Prenatal metal exposure is associated with DNAm, including DMRs annotated to genes involved in neurodevelopment. Future research is needed to determine if DNAm partially explains the relationship between prenatal metal exposures and health outcomes.

Identification of 34 genes conferring genetic and pharmacological risk for the comorbidity of schizophrenia and smoking behaviors

The prevalence of smoking is significantly higher in persons with schizophrenia (SCZ) than in the general population. However, the biological mechanisms of the comorbidity of smoking and SCZ are largely unknown. This study aimed to reveal shared biological pathways for the two diseases by analyzing data from two genome-wide association studies with a total sample size of 153,898. With pathway-based analysis, we first discovered 18 significantly enriched pathways shared by SCZ and smoking, which were classified into five groups: postsynaptic density, cadherin binding, dendritic spine, long-term depression, and axon guidance. Then, by using an integrative analysis of genetic, epigenetic, and expression data, we found not only 34 critical genes (e.g., PRKCZ, ARHGEF3, and CDKN1A) but also various risk-associated SNPs in these genes, which convey susceptibility to the comorbidity of the two disorders. Finally, using both in vivo and in vitro data, we demonstrated that the expression profiles of the 34 genes were significantly altered by multiple psychotropic drugs. Together, this multi-omics study not only reveals target genes for new drugs to treat SCZ but also reveals new insights into the shared genetic vulnerabilities of SCZ and smoking behaviors.

Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder

BACKGROUND

Epigenetic mechanisms have been suggested to play a role in the development of post-traumatic stress disorder (PTSD). Here, blood-derived DNA methylation data (HumanMethylation450 BeadChip) collected prior to and following combat exposure in three cohorts of male military members were analyzed to assess whether DNA methylation profiles are associated with the development of PTSD. A total of 123 PTSD cases and 143 trauma-exposed controls were included in the analyses. The Psychiatric Genomics Consortium (PGC) PTSD EWAS QC pipeline was used on all cohorts, and results were combined using a sample size weighted meta-analysis in a two-stage design. In stage one, we jointly analyzed data of two new cohorts (N = 126 and 78) for gene discovery, and sought to replicate significant findings in a third, previously published cohort (N = 62) to assess the robustness of our results. In stage 2, we aimed at maximizing power for gene discovery by combining all three cohorts in a meta-analysis.

RESULTS

Stage 1 analyses identified four CpG sites in which, conditional on pre-deployment DNA methylation, post-deployment DNA methylation was significantly associated with PTSD status after epigenome-wide adjustment for multiple comparisons. The most significant (intergenic) CpG cg05656210 (p = 1.0 × 10-08) was located on 5q31 and significantly replicated in the third cohort. In addition, 19 differentially methylated regions (DMRs) were identified, but failed replication. Stage 2 analyses identified three epigenome-wide significant CpGs, the intergenic CpG cg05656210 and two additional CpGs located in MAD1L1 (cg12169700) and HEXDC (cg20756026). Interestingly, cg12169700 had an underlying single nucleotide polymorphism (SNP) which was located within the same LD block as a recently identified PTSD-associated SNP in MAD1L1. Stage 2 analyses further identified 12 significant differential methylated regions (DMRs), 1 of which was located in MAD1L1 and 4 were situated in the human leukocyte antigen (HLA) region.

CONCLUSIONS

This study suggests that the development of combat-related PTSD is associated with distinct methylation patterns in several genomic positions and regions. Our most prominent findings suggest the involvement of the immune system through the HLA region and HEXDC, and MAD1L1 which was previously associated with PTSD.

Association of MAD1L1 polymorphism (rs871925) with prenatal famine exposure and schizophrenia in a Chinese population: A case-control study

Schizophrenia (SCZ) is a complex, frequently disabling psychiatric disorder. Prenatal exposure to famine, an environmental factor, plays a significant role in the cause of SCZ. We used DNA methylation related sites to analyze their association with prenatal famine exposure and SCZ risk in a Northeast Han Chinese population. A total of 967 subjects (446 patients with SCZ/521 health controls) were recruited. Five single-nucleotide polymorphisms (rs2300149 in ITIH1, rs2675956 in NGEF, rs3758543 in NT5C2, rs7003288 in NA, and rs871925 in MAD1L1) were selected and genotyped. Genotype distribution and allele frequency analysis indicated that rs871925 was significantly associated with SCZ. We also found a significant association between prenatal exposure to famine and rs871925 in the recessive model in the health control group. The generalized multifactor dimensionality reduction analysis suggested a five-locus interaction model association with the risk of developing SCZ. Our data suggested that MAD1L1 rs871925 was associated with prenatal famine exposure and SCZ susceptibility in a Northeast Han Chinese population.

Replicated associations of FADS1, MAD1L1, and a rare variant at 10q26.13 with bipolar disorder in Chinese population

Genetic analyses of psychiatric illnesses, such as bipolar disorder (BPD), have revealed essential information regarding the underlying pathological mechanisms. While such studies in populations of European ancestry have achieved prominent success, understanding the genetic risk factors of these illnesses (especially BPD) in Chinese population remains an urgent task. Given the lack of genome-wide association study (GWAS) of BPD in Chinese population from Mainland China, replicating the previously reported GWAS hits in distinct populations will provide valuable information for future GWAS analysis in Han Chinese. In the present study, we have recruited 1146 BPD cases and 1956 controls from Mainland China for genetic analyses, as well as 65 Han Chinese brain amygdala tissues for mRNA expression analyses. Using this clinical sample, one of the largest Han Chinese BPD samples till now, we have conducted replication analyses of 21 single nucleotide polymorphisms (SNPs) extracted from previous GWAS of distinct populations. Among the 21 tested SNPs, 16 showed the same direction of allelic effects in our samples compared with previous studies; 6 SNPs achieved nominal significance (p < 0.05) at one-tailed test, and 2 additional SNPs showed marginal significance (p < 0.10). Aside from replicating previously reported BPD risk SNPs, we herein also report several intriguing findings: (1) the SNP rs174576 was associated with BPD in our Chinese sample and in the overall global meta-analysis, and was significantly correlated with FADS1 mRNA in diverse public RNA-seq datasets as well as our in house collected Chinese amygdala samples; (2) two (partially) independent SNPs in MAD1L1 were both significantly associated with BPD in our Chinese sample, which was also supported by haplotype analysis; (3) a rare SNP rs78089757 in 10q26.13 region was a genome-wide significant variant for BPD in East Asians, and this SNP was near monomorphic in Europeans. In sum, these results confirmed several significant BPD risk genes. We hope this Chinese BPD case-control sample and the current brain amygdala tissues (with continuous increasing sample size in the near future) will provide helpful resources in elucidating the genetic and molecular basis of BPD in this major world population.

Regional enrichment analyses on genetic profiles for schizophrenia and bipolar disorder

Both schizophrenia (SZ) and bipolar disorder (BD) are highly heritable psychiatric disorders. The significant genomic risk loci are of great importance but with no guarantee of known functional impact and they cannot totally explain the genetic inheritance. In this study we present regional enrichment analyses across the genome, aiming to strike a balance between individual risk loci and integrated regional effects. Chromosomes were partitioned into 2 million base-pair regions (indicated by an underscore sign in the cytogenetic bands) on which enrichment tests are performed. In the discovery phase, we leverage the Psychiatric Genomics Consortium SZ and BD initial association test results for European Ancestry (EA) population and dbGAP SNP data for African Ancestry (AA) population. 78 and 48 regions show significantly enriched associations with SZ and BD respectively in the EA population, and nine are in common including MHC, 3p21.1, 7p22.3_2, 2q32.3_2, 8q24.3_4, and 19q13.33_1. The most unique SZ associated region is 1p21.3_3, while the most unique BD associated region is 6q25.2_1. For the AA population fewer regions are discovered with only 10% overlapping with that of EA population. A replication test using Wellcome Trust Case Control Consortium data for EA population verified 9% of the SZ enriched regions and 40% of the BD enriched regions. In summary, we showed that regional enrichment analyses produce reliable genetic association profiles using about one tenth of samples compared to single base-pair genome wide association approach. The identified association regions will be useful for further genetic functional studies.

Progress in genome-wide association studies of schizophrenia in Han Chinese populations

Since 2006, genome-wide association studies of schizophrenia have led to the identification of numerous novel risk loci for this disease. However, there remains a geographical imbalance in genome-wide association studies, which to date have primarily focused on Western populations. During the last 6 years, genome-wide association studies in Han Chinese populations have identified both the sharing of susceptible loci across ethnicities and genes unique to Han Chinese populations. Here, we review recent progress in genome-wide association studies of schizophrenia in Han Chinese populations. Researchers have identified and replicated the sharing of susceptible genes, such as within the major histocompatibility complex, microRNA 137 (MIR137), zinc finger protein 804A (ZNF804A), vaccinia related kinase 2 (VRK2), and arsenite methyltransferase (AS3MT), across both European and East Asian populations. Several copy number variations identified in European populations have also been validated in the Han Chinese, including duplications at 16p11.2, 15q11.2-13.1, 7q11.23, and VIPR2 and deletions at 22q11.2, 1q21.1-q21.2, and NRXN1. However, these studies have identified some potential confounding factors, such as genetic heterogeneity and the effects of natural selection on tetraspanin 18 (TSPAN18) or zinc finger protein 323 (ZNF323), which may explain the population differences in genome-wide association studies. In the future, genome-wide association studies in Han Chinese populations should include meta-analyzes or mega-analyses with enlarged sample sizes across populations, deep sequencing, precision medicine treatment, and functional exploration of the risk genes for schizophrenia.