Association Study of 167 Candidate Genes for Schizophrenia Selected by a Multi-Domain Evidence-Based Prioritization Algorithm and Neurodevelopmental Hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this
ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

Genetic variation in the 6p22.3 gene DTNBP1, the human ortholog of the mouse dysbindin gene, is associated with schizophrenia

Prior evidence has supported the existence of multiple susceptibility genes for schizophrenia. Multipoint linkage analysis of the 270 Irish high-density pedigrees that we have studied, as well as results from several other samples, suggest that at least one such gene is located in region 6p24-21. In the present study, family-based association analysis of 36 simple sequence-length-polymorphism markers and of 17 SNP markers implicated two regions, separated by approximately 7 Mb. The first region, and the focus of this report, is 6p22.3. In this region, single-nucleotide polymorphisms within the 140-kb gene DTNBP1 (dystrobrevin-binding protein 1, or dysbindin) are strongly associated with schizophrenia. Uncorrected, empirical P values produced by the program TRANSMIT were significant (P

Identification of a high-risk haplotype for the dystrobrevin binding protein 1 (DTNBP1) gene in the Irish study of high-density schizophrenia families

A recent report showed significant associations between several SNPs in a previously unknown EST cluster with schizophrenia. (1). The cluster was identified as the human dystrobrevin binding protein 1 gene (DTNBP1) by sequence database comparisons and homology with mouse DTNBP1. (2). However, the linkage disequilibrium (LD) among the SNPs in DTNBP1 as well as the pattern of significant SNP-schizophrenia association was complex. This raised several questions such as the number of susceptibility alleles that may be involved and the size of the region where the actual disease mutation(s) could be located. To address these questions, we performed different single-marker tests on the 12 previously studied and 2 new SNPs in DTNBP1 that were re-scored using an improved procedure, and performed a variety of haplotype analyses. The sample consisted of 268 Irish multiplex families selected for high density of schizophrenia. Results suggested a simple structure where the LD in the target region could be explained by 6 haplotypes that together accounted for 96% of haplotype diversity in the whole sample. From these six, a single high-risk haplotype was identified that showed a significant association with schizophrenia and explained the pattern of significant findings in the analyses with individual markers. This haplotype was 30 kb long, had a large effect, could be measured with two tag SNPs only, had a frequency of 6% in our sample, seemed to be of relatively recent origin in evolutionary terms, and was equally distributed over Ireland. Implications of these findings for follow-up and replication studies are discussed.

Variants in the catechol-o-methyltransferase (COMT) gene are associated with schizophrenia in Irish high-density families

The enzyme catechol-o-methyltransferase (COMT) transfers a methyl group from adenosylmethionine to catecholamines including the neurotransmitters dopamine, epinephrine and norepinephrine. This methylation results in the degradation of catecholamines. The involvement of the COMT gene in the metabolic pathway of these neurotransmitters has made it an attractive candidate gene for many psychiatric disorders. In this article, we reported our study of association of COMT with schizophrenia in Irish families with a high density of schizophrenia. Three single nucleotide polymorphisms (SNPs) were genotyped for the 274 such families and within-family transmission disequilibrium tests were performed. SNP rs4680, which is the functional Val/Met polymorphism, showed modest association with the disease by the TRANSMIT, FBAT and PDT programs, while the other two SNPs were negative. These SNPs showed lower level of LDs with each other in the Irish subjects than in Ashkenazi Jews. Haplotype analysis indicated that a haplotype, haplotype A-G-A for SNPs rs737865-rs4680-rs165599, was preferentially transmitted to the affected subjects. This was different from the reported G-G-G haplotype found in Ashkenazi Jews, but both haplotypes shared the Val allele. We concluded that COMT gene is associated with schizophrenia and carries a small but significant risk to the susceptibility in the Irish subjects.

Regulator of G-protein signaling 4 (RGS4) gene is associated with schizophrenia in Irish high density families

The regulator of the G-protein signaling 4 (RGS4) gene was shown to have a different expression pattern in schizophrenia patients in a microarray study. A family-based study subsequently implicated the association of this gene with schizophrenia. We replicated the study with our sample from the Irish Study of High Density Schizophrenia Families (ISHDSF). Single marker transmission disequilibrium tests (TDT) for the four core SNPs showed modest association for SNP 18 (using a narrow diagnostic approach with FBAT P = 0.044; with PDT P = 0.0073) and a trend for SNP 4 (with FBAT P = 0.1098; with PDT P = 0.0249). For SNP 1 and 7, alleles overtransmitted to affected subjects were the same as previously reported. Haplotype analyses suggested that haplotype G-G-G for SNP1-4-18, which is the most abundant haplotype (42.3%) in the Irish families, was associated with the disease (narrow diagnosis, FBAT P = 0.0061, PDT P = 0.0498). This was the same haplotype implicated in the original study. While P values were not corrected for multiple testing because of the clear prior hypothesis, these results could be interpreted as supporting evidence for the association between RGS4 and schizophrenia.

Relationship between a high-risk haplotype in the DTNBP1 (dysbindin) gene and clinical features of schizophrenia

The purpose of this study was to determine whether a haplotype in the dystrobrevin binding protein 1 (DTNBP1) gene previously associated with schizophrenia not only increases the susceptibility to psychotic illness but also to a more or less clinically specific form of psychotic illness.

Association between the 5q31.1 gene neurogenin1 and schizophrenia

Multiple lines of evidence suggest that schizophrenia results from aberrant neurodevelopment. The neurogenin1 gene (neurog1) consists of a single 1,666 bp exon that encodes a basic helix-loop-helix (bHLH) transcription factor that causes neuronal differentiation and induces cortical and glutamatergic differentiation programs. Because of its function and its location in 5q31.1, which has been linked to schizophrenia in multiple samples, we tested it for association with the disorder. We sequenced neurog1 in 25 affected subjects from the Irish Study of High-Density Schizophrenia Families. We observed a 5'-UTR SNP at position -60, already present in databases as rs8192558, and tested it along with rs2344485, rs8192559, and rs2344484. Narrow, intermediate, and broad diagnostic definitions were used. The major alleles of rs8192558 and rs2344484 were over-transmitted to affected subjects using both Pedigree Disequilibrium Test (PDT) (0.01 <or = P <or = 0.06) and FBAT (0.02 <or = P <or = 0.07). A haplotype consisting of the major alleles of all four SNPs was significantly over-transmitted in FBAT to the broad definition (P = 0.049), with trend significance to the narrow and intermediate definitions, and with trend significance in PDT. In confirmatory tests using 657 cases and 411 controls, this haplotype was slightly but not significantly over-represented in cases (81% vs. 77%, P = 0.21). These results, along with a priori evidence for the involvement of neurog1 in neurodevelopment, suggest that variants in neurog1 might have a small effect on susceptibility to schizophrenia. This gene should be tested in additional and larger samples.

Interleukin 3 and schizophrenia:the impact of sex and family history

Chromosome 5q21-33 has been implicated in harboring risk genes for schizophrenia. In this paper, we report evidence that multiple single nucleotide polymorphisms in and around interleukin 3 (IL3) are associated with the disease in the Irish Study of High-Density Schizophrenia Families (ISHDSF), the Irish Case-Control Study of Schizophrenia (ICCSS) and the Irish Trio Study of Schizophrenia (ITRIO). The associations are sex-specific and depend on the family history (FH) of schizophrenia. In all three samples, rs31400 shows female-specific and FH-dependent associations (P=0.0062, 0.0647 and 0.0284 for the ISHDSF, ICCSS and ITRIO, respectively). Several markers have similar associations in one or two of the three samples. In haplotype analyses, identical risk and protective haplotypes are identified in the ISHDSF and ITRIO samples in several multimarker combinations. For ICCSS, the same haplotypes are implicated; however, the risk haplotypes observed in the family samples become protective. Several significant markers, rs440970, rs31400 and rs2069803, are located in and around known estrogen response elements, promoter and enhancer of the IL3 gene. They may explain the sex-specific associations and be functional for the expression of IL3 gene.

Vitamin D receptor gene variants and esophageal adenocarcinoma risk:A population-based case-control study

Purpose: Polymorphisms in the vitamin D receptor (VDR) gene may be of etiological importance in determining cancer risk. The aim of this study was to assess the association between common VDR gene polymorphisms and esophageal adenocarcinoma (EAC) risk in an all-Ireland population-based case-control study. Methods: EAC cases and frequency-matched controls by age and gender recruited between March 2002 and December 2004 throughout Ireland were included. Participants were interviewed, and a blood sample collected for DNA extraction. Twenty-seven single nucleotide polymorphisms in the VDR gene were genotyped using Sequenom or TaqMan assays while the poly(A) microsatellite was genotyped by fluorescent fragment analysis. Unconditional logistic regression was applied to assess the association between VDR polymorphisms and EAC risk. Results: A total of 224 cases of EAC and 256 controls were involved in analyses. After adjustment for potential confounders, TT homozygotes at rs2238139 and rs2107301 had significantly reduced risks of EAC compared with CC homozygotes. In contrast, SS alleles of the poly(A) microsatellite had significantly elevated risks of EAC compared with SL/LL alleles. However, following permutation analyses to adjust for multiple comparisons, no significant associations were observed between any VDR gene polymorphism and EAC risk. Conclusions: VDR gene polymorphisms were not significantly associated with EAC development in this Irish population. Confirmation is required from larger studies. © Springer Science+Business Media, LLC 2011.

Genetic examination of SETD7 and SUV39H1/H2 methyltransferases and the risk of diabetes complications in patients with type 1 diabetes

Hyperglycemia plays a pivotal role in the development and progression of vascular complications, which are the major sources of morbidity and mortality in diabetes. Furthermore, these vascular complications often persist and progress despite improved glucose control, possibly as a result of prior episodes of hyperglycemia. Epigenetic modifications mediated by histone methyltransferases are associated with gene-activating events that promote enhanced expression of key proinflammatory molecules implicated in vascular injury. In this study, we investigated genetic polymorphisms of the SETD7, SUV39H1, and SUV39H2 methyltransferases as predictors of risk for micro- and macrovascular complications in type 1 diabetes.

The role of adiposity in cardiometabolic traits:a Mendelian randomization analysis

The association between adiposity and cardiometabolic traits is well known from epidemiological studies. Whilst the causal relationship is clear for some of these traits, for others it is not. We aimed to determine whether adiposity is causally related to various cardiometabolic traits using the Mendelian randomization approach.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Large-scale association analysis identifies new risk loci for coronary artery disease

Coronary artery disease (CAD) is the commonest cause of death. Here, we report an association analysis in 63,746 CAD cases and 130,681 controls identifying 15 loci reaching genome-wide significance, taking the number of susceptibility loci for CAD to 46, and a further 104 independent variants (r(2)

Genome-wide significant associations in schizophrenia to ITIH3/4, CACNA1C and SDCCAG8, and extensive replication of associations reported by the Schizophrenia PGC

The Schizophrenia Psychiatric Genome-Wide Association Study Consortium (PGC) highlighted 81 single-nucleotide polymorphisms (SNPs) with moderate evidence for association to schizophrenia. After follow-up in independent samples, seven loci attained genome-wide significance (GWS), but multi-locus tests suggested some SNPs that did not do so represented true associations. We tested 78 of the 81 SNPs in 2640 individuals with a clinical diagnosis of schizophrenia attending a clozapine clinic (CLOZUK), 2504 cases with a research diagnosis of bipolar disorder, and 2878 controls. In CLOZUK, we obtained significant replication to the PGC-associated allele for no fewer than 37 (47%) of the SNPs, including many prior GWS major histocompatibility complex (MHC) SNPs as well as 3/6 non-MHC SNPs for which we had data that were reported as GWS by the PGC. After combining the new schizophrenia data with those of the PGC, variants at three loci (ITIH3/4, CACNA1C and SDCCAG8) that had not previously been GWS in schizophrenia attained that level of support. In bipolar disorder, we also obtained significant evidence for association for 21% of the alleles that had been associated with schizophrenia in the PGC. Our study independently confirms association to three loci previously reported to be GWS in schizophrenia, and identifies the first GWS evidence in schizophrenia for a further three loci. Given the number of independent replications and the power of our sample, we estimate 98% (confidence interval (CI) 78-100%) of the original set of 78 SNPs represent true associations. We also provide strong evidence for overlap in genetic risk between schizophrenia and bipolar disorder.

Polymorphisms near TBX5 and GDF7 are associated with increased risk for Barrett's esophagus

BACKGROUND & AIMS: Barrett's esophagus (BE) increases the risk of esophageal adenocarcinoma (EAC). We found the risk to be BE has been associated with single nucleotide polymorphisms (SNPs) on chromosome 6p21 (within the HLA region) and on 16q23, where the closest protein-coding gene is FOXF1. Subsequently, the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON) identified risk loci for BE and esophageal adenocarcinoma near CRTC1 and BARX1, and within 100 kb of FOXP1. We aimed to identify further SNPs that increased BE risk and to validate previously reported associations.

METHODS: We performed a genome-wide association study (GWAS) to identify variants associated with BE and further analyzed promising variants identified by BEACON by genotyping 10,158 patients with BE and 21,062 controls.

RESULTS: We identified 2 SNPs not previously associated with BE: rs3072 (2p24.1; odds ratio [OR] = 1.14; 95% CI: 1.09-1.18; P = 1.8 × 10(-11)) and rs2701108 (12q24.21; OR = 0.90; 95% CI: 0.86-0.93; P = 7.5 × 10(-9)). The closest protein-coding genes were respectively GDF7 (rs3072), which encodes a ligand in the bone morphogenetic protein pathway, and TBX5 (rs2701108), which encodes a transcription factor that regulates esophageal and cardiac development. Our data also supported in BE cases 3 risk SNPs identified by BEACON (rs2687201, rs11789015, and rs10423674). Meta-analysis of all data identified another SNP associated with BE and esophageal adenocarcinoma: rs3784262, within ALDH1A2 (OR = 0.90; 95% CI: 0.87-0.93; P = 3.72 × 10(-9)).

CONCLUSIONS: We identified 2 loci associated with risk of BE and provided data to support a further locus. The genes we found to be associated with risk for BE encode transcription factors involved in thoracic, diaphragmatic, and esophageal development or proteins involved in the inflammatory response.