Next-generation sequencing in clinical diagnostics is providing valuable genomic variant data, which can be used to support healthcare decisions. In silico tools to predict pathogenicity are crucial to assess such variants and we have evaluated a new tool, Combined Annotation Dependent Depletion (CADD), and its classification of gene variants in Lynch syndrome by using a set of 2,210 DNA mismatch repair gene variants. These had already been classified by experts from InSiGHT's Variant Interpretation Committee. Overall, we found CADD scores do predict pathogenicity (Spearman's ρ = 0.595, P < 0.001). However, we discovered 31 major discrepancies between the InSiGHT classification and the CADD scores; these were explained in favor of the expert classification using population allele frequencies, cosegregation analyses, disease association studies, or a second-tier test. Of 751 variants that could not be clinically classified by InSiGHT, CADD indicated that 47 variants were worth further study to confirm their putative pathogenicity. We demonstrate CADD is valuable in prioritizing variants in clinically relevant genes for further assessment by expert classification teams.

Mutations create variation in the population, fuel evolution and cause genetic diseases. Current knowledge about de novo mutations is incomplete and mostly indirect. Here we analyze 11,020 de novo mutations from the whole genomes of 250 families. We show that de novo mutations in the offspring of older fathers are not only more numerous but also occur more frequently in early-replicating, genic regions. Functional regions exhibit higher mutation rates due to CpG dinucleotides and show signatures of transcription-coupled repair, whereas mutation clusters with a unique signature point to a new mutational mechanism. Mutation and recombination rates independently associate with nucleotide diversity, and regional variation in human-chimpanzee divergence is only partly explained by heterogeneity in mutation rate. Finally, we provide a genome-wide mutation rate map for medical and population genetics applications. Our results provide new insights and refine long-standing hypotheses about human mutagenesis.

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of The Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value <6.61 × 10(-4)), including a rare missense variant in ABCA6 (rs77542162, p.Cys1359Arg, frequency 0.034), which is predicted to be deleterious. The frequency of this ABCA6 variant is 3.65-fold increased in the Dutch and its effect (βLDL-C=0.135, βTC=0.140) is estimated to be very similar to those observed for single variants in well-known lipid genes, such as LDLR.

Many cancer-associated somatic copy number alterations (SCNAs) are known. Currently, one of the challenges is to identify the molecular downstream effects of these variants. Although several SCNAs are known to change gene expression levels, it is not clear whether each individual SCNA affects gene expression. We reanalyzed 77,840 expression profiles and observed a limited set of 'transcriptional components' that describe well-known biology, explain the vast majority of variation in gene expression and enable us to predict the biological function of genes. On correcting expression profiles for these components, we observed that the residual expression levels (in 'functional genomic mRNA' profiling) correlated strongly with copy number. DNA copy number correlated positively with expression levels for 99% of all abundantly expressed human genes, indicating global gene dosage sensitivity. By applying this method to 16,172 patient-derived tumor samples, we replicated many loci with aberrant copy numbers and identified recurrently disrupted genes in genomically unstable cancers.

Studies in European populations have contributed to a better understanding of the genetics of complex diseases, for example, in coeliac disease (CeD), studies of over 23 000 European samples have reported association to the HLA locus and another 39 loci. However, these associations have not been evaluated in detail in other ethnicities. We sought to better understand how disease-associated loci that have been mapped in Europeans translate to a disease risk for a population with a different ethnic background. We therefore performed a validation of European risk loci for CeD in 497 cases and 736 controls of north Indian origin. Using a dense-genotyping platform (Immunochip), we confirmed the strong association to the HLA region (rs2854275, P=8.2 × 10(-49)). Three loci showed suggestive association (rs4948256, P=9.3 × 10(-7), rs4758538, P=8.6 × 10(-5) and rs17080877, P=2.7 × 10(-5)). We directly replicated five previously reported European variants (P<0.05; mapping to loci harbouring FASLG/TNFSF18, SCHIP1/IL12A, PFKFB3/PRKCQ, ZMIZ1 and ICOSLG). Using a transferability test, we further confirmed association at PFKFB3/PRKCQ (rs2387397, P=2.8 × 10(-4)) and PTPRK/THEMIS (rs55743914, P=3.4 × 10(-4)). The north Indian population has a higher degree of consanguinity than Europeans and we therefore explored the role of recessively acting variants, which replicated the HLA locus (rs9271850, P=3.7 × 10(-23)) and suggested a role of additional four loci. To our knowledge, this is the first replication study of CeD variants in a non-European population.

To perform a genome-wide association study (GWAS) using the Immunochip array in 3,420 cases of ischemic stroke and 6,821 controls, followed by a meta-analysis with data from more than 14,000 additional ischemic stroke cases.

