Cohort Profile: The Nijmegen Biomedical Study (NBS)

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

BACKGROUND

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

METHODS

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis.

FINDINGS

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

INTERPRETATION

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

FUNDING

Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Genetic relationship between the immune system and autism

Autism spectrum disorder (ASD) is a common and complex neurodevelopmental condition. The pathophysiology of ASD is poorly defined; however, it includes a strong genetic component and there is increasing evidence to support a role of immune dysregulation. Nonetheless, it is unclear which immune phenotypes link to ASD through genetics. Hence, we investigated the genetic correlation between ASD and diverse classes of immune conditions and markers; and if these immune-related genetic factors link to specific autistic-like traits in the population. We estimated global and local genetic correlations between ASD (n = 55,420) and 11 immune phenotypes (n = 14,256-755,406) using genome-wide association study summary statistics. Subsequently, polygenic scores (PGS) for these immune phenotypes were calculated in a population-based sample (n = 2487) and associated to five autistic-like traits (i.e., attention to detail, childhood behaviour, imagination, rigidity, social skills), and a total autistic-like traits score. Sex-stratified PGS analyses were also performed. At the genome-wide level, ASD was positively correlated with allergic diseases (ALG), and negatively correlated with lymphocyte count, rheumatoid arthritis (RA), and systemic lupus erythematosus (SLE) (FDR-p = 0.01-0.02). At the local genetic level, ASD was correlated with RA, C-reactive protein, and granulocytes and lymphocyte counts (p = 5.8 × 10-6-0.002). In the general population sample, increased genetic liability for SLE, RA, ALG, and lymphocyte levels, captured by PGS, was associated with the total autistic score and with rigidity and childhood behaviour (FDR-p = 0.03). In conclusion, we demonstrated a genetic relationship between ASD and immunity that depends on the type of immune phenotype considered; some increase likelihood whereas others may potentially help build resilience. Also, this relationship may be restricted to specific genetic loci and link to specific autistic dimensions (e.g., rigidity).

Association of Risk Variants in the CFH Gene With Elevated Levels of Coagulation and Complement Factors in Idiopathic Multifocal Choroiditis

Importance

Idiopathic multifocal choroiditis (MFC) is poorly understood, thereby hindering optimal treatment and monitoring of patients.

Objective

To identify the genes and pathways associated with idiopathic MFC.

Design, Setting, and Participants

This was a case-control genome-wide association study (GWAS) and protein study of blood plasma samples conducted from March 2006 to February 2022. This was a multicenter study involving 6 Dutch universities. Participants were grouped into 2 cohorts: cohort 1 consisted of Dutch patients with idiopathic MFC and controls, and cohort 2 consisted of patients with MFC and controls. Plasma samples from patients with idiopathic MFC who had not received treatment were subjected to targeted proteomics. Idiopathic MFC was diagnosed according to the Standardization of Uveitis Nomenclature (SUN) Working Group guidelines for punctate inner choroidopathy and multifocal choroiditis with panuveitis. Data were analyzed from July 2021 to October 2022.

Main outcomes and measures

Genetic variants associated with idiopathic MFC and risk variants associated with plasma protein concentrations in patients.

Results

This study included a total of 4437 participants in cohort 1 (170 [3.8%] Dutch patients with idiopathic MFC and 4267 [96.2%] controls; mean [SD] age, 55 [18] years; 2443 female [55%]) and 1344 participants in cohort 2 (52 [3.9%] patients with MFC and 1292 [96.1%] controls; 737 male [55%]). The primary GWAS association mapped to the CFH gene with genome-wide significance (lead variant the A allele of rs7535263; odds ratio [OR], 0.52; 95% CI, 0.41-0.64; P = 9.3 × 10-9). There was no genome-wide significant association with classical human leukocyte antigen (HLA) alleles (lead classical allele, HLA-A*31:01; P = .002). The association with rs7535263 showed consistent direction of effect in an independent cohort of 52 cases and 1292 control samples (combined meta-analysis OR, 0.58; 95% CI, 0.38-0.77; P = 3.0 × 10-8). In proteomic analysis of 87 patients, the risk allele G of rs7535263 in the CFH gene was strongly associated with increased plasma concentrations of factor H-related (FHR) proteins (eg, FHR-2, likelihood ratio test, adjusted P = 1.1 × 10-3) and proteins involved in platelet activation and the complement cascade.

