Systematic discovery of gene fusions in pediatric cancer by integrating RNA-seq and WGS

BACKGROUND

Gene fusions are important cancer drivers in pediatric cancer and their accurate detection is essential for diagnosis and treatment. Clinical decision-making requires high confidence and precision of detection. Recent developments show RNA sequencing (RNA-seq) is promising for genome-wide detection of fusion products but hindered by many false positives that require extensive manual curation and impede discovery of pathogenic fusions.

METHODS

We developed Fusion-sq to overcome existing disadvantages of detecting gene fusions. Fusion-sq integrates and "fuses" evidence from RNA-seq and whole genome sequencing (WGS) using intron-exon gene structure to identify tumor-specific protein coding gene fusions. Fusion-sq was then applied to the data generated from a pediatric pan-cancer cohort of 128 patients by WGS and RNA sequencing.

RESULTS

In a pediatric pan-cancer cohort of 128 patients, we identified 155 high confidence tumor-specific gene fusions and their underlying structural variants (SVs). This includes all clinically relevant fusions known to be present in this cohort (30 patients). Fusion-sq distinguishes healthy-occurring from tumor-specific fusions and resolves fusions in amplified regions and copy number unstable genomes. A high gene fusion burden is associated with copy number instability. We identified 27 potentially pathogenic fusions involving oncogenes or tumor-suppressor genes characterized by underlying SVs, in some cases leading to expression changes indicative of activating or disruptive effects.

CONCLUSIONS

Our results indicate how clinically relevant and potentially pathogenic gene fusions can be identified and their functional effects investigated by combining WGS and RNA-seq. Integrating RNA fusion predictions with underlying SVs advances fusion detection beyond extensive manual filtering. Taken together, we developed a method for identifying candidate gene fusions that is suitable for precision oncology applications. Our method provides multi-omics evidence for assessing the pathogenicity of tumor-specific gene fusions for future clinical decision making.

Molecular Characterization Reveals Subclasses of 1q Gain in Intermediate Risk Wilms Tumors

Simple Summary

Chromosomal alterations and other structural variants have been recurrently identified in Wilms tumors (WT) and are promising biomarkers for risk stratification. Chromosome 1q gain occurs in one in three WTs and is associated with poor prognosis, but its impact on tumor biology remains unknown. Here, we investigated the mutational mechanisms and functional effects of chromosomal alterations in WTs, and in particular 1q gain. We identified subgroups of tumors with typical activated biological processes: muscle differentiation, immune system, kidney development and proliferation. Combining these subgroups with genomic data showed that tumors with 1q gain occur in all subgroups and can be associated with different functional effects. Also, 1q gain tumors differ in mutational mechanisms and co-occurring tumor-specific mutations. In conclusion, we identified subgroups of tumors with 1q gain and therefore propose that incorporating expression data in risk stratification could improve the clinical utility of 1q gain.

Abstract

Chromosomal alterations have recurrently been identified in Wilms tumors (WTs) and some are associated with poor prognosis. Gain of 1q (1q+) is of special interest given its high prevalence and is currently actively studied for its prognostic value. However, the underlying mutational mechanisms and functional effects remain unknown. In a national unbiased cohort of 30 primary WTs, we integrated somatic SNVs, CNs and SVs with expression data and distinguished four clusters characterized by affected biological processes: muscle differentiation, immune system, kidney development and proliferation. Combined genome-wide CN and SV profiles showed that tumors profoundly differ in both their types of 1q+ and genomic stability and can be grouped into WTs with co-occurring 1p−/1q+, multiple chromosomal gains or CN neutral tumors. We identified 1q+ in eight tumors that differ in mutational mechanisms, subsequent rearrangements and genomic contexts. Moreover, 1q+ tumors were present in all four expression clusters reflecting activation of various biological processes, and individual tumors overexpress different genes on 1q. In conclusion, by integrating CNs, SVs and gene expression, we identified subgroups of 1q+ tumors reflecting differences in the functional effect of 1q gain, indicating that expression data is likely needed for further risk stratification of 1q+ WTs.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing

PURPOSE

Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events.

METHODS

We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing.

RESULTS

We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents.

CONCLUSION

We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

Amplicon-Based Targeted Next-Generation Sequencing of Formalin-Fixed, Paraffin-Embedded Tissue

Next-generation sequencing (NGS) is rapidly becoming the method of choice for mutation analysis in both research and diagnostics. The benefit of targeted NGS compared to whole-genome and whole-exome sequencing is that smaller amounts of input material can be used as well as qualitatively suboptimal tissue samples, like formalin-fixed, paraffin-embedded archival tissue.Here, we describe the protocol for targeted next-generation sequencing using the Ion Torrent PGM platform in combination with Ion Ampliseq NGS gene panels for formalin-fixed, paraffin-embedded tissues. Both the manual and the automated workflow are described as well as the bioinformatics for data analysis.

E-cadherin loss induces targetable autocrine activation of growth factor signalling in lobular breast cancer

Despite the fact that loss of E-cadherin is causal to the development and progression of invasive lobular carcinoma (ILC), options to treat this major breast cancer subtype are limited if tumours develop resistance to anti-oestrogen treatment regimens. This study aimed to identify clinically targetable pathways that are aberrantly active downstream of E-cadherin loss in ILC. Using a combination of reverse-phase protein array (RPPA) analyses, mRNA sequencing, conditioned medium growth assays and CRISPR/Cas9-based knock-out experiments, we demonstrate that E-cadherin loss causes increased responsiveness to autocrine growth factor receptor (GFR)-dependent activation of phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K)/Akt signalling. Autocrine activation of GFR signalling and its downstream PI3K/Akt hub was independent of oncogenic mutations in PIK3CA, AKT1 or PTEN. Analyses of human ILC samples confirmed growth factor production and pathway activity. Pharmacological inhibition of Akt using AZD5363 or MK2206 resulted in robust inhibition of cell growth and survival of ILC cells, and impeded tumour growth in a mouse ILC model. Because E-cadherin loss evokes hypersensitisation of PI3K/Akt activation independent of oncogenic mutations in this pathway, we propose clinical intervention of PI3K/Akt in ILC based on functional E-cadherin inactivation, irrespective of activating pathway mutations.

Genetic vulnerability to DUSP22 promoter hypermethylation is involved in the relation between in utero famine exposure and schizophrenia

Epigenetic changes may account for the doubled risk to develop schizophrenia in individuals exposed to famine in utero. We therefore investigated DNA methylation in a unique sample of patients and healthy individuals conceived during the great famine in China. Subsequently, we examined two case-control samples without famine exposure in whole blood and brain tissue. To shed light on the causality of the relation between famine exposure and DNA methylation, we exposed human fibroblasts to nutritional deprivation. In the famine-exposed schizophrenia patients, we found significant hypermethylation of the dual specificity phosphatase 22 (DUSP22) gene promoter (Chr6:291687-293285) (N = 153, p = 0.01). In this sample, DUSP22 methylation was also significantly higher in patients independent of famine exposure (p = 0.025), suggesting that hypermethylation of DUSP22 is also more generally involved in schizophrenia risk. Similarly, DUSP22 methylation was also higher in two separate case-control samples not exposed to famine using DNA from whole blood (N = 64, p = 0.03) and postmortem brains (N = 214, p = 0.007). DUSP22 methylation showed strong genetic regulation across chromosomes by a region on chromosome 16 which was consistent with new 3D genome interaction data. The presence of a direct link between famine and DUSP22 transcription was supported by data from cultured human fibroblasts that showed increased methylation (p = 0.048) and expression (p = 0.019) in response to nutritional deprivation (N = 10). These results highlight an epigenetic locus that is genetically regulated across chromosomes and that is involved in the response to early-life exposure to famine and that is relevant for a major psychiatric disorder.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

