Genetic insights into resting heart rate and its role in cardiovascular disease

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for IgA nephropathy

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. Here we performed a genome-wide association study of 10,146 kidney-biopsy-diagnosed IgAN cases and 28,751 controls across 17 international cohorts. We defined 30 genome-wide significant risk loci explaining 11% of disease risk. A total of 16 loci were new, including TNFSF4/TNFSF18, REL, CD28, PF4V1, LY86, LYN, ANXA3, TNFSF8/TNFSF15, REEP3, ZMIZ1, OVOL1/RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The risk loci were enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. We also observed a positive genetic correlation between IgAN and serum IgA levels. High polygenic score for IgAN was associated with earlier onset of kidney failure. In a comprehensive functional annotation analysis of candidate causal genes, we observed convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

Dissecting the genetic heterogeneity of gastric cancer

BACKGROUND

Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture.

METHODS

We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO.

FINDINGS

Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level.

INTERPRETATION

Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO.

FUNDING

German Research Foundation (DFG).

Modulation of gut microbiota: The effects of a fruits and vegetables supplement

The consumption of an optimal amount of fruits and vegetables is known to improve physical fitness and physiological body functions. Healthy eating habits, including intake of fruits and vegetables, can modify gut microbiota. This study aimed to demonstrate the effectiveness of a formulated fruit and vegetable supplement (FVS) in modulating the antioxidant capacity and the gut microbiota composition. We enrolled 30 healthy volunteer subjects, matched for age, gender, BMI, and smoking habits, and randomized them into the FVS and the placebo (PLA) groups. Among the serum vitamins, the folic acid level was significantly higher (p = 0.001) in the FVS group than in the PLA group, whereas the vitamin B2 level was significantly higher in the PLA group than in the FVS group (p = 0.028). The antioxidant capacity, measured by using the oxygen radical absorbance capacity (ORAC) method, was also slightly higher in the FVS group than in the PLA group but did not reach statistical significance. The dietary intake, assessed by 24-h recalls, did not show any significant changes after the supplementation in both the groups. The gut microbiome composition, measured by 16S rDNA sequencing, showed no difference in both alpha and beta diversities, whereas the LEfse analysis revealed a microbial shift after the treatment, with a decreased abundance of the genus Ruminococcus from the Lachnospiraceae family (p = 0.009), and the unclassified genus from the family Erysipelotrichaceae (UC36, p = 0.003) in the FVS group compared with the PLA group (confirmed by SIAMCAT analysis, AUC = 74.1%). With a minor effect, the genus Faecalibacterium and unclassified genus and family from the order Lactobacillales (UC31) were also increased in the FVS group compared with the PLA group (p = 0.0474, p = 0.0352, respectively). SCFA measurement by gas chromatography–mass spectrometry showed an increased level of 2-methylbutyrate in the FVS group compared with the PLA group (p = 0.0385). Finally, the Spearman correlation analysis showed that in the FVS group, the genus Faecalibacterium positively correlated with 2-methyl butyrate (p = 0.040). In the PLA group, none of the significant bacteria correlated with either SCFA or serum biomarkers. The network analysis confirmed the positive correlation between genus Faecalibacterium and 2-methyl butyrate. We can conclude that the FVS in healthy individuals modified the gut microbiota composition and metabolites, and it can potentially contribute to reduce the pro-inflammatory response along with the antioxidant capacity.

Identification of risk loci for primary aldosteronism in genome-wide association studies

Primary aldosteronism affects up to 10% of hypertensive patients and is responsible for treatment resistance and increased cardiovascular risk. Here we perform a genome-wide association study in a discovery cohort of 562 cases and 950 controls and identify three main loci on chromosomes 1, 13 and X; associations on chromosome 1 and 13 are replicated in a second cohort and confirmed by a meta-analysis involving 1162 cases and 3296 controls. The association on chromosome 13 is specific to men and stronger in bilateral adrenal hyperplasia than aldosterone producing adenoma. Candidate genes located within the two loci, CASZ1 and RXFP2, are expressed in human and mouse adrenals in different cell clusters. Their overexpression in adrenocortical cells suppresses mineralocorticoid output under basal and stimulated conditions, without affecting cortisol biosynthesis. Our study identifies the first risk loci for primary aldosteronism and highlights new mechanisms for the development of aldosterone excess.

Detection of primary aldosteronism, the most common form of secondary arterial hypertension, is essential for targeted management and prevention of cardiovascular complications. Here, the authors identify genetic loci associated with primary aldosteronism, suggesting new mechanisms of disease.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

A multi-step genomic approach prioritized TBKBP1 gene as relevant for multiple sclerosis susceptibility

Background

Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the “missing heritability” phenomenon.

Objective

To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease.

Methods

We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed.

Results

Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant.

Conclusions

No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00415-022-11109-8.

Association of colorectal cancer with genetic and epigenetic variation in PEAR1—A population-based cohort study

Platelet Endothelial Aggregation Receptor 1 (PEAR1) modulates angiogenesis and platelet contact-induced activation, which play a role in the pathogenesis of colorectal cancer. We therefore tested the association of incident colorectal cancer and genetic and epigenetic variability in PEAR1 among 2532 randomly recruited participants enrolled in the family-based Flemish Study on Environment, Genes and Health Outcomes (51.2% women; mean age 44.8 years). All underwent genotyping of rs12566888 located in intron 1 of the PEAR1 gene; in 926 participants, methylation at 16 CpG sites in the PEAR1 promoter was also assessed. Over 18.1 years (median), 49 colorectal cancers occurred, all in different pedigrees. While accounting for clustering of risk factors within families and adjusting for sex, age, body mass index, the total-to-HDL cholesterol ratio, serum creatinine, plasma glucose, smoking and drinking, use of antiplatelet and nonsteroidal anti-inflammatory drug, the hazard ratio of colorectal cancer contrasting minor-allele (T) carriers vs. major-allele (GG) homozygotes was 2.17 (95% confidence interval, 1.18–3.99; P = 0.013). Bootstrapped analyses, from which we randomly excluded from two to nine cancer cases, provided confirmatory results. In participants with methylation data, we applied partial least square discriminant analysis (PLS-DA) and identified two methylation sites associated with higher colorectal cancer risk and two with lower risk. In-silico analysis suggested that methylation of the PEAR1 promoter at these four sites might affect binding of transcription factors p53, PAX5, and E2F-1, thereby modulating gene expression. In conclusion, our findings suggest that genetic and epigenetic variation in PEAR1 modulates the risk of colorectal cancer in white Flemish. To what extent, environmental factors as exemplified by our methylation data, interact with genetic predisposition and modulate penetrance of colorectal cancer risk is unknown.

