Cas9-directed long-read sequencing to resolve optical genome mapping findings in leukemia diagnostics.

Leukemias are genetically heterogeneous and diagnostics therefore includes various standard-of-care (SOC) techniques, including karyotyping, SNP-array and FISH. Optical genome mapping (OGM) may replace these as it detects different types of structural aberrations simultaneously and additionally detects much smaller aberrations (500 bp vs 5-10 Mb with karyotyping). However, its resolution may still be too low to define clinical relevance of aberrations when they are located between two OGM labels or when labels are not distinct enough. Here, we test the potential of Cas9-directed long-read sequencing (LRS) as an additional technique to resolve such potentially relevant new findings. From an internal Bionano implementation study we selected ten OGM calls that could not be validated with SOC methods. Per variant we designed crRNAs for Cas9 enrichment, prepared libraries and sequenced them on a MinION/GridION device. We could confirm all aberrations and, importantly, the actual breakpoints of the OGM calls were located between 0.2 and 5.5 kb of the OGM-estimated breakpoints, confirming the high reliability of OGM. Furthermore, we show examples of redefinition of aberrations between labels that enable judgment of clinical relevance. Our results suggest that Cas9-directed LRS can be a relevant and flexible secondary technique in diagnostic workflows including OGM.

Targeted RNA-Sequencing Enables Detection of Relevant Translocations and Single Nucleotide Variants and Provides a Method for Classification of Hematological Malignancies-RANKING.

BACKGROUND: Patients with hematological malignancies (HMs) carry a wide range of chromosomal and molecular abnormalities that impact their prognosis and treatment. Since no current technique can detect all relevant abnormalities, technique(s) are chosen depending on the reason for referral, and abnormalities can be missed. We tested targeted transcriptome sequencing as a single platform to detect all relevant abnormalities and compared it to current techniques. MATERIAL AND METHODS: We performed RNA-sequencing of 1385 genes (TruSight RNA Pan-Cancer, Illumina) in bone marrow from 136 patients with a primary diagnosis of HM. We then applied machine learning to expression profile data to perform leukemia classification, a method we named RANKING. Gene fusions for all the genes in the panel were detected, and overexpression of the genes EVI1, CCND1, and BCL2 was quantified. Single nucleotide variants/indels were analyzed in acute myeloid leukemia (AML), myelodysplastic syndrome and patients with acute lymphoblastic leukemia (ALL) using a virtual myeloid (54 genes) or lymphoid panel (72 genes). RESULTS: RANKING correctly predicted the leukemia classification of all AML and ALL samples and improved classification in 3 patients. Compared to current methods, only one variant was missed, c.2447A>T in KIT (RT-PCR at 10-4), and BCL2 overexpression was not seen due to a t(14; 18)(q32; q21) in 2% of the cells. Our RNA-sequencing method also identified 6 additional fusion genes and overexpression of CCND1 due to a t(11; 14)(q13; q32) in 2 samples. CONCLUSIONS: Our combination of targeted RNA-sequencing and data analysis workflow can improve the detection of relevant variants, and expression patterns can assist in establishing HM classification.

Detection of Fusion Genes to Determine Minimal Residual Disease in Leukemia Using Next-Generation Sequencing.

BACKGROUND: Measuring minimal residual disease (MRD), the persistence of leukemic cells after treatment, is important for monitoring leukemia recurrence. The current methods for monitoring MRD are flow cytometry, to assess aberrant immune phenotypes, and digital droplet PCR (ddPCR), to target genetic aberrations such as single-nucleotide variants and gene fusions. We present the performance of an RNA-based next-generation sequencing (NGS) method for MRD gene fusion detection compared with ddPCR. This method may have advantages, including the capacity to analyze different genetic aberrations and patients in 1 experiment. In particular, detection at the RNA level may be highly sensitive if the genetic aberration is highly expressed. METHODS: We designed a probe-based NGS panel targeting the breakpoints of 11 fusion genes previously identified in clinical patients and 2 fusion genes present in cell lines. Blocking probes were added to prevent nonspecific enrichment. Each patient RNA sample was diluted in background RNA, depleted for rRNA and globin mRNA, converted to cDNA, and prepared for sequencing. Unique sequence reads, identified by unique molecular identifiers, were aligned directly to reference transcripts. The same patient and cell-line samples were also analyzed with ddPCR for direct comparison. RESULTS: Our NGS method reached a maximum sensitivity of 1 aberrant cell in 10 000 cells and was mostly within a factor of 10 compared with ddPCR. CONCLUSIONS: Our detection limit was below the threshold of 1:1000 recommended by European Leukemia Net. Further optimizations are easy to implement and are expected to boost the sensitivity of our method to diagnostically obtained ddPCR thresholds.

