Polygenic interactions with environmental adversity in the aetiology of major depressive disorder.

BACKGROUND: Major depressive disorder (MDD) is a common and disabling condition with well-established heritability and environmental risk factors. Gene-environment interaction studies in MDD have typically investigated candidate genes, though the disorder is known to be highly polygenic. This study aims to test for interaction between polygenic risk and stressful life events (SLEs) or childhood trauma (CT) in the aetiology of MDD. METHOD: The RADIANT UK sample consists of 1605 MDD cases and 1064 controls with SLE data, and a subset of 240 cases and 272 controls with CT data. Polygenic risk scores (PRS) were constructed using results from a mega-analysis on MDD by the Psychiatric Genomics Consortium. PRS and environmental factors were tested for association with case/control status and for interaction between them. RESULTS: PRS significantly predicted depression, explaining 1.1% of variance in phenotype (p = 1.9 × 10(-6)). SLEs and CT were also associated with MDD status (p = 2.19 × 10(-4) and p = 5.12 × 10(-20), respectively). No interactions were found between PRS and SLEs. Significant PRSxCT interactions were found (p = 0.002), but showed an inverse association with MDD status, as cases who experienced more severe CT tended to have a lower PRS than other cases or controls. This relationship between PRS and CT was not observed in independent replication samples. CONCLUSIONS: CT is a strong risk factor for MDD but may have greater effect in individuals with lower genetic liability for the disorder. Including environmental risk along with genetics is important in studying the aetiology of MDD and PRS provide a useful approach to investigating gene-environment interactions in complex traits.

The association between lower educational attainment and depression owing to shared genetic effects? Results in ~25,000 subjects.

An association between lower educational attainment (EA) and an increased risk for depression has been confirmed in various western countries. This study examines whether pleiotropic genetic effects contribute to this association. Therefore, data were analyzed from a total of 9662 major depressive disorder (MDD) cases and 14,949 controls (with no lifetime MDD diagnosis) from the Psychiatric Genomics Consortium with additional Dutch and Estonian data. The association of EA and MDD was assessed with logistic regression in 15,138 individuals indicating a significantly negative association in our sample with an odds ratio for MDD 0.78 (0.75-0.82) per standard deviation increase in EA. With data of 884,105 autosomal common single-nucleotide polymorphisms (SNPs), three methods were applied to test for pleiotropy between MDD and EA: (i) genetic profile risk scores (GPRS) derived from training data for EA (independent meta-analysis on ~120,000 subjects) and MDD (using a 10-fold leave-one-out procedure in the current sample), (ii) bivariate genomic-relationship-matrix restricted maximum likelihood (GREML) and (iii) SNP effect concordance analysis (SECA). With these methods, we found (i) that the EA-GPRS did not predict MDD status, and MDD-GPRS did not predict EA, (ii) a weak negative genetic correlation with bivariate GREML analyses, but this correlation was not consistently significant, (iii) no evidence for concordance of MDD and EA SNP effects with SECA analysis. To conclude, our study confirms an association of lower EA and MDD risk, but this association was not because of measurable pleiotropic genetic effects, which suggests that environmental factors could be involved, for example, socioeconomic status.

Type I interferon signaling genes in recurrent major depression: increased expression detected by whole-blood RNA sequencing.

A study of genome-wide gene expression in major depressive disorder (MDD) was undertaken in a large population-based sample to determine whether altered expression levels of genes and pathways could provide insights into biological mechanisms that are relevant to this disorder. Gene expression studies have the potential to detect changes that may be because of differences in common or rare genomic sequence variation, environmental factors or their interaction. We recruited a European ancestry sample of 463 individuals with recurrent MDD and 459 controls, obtained self-report and semi-structured interview data about psychiatric and medical history and other environmental variables, sequenced RNA from whole blood and genotyped a genome-wide panel of common single-nucleotide polymorphisms. We used analytical methods to identify MDD-related genes and pathways using all of these sources of information. In analyses of association between MDD and expression levels of 13 857 single autosomal genes, accounting for multiple technical, physiological and environmental covariates, a significant excess of low P-values was observed, but there was no significant single-gene association after genome-wide correction. Pathway-based analyses of expression data detected significant association of MDD with increased expression of genes in the interferon α/ß signaling pathway. This finding could not be explained by potentially confounding diseases and medications (including antidepressants) or by computationally estimated proportions of white blood cell types. Although cause-effect relationships cannot be determined from these data, the results support the hypothesis that altered immune signaling has a role in the pathogenesis, manifestation, and/or the persistence and progression of MDD.

