Resequencing Three Candidate Genes for Major Depressive Disorder in a Dutch Cohort

Major depressive disorder (MDD) is a psychiatric disorder, characterized by periods of low mood of more than two weeks, loss of interest in normally enjoyable activities and behavioral changes. MDD is a complex disorder and does not have a single genetic cause. In 2009 a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. Many of the top signals of this GWAS mapped to a region spanning the gene PCLO, and the non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became genome wide significant after post-hoc analysis. We performed resequencing of PCLO, GRM7, and SLC6A4 in 50 control samples from the GAIN-MDD cohort, to detect new genomic variants. Subsequently, we genotyped these variants in the entire GAIN-MDD cohort and performed association analysis to investigate if rs2522833 is the causal variant or simply in linkage disequilibrium with a more associated variant. GRM7 and SLC6A4 are both candidate genes for MDD from literature. We aimed to gather more evidence that rs2522833 is indeed the causal variant in the GAIN-MDD cohort or to find a previously undetected common variant in either PCLO, GRM7, or SLC6A4 with a higher association in this cohort. After next generation sequencing and association analysis we excluded the possibility of an undetected common variant to be more associated. For neither PCLO nor GRM7 we found a more associated variant. For SLC6A4, we found a new SNP that showed a lower P-value (P = 0.07) than in the GAIN-MDD GWAS (P = 0.09). However, no evidence for genome-wide significance was found. Although we did not take into account rare variants, we conclude that our results provide further support for the hypothesis that the non-synonymous coding SNP rs2522833 in the PCLO gene is indeed likely to be the causal variant in the GAIN-MDD cohort.     

Identification of a novel prostate cancer susceptibility variant in the KLK3 gene transcript

Genome-wide association studies (GWAS) have identified more than 30 prostate cancer (PrCa) susceptibility loci. One of these (rs2735839) is located close to a plausible candidate susceptibility gene, KLK3, which encodes prostate-specific antigen (PSA). PSA is widely used as a biomarker for PrCa detection and disease monitoring. To refine the association between PrCa and variants in this region, we used genotyping data from a two-stage GWAS using samples from the UK and Australia, and the Cancer Genetic Markers of Susceptibility (CGEMS) study. Genotypes were imputed for 197 and 312 single nucleotide polymorphisms (SNPs) from HapMap2 and the 1000 Genome Project, respectively. The most significant association with PrCa was with a previously unidentified SNP, rs17632542 (combined P?=?3.9?×?10(-22)). This association was confirmed by direct genotyping in three stages of the UK/Australian GWAS, involving 10,405 cases and 10,681 controls (combined P?=?1.9?×?10(-34)). rs17632542 is also shown to be associated with PSA levels and it is a non-synonymous coding SNP (Ile179Thr) in KLK3. Using molecular dynamic simulation, we showed evidence that this variant has the potential to introduce alterations in the protein or affect RNA splicing. We propose that rs17632542 may directly influence PrCa risk.     

Comprehensive SNP Scan of DNA Repair and DNA Damage Response Genes Reveal Multiple Susceptibility Loci Conferring Risk to Tobacco Associated Leukoplakia and Oral Cancer

Polymorphic variants of DNA repair and damage response genes play major role in carcinogenesis. These variants are suspected as predisposition factors to Oral Squamous Cell Carcinoma (OSCC). For identification of susceptible variants affecting OSCC development in Indian population, the “maximally informative” method of SNP selection from HapMap data to non-HapMap populations was applied. Three hundred twenty-five SNPs from 11 key genes involved in double strand break repair, mismatch repair and DNA damage response pathways were genotyped on a total of 373 OSCC, 253 leukoplakia and 535 unrelated control individuals. The significantly associated SNPs were validated in an additional cohort of 144 OSCC patients and 160 controls. The rs12515548 of MSH3 showed significant association with OSCC both in the discovery and validation phases (discovery P-value: 1.43E-05, replication P-value: 4.84E-03). Two SNPs (rs12360870 of MRE11A, P-value: 2.37E-07 and rs7003908 of PRKDC, P-value: 7.99E-05) were found to be significantly associated only with leukoplakia. Stratification of subjects based on amount of tobacco consumption identified SNPs that were associated with either high or low tobacco exposed group. The study reveals a synergism between associated SNPs and lifestyle factors in predisposition to OSCC and leukoplakia.     