A homozygous variant in growth and differentiation factor 2 (GDF2) may cause lymphatic dysplasia with hydrothorax and nonimmune hydrops fetalis.

The etiology of nonimmune hydrops fetalis is extensive and includes genetic disorders. We describe a term-born female neonate with late onset extensive nonimmune hydrops, that is, polyhydramnios, edema, and congenital bilateral chylothorax. This newborn was successfully treated with repetitive thoracocentesis, total parenteral feeding, octreotide intravenously and finally surgical pleurodesis and corticosteroids. A genetic cause seemed plausible as the maternal history revealed a fatal nonimmune hydrops fetalis. A homozygous truncating variant in GDF2 (c.451C>T, p.(Arg151*)) was detected with exome sequencing. Genetic analysis of tissue obtained from the deceased fetal sibling revealed the same homozygous variant. The parents and two healthy siblings were heterozygous for the GDF2 variant. Skin and lung biopsies in the index patient, as well as the revised lung biopsy of the deceased fetal sibling, showed lymphatic dysplasia and lymphangiectasia. To the best of our knowledge, this is the first report of an association between a homozygous variant in GDF2 with lymphatic dysplasia, hydrothorax and nonimmune hydrops fetalis.

A prospective study on rapid exome sequencing as a diagnostic test for multiple congenital anomalies on fetal ultrasound.

OBJECTIVE: Conventional genetic tests (quantitative fluorescent-PCR [QF-PCR] and single nucleotide polymorphism-array) only diagnose ~40% of fetuses showing ultrasound abnormalities. Rapid exome sequencing (rES) may improve this diagnostic yield, but includes challenges such as uncertainties in fetal phenotyping, variant interpretation, incidental unsolicited findings, and rapid turnaround times. In this study, we implemented rES in prenatal care to increase diagnostic yield. METHODS: We prospectively studied 55 fetuses. Inclusion criteria were: (a) two or more independent major fetal anomalies, (b) hydrops fetalis or bilateral renal cysts alone, or (c) one major fetal anomaly and a first-degree relative with the same anomaly. In addition to conventional genetic tests, we performed trio rES analysis using a custom virtual gene panel of ~3850 Online Mendelian Inheritance in Man (OMIM) genes. RESULTS: We established a genetic rES-based diagnosis in 8 out of 23 fetuses (35%) without QF-PCR or array abnormalities. Diagnoses included MIRAGE (SAMD9), Zellweger (PEX1), Walker-Warburg (POMGNT1), Noonan (PTNP11), Kabuki (KMT2D), and CHARGE (CHD7) syndrome and two cases of Osteogenesis Imperfecta type 2 (COL1A1). In six cases, rES diagnosis aided perinatal management. The median turnaround time was 14 (range 8-20) days. CONCLUSION: Implementing rES as a routine test in the prenatal setting is challenging but technically feasible, with a promising diagnostic yield and significant clinical relevance.

Occupational exposure to gases/fumes and mineral dust affect DNA methylation levels of genes regulating expression.

