The phenotypic effect of some single nucleotide polymorphisms (SNPs) depends on their parental origin. We present a novel approach to detect parent-of-origin effects (POEs) in genome-wide genotype data of unrelated individuals. The method exploits increased phenotypic variance in the heterozygous genotype group relative to the homozygous groups. We applied the method to >56,000 unrelated individuals to search for POEs influencing body mass index (BMI). Six lead SNPs were carried forward for replication in five family-based studies (of ∼4,000 trios). Two SNPs replicated: the paternal rs2471083-C allele (located near the imprinted KCNK9 gene) and the paternal rs3091869-T allele (located near the SLC2A10 gene) increased BMI equally (beta = 0.11 (SD), P<0.0027) compared to the respective maternal alleles. Real-time PCR experiments of lymphoblastoid cell lines from the CEPH families showed that expression of both genes was dependent on parental origin of the SNPs alleles (P<0.01). Our scheme opens new opportunities to exploit GWAS data of unrelated individuals to identify POEs and demonstrates that they play an important role in adult obesity. Large genetic association studies have revealed many genetic factors influencing common traits, such as body mass index (BMI). These studies assume that the effect of genetic variants is the same regardless of whether they are inherited from the mother or the father. In our study, we have developed a new approach that allows us to investigate variants whose impact depends on their parental origin (parent-of-origin effects), in unrelated samples when the parental origin cannot be inferred. This is feasible because at genetic markers at which such effects occur there is increased variability of the trait among individuals who inherited different genetic codes from their mother and their father compared to individuals who inherited the same genetic code from both parents. We applied this methodology to discover genetic markers with parent-of-origin effects (POEs) on BMI. This resulted in six candidate markers showing strong POE association. We then attempted to replicate the POE effects of these markers in family studies (where one can infer the parental origin of the inherited variants). Two of our candidates showed significant association in the family studies, the paternal and maternal effects of these markers were in the opposite direction.

DNA methylation profiles are responsive to environmental stimuli and metabolic shifts. This makes DNA methylation a potential biomarker of environmental-related and lifestyle-driven diseases of adulthood. Therefore, we investigated if white blood cells’ (WBCs) DNA methylation profiles are associated with myocardial infarction (MI) occurrence. Whole-genome DNA methylation was investigated by microarray analysis in 292 MI cases and 292 matched controls from the large prospective Italian European Prospective Investigation into Cancer and Nutrition (EPIC) cohort (EPICOR study). Significant signals (false discovery rate (FDR) adjusted P < 0.05) were replicated by mass spectrometry in 317 MI cases and 262 controls from the Dutch EPIC cohort (EPIC-NL). Long interspersed nuclear element-1 (LINE-1) methylation profiles were also evaluated in both groups. A differentially methylated region (DMR) within the zinc finger and BTB domain-containing protein 12 (ZBTB12) gene body and LINE-1 hypomethylation were identified in EPICOR MI cases and replicated in the EPIC-NL sample (ZBTB12-DMR meta-analysis: effect size ± se = −0.016 ± 0.003, 95 % CI = −0.021;−0.011, P = 7.54 × 10−10; LINE-1 methylation meta-analysis: effect size ± se = −0.161 ± 0.040, 95 % CI = −0.239;−0.082, P = 6.01 × 10−5). Moreover, cases with shorter time to disease had more pronounced ZBTB12-DMR hypomethylation (meta-analysis, men: effect size ± se = −0.0059 ± 0.0017, P TREND = 5.0 × 10−4; women: effect size ± se = −0.0053 ± 0.0019, P TREND = 6.5 × 10−3) and LINE-1 hypomethylation (meta-analysis, men: effect size ± se = −0.0010 ± 0.0003, P TREND = 1.6 × 10−3; women: effect size ± se = −0.0008 ± 0.0004, P TREND = 0.026) than MI cases with longer time to disease. In the EPIC-NL replication panel, DNA methylation profiles improved case-control discrimination and reclassification when compared with traditional MI risk factors only (net reclassification improvement (95 % CI) between 0.23 (0.02–0.43), P = 0.034, and 0.89 (0.64–1.14), P < 1 × 10−5). Our data suggest that specific methylation profiles can be detected in WBCs, in a preclinical condition, several years before the occurrence of MI, providing an independent signature of cardiovascular risk. We showed that prediction accuracy can be improved when DNA methylation is taken into account together with traditional MI risk factors, although further confirmation on a larger sample is warranted. Our findings support the potential use of DNA methylation patterns in peripheral blood white cells as promising early biomarkers of MI.