Genome-wide association analysis of circulating vitamin D levels in children with asthma

Vitamin D deficiency is becoming more apparent in many populations. Genetic factors may play a role in the maintenance of vitamin D levels. The objective of this study was to perform a genome-wide analysis (GWAS) of vitamin D levels, including replication of prior GWAS results. We measured 25-hydroxyvitamin D (25(OH)D) levels in serum collected at the time of enrollment and at year 4 in 572 Caucasian children with asthma, who were part of a multi-center clinical trial, the Childhood Asthma Management Program. Replication was performed in a second cohort of 592 asthmatics from Costa Rica and a third cohort of 516 Puerto Rican asthmatics. In addition, we attempted replication of three SNPs that were previously identified in a large GWAS of Caucasian individuals. The setting included data from a clinical trial of childhood asthmatics and two cohorts of asthmatics recruited for genetic studies of asthma. The main outcome measure was circulating 25(OH)D levels. The 25(OH)D levels at the two time-points were only modestly correlated with each other (intraclass correlation coefficient?=?0.33) in the CAMP population. We identified SNPs that were nominally associated with 25(OH)D levels at two time-points in CAMP, and replicated four SNPs in the Costa Rican cohort: rs11002969, rs163221, rs1678849, and rs4864976. However, these SNPs were not significantly associated with 25(OH)D levels in a third population of Puerto Rican asthmatics. We were able to replicate the SNP with the strongest effect, previously reported in a large GWAS: rs2282679 (GC), and we were able to replicate another SNP, rs10741657 (CYP2R1), to a lesser degree. We were able to replicate two of three prior significant findings in a GWAS of 25(OH)D levels. Other SNPs may be additionally associated with 25(OH)D levels in certain populations.     

Identification of sialyltransferase 8B as a generalized susceptibility gene for psychotic and mood disorders on chromosome 15q25-26

We previously identified a significant bipolar spectrum disorder linkage peak on 15q25-26 using 35 extended families with a broad clinical phenotype, including bipolar disorder (types I and II), recurrent unipolar depression and schizoaffective disorder. However, the specific gene(s) contributing to this signal had not been identified. By a fine mapping association study in an Australian case-control cohort (n?=?385), we find that the sialyltransferase 8B (ST8SIA2) gene, coding for an enzyme that glycosylates proteins involved in neuronal plasticity which has previously shown association to both schizophrenia and autism, is associated with increased risk to bipolar spectrum disorder. Nominal single point association was observed with SNPs in ST8SIA2 (rs4586379, P?=?0.0043; rs2168351, P?=?0.0045), and a specific risk haplotype was identified (frequency: bipolar vs controls?=?0.41 vs 0.31; χ(2)?=?6.46, P?=?0.011, OR?=?1.47). Over-representation of the specific risk haplotype was also observed in an Australian schizophrenia case-control cohort (n?=?256) (χ(2)?=?8.41, P?=?0.004, OR?=?1.82). Using GWAS data from the NIMH bipolar disorder (n?=?2055) and NIMH schizophrenia (n?=?2550) cohorts, the equivalent haplotype was significantly over-represented in bipolar disorder (χ(2)?=?5.91, P?=?0.015, OR?=?1.29), with the same direction of effect in schizophrenia, albeit non-significant (χ(2)?=?2.3, P?=?0.129, OR?=?1.09). We demonstrate marked down-regulation of ST8SIA2 gene expression across human brain development and show a significant haplotype×diagnosis effect on ST8SIA2 mRNA levels in adult cortex (ANOVA: F(1,87)?=?6.031, P?=?0.016). These findings suggest that variation the ST8SIA2 gene is associated with increased risk to mental illness, acting to restrict neuronal plasticity and disrupt early neuronal network formation, rendering the developing and adult brain more vulnerable to secondary genetic or environmental insults.     

Functional analysis of HapMap SNPs