Many workers are daily exposed to occupational agents like gases/fumes, mineral dust or biological dust, which could induce adverse health effects. Epigenetic mechanisms, such as DNA methylation, have been suggested to play a role. We therefore aimed to identify differentially methylated regions (DMRs) upon occupational exposures in never-smokers and investigated if these DMRs associated with gene expression levels. To determine the effects of occupational exposures independent of smoking, 903 never-smokers of the LifeLines cohort study were included. We performed three genome-wide methylation analyses (Illumina 450 K), one per occupational exposure being gases/fumes, mineral dust and biological dust, using robust linear regression adjusted for appropriate confounders. DMRs were identified using comb-p in Python. Results were validated in the Rotterdam Study (233 never-smokers) and methylation-expression associations were assessed using Biobank-based Integrative Omics Study data (n = 2802). Of the total 21 significant DMRs, 14 DMRs were associated with gases/fumes and 7 with mineral dust. Three of these DMRs were associated with both exposures (RPLP1 and LINC02169 (2×)) and 11 DMRs were located within transcript start sites of gene expression regulating genes. We replicated two DMRs with gases/fumes (VTRNA2-1 and GNAS) and one with mineral dust (CCDC144NL). In addition, nine gases/fumes DMRs and six mineral dust DMRs significantly associated with gene expression levels. Our data suggest that occupational exposures may induce differential methylation of gene expression regulating genes and thereby may induce adverse health effects. Given the millions of workers that are exposed daily to occupational exposures, further studies on this epigenetic mechanism and health outcomes are warranted.

Limited overlap in significant hits between genome-wide association studies on two airflow obstruction definitions in the same population.

BACKGROUND: Airflow obstruction is a hallmark of chronic obstructive pulmonary disease (COPD), and is defined as either the ratio between forced expiratory volume in one second and forced vital capacity (FEV1/FVC) < 70% or < lower limit of normal (LLN). This study aimed to assess the overlap between genome-wide association studies (GWAS) on airflow obstruction using these two definitions in the same population stratified by smoking. METHODS: GWASes were performed in the LifeLines Cohort Study for both airflow obstruction definitions in never-smokers (NS = 5071) and ever-smokers (ES = 4855). The FEV1/FVC < 70% models were adjusted for sex, age, and height; FEV1/FVC < LLN models were not adjusted. Ever-smokers models were additionally adjusted for pack-years and current-smoking. The overlap in significantly associated SNPs between the two definitions and never/ever-smokers was assessed using several p-value thresholds. To quantify the agreement, the Pearson correlation coefficient was calculated between the p-values and ORs. Replication was performed in the Vlagtwedde-Vlaardingen study (NS = 432, ES = 823). The overlapping SNPs with p < 10- 4 were validated in the Vlagtwedde-Vlaardingen and Rotterdam Study cohorts (NS = 1966, ES = 3134) and analysed for expression quantitative trait loci (eQTL) in lung tissue (n = 1087). RESULTS: In the LifeLines cohort, 96% and 93% of the never- and ever-smokers were classified concordantly based on the two definitions. 26 and 29% of the investigated SNPs were overlapping at p < 0.05 in never- and ever-smokers, respectively. At p < 10- 4 the overlap was 4% and 6% respectively, which could be change findings as shown by simulation studies. The effect estimates of the SNPs of the two definitions correlated strongly, but the p-values showed more variation and correlated only moderately. Similar observations were made in the Vlagtwedde-Vlaardingen study. Two overlapping SNPs in never-smokers (NFYC and FABP7) had the same direction of effect in the validation cohorts and the NFYC SNP was an eQTL for NFYC-AS1. NFYC is a transcription factor that binds to several known COPD genes, and FABP7 may be involved in abnormal pulmonary development. CONCLUSIONS: The definition of airflow obstruction and the population under study may be important determinants of which SNPs are associated with airflow obstruction. The genes FABP7 and NFYC(-AS1) could play a role in airflow obstruction in never-smokers specifically.

From blood to lung tissue: effect of cigarette smoke on DNA methylation and lung function.