Genotyping serotonin transporter polymorphisms 5-HTTLPR and rs25531 in European- and African-American subjects from the National Institute of Mental Health's Collaborative Center for Genomic Studies.

A number of studies have suggested DNA sequence variability in the serotonin transporter gene (SLC6A4) between European-American (EA) and African-American (AA) populations, which could be clinically important, given the central role SLC6A4 has in serotonin transmission. However, these studies have had relatively small samples, used self-reported measures of race, and have only tested the promoter-linked polymorphism 5-HTTLPR. Here we genotype 5-HTTLPR and rs25531, a neighboring functional polymorphism, in 954 AA and 2622EA subjects from a National Institute of Mental Health repository sample. Genotyping was performed using fragment analysis by capillary electrophoresis. AA, as compared with EA, groups had lower frequencies of the S allele (0.25 vs 0.43) and SS genotype (0.06 vs 0.19) at 5-HTTLPR, and higher rates of the G allele at rs25531 (0.21 vs 0.075). A rare xL variant at 5-HTTLPR was also more common among AAs (0.017 vs 0.008). When the polymorphisms were redefined into a high- and low-transcription haplotypes, the AA group showed significantly fewer low-transcription variants (χ(2)=4.8, P=0.03). No genotypes were associated with major depression, any anxiety disorder, or neuroticism in either EA or AA populations. This is the largest study to show SLC6A4 genotype differences between EA and AA populations, and the first to include rs25531. Lack of associations with clinical outcomes may reflect untested moderating environmental influences.

Schizophrenia susceptibility alleles are enriched for alleles that affect gene expression in adult human brain.

It is widely thought that alleles that influence susceptibility to common diseases, including schizophrenia, will frequently do so through effects on gene expression. As only a small proportion of the genetic variance for schizophrenia has been attributed to specific loci, this remains an unproven hypothesis. The International Schizophrenia Consortium (ISC) recently reported a substantial polygenic contribution to that disorder, and that schizophrenia risk alleles are enriched among single-nucleotide polymorphisms (SNPs) selected for marginal evidence for association (P<0.5) from genome-wide association studies (GWAS). It follows that if schizophrenia susceptibility alleles are enriched for those that affect gene expression, those marginally associated SNPs, which are also expression quantitative trait loci (eQTLs), should carry more true association signals compared with SNPs that are not marginally associated. To test this, we identified marginally associated (P<0.5) SNPs from two of the largest available schizophrenia GWAS data sets. We assigned eQTL status to those SNPs based upon an eQTL data set derived from adult human brain. Using the polygenic score method of analysis reported by the ISC, we observed and replicated the observation that higher probability cis-eQTLs predicted schizophrenia better than those with a lower probability for being a cis-eQTL. Our data support the hypothesis that alleles conferring risk of schizophrenia are enriched among those that affect gene expression. Moreover, our data show that notwithstanding the likely developmental origin of schizophrenia, studies of adult brain tissue can, in principle, allow relevant susceptibility eQTLs to be identified.

GWA study data mining and independent replication identify cardiomyopathy-associated 5 (CMYA5) as a risk gene for schizophrenia.

