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.

Common variant at 16p11.2 conferring risk of psychosis

Epidemiological and genetic data support the notion that schizophrenia and bipolar disorder share genetic risk factors. In our previous genome-wide association (GWA) study, meta-analysis and follow-up (totaling as many as 18,206 cases and 42,536 controls), we identified four loci showing genome-wide significant association with schizophrenia. Here we consider a mixed schizophrenia and bipolar disorder (psychosis) phenotype (addition of 7,469 bipolar disorder cases, 1,535 schizophrenia cases, 333 other psychosis cases, 808 unaffected family members and 46,160 controls). Combined analysis reveals a novel variant at 16p11.2 showing genome-wide significant association (rs4583255[T], OR = 1.08, P = 6.6 × 10−11). The new variant is located within a 593 kb region that substantially increases risk of psychosis when duplicated. In line with the association of the duplication with reduced body mass index (BMI), rs4583255[T] is also associated with lower BMI (P = 0.0039 in the public GIANT consortium dataset; P = 0.00047 in 22,651 additional Icelanders).

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.

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.

Catechol-O-methyltransferase and the clinical features of psychosis

A functional polymorphism (Val-158-Met) at the Catechol-O-methyltransferase (COMT) locus has been identified as a potential etiological factor in schizophrenia. Yet the association has not been convincingly replicated across independent samples. We hypothesized that phenotypic heterogeneity might be diluting the COMT effect. To clarify the putative association, we performed an exploratory analysis to test for association between COMT and five psychosis symptom scales. These were derived through factor analysis of the Operational Criteria Checklist for Psychiatric Illness. Our sample was the Irish Study of High Density Schizophrenia Families, a large collection consisting of 268 multiplex families. This sample has previously shown a small but significant effect of the COMT Val allele in conferring risk for schizophrenia. We tested for preferential transmission of COMT alleles from parent to affected offspring (n = 749) for each of the five factor-derived scales (negative symptoms, delusions, hallucinations, mania, and depression). Significant overtransmission of the Val allele was found for mania (P <0.05) and depression (P = 0.01) scales. Examination of odds ratios (ORs) revealed a heterogeneous effect of COMT, whereby it had no effect on Negative Symptoms, but largest impact on Depression (OR = 1.4). These results suggest a modest affective vulnerability conferred by this allele in psychosis, but will require replication.

Haplotypes spanning SPEC2, PDZ-GEF2 and ACSL6 genes are associated with schizophrenia

Chromosome 5q22-33 is a region where studies have repeatedly found evidence for linkage to schizophrenia. In this report, we took a stepwise approach to systematically map this region in the Irish Study of High Density Schizophrenia Families (ISHDSF, 267 families, 1337 subjects) sample. We typed 289 SNPs in the critical interval of 8 million basepairs and found a 758 kb interval coding for the SPEC2/PDZ-GEF2/ACSL6 genes to be associated with the disease. Using sex and genotype-conditioned transmission disequilibrium test analyses, we found that 19 of the 24 typed markers were associated with the disease and the associations were sex-specific. We replicated these findings with an Irish case-control sample (657 cases and 414 controls), an Irish parent-proband trio sample (187 families, 564 subjects), a German nuclear family sample (211 families, 751 subjects) and a Pittsburgh nuclear family sample (247 families, 729 subjects). In all four samples, we replicated the sex-specific associations at the levels of both individual markers and haplotypes using sex- and genotype-conditioned analyses. Three risk haplotypes were identified in the five samples, and each haplotype was found in at least two samples. Consistent with the discovery of multiple estrogen-response elements in this region, our data showed that the impact of these haplotypes on risk for schizophrenia differed in males and females. From these data, we concluded that haplotypes underlying the SPEC2/PDZ-GEF2/ACSL6 region are associated with schizophrenia. However, due to the extended high LD in this region, we were unable to distinguish whether the association signals came from one or more of these genes.

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.

Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age