Genome-wide association analyses identify new Brugada syndrome risk loci and highlight a new mechanism of sodium channel regulation in disease susceptibility

Brugada syndrome (BrS) is a cardiac arrhythmia disorder associated with sudden death in young adults. With the exception of SCN5A, encoding the cardiac sodium channel Na_V1.5, susceptibility genes remain largely unknown. Here we performed a genome-wide association meta-analysis comprising 2,820 unrelated cases with BrS and 10,001 controls, and identified 21 association signals at 12 loci (10 new). Single nucleotide polymorphism (SNP)-heritability estimates indicate a strong polygenic influence. Polygenic risk score analyses based on the 21 susceptibility variants demonstrate varying cumulative contribution of common risk alleles among different patient subgroups, as well as genetic associations with cardiac electrical traits and disorders in the general population. The predominance of cardiac transcription factor loci indicates that transcriptional regulation is a key feature of BrS pathogenesis. Furthermore, functional studies conducted on MAPRE2, encoding the microtubule plus-end binding protein EB2, point to microtubule-related trafficking effects on Na_V1.5 expression as a new underlying molecular mechanism. Taken together, these findings broaden our understanding of the genetic architecture of BrS and provide new insights into its molecular underpinnings.

Dissecting the Polygenic Basis of Primary Hypertension: Identification of Key Pathway-Specific Components

Introduction and Objectives

Genome-wide association studies have identified a high number of genetic loci associated with hypertension suggesting the presence of an underlying polygenic architecture. In this study, we aimed to dissect the polygenic component of primary hypertension searching also for pathway-specific components.

Methods

The polygenic risk score (PRS) models, based on the UK biobank genetic signals for hypertension status, were obtained on a target Italian case/control cohort including 561 cases and 731 hyper-normal controls from HYPERGENES, and were then applied to an independent validation cohort composed by multi-countries European-based samples including 1,284 cases and 960 hyper-normal controls.

Results

The resulting genome-wide PRS was capable of stratifying the individuals for hypertension risk by comparing between individuals in the last PRS decile and the median decile: we observed an odds ratio (OR) of 3.62, CI = [2.01, 6.32] (P = 9.01E-07) and 3.22, 95% CI = [2.06, 5.10] (P = 6.47E-08) in the target and validation cohorts, respectively. The relatively high case/control ORs across PRS quantiles corroborates the presence of strong polygenic components which could be driven by an enrichment of risk alleles within the cases but also by potential enrichment of protective alleles in the old normotensive controls. Moreover, novel pathway-specific PRS revealed an enrichment of the polygenic signal attributable to specific biological pathways. Among those the most significantly associated with hypertension status was the calcium signaling pathway together with other mainly related such as the phosphatidylinositol/inositol phosphate pathways.

Conclusions

The development of pathway-specific PRS could prioritize biological mechanisms, according to their contribution to the genetic susceptibility, whose regulations might be a potential pharmacological preventive target.

Antihypertensive treatment guided by genetics: PEARL-HT, the randomized proof-of-concept trial comparing rostafuroxin with losartan

We compared a standard antihypertensive losartan treatment with a pharmacogenomics-guided rostafuroxin treatment in never-treated Caucasian and Chinese patients with primary hypertension. Rostafuroxin is a digitoxigenin derivative that selectively disrupts the binding to the cSrc-SH2 domain of mutant α-adducin and of the ouabain-activated Na-K pump at 10^-11 M. Of 902 patients screened, 172 were enrolled in Italy and 107 in Taiwan. After stratification for country and genetic background, patients were randomized to rostafuroxin or losartan, being the difference in the fall in office systolic blood pressure (OSBP) after 2-month treatment the primary endpoint. Three pharmacogenomic profiles (P) were examined, considering: P1, adding to the gene variants included in the subsequent P2, the variants detected by post-hoc analysis of a previous trial; P2, variants of genes encoding enzymes for endogenous ouabain (EO) synthesis (LSS and HSD3B1), EO transport (MDR1/ABCB1), adducin (ADD1 and ADD3); P3, variants of the LSS gene only. In Caucasians, the group differences (rostafuroxin 50 μg minus losartan 50 mg in OSBP mmHg) were significant both in P2 adjusted for genetic heterogeneity (P2a) and P3 LSS rs2254524 AA [9.8 (0.6-19.0), P = 0.038 and 13.4 (25.4-2.5), P = 0.031, respectively]. In human H295R cells transfected with LSS A and LSS C variants, the EO production was greater in the former (P = 0.038); this difference was abolished by rostafuroxin at 10^-11 M. Chinese patients had a similar drop in OSBP to Caucasians with losartan but no change in OSBP with rostafuroxin. These results show that genetics may guide drug treatment for primary hypertension in Caucasians.

Genome sequencing analysis identifies new loci associated with Lewy body dementia and provides insights into its genetic architecture

The genetic basis of Lewy body dementia (LBD) is not well understood. Here, we performed whole-genome sequencing in large cohorts of LBD cases and neurologically healthy controls to study the genetic architecture of this understudied form of dementia, and to generate a resource for the scientific community. Genome-wide association analysis identified five independent risk loci, whereas genome-wide gene-aggregation tests implicated mutations in the gene GBA. Genetic risk scores demonstrate that LBD shares risk profiles and pathways with Alzheimer's disease and Parkinson's disease, providing a deeper molecular understanding of the complex genetic architecture of this age-related neurodegenerative condition.