Maintenance treatment with capecitabine and bevacizumab versus observation in metastatic colorectal cancer: updated results and molecular subgroup analyses of the phase 3 CAIRO3 study

Background

The phase 3 CAIRO3 study showed that capecitabine plus bevacizumab (CAP-B) maintenance treatment after six cycles capecitabine, oxaliplatin, and bevacizumab (CAPOX-B) in metastatic colorectal cancer (mCRC) patients is effective, without compromising quality of life. In this post hoc analysis with updated follow-up and data regarding sidedness, we defined subgroups according to RAS/BRAF mutation status and mismatch repair (MMR) status, and investigated their influence on treatment efficacy.

Patients and methods

A total of 558 patients with previously untreated mCRC and stable disease or better after six cycles CAPOX-B induction treatment were randomised to either CAP-B maintenance treatment (n = 279) or observation (n = 279). Upon first progression, patients were to receive CAPOX-B reintroduction until second progression (PFS2, primary end point). We centrally assessed RAS/BRAF mutation status and MMR status, or used local results if central assessment was not possible. Intention-to-treat stratified Cox models adjusted for baseline covariables were used to examine whether treatment efficacy was modified by RAS/BRAF mutation status.

Results

RAS, BRAF mutations, and MMR deficiency were detected in 240/420 (58%), 36/381 (9%), and 4/279 (1%) patients, respectively. At a median follow-up of 87 months (IQR 69-97), all mutational subgroups showed significant improvement from maintenance treatment for the primary end point PFS2 [RAS/BRAF wild-type: hazard ratio (HR) 0.57 (95% CI 0.39-0.84); RAS-mutant: HR 0.74 (0.55-0.98); V600EBRAF-mutant: HR 0.28 (0.12-0.64)] and secondary end points, except for the RAS-mutant subgroup regarding overall survival. Adjustment for sidedness instead of primary tumour location yielded comparable results. Although right-sided tumours were associated with inferior prognosis, both patients with right- and left-sided tumours showed significant benefit from maintenance treatment.

Conclusions

CAP-B maintenance treatment after six cycles CAPOX-B is effective in first-line treatment of mCRC across all mutational subgroups. The benefit of maintenance treatment was most pronounced in patients with RAS/BRAF wild-type and V600EBRAF-mutant tumours.

ClinicalTrials.gov number

NCT00442637.

Mutational analysis using Sanger and next generation sequencing in sporadic spindle cell hemangiomas: A study of 19 cases

Spindle cell hemangioma (SCH) is a distinct vascular soft-tissue lesion characterized by cavernous blood vessels and a spindle cell component mainly occurring in the distal extremities of young adults. The majority of cases harbor heterozygous mutations in IDH1/2 sporadically or rarely in association with Maffucci syndrome. However, based on mosaicism and accordingly a low percentage of lesional cells harboring a mutant allele, detection can be challenging. We tested 19 sporadic SCHs by Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA), conventional next generation sequencing (NGS), and NGS using a single molecule molecular inversion probes (smMIP)-based library preparation to compare their diagnostic value. Out of 10 cases tested by Sanger sequencing and 2 analyzed using MLPA, 4 and 1, respectively, revealed a mutation in IDH1 (p.R132C). The 7 remaining negative cases and additional 6 cases were investigated using smMIP/NGS, showing hot spot mutations in IDH1 (p.R132C) (8 cases) and IDH2 (3 cases; twice p.R172S and once p.R172G, respectively). One case was negative. Owing to insufficient DNA quality and insufficient coverage, 2 cases were excluded. In total, in 16 out of 17 cases successfully tested, an IDH1/2 mutation was found. Given that IDH1/2 mutations were absent in 161 other vascular lesions tested by smMIP/NGS, the mutation can be considered as highly specific for SCH.

Targeted next-generation sequencing of commonly mutated genes in esophageal adenocarcinoma patients with long-term survival

Survival of patients with esophageal adenocarcinoma remains poor and individual differences in prognosis remain unexplained. This study investigated whether gene mutations can explain why patients with high-risk (pT3-4, pN+) esophageal adenocarcinoma survive past 5 years after esophagectomy. Six long-term survivors (LTS) (≥5 years survival without recurrence) and six short-term survivors (STS) (<2 years survival due to recurrence) who underwent resection without neoadjuvant therapy for high-risk esophageal adenocarcinoma were included. Targeted next-generation sequencing of 16 genes related to esophageal adenocarcinoma was performed. Mutations were compared between the LTS and STS and described in comparison with literature. A total of 48 mutations in 10 genes were identified. In the LTS, the median number of mutated genes per sample was 5 (range: 0-5) and the samples together harbored 22 mutations in 8 genes: APC (n = 1), CDH11 (n = 2), CDKN2A (n = 2), FAT4 (n = 5), KRAS (n = 1), PTPRD (n = 1), TLR4 (n = 8), and TP53 (n = 2). The median number of mutated genes per sample in the STS was 4 (range: 1-8) and in total 26 mutations were found in six genes: CDH11 (n = 5), FAT4 (n = 7), SMAD4 (n = 1), SMARCA4 (n = 1), TLR4 (n = 7), and TP53 (n = 5). CDH11, CDKN2A, FAT4, TLR4, and TP53 were mutated in at least 2 LTS or STS, exceeding mutation rates in literature. Mutations across the LTS and STS were found in 10 of the 16 genes. The results warrant future studies to investigate a larger range of genes in a larger sample size. This may result in a panel with prognostic genes, to predict individual prognosis and to select effective individualized therapy for patients with esophageal adenocarcinoma.

Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa

OBJECTIVE

The authors conducted a genome-wide association study of anorexia nervosa and calculated genetic correlations with a series of psychiatric, educational, and metabolic phenotypes.

METHOD

Following uniform quality control and imputation procedures using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, the authors performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression was used to calculate genome-wide common variant heritability (single-nucleotide polymorphism [SNP]-based heritability [h²_SNP]), partitioned heritability, and genetic correlations (r_g) between anorexia nervosa and 159 other phenotypes.

RESULTS

Results were obtained for 10,641,224 SNPs and insertion-deletion variants with minor allele frequencies >1% and imputation quality scores >0.6. The h²_SNP of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. The authors identified one genome-wide significant locus on chromosome 12 (rs4622308) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high-density lipoprotein cholesterol, and significant negative genetic correlations were observed between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes.

CONCLUSIONS

Anorexia nervosa is a complex heritable phenotype for which this study has uncovered the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. The study results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

Peripheral blood gene expression profiles linked to monoamine metabolite levels in cerebrospinal fluid

