Judith Hendriks

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.

Effect of mammographic breast density on breast cancer screening performance: a study in Nijmegen, the Netherlands.

STUDY OBJECTIVE: To study the implications of breast density on mammographic screening performance. DESIGN: Screening outcomes of women with dense breast patterns were compared with those of women with lucent breast patterns (dense > 25% densities, lucent < or = 25% densities); the women were screened in different periods (before/after improvement of the mammographic technique in 1982). SETTING: Nijmegen, the Netherlands, 1977-1994. PARTICIPANTS: Between 1977 and 1994, 73,525 repeat screenings were performed in 19,152 participants (aged 50-69 years) in the Nijmegen breast cancer screening programme (repeat screenings were defined as mammographic examinations that were preceded by an examination in the previous screening round). Participants were screened biennially with mammography. There were 258 screen detected and 145 interval cancers. MAIN RESULTS: Before 1982 (rounds 2-4) the predictive value of a positive screening test (PV+) was lower in women with dense breasts than in those with lucent breasts (dense 29% v lucent 52%, p = 0.003). Also, the ratio of screen detected cancers to the total number of screen detected plus interval cancers (as a proxy for sensitivity) was lower in this group (based on a one year interval: dense 63% v lucent 92%, p = 0.001 and based on a two year interval: dense 41% v lucent 68%, p = 0.002). Moreover, the survival rate was less favourable for those with dense breasts (p = 0.07). In rounds 5-10, there were no important differences with respect to PV+ (dense 66% v lucent 62%, p = 0.57) or survival (p = 0.48). Moreover, sensitivity based on a one year interval was nearly as high in women with dense breasts as in those with lucent breasts (85% v 86%, p = 0.75). However, based on a two year interval sensitivity was lower (dense 59% v lucent 72%, p = 0.04). CONCLUSIONS: In the early screening years (rounds 2-4) high breast density had an unfavourable effect on screening performance. Nowadays, the situation has improved with respect to PV+, survival and detecting tumours in dense breasts with a lead time of up to one year, but little improvement has occurred in the detection of tumours with a lead time greater than one year.

Changes in mammographic breast density and concomitant changes in breast cancer risk

Among participants of the biennial Nijmegen breast cancer screening programme, we examined whether diminution of mammographic breast density lowered breast cancer risk. Post-menopausal breast cancer cases (n = 108), who had to have participated in all the five screening rounds prior to their diagnosis, were matched to 400 controls on year of birth and screening history. Controls had to be free of breast cancer at the time of the cases diagnosis. Changes in breast density were measured over a 10-year period, by a fully computerized method. Women in whom 5-25% or >25% of the breast was composed of fibro-glandular density showed a threefold increased 10-year risk compared to women with <5% density. In women with 5-25% density initially, we observed a trend of decreasing risk with diminishing density: when women with <5% density throughout the whole period formed the reference category, the odds ratio (OR) for those who decreased from 5-25% to <5% density was 1.9 [95% confidence interval (CI) = 0.6-6.1] in contrast to the OR of 5.7 (95% CI = 2.2-15.2) for those with persisting 5-25% density. In women who increased from 5-25% density to >25% density the OR was 6.9 (95% CI = 2.1-22.9). In women with >25% density initially, diminishing density was not clearly associated with lowering risk, which may be partly explained by the low number of women who decreased to <5% (n = 12). Due to the limited size of the study these results have to be interpreted with caution. Although the results are not conclusive, they could indicate a trend of decreasing risk with diminishing breast density. Should this effect be real, it may have great implications for the primary prevention of breast cancer or for the identification of high-risk groups who would benefit by more frequent screening. Therefore, large-scale, long-term follow-up studies on the effects of changes in breast density are needed.

Common Genetic Variations in CCK, Leptin, and Leptin Receptor Genes Are Associated With Specific Human Eating Patterns

Obesity has a heritable component; however, the heterogeneity of obesity complicates dissection of its genetic background. In this study, we therefore focused on eating patterns as specific traits within obesity. These traits have a heritable component; genes associated with a specific eating pattern have not yet been reported at the population level. In this study, we determined whether genetic variations in cholecystokinin (CCK) and leptin genes underlie specific eating patterns. We selected obese individuals showing extreme snacking behavior or use of excessive portion sizes from a large population-based sample (n = 17,357) from the Prospect-EPIC (European Prospective Study into Cancer and Nutrition) study. Using allele-specific PCRs, we tested several single nucleotide polymorphisms in the candidate genes and performed haplotype analysis. Obese carriers of common allelic variations in leptin or the leptin receptor gene had an increased risk to display extreme snacking behavior. In contrast, obese carriers of common allelic variations in CCK had an increased risk to eating increased meal sizes. In conclusion, we identified common allelic variants specifically associated with distinctly different eating patterns, namely extreme snacking behavior or excessive portion size.