We conducted data-mining analyses using the Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) and molecular genetics of schizophrenia genome-wide association study supported by the genetic association information network (MGS-GAIN) schizophrenia data sets and performed bioinformatic prioritization for all the markers with P-values ≤0.05 in both data sets. In this process, we found that in the CMYA5 gene, there were two non-synonymous markers, rs3828611 and rs10043986, showing nominal significance in both the CATIE and MGS-GAIN samples. In a combined analysis of both the CATIE and MGS-GAIN samples, rs4704591 was identified as the most significant marker in the gene. Linkage disequilibrium analyses indicated that these markers were in low LD (3 828 611-rs10043986, r(2)=0.008; rs10043986-rs4704591, r(2)=0.204). In addition, CMYA5 was reported to be physically interacting with the DTNBP1 gene, a promising candidate for schizophrenia, suggesting that CMYA5 may be involved in the same biological pathway and process. On the basis of this information, we performed replication studies for these three single-nucleotide polymorphisms. The rs3828611 was found to have conflicting results in our Irish samples and was dropped out without further investigation. The other two markers were verified in 23 other independent data sets. In a meta-analysis of all 23 replication samples (family samples, 912 families with 4160 subjects; case-control samples, 11 380 cases and 15 021 controls), we found that both markers are significantly associated with schizophrenia (rs10043986, odds ratio (OR)=1.11, 95% confidence interval (CI)=1.04-1.18, P=8.2 × 10(-4) and rs4704591, OR=1.07, 95% CI=1.03-1.11, P=3.0 × 10(-4)). The results were also significant for the 22 Caucasian replication samples (rs10043986, OR=1.11, 95% CI=1.03-1.17, P=0.0026 and rs4704591, OR=1.07, 95% CI=1.02-1.11, P=0.0015). Furthermore, haplotype conditioned analyses indicated that the association signals observed at these two markers are independent. On the basis of these results, we concluded that CMYA5 is associated with schizophrenia and further investigation of the gene is warranted.

Genome-wide association study of recurrent early-onset major depressive disorder.

A genome-wide association study was carried out in 1020 case subjects with recurrent early-onset major depressive disorder (MDD) (onset before age 31) and 1636 control subjects screened to exclude lifetime MDD. Subjects were genotyped with the Affymetrix 6.0 platform. After extensive quality control procedures, 671 424 autosomal single nucleotide polymorphisms (SNPs) and 25 068 X chromosome SNPs with minor allele frequency greater than 1% were available for analysis. An additional 1 892 186 HapMap II SNPs were analyzed based on imputed genotypic data. Single-SNP logistic regression trend tests were computed, with correction for ancestry-informative principal component scores. No genome-wide significant evidence for association was observed, assuming that nominal P<5 × 10(-8) approximates a 5% genome-wide significance threshold. The strongest evidence for association was observed on chromosome 18q22.1 (rs17077540, P=1.83 × 10(-7)) in a region that has produced some evidence for linkage to bipolar-I or -II disorder in several studies, within an mRNA detected in human brain tissue (BC053410) and approximately 75 kb upstream of DSEL. Comparing these results with those of a meta-analysis of three MDD GWAS data sets reported in a companion article, we note that among the strongest signals observed in the GenRED sample, the meta-analysis provided the greatest support (although not at a genome-wide significant level) for association of MDD to SNPs within SP4, a brain-specific transcription factor. Larger samples will be required to confirm the hypothesis of association between MDD (and particularly the recurrent early-onset subtype) and common SNPs.

Novel loci for major depression identified by genome-wide association study of Sequenced Treatment Alternatives to Relieve Depression and meta-analysis of three studies.

We report a genome-wide association study (GWAS) of major depressive disorder (MDD) in 1221 cases from the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study and 1636 screened controls. No genome-wide evidence for association was detected. We also carried out a meta-analysis of three European-ancestry MDD GWAS data sets: STAR*D, Genetics of Recurrent Early-onset Depression and the publicly available Genetic Association Information Network-MDD data set. These data sets, totaling 3957 cases and 3428 controls, were genotyped using four different platforms (Affymetrix 6.0, 5.0 and 500 K, and Perlegen). For each of 2.4 million HapMap II single-nucleotide polymorphisms (SNPs), using genotyped data where available and imputed data otherwise, single-SNP association tests were carried out in each sample with correction for ancestry-informative principal components. The strongest evidence for association in the meta-analysis was observed for intronic SNPs in ATP6V1B2 (P=6.78 x 10⁻7), SP4 (P=7.68 x 10⁻7) and GRM7 (P=1.11 x 10⁻6). Additional exploratory analyses were carried out for a narrower phenotype (recurrent MDD with onset before age 31, N=2191 cases), and separately for males and females. Several of the best findings were supported primarily by evidence from narrow cases or from either males or females. On the basis of previous biological evidence, we consider GRM7 a strong MDD candidate gene. Larger samples will be required to determine whether any common SNPs are significantly associated with MDD.