The blood-brain barrier separates circulating blood from the central nervous system (CNS). The scope of this barrier is not fully understood which limits our ability to relate biological measurements from peripheral to central phenotypes. For example, it is unknown to what extent gene expression levels in peripheral blood are reflective of CNS metabolism. In this study, we examine links between central monoamine metabolite levels and whole-blood gene expression to better understand the connection between peripheral systems and the CNS. To that end, we correlated the prime monoamine metabolites in cerebrospinal fluid (CSF) with whole-genome gene expression microarray data from blood (N=240 human subjects). We additionally applied gene-enrichment analysis and weighted gene co-expression network analyses (WGCNA) to identify modules of co-expressed genes in blood that may be involved with monoamine metabolite levels in CSF. Transcript levels of two genes were significantly associated with CSF serotonin metabolite levels after Bonferroni correction for multiple testing: THAP7 (P=2.8 × 10^-8, β=0.08) and DDX6 (P=2.9 × 10^-7, β=0.07). Differentially expressed genes were significantly enriched for genes expressed in the brain tissue (P=6.0 × 10^-52). WGCNA revealed significant correlations between serotonin metabolism and hub genes with known functions in serotonin metabolism, for example, HTR2A and COMT. We conclude that gene expression levels in whole blood are associated with monoamine metabolite levels in the human CSF. Our results, including the strong enrichment of brain-expressed genes, illustrate that gene expression profiles in peripheral blood can be relevant for quantitative metabolic phenotypes in the CNS.

Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase.

Clinical and genetic analysis of a family with two rare reflex epilepsies

PURPOSE

To determine clinical phenotypes, evolution and genetic background of a large family with a combination of two unusual forms of reflex epilepsies.

METHOD

Phenotyping was performed in eighteen family members (10 F, 8 M) including standardized EEG recordings with intermittent photic stimulation (IPS). Genetic analyses (linkage scans, Whole Exome Sequencing (WES) and Functional studies) were performed using photoparoxysmal EEG responses (PPRs) as affection status.

RESULTS

The proband suffered from speaking induced jaw-jerks and increasing limb jerks evoked by flickering sunlight since about 50 years of age. Three of her family members had the same phenotype. Generalized PPRs were found in seven members (six above 50 years of age) with myoclonus during the PPR. Evolution was typical: Sensitivity to lights with migraine-like complaints around adolescence, followed by jerks evoked by lights and spontaneously with dropping of objects, and strong increase of light sensitivity and onset of talking induced jaw jerks around 50 years. Linkage analysis showed suggestive evidence for linkage to four genomic regions. All photosensitive family members shared a heterozygous R129C mutation in the SCNM1 gene that regulates splicing of voltage gated ion channels. Mutation screening of 134 unrelated PPR patients and 95 healthy controls, did not replicate these findings.

CONCLUSION

This family presents a combination of two rare reflex epilepsies. Genetic analysis favors four genomic regions and points to a shared SCNM1 mutation that was not replicated in a general cohort of photosensitive subjects. Further genetic studies in families with similar combination of features are warranted.

Cross-disorder genome-wide analyses suggest a complex genetic relationship between Tourette's syndrome and OCD

OBJECTIVE

Obsessive-compulsive disorder (OCD) and Tourette's syndrome are highly heritable neurodevelopmental disorders that are thought to share genetic risk factors. However, the identification of definitive susceptibility genes for these etiologically complex disorders remains elusive. The authors report a combined genome-wide association study (GWAS) of Tourette's syndrome and OCD.

METHOD

The authors conducted a GWAS in 2,723 cases (1,310 with OCD, 834 with Tourette's syndrome, 579 with OCD plus Tourette's syndrome/chronic tics), 5,667 ancestry-matched controls, and 290 OCD parent-child trios. GWAS summary statistics were examined for enrichment of functional variants associated with gene expression levels in brain regions. Polygenic score analyses were conducted to investigate the genetic architecture within and across the two disorders.

RESULTS

Although no individual single-nucleotide polymorphisms (SNPs) achieved genome-wide significance, the GWAS signals were enriched for SNPs strongly associated with variations in brain gene expression levels (expression quantitative loci, or eQTLs), suggesting the presence of true functional variants that contribute to risk of these disorders. Polygenic score analyses identified a significant polygenic component for OCD (p=2×10(-4)), predicting 3.2% of the phenotypic variance in an independent data set. In contrast, Tourette's syndrome had a smaller, nonsignificant polygenic component, predicting only 0.6% of the phenotypic variance (p=0.06). No significant polygenic signal was detected across the two disorders, although the sample is likely underpowered to detect a modest shared signal. Furthermore, the OCD polygenic signal was significantly attenuated when cases with both OCD and co-occurring Tourette's syndrome/chronic tics were included in the analysis (p=0.01).

CONCLUSIONS

Previous work has shown that Tourette's syndrome and OCD have some degree of shared genetic variation. However, the data from this study suggest that there are also distinct components to the genetic architectures of these two disorders. Furthermore, OCD with co-occurring Tourette's syndrome/chronic tics may have different underlying genetic susceptibility compared with OCD alone.

Identification of schizophrenia-associated loci by combining DNA methylation and gene expression data from whole blood

Emerging evidence suggests that schizophrenia (SZ) susceptibility involves variation at genetic, epigenetic and transcriptome levels. We describe an integrated approach that leverages DNA methylation and gene expression data to prioritize genetic variation involved in disease. DNA methylation levels were obtained from whole blood of 260 SZ patients and 250 unaffected controls of which a subset with gene expression data was available. By assessing DNA methylation and gene expression in cases and controls, we identified 432 CpG sites with differential methylation levels that are associated with differential gene expression. We hypothesized that genetic factors involved in these methylation levels may be associated with the genetic risk of SZ susceptibility. To test this hypothesis, we used results from the Psychiatric Genomics Consortium SZ genome-wide association study (GWAS). We observe an enrichment of SZ-associated SNPs in the mQTLs of which the associated CpG site is also correlated with differential gene expression in SZ. While this enrichment was already apparent when using nominal significant thresholds, enrichment was even more pronounced when applying more stringent significance levels. One locus, previously identified as susceptibility locus in a SZ GWAS, involves SNP rs11191514:C>T, which regulates DNA methylation of calcium homeostasis modulator 1 that is also associated with differential gene expression in patients. Overall, our results suggest that epigenetic variation plays an important role in SZ susceptibility and that the integration of analyses of genetic, epigenetic and gene expression profiles may be a biologically meaningful approach for identifying disease susceptibility loci, even when using whole blood data in studies of brain-related disorders.

Investigation of the involvement of MIR185 and its target genes in the development of schizophrenia

BACKGROUND

Schizophrenia is a complex neuropsychiatric disorder of unclear etiology. The strongest known genetic risk factor is the 22q11.2 microdeletion. Research has yet to confirm which genes within the deletion region are implicated in schizophrenia. The minimal 1.5 megabase deletion contains MIR185, which encodes microRNA 185.

METHODS

We determined miR-185 expression in embryonic and adult mouse brains. Common and rare variants at this locus were then investigated using a human genetics approach. First, we performed gene-based analyses for MIR185 common variants and target genes using Psychiatric Genomics Consortium genome-wide association data. Second, MIR185 was resequenced in German patients (n = 1000) and controls (n = 500). We followed up promising variants by genotyping an additional European sample (patients, n = 3598; controls, n = 4082).