BACKGROUND: Genetic and environmental factors play a role in the development of COPD. The epigenome, and more specifically DNA methylation, is recognized as important link between these factors. We postulate that DNA methylation is one of the routes by which cigarette smoke influences the development of COPD. In this study, we aim to identify CpG-sites that are associated with cigarette smoke exposure and lung function levels in whole blood and validate these CpG-sites in lung tissue. METHODS: The association between pack years and DNA methylation was studied genome-wide in 658 current smokers with >5 pack years using robust linear regression analysis. Using mediation analysis, we subsequently selected the CpG-sites that were also associated with lung function levels. Significant CpG-sites were validated in lung tissue with pyrosequencing and expression quantitative trait methylation (eQTM) analysis was performed to investigate the association between DNA methylation and gene expression. RESULTS: 15 CpG-sites were significantly associated with pack years and 10 of these were additionally associated with lung function levels. We validated 5 CpG-sites in lung tissue and found several associations between DNA methylation and gene expression. CONCLUSION: This study is the first to validate a panel of CpG-sites that are associated with cigarette smoking and lung function levels in whole blood in the tissue of interest: lung tissue.

Occupational exposure to pesticides is associated with differential DNA methylation.

OBJECTIVES: Occupational pesticide exposure is associated with a wide range of diseases, including lung diseases, but it is largely unknown how pesticides influence airway disease pathogenesis. A potential mechanism might be through epigenetic mechanisms, like DNA methylation. Therefore, we assessed associations between occupational exposure to pesticides and genome-wide DNA methylation sites. METHODS: 1561 subjects of LifeLines were included with either no (n=1392), low (n=108) or high (n=61) exposure to any type of pesticides (estimated based on current or last held job). Blood DNA methylation levels were measured using Illumina 450K arrays. Associations between pesticide exposure and 420 938 methylation sites (CpGs) were assessed using robust linear regression adjusted for appropriate confounders. In addition, we performed genome-wide stratified and interaction analyses by gender, smoking and airway obstruction status, and assessed associations between gene expression and methylation for genome-wide significant CpGs (n=2802). RESULTS: In total for all analyses, high pesticide exposure was genome-wide significantly (false discovery rate P<0.05) associated with differential DNA methylation of 31 CpGs annotated to 29 genes. Twenty of these CpGs were found in subjects with airway obstruction. Several of the identified genes, for example, RYR1, ALLC, PTPRN2, LRRC3B, PAX2 and VTRNA2-1, are genes previously linked to either pesticide exposure or lung-related diseases. Seven out of 31 CpGs were associated with gene expression levels. CONCLUSIONS: We show for the first time that occupational exposure to pesticides is genome-wide associated with differential DNA methylation. Further research should reveal whether this differential methylation plays a role in the airway disease pathogenesis induced by pesticides.

Understanding the role of the chromosome 15q25.1 in COPD through epigenetics and transcriptomics.

Chronic obstructive pulmonary disease (COPD) is a major health burden in adults and cigarette smoking is considered the most important environmental risk factor of COPD. Chromosome 15q25.1 locus is associated with both COPD and smoking. Our study aims at understanding the mechanism underlying the association of chromosome 15q25.1 with COPD through epigenetic and transcriptional variation in a population-based setting. To assess if COPD-associated variants in 15q25.1 are methylation quantitative trait loci, epigenome-wide association analysis of four genetic variants, previously associated with COPD (P < 5 × 10-8) in the 15q25.1 locus (rs12914385:C>T-CHRNA3, rs8034191:T>C-HYKK, rs13180:C>T-IREB2 and rs8042238:C>T-IREB2), was performed in the Rotterdam study (n = 1489). All four variants were significantly associated (P < 1.4 × 10-6) with blood DNA methylation of IREB2, CHRNA3 and PSMA4, of which two, including IREB2 and PSMA4, were also differentially methylated in COPD cases and controls (P < 0.04). Further additive and multiplicative effects of smoking were evaluated and no significant effect was observed. To evaluate if these four genetic variants are expression quantitative trait loci, transcriptome-wide association analysis was performed in 1087 lung samples. All four variants were also significantly associated with differential expression of the IREB2 3'UTR in lung tissues (P < 5.4 × 10-95). We conclude that regulatory mechanisms affecting the expression of IREB2 gene, such as DNA methylation, may explain the association between genetic variants in chromosome 15q25.1 and COPD, largely independent of smoking.

