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. Molecular Psychiatry (2009) 14, 774-785; doi:10.1038/mp.2008.135; published online 30 December 2008

Genome-wide association identifies a common variant in the reelin gene that increases the risk of schizophrenia only in women

Sex differences in schizophrenia are well known, but their genetic basis has not been identified. We performed a genome-wide association scan for schizophrenia in an Ashkenazi Jewish population using DNA pooling. We found a female-specific association with rs7341475, a SNP in the fourth intron of the reelin ( RELN) gene (p = 2.9 x 10(-5) in women), with a significant gene-sex effect (p = 1.8 x 10(-4)). We studied rs7341475 in four additional populations, totaling 2,274 cases and 4,401 controls. A significant effect was observed only in women, replicating the initial result (p = 2.1 x 10(-3) in women; p = 4.2 x 10(-3) for gene-sex interaction). Based on all populations the estimated relative risk of women carrying the common genotype is 1.58 (p = 8.8 x 10(-7); p = 1.6 x 10(-5) for gene-sex interaction). The female-specific association between RELN and schizophrenia is one of the few examples of a replicated sex-specific genetic association in any disease.

Significant correlation in linkage signals from genome-wide scans of schizophrenia and schizotypy

Prior family and adoption studies have suggested a genetic relationship between schizophrenia and schizotypy. However, this has never been verified using linkage methods. We therefore attempted to test for a correlation in linkage signals from genome-wide scans of schizophrenia and schizotypy. The Irish study of high-density schizophrenia families comprises 270 families with at least two members with schizophrenia or poor-outcome schizoaffective disorder (n = 637). Non-psychotic relatives were assessed using the structured interview for schizotypy (n = 746). A 10-cM multipoint, non-parametric, autosomal genomewide scan of schizophrenia was performed in Merlin. A scan of a quantitative trait comprising ratings of DSM-III-R criteria for schizotypal personality disorder in non-psychotic relatives was also performed. Schizotypy logarithm of the odds (LOD) scores were regressed onto schizophrenia LOD scores at all loci, with adjustment for spatial autocorrelation. To assess empirical significance, this was also carried out using 1000 null scans of schizotypy. The number of jointly linked loci in the real data was compared to distribution of jointly linked loci in the null scans. No markers were suggestively linked to schizotypy based on strict Lander Kruglyak criteria. Schizotypy LODs predicted schizophrenia LODs above chance expectation genome wide (empirical P = 0.04). Two and four loci yielded nonparametric LOD (NPLs) > 1.0 and > 0.75, respectively, for both schizophrenia and schizotypy (genome-wide empirical P = 0.04 and 0.02, respectively). These results suggest that at least a subset of schizophrenia susceptibility genes also affects schizotypy in non-psychotic relatives. Power may therefore be increased in molecular genetic studies of schizophrenia if they incorporate measures of schizotypy in non-psychotic relatives.

Genome-wide association of early-onset myocardial infarction with single nucleotide polymorphisms and copy number variants

We conducted a genome-wide association study testing single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) for association with early-onset myocardial infarction in 2,967 cases and 3,075 controls. We carried out replication in an independent sample with an effective sample size of up to 19,492. SNPs at nine loci reached genome-wide significance: three are newly identified (21q22 near MRPS6-SLC5A3-KCNE2, 6p24 in PHACTR1 and 2q33 in WDR12) and six replicated prior observations1-4 (9p21, 1p13 near CELSR2-PSRC1-SORT1, 10q11 near CXCL12, 1q41 in MIA3, 19p13 near LDLR and 1p32 near PCSK9). We tested 554 common copy number polymorphisms (>1% allele frequency) and none met the pre-specified threshold for replication (P < 10-3). We identified 8,065 rare CNVs but did not detect a greater CNV burden in cases compared to controls, in genes compared to the genome as a whole, or at any individual locus. SNPs at nine loci were reproducibly associated with myocardial infarction, but tests of common and rare CNVs failed to identify additional associations with myocardial infarction risk.

Genome-wide association study identifies five new schizophrenia loci