RESULTS

In situ hybridization in mice revealed miR-185 expression in brain regions implicated in schizophrenia. Gene-based tests revealed association between common variants in 3 MIR185 target genes (ATAT1, SH3PXD2A, NTRK3) and schizophrenia. Further analyses in mice revealed overlapping expression patterns for these target genes and miR-185. Resequencing identified 2 rare patient-specific novel variants flanking MIR185. However, follow-up genotyping provided no further evidence of their involvement in schizophrenia.

LIMITATIONS

Power to detect rare variant associations was limited.

CONCLUSION

Human genetic analyses generated no evidence of the involvement of MIR185 in schizophrenia. However, the expression patterns of miR-185 and its target genes in mice, and the genetic association results for the 3 target genes, suggest that further research into the involvement of miR-185 and its downstream pathways in schizophrenia is warranted.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Vitamin B-6 vitamers in human plasma and cerebrospinal fluid

BACKGROUND

Vitamin B-6 comprises a group of 6 interrelated vitamers and is essential for numerous physiologic processes, including brain functioning. Genetic disorders disrupting vitamin B-6 metabolism have severe clinical consequences.

OBJECTIVE

To adequately diagnose known and novel disorders in vitamin B-6 metabolism, a reference set is required containing information on all vitamin B-6 vitamers in plasma and cerebrospinal fluid (CSF).

DESIGN

Concentrations of vitamin B-6 vitamers in the plasma and CSF of 533 adult subjects were measured by ultra high-performance liquid chromatography-tandem mass spectrometry.

RESULTS

The relative vitamin B-6 vitamer composition of plasma [pyridoxal phosphate (PLP) > pyridoxic acid (PA) > pyridoxal] differed from that of CSF (pyridoxal > PLP > PA > pyridoxamine). Sex influenced vitamin B-6 vitamer concentrations in plasma and CSF and should therefore be taken into account when interpreting vitamin B-6 vitamer concentrations. The strict ratios and strong correlations between vitamin B-6 vitamers point to a tight regulation of vitamin B-6 vitamer concentrations in blood and CSF. Given the unique design of this study, with simultaneously withdrawn blood and CSF from a large number of subjects, reliable CSF:plasma ratios and correlations of vitamin B-6 vitamers could be established.

CONCLUSIONS

We provide an extensive reference set of vitamin B-6 vitamer concentrations in plasma and CSF. In addition to providing insight on the regulation of individual vitamers and their intercompartmental distribution, we anticipate that these data will prove to be a valuable reference set for the diagnosis and treatment of conditions associated with altered vitamin B-6 metabolism.

Seasonal changes in gene expression represent cell-type composition in whole blood

Seasonal patterns in behavior and biological parameters are widespread. Here, we examined seasonal changes in whole blood gene expression profiles of 233 healthy subjects. Using weighted gene co-expression network analysis, we identified three co-expression modules showing circannual patterns. Enrichment analysis suggested that this signal stems primarily from red blood cells and blood platelets. Indeed, a large clinical database with 51 142 observations of blood cell counts over 3 years confirmed a corresponding seasonal pattern of counts of red blood cells, reticulocytes and platelets. We found no direct evidence that these changes are linked to genes known to be key players in regulating immune function or circadian rhythm. It is likely, however, that these seasonal changes in cell counts and gene expression profiles in whole blood represent biological and clinical relevant phenomena. Moreover, our findings highlight possible confounding factors relevant to the study of gene expression profiles in subjects collected at geographical locations with disparaging seasonality patterns.

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 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 7469 bipolar disorder cases, 1535 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]; odds ratio=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 data set; P=0.00047 in 22 651 additional Icelanders).

Partitioning the heritability of Tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.

D-amino acid aberrations in cerebrospinal fluid and plasma of smokers

The glutamatergic neurotransmission system and the N-methyl-D-aspartate receptor (NMDAR) have been implicated in smoking and alcohol consumption behavior. Preclinical studies have demonstrated that nicotine and ethanol influence NMDAR functionality, which may have a role in tendencies to consume these substances. Nonetheless, little is known about concentrations of NMDAR coagonists in the cerebrospinal fluid (CSF) and plasma of individuals who smoke or consume alcohol. Glycine and L- and D-stereoisomers of alanine, serine, and proline were therefore measured using ultra-high-performance liquid chromatography-tandem mass spectrometry in 403 healthy subjects. Nicotine and alcohol consumption were quantified using questionnaires. Possible differences in NMDAR coagonist concentrations in plasma and CSF were investigated using ANCOVA with age, body mass index, and storage duration as covariates. The significance threshold was Bonferroni corrected (α=0.00625). Compared with non-smokers, smokers displayed lower levels of D-proline in plasma (p=0.0027, Cohen's d=-0.41) and D-proline in CSF (p=0.0026, Cohen's d=-0.43). D-Serine in CSF was higher in smokers than in non-smokers (p=0.0052, Cohen's d=0.41). After subdividing participants based on smoking quantity, dose-dependent decreases were demonstrated in smokers for D-proline in plasma (F=5.65, p=0.0039) and D-proline in CSF (F=5.20, p=0.0060). No differences in NMDAR coagonist levels between alcohol consumption groups were detected. To our knowledge, this is the first report to implicate D-amino acids in smoking behavior of humans. Whether such concentration differences lie at the root of or result from smoking habits may be addressed in prospective studies.

Genome-wide association study of obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a common, debilitating neuropsychiatric illness with complex genetic etiology. The International OCD Foundation Genetics Collaborative (IOCDF-GC) is a multi-national collaboration established to discover the genetic variation predisposing to OCD. A set of individuals affected with DSM-IV OCD, a subset of their parents, and unselected controls, were genotyped with several different Illumina SNP microarrays. After extensive data cleaning, 1465 cases, 5557 ancestry-matched controls and 400 complete trios remained, with a common set of 469,410 autosomal and 9657 X-chromosome single nucleotide polymorphisms (SNPs). Ancestry-stratified case-control association analyses were conducted for three genetically-defined subpopulations and combined in two meta-analyses, with and without the trio-based analysis. In the case-control analysis, the lowest two P-values were located within DLGAP1 (P=2.49 × 10(-6) and P=3.44 × 10(-6)), a member of the neuronal postsynaptic density complex. In the trio analysis, rs6131295, near BTBD3, exceeded the genome-wide significance threshold with a P-value=3.84 × 10(-8). However, when trios were meta-analyzed with the case-control samples, the P-value for this variant was 3.62 × 10(-5), losing genome-wide significance. Although no SNPs were identified to be associated with OCD at a genome-wide significant level in the combined trio-case-control sample, a significant enrichment of methylation QTLs (P<0.001) and frontal lobe expression quantitative trait loci (eQTLs) (P=0.001) was observed within the top-ranked SNPs (P<0.01) from the trio-case-control analysis, suggesting these top signals may have a broad role in gene expression in the brain, and possibly in the etiology of OCD.

BDNF Val66Met homozygosity does not influence plasma BDNF levels in healthy human subjects

A putative pathway by which the BDNF Val66Met polymorphism (rs6265) leads to aberrant phenotypes is its influence on plasma BDNF. Research into the impact of rs6265 on plasma BDNF has given rise to conflicting results. Moreover, most such studies have compared Met-carriers with Val-homozygous subjects. We therefore genotyped subjects from a population-based cohort (the Utrecht Health Project, N=2743) and assessed whether plasma BDNF differs between rs6265 homozygous groups. We maximized the number of Met-homozygous subjects in whom we measured plasma BDNF, resulting in plasma BDNF being available for 19 Met-homozygous and 42 matched Val-homozygous subjects. Mean concentrations (S.D.) were 1963.1 (750.1) and 2133.2 pg/ml (1164.3) for the Val/Val and Met/Met groups, respectively. Using ANOVA, no differences in plasma BDNF between the two groups were detected. In conclusion, these results add to a growing body of evidence indicating that allelic variation at rs6265 does not have medium to large effects on plasma BDNF concentrations.

Genetic analysis of DNA methylation and gene expression levels in whole blood of healthy human subjects

BACKGROUND

The predominant model for regulation of gene expression through DNA methylation is an inverse association in which increased methylation results in decreased gene expression levels. However, recent studies suggest that the relationship between genetic variation, DNA methylation and expression is more complex.

RESULTS

Systems genetic approaches for examining relationships between gene expression and methylation array data were used to find both negative and positive associations between these levels. A weighted correlation network analysis revealed that i) both transcriptome and methylome are organized in modules, ii) co-expression modules are generally not preserved in the methylation data and vice-versa, and iii) highly significant correlations exist between co-expression and co-methylation modules, suggesting the existence of factors that affect expression and methylation of different modules (i.e., trans effects at the level of modules). We observed that methylation probes associated with expression in cis were more likely to be located outside CpG islands, whereas specificity for CpG island shores was present when methylation, associated with expression, was under local genetic control. A structural equation model based analysis found strong support in particular for a traditional causal model in which gene expression is regulated by genetic variation via DNA methylation instead of gene expression affecting DNA methylation levels.