Abstract P118: Antihypertensive Treatment Guided By Genetics: Pearl-ht, The Randomized Proof-of-concept Trial Comparing Rostafuroxin With Losartan

Primary hypertension contributes to more than 50% of all the worldwide cardiovascular deaths. Current antihypertensive therapy is still targeting the physiological mechanisms regulating blood pressure, but not the “causal” ones. So far, genetic studies have been unable to prove gene causation in human primary hypertension. In this study we compared a standard antihypertensive losartan treatment with a pharmacogenomics-guided rostafuroxin treatment in never-treated Caucasian and Chinese patients with primary hypertension (PEARL-HT trial registration: EudraCT: 2010-022073-34, ClinicalTrials.gov: NCT01320397). Rostafuroxin is a digitoxigenin derivative that selectively disrupts the binding to the cSrc-SH2 domain of mutant α-adducin and of the ouabain-activated Na-K pump at 10 -11 M. Of 902 patients screened, 172 were enrolled in Italy and 107 in Taiwan. After stratification for country and genetic background, patients were randomized to rostafuroxin or losartan, being the difference in the fall in Office Systolic Blood Pressure (OSBP) after two-month treatment the primary endpoint. Three pharmacogenomic profiles (P) were examined, considering: P1, adding to the gene variants included in the subsequent P2, the variants detected by post-hoc analysis of a previous trial; P2, variants of genes encoding enzymes for endogenous ouabain (EO) synthesis ( LSS and HSD3B1 ), EO transport ( MDR1/ABCB1 ), adducin ( ADD1 and ADD3 ); P3, variants of the LSS gene only. In Caucasians, the group differences (rostafuroxin 50 μg minus losartan 50 mg in OSBP mmHg) were significant both in P2 adjusted for genetic heterogeneity (P2a) and P3 LSS rs2254524 AA [9.8 (0.6-19.0), P =0.038 and 13.4 (25.4-2.5), P =0.031, respectively]. In human H295R cells transfected with LSS A and LSS C variants, the EO production was greater in the former ( P =0.038); this difference was abolished by rostafuroxin at 10 -11 M. Chinese patients had a similar drop in OSBP to Caucasians with losartan but no change in OSBP with rostafuroxin. These results show that the selective blockade of mutant adducin and endogenous ouabain pressor mechanisms by rostafuroxin improves the therapy of primary hypertension in in Caucasians carriers of the related gene variants.

Abstract P122: Lanosterol Synthase (LSS) Gene As Predictor Of Kidney Dysfunction In Hypertensive Patients

Recently we have proposed LSS as genotype-based risk stratification to predict accelerated eGFR decay in a prospective cohort study of essential HYP patients. We extract clinical data of a general population from HYPERGENES consortium. We also collect genetic data for LSS polymorphism that was already demonstrated as involved in kidney damage. A cohort of 3137 patients were selected. At different age, as expected, HYP status have a deep impact on eGFR. Indeed in young (age < 50ys) eGFR is higher in HYP vs control (86.2 ± 19.3 vs 82.7 ± 16.3 ml/min; p=0.04); vice-versa in elder (age > 65ys) HYP have a reduction in eGFR (55.7 ± 12.9 vs 71.0 ± 12.9 ml/min; p<0.001) and higher incidence of severe CKD (2.7% vs 0.1%; p=0.001). We do not observe a direct influence of LSS polymorphism of eGFR. When LSS is considered according to HYP status and age class it is possible to observe a preservation of age-associated eGFR reduction in normotensive patients (eGFR 91.4 vs 90.1 vs 81.8 at age classes for LSS risk allele); on the other side, LSS seems to enhance eGFR reduction associate to HYP status (eGFR 84.6 vs 75.6 vs 48.0 for LSS risk allele; p=0.032 after correction for covariates; see table for details and figure 1). Furthermore, our genomics lab just developed a LSS V642L Knock-in mouse model by CRISPR-Cas9 technique. The sequencing analysis confirmed the correct insertion of 643L mutation for both alleles in mice gene. The 642L LSS Knock-in mice at 1 month of age were viable, healthy and indistinguishable from the WT counterpart

Conclusion: Patients caring the risk allele of this specific LSS polymorphism seems to express hyper-filtration if compared to their counterparts.

Transethnic Genome-Wide Association Study Provides Insights in the Genetic Architecture and Heritability of Long QT Syndrome

Supplemental Digital Content is available in the text.

Long QT syndrome (LQTS) is a rare genetic disorder and a major preventable cause of sudden cardiac death in the young. A causal rare genetic variant with large effect size is identified in up to 80% of probands (genotype positive) and cascade family screening shows incomplete penetrance of genetic variants. Furthermore, a proportion of cases meeting diagnostic criteria for LQTS remain genetically elusive despite genetic testing of established genes (genotype negative). These observations raise the possibility that common genetic variants with small effect size contribute to the clinical picture of LQTS. This study aimed to characterize and quantify the contribution of common genetic variation to LQTS disease susceptibility.

Klotho Gene in Human Salt-Sensitive Hypertension

BACKGROUND AND OBJECTIVES

Hypertension is a common aging-related disorder. Salt intake is one of the main environmental factors contributing to the development of hypertension. Transgenic mice with one-half Klotho deficiency displayed a spontaneous BP increase and salt-sensitive hypertension in response to high sodium intake. Usually circulating levels of α-Klotho decrease with age, and this reduction may be stronger in patients with several aging-related diseases. This study aimed at exploring the association of Klotho with salt sensitivity in humans.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

The role of Klotho polymorphisms and α-Klotho serum levels was evaluated in patients with hypertension who were treatment naive and underwent an acute salt-sensitivity test (discovery n=673, intravenous 2 L of 0.9% saline in 2 hours). Salt sensitivity was defined as a mean BP increase of >4 mm Hg at the end of the infusion. A total of 32 single nucleotide polymorphisms in the Klotho gene (KL), previously identified with a genome-wide association study, were used in the genetic analysis and studied for a pressure-natriuresis relationship.

RESULTS