The Copper Metabolism MURR1 domain protein 1 (COMMD1) is a protein involved in multiple cellular pathways, including copper homeostasis, NF-κB and hypoxia signalling. Acting as a scaffold protein, COMMD1 mediates the levels, stability and proteolysis of its substrates (e.g. the copper-transporters ATP7B and ATP7A, RELA and HIF-1α). Recently, we established an interaction between the Cu/Zn superoxide dismutase 1 (SOD1) and COMMD1, resulting in a decreased maturation and activation of SOD1. Mutations in SOD1, associated with the progressive neurodegenerative disorder Amyotrophic Lateral Sclerosis (ALS), cause misfolding and aggregation of the mutant SOD1 (mSOD1) protein. Here, we identify COMMD1 as a novel regulator of misfolded protein aggregation as it enhances the formation of mSOD1 aggregates upon binding. Interestingly, COMMD1 co-localizes to the sites of mSOD1 inclusions and forms high molecular weight complexes in the presence of mSOD1. The effect of COMMD1 on protein aggregation is client-specific as, in contrast to mSOD1, COMMD1 decreases the abundance of mutant Parkin inclusions, associated with Parkinson's disease. Aggregation of a polyglutamine-expanded Huntingtin, causative of Huntington's disease, appears unaltered by COMMD1. Altogether, this study offers new research directions to expand our current knowledge on the mechanisms underlying aggregation disease pathologies.

A gluten-free diet (GFD) is the most commonly adopted special diet worldwide. It is an effective treatment for coeliac disease and is also often followed by individuals to alleviate gastrointestinal complaints. It is known there is an important link between diet and the gut microbiome, but it is largely unknown how a switch to a GFD affects the human gut microbiome.

Effective immunity requires a complex network of cellular and humoral components that interact with each other and are influenced by different environmental and host factors. We used a systems biology approach to comprehensively assess the impact of environmental and genetic factors on immune cell populations in peripheral blood, including associations with immunoglobulin concentrations, from ∼500 healthy volunteers from the Human Functional Genomics Project. Genetic heritability estimation showed that variations in T cell numbers are more strongly driven by genetic factors, while B cell counts are more environmentally influenced. Quantitative trait loci (QTL) mapping identified eight independent genomic loci associated with leukocyte count variation, including four associations with T and B cell subtypes. The QTLs identified were enriched among genome-wide association study (GWAS) SNPs reported to increase susceptibility to immune-mediated diseases. Our systems approach provides insights into cellular and humoral immune trait variability in humans.