The genetic contribution to the variation in human lifespan is approximately 25%.  Despite the large number of identified disease-susceptibility loci, it is not known which loci influence population mortality.  We performed a genome-wide association meta-analysis of 7729 long-lived individuals of European descent (≥ 85 years) and 16121 younger controls (< 65 years) followed by replication in an additional set of 13060 long-lived individuals and 61156 controls. In addition, we performed a subset analysis in cases ≥ 90 years. We observed genome-wide significant association with longevity, as reflected by survival to ages beyond 90 years, at a novel locus, rs2149954, on chromosome 5q33.3 (OR = 1.10, P =1.74 x 10-8). We also confirmed association of rs4420638 on chromosome 19q13.32 (OR = 0.72, P = 3.40 x 10-36), representing the TOMM40/APOE/APOC1 locus. In a prospective meta-analysis (n = 34103) the minor allele of rs2149954 (T) on chromosome 5q33.3 associates with increased survival (HR = 0.95, P = 0.003). This allele has previously been reported to associate with low blood pressure in middle age. Interestingly, the minor allele (T) associates with decreased cardiovascular mortality risk, independent of blood pressure. We report on the first GWAS-identified longevity locus on chromosome 5q33.3 influencing survival in the general European population. The minor allele of this locus associates with low blood pressure in middle age, although the contribution of this allele to survival may be less dependent on blood pressure. Hence, the pleiotropic mechanisms by which this intragenic variation contributes to lifespan regulation have to be elucidated.

Uranium halide complexes in ionic liquids: an electrochemical and structural study

The electrochemistry of the salts, [emim](2)[UBr6] and [emim](2)[UO2Br4] ([emim] = 1-ethyl-3-methylimidazolium), has been investigated in both a basic and an acidic bromoaluminate(III) ionic liquid. In the basic ionic liquid, the hexabromo salt undergoes a one-electron reversible reduction process at a stationary glassy carbon disc electrode, while the tetrabromodioxo salt was reduced to a uranium(IV) species by an irreversible two-electron process with the simultaneous transfer of oxide to the ionic liquid. On the other hand, dissolution of either of the salts in an acidic bromoaluminate( III) ionic liquid resulted in the formation of the same electroactive species. The solid state structures of the uranium chloride salts, [emim](2)[UCl6] and [emim](2)[UO2Cl4], have previously been reported, but have now been re-evaluated using a new statistical model developed in our group, to determine the presence or absence of weak hydrogen bonding interactions in the crystalline state.

Transatlantic distribution of the Alaskan White River Ash

Volcanic ash layers preserved within the geologic record represent precise time markers that correlate disparate depositional environments and enable the investigation of synchronous and/or asynchronous behaviors in Earth system and archaeological sciences. However, it is generally assumed that only exceptionally powerful events, such as supereruptions (≥450 km3 of ejecta as dense-rock equivalent; recurrence interval of ∼105 yr), distribute ash broadly enough to have an impact on human society, or allow us to address geologic, climatic, and cultural questions on an intercontinental scale. Here we use geochemical, age, and morphological evidence to show that the Alaskan White River Ash (eastern lobe; A.D. 833–850) correlates to the “AD860B” ash (A.D. 846–848) found in Greenland and northern Europe. These occurrences represent the distribution of an ash over 7000 km, linking marine, terrestrial, and ice-core records. Our results indicate that tephra from more moderate-size eruptions, with recurrence intervals of ∼100 yr, can have substantially greater distributions than previously thought, with direct implications for volcanic dispersal studies, correlation of widely distributed proxy records, and volcanic hazard assessment.

A worldwide survey of chronic cough: A manifestation of enhanced somatosensory response:A manifestation of enhanced somatosensory response

Reports from individual centres suggest a preponderance of females with chronic cough. Females also have heightened cough reflex sensitivity. Here we have reviewed the age and sex of unselected referrals to 11 cough clinics. To investigate the cause of any observed sex dimorphism, functional magnetic resonance imaging of putative cough centres was analysed in normal volunteers. The demographic profile of consecutive patients presenting with chronic cough was evaluated. Cough challenge with capsaicin was undertaken in normal volunteers to construct a concentration-response curve. Subsequent functional magnetic resonance imaging during repeated inhalation of sub-tussive concentrations of capsaicin observed areas of activation within the brain and differences in the sexes identified. Of the 10 032 patients presenting with chronic cough, two-thirds (6591) were female (mean age 55 years). The patient profile was largely uniform across centres. The most common age for presentation was 60-69 years. The maximum tolerable dose of inhaled capsaicin was lower in females; however, a significantly greater activation of the somatosensory cortex was observed. Patients presenting with chronic cough from diverse racial and geographic backgrounds have a strikingly homogeneous demographic profile, suggesting a distinct clinical entity. The preponderance of females may be explained by sex-related differences in the central processing of cough sensation.

A solution concept for network games: The role of multilateral interactions

We propose an allocation rule that takes into account the importance of both players and their links and characterize it for a fixed network. Our characterization is along the lines of the characterization of the Position value for Network games by van den Nouweland and Slikker (2012). The allocation rule so defined admits multilateral interactions among the players through their links which distinguishes it from the other existing rules. Next, we extend our allocation rule to flexible networks à la Jackson (2005).

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.