Of the patients with hypertension, 35% were classified as salt sensitive. The most relevant polymorphism associated with pressure natriuresis was the common missense single nucleotide polymorphism rs9536314, and the GG and GT genotypes were more represented among patients who were salt sensitive (P=0.001). Those carrying the G allele showed a less steep pressure-natriuresis relationship, meaning that a significant increase in mean BP was needed to excrete the same quantity of salt compared with patients who were salt resistant. KL rs9536314 also replicated the pressure-natriuresis association in an independent replication cohort (n=193) and in the combined analysis (n=866). There was an inverse relationship between circulating Klotho and mean BP changes after the saline infusion (r=-0.14, P=0.03). Moreover, circulating α-Klotho was directly related to kidney function at baseline eGFR (r=0.22, P<0.001).

CONCLUSIONS

KL rs9536314 is associated with salt-sensitive hypertension in patients with hypertension who are treatment naive. Moreover, circulating α-Klotho levels were mainly related to diastolic BP changes at the end of a salt load and to eGFR as an expression of kidney aging.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

P1576Association of colorectal cancer with genetic and epigenetic variation in PEAR1 - a population-based cohort study

Background

Platelet Endothelial Aggregation Receptor 1 (PEAR1) modulates angiogenesis and platelet contact-induced activation, which play a role in the pathogenesis of colorectal cancer.

Purpose

To study the association of colorectal cancer with genetic and epigenetic variation in PEAR1.

Methods

Among 2532 randomly recruited participants enrolled in the family-based Flemish Study on Environment, Genes and Health Outcomes (51.2% women; mean age 44.8 years), we recorded the incidence of colorectal cancer and genotyped SNP rs12566888 located in intron 1 of the PEAR1 gene. In 929 participants, we also measured the methylation at 16 CpG sites in the PEAR1 promoter. In multivariable-adjusted analyses, we contrasted the risk of colorectal cancer in minor-allele (T) carriers vs. major allele (GG) homozygotes. We applied partial least square discriminant analysis (PLS-DA) to identify methylation sites associated with colorectal cancer.

Results

Over a median follow-up of 18.1 years, 49 patients developed colorectal cancer. While accounting for clustering within families and adjusting for sex, age, body mass index, the total-to-HDL cholesterol ratio, serum creatinine, plasma glucose, smoking and drinking, use of antiplatelet and nonsteroidal anti-inflammatory drug, the hazard ratio contrasting minor allele carriers vs. major allele homozygotes was 2.17 (95% confidence interval, 1.18–3.99; P=0.013). Bootstrapped analyses, from which we randomly excluded from two to nine cancer cases, provided confirmatory results. PLS-DA identified two methylation sites in the PEAR1 promoter associated with higher colorectal cancer risk and two with lower risk. In-silico analysis suggested that methylation of the PEAR1 promoter at these four sites affects binding of the transcription factors p53, PAX5, and E2F-1, thereby modulating gene expression.

Potential pathways

Conclusions

Our findings suggest that genetic and epigenetic variation in PEAR1 modulates the risk of colorectal cancer in white Flemish.

Acknowledgement/Funding

The European Union, the European Research Council, the European Research Area Net for Cardiovascular Diseases, and the Research Foundation Flanders.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

A catalog of genetic loci associated with kidney function from analyses of a million individuals

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Author Correction: The copy number variation landscape of congenital anomalies of the kidney and urinary tract

In the version of this article initially published, affiliation 38 incorrectly read "ICNU-Nephrology and Urology Department, Barcelona, Spain"; "Renal Division, Hospital Clinic, IDIBAPS, University of Barcelona, Barcelona, Spain" is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

The copy number variation landscape of congenital anomalies of the kidney and urinary tract

Congenital anomalies of the kidney and urinary tract (CAKUT) are a major cause of pediatric kidney failure. We performed a genome-wide analysis of copy number variants (CNVs) in 2,824 cases and 21,498 controls. Affected individuals carried a significant burden of rare exonic (i.e. affecting coding regions) CNVs and were enriched for known genomic disorders (GD). Kidney anomaly (KA) cases were most enriched for exonic CNVs, encompassing GD-CNVs and novel deletions; obstructive uropathy (OU) had a lower CNV burden and an intermediate prevalence of GD-CNVs; vesicoureteral reflux (VUR) had the fewest GD-CNVs but was enriched for novel exonic CNVs, particularly duplications. Six loci (1q21, 4p16.1-p16.3, 16p11.2, 16p13.11, 17q12, and 22q11.2) accounted for 65% of patients with GD-CNVs. Deletions at 17q12, 4p16.1-p16.3, and 22q11.2 were specific for KA; the 16p11.2 locus showed extensive pleiotropy. Using a multidisciplinary approach, we identified TBX6 as a driver for the CAKUT subphenotypes in the 16p11.2 microdeletion syndrome.

Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis Risk

Multiple sclerosis is a complex neurological disease, with ∼20% of risk heritability attributable to common genetic variants, including >230 identified by genome-wide association studies. Multiple strands of evidence suggest that much of the remaining heritability is also due to additive effects of common variants rather than epistasis between these variants or mutations exclusive to individual families. Here, we show in 68,379 cases and controls that up to 5% of this heritability is explained by low-frequency variation in gene coding sequence. We identify four novel genes driving MS risk independently of common-variant signals, highlighting key pathogenic roles for regulatory T cell homeostasis and regulation, IFNγ biology, and NFκB signaling. As low-frequency variants do not show substantial linkage disequilibrium with other variants, and as coding variants are more interpretable and experimentally tractable than non-coding variation, our discoveries constitute a rich resource for dissecting the pathobiology of MS.

Claudin-14 Gene Polymorphisms and Urine Calcium Excretion

BACKGROUND AND OBJECTIVES

Claudin-16 and -19 are proteins forming pores for the paracellular reabsorption of divalent cations in the ascending limb of Henle loop; conversely, claudin-14 decreases ion permeability of these pores. Single-nucleotide polymorphisms in gene coding for claudin-14 were associated with kidney stones and calcium excretion. This study aimed to explore the association of claudin-14, claudin-16, and claudin-19 single-nucleotide polymorphisms with calcium excretion.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