CONCLUSIONS

Our results provide new insights into the complex mechanisms between genetic markers, epigenetic mechanisms and gene expression. We find strong support for the classical model of genetic variants regulating methylation, which in turn regulates gene expression. Moreover we show that, although the methylation and expression modules differ, they are highly correlated.

Variance-Component Analysis of Obesity in Type 2 Diabetes Confirms Loci on Chromosomes 1q and 11q

To study genetic loci influencing obesity in nuclear families with type 2 diabetes, we performed a genome-wide screen with 325 microsatellite markers that had an average spacing of 11 cM and a mean heterozygosity of approximately 75% covering all 22 autosomes. Genotype data were obtained from 562 individuals from 178 families from the Breda Study Cohort. These families were determined to have at least two members with type 2 diabetes. As a measure of obesity, the BMI of each diabetes patient was determined. The genotypes were analyzed using variance components (VCs) analysis implemented in GENEHUNTER 2 to determine quantitative trait loci influencing BMI. The VC analysis revealed two genomic regions showing VC logarithm of odds (LOD) scores > or =1.0 on chromosome 1 and chromosome 11. The regions of interest on both chromosomes were further investigated by fine-mapping with additional markers, resulting in a VC LOD score of 1.5 on chromosome 1q and a VC LOD of 2.4 on chromosome 11q. The locus on chromosome 1 has been implicated previously in diabetes. The locus on chromosome 11 has been implicated previously in diabetes and obesity. Our study to determine linkage for BMI confirms the presence of quantitative trait loci influencing obesity in subjects with type 2 diabetes on chromosomes 1q31-q42 and 11q14-q24.

Common inversion polymorphism at 17q21.31 affects expression of multiple genes in tissue-specific manner

Background

Chromosome 17q21.31 contains a common inversion polymorphism of approximately 900 kb in populations with European ancestry. Two divergent MAPT haplotypes, H1 and H2 are described with distinct linkage disequilibrium patterns across the region reflecting the inversion status at this locus. The MAPT H1 haplotype has been associated with progressive supranuclear palsy, corticobasal degeneration, Parkinson’s disease and Alzheimer’s disease, while the H2 is linked to recurrent deletion events associated with the 17q21.31 microdeletion syndrome, a disease characterized by developmental delay and learning disability.

Results

In this study, we investigate the effect of the inversion on the expression of genes in the 17q21.31 region. We find the expression of several genes in and at the borders of the inversion to be affected; specific either to whole blood or different regions of the human brain. The H1 haplotype was found to be associated with an increased expression of LRRC37A4, PLEKH1M and MAPT. In contrast, a decreased expression of MGC57346, LRRC37A and CRHR1 was associated with H1.

Conclusions

Studies thus far have focused on the expression of MAPT in the inversion region. However, our results show that the inversion status affects expression of other genes in the 17q21.31 region as well. Given the link between the inversion status and different neurological diseases, these genes may also be involved in disease pathology, possibly in a tissue-specific manner.

Association between genetic variation in a region on chromosome 11 and schizophrenia in large samples from Europe

Recent molecular studies have implicated common alleles of small to moderate effect and rare alleles with larger effect sizes in the genetic architecture of schizophrenia (SCZ). It is expected that the reliable detection of risk variants with very small effect sizes can only be achieved through the recruitment of very large samples of patients and controls (that is tens of thousands), or large, potentially more homogeneous samples that have been recruited from confined geographical areas using identical diagnostic criteria. Applying the latter strategy, we performed a genome-wide association study (GWAS) of 1169 clinically well characterized and ethnically homogeneous SCZ patients from a confined area of Western Europe (464 from Germany, 705 from The Netherlands) and 3714 ethnically matched controls (1272 and 2442, respectively). In a subsequent follow-up study of our top GWAS results, we included an additional 2569 SCZ patients and 4088 controls (from Germany, The Netherlands and Denmark). Genetic variation in a region on chromosome 11 that contains the candidate genes AMBRA1, DGKZ, CHRM4 and MDK was significantly associated with SCZ in the combined sample (n=11 540; P=3.89 × 10(-9), odds ratio (OR)=1.25). This finding was replicated in 23 206 independent samples of European ancestry (P=0.0029, OR=1.11). In a subsequent imaging genetics study, healthy carriers of the risk allele exhibited altered activation in the cingulate cortex during a cognitive control task. The area of interest is a critical interface between emotion regulation and cognition that is structurally and functionally abnormal in SCZ and bipolar disorder.

A gene co-expression network in whole blood of schizophrenia patients is independent of antipsychotic-use and enriched for brain-expressed genes

Despite large-scale genome-wide association studies (GWAS), the underlying genes for schizophrenia are largely unknown. Additional approaches are therefore required to identify the genetic background of this disorder. Here we report findings from a large gene expression study in peripheral blood of schizophrenia patients and controls. We applied a systems biology approach to genome-wide expression data from whole blood of 92 medicated and 29 antipsychotic-free schizophrenia patients and 118 healthy controls. We show that gene expression profiling in whole blood can identify twelve large gene co-expression modules associated with schizophrenia. Several of these disease related modules are likely to reflect expression changes due to antipsychotic medication. However, two of the disease modules could be replicated in an independent second data set involving antipsychotic-free patients and controls. One of these robustly defined disease modules is significantly enriched with brain-expressed genes and with genetic variants that were implicated in a GWAS study, which could imply a causal role in schizophrenia etiology. The most highly connected intramodular hub gene in this module (ABCF1), is located in, and regulated by the major histocompatibility (MHC) complex, which is intriguing in light of the fact that common allelic variants from the MHC region have been implicated in schizophrenia. This suggests that the MHC increases schizophrenia susceptibility via altered gene expression of regulatory genes in this network.