COPD GWAS variant at 19q13.2 in relation with DNA methylation and gene expression.

Chronic obstructive pulmonary disease (COPD) is among the major health burdens in adults. While cigarette smoking is the leading risk factor, a growing number of genetic variations have been discovered to influence disease susceptibility. Epigenetic modifications may mediate the response of the genome to smoking and regulate gene expression. Chromosome 19q13.2 region is associated with both smoking and COPD, yet its functional role is unclear. Our study aimed to determine whether rs7937 (RAB4B, EGLN2), a top genetic variant in 19q13.2 region identified in genome-wide association studies of COPD, is associated with differential DNA methylation in blood (N = 1490) and gene expression in blood (N = 721) and lungs (N = 1087). We combined genetic and epigenetic data from the Rotterdam Study (RS) to perform the epigenome-wide association analysis of rs7937. Further, we used genetic and transcriptomic data from blood (RS) and from lung tissue (Lung expression quantitative trait loci mapping study), to perform the transcriptome-wide association study of rs7937. Rs7937 was significantly (FDR < 0.05) and consistently associated with differential DNA methylation in blood at 4 CpG sites in cis, independent of smoking. One methylation site (cg11298343-EGLN2) was also associated with COPD (P = 0.001). Additionally, rs7937 was associated with gene expression levels in blood in cis (EGLN2), 42% mediated through cg11298343, and in lung tissue, in cis and trans (NUMBL, EGLN2, DNMT3A, LOC101929709 and PAK2). Our results suggest that changes of DNA methylation and gene expression may be intermediate steps between genetic variants and COPD, but further causal studies in lung tissue should confirm this hypothesis.

The emerging landscape of dynamic DNA methylation in early childhood.

BACKGROUND: DNA methylation has been found to associate with disease, aging and environmental exposure, but it is unknown how genome, environment and disease influence DNA methylation dynamics in childhood. RESULTS: By analysing 538 paired DNA blood samples from children at birth and at 4-5 years old and 726 paired samples from children at 4 and 8 years old from four European birth cohorts using the Illumina Infinium Human Methylation 450 k chip, we have identified 14,150 consistent age-differential methylation sites (a-DMSs) at epigenome-wide significance of p < 1.14 × 10-7. Genes with an increase in age-differential methylation were enriched in pathways related to 'development', and were more often located in bivalent transcription start site (TSS) regions, which can silence or activate expression of developmental genes. Genes with a decrease in age-differential methylation were involved in cell signalling, and enriched on H3K27ac, which can predict developmental state. Maternal smoking tended to decrease methylation levels at the identified da-DMSs. We also found 101 a-DMSs (0.71%) that were regulated by genetic variants using cis-differential Methylation Quantitative Trait Locus (cis-dMeQTL) mapping. Moreover, a-DMS-associated genes during early development were significantly more likely to be linked with disease. CONCLUSION: Our study provides new insights into the dynamic epigenetic landscape of the first 8 years of life.

Genome-wide association study on the FEV1/FVC ratio in never-smokers identifies HHIP and FAM13A.