Approximately 1.5 billion people worldwide are overweight or affected by obesity, and are at risk of developing type 2 diabetes, cardiovascular disease and related metabolic and inflammatory disturbances. Although the mechanisms linking adiposity to associated clinical conditions are poorly understood, recent studies suggest that adiposity may influence DNA methylation, a key regulator of gene expression and molecular phenotype. Here we use epigenome-wide association to show that body mass index (BMI; a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci with P < 1 × 10-7, range P = 9.2 × 10-8 to 6.0 × 10-46; n = 10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find that methylation loci are enriched for functional genomic features in multiple tissues (P < 0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P < 9.0 × 10-6, range P = 5.5 × 10-6 to 6.1 × 10-35, n = 1,785 samples). The methylation loci identify genes involved in lipid and lipoprotein metabolism, substrate transport and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future development of type 2 diabetes (relative risk per 1 standard deviation increase in methylation risk score: 2.3 (2.07-2.56); P = 1.1 × 10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type 2 diabetes and other adverse clinical consequences of obesity.

Using the Immunochip custom SNP array, which was designed for dense genotyping of 186 loci identified through genome-wide association studies (GWAS), we analyzed 11,475 individuals with rheumatoid arthritis (cases) of European ancestry and 15,870 controls for 129,464 markers. We combined these data in a meta-analysis with GWAS data from additional independent cases (n = 2,363) and controls (n = 17,872). We identified 14 new susceptibility loci, 9 of which were associated with rheumatoid arthritis overall and five of which were specifically associated with disease that was positive for anticitrullinated peptide antibodies, bringing the number of confirmed rheumatoid arthritis risk loci in individuals of European ancestry to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at 6 loci and identified association to low-frequency variants at 4 loci. Bioinformatic analyses generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.

In addition to HLA genetic incompatibility, non-HLA difference between donor and recipients of transplantation leading to allograft rejection are now becoming evident. We aimed to create a unique genome-wide platform to facilitate genomic research studies in transplant-related studies. We designed a genome-wide genotyping tool based on the most recent human genomic reference datasets, and included customization for known and potentially relevant metabolic and pharmacological loci relevant to transplantation.

Candida albicans is the most common human fungal pathogen causing mucosal and systemic infections. However, human antifungal immunity remains poorly defined. Here by integrating transcriptional analysis and functional genomics, we identified Candida-specific host defence mechanisms in humans. Candida induced significant expression of genes from the type I interferon pathway in human peripheral blood mononuclear cells. This unexpectedly prominent role of type I interferon pathway in anti-Candida host defence was supported by additional evidence. Polymorphisms in type I interferon genes modulated Candida-induced cytokine production and were correlated with susceptibility to systemic candidiasis. In in vitro experiments, type I interferons skewed Candida-induced inflammation from a Th17 response towards a Th1 response. Patients with chronic mucocutaneous candidiasis displayed defective expression of genes in the type I interferon pathway. These findings indicate that the type I interferon pathway is a main signature of Candida-induced inflammation and has a crucial role in anti-Candida host defence in humans.

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV(1)), and its ratio to forced vital capacity (FEV(1)/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV(1) and FEV(1)/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P(JMA = )5.00×10(-11)), HLA-DQB1 and HLA-DQA2 (smallest P(JMA = )4.35×10(-9)), and KCNJ2 and SOX9 (smallest P(JMA = )1.28×10(-8)) were associated with FEV(1)/FVC or FEV(1) in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

Genome-wide association studies in Japanese populations recently identified common variants in the KCNQ1 gene to be associated with type 2 diabetes. We examined the association of these variants within KCNQ1 with type 2 diabetes in a Dutch population, investigated their effects on insulin secretion and metabolic traits and on the risk of developing complications in type 2 diabetes patients.

There is no consensus on the preferred method for defining the non-inferiority margin in non-inferiority trials, and previous studies showed that the rationale for its choice is often not reported. This study investigated how the non-inferiority margin is defined in the published literature, and whether its reporting has changed over time.

Tissue-resident lymphocytes play a key role in immune surveillance, but it remains unclear how these inherently stable cell populations respond to chronic inflammation. In the setting of celiac disease (CeD), where exposure to dietary antigen can be controlled, gluten-induced inflammation triggered a profound depletion of naturally occurring Vγ4+/Vδ1+ intraepithelial lymphocytes (IELs) with innate cytolytic properties and specificity for the butyrophilin-like (BTNL) molecules BTNL3/BTNL8. Creation of a new niche with reduced expression of BTNL8 and loss of Vγ4+/Vδ1+ IELs was accompanied by the expansion of gluten-sensitive, interferon-γ-producing Vδ1+ IELs bearing T cell receptors (TCRs) with a shared non-germline-encoded motif that failed to recognize BTNL3/BTNL8. Exclusion of dietary gluten restored BTNL8 expression but was insufficient to reconstitute the physiological Vγ4+/Vδ1+ subset among TCRγδ+ IELs. Collectively, these data show that chronic inflammation permanently reconfigures the tissue-resident TCRγδ+ IEL compartment in CeD. VIDEO ABSTRACT.

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10-10), among which transcription factors were overrepresented (Fisher's P = 3.3 × 10-7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser.

Due to limited sepsis patient cohort size and extreme heterogeneity, only one significant locus and suggestive associations at several independent loci were implicated by three genome-wide association studies. However, genes from such suggestive loci may also provide crucial information to unravel genetic mechanisms that determine sepsis heterogeneity. Therefore, in this study, we made use of integrative approaches to prioritize genes and pathways affected by sepsis associated genetic variants. By integrating expression quantitative trait loci (eQTL) results from the largest whole-blood eQTL database, cytokine QTLs from pathogen-stimulated peripheral blood mononuclear cells (PBMCs), publicly available blood transcriptome data from pneumoniae-derived sepsis patients, and transcriptome data from pathogen-stimulated PBMCs, we identified 55 potential genes affected by 39 independent loci. By performing pathway enrichment analysis at these loci we found enrichment of genes for adherences-junction pathway. Finally, we investigated the functional role of the only one GWAS significant SNP rs4957796 on sepsis survival in altering transcription factor binding affinity in monocytes and endothelial cells. We also found that transient deficiency of FER and MAN2A1 affect endothelial response to stimulation, indicating that both FER and MAN2A1 could be the causal genes at this locus. Taken together, our study suggests that in addition to immune pathways, genetic variants may also affect non-immune related pathways.

The risk of peripheral artery disease (PAD) in patients with diabetes mellitus (DM) and coronary artery disease (CAD) is an important and inadequately addressed issue. Our aim is to examine the impact of DM on risk of PAD in patients with different degrees of CAD characterized by coronary angiography (CAG).