Season of sampling and season of birth influence serotonin metabolite levels in human cerebrospinal fluid

Background

Animal studies have revealed seasonal patterns in cerebrospinal fluid (CSF) monoamine (MA) turnover. In humans, no study had systematically assessed seasonal patterns in CSF MA turnover in a large set of healthy adults.

Methodology/Principal Findings

Standardized amounts of CSF were prospectively collected from 223 healthy individuals undergoing spinal anesthesia for minor surgical procedures. The metabolites of serotonin (5-hydroxyindoleacetic acid, 5-HIAA), dopamine (homovanillic acid, HVA) and norepinephrine (3-methoxy-4-hydroxyphenylglycol, MPHG) were measured using high performance liquid chromatography (HPLC). Concentration measurements by sampling and birth dates were modeled using a non-linear quantile cosine function and locally weighted scatterplot smoothing (LOESS, span = 0.75). The cosine model showed a unimodal season of sampling 5-HIAA zenith in April and a nadir in October (p-value of the amplitude of the cosine = 0.00050), with predicted maximum (PC_max) and minimum (PC_min) concentrations of 173 and 108 nmol/L, respectively, implying a 60% increase from trough to peak. Season of birth showed a unimodal 5-HIAA zenith in May and a nadir in November (p = 0.00339; PC_max = 172 and PC_min = 126). The non-parametric LOESS showed a similar pattern to the cosine in both season of sampling and season of birth models, validating the cosine model. A final model including both sampling and birth months demonstrated that both sampling and birth seasons were independent predictors of 5-HIAA concentrations.

Conclusion

In subjects without mental illness, 5-HT turnover shows circannual variation by season of sampling as well as season of birth, with peaks in spring and troughs in fall.

Genome-wide analysis shows increased frequency of copy number variation deletions in Dutch schizophrenia patients

BACKGROUND

Since 2008, multiple studies have reported on copy number variations (CNVs) in schizophrenia. However, many regions are unique events with minimal overlap between studies. This makes it difficult to gain a comprehensive overview of all CNVs involved in the etiology of schizophrenia. We performed a systematic CNV study on the basis of a homogeneous genome-wide dataset aiming at all CNVs ≥ 50 kilobase pair. We complemented this analysis with a review of cytogenetic and chromosomal abnormalities for schizophrenia reported in the literature with the purpose of combining classical genetic findings and our current understanding of genomic variation.

METHODS

We investigated 834 Dutch schizophrenia patients and 672 Dutch control subjects. The CNVs were included if they were detected by QuantiSNP (http://www.well.ox.ac.uk/QuantiSNP/) as well as PennCNV (http://www.neurogenome.org/cnv/penncnv/) and contain known protein coding genes. The integrated identification of CNV regions and cytogenetic loci indicates regions of interest (cytogenetic regions of interest [CROIs]).

RESULTS

In total, 2437 CNVs were identified with an average number of 2.1 CNVs/subject for both cases and control subjects. We observed significantly more deletions but not duplications in schizophrenia cases versus control subjects. The CNVs identified coincide with loci previously reported in the literature, confirming well-established schizophrenia CROIs 1q42 and 22q11.2 as well as indicating a potentially novel CROI on chromosome 5q35.1.

CONCLUSIONS

Chromosomal deletions are more prevalent in schizophrenia patients than in healthy subjects and therefore confer a risk factor for pathogenicity. The combination of our CNV data with previously reported cytogenetic abnormalities in schizophrenia provides an overview of potentially interesting regions for positional candidate genes.

Common variants at VRK2 and TCF4 conferring risk of schizophrenia

Common sequence variants have recently joined rare structural polymorphisms as genetic factors with strong evidence for association with schizophrenia. Here we extend our previous genome-wide association study and meta-analysis (totalling 7 946 cases and 19 036 controls) by examining an expanded set of variants using an enlarged follow-up sample (up to 10 260 cases and 23 500 controls). In addition to previously reported alleles in the major histocompatibility complex region, near neurogranin (NRGN) and in an intron of transcription factor 4 (TCF4), we find two novel variants showing genome-wide significant association: rs2312147[C], upstream of vaccinia-related kinase 2 (VRK2) [odds ratio (OR) = 1.09, P = 1.9 × 10(-9)] and rs4309482[A], between coiled-coiled domain containing 68 (CCDC68) and TCF4, about 400 kb from the previously described risk allele, but not accounted for by its association (OR = 1.09, P = 7.8 × 10(-9)).

Copy number variations of chromosome 16p13.1 region associated with schizophrenia

Deletions and reciprocal duplications of the chromosome 16p13.1 region have recently been reported in several cases of autism and mental retardation (MR). As genomic copy number variants found in these two disorders may also associate with schizophrenia, we examined 4345 schizophrenia patients and 35,079 controls from 8 European populations for duplications and deletions at the 16p13.1 locus, using microarray data. We found a threefold excess of duplications and deletions in schizophrenia cases compared with controls, with duplications present in 0.30% of cases versus 0.09% of controls (P=0.007) and deletions in 0.12 % of cases and 0.04% of controls (P>0.05). The region can be divided into three intervals defined by flanking low copy repeats. Duplications spanning intervals I and II showed the most significant (P = 0.00010) association with schizophrenia. The age of onset in duplication and deletion carriers among cases ranged from 12 to 35 years, and the majority were males with a family history of psychiatric disorders. In a single Icelandic family, a duplication spanning intervals I and II was present in two cases of schizophrenia, and individual cases of alcoholism, attention deficit hyperactivity disorder and dyslexia. Candidate genes in the region include NTAN1 and NDE1. We conclude that duplications and perhaps also deletions of chromosome 16p13.1, previously reported to be associated with autism and MR, also confer risk of schizophrenia.

A large genome scan for rare CNVs in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease selectively affecting motor neurons in the brain and spinal cord. Recent genome-wide association studies (GWASs) have identified several common variants which increase disease susceptibility. In contrast, rare copy-number variants (CNVs), which have been associated with several neuropsychiatric traits, have not been studied for ALS in well-powered study populations. To examine the role of rare CNVs in ALS susceptibility, we conducted a CNV association study including over 19,000 individuals. In a genome-wide screen of 1875 cases and 8731 controls, we did not find evidence for a difference in global CNV burden between cases and controls. In our association analyses, we identified two loci that met our criteria for follow-up: the DPP6 locus (OR = 3.59, P = 6.6 × 10(-3)), which has already been implicated in ALS pathogenesis, and the 15q11.2 locus, containing NIPA1 (OR = 12.46, P = 9.3 × 10(-5)), the gene causing hereditary spastic paraparesis type 6 (HSP 6). We tested these loci in a replication cohort of 2559 cases and 5887 controls. Again, results were suggestive of association, but did not meet our criteria for independent replication: DPP6 locus: OR = 1.92, P = 0.097, pooled results: OR = 2.64, P = 1.4 × 10(-3); NIPA1: OR = 3.23, P = 0.041, pooled results: OR = 6.20, P = 2.2 × 10(-5)). Our results highlight DPP6 and NIPA1 as candidates for more in-depth studies. Unlike other complex neurological and psychiatric traits, rare CNVs with high effect size do not play a major role in ALS pathogenesis.