We performed a retrospective observational study of 393 patients with hypertension who were naïve to antihypertensive drugs, in whom we measured 24-hour urine calcium excretion; history of kidney stones was ascertained by interview; 370 of these patients underwent an intravenous 0.9% sodium chloride infusion (2 L in 2 hours) to evaluate the response of calcium excretion in three different 2-hour urine samples collected before, during, and after saline infusion. Genotypes of claudin-14, claudin-16, and claudin-19 were obtained from data of a previous genome-wide association study in the same patients.

RESULTS

Thirty-one single-nucleotide polymorphisms of the 3' region of the claudin-14 gene were significantly associated with 24-hour calcium excretion and calcium excretion after saline infusion. The most significant associated single-nucleotide polymorphism was rs219755 (24-hour calcium excretion in GG, 225±124 mg/24 hours; 24-hour calcium excretion in GA, 194±100 mg/24 hours; 24-hour calcium excretion in AA, 124±73 mg/24 hours; P<0.001; calcium excretion during saline infusion in GG, 30±21 mg/2 hours; calcium excretion during saline infusion in GA, 29±18 mg/2 hours; calcium excretion during saline infusion in AA, 17±11 mg/2 hours; P=0.03). No significant associations were found among claudin-16 and claudin-19 single-nucleotide polymorphisms and calcium excretion and between claudin-14, claudin-16, and claudin-19 single-nucleotide polymorphisms and stones. Bioinformatic analysis showed that one single-nucleotide polymorphism at claudin-14 among those associated with calcium excretion may potentially influence splicing of transcript.

CONCLUSIONS

Claudin-14 genotype at the 3' region is associated with calcium excretion in 24-hour urine and after the calciuretic stimulus of saline infusion.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

The risk of nephrolithiasis is causally related to inactive matrix Gla protein, a marker of vitamin K status: a Mendelian randomization study in a Flemish population

Background

Vitamin K (VK)-dependent γ-glutamate carboxylation and serine phosphorylation activate matrix Gla protein (MGP) to a potent locally acting inhibitor of calcification. Nephrolithiasis represents a process of unwanted calcification associated with substantial mortality and high recurrence rates. We hypothesized that the risk of nephrolithiasis increases with VK shortage, as exemplified by higher plasma levels of desphospho-uncarboxylated MGP (dp-ucMGP).

Methods

In 1748 randomly recruited Flemish individuals (51.1% women; mean age 46.8 years), we determined dp-ucMGP and the prevalence of nephrolithiasis at baseline (April 1996-February 2015) and its incidence during follow-up until March 2016. We estimated the multivariable-adjusted relative risk associated with the doubling of dp-ucMGP, using logistic or Cox regression. We did a Mendelian randomization analysis using four MGP genotypes as instrumental variables.

Results

With adjustments applied for sex, age and 24-h urinary volume and calcium excretion, the odds of having prevalent nephrolithiasis [n = 144 (8.2%)] associated with dp-ucMGP was 1.31 [95% confidence interval (CI) 1.04-1.64; P = 0.022]. dp-ucMGP levels were associated (P ≤ 0.001) with MGP variants rs2098435, rs4236 and rs2430692. In the Mendelian analysis, the causal odds ratio was 3.82 (95% CI 1.15-12.7; P = 0.029). The incidence of nephrolithiasis over 12.0 years (median) was 37 cases (0.2%). With similar adjustments as before, the hazard ratio in relation to dp-ucMGP was 2.48 (95% CI 1.71-3.61; P < 0.001). Additional adjustment for a nephrolithiasis propensity score produced consistent results.

Conclusion

Higher levels of inactive dp-ucMGP may be causally associated with the risk of nephrolithiasis. Whether or not VK deficiency plays a role in these observations remains to be firmly established.

A novel network analysis approach reveals DNA damage, oxidative stress and calcium/cAMP homeostasis-associated biomarkers in frontotemporal dementia

Frontotemporal Dementia (FTD) is the form of neurodegenerative dementia with the highest prevalence after Alzheimer’s disease, equally distributed in men and women. It includes several variants, generally characterized by behavioural instability and language impairments. Although few mendelian genes (MAPT, GRN, and C9orf72) have been associated to the FTD phenotype, in most cases there is only evidence of multiple risk loci with relatively small effect size. To date, there are no comprehensive studies describing FTD at molecular level, highlighting possible genetic interactions and signalling pathways at the origin FTD-associated neurodegeneration. In this study, we designed a broad FTD genetic interaction map of the Italian population, through a novel network-based approach modelled on the concepts of disease-relevance and interaction perturbation, combining Steiner tree search and Structural Equation Model (SEM) analysis. Our results show a strong connection between Calcium/cAMP metabolism, oxidative stress-induced Serine/Threonine kinases activation, and postsynaptic membrane potentiation, suggesting a possible combination of neuronal damage and loss of neuroprotection, leading to cell death. In our model, Calcium/cAMP homeostasis and energetic metabolism impairments are primary causes of loss of neuroprotection and neural cell damage, respectively. Secondly, the altered postsynaptic membrane potentiation, due to the activation of stress-induced Serine/Threonine kinases, leads to neurodegeneration. Our study investigates the molecular underpinnings of these processes, evidencing key genes and gene interactions that may account for a significant fraction of unexplained FTD aetiology. We emphasized the key molecular actors in these processes, proposing them as novel FTD biomarkers that could be crucial for further epidemiological and molecular studies.

Genetic loci associated with heart rate variability and their effects on cardiac disease risk

Reduced cardiac vagal control reflected in low heart rate variability (HRV) is associated with greater risks for cardiac morbidity and mortality. In two-stage meta-analyses of genome-wide association studies for three HRV traits in up to 53,174 individuals of European ancestry, we detect 17 genome-wide significant SNPs in eight loci. HRV SNPs tag non-synonymous SNPs (in NDUFA11 and KIAA1755), expression quantitative trait loci (eQTLs) (influencing GNG11, RGS6 and NEO1), or are located in genes preferentially expressed in the sinoatrial node (GNG11, RGS6 and HCN4). Genetic risk scores account for 0.9 to 2.6% of the HRV variance. Significant genetic correlation is found for HRV with heart rate (-0.74<r_g<-0.55) and blood pressure (-0.35<r_g<-0.20). These findings provide clinically relevant biological insight into heritable variation in vagal heart rhythm regulation, with a key role for genetic variants (GNG11, RGS6) that influence G-protein heterotrimer action in GIRK-channel induced pacemaker membrane hyperpolarization.