Genome-wide association studies (GWAS) have successfully identified many genetic variants associated with complex diseases and traits. However, functional consequence of genetic variants studied in GWAS is not yet fully investigated, which would hinder the application of GWAS. We therefore performed a systematic functional analysis of HapMap SNPs, which have been most commonly used as the reference panel for GWAS. Our study highlights several characteristics of HapMap SNPs and identifies subsets of genetic variants with interesting functional implication. The results show that HapMap SNPs have good coverage within RefSeq genes, especially within known disease-related genes. On the other hand, only a small percentage of SNPs are non-synonymous SNPs while many SNPs are actually located at gene deserts. Moreover, many functionally important variants are not yet still interrogated. A redesigned SNP reference panel with additional functionally important variants would be useful to identify disease-causal variants in the future genome-wide studies.     

Exonic Re-Sequencing of the Chromosome 2q24.3 Parkinson’s Disease Locus

Genome-wide association studies (GWAS) in Parkinson’s disease (PD) have identified over 20 genomic regions associated with disease risk. Many of these loci include several candidate genes making it difficult to pinpoint the causal gene. The locus on chromosome 2q24.3 encompasses three genes: B3GALT1, STK39, and CERS6. In order to identify if the causal variants are simple missense changes, we sequenced all 31 exons of these three genes in 187 patients with PD. We identified 13 exonic variants including four non-synonymous and three insertion/deletion variants (indels). These non-synonymous variants and rs2102808, the GWAS tag SNP, were genotyped in three independent series consisting of a total of 1976 patients and 1596 controls. Our results show that the seven identified 2q24.3 coding variants are not independently responsible for the GWAS association signal at the locus; however, there is a haplotype, which contains both rs2102808 and a STK39 exon 1 6bp indel variant, that is significantly associated with PD risk (Odds Ratio [OR] = 1.35, 95% CI: 1.11–1.64, P = 0.003). This haplotype is more associated than each of the two variants independently (OR = 1.23, P = 0.005 and 1.10, P = 0.10, respectively). Our findings suggest that the risk variant is likely located in a non-coding region. Additional sequencing of the locus including promoter and regulatory regions will be needed to pinpoint the association at this locus that leads to an increased risk to PD.     

An analysis of single nucleotide polymorphisms of 125 DNA repair genes in the Texas genome-wide association study of lung cancer with a replication for the XRCC4 SNPs

DNA repair genes are important for maintaining genomic stability and limiting carcinogenesis. We analyzed all single nucleotide polymorphisms (SNPs) of 125 DNA repair genes covered by the Illumina HumanHap300 (v1.1) BeadChips in a previously conducted genome-wide association study (GWAS) of 1154 lung cancer cases and 1137 controls and replicated the top-hits of XRCC4 SNPs in an independent set of 597 cases and 611 controls in Texas populations. We found that six of 20 XRCC4 SNPs were associated with a decreased risk of lung cancer with a P-value of 0.01 or lower in the discovery dataset, of which the most significant SNP was rs10040363 (P for allelic test=4.89 x 10??). Moreover, the data in this region allowed us to impute a potentially functional SNP rs2075685 (imputed P for allelic test=1.3 x 10?3). A luciferase reporter assay demonstrated that the rs2075685G>T change in the XRCC4 promoter increased expression of the gene. In the replication study of rs10040363, rs1478486, rs9293329, and rs2075685, however, only rs10040363 achieved a borderline association with a decreased risk of lung cancer in a dominant model (adjusted OR=0.80, 95% CI=0.62-1.03 and P=0.079). In the final combined analysis of both the Texas GWAS discovery and replication datasets, the strength of the association was increased for rs10040363 (adjusted OR=0.77, 95% CI=0.66-0.89, P(dominant)=5 x 10?? and P for trend=5 x 10??) and rs1478486 (adjusted OR=0.82, 95% CI=0.71-0.94, P(dominant)=6 x 10?3 and P for trend=3.5 x 10?3). Finally, we conducted a meta-analysis of these XRCC4 SNPs with available data from published GWA studies of lung cancer with a total of 12,312 cases and 47,921 controls, in which none of these XRCC4 SNPs was associated with lung cancer risk. It appeared that rs2075685, although associated with increased expression of a reporter gene and lung cancer risk in the Texas populations, did not have an effect on lung cancer risk in other populations. This study underscores the importance of replication using published data in larger populations.     