Worldwide, hypertension still represents a serious health burden with nine million people dying as a consequence of hypertension-related complications. Essential hypertension is a complex trait supported by multifactorial genetic inheritance together with environmental factors. The heritability of blood pressure (BP) is estimated to be 30–50%. A great effort was made to find genetic variants affecting BP levels through Genome-Wide Association Studies (GWAS). This approach relies on the “common disease–common variant” hypothesis and led to the identification of multiple genetic variants which explain, in aggregate, only 2–3% of the genetic variance of hypertension. Part of the missing genetic information could be caused by variants too rare to be detected by GWAS. The use of exome chips and Next-Generation Sequencing facilitated the discovery of causative variants. Here, we report the advances in the detection of novel rare variants, genes, and/or pathways through the most promising approaches, and the recent statistical tests that have emerged to handle rare variants. We also discuss the need to further support rare novel variants with replication studies within larger consortia and with deeper functional studies to better understand how new genes might improve patient care and the stratification of the response to antihypertensive treatments.

Asbestos exposure is the main risk factor for malignant pleural mesothelioma (MPM), a rare aggressive tumor. Nevertheless, only 5–17% of those exposed to asbestos develop MPM, suggesting the involvement of other environmental and genetic risk factors. To identify the genetic risk factors that may contribute to the development of MPM, we conducted a genome-wide association study (GWAS; 370,000 genotyped SNPs, 5 million imputed SNPs) in Italy, among 407 MPM cases and 389 controls with a complete history of asbestos exposure. A replication study was also undertaken and included 428 MPM cases and 1269 controls from Australia. Although no single marker reached the genome-wide significance threshold, several associations were supported by haplotype-, chromosomal region-, gene- and gene-ontology process-based analyses. Most of these SNPs were located in regions reported to harbor aberrant alterations in mesothelioma (SLC7A14, THRB, CEBP350, ADAMTS2, ETV1, PVT1 and MMP14 genes), causing at most a 2–3-fold increase in MPM risk. The Australian replication study showed significant associations in five of these chromosomal regions (3q26.2, 4q32.1, 7p22.2, 14q11.2, 15q14). Multivariate analysis suggested an independent contribution of 10 genetic variants, with an Area Under the ROC Curve (AUC) of 0.76 when only exposure and covariates were included in the model, and of 0.86 when the genetic component was also included, with a substantial increase of asbestos exposure risk estimation (odds ratio, OR: 45.28, 95% confidence interval, CI: 21.52–95.28). These results showed that genetic risk factors may play an additional role in the development of MPM, and that these should be taken into account to better estimate individual MPM risk in individuals who have been exposed to asbestos.

Genome-wide association studies (GWAS) of urinary bladder cancer (UBC) have yielded common variants at 12 loci that associate with risk of the disease. We report here the results of a GWAS of UBC including 1670 UBC cases and 90 180 controls, followed by replication analysis in additional 5266 UBC cases and 10 456 controls. We tested a dataset containing 34.2 million variants, generated by imputation based on whole-genome sequencing of 2230 Icelanders. Several correlated variants at 20p12, represented by rs62185668, show genome-wide significant association with UBC after combining discovery and replication results (OR = 1.19, P = 1.5 × 10⁻¹¹ for rs62185668-A, minor allele frequency = 23.6%). The variants are located in a non-coding region approximately 300 kb upstream from the JAG1 gene, an important component of the Notch signaling pathways that may be oncogenic or tumor suppressive in several forms of cancer. Our results add to the growing number of UBC risk variants discovered through GWAS.