Conclusions and relevance

Results suggest that CFH gene variants increase systemic concentrations of key factors of the complement and coagulation cascades, thereby conferring susceptibility to idiopathic MFC. These findings suggest that the complement and coagulation pathways may be key targets for the treatment of idiopathic MFC.

Novel (sulfated) thyroid hormone transporters in the solute carrier 22 family

Objective

Thyroid hormone (TH) transport represents a critical first step in governing intracellular TH regulation. It is still unknown whether the full repertoire of TH transporters has been identified. Members of the solute carrier (SLC) 22 family have substrates in common with the known TH transporters of the organic anion-transporting peptide family. Therefore, we screened the SLC22 family for TH transporters.

Methods

Uptake of 1 nM of iodothyronines or sulfated iodothyronines in COS1 cells expressing SLC22 proteins was performed.

Results

We first tested 25 mouse (m) SLC22 proteins for TH uptake and found that the majority of the organic anion transporter (OAT) clade were capable of 3,3',5-triiodothyronine and/or thyroxine (T4) transport. Based on phylogenetic tree analysis of the mouse and human (h) SLC22 family, we selected eight hSLC22s that grouped with the newly identified mouse TH transporters. Of these, four tested positive for uptake of one or more substrates, particularly hSLC22A11 showed robust (3-fold over control) uptake of T4. Uptake of sulfated iodothyronines was strongly (up to 17-fold) induced by some SLC22s, most notably SLC22A8, hSLC22A9, mSLC22A27 and mSLC22A29. Finally, the zebrafish orthologues of SLC22A6/8 drOatx and drSlc22a6l also transported almost all (sulfated) iodothyronines tested. The OAT inhibitors lesinurad and probenecid inhibited most SLC22 proteins.

Conclusions

Our results demonstrated that members of the OAT clade of the SLC22 family constitute a novel, evolutionary conserved group of transporters for (sulfated) iodothyronines. Future studies should reveal the relevance of these transporters in TH homeostasis and physiology.

Prevalence and Therapeutic Implications of Clonal Hematopoiesis of Indeterminate Potential in Young Patients With Stroke

BACKGROUND

Undetermined stroke etiology hampers optimal secondary prevention in a large proportion of young patients. We explored whether genetic screening for clonal hematopoiesis of indetermined potential (CHIP), a novel risk factor for stroke, could identify patients with myeloid precursor lesions or covert myeloid neoplasm requiring specific treatment.

METHODS

We performed targeted sequencing on 56 genes recurrently mutated in hematologic neoplasms in a prospective cohort of patients with acute brain ischemia between 18 and 60 years. CHIP prevalence was compared with age-matched healthy controls from the Nijmegen Biomedical Study (n=1604) and the UK Biobank (n=101 678). Patients with suspicion of high-risk CHIP or myeloid neoplasm were invited for further hematologic evaluation.

RESULTS

We included 248 consecutive patients (39% women) of whom 176 (71%) had cryptogenic stroke etiology. Fifty-one (21%) patients had CHIP, 3-fold more than in the general population (7.7% versus 2.6% for the Nijmegen Biomedical Study and 11.9% versus 4.1% for UK Biobank; P<0.001 for both). Patients with CHIP were older (median [interquartile range], 53 [50-59] versus 51 [41-56] years; P<0.001), had higher carotid intima-media thickness (0.68 [0.58-0.80] versus 0.59 [0.51-0.73] mm; P=0.009), and had higher burden of atherosclerosis (29.4% versus 16.7%; P=0.04). We invited 11 patients (4.4%) for further hematologic assessment, which in 7 led to the diagnosis of high-risk CHIP and in 2 to the new diagnosis of a myeloproliferative neoplasm with indication for cytoreductive therapy.

CONCLUSIONS

Using genetic screening for myeloid disorders in patients with stroke of predominantly undetermined etiology, we found a 3-fold higher CHIP prevalence than in the general population. We identified high-risk CHIP and previously covert myeloproliferative neoplasms as potential stroke etiologies in 4.4% and 1% of patients, respectively. Our findings demonstrate the diagnostic and therapeutic yield of genetic screening in young patients with stroke. Future studies should investigate the role of CHIP for stroke recurrence and optimal secondary prevention.