Meta-analysis of 32 genome-wide linkage studies of schizophrenia.

A genome scan meta-analysis (GSMA) was carried out on 32 independent genome-wide linkage scan analyses that included 3255 pedigrees with 7413 genotyped cases affected with schizophrenia (SCZ) or related disorders. The primary GSMA divided the autosomes into 120 bins, rank-ordered the bins within each study according to the most positive linkage result in each bin, summed these ranks (weighted for study size) for each bin across studies and determined the empirical probability of a given summed rank (P(SR)) by simulation. Suggestive evidence for linkage was observed in two single bins, on chromosomes 5q (142-168 Mb) and 2q (103-134 Mb). Genome-wide evidence for linkage was detected on chromosome 2q (119-152 Mb) when bin boundaries were shifted to the middle of the previous bins. The primary analysis met empirical criteria for 'aggregate' genome-wide significance, indicating that some or all of 10 bins are likely to contain loci linked to SCZ, including regions of chromosomes 1, 2q, 3q, 4q, 5q, 8p and 10q. In a secondary analysis of 22 studies of European-ancestry samples, suggestive evidence for linkage was observed on chromosome 8p (16-33 Mb). Although the newer genome-wide association methodology has greater power to detect weak associations to single common DNA sequence variants, linkage analysis can detect diverse genetic effects that segregate in families, including multiple rare variants within one locus or several weakly associated loci in the same region. Therefore, the regions supported by this meta-analysis deserve close attention in future studies.

Genomewide linkage scan of schizophrenia in a large multicenter pedigree sample using single nucleotide polymorphisms.

A genomewide linkage scan was carried out in eight clinical samples of informative schizophrenia families. After all quality control checks, the analysis of 707 European-ancestry families included 1615 affected and 1602 unaffected genotyped individuals, and the analysis of all 807 families included 1900 affected and 1839 unaffected individuals. Multipoint linkage analysis with correction for marker-marker linkage disequilibrium was carried out with 5861 single nucleotide polymorphisms (SNPs; Illumina version 4.0 linkage map). Suggestive evidence for linkage (European families) was observed on chromosomes 8p21, 8q24.1, 9q34 and 12q24.1 in nonparametric and/or parametric analyses. In a logistic regression allele-sharing analysis of linkage allowing for intersite heterogeneity, genomewide significant evidence for linkage was observed on chromosome 10p12. Significant heterogeneity was also observed on chromosome 22q11.1. Evidence for linkage across family sets and analyses was most consistent on chromosome 8p21, with a one-LOD support interval that does not include the candidate gene NRG1, suggesting that one or more other susceptibility loci might exist in the region. In this era of genomewide association and deep resequencing studies, consensus linkage regions deserve continued attention, given that linkage signals can be produced by many types of genomic variation, including any combination of multiple common or rare SNPs or copy number variants in a region.

Analysis of 10 independent samples provides evidence for association between schizophrenia and a SNP flanking fibroblast growth factor receptor 2.