Genetic association study of exfoliation syndrome identifies a protective rare variant at LOXL1 and five new susceptibility loci

Exfoliation syndrome (XFS) is the most common known risk factor for secondary glaucoma and a major cause of blindness worldwide. Variants in two genes, LOXL1 and CACNA1A, have previously been associated with XFS. To further elucidate the genetic basis of XFS, we collected a global sample of XFS cases to refine the association at LOXL1, which previously showed inconsistent results across populations, and to identify new variants associated with XFS. We identified a rare protective allele at LOXL1 (p.Phe407, odds ratio (OR) = 25, P = 2.9 × 10^-14) through deep resequencing of XFS cases and controls from nine countries. A genome-wide association study (GWAS) of XFS cases and controls from 24 countries followed by replication in 18 countries identified seven genome-wide significant loci (P < 5 × 10^-8). We identified association signals at 13q12 (POMP), 11q23.3 (TMEM136), 6p21 (AGPAT1), 3p24 (RBMS3) and 5q23 (near SEMA6A). These findings provide biological insights into the pathology of XFS and highlight a potential role for naturally occurring rare LOXL1 variants in disease biology.

Pseudoexfoliation syndrome-associated genetic variants affect transcription factor binding and alternative splicing of LOXL1

Although lysyl oxidase-like 1 (LOXL1) is known as the principal genetic risk factor for pseudoexfoliation (PEX) syndrome, a major cause of glaucoma and cardiovascular complications, no functional variants have been identified to date. Here, we conduct a genome-wide association scan on 771 German PEX patients and 1,350 controls, followed by independent testing of associated variants in Italian and Japanese data sets. We focus on a 3.5-kb four-component polymorphic locus positioned spanning introns 1 and 2 of LOXL1 with enhancer-like chromatin features. We find that the rs11638944:C>G transversion exerts a cis-acting effect on the expression levels of LOXL1, mediated by differential binding of the transcription factor RXRα (retinoid X receptor alpha) and by modulating alternative splicing of LOXL1, eventually leading to reduced levels of LOXL1 mRNA in cells and tissues of risk allele carriers. These findings uncover a functional mechanism by which common noncoding variants influence LOXL1 expression.

LOXL1 is a genetic risk factor for pseudoexfoliation syndrome of the eye but a causal variant has not been identified. Here, Pasutto et al., find intronic LOXL1 risk variants influence transcription factor binding and alternative splicing of LOXL1 in affected tissues reducing levels of LOXL1 mRNA.

PEAR1 is not a major susceptibility gene for cardiovascular disease in a Flemish population

Background

Platelet Endothelial Aggregation Receptor 1 (PEAR1), a membrane protein highly expressed in platelets and endothelial cells, plays a role in platelet contact-induced activation, sustained platelet aggregation and endothelial function. Previous reports implicate PEAR1 rs12041331 as a variant influencing risk in patients with coronary heart disease. We investigated whether genetic variation in PEAR1 predicts cardiovascular outcome in a white population.

Methods

In 1938 participants enrolled in the Flemish Study on Environment, Genes and Health Outcomes (51.3% women; mean age 43.6 years), we genotyped 9 tagging SNPs in PEAR1, measured baseline cardiovascular risk factors, and recorded Cardiovascular disease incidence. For SNPs, we contrasted cardiovascular disease incidence of minor-allele heterozygotes and homozygotes (variant) vs. major-allele homozygotes (reference) and for haplotypes carriers vs. non-carriers. In adjusted analyses, we accounted for family clusters and baseline covariables, including sex, age, body mass index, mean arterial pressure, the total-to-HDL cholesterol ratio, smoking and drinking, antihypertensive drug treatment, and history of cardiovascular disease and diabetes mellitus.

Results

Over a median follow-up of 15.3 years, 238 died and 181 experienced a major cardiovascular endpoint. The multivariable-adjusted hazard ratios of eight PEAR1 SNPs, including rs12566888, ranged from 0.87 to 1.07 (P ≥0.35) and from 0.78 to 1.30 (P ≥0.15), respectively. The hazard ratios of three haplotypes with frequency ≥10% ranged from 0.93 to 1.11 (P ≥0.49) for mortality and from 0.84 to 1.03 (P ≥0.29) for a cardiovascular complications. These results were not influenced by intake of antiplatelet drugs, nonsteroidal anti-inflammatory drugs, or both (P-values for interaction ≥ 0.056).

Conclusions

In a White population, we could not replicate previous reports from experimental studies or obtained in patients suggesting that PEAR1 might be a susceptibility gene for cardiovascular complications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-017-0411-x) contains supplementary material, which is available to authorized users.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genetic Drivers of Kidney Defects in the DiGeorge Syndrome

BACKGROUND

The DiGeorge syndrome, the most common of the microdeletion syndromes, affects multiple organs, including the heart, the nervous system, and the kidney. It is caused by deletions on chromosome 22q11.2; the genetic driver of the kidney defects is unknown.

METHODS

We conducted a genomewide search for structural variants in two cohorts: 2080 patients with congenital kidney and urinary tract anomalies and 22,094 controls. We performed exome and targeted resequencing in samples obtained from 586 additional patients with congenital kidney anomalies. We also carried out functional studies using zebrafish and mice.

RESULTS

We identified heterozygous deletions of 22q11.2 in 1.1% of the patients with congenital kidney anomalies and in 0.01% of population controls (odds ratio, 81.5; P=4.5×10^-14). We localized the main drivers of renal disease in the DiGeorge syndrome to a 370-kb region containing nine genes. In zebrafish embryos, an induced loss of function in snap29, aifm3, and crkl resulted in renal defects; the loss of crkl alone was sufficient to induce defects. Five of 586 patients with congenital urinary anomalies had newly identified, heterozygous protein-altering variants, including a premature termination codon, in CRKL. The inactivation of Crkl in the mouse model induced developmental defects similar to those observed in patients with congenital urinary anomalies.