BACKGROUND: Although a striking proportion (25% to 45%) of patients with chronic obstructive pulmonary disease are never-smokers, most genetic susceptibility studies have not focused on this group exclusively. OBJECTIVE: The aim of this study was to identify common genetic variants associated with FEV1 and its ratio to forced vital capacity (FVC) in never-smokers. METHODS: Genome-wide association studies were performed in 5070 never-smokers of the identification cohort LifeLines, and results (P < 10-5) were verified by using a meta-analysis of the Vlagtwedde-Vlaardingen study and the Rotterdam Study I-III (total n = 1966). Furthermore, we aimed to assess the effects of the replicated variants in more detail by performing genetic risk score, expression quantitative trait loci, and variant*ever-smoking interaction analyses. RESULTS: We identified associations between the FEV1/FVC ratio and 5 common genetic variants in the identification cohort, and 2 of these associations were replicated. The 2 variants annotated to the genes hedgehog interacting protein (HHIP) and family with sequence similarity 13 member A (FAM13A) were shown to have an additive effect on FEV1/FVC levels in the genetic risk score analysis; were associated with gene expression of HHIP and FAM13A in lung tissue, respectively; and were genome-wide significant in a meta-analysis including both identification and 4 verification cohorts (P < 2.19 × 10-7). Finally, we did not identify significant interactions between the variants and ever smoking. Results of the FEV1 identification analysis were not replicated. CONCLUSION: The genes HHIP and FAM13A confer a risk for airway obstruction in general that is not driven exclusively by cigarette smoking, which is the main risk factor for chronic obstructive pulmonary disease.

Susceptibility to chronic mucus hypersecretion, a genome wide association study.

BACKGROUND: Chronic mucus hypersecretion (CMH) is associated with an increased frequency of respiratory infections, excess lung function decline, and increased hospitalisation and mortality rates in the general population. It is associated with smoking, but it is unknown why only a minority of smokers develops CMH. A plausible explanation for this phenomenon is a predisposing genetic constitution. Therefore, we performed a genome wide association (GWA) study of CMH in Caucasian populations. METHODS: GWA analysis was performed in the NELSON-study using the Illumina 610 array, followed by replication and meta-analysis in 11 additional cohorts. In total 2,704 subjects with, and 7,624 subjects without CMH were included, all current or former heavy smokers (≥20 pack-years). Additional studies were performed to test the functional relevance of the most significant single nucleotide polymorphism (SNP). RESULTS: A strong association with CMH, consistent across all cohorts, was observed with rs6577641 (p = 4.25×10(-6), OR = 1.17), located in intron 9 of the special AT-rich sequence-binding protein 1 locus (SATB1) on chromosome 3. The risk allele (G) was associated with higher mRNA expression of SATB1 (4.3×10(-9)) in lung tissue. Presence of CMH was associated with increased SATB1 mRNA expression in bronchial biopsies from COPD patients. SATB1 expression was induced during differentiation of primary human bronchial epithelial cells in culture. CONCLUSIONS: Our findings, that SNP rs6577641 is associated with CMH in multiple cohorts and is a cis-eQTL for SATB1, together with our additional observation that SATB1 expression increases during epithelial differentiation provide suggestive evidence that SATB1 is a gene that affects CMH.

Correlation of genetic risk and messenger RNA expression in a Th17/IL23 pathway analysis in inflammatory bowel disease.

BACKGROUND: The Th17/IL23 pathway has both genetically and biologically been implicated in the pathogenesis of the inflammatory bowel diseases (IBD), Crohn's disease, and ulcerative colitis. So far, it is unknown whether and how associated risk variants affect expression of the genes encoding for Th17/IL23 pathway proteins. METHODS: Ten IBD-associated SNPs residing near Th17/IL23 genes were used to construct a genetic risk model in 753 Dutch IBD cases and 1045 controls. In an independent cohort of 40 Crohn's disease, 40 ulcerative colitis, and 40 controls, the genetic risk load and presence of IBD were correlated to quantitative PCR-generated messenger RNA (mRNA) expression of 9 representative Th17/IL23 genes in both unstimulated and PMA/CaLo stimulated peripheral blood mononuclear cells. In 1240 individuals with various immunological diseases with whole genome genotype and mRNA-expression data, we also assessed correlation between genetic risk load and differential mRNA expression and sought for SNPs affecting expression of all currently known Th17/IL23 pathway genes (cis-expression quantitative trait locus). RESULTS: The presence of IBD, but not the genetic risk load, was correlated to differential mRNA expression for IL6 in unstimulated peripheral blood mononuclear cells and to IL23A and RORC in response to stimulation. The cis-expression quantitative trait locus analysis showed little evidence for correlation between genetic risk load and mRNA expression of Th17/IL23 genes, because we identified for only 2 of 22 Th17/IL23 genes a cis-expression quantitative trait locus single nucleotide polymorphism that is also associated to IBD (STAT3 and CCR6). CONCLUSIONS: Our results suggest that only the presence of IBD and not the genetic risk load alters mRNA expression levels of IBD-associated Th17/IL23 genes.