We examined the role of common genetic variation in schizophrenia in a genome-wide association study of substantial size: a stage 1 discovery sample of 21,856 individuals of European ancestry and a stage 2 replication sample of 29,839 independent subjects. The combined stage 1 and 2 analysis yielded genome-wide significant associations with schizophrenia for seven loci, five of which are new (1p21.3, 2q32.3, 8p23.2, 8q21.3 and 10q24.32-q24.33) and two of which have been previously implicated (6p21.32-p22.1 and 18q21.2). The strongest new finding (P = 1.6 x 10(-11)) was with rs1625579 within an intron of a putative primary transcript for MIR137 (microRNA 137), a known regulator of neuronal development. Four other schizophrenia loci achieving genome-wide significance contain predicted targets of MIR137, suggesting MIR137-mediated dysregulation as a previously unknown etiologic mechanism in schizophrenia. In a joint analysis with a bipolar disorder sample (16,374 affected individuals and 14,044 controls), three loci reached genome-wide significance: CACNA1C (rs4765905, P = 7.0 x 10(-9)), ANK3 (rs10994359, P = 2.5 x 10(-8)) and the ITIH3-ITIH4 region (rs2239547, P = 7.8 x 10(-9)).

Characterisation of Genome-Wide PLZF/RARA Target Genes

The PLZF/RARA fusion protein generated by the t(11;17)(q23;q21) translocation in acute promyelocytic leukaemia (APL) is believed to act as an oncogenic transcriptional regulator recruiting epigenetic factors to genes important for its transforming potential. However, molecular mechanisms associated with PLZF/RARA-dependent leukaemogenesis still remain unclear. We searched for specific PLZF/RARA target genes by ChIP-on-chip in the haematopoietic cell line U937 conditionally expressing PLZF/RARA. By comparing bound regions found in U937 cells expressing endogenous PLZF with PLZF/RARA-induced U937 cells, we isolated specific PLZF/RARA target gene promoters. We next analysed gene expression profiles of our identified target genes in PLZF/RARA APL patients and analysed DNA sequences and epigenetic modification at PLZF/RARA binding sites. We identify 413 specific PLZF/RARA target genes including a number encoding transcription factors involved in the regulation of haematopoiesis. Among these genes, 22 were significantly down regulated in primary PLZF/RARA APL cells. In addition, repressed PLZF/RARA target genes were associated with increased levels of H3K27me3 and decreased levels of H3K9K14ac. Finally, sequence analysis of PLZF/RARA bound sequences reveals the presence of both consensus and degenerated RAREs as well as enrichment for tissue-specific transcription factor motifs, highlighting the complexity of targeting fusion protein to chromatin. Our study suggests that PLZF/RARA directly targets genes important for haematopoietic development and supports the notion that PLZF/RARA acts mainly as an epigenetic regulator of its direct target genes.

Genome-wide analysis of p63 binding sites identifies AP-2 factors as co-regulators of epidermal differentiation

The p63 transcription factor (TP63) is critical in development, growth and differentiation of stratifying epithelia. This is highlighted by the severity of congenital abnormalities caused by TP63 mutations in humans, the dramatic phenotypes in knockout mice and de-regulation of TP63 expression in neoplasia altering the tumour suppressive roles of the TP53 family. In order to define the normal role played by TP63 and provide the basis for better understanding how this network is perturbed in disease, we used chromatin immunoprecipitation combined with massively parallel sequencing (ChIP-seq) to identify >7500 high-confidence TP63-binding regions across the entire genome, in primary human neonatal foreskin keratinocytes (HFKs). Using integrative strategies, we demonstrate that only a subset of these sites are bound by TP53 in response to DNA damage. We identify a role for TP63 in transcriptional regulation of multiple genes genetically linked to cleft palate and identify AP-2alpha (TFAP2A) as a co-regulator of a subset of these genes. We further demonstrate that AP-2gamma (TFAP2C) can bind a subset of these regions and that acute depletion of either TFAP2A or TFAP2C alone is sufficient to reduce terminal differentiation of organotypic epidermal skin equivalents, indicating overlapping physiological functions with TP63.