The PIP5K2A and RGS4 genes are differentially associated with deficit and non-deficit schizophrenia

Several putative schizophrenia susceptibility genes have recently been reported, but it is not clear whether these genes are associated with schizophrenia in general or with specific disease subtypes. In a previous study, we found an association of the neuregulin 1 (NRG1) gene with non‐deficit schizophrenia only. We now report an association study of four schizophrenia candidate genes in patients with and without deficit schizophrenia, which is characterized by severe and enduring negative symptoms. Single‐nucleotide polymorphisms (SNPs) were genotyped in the DTNBP1 (dysbindin), G72/G30 and RGS4 genes, and the relatively unknown PIP5K2A gene, which is located in a region of linkage with both schizophrenia and bipolar disorder. The sample consisted of 273 Dutch schizophrenia patients, 146 of whom were diagnosed with deficit schizophrenia and 580 controls. The strongest evidence for association was found for the A‐allele of SNP rs10828317 in the PIP5K2A gene, which was associated with both clinical subtypes (P = 0.0004 in the entire group; non‐deficit P = 0.016, deficit P = 0.002). Interestingly, this SNP leads to a change in protein composition. In RGS4, the G‐allele of the previously reported SNP RGS4‐1 (single and as part of haplotypes with SNP RGS4‐18) was associated with non‐deficit schizophrenia (P = 0.03) but not with deficit schizophrenia (P = 0.79). SNPs in the DTNBP1 and G72/G30 genes were not significantly associated in any group. In conclusion, our data provide further evidence that specific genes may be involved in different schizophrenia subtypes and suggest that the PIP5K2A gene deserves further study as a general susceptibility gene for schizophrenia.

Short communication: Breast parenchymal patterns and their changes with age

In studies on mass screening, it has often been reported that tumours in breasts with dense parenchyma are difficult to detect and may have a more advanced stage at diagnosis. Shorter rescreening intervals have been suggested for these women but, before recommending such a strategy, it is important to investigate how often dense breast parenchyma (P2 and DY patterns according to Wolfe) is actually present in a screening population and to what extent these patterns change with age. The prevalence of dense breast parenchyma (P2 and DY) in our study population was 33% at first examination (n = 2581), which is fairly low compared with other screening populations. Its presence was strongly, inversely age-dependent. Breast patterns of 1177 women, aged 35-85 years, were followed for 12 years. In 39% (182/461) of the women with a P2 or DY pattern at their first examination, regression to a lucent pattern (N1 and P1) occurred over the years. The majority of these women were assumed to have reached menopause in the follow-up period. These findings support the hypothesis that the presence of dense breast parenchyma is related to the reproductive period and indicate that shortening the rescreening intervals would be most effective in pre-menopausal age groups.

Parity and mammographic breast density in relation to breast cancer risk: indication of interaction

We examined whether the harmful influence of nulliparity on breast cancer risk could be mediated by high mammographic density. Another possibility is that mammographic density and nulliparity act independently or perhaps synergistically on breast cancer risk. Our study population consisted of 129 cases and 517 controls who had been participants in the Nijmegen breast cancer screening programme for 10 years. Breast density was classified with a fully automated technique on digitized mammograms from the screening examination 10 years before diagnosis. Classification was based on the proportion of the breast that was composed of high density: < 5%, 5-25% or > 25%. Data on parity and potential confounders were obtained using a questionnaire, administered at the same examination. We found that nulliparae with low breast density (< 5%) were not at increased risk compared to parous women with low density: OR 1.1 (95% CI 0.2-5.8). Parous women with < 5% density formed the reference category throughout all analyses. The risks for parous women with 5-25% or > 25% density were 2.7 (95% CI 1.3-5.6) and 3.6 (95% CI 1.7-7.7) fold increased, respectively. However, when both factors were present (nulliparity and > or = 5% density), breast cancer risk was 7.1 times higher (95% CI 3.2-15.9). This could indicate that nulliparity and high breast density might work synergistically and that breast density is not just an explanatory factor in the influence of nulliparity on breast cancer risk. It is hypothesized that high breast density (reflecting fibro-glandular tissue with increased epithelial cell proliferation) is more susceptible to carcinogenic effects in the undifferentiated epithelial breast tissue of nulliparae than in the differentiated tissue of parous women. Since there were few data, no firm conclusions can be drawn. If these findings can be confirmed in a larger study population, however, they may have important implications for the prevention and early detection of breast cancer.