ADAM33 gene polymorphisms and mortality. A prospective cohort study.

The ADAM33 gene is associated with the pathophysiology of Chronic Obstructive Pulmonary Disease (COPD) and atherosclerosis. In this study we investigated all-cause, COPD and cardiovascular mortality, in relation to single nucleotide polymorphisms (SNPs) in ADAM33 (Q_1, S_1, S_2, T_1 and T_2) that were genotyped in 1,390 subjects from the Vlagtwedde/Vlaardingen cohort. Participants were examined at entry in 1989/1990 and followed up till evaluation of the vital status on December 31(st), 2008. Using Cox proportional hazards regression we estimated the risk of the SNPs in relation to mortality, adjusting for gender, age, FEV1, height, place of residence and packyears of smoking. Additionally, we performed stratified analyses according to gender and smoking habits. After 18 years, 284 (20.4%) subjects had died (107 due to cardiovascular disease and 20 due to COPD). Individuals homozygous for the minor allele of SNP T_2 had an increased risk of all-cause and cardiovascular mortality compared to wild types: hazard ratio 3.6 (95% confidence interval 2.0 to 6.7) and 3.4 (1.2 to 9.5) respectively. Individuals homozygous for the minor allele of S_1, S_2, T_2 or Q_1 had a significantly increased risk of COPD mortality. In stratified analyses the risk of all-cause mortality associated with SNP T_2 did not change: females 3.5 (1.5 to 8.3), males 3.1 (1.2 to 7.6), never smokers 3.8 (0.9 to 16.3), ever smokers 3.6 (1.8 to 7.2). This study shows for the first time that ADAM33 is a pleiotropic gene that is associated with all-cause, COPD and cardiovascular mortality, independent of potential confounders.

Limited evidence for parent-of-origin effects in inflammatory bowel disease associated loci.

BACKGROUND: Genome-wide association studies of two main forms of inflammatory bowel diseases (IBD), Crohn's disease (CD) and ulcerative colitis (UC), have identified 99 susceptibility loci, but these explain only 23% of the genetic risk. Part of the 'hidden heritability' could be in transmissible genetic effects in which mRNA expression in the offspring depends on the parental origin of the allele (genomic imprinting), since children whose mothers have CD are more often affected than children with affected fathers. We analyzed parent-of-origin (POO) effects in Dutch and Indian cohorts of IBD patients. METHODS: We selected 28 genetic loci associated with both CD and UC, and tested them for POO effects in 181 Dutch IBD case-parent trios. Three susceptibility variants in NOD2 were tested in 111 CD trios and a significant finding was re-evaluated in 598 German trios. The UC-associated gene, BTNL2, reportedly imprinted, was tested in 70 Dutch UC trios. Finally, we used 62 independent Indian UC trios to test POO effects of five established Indian UC risk loci. RESULTS: We identified POO effects for NOD2 (L1007fs; OR = 21.0, P-value = 0.013) for CD; these results could not be replicated in an independent cohort (OR = 0.97, P-value = 0.95). A POO effect in IBD was observed for IL12B (OR = 3.2, P-value = 0.019) and PRDM1 (OR = 5.6, P-value = 0.04). In the Indian trios the IL10 locus showed a POO effect (OR = 0.2, P-value = 0.03). CONCLUSIONS: Little is known about the effect of genomic imprinting in complex diseases such as IBD. We present limited evidence for POO effects for the tested IBD loci. POO effects explain part of the hidden heritability for complex genetic diseases but need to be investigated further.