SNP in the genome-wide association study hotspot on chromosome 9p21 confers susceptibility to diabetic nephropathy in type 1 diabetes

AIMS/HYPOTHESIS: Parental type 2 diabetes mellitus increases the risk of diabetic nephropathy in offspring with type 1 diabetes mellitus. Several single nucleotide polymorphisms (SNPs) that predispose to type 2 diabetes mellitus have recently been identified. It is, however, not known whether such SNPs also confer susceptibility to diabetic nephropathy in patients with type 1 diabetes mellitus. METHODS: We genotyped nine SNPs associated with type 2 diabetes mellitus in genome-wide association studies in the Finnish population, and tested for their association with diabetic nephropathy as well as with severe retinopathy and cardiovascular disease in 2,963 patients with type 1 diabetes mellitus. Replication of significant SNPs was sought in 2,980 patients from three other cohorts. RESULTS: In the discovery cohort, rs10811661 near gene CDKN2A/B was associated with diabetic nephropathy. The association remained after robust Bonferroni correction for the total number of tests performed in this study (OR 1.33 [95% CI 1.14, 1.56], p?=?0.00045, p (36tests)?=?0.016). In the meta-analysis, the combined result for diabetic nephropathy was significant, with a fixed effects p value of 0.011 (OR 1.15 [95% CI 1.02, 1.29]). The association was particularly strong when patients with end-stage renal disease were compared with controls (OR 1.35 [95% CI 1.13, 1.60], p?=?0.00038). The same SNP was also associated with severe retinopathy (OR 1.37 [95% CI 1.10, 1.69] p?=?0.0040), but the association did not remain after Bonferroni correction (p (36tests)?=?0.14). None of the other selected SNPs was associated with nephropathy, severe retinopathy or cardiovascular disease. CONCLUSIONS/INTERPRETATION: A SNP predisposing to type 2 diabetes mellitus, rs10811661 near CDKN2A/B, is associated with diabetic nephropathy in patients with type 1 diabetes mellitus.

Heritability and Genome-Wide Association Study to Assess Genetic Differences between Advanced Age-Related Macular Degeneration Subtypes

PURPOSE: To investigate whether the 2 subtypes of advanced age-related macular degeneration (AMD), choroidal neovascularization (CNV), and geographic atrophy (GA) segregate separately in families and to identify which genetic variants are associated with these 2 subtypes. DESIGN: Sibling correlation study and genome-wide association study (GWAS). PARTICIPANTS: For the sibling correlation study, 209 sibling pairs with advanced AMD were included. For the GWAS, 2594 participants with advanced AMD subtypes and 4134 controls were included. Replication cohorts included 5383 advanced AMD participants and 15 240 controls. METHODS: Participants had the AMD grade assigned based on fundus photography, examination, or both. To determine heritability of advanced AMD subtypes, a sibling correlation study was performed. For the GWAS, genome-wide genotyping was conducted and 6 036 699 single nucleotide polymorphisms (SNPs) were imputed. Then, the SNPs were analyzed with a generalized linear model controlling for genotyping platform and genetic ancestry. The most significant associations were evaluated in independent cohorts. MAIN OUTCOME MEASURES: Concordance of advanced AMD subtypes in sibling pairs and associations between SNPs with GA and CNV advanced AMD subtypes. RESULTS: The difference between the observed and expected proportion of siblings concordant for the same subtype of advanced AMD was different to a statistically significant degree (P = 4.2×10(-5)), meaning that in siblings of probands with CNV or GA, the same advanced subtype is more likely to develop. In the analysis comparing participants with CNV to those with GA, a statistically significant association was observed at the ARMS2/HTRA1 locus (rs10490924; odds ratio [OR], 1.47; P = 4.3×10(-9)), which was confirmed in the replication samples (OR, 1.38; P = 7.4×10(-14) for combined discovery and replication analysis). CONCLUSIONS: Whether CNV versus GA develops in a patient with AMD is determined in part by genetic variation. In this large GWAS meta-analysis and replication analysis, the ARMS2/HTRA1 locus confers increased risk for both advanced AMD subtypes, but imparts greater risk for CNV than for GA. This locus explains a small proportion of the excess sibling correlation for advanced AMD subtype. Other loci were detected with suggestive associations that differ for advanced AMD subtypes and deserve follow-up in additional studies. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