We and others have previously reported linkage to schizophrenia on chromosome 10q25-q26 but, to date, a susceptibility gene in the region has not been identified. We examined data from 3606 single-nucleotide polymorphisms (SNPs) mapping to 10q25-q26 that had been typed in a genome-wide association study (GWAS) of schizophrenia (479 UK cases/2937 controls). SNPs with P<0.01 (n=40) were genotyped in an additional 163 UK cases and those markers that remained nominally significant at P<0.01 (n=22) were genotyped in replication samples from Ireland, Germany and Bulgaria consisting of a total of 1664 cases with schizophrenia and 3541 controls. Only one SNP, rs17101921, was nominally significant after meta-analyses across the replication samples and this was genotyped in an additional six samples from the United States/Australia, Germany, China, Japan, Israel and Sweden (n=5142 cases/6561 controls). Across all replication samples, the allele at rs17101921 that was associated in the GWAS showed evidence for association independent of the original data (OR 1.17 (95% CI 1.06-1.29), P=0.0009). The SNP maps 85 kb from the nearest gene encoding fibroblast growth factor receptor 2 (FGFR2) making this a potential susceptibility gene for schizophrenia.

Association analyses of the neuregulin 1 gene with schizophrenia and manic psychosis in a Hispanic population.

OBJECTIVE: This study used the population of the Central Valley of Costa Rica (CVCR) and phenotyping strategies alternative to DSMIV classifications to investigate the association of neuregulin 1 with schizophrenia. METHOD: Using 134 family trios with a history of psychosis, we genotyped six of the seven markers originally identified to be associated with schizophrenia in Iceland. RESULTS: The neuregulin Icelandic haplotype was not associated with schizophrenia in the CVCR population. However, a novel haplotype was found to be overrepresented in subjects with functional psychosis (global P-value > 0.05). Stratification of the sample by history of mania suggests that this haplotype may be preferentially over-transmitted to persons with a history of manic psychosis. CONCLUSION: These results suggest that the neuregulin 1 gene is unlikely to play a major role in predisposing to schizophrenia in the CVCR. Further studies in the CVCR and other Latin American populations should be performed in order to corroborate these findings.

Multicenter linkage study of schizophrenia loci on chromosome 22q.

The hypothesis of the existence of one or more schizophrenia susceptibility loci on chromosome 22q is supported by reports of genetic linkage and association, meta-analyses of linkage, and the observation of elevated risk for psychosis in people with velocardiofacial syndrome, caused by 22q11 microdeletions. We tested this hypothesis by evaluating 10 microsatellite markers spanning 22q in a multicenter sample of 779 pedigrees. We also incorporated age at onset and sex into the analysis as covariates. No significant evidence for linkage to schizophrenia or for linkage associated with earlier age at onset, gender, or heterogeneity across sites was observed. We interpret these findings to mean that the population-wide effects of putative 22q schizophrenia susceptibility loci are too weak to detect with linkage analysis even in large samples.

Examination of IMPA1 and IMPA2 genes in manic-depressive patients: association between IMPA2 promoter polymorphisms and bipolar disorder.

Manic-depressive (bipolar) illness is a serious psychiatric disorder with a strong genetic predisposition. The disorder is likely to be multifactorial and etiologically complex, and the causes of genetic susceptibility have been difficult to unveil. Lithium therapy is a widely used pharmacological treatment of manic-depressive illness, which both stabilizes the ongoing episodes and prevents relapses. A putative target of lithium treatment has been the inhibition of the myo-inositol monophosphatase (IMPase) enzyme, which dephosphorylates myo-inositol monophosphate in the phosphatidylinositol signaling system. Two genes encoding human IMPases have so far been isolated, namely myo-inositol monophosphatase 1 (IMPA1) on chromosome 8q21.13-21.3 and myo-inositol monophosphatase 2 (IMPA2) on chromosome 18p11.2. In the present study, we have scanned for DNA variants in the human IMPA1 and IMPA2 genes in a pilot sample of Norwegian manic-depressive patients, followed by examination of selected polymorphisms and haplotypes in a family-based bipolar sample of Palestinian Arab proband-parent trios. Intriguingly, two frequent single-nucleotide polymorphisms (-461C>T and -207T>C) in the IMPA2 promoter sequence and their corresponding haplotypes showed transmission disequilibrium in the Palestinian Arab trios. No association was found between the IMPA1 polymorphisms and bipolar disorder, neither with respect to disease susceptibility nor with variation in lithium treatment response. The association between manic-depressive illness and IMPA2 variants supports several reports on the linkage of bipolar disorder to chromosome 18p11.2, and sustains the possible role of IMPA2 as a susceptibility gene in bipolar disorder.