IL12B and IL23R gene SNPs in Japanese psoriasis

Psoriasis is a common human skin disease whereby abnormal production of inflammatory mediators is believed to play an important role in its pathogenesis. The IL12B gene, which encodes the shared IL-12p40 subunit in two cytokines, IL-12 and IL-23, and the IL23R gene, which encodes a subunit of the receptor for IL-23, were identified as psoriasis-susceptibility genetic factors in recent candidate gene and genome-wide association studies of Chinese and Europeans. Since there are significant differences in the incidence of psoriasis between Europeans and Japanese suggesting a genetic ethnic effect, we examined the association of IL12B and IL23R gene polymorphisms with psoriasis in a cohort of Japanese. In this study, we genotyped two SNPs (rs3212227 and rs6887695) in the IL12B gene and one SNP (rs11209026) in the IL23R gene using 560 Japanese psoriasis cases and 560 controls and compared our results with those previously published for Europeans and Asians. Our study showed significant associations between psoriasis and both IL12B gene SNPs, rs3212227 (odds ratio (OR)?=?1.35, P?=?4.94E-04) and rs6887695 (OR?=?1.32, P?=?2.00E-03), but no significant association between psoriasis and the IL23R SNP, rs11209026. Furthermore, a significant haplotype association was found for the IL12B gene protective haplotype C-C (OR?=?0.71, P?=?1.84E-04) in Japanese, as previously elucidated in the studies of European ancestry.     

EphB2 SNPs and sporadic prostate cancer risk in African American men

The EphB2 gene has been implicated as a tumor suppressor gene somatically altered in both prostate cancer (PC) and colorectal cancer. We have previously shown an association between an EphB2 germline nonsense variant and risk of familial prostate cancer among African American Men (AAM). Here we set out to test the hypothesis that common variation within the EphB2 locus is associated with increased risk of sporadic PC in AAM. We genotyped a set of 341 single nucleotide polymorphisms (SNPs) encompassing the EphB2 locus, including known and novel coding and noncoding variants, in 490 AA sporadic PC cases and 567 matched controls. Single marker-based logistical regression analyses revealed seven EphB2 SNPs showing statistically significant association with prostate cancer risk in our population. The most significant association was achieved for a novel synonymous coding SNP, TGen-624, (Odds Ratio (OR) =?0.22; 95% Confidence Interval (CI) 0.08-0.66, p?=?1×10(-5)). Two other SNPs also show significant associations toward a protective effect rs10465543 and rs12090415 (p?=?1×10(-4)), OR?=?0.49 and 0.7, respectively. Two additional SNPs revealed trends towards an increase in risk of prostate cancer, rs4612601 and rs4263970 (p?=?0.001), OR?=?1.35 and 1.31, respectively. Furthermore, haplotype analysis revealed low levels of linkage disequilibrium within the region, with two blocks being associated with prostate cancer risk among our population. These data suggest that genetic variation at the EphB2 locus may increase risk of sporadic PC among AAM.     

A Germline Polymorphism of DNA Polymerase Beta Induces Genomic Instability and Cellular Transformation

Several germline single nucleotide polymorphisms (SNPs) have been identified in the POLB gene, but little is known about their cellular and biochemical impact. DNA Polymerase β (Pol β), encoded by the POLB gene, is the main gap-filling polymerase involved in base excision repair (BER), a pathway that protects the genome from the consequences of oxidative DNA damage. In this study we tested the hypothesis that expression of the POLB germline coding SNP (rs3136797) in mammalian cells could induce a cancerous phenotype. Expression of this SNP in both human and mouse cells induced double-strand breaks, chromosomal aberrations, and cellular transformation. Following treatment with an alkylating agent, cells expressing this coding SNP accumulated BER intermediate substrates, including single-strand and double-strand breaks. The rs3136797 SNP encodes the P242R variant Pol β protein and biochemical analysis showed that P242R protein had a slower catalytic rate than WT, although P242R binds DNA similarly to WT. Our results suggest that people who carry the rs3136797 germline SNP may be at an increased risk for cancer susceptibility.     