Genome-Wide Association Study of Multiplex Schizophrenia Pedigrees

OBJECTIVE The authors used a genome-wide association study (GWAS) of multiply affected families to investigate the association of schizophrenia to common single-nucleotide polymorphisms (SNPs) and rare copy number variants (CNVs). METHOD The family sample included 2,461 individuals from 631 pedigrees (581 in the primary European-ancestry analyses). Association was tested for single SNPs and genetic pathways. Polygenic scores based on family study results were used to predict case-control status in the Schizophrenia Psychiatric GWAS Consortium (PGC) data set, and consistency of direction of effect with the family study was determined for top SNPs in the PGC GWAS analysis. Within-family segregation was examined for schizophrenia-associated rare CNVs. RESULTS No genome-wide significant associations were observed for single SNPs or for pathways. PGC case and control subjects had significantly different genome-wide polygenic scores (computed by weighting their genotypes by log-odds ratios from the family study) (best p=10-17, explaining 0.4% of the variance). Family study and PGC analyses had consistent directions for 37 of the 58 independent best PGC SNPs (p=0.024). The overall frequency of CNVs in regions with reported associations with schizophrenia (chromosomes 1q21.1, 15q13.3, 16p11.2, and 22q11.2 and the neurexin-1 gene [NRXN1]) was similar to previous case-control studies. NRXN1 deletions and 16p11.2 duplications (both of which were transmitted from parents) and 22q11.2 deletions (de novo in four cases) did not segregate with schizophrenia in families. CONCLUSIONS Many common SNPs are likely to contribute to schizophrenia risk, with substantial overlap in genetic risk factors between multiply affected families and cases in large case-control studies. Our findings are consistent with a role for specific CNVs in disease pathogenesis, but the partial segregation of some CNVs with schizophrenia suggests that researchers should exercise caution in using them for predictive genetic testing until their effects in diverse populations have been fully studied.

Two non-synonymous markers in PTPN21, identified by genome-wide association study data-mining and replication, are associated with schizophrenia

We conducted data-mining analyses of genome wide association (GWA) studies of the CATIE and MGS-GAIN datasets, and found 13 markers in the two physically linked genes, PTPN21 and EML5, showing nominally significant association with schizophrenia. Linkage disequilibrium (LD) analysis indicated that all 7 markers from PTPN21 shared high LD (r(2)>0.8), including rs2274736 and rs2401751, the two non-synonymous markers with the most significant association signals (rs2401751, P=1.10 × 10(-3) and rs2274736, P=1.21 × 10(-3)). In a meta-analysis of all 13 replication datasets with a total of 13,940 subjects, we found that the two non-synonymous markers are significantly associated with schizophrenia (rs2274736, OR=0.92, 95% CI: 0.86-0.97, P=5.45 × 10(-3) and rs2401751, OR=0.92, 95% CI: 0.86-0.97, P=5.29 × 10(-3)). One SNP (rs7147796) in EML5 is also significantly associated with the disease (OR=1.08, 95% CI: 1.02-1.14, P=6.43 × 10(-3)). These 3 markers remain significant after Bonferroni correction. Furthermore, haplotype conditioned analyses indicated that the association signals observed between rs2274736/rs2401751 and rs7147796 are statistically independent. Given the results that 2 non-synonymous markers in PTPN21 are associated with schizophrenia, further investigation of this locus is warranted.

Genome-wide Association Study Implicates HLA-C*01:02 as a Risk Factor at the Major Histocompatibility Complex Locus in Schizophrenia