Polymorphisms in the 5'-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression.

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 5' segment, spanning common DNA sequence variations, T-261G, A-161T, and -182INS/DEL-181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype -261G_-182INS-181_A-161 enhanced transcriptional activity 2.3-fold compared with the haplotype T-261_-182INS-181_A-161. Conversely, -161T reversed this, and the net effect when -261G and -161T were in the same haplotype (-261G_-182INS-181_-161T) was equivalent to the major haplotype (T-261_-182INS-181_A-161). Electrophoretic mobility shift experiments showed that -261G and -161T modify the binding of transcription factors (TFs): -261G generates a new AP2 binding site, while alleles A-161 and -161T exhibit different binding characteristics to AP1. T-261G and A-161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.

Simulation studies of detection of a complex disease in a partially isolated population.

Simulation studies were undertaken with POPGEN, a new population simulation program, to explore strategies for detecting loci underlying rare and common disorders in a small population that has been partially isolated for 10 generations. Haplotype-sharing analysis (HSA) and non-parametric linkage analysis (NPL) were applied to the simulated haplotype and pedigree data for 100 cases, 100 controls, and an average of 28 multiplex pedigrees from cases' families, for a 2-5 cM map of markers. When identity by descent (IBD) status was known (using unique founder marker allele designations assigned during simulation), a linkage disequilibrium (LD) signal could be detected under disease-generating models predicting relative risk to sibs of 11.8 (high-RR) or 2.67 (mod-RR). Detection was more difficult when marker alleles were down-coded to resemble microsatellites (heterozygosities 0.75-0.80). False-positive peaks on nondisease chromosomes were uncommon. NPL analysis was more powerful than HSA at this marker density using down-coded alleles and assuming availability of all affected relatives. LD mapping of common disorders is likely to require denser maps of highly polymorphic markers to approximate full IBD information. LD and linkage mapping provide independent information, and strategies that combine these two methods could be useful in studies of small isolated populations.

Genetic diversity of the human serotonin receptor 1B (HTR1B) gene.

We systematically and comprehensively investigated polymorphisms of the HTR1B gene as well as their linkage disequilibrium and ancestral relationships. We have detected the following polymorphisms in our sample via denaturing gradient gel electrophoresis, database comparisons, and/or previously published assays: G-511T, T-261G, -182INS/DEL-181, A-161T, C129T, T371G, T655C, C705T, G861C, A1099G, G1120A, and A1180G. The results of the intermarker analyses showed strong linkage disequilibrium between the C129T and the G861C polymorphisms and revealed four common haplotypes: ancestral (via chimpanzee comparisons), 129T/861C, -161T, and -182DEL-181. The results of association tests with schizophrenia were negative, although A-161T had a nominal P = 0.04 via ASPEX/sib_tdt. The expressed missense substitutions, Phe124Cys, Phe219Leu, Ile367Val, and Glu374Lys, could potentially affect ligand binding or interaction with G proteins and thus modify drug response in carriers of these variants. On average, the human cSNPs and differences among other primates clustered in the more thermodynamically unstable regions of the mRNA, which suggests that the evolutionary survival of nucleotide sequence variation may be influenced by the mRNA structure of this gene.

Haplotype sharing tests of linkage disequilibrium in a Hutterite asthma data set.

The Genetic Analysis Workshop 12 genome scan data set for "strict" asthma in a Hutterite population was analyzed using haplotype sharing analysis (HSA), which tests for differences in mean length of haplotype sharing around each marker for pairs of chromosomes in cases versus controls. The regions of chromosome 1 and 8 where evidence for linkage was observed in published analyses were negative by HSA. HSA yielded positive results on chromosomes 7, 12, 16, 18, and 21 (p = 0.003 on 21q). Although there are reports of support for linkage to asthma in some of these regions, it is not known whether any represent true positives. Further study is needed of the possible role of length-based measures of linkage disequilibrium in recent population isolates.