Genome-wide association study in patients with posterior urethral valves

Congenital lower urinary tract obstructions (LUTO) are most often caused by posterior urethral valves (PUV), a male limited anatomical obstruction of the urethra affecting 1 in 4,000 male live births. Little is known about the genetic background of PUV. Here, we report the largest genome-wide association study (GWAS) for PUV in 4 cohorts of patients and controls. The final meta-analysis included 756 patients and 4,823 ethnicity matched controls and comprised 5,754,208 variants that were genotyped or imputed and passed quality control in all 4 cohorts. No genome-wide significant locus was identified, but 33 variants showed suggestive significance (P < 1 × 10^-5). When considering only loci with multiple variants residing within < 10 kB of each other showing suggestive significance and with the same effect direction in all 4 cohorts, 3 loci comprising a total of 9 variants remained. These loci resided on chromosomes 13, 16, and 20. The present GWAS and meta-analysis is the largest genetic study on PUV performed to date. The fact that no genome-wide significant locus was identified, can be explained by lack of power or may indicate that common variants do not play a major role in the etiology of PUV. Nevertheless, future studies are warranted to replicate and validate the 3 loci that yielded suggestive associations.

Recommendations for a Dutch Sustainable Biobanking Environment

Biobanks and their collections are considered essential for contemporary biomedical research and a critical resource toward personalized medicine. However, they need to operate in a sustainable manner to prevent research waste and maximize impact. Sustainability is the capacity of a biobank to remain operative, effective, and competitive over its expected lifetime. This remains a challenge given a biobank's position at the interplay of ethical, societal, scientific, and commercial values and the difficulties in finding continuous funding. In the end, biobanks are responsible for their own sustainability. Still, biobanks also depend on their surrounding environment, which contains overarching legislative, policy, financial, and other factors that can either impede or promote sustainability. The Biobanking and Biomolecular Research Infrastructure for The Netherlands (BBMRI.nl) has worked on improving the national environment for sustainable biobanking. In this article, we present the final outcomes of this BBMRI.nl project. First, we summarize the current overarching challenges of the Dutch biobanking landscape. These challenges were gathered during workshops and focus groups with Dutch biobanks and their users, for which the full results are described in separate reports. The main overarching challenges relate to sample and data quality, funding, use and reuse, findability and accessibility, and the general image of biobanks. Second, we propose a package of recommendations—across nine themes—toward creating overarching conditions that stimulate and enable sustainable biobanking. These recommendations serve as a guideline for the Dutch biobanking community and their stakeholders to jointly work toward practical implementation and a better biobanking environment. There are undoubtedly parallels between the Dutch situation and the challenges found in other countries. We hope that sharing our project's approach, outcomes, and recommendations will support other countries in their efforts toward sustainable biobanking.

Sex-Specific Regulation of Inflammation and Metabolic Syndrome in Obesity

Supplemental Digital Content is available in the text.

Objective:

Metabolic dysregulation and inflammation are important consequences of obesity and impact susceptibility to cardiovascular disease. Anti-inflammatory therapy in cardiovascular disease is being developed under the assumption that inflammatory pathways are identical in women and men, but it is not known if this is indeed the case. In this study, we assessed the sex-specific relation between inflammation and metabolic dysregulation in obesity.

Approach and Results:

Three hundred two individuals were included, half with a BMI 27 to 30 kg/m² and half with a BMI>30 kg/m², 45% were women. The presence of metabolic syndrome was assessed according to the National Cholesterol Education Program-ATPIII criteria, and inflammation was studied using circulating markers of inflammation, cell counts, and ex vivo cytokine production capacity of isolated immune cells. Additionally, lipidomic and metabolomic data were gathered, and subcutaneous fat biopsies were histologically assessed.

Metabolic syndrome is associated with an increased inflammatory profile that profoundly differs between women and men: women with metabolic syndrome show a lower concentration of the anti-inflammatory adiponectin, whereas men show increased levels of several pro-inflammatory markers such as IL (interleukin)-6 and leptin. Adipose tissue inflammation showed similar sex-specific associations with these markers. Peripheral blood mononuclear cells isolated from men, but not women, with metabolic syndrome display enhanced cytokine production capacity.