Association study of NRG1, DTNBP1, RGS4, G72/G30, and PIP5K2A with schizophrenia and symptom severity in a Hungarian sample

Genetic association studies have yielded extensive but frequently inconclusive data about genetic risk factors for schizophrenia. Clinical and genetic heterogeneity are possible factors explaining the inconsistent findings. The objective of this study was to test the association of commonly incriminated candidate genes with two clinically divergent subgroups, non-deficit (SZ-ND) and deficit-schizophrenia (SZ-D), and symptom severity, in order to test for replication of previously reported results. A homogeneous sample of 280 schizophrenia patients and 230 healthy controls of Hungarian, Caucasian descent were genotyped for polymorphisms in schizophrenia candidate genes NRG1, DTNBP1, RGS4, G72/G30, and PIP5K2A. Patients were divided into the diagnostic subgroups of SZ-ND and SZ-D using the Schedule for Deficit Syndrome (SDS), and assessed clinically by the Positive and Negative Symptom Scale (PANSS). SNP8NRG241930 in NRG1 and rs1011313 in DTNBP1 were associated with SZ-ND (P = 0.04 and 0.03, respectively). Polymorphisms in RGS4, G72/G30, and PIP5K2A were neither associated with SZ-ND nor with SZ-D. SNP8NRG241930 showed association with the PANSS cognitive and hostility/excitability factors, rs1011313 with the negative factor and SDS total score, and rs10917670 in RGS4 was associated with the depression factor. Although these results replicate earlier findings about the genetic background of SZ-ND and SZ-D only partially, our data seem to confirm previously reported association of NRG1 with schizophrenia without prominent negative symptoms. It was possible to detect associations of small-to-medium effect size between the investigated candidate genes and symptom severity. Such studies have the potential to unravel the possible connection between genetic and clinical heterogeneity in schizophrenia.

A large replication study and meta-analysis in European samples provides further support for association of AHI1 markers with schizophrenia

The Abelson helper integration site 1 (AHI1) gene locus on chromosome 6q23 is among a group of candidate loci for schizophrenia susceptibility that were initially identified by linkage followed by linkage disequilibrium mapping, and subsequent replication of the association in an independent sample. Here, we present results of a replication study of AHI1 locus markers, previously implicated in schizophrenia, in a large European sample (in total 3907 affected and 7429 controls). Furthermore, we perform a meta-analysis of the implicated markers in 4496 affected and 18,920 controls. Both the replication study of new samples and the meta-analysis show evidence for significant overrepresentation of all tested alleles in patients compared with controls (meta-analysis; P = 8.2 x 10(-5)-1.7 x 10(-3), common OR = 1.09-1.11). The region contains two genes, AHI1 and C6orf217, and both genes-as well as the neighbouring phosphodiesterase 7B (PDE7B)-may be considered candidates for involvement in the genetic aetiology of schizophrenia.

Gene expression profiling in C57BL/6J and A/J mouse inbred strains reveals gene networks specific for brain regions independent of genetic background

BACKGROUND

We performed gene expression profiling of the amygdala and hippocampus taken from inbred mouse strains C57BL/6J and A/J. The selected brain areas are implicated in neurobehavioral traits while these mouse strains are known to differ widely in behavior. Consequently, we hypothesized that comparing gene expression profiles for specific brain regions in these strains might provide insight into the molecular mechanisms of human neuropsychiatric traits. We performed a whole-genome gene expression experiment and applied a systems biology approach using weighted gene co-expression network analysis.

RESULTS

We were able to identify modules of co-expressed genes that distinguish a strain or brain region. Analysis of the networks that are most informative for hippocampus and amygdala revealed enrichment in neurologically, genetically and psychologically related pathways. Close examination of the strain-specific gene expression profiles, however, revealed no functional relevance but a significant enrichment of single nucleotide polymorphisms in the probe sequences used for array hybridization. This artifact was not observed for the modules of co-expressed genes that distinguish amygdala and hippocampus.

CONCLUSIONS

The brain-region specific modules were found to be independent of genetic background and are therefore likely to represent biologically relevant molecular networks that can be studied to complement our knowledge about pathways in neuropsychiatric disease.

Genome-wide association study identifies 19p13.3 (UNC13A) and 9p21.2 as susceptibility loci for sporadic amyotrophic lateral sclerosis

We conducted a genome-wide association study among 2,323 individuals with sporadic amyotrophic lateral sclerosis (ALS) and 9,013 control subjects and evaluated all SNPs with P < 1.0 x 10(-4) in a second, independent cohort of 2,532 affected individuals and 5,940 controls. Analysis of the genome-wide data revealed genome-wide significance for one SNP, rs12608932, with P = 1.30 x 10(-9). This SNP showed robust replication in the second cohort (P = 1.86 x 10(-6)), and a combined analysis over the two stages yielded P = 2.53 x 10(-14). The rs12608932 SNP is located at 19p13.3 and maps to a haplotype block within the boundaries of UNC13A, which regulates the release of neurotransmitters such as glutamate at neuromuscular synapses. Follow-up of additional SNPs showed genome-wide significance for two further SNPs (rs2814707, with P = 7.45 x 10(-9), and rs3849942, with P = 1.01 x 10(-8)) in the combined analysis of both stages. These SNPs are located at chromosome 9p21.2, in a linkage region for familial ALS with frontotemporal dementia found previously in several large pedigrees.

Common variants conferring risk of schizophrenia

Schizophrenia is a complex disorder, caused by both genetic and environmental factors and their interactions. Research on pathogenesis has traditionally focused on neurotransmitter systems in the brain, particularly those involving dopamine. Schizophrenia has been considered a separate disease for over a century, but in the absence of clear biological markers, diagnosis has historically been based on signs and symptoms. A fundamental message emerging from genome-wide association studies of copy number variations (CNVs) associated with the disease is that its genetic basis does not necessarily conform to classical nosological disease boundaries. Certain CNVs confer not only high relative risk of schizophrenia but also of other psychiatric disorders. The structural variations associated with schizophrenia can involve several genes and the phenotypic syndromes, or the 'genomic disorders', have not yet been characterized. Single nucleotide polymorphism (SNP)-based genome-wide association studies with the potential to implicate individual genes in complex diseases may reveal underlying biological pathways. Here we combined SNP data from several large genome-wide scans and followed up the most significant association signals. We found significant association with several markers spanning the major histocompatibility complex (MHC) region on chromosome 6p21.3-22.1, a marker located upstream of the neurogranin gene (NRGN) on 11q24.2 and a marker in intron four of transcription factor 4 (TCF4) on 18q21.2. Our findings implicating the MHC region are consistent with an immune component to schizophrenia risk, whereas the association with NRGN and TCF4 points to perturbation of pathways involved in brain development, memory and cognition.

Recurrent CNVs disrupt three candidate genes in schizophrenia patients