Genome-wide association analysis in asthma subjects identifies SPATS2L as a novel bronchodilator response gene

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV(1)) before and after the administration of a short-acting β(2)-agonist, the most common rescue medications used for the treatment of asthma. BDR also serves as a test of β(2)-agonist efficacy. BDR is a complex trait that is partly under genetic control. A genome-wide association study (GWAS) of BDR, quantified as percent change in baseline FEV(1) after administration of a β(2)-agonist, was performed with 1,644 non-Hispanic white asthmatic subjects from six drug clinical trials: CAMP, LOCCS, LODO, a medication trial conducted by Sepracor, CARE, and ACRN. Data for 469,884 single-nucleotide polymorphisms (SNPs) were used to measure the association of SNPs with BDR using a linear regression model, while adjusting for age, sex, and height. Replication of primary P-values was attempted in 501 white subjects from SARP and 550 white subjects from DAG. Experimental evidence supporting the top gene was obtained via siRNA knockdown and Western blotting analyses. The lowest overall combined P-value was 9.7E-07 for SNP rs295137, near the SPATS2L gene. Among subjects in the primary analysis, those with rs295137 TT genotype had a median BDR of 16.0 (IQR = [6.2, 32.4]), while those with CC or TC genotypes had a median BDR of 10.9 (IQR = [5.0, 22.2]). SPATS2L mRNA knockdown resulted in increased β(2)-adrenergic receptor levels. Our results suggest that SPATS2L may be an important regulator of β(2)-adrenergic receptor down-regulation and that there is promise in gaining a better understanding of the biological mechanisms of differential response to β(2)-agonists through GWAS.     

Combinations of Susceptibility Genes Are Associated with Higher Risk for Multiple Sclerosis and Imply Disease Course Specificity

Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system that predominantly affects young adults. The genetic contributions to this multifactorial disease were underscored by a genome wide association study (GWAS) conducted by the International Multiple Sclerosis Genetic Consortium in a multinational cohort prompting the discovery of 57 non-MHC MS-associated common genetic variants. Hitherto, few of these newly reported variants have been replicated in larger independent patient cohorts. We genotyped a cohort of 1033 MS patients and 644 healthy controls with a consistent genetic background for the 57 non-MHC variants reported to be associated with MS by the first large GWAS as well as the HLA DRB1*1501 tagging SNP rs3135388. We robustly replicated three of the 57 non-MHC reported MS-associated single nucleotide polymorphisms (SNPs). In addition, our study revealed several genotype-genotype combinations with an evidently higher degree of disease association than the genotypes of the single SNPs. We further correlated well-defined clinical phenotypes, i.e. ataxia, visual impairment due to optic neuritis and paresis with single SNPs and genotype combinations, and identified several associations. The results may open new avenues for clinical implications of the MS associated genetic variants reported from large GWAS.     

Analysis of single nucleotide polymorphism rs415430 in the WNT3 gene in the Russian population with the Parkinson disease

The Parkinson disease (PD) is the second most common progressive neurodegenerative disorder that arises due to degeneration of dopaminergic neurons. The causes of this disease are still unknown, but a number of genes involved in pathogenesis of familial and sporadic forms of PD has been identified. According to recent data of genome wide association studies (GWAS), single nucleotide polymorphisms (SNPs) in these genes (including MAPT locus) may play an important role in the development of PD. Therefore, we analyzed distribution of genotype frequencies of SNP rs415430 in the WNT3 gene in the Russian patients with sporadic PD and in the Russian population controls (OR = 0.84, Confidence Interval (95% CI) 0.58-1.23, p = 0.39). It was concluded that SNP rs415430 in the WNT3 gene was not associated with the risk of development of PD.     

A replication study of GWAS-derived lipid genes in Asian Indians: the chromosomal region 11q23.3 harbors loci contributing to triglycerides