Conclusions:

We identified sex-specific pathways that influence inflammation in obesity. Excessive production of proinflammatory cytokines was observed in men with metabolic syndrome. In contrast, women typically showed reduced levels of the anti-inflammatory adipokine adiponectin. These different mechanisms of inflammatory dysregulation between women and men with obesity argue for sex-specific therapeutic strategies.

Measurement and genetic architecture of lifetime depression in the Netherlands as assessed by LIDAS (Lifetime Depression Assessment Self-report)

BACKGROUND

Major depressive disorder (MDD) is a common mood disorder, with a heritability of around 34%. Molecular genetic studies made significant progress and identified genetic markers associated with the risk of MDD; however, progress is slowed down by substantial heterogeneity as MDD is assessed differently across international cohorts. Here, we used a standardized online approach to measure MDD in multiple cohorts in the Netherlands and evaluated whether this approach can be used in epidemiological and genetic association studies of depression.

METHODS

Within the Biobank Netherlands Internet Collaboration (BIONIC) project, we collected MDD data in eight cohorts involving 31 936 participants, using the online Lifetime Depression Assessment Self-report (LIDAS), and estimated the prevalence of current and lifetime MDD in 22 623 unrelated individuals. In a large Netherlands Twin Register (NTR) twin-family dataset (n ≈ 18 000), we estimated the heritability of MDD, and the prediction of MDD in a subset (n = 4782) through Polygenic Risk Score (PRS).

RESULTS

Estimates of current and lifetime MDD prevalence were 6.7% and 18.1%, respectively, in line with population estimates based on validated psychiatric interviews. In the NTR heritability estimates were 0.34/0.30 (s.e. = 0.02/0.02) for current/lifetime MDD, respectively, showing that the LIDAS gives similar heritability rates for MDD as reported in the literature. The PRS predicted risk of MDD (OR 1.23, 95% CI 1.15-1.32, R2 = 1.47%).

CONCLUSIONS

By assessing MDD status in the Netherlands using the LIDAS instrument, we were able to confirm previously reported MDD prevalence and heritability estimates, which suggests that this instrument can be used in epidemiological and genetic association studies of depression.