BACKGROUND: We performed a genome-wide association study (GWAS) to identify common risk variants for schizophrenia. METHODS: The discovery scan included 1606 patients and 1794 controls from Ireland, using 6,212,339 directly genotyped or imputed single nucleotide polymorphisms (SNPs). A subset of this sample (270 cases and 860 controls) was subsequently included in the Psychiatric GWAS Consortium-schizophrenia GWAS meta-analysis. RESULTS: One hundred eight SNPs were taken forward for replication in an independent sample of 13,195 cases and 31,021 control subjects. The most significant associations in discovery, corrected for genomic inflation, were (rs204999, p combined = 1.34 × 10(-9) and in combined samples (rs2523722 p combined = 2.88 × 10(-16)) mapped to the major histocompatibility complex (MHC) region. We imputed classical human leukocyte antigen (HLA) alleles at the locus; the most significant finding was with HLA-C*01:02. This association was distinct from the top SNP signal. The HLA alleles DRB1*03:01 and B*08:01 were protective, replicating a previous study. CONCLUSIONS: This study provides further support for involvement of MHC class I molecules in schizophrenia. We found evidence of association with previously reported risk alleles at the TCF4, VRK2, and ZNF804A loci.

Genome-wide profiling of methylation identifies novel targets with aberrant hyper-methylation and reduced expression in low-risk myelodysplastic syndromes

Gene expression profiling signatures may be used to classify the subtypes of Myelodysplastic syndrome (MDS) patients. However, there are few reports on the global methylation status in MDS. The integration of genome-wide epigenetic regulatory marks with gene expression levels would provide additional information regarding the biological differences between MDS and healthy controls. Gene expression and methylation status were measured using high-density microarrays. A total of 552 differentially methylated CpG loci were identified as being present in low-risk MDS; hypermethylated genes were more frequent than hypomethylated genes. In addition, mRNA expression profiling identified 1005 genes that significantly differed between low-risk MDS and the control group. Integrative analysis of the epigenetic and expression profiles revealed that 66.7% of the hypermethylated genes were underexpressed in low-risk MDS cases. Gene network analysis revealed molecular mechanisms associated with the low-risk MDS group, including altered apoptosis pathways. The two key apoptotic genes BCL2 and ETS1 were identified as silenced genes. In addition, the immune response and micro RNA biogenesis were affected by the hypermethylation and underexpression of IL27RA and DICER1. Our integrative analysis revealed that aberrant epigenetic regulation is a hallmark of low-risk MDS patients and could have a central role in these diseases.

Genome-wide association study to identify genetic determinants of severe asthma

BACKGROUND: The genetic basis for developing asthma has been extensively studied. However, association studies to date have mostly focused on mild to moderate disease and genetic risk factors for severe asthma remain unclear. OBJECTIVE: To identify common genetic variants affecting susceptibility to severe asthma. METHODS: A genome-wide association study was undertaken in 933 European ancestry individuals with severe asthma based on Global Initiative for Asthma (GINA) criteria 3 or above and 3346 clean controls. After standard quality control measures, the association of 480?889 genotyped single nucleotide polymorphisms (SNPs) was tested. To improve the resolution of the association signals identified, non-genotyped SNPs were imputed in these regions using a dense reference panel of SNP genotypes from the 1000 Genomes Project. Then replication of SNPs of interest was undertaken in a further 231 cases and 1345 controls and a meta-analysis was performed to combine the results across studies. RESULTS: An association was confirmed in subjects with severe asthma of loci previously identified for association with mild to moderate asthma. The strongest evidence was seen for the ORMDL3/GSDMB locus on chromosome 17q12-21 (rs4794820, p=1.03×10((-8)) following meta-analysis) meeting genome-wide significance. Strong evidence was also found for the IL1RL1/IL18R1 locus on 2q12 (rs9807989, p=5.59×10((-8)) following meta-analysis) just below this threshold. No novel loci for susceptibility to severe asthma met strict criteria for genome-wide significance. CONCLUSIONS: The largest genome-wide association study of severe asthma to date was carried out and strong evidence found for the association of two previously identified asthma susceptibility loci in patients with severe disease. A number of novel regions with suggestive evidence were also identified warranting further study.