Recent genome-wide association scans (GWAS) and meta-analysis studies on European populations have identified many genes previously implicated in lipid regulation. Validation of these loci on different global populations is important in determining their clinical relevance, particularly for development of novel drug targets for treating and preventing diabetic dyslipidemia and coronary artery disease (CAD). In an attempt to replicate GWAS findings on a non-European sample, we examined the role of six of these loci (CELSR2-PSRC1-SORT1 rs599839; CDKN2A-2B rs1333049; BUD13-ZNF259 rs964184; ZNF259 rs12286037; CETP rs3764261; APOE-C1-C4-C2 rs4420638) in our Asian Indian cohort from the Sikh Diabetes Study (SDS) comprising 3,781 individuals (2,902 from Punjab and 879 from the US). Two of the six SNPs examined showed convincing replication in these populations of Asian Indian origin. Our study confirmed a strong association of CETP rs3764261 with high-density lipoprotein cholesterol (HDL-C) (p?=?2.03×10(-26)). Our results also showed significant associations of two GWAS SNPs (rs964184 and rs12286037) from BUD13-ZNF259 near the APOA5-A4-C3-A1 genes with triglyceride (TG) levels in this Asian Indian cohort (rs964184: p?=?1.74×10(-17); rs12286037: p?=?1.58×10(-2)). We further explored 45 SNPs in a ～195 kb region within the chromosomal region 11q23.3 (encompassing the BUD13-ZNF259, APOA5-A4-C3-A1, and SIK3 genes) in 8,530 Asian Indians from the London Life Sciences Population (LOLIPOP) (UK) and SDS cohorts. Five more SNPs revealed significant associations with TG in both cohorts individually as well as in a joint meta-analysis. However, the strongest signal for TG remained with BUD13-ZNF259 (rs964184: p?=?1.06×10(-39)). Future targeted deep sequencing and functional studies should enhance our understanding of the clinical relevance of these genes in dyslipidemia and hypertriglyceridemia (HTG) and, consequently, diabetes and CAD.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

A haplotype of glycogen synthase kinase 3β is associated with early onset of unipolar major depression

Recent findings suggest that glycogen synthase kinase 3β (GSK3β) may play a role in the pathophysiology and treatment of mood disorders. Various genetic studies have shown the association of GSK3β polymorphisms with different mood disorder phenotypes. We hypothesized that genetic variants in the GSK3β gene could partially underlie the susceptibility to mood disorders. We performed a genetic case–control study of 440 psychiatrically screened control subjects and 445 mood disorder patients [256 unipolar major depressive disorder (MDD) and 189 bipolar disorder (BD)]. We genotyped a set of 11 single nucleotide polymorphisms (SNPs) and determined the relative frequency of a known copy number variant (CNV) overlapping the GSK3β by quantitative real‐time polymerase chain reaction (PCR). We found no evidence of association with MDD or BD diagnosis, and we further investigated the age at onset (AAO) of the disorder and severity of depressive index episode. We found that rs334555, located in intron 1 of GSK3β, was nominally associated with an earlier AAO of the disease in MDD (P = 0.001). We also identified a haplotype containing three SNPs (rs334555, rs119258668 and rs11927974) associated with AAO of the disorder (permutated P = 0.0025). We detected variability for the CNV, but we could not detect differences between patients and controls for any of the explored phenotypes. This study presents further evidence of the contribution of GSK3β to mood disorders, implicating a specific SNP and a haplotype with an earlier onset of the disorder in a group of well‐characterized patients with unipolar MDD. Further replication studies in patients with the same phenotypic characteristics should confirm the results reported here.     

Association of polymorphisms in the Chr18q11.2 locus with tuberculosis in Chinese population