Identification of ADHD risk genes in extended pedigrees by combining linkage analysis and whole-exome sequencing

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (N_total = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.

Non-Syndromic Cleft Lip with or without Cleft Palate: Genome-Wide Association Study in Europeans Identifies a Suggestive Risk Locus at 16p12.1 and Supports SH3PXD2A as a Clefting Susceptibility Gene

Non-syndromic cleft lip with or without cleft palate (nsCL/P) ranks among the most common human congenital malformations, and has a multifactorial background in which both exogenous and genetic risk factors act in concert. The present report describes a genome-wide association study (GWAS) involving a total of 285 nsCL/P patients and 1212 controls from the Netherlands and Belgium. Twenty of the 40 previously reported nsC/LP susceptibility loci were replicated, which underlined the validity of this sample. SNV-based analysis of the data identified an as yet unreported suggestive locus at chromosome 16p12.1 (p-value of the lead SNV: 4.17 × 10^-7). This association was replicated in two of three patient/control replication series (Central European and Yemeni). Gene analysis of the GWAS data prioritized SH3PXD2A at chromosome 10q24.33 as a candidate gene for nsCL/P. To date, support for this gene as a cleft gene has been restricted to data from zebrafish and a knockout mouse model. The present GWAS was the first to implicate SH3PXD2A in non-syndromic cleft formation in humans. In summary, although performed in a relatively small sample, the present GWAS generated novel insights into nsCL/P etiology.

Role of the Complement System in Chronic Central Serous Chorioretinopathy: A Genome-Wide Association Study

Importance

To date, several targeted genetic studies on chronic central serous chorioretinopathy (cCSC) have been performed; however, unbiased genome-wide studies into the genetics of cCSC have not been reported. To discover new genetic loci associated with cCSC and to better understand the causative mechanism of this disease, we performed a genome-wide association study (GWAS) on patients with cCSC.

Objective

To discover new genetic loci and pathways associated with cCSC and to predict the association of genetic variants with gene expression in patients with cCSC.

Design, Setting, and Participants

This case-control GWAS was completed in the general community, 3 referral university medical centers, and outpatient care on Europeans individuals with cCSC and population-based control participants. Genotype data was collected from May 2013 to August 2017, and data analysis occurred from August 2017 to November 2017.

Main Outcomes and Measures

Associations of single-nucleotide polymorphisms, haplotypes, genetic pathways, and predicted gene expression with cCSC.

Results

A total of 521 patients with cCSC (median age, 51 years; interquartile range [IQR], 44-59 years; 420 [80.6%] male) and 3577 European population-based control participants (median age, 52 years; IQR, 37-71 years; 1630 [45.6%] male) were included. One locus on chromosome 1 at the complement factor H (CFH) gene reached genome-wide significance and was associated with an increased risk of cCSC (rs1329428; odds ratio [OR], 1.57 [95% CI, 1.38-1.80]; P = 3.12 × 10-11). The CFH haplotypes H1 and H3 were protective for cCSC (H1: OR, 0.64 [95% CI, 0.53-0.77]; P = 2.18 × 10-6; H3: OR, 0.54 [95% CI, 0.42-0.70]; P = 2.49 × 10-6), whereas haplotypes H2, H4, H5, and the aggregate of rare CFH haplotypes conferred increased risk (H2: OR, 1.57 [95% CI, 1.30-1.89]; P = 2.18 × 10-6; H4: OR, 1.43 [95% CI, 1.13-1.80]; P = 2.49 × 10-3; H5: OR, 1.80 [95% CI, 1.36-2.39]; P = 4.61 × 10-5; rare haplotypes: OR, 1.99 [95% CI, 1.43-2.77]; P = 4.59 × 10-5). Pathway analyses showed involvement of the complement cascade and alternative open reading frame (ARF) pathway in cCSC. Using PrediXcan, we identified changes in predicted expression of complement genes CFH, complement factor H related 1 (CFHR1), complement factor related 4 (CFHR4), and membrane cofactor protein (MCP/CD46). Additionally, the potassium sodium-activated channel subfamily T member 2 (KCNT2) and tumor necrosis factor receptor superfamily member 10a (TNFRSF10A) genes were differentially expressed in patients with cCSC.

Conclusions and Relevance

In this GWAS on cCSC, we identified a locus on chromosome 1 at the CFH gene that was significantly associated with cCSC, and we report protective and risk-conferring haplotypes in this gene. Pathway analyses were enriched for complement genes, and gene expression analysis suggests a role for CFH, CFHR1, CFHR4, CD46, KCNT2, and TNFRSF10A in the disease. Taken together, these results underscore the potential importance of the complement pathway in the causative mechanisms of cCSC.

Cigarette Smoking and the Risk of Cutaneous Melanoma: A Case-Control Study

OBJECTIVE

Previous studies regarding cigarette smoking causing a lower risk of melanoma are inconclusive. Here, we re-examined melanoma risk in relation to cigarette smoking in a large, case-control study.

METHODS

In total 1,157 patients with melanoma diagnosed between 2003 and 2011 in the Netherlands and 5,595 controls from the Nijmegen Biomedical Study were included. Information concerning smoking habits and known risk factors for melanoma were obtained through self-administered questionnaires. Logistic regression analyses stratified by gender were performed to study the risk of cigarette smoking on melanoma risk, adjusted for age, marital status, highest level of education, skin type, sun vacation, use of solarium, time spent outdoors, and sun protective measures.

RESULTS

Among men, current and former smokers did not have a higher risk of melanoma compared to never smokers: adjusted odds ratio (OR) = 0.56 (95% confidence interval [CI]: 0.40-0.79) and adjusted OR = 0.50 (95% CI: 0.39-0.64), respectively. With an increasing number of years smoked the risk of melanoma decreased: <20 years: OR = 0.61 (95% CI: 0.46-0.80); 21-40 years: OR = 0.50 (95% CI: 0.37-0.68); >40 years: OR = 0.26 (95% CI: 0.15-0.44). No clear trend was found for the number of cigarettes smoked. Results for females were less clear and not statistically significant (current smoker: adjusted OR = 0.96, 95% CI: 0.74-1.26, former smoker: adjusted OR = 0.89, 95% CI: 0.73-1.08).

CONCLUSION

This study shows a strong inverse association between cigarette smoking and melanoma risk in men. Fundamental laboratory research is necessary to investigate the biological relation between smoking cigarettes and melanoma.

Exome chip association study excluded the involvement of rare coding variants with large effect sizes in the etiology of anorectal malformations

INTRODUCTION

Anorectal malformations (ARM) are rare congenital malformations, resulting from disturbed hindgut development. A genetic etiology has been suggested, but evidence for the involvement of specific genes is scarce. We evaluated the contribution of rare and low-frequency coding variants in ARM etiology, assuming a multifactorial model.

METHODS

We analyzed 568 Caucasian ARM patients and 1,860 population-based controls using the Illumina HumanExome Beadchip array, which contains >240,000 rare and low-frequency coding variants. GenomeStudio clustering and calling was followed by re-calling of 'no-calls' using zCall for patients and controls simultaneously. Single variant and gene-based analyses were performed to identify statistically significant associations, applying Bonferroni correction. Following an extra quality control step, candidate variants were selected for validation using Sanger sequencing.

RESULTS

When we applied a MAF of ≥1.0%, no variants or genes showed statistically significant associations with ARM. Using a MAF cut-off at 0.4%, 13 variants initially reached statistical significance, but had to be discarded upon further inspection: ten variants represented calling errors of the software, while the minor alleles of the remaining three variants were not confirmed by Sanger sequencing.

CONCLUSION

Our results show that rare and low-frequency coding variants with large effect sizes, present on the exome chip do not contribute to ARM etiology.

Donor InSight: characteristics and representativeness of a Dutch cohort study on blood and plasma donors

Background and Objectives

More insight into donor health and behaviour may contribute to more efficient and focused strategies regarding donor care and management. Donor InSight (DIS) is a Dutch cohort study of blood and plasma donors. We aimed to outline the objectives and methods of DIS, describe the cohort, and compare it to the active Dutch donor population.

Materials and Methods

In 2007‐2009 (DIS‐I, n = 31 338) and 2012‐2013 (DIS‐II, 34 826, of whom 22 132 also participated in DIS‐I) questionnaire data on demographics, donation, lifestyle, family composition, health and disease were collected. A second follow‐up (DIS‐III, n = 3046), including donors with differing haemoglobin trajectories, was completed in 2015‐2016. DIS‐III includes data on genetic determinants, iron and red cell indices. Representativeness of the DIS‐I sample for the entire Dutch donor population was assessed by comparing characteristics of both.

Results

Donor InSight was initially set up because of a need for more detailed information and evidence as a basis for decision‐making in blood banks. DIS‐I sample is comparable to the total Dutch donor population in terms of age, body mass index, haemoglobin level, blood pressure, blood type and donation behaviour.

Conclusion

Donor InSight is a cohort study representative of the Dutch donor population. It provides evidence to support evidence‐based decision making.

Association of BRCA2 K3326* With Small Cell Lung Cancer and Squamous Cell Cancer of the Skin

Background

Most pathogenic mutations in the BRCA2 gene carry a high risk of hereditary breast and ovarian cancer (HBOC). However, a stop-gain mutation, K3326* (rs11571833), confers risk of lung cancer and cancers of the upper-aero-digestive tract but only a modest risk of breast or ovarian cancer. The Icelandic population provides an opportunity for comprehensive characterization of the cancer risk profiles of K3326* and HBOC mutations because a single mutation, BRCA2 999del5, is responsible for almost all BRCA2-related HBOC in the population.

Methods

Genotype information on 43 641 cancer patients and 370 971 control subjects from Iceland, the Netherlands, and the United States was used to assess the cancer risk profiles of K3326* and BRCA2 999del5. BRCA2 expression was assessed using RNAseq data from blood (n = 2233), as well as 52 tissues reported in the GTEx database.

Results

The cancer risks associated with K3326* are fundamentally different from those associated with 999del5. We report for the first time an association between K3326* and small cell lung cancer (odds ratio [OR] = 2.06, 95% confidence interval [CI] = 1.35 to 3.16) and squamous cell carcinoma of the skin (OR = 1.69, 95% CI = 1.26 to 2.26). Individuals homozygous for K3326* reach old age and have children. Unlike BRCA2 999del5, the K3326* allele does not affect the level of BRCA2 transcripts, and the allele is expressed to the same extent as the wild-type allele.

Conclusions

K3326* associates primarily with cancers that have strong environmental genotoxic risk factors. Expression of the K3326* allele suggests that a variant protein may be made that retains the DNA repair capabilities important to hormone-responsive tissues but may be less efficient in responding to genotoxic stress.

Autism spectrum disorders and autistic traits share genetics and biology

Autism spectrum disorders (ASDs) and autistic traits in the general population may share genetic susceptibility factors. In this study, we investigated such potential overlap based on common genetic variants. We developed and validated a self-report questionnaire of autistic traits in adults. We then conducted genome-wide association studies (GWASs) of six trait scores derived from the questionnaire through exploratory factor analysis in 1981 adults from the general population. Using the results from the Psychiatric Genomics Consortium GWAS of ASDs, we observed genetic sharing between ASDs and the autistic traits 'childhood behavior', 'rigidity' and 'attention to detail'. Gene-set analysis subsequently identified 'rigidity' to be significantly associated with a network of ASD gene-encoded proteins that regulates neurite outgrowth. Gene-wide association with the well-established ASD gene MET reached significance. Taken together, our findings provide evidence for an overlapping genetic and biological etiology underlying ASDs and autistic population traits, which suggests that genetic studies in the general population may yield novel ASD genes.

Ultra-sensitive Sequencing Identifies High Prevalence of Clonal Hematopoiesis-Associated Mutations throughout Adult Life

Clonal hematopoiesis results from somatic mutations in hematopoietic stem cells, which give an advantage to mutant cells, driving their clonal expansion and potentially leading to leukemia. The acquisition of clonal hematopoiesis-driver mutations (CHDMs) occurs with normal aging and these mutations have been detected in more than 10% of individuals ≥65 years. We aimed to examine the prevalence and characteristics of CHDMs throughout adult life. We developed a targeted re-sequencing assay combining high-throughput with ultra-high sensitivity based on single-molecule molecular inversion probes (smMIPs). Using smMIPs, we screened more than 100 loci for CHDMs in more than 2,000 blood DNA samples from population controls between 20 and 69 years of age. Loci screened included 40 regions known to drive clonal hematopoiesis when mutated and 64 novel candidate loci. We identified 224 somatic mutations throughout our cohort, of which 216 were coding mutations in known driver genes (DNMT3A, JAK2, GNAS, TET2, and ASXL1), including 196 point mutations and 20 indels. Our assay's improved sensitivity allowed us to detect mutations with variant allele frequencies as low as 0.001. CHDMs were identified in more than 20% of individuals 60 to 69 years of age and in 3% of individuals 20 to 29 years of age, approximately double the previously reported prevalence despite screening a limited set of loci. Our findings support the occurrence of clonal hematopoiesis-associated mutations as a widespread mechanism linked with aging, suggesting that mosaicism as a result of clonal evolution of cells harboring somatic mutations is a universal mechanism occurring at all ages in healthy humans.

A genome-wide association study yields five novel thyroid cancer risk loci

The great majority of thyroid cancers are of the non-medullary type. Here we report findings from a genome-wide association study of non-medullary thyroid cancer, including in total 3,001 patients and 287,550 controls from five study groups of European descent. Our results yield five novel loci (all with Pcombined<3 × 10-8): 1q42.2 (rs12129938 in PCNXL2), 3q26.2 (rs6793295 a missense mutation in LRCC34 near TERC), 5q22.1 (rs73227498 between NREP and EPB41L4A), 10q24.33 (rs7902587 near OBFC1), and two independently associated variants at 15q22.33 (rs2289261 and rs56062135; both in SMAD3). We also confirm recently published association results from a Chinese study of a variant on 5p15.33 (rs2736100 near the TERT gene) and present a stronger association result for a moderately correlated variant (rs10069690; OR=1.20, P=3.2 × 10-7) based on our study of individuals of European ancestry. In combination, these results raise several opportunities for future studies of the pathogenesis of thyroid cancer.