Schizophrenia is a severe psychiatric disease with complex etiology, affecting approximately 1% of the general population. Most genetics studies so far have focused on disease association with common genetic variation, such as single-nucleotide polymorphisms (SNPs), but it has recently become apparent that large-scale genomic copy-number variants (CNVs) are involved in disease development as well. To assess the role of rare CNVs in schizophrenia, we screened 54 patients with deficit schizophrenia using Affymetrix's GeneChip 250K SNP arrays. We identified 90 CNVs in total, 77 of which have been reported previously in unaffected control cohorts. Among the genes disrupted by the remaining rare CNVs are MYT1L, CTNND2, NRXN1, and ASTN2, genes that play an important role in neuronal functioning but--except for NRXN1--have not been associated with schizophrenia before. We studied the occurrence of CNVs at these four loci in an additional cohort of 752 patients and 706 normal controls from The Netherlands. We identified eight additional CNVs, of which the four that affect coding sequences were found only in the patient cohort. Our study supports a role for rare CNVs in schizophrenia susceptibility and identifies at least three candidate genes for this complex disorder.

Neutrophil recruitment and barrier impairment in celiac disease: a genomic study

BACKGROUND & AIMS

Celiac disease is an enteropathy featuring villous atrophy, crypt hyperplasia, and lymphocytosis. Tissue remodeling is driven by an inflammatory reaction to gluten in genetically susceptible individuals. The adaptive pathway is considered the major immune response but recent evidence has indicated the involvement of innate immunity as well. To assess the contribution of either immune response we performed global gene expression profiling of the regenerating mucosa.

METHODS

Microarray hybridizations were performed with biopsy samples from 13 untreated patients, 31 patients on a gluten-free diet in various stages of remission, and 21 controls. Additional data were generated using low-density array and conventional quantitative reverse-transcription polymerase chain reaction, and immunohistochemistry.

RESULTS

A total of 108 differentially expressed immune-related genes were identified (50 innate, 43 adaptive, 9 both innate/adaptive, and 6 immunoregulatory). Expression levels showed a gradual change as opposed to the discrete histological transitions. In addition to details provided on the adaptive and innate immune pathways used, we observed a chronic recruitment of activated neutrophils. Neutrophil involvement was unabated in otherwise completely normalized remission patients.

CONCLUSIONS

We observed a contribution of both the innate and adaptive immune response in celiac disease pathogenesis. The discrepancy between the histological classification and the observed incremental change in immune-gene expression may have consequences for current diagnostic inclusion criteria. Enhanced neutrophil infiltration in both active and remission patients points to a genetic impairment of the intestinal barrier that may contribute to the cause rather than the consequence of celiac disease.

Association of variants of transcription factor 7-like 2 (TCF7L2) with susceptibility to type 2 diabetes in the Dutch Breda cohort

AIM/HYPOTHESIS

A strong association between susceptibility to type 2 diabetes and common variants of transcription factor 7-like 2 (TCF7L2), encoding an enteroendocrine transcription factor involved in glucose homeostasis, has been reported in three different populations (Iceland, Denmark and USA) by Grant et al. We aimed to replicate these findings in a Dutch cohort.

METHODS

We analysed the genotypes of two intronic single nucleotide polymorphisms (SNPs) in TCF7L2 gene in 502 unrelated type 2 diabetes patients and in a set of healthy controls (n = 920). The two SNPs showed almost complete linkage disequilibrium (D' = 0.91).

RESULTS

We were able to replicate the previously reported association in our Breda cohort. The minor alleles of both variants were significantly over-represented in cases (odds ratio [OR] 1.29, 95% CI 1.09-1.52, [Formula: see text] for rs12255372; OR 1.41, 95% CI 1.19-1.66, [Formula: see text] for rs7903146). In addition, TCF7L2 haplotypes were analysed for association with the disease. The analysis of haplotypes did not reveal any strong association beyond that expected from analysing individual SNPs. The TT haplotype carrying the minor alleles was more frequent among cases (OR 1.38, [Formula: see text]).

CONCLUSIONS/INTERPRETATION

Our data strongly confirm that variants of the TCF7L2 gene contribute to the risk of type 2 diabetes. The population-attributable risk from this factor in the Dutch type 2 diabetes population is 10%.

Myosin IXB variant increases the risk of celiac disease and points toward a primary intestinal barrier defect

Celiac disease is probably the best-understood immune-related disorder. The disease presents in the small intestine and results from the interplay between multiple genes and gluten, the triggering environmental factor. Although HLA class II genes explain 40% of the heritable risk, non-HLA genes accounting for most of the familial clustering have not yet been identified. Here we report significant and replicable association (P = 2.1 x 10(-6)) to a common variant located in intron 28 of the gene myosin IXB (MYO9B), which encodes an unconventional myosin molecule that has a role in actin remodeling of epithelial enterocytes. Individuals homozygous with respect to the at-risk allele have a 2.3-times higher risk of celiac disease (P = 1.55 x 10(-5)). This result is suggestive of a primary impairment of the intestinal barrier in the etiology of celiac disease, which may explain why immunogenic gluten peptides are able to pass through the epithelial barrier.

Gene expression profiling and phenotype analyses ofS. cerevisiaein response to changing copper reveals six genes with new roles in copper and iron metabolism

Exhaustive microarray time course analyses of Saccharomyces cerevisiae during copper starvation and copper excess reveal new aspects of metal-induced gene regulation. Aside from identifying targets of established copper- and iron-responsive transcription factors, we find that genes encoding mitochondrial proteins are downregulated and that copper-independent iron transport genes are preferentially upregulated, both during prolonged copper deprivation. The experiments also suggest the presence of a small regulatory iron pool that links copper and iron responses. One hundred twenty-eight genes with putative roles in metal metabolism were further investigated by several systematic phenotype screens. Of the novel phenotypes uncovered, hsp12-Δ and arn1-Δ display increased sensitivity to copper, cyc1-Δ and crr1-Δ show resistance to high copper, vma13-Δ exhibits increased sensitivity to iron deprivation, and pep12-Δ results in reduced growth in high copper and low iron. Besides revealing new components of eukaryotic metal trafficking pathways, the results underscore the previously determined intimate links between iron and copper metabolism and mitochondrial and vacuolar function in metal trafficking. The analyses further suggest that copper starvation can specifically lead to downregulation of respiratory function to preserve iron and copper for other cellular processes.

Genome-wide screen in obese pedigrees with type 2 diabetes mellitus from a defined Dutch population

A genome scan was performed in obese type 2 diabetes mellitus pedigrees to identify susceptibility loci involved in obesity-driven type 2 diabetes mellitus. We studied the 20% most obese diabetes pedigrees from a confined Dutch population from around the town of Breda. Previously we, and others, have already shown that a susceptibility locus influencing obesity in diabetes may reside on chromosome 18p11. We now report evidence to also suggest linkage for type 2 diabetes in these obese pedigrees on chromosome regions 11p (genome-wide P-value </= 0.061) and 12q (genome-wide P-value </= 0.029), thereby confirming previous findings from corresponding regions. The linkage found in the Breda Cohort of type 2 diabetes patients is influenced by obesity. This supports the notion that a genetic predisposition to obesity is probably intertwined with a genetic predisposition to type 2 diabetes. Further efforts should address the question of how, on a genetic level, these two factors interact.