A GWAS study has reported that two single nucleotide polymorphisms (SNPs) were associated with predisposition to tuberculosis (TB) in African populations. These two loci represented the long-waited GWAS hits for TB susceptibility. To determine whether these two SNPs are associated with TB in Chinese population, we attempted an replication in a cohort of over one thousand Chinese TB patients and 1,280 healthy controls using melting temperature shift allele-specific genotyping analysis. We found that only SNP rs4331426 was significantly associated with TB in Chinese population (p = 0.011). However, the effect was opposite. The G allele of the SNP in Chinese population is a protective allele (OR = 0.62, 95 % CI 0.44–0.87), while it was the risk allele for African population (OR = 1.19, 95 % CI 1.12–1.26). No significance was found for SNP rs2335704. The results provided an independent support for a role in susceptibility to TB for SNP rs4331426. However, it also indicated that direct predisposition element to TB and the association effects may vary across ethnic groups.     

A Two-Stage Association Study Suggests BRAP as a Susceptibility Gene for Schizophrenia

Schizophrenia (SZ) is a neurodevelopmental disorder in which altered immune function typically plays an important role in mediating the effect of environmental insults and regulation of inflammation. The breast cancer suppressor protein associated protein (BRAP) is suggested to exert vital effects in neurodevelopment by modulating the mitogen-activated protein kinase cascade and inflammation signaling. To explore the possible role of BRAP in SZ, we conducted a two-stage study to examine the association of BRAP polymorphisms with SZ in the Han Chinese population. In stage one, we screened SNPs in BRAP from our GWAS data, which detected three associated SNPs, with rs3782886 being the most significant one (P  =  2.31E-6, OR  =  0.67). In stage two, we validated these three SNPs in an independently collected population including 1957 patients and 1509 controls, supporting the association of rs3782886 with SZ (P  =  1.43E-6, OR  =  0.73). Furthermore, cis-eQTL analysis indicates that rs3782886 genotypes are associated with mRNA levels of aldehyde dehydrogenase 2 family (ALDH2) (P  =  0.0039) and myosin regulatory light chain 2 (MYL2) (P < 1.0E-4). Our data suggest that the BRAP gene may confer vulnerability for SZ in Han Chinese population, adding further evidence for the involvement of developmental and/or neuroinflammatory cascades in the illness.     

Placental Genome and Maternal-Placental Genetic Interactions: A Genome-Wide and Candidate Gene Association Study of Placental Abruption

While available evidence supports the role of genetics in the pathogenesis of placental abruption (PA), PA-related placental genome variations and maternal-placental genetic interactions have not been investigated. Maternal blood and placental samples collected from participants in the Peruvian Abruptio Placentae Epidemiology study were genotyped using Illumina’s Cardio-Metabochip platform. We examined 118,782 genome-wide SNPs and 333 SNPs in 32 candidate genes from mitochondrial biogenesis and oxidative phosphorylation pathways in placental DNA from 280 PA cases and 244 controls. We assessed maternal-placental interactions in the candidate gene SNPS and two imprinted regions (IGF2/H19 and C19MC). Univariate and penalized logistic regression models were fit to estimate odds ratios. We examined the combined effect of multiple SNPs on PA risk using weighted genetic risk scores (WGRS) with repeated ten-fold cross-validations. A multinomial model was used to investigate maternal-placental genetic interactions. In placental genome-wide and candidate gene analyses, no SNP was significant after false discovery rate correction. The top genome-wide association study (GWAS) hits were rs544201, rs1484464 (CTNNA2), rs4149570 (TNFRSF1A) and rs13055470 (ZNRF3) (p-values: 1.11e-05 to 3.54e-05). The top 200 SNPs of the GWAS overrepresented genes involved in cell cycle, growth and proliferation. The top candidate gene hits were rs16949118 (COX10) and rs7609948 (THRB) (p-values: 6.00e-03 and 8.19e-03). Participants in the highest quartile of WGRS based on cross-validations using SNPs selected from the GWAS and candidate gene analyses had a 8.40-fold (95% CI: 5.8–12.56) and a 4.46-fold (95% CI: 2.94–6.72) higher odds of PA compared to participants in the lowest quartile. We found maternal-placental genetic interactions on PA risk for two SNPs in PPARG (chr3∶12313450 and chr3∶12412978) and maternal imprinting effects for multiple SNPs in the C19MC and IGF2/H19 regions. Variations in the placental genome and interactions between maternal-placental genetic variations may contribute to PA risk. Larger studies may help advance our understanding of PA pathogenesis.