Analysis of SNPs with an effect on gene expression identifies UBE2L3 and BCL3 as potential new risk genes for Crohn's disease.

Genome-wide association studies (GWAS) for Crohn's disease (CD) have identified loci explaining approximately 20% of the total genetic risk of CD. Part of the other genetic risk loci is probably partly hidden among signals discarded by the multiple testing correction needed in the analysis of GWAS data. Strategies for finding these hidden loci require large replication cohorts and are costly to perform. We adopted a strategy of selecting SNPs for follow-up that showed a correlation to gene expression [cis-expression quantitative trait loci (eQTLs)] since these have been shown more likely to be trait-associated. First we show that there is an overrepresentation of cis-eQTLs in the known CD-associated loci. Then SNPs were selected for follow-up by screening the top 500 SNP hits from a CD GWAS data set. We identified 10 cis-eQTL SNPs. These 10 SNPs were tested for association with CD in two independent cohorts of Dutch CD patients (1539) and healthy controls (2648). In a combined analysis, we identified two cis-eQTL SNPs that were associated with CD rs2298428 in UBE2L3 (P=5.22x10(-5)) and rs2927488 in BCL3 (P=2.94x10(-4)). After adding additional publicly available data from a previously reported meta-analysis, the association with rs2298428 almost reached genome-wide significance (P=2.40x10(-7)) and the association with rs2927488 was corroborated (P=6.46x10(-4)). We have identified UBE2L3 and BCL3 as likely novel risk genes for CD. UBE2L3 is also associated with other immune-mediated diseases. These results show that eQTL-based pre-selection for follow-up is a useful approach for identifying risk loci from a moderately sized GWAS.

Detecting shared pathogenesis from the shared genetics of immune-related diseases.

Recent genetic studies have revealed shared immunological mechanisms in several immune-related disorders that further our understanding of the development and concomitance of these diseases. Our Review focuses on these shared aspects, using the novel findings of recently performed genome-wide association studies and non-synonymous SNP scans as a starting point. We discuss how identifying new genes that are associated with more than one autoimmune or chronic inflammatory disorder could explain the genetic basis of the shared pathogenesis of immune-related diseases. This analysis helps to highlight the key molecular pathways that are involved in these disorders and the potential roles of novel genes in immune-related diseases.

Genetic analysis of innate immunity in Crohn's disease and ulcerative colitis identifies two susceptibility loci harboring CARD9 and IL18RAP.

The two main phenotypes of inflammatory bowel disease (IBD)--Crohn's disease (CD) and ulcerative colitis (UC)--are chronic intestinal inflammatory disorders with a complex genetic background. Using a three-stage design, we performed a functional candidate-gene analysis of innate immune pathway in IBD. In phase I, we typed 354 SNPs from 85 innate immunity genes in 520 Dutch IBD patients (284 CD, 236 UC) and 808 controls. In phase II, ten autosomal SNPs showing association at p < 0.006 in phase I were replicated in a second cohort of 545 IBD patients (326 CD, 219 UC) and 360 controls. In phase III, four SNPs with p < 0.01 in the combined phase I and phase II analysis were genotyped in an additional 786 IBD samples (452 CD, 334 UC) and 768 independent controls. Joint analysis of 1851 IBD patients (1062 CD, 789 UC) and 1936 controls demonstrated strong association to the IL18RAP rs917997 SNP for both CD and UC (p(IBD) 1.9 x 10(-8); OR 1.35). Association in CD is independently supported by the Crohn's disease dataset of the Wellcome Trust Case Control Consortium (imputed SNP rs917997, p = 9.19 x 10(-4)). In addition, an association of the CARD9 rs10870077 SNP to CD and UC was observed (p(IBD) = 3.25 x 10(-5); OR 1.21). Both genes are located in extended haplotype blocks on 2q11-2q12 and 9q34.3, respectively. Our results indicate two IBD loci and further support the importance of the innate immune system in the predisposition to both CD and UC.