Interspecies Trait Genetics Reveals Association of Adcy8 with Mouse Avoidance Behavior and a Human Mood Disorder

BACKGROUND Identifying susceptibility genes for endophenotypes by studying analogous behaviors across species is an important strategy for understanding the pathophysiology underlying psychiatric disorders. This approach provides novel biological pathways plus validated animal models critical for selective drug development. One such endophenotype is avoidance behavior. METHODS In the present study, novel automated registration methods for longitudinal behavioral assessment in home cages are used to screen a panel of recently generated mouse chromosome substitution strains that are very powerful in quantitative trait loci (QTL) detection of complex traits. In this way, we identified chromosomes regulating avoidance behavior (increased sheltering preference) independent of motor activity levels (horizontal distance moved). Genetic information from the mouse QTL-interval was integrated with that from the homologous human linkage region for a mood disorder. RESULTS We genetically mapped a QTL for avoidance behavior on mouse chromosome 15, homologous with a human genome region (8q24) linked to bipolar disorder. Integrating the syntenic mouse QTL-interval with genotypes of 1868 BPD cases versus 14,311 control subjects revealed two associated genes (ADCY8 and KCNQ3). Adenylyl cyclase 8 (Adcy8) was differentially expressed in specific brain regions of mouse strains that differ in avoidance behavior levels. Finally, we showed that chronic infusion of the human mood stabilizer carbamazepine (that acts via adenylyl cyclase activity) significantly reduced mouse avoidance behavior, providing a further link between human mood disorders and this mouse home cage behavior. CONCLUSIONS Our data suggest that Adcy8 might encode a translational behavioral endophenotype of bipolar disorder.

A common variant in DRD3 receptor is associated with autism spectrum disorder.

BACKGROUND The presence of specific and common genetic etiologies for autism spectrum disorders (ASD) and attention-deficit/hyperactivity disorder (ADHD) was investigated for 132 candidate genes in a two-stage design-association study. METHODS 1,536 single nucleotide polymorphisms (SNPs) covering these candidate genes were tested in ASD (n = 144) and ADHD (n = 110) patients and control subjects (n = 404) from The Netherlands. A second stage was performed with those SNPs from Stage I reaching a significance threshold for association of p < .01 in an independent sample of ASD patients (n = 128) and controls (n = 124) from the United Kingdom and a Dutch ADHD (n = 150) and control (n = 149) sample. RESULTS No shared association was found between ASD and ADHD. However, in the first and second ASD samples and in a joint statistical analysis, a significant association between SNP rs167771 located in the DRD3 gene was found (joint analysis uncorrected: p = 3.11 x 10(-6); corrected for multiple testing and potential stratification: p = .00162). CONCLUSIONS The DRD3 gene is related to stereotyped behavior, liability to side effects of antipsychotic medication, and movement disorders and may therefore have important clinical implications for ASD.

High-resolution genetic mapping of mammalian motor activity levels in mice

The generation of motor activity levels is under tight neural control to execute essential behaviors, such as movement toward food or for social interaction. To identify novel neurobiological mechanisms underlying motor activity levels, we studied a panel of chromosome substitution (CS) strains derived from mice with high (C57BL/6J strain) or low motor activity levels (A/J strain) using automated home cage behavioral registration. In this study, we genetically mapped the expression of baseline motor activity levels (horizontal distance moved) to mouse chromosome 1. Further genetic mapping of this trait revealed an 8.3‐Mb quantitative trait locus (QTL) interval. This locus is distinct from the QTL interval for open‐field anxiety‐related motor behavior on this chromosome. By data mining, an existing phenotypic and genotypic data set of 2445 genetically heterogeneous mice (http://gscan.well.ox.ac.uk/), we confirmed linkage to the peak marker at 79 970 253 bp and refined the QTL to a 312‐kb interval containing a single gene (A830043J08Rik). Sequence analysis showed a nucleotide deletion in the 3′ untranslated region of the Riken gene. Genome‐wide microarray gene expression profiling in brains of discordant F2 individuals from CS strain 1 showed a significant upregulation of Epha4 in low‐active F2 individuals. Inclusion of a genetic marker for Epha4 confirmed that this gene is located outside of the QTL interval. Both Epha4 and A830043J08Rik are expressed in brain motor circuits, and similar to Epha4 mutants, we found linkage between reduced motor neurons number and A/J chromosome 1. Our findings provide a novel QTL and a potential downstream target underlying motor circuitry development and the expression of physical activity levels.