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.

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.

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.

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.

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.

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.

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.

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.

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%.

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.

A genome-wide scan in type 2 diabetes mellitus provides independent replication of a susceptibility locus on 18p11 and suggests the existence of novel Loci on 2q12 and 19q13

A genome-wide scan was performed, using nonparametric linkage analyses, to find susceptibility loci for type 2 diabetes mellitus in the Dutch population. We studied 178 families from The Netherlands, who constituted 312 affected sibling pairs. The first stage of the genome scan consisted of 270 DNA markers, with an average intermarker spacing of 13 cM. Because obesity and type 2 diabetes mellitus are interrelated, the data set was stratified for the subphenotype body mass index, corrected for age and gender. This resulted in a suggestive maximum multipoint LOD score of 2.3 (single-point P value, 9.7 x 10(-4); genome-wide P value, 0.028) for the most obese 20% pedigrees of the data set, between marker loci D18S471 and D18S843. In the lowest 80% obese pedigrees, two interesting loci on chromosome 2 and 19 were found, with LOD scores of 1.5 and 1.3. We provide independent evidence that the chromosome 18p11 locus, reported earlier from a Finnish/Swedish population, is of definite interest for type 2 diabetes mellitus in connection with obesity. Subsequently, our results indicate that two novel loci may reside on chromosomes 2 and 19, with minor effects involved in the development of type 2 diabetes mellitus in the Dutch population.

Genetic variation in ICF syndrome: evidence for genetic heterogeneity

ICF syndrome is a rare autosomal recessive immunoglobulin deficiency, sometimes combined with defective cellular immunity. Other features that are frequently observed in ICF syndrome patients include facial dysmorphism, developmental delay, and recurrent infections. The most diagnostic feature of ICF syndrome is the branching of chromosomes 1, 9, and 16 due to pericentromeric instability. Positional candidate cloning recently discovered the de novo DNA methyltransferase 3B (DNMT3B) as the responsible gene by identifying seven different mutations in nine ICF patients. DNMT3B specifically methylates repeat sequences adjacent to the centromeres of chromosome 1, 9, and 16. Our panel of 14 ICF patients was subjected to mutation analysis in the DNMT3B gene. Mutations in DNMT3B were discovered in only nine of our 14 ICF patients. Moreover, two ICF patients from consanguineous families who did not show autozygosity (i.e. homozygosity by descent) for the DNMT3B locus did not reveal DNMT3B mutations, suggesting genetic heterogeneity for this disease. Mutation analysis revealed 11 different mutations, including seven novel ones: eight different missense mutations, two different nonsense mutations, and a splice-site mutation leading to the insertion of three aa's. The missense mutations occurred in or near the catalytic domain of DNMT3B protein, indicating a possible interference with the normal functioning of the enzyme. However, none of the ICF patients was homozygous for a nonsense allele, suggesting that absence of this enzyme is not compatible with life. Compound heterozygosity for a missense and a nonsense mutation did not seem to correlate with a more severe phenotype.

Localization of the ICF syndrome to chromosome 20 by homozygosity mapping

Immunodeficiency in association with centromere instability of chromosomes 1, 9, and 16 and facial anomalies (ICF syndrome) is a rare autosomal recessive disorder. ICF patients show marked hypomethylation of their DNA; undermethylation of classical satellites II and III is thought to be associated with the centromere instability. We used DNA from three consanguineous families with a total of four ICF patients and performed a total genome screen, to localize the ICF syndrome gene by homozygosity mapping. One chromosomal region (20q11-q13) was consistently found to be homozygous in ICF patients, whereas all healthy sibs showed a heterozygous pattern. Comparison of the regions of homozygosity in the four ICF patients localized the ICF locus to a 9-cM region between the markers D20S477 and D20S850. Analysis of more families will be required, to refine the map location further. Isolation of the gene associated with the ICF syndrome not only will give insight into the etiology of the ICF syndrome but will also broaden our understanding of DNA methylation processes.