CONCLUSIONS

We identified a recurrent 370-kb deletion at the 22q11.2 locus as a driver of kidney defects in the DiGeorge syndrome and in sporadic congenital kidney and urinary tract anomalies. Of the nine genes at this locus, SNAP29, AIFM3, and CRKL appear to be critical to the phenotype, with haploinsufficiency of CRKL emerging as the main genetic driver. (Funded by the National Institutes of Health and others.).

SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium (n_Stage1: 111,666; n_Stage2: 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea (PPM1J, EDEM3, ACP1, SPEG, EYA4, CYP1A1, and ATXN2L; P_Stage1<3.7×10^-7), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 (P=5.4×10^-8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2-knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Genome-Wide and Gene-Based Meta-Analyses Identify Novel Loci Influencing Blood Pressure Response to Hydrochlorothiazide

This study aimed to identify novel loci influencing the antihypertensive response to hydrochlorothiazide monotherapy. A genome-wide meta-analysis of blood pressure (BP) response to hydrochlorothiazide was performed in 1739 white hypertensives from 6 clinical trials within the International Consortium for Antihypertensive Pharmacogenomics Studies, making it the largest study to date of its kind. No signals reached genome-wide significance (P<5×10^-8), and the suggestive regions (P<10^-5) were cross-validated in 2 black cohorts treated with hydrochlorothiazide. In addition, a gene-based analysis was performed on candidate genes with previous evidence of involvement in diuretic response, in BP regulation, or in hypertension susceptibility. Using the genome-wide meta-analysis approach, with validation in blacks, we identified 2 suggestive regulatory regions linked to gap junction protein α1 gene (GJA1) and forkhead box A1 gene (FOXA1), relevant for cardiovascular and kidney function. With the gene-based approach, we identified hydroxy-delta-5-steroid dehydrogenase, 3 β- and steroid δ-isomerase 1 gene (HSD3B1) as significantly associated with BP response (P<2.28×10^-4 ). HSD3B1 encodes the 3β-hydroxysteroid dehydrogenase enzyme and plays a crucial role in the biosynthesis of aldosterone and endogenous ouabain. By amassing all of the available pharmacogenomic studies of BP response to hydrochlorothiazide, and using 2 different analytic approaches, we identified 3 novel loci influencing BP response to hydrochlorothiazide. The gene-based analysis, never before applied to pharmacogenomics of antihypertensive drugs to our knowledge, provided a powerful strategy to identify a locus of interest, which was not identified in the genome-wide meta-analysis because of high allelic heterogeneity. These data pave the way for future investigations on new pathways and drug targets to enhance the current understanding of personalized antihypertensive treatment.

PHACTR1 Is a Genetic Susceptibility Locus for Fibromuscular Dysplasia Supporting Its Complex Genetic Pattern of Inheritance

Fibromuscular dysplasia (FMD) is a nonatherosclerotic vascular disease leading to stenosis, dissection and aneurysm affecting mainly the renal and cerebrovascular arteries. FMD is often an underdiagnosed cause of hypertension and stroke, has higher prevalence in females (~80%) but its pathophysiology is unclear. We analyzed ~26K common variants (MAF>0.05) generated by exome-chip arrays in 249 FMD patients and 689 controls. We replicated 13 loci (P<10-4) in 402 cases and 2,537 controls and confirmed an association between FMD and a variant in the phosphatase and actin regulator 1 gene (PHACTR1). Three additional case control cohorts including 512 cases and 669 replicated this result and overall reached the genomic level of significance (OR = 1.39, P = 7.4×10-10, 1,154 cases and 3,895 controls). The top variant, rs9349379, is intronic to PHACTR1, a risk locus for coronary artery disease, migraine, and cervical artery dissection. The analyses of geometrical parameters of carotids from ~2,500 healthy volunteers indicate higher intima media thickness (P = 1.97×10-4) and wall to lumen ratio (P = 0.002) in rs9349379-A carriers, suggesting indices of carotid hypertrophy previously described in carotids of FMD patients. Immunohistochemistry detected PHACTR1 in endothelium and smooth muscle cells of FMD and normal human carotids. The expression of PHACTR1 by genotypes in primary human fibroblasts showed higher expression in rs9349379-A carriers (N = 86, P = 0.003). Phactr1 knockdown in zebrafish resulted in dilated vessels indicating subtle impaired vascular development. We report the first susceptibility locus for FMD and provide evidence for a complex genetic pattern of inheritance and indices of shared pathophysiology between FMD and other cardiovascular and neurovascular diseases.

A candidate gene study identifies a haplotype of CD2 as novel susceptibility factor for systemic sclerosis

OBJECTIVES

Systemic sclerosis (SSc) is a rare autoimmune disease (AID) with a complex genetic etiology. Evidence for a shared pathogenesis across AIDs is given by the well-known pleiotropism of autoimmune genes. Recently, several unbiased approaches have identified an association between polymorphisms of the CD2 gene, and rheumatoid arthritis (RA) susceptibility. The objective of this study was to investigate whether CD2 polymorphisms are associated with SSc.

METHODS

Two SNPs of CD2, rs624988 and rs798036, were genotyped in a total of 1,786 SSc patients and 2,360 healthy individuals from two European populations (France and Italy). Meta-analyses were performed to assess whether an association exists between CD2 polymorphisms or haplotypes and SSc or its main subtypes.

RESULTS

The combined analyses revealed an association between the rs624988 A allele and SSc susceptibility: padj=0.023, OR=1.14 (95%CI 1.04-1.25). Single marker analysis did not reveal any association between rs798036 and SSc. Haplotype analysis identified that the A-T haplotype, previously described in RA, was associated with higher susceptibility for SSc (padj=0.029, OR=1.14, 95%CI 1.04-1.25) and with the positive anti-centromere antibody sub-group of SSc patients (padj=0.009, OR=1.19 95%CI 1.07-1.32). Genotype-mRNA expression correlations revealed that the CD2 risk haplotype was associated with decreased CD2 mRNA expression in SSc patients.

CONCLUSIONS

Our study establishes CD2 as a new susceptibility factor for SSc, in a European Caucasian population, confirming the sharing of autoimmune risk factors by SSc and RA.

The burden of multiple sclerosis variants in continental Italians and Sardinians

BACKGROUND

Recent studies identified > 100 non-HLA (human leukocyte antigen) multiple sclerosis (MS) susceptibility variants in Northern European populations, but their role in Southern Europeans is largely unexplored.

OBJECTIVE

We aimed to investigate the cumulative impact of those variants in two Mediterranean populations: Continental Italians and Sardinians.

METHODS

We calculated four weighted Genetic Risk Scores (wGRS), using up to 102 non-HLA MS risk variants and 5 HLA MS susceptibility markers in 1691 patients and 2194 controls from continental Italy; and 2861 patients and 3034 controls from Sardinia. We then assessed the differences between populations using Nagelkerke's R(2) and the area under the Receiver Operating Characteristic (ROC) curves.

RESULTS

As expected, the genetic burden (mean wGRS value) was significantly higher in MS patients than in controls, in both populations. Of note, the burden was significantly higher in Sardinians. Conversely, the proportion of variability explained and the predictive power were significantly higher in continental Italians. Notably, within the Sardinian patients, we also observed a significantly higher burden of non-HLA variants in individuals who do not carry HLA risk alleles.

CONCLUSIONS

The observed differences in MS genetic burden between the two Mediterranean populations highlight the need for more genetic studies in South Europeans, to further expand the knowledge of MS genetics.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Interaction between polyphenols intake and PON1 gene variants on markers of cardiovascular disease: a nutrigenetic observational study

Background

Paraoxonase 1 (PON1) gene polymorphisms and polyphenols intake have been reported independently associated to lipid profile and susceptibility to atherosclerosis and cardiovascular disease. However, the interaction between these factors remains to be investigated. We performed an observational nutrigenetic study to examine whether the interaction between polyphenols and anthocyanins intake and PON1 genetic variants can modulate biomarkers of cardiovascular health in an Italian healthy population.

Methods

We recruited 443 healthy volunteers who participated in the EC funded ATHENA project (AnThocyanin and polyphenols bioactive for Health Enhancement through Nutritional Advancement). Data collection included detailed demographic, clinical, dietary, lifestyle, biochemical and genetic data. Polyphenols and anthocyanins intake was measured by 24 h dietary recall repeated three times a year in order to get seasonal variations. We tested the interaction between 18 independent tagging SNPs in PON1 gene and polyphenols intake on HDL, LDL, cholesterol, triglycerides and atherogenic index of plasma.

Results

Without considering the genetic background, we could not observe significant differences in the lipid profile between high and low polyphenols and anthocyanins intake. Using a nutrigenetic approach, we identified protective genotypes in four independent polymorphisms that, at Bonferroni level (p ≤ 0.0028), present a significant association with increased HDL level under high polyphenols and anthocyanins intake, compared to risk genotypes (rs854549, Beta = 4.7 per C allele; rs854552, Beta = 5.6 per C allele; rs854571, Beta = 3.92 per T allele; rs854572, Beta = 3.94 per C allele).

Conclusions

We highlight the protective role of genetic variants in PON1 towards cardiovascular risk under high polyphenols and anthocyanins consumption. PON1 variants could represent novel biomarkers to stratify individuals who might benefit from targeted dietary recommendation for health promotion and strategies of preventive medicine.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-016-0941-6) contains supplementary material, which is available to authorized users.

Genome-wide association study identifies 74 loci associated with educational attainment

Educational attainment is strongly influenced by social and other environmental factors, but genetic factors are estimated to account for at least 20% of the variation across individuals. Here we report the results of a genome-wide association study (GWAS) for educational attainment that extends our earlier discovery sample of 101,069 individuals to 293,723 individuals, and a replication study in an independent sample of 111,349 individuals from the UK Biobank. We identify 74 genome-wide significant loci associated with the number of years of schooling completed. Single-nucleotide polymorphisms associated with educational attainment are disproportionately found in genomic regions regulating gene expression in the fetal brain. Candidate genes are preferentially expressed in neural tissue, especially during the prenatal period, and enriched for biological pathways involved in neural development. Our findings demonstrate that, even for a behavioural phenotype that is mostly environmentally determined, a well-powered GWAS identifies replicable associated genetic variants that suggest biologically relevant pathways. Because educational attainment is measured in large numbers of individuals, it will continue to be useful as a proxy phenotype in efforts to characterize the genetic influences of related phenotypes, including cognition and neuropsychiatric diseases.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

A genome-wide screening and SNPs-to-genes approach to identify novel genetic risk factors associated with frontotemporal dementia

Frontotemporal dementia (FTD) is the second most prevalent form of early onset dementia after Alzheimer's disease (AD). We performed a case-control association study in an Italian FTD cohort (n = 530) followed by the novel single nucleotide polymorphisms (SNPs)-to-genes approach and functional annotation analysis. We identified 2 novel potential loci for FTD. Suggestive SNPs reached p-values ∼10⁻⁷ and odds ratio > 2.5 (2p16.3) and 1.5 (17q25.3). Suggestive alleles at 17q25.3 identified a disease-associated haplotype causing decreased expression of –cis genes such as RFNG and AATK involved in neuronal genesis and differentiation and axon outgrowth, respectively. We replicated this locus through the SNPs-to-genes approach. Our functional annotation analysis indicated significant enrichment for functions of the brain (neuronal genesis, differentiation, and maturation), the synapse (neurotransmission and synapse plasticity), and elements of the immune system, the latter supporting our recent international FTD–genome-wide association study. This is the largest genome-wide study in Italian FTD to date. Although our results are not conclusive, we set the basis for future replication studies and identification of susceptible molecular mechanisms involved in FTD pathogenesis.