ProbABEL package for genome-wide association analysis of imputed data

Common Polymorphisms in IFI16 and AIM2 Genes Are Associated With Periodontal Disease

BACKGROUND

The single nucleotide polymorphism (SNP) context of a previously identified periodontitis-associated locus is investigated, and its association with microbial, biologic, and periodontal disease clinical parameters is examined.

METHODS

A 200-kb spanning region of 1q12 previously highlighted in a genome-wide association scan among 4,766 European American individuals (SNP rs1633266) was annotated. Two haplotype blocks were selected. Association of these polymorphisms with data on microbial plaque composition, gingival crevicular fluid (GCF)-interleukin (IL)-1β levels, and clinical parameters of periodontal disease were examined. Descriptive analysis of IFI16 and AIM2 protein expression in gingival tissues from healthy individuals (n = 2) and individuals with chronic periodontitis (n = 2) was done via immunohistochemistry.

RESULTS

The highlighted locus is a 100-kb region containing the interferon γ-inducible protein 16 (IFI16) and absent in melanoma 2 (AIM2) genes. Two haplotype blocks, rs6940 and rs1057028, were significantly associated with increased extent bleeding on probing and levels of microorganisms Porphyromonas gingivalis, Tannerella forsythia, and Campylobacter rectus (P ≤0.05). Haplotype block rs1057028 was also significantly associated with pathogens Fusobacterium nucleatum and Aggregatibacter actinomycetemcomitans, increased GCF-IL-1β levels, and extent of probing depth ≥4 mm (P ≤0.05). Prevalence of severe periodontitis (biofilm-gingival interface P3 classification) was positively associated with haplotype block rs1057028. Similar trends were observed for haplotype block rs1057028. IFI16 and AIM2 protein expression was observed in multiple cell types of gingival tissues, including inflammatory cells.

CONCLUSION

This study found IFI16 and AIM2 SNPs associated with higher levels of periodontal microorganisms and an increased percentage of periodontal disease clinical parameters, suggesting the need for functional studies and additional fine-mapping of variants in the 1q12-locus.

Robustness of meta-analyses in finding gene × environment interactions

Meta-analyses that synthesize statistical evidence across studies have become important analytical tools for genetic studies. Inspired by the success of genome-wide association studies of the genetic main effect, researchers are searching for gene × environment interactions. Confounders are routinely included in the genome-wide gene × environment interaction analysis as covariates; however, this does not control for any confounding effects on the results if covariate × environment interactions are present. We carried out simulation studies to evaluate the robustness to the covariate × environment confounder for meta-regression and joint meta-analysis, which are two commonly used meta-analysis methods for testing the gene × environment interaction or the genetic main effect and interaction jointly. Here we show that meta-regression is robust to the covariate × environment confounder while joint meta-analysis is subject to the confounding effect with inflated type I error rates. Given vast sample sizes employed in genome-wide gene × environment interaction studies, non-significant covariate × environment interactions at the study level could substantially elevate the type I error rate at the consortium level. When covariate × environment confounders are present, type I errors can be controlled in joint meta-analysis by including the covariate × environment terms in the analysis at the study level. Alternatively, meta-regression can be applied, which is robust to potential covariate × environment confounders.

Screening for interaction effects in gene expression data

Expression quantitative trait (eQTL) studies are a powerful tool for identifying genetic variants that affect levels of messenger RNA. Since gene expression is controlled by a complex network of gene-regulating factors, one way to identify these factors is to search for interaction effects between genetic variants and mRNA levels of transcription factors (TFs) and their respective target genes. However, identification of interaction effects in gene expression data pose a variety of methodological challenges, and it has become clear that such analyses should be conducted and interpreted with caution. Investigating the validity and interpretability of several interaction tests when screening for eQTL SNPs whose effect on the target gene expression is modified by the expression level of a transcription factor, we characterized two important methodological issues. First, we stress the scale-dependency of interaction effects and highlight that commonly applied transformation of gene expression data can induce or remove interactions, making interpretation of results more challenging. We then demonstrate that, in the setting of moderate to strong interaction effects on the order of what may be reasonably expected for eQTL studies, standard interaction screening can be biased due to heteroscedasticity induced by true interactions. Using simulation and real data analysis, we outline a set of reasonable minimum conditions and sample size requirements for reliable detection of variant-by-environment and variant-by-TF interactions using the heteroscedasticity consistent covariance-based approach.

A Review of the Genetics of Hypertension with a Focus on Gene-Environment Interactions

PURPOSE OF REVIEW

Here, we discuss the interpretation and modeling of gene-environment interactions in hypertension-related phenotypes, with a focus on the necessary assumptions and possible challenges.

RECENT FINDINGS

Recently, small cohort studies have discovered several novel genetic variants associated with hypertension-related phenotypes through modeling gene-environment interactions. Several consortia-based meta-analytic efforts have uncovered many novel genetic variants in hypertension without modeling interaction terms, giving promise to future meta-analytic efforts that incorporate gene-environment interactions. Heritability studies and genome-wide association studies have established that hypertension, a prevalent cardiovascular disease, has a genetic component that may be modulated by the environment (such as lifestyle factors). This review includes a discussion of known genetic associations for hypertension/blood pressure, including those resulting from the incorporation of gene-environmental interaction modeling.

Genome-Wide Association Studies of Multiple Keratinocyte Cancers

There is strong evidence for a role of environmental risk factors involved in susceptibility to develop multiple keratinocyte cancers (mKCs), but whether genes are also involved in mKCs susceptibility has not been thoroughly investigated. We investigated whether single nucleotide polymorphisms (SNPs) are associated with susceptibility for mKCs. A genome-wide association study (GWAS) of 1,666 cases with mKCs and 1,950 cases with single KC (sKCs; controls) from Harvard cohorts (the Nurses' Health Study [NHS], NHS II, and the Health Professionals Follow-Up Study) and the Framingham Heart Study was carried-out using over 8 million SNPs (stage-1). We sought to replicate the most significant statistical associations (p-value≤ 5.5x10-6) in an independent cohort of 574 mKCs and 872 sKCs from the Rotterdam Study. In the discovery stage, 40 SNPs with suggestive associations (p-value ≤5.5x10-6) were identified, with eight independent SNPs tagging all 40 SNPs. The most significant SNP was located at chromosome 9 (rs7468390; p-value = 3.92x10-7). In stage-2, none of these SNPs replicated and only two of them were associated with mKCs in the same direction in the combined meta-analysis. We tested the associations for 19 previously reported basal cell carcinoma-related SNPs (candidate gene association analysis), and found that rs1805007 (MC1R locus) was significantly associated with risk of mKCs (p-value = 2.80x10-4). Although the suggestive SNPs with susceptibility for mKCs were not replicated, we found that previously identified BCC variants may also be associated with mKC, which the most significant association (rs1805007) located at the MC1R gene.

A genome-wide analysis of gene-caffeine consumption interaction on basal cell carcinoma

Animal models have suggested that oral or topical administration of caffeine could inhibit ultraviolet-induced carcinogenesis via the ataxia telangiectasia and rad3 (ATR)-related apoptosis. Previous epidemiological studies have demonstrated that increased caffeine consumption is associated with reduced risk of basal cell carcinoma (BCC). To identify common genetic markers that may modify this association, we tested gene-caffeine intake interaction on BCC risk in a genome-wide analysis. We included 3383 BCC cases and 8528 controls of European ancestry from the Nurses' Health Study and Health Professionals Follow-up Study. Single nucleotide polymorphism (SNP) rs142310826 near the NEIL3 gene showed a genome-wide significant interaction with caffeine consumption (P = 1.78 × 10-8 for interaction) on BCC risk. There was no gender difference for this interaction (P = 0.64 for heterogeneity). NEIL3, a gene belonging to the base excision DNA repair pathway, encodes a DNA glycosylase that recognizes and removes lesions produced by oxidative stress. In addition, we identified several loci with P value for interaction <5 × 10-7 in gender-specific analyses (P for heterogeneity between genders < 0.001) including those mapping to the genes LRRTM4, ATF3 and DCLRE1C in women and POTEA in men. Finally, we tested the associations between caffeine consumption-related SNPs reported by previous genome-wide association studies and risk of BCC, both individually and jointly, but found no significant association. In sum, we identified a DNA repair gene that could be involved in caffeine-mediated skin tumor inhibition. Further studies are warranted to confirm these findings.

A comprehensive genome-wide analysis of melanoma Breslow thickness identifies interaction between CDC42 and SCIN genetic variants

Breslow thickness (BT) is a major prognostic factor of cutaneous melanoma (CM), the most fatal skin cancer. The genetic component of BT has only been explored by candidate gene studies with inconsistent results. Our objective was to uncover the genetic factors underlying BT using an hypothesis-free genome-wide approach. Our analysis strategy integrated a genome-wide association study (GWAS) of single nucleotide polymorphisms (SNPs) for BT followed by pathway analysis of GWAS outcomes using the gene-set enrichment analysis (GSEA) method and epistasis analysis within BT-associated pathways. This strategy was applied to two large CM datasets with Hapmap3-imputed SNP data: the French MELARISK study for discovery (966 cases) and the MD Anderson Cancer Center study (1,546 cases) for replication. While no marginal effect of individual SNPs was revealed through GWAS, three pathways, defined by gene ontology (GO) categories were significantly enriched in genes associated with BT (false discovery rate ≤5% in both studies): hormone activity, cytokine activity and myeloid cell differentiation. Epistasis analysis, within each significant GO, identified a statistically significant interaction between CDC42 and SCIN SNPs (pmeta-int =2.2 × 10(-6) , which met the overall multiple-testing corrected threshold of 2.5 × 10(-6) ). These two SNPs (and proxies) are strongly associated with CDC42 and SCIN gene expression levels and map to regulatory elements in skin cells. This interaction has important biological relevance since CDC42 and SCIN proteins have opposite effects in actin cytoskeleton organization and dynamics, a key mechanism underlying melanoma cell migration and invasion.

Ascertainment bias from imputation methods evaluation in wheat

BACKGROUND

Whole-genome genotyping techniques like Genotyping-by-sequencing (GBS) are being used for genetic studies such as Genome-Wide Association (GWAS) and Genomewide Selection (GS), where different strategies for imputation have been developed. Nevertheless, imputation error may lead to poor performance (i.e. smaller power or higher false positive rate) when complete data is not required as it is for GWAS, and each marker is taken at a time. The aim of this study was to compare the performance of GWAS analysis for Quantitative Trait Loci (QTL) of major and minor effect using different imputation methods when no reference panel is available in a wheat GBS panel.

RESULTS

In this study, we compared the power and false positive rate of dissecting quantitative traits for imputed and not-imputed marker score matrices in: (1) a complete molecular marker barley panel array, and (2) a GBS wheat panel with missing data. We found that there is an ascertainment bias in imputation method comparisons. Simulating over a complete matrix and creating missing data at random proved that imputation methods have a poorer performance. Furthermore, we found that when QTL were simulated with imputed data, the imputation methods performed better than the not-imputed ones. On the other hand, when QTL were simulated with not-imputed data, the not-imputed method and one of the imputation methods performed better for dissecting quantitative traits. Moreover, larger differences between imputation methods were detected for QTL of major effect than QTL of minor effect. We also compared the different marker score matrices for GWAS analysis in a real wheat phenotype dataset, and we found minimal differences indicating that imputation did not improve the GWAS performance when a reference panel was not available.

CONCLUSIONS

Poorer performance was found in GWAS analysis when an imputed marker score matrix was used, no reference panel is available, in a wheat GBS panel.

An Adaptive Fisher's Combination Method for Joint Analysis of Multiple Phenotypes in Association Studies

Currently, the analyses of most genome-wide association studies (GWAS) have been performed on a single phenotype. There is increasing evidence showing that pleiotropy is a widespread phenomenon in complex diseases. Therefore, using only one single phenotype may lose statistical power to identify the underlying genetic mechanism. There is an increasing need to develop and apply powerful statistical tests to detect association between multiple phenotypes and a genetic variant. In this paper, we develop an Adaptive Fisher’s Combination (AFC) method for joint analysis of multiple phenotypes in association studies. The AFC method combines p-values obtained in standard univariate GWAS by using the optimal number of p-values which is determined by the data. We perform extensive simulations to evaluate the performance of the AFC method and compare the power of our method with the powers of TATES, Tippett’s method, Fisher’s combination test, MANOVA, MultiPhen, and SUMSCORE. Our simulation studies show that the proposed method has correct type I error rates and is either the most powerful test or comparable with the most powerful test. Finally, we illustrate our proposed methodology by analyzing whole-genome genotyping data from a lung function study.

Genetic variants in RBFOX3 are associated with sleep latency

Time to fall asleep (sleep latency) is a major determinant of sleep quality. Chronic, long sleep latency is a major characteristic of sleep-onset insomnia and/or delayed sleep phase syndrome. In this study we aimed to discover common polymorphisms that contribute to the genetics of sleep latency. We performed a meta-analysis of genome-wide association studies (GWAS) including 2 572 737 single nucleotide polymorphisms (SNPs) established in seven European cohorts including 4242 individuals. We found a cluster of three highly correlated variants (rs9900428, rs9907432 and rs7211029) in the RNA-binding protein fox-1 homolog 3 gene (RBFOX3) associated with sleep latency (P-values=5.77 × 10(-08), 6.59 × 10(-)(08) and 9.17 × 10(-)(08)). These SNPs were replicated in up to 12 independent populations including 30 377 individuals (P-values=1.5 × 10(-)(02), 7.0 × 10(-)(03) and 2.5 × 10(-)(03); combined meta-analysis P-values=5.5 × 10(-07), 5.4 × 10(-07) and 1.0 × 10(-07)). A functional prediction of RBFOX3 based on co-expression with other genes shows that this gene is predominantly expressed in brain (P-value=1.4 × 10(-316)) and the central nervous system (P-value=7.5 × 10(-)(321)). The predicted function of RBFOX3 based on co-expression analysis with other genes shows that this gene is significantly involved in the release cycle of neurotransmitters including gamma-aminobutyric acid and various monoamines (P-values<2.9 × 10(-11)) that are crucial in triggering the onset of sleep. To conclude, in this first large-scale GWAS of sleep latency we report a novel association of variants in RBFOX3 gene. Further, a functional prediction of RBFOX3 supports the involvement of RBFOX3 with sleep latency.

Testing the link between genome size and growth rate in maize

Little is known about the factors driving within species Genome Size (GS) variation. GS may be shaped indirectly by natural selection on development and adaptative traits. Because GS variation is particularly pronounced in maize, we have sampled 83 maize inbred lines from three well described genetic groups adapted to contrasted climate conditions: inbreds of tropical origin, Flint inbreds grown in temperate climates, and Dent inbreds distributed in the Corn Belt. As a proxy for growth rate, we measured the Leaf Elongation Rate maximum during nighttime (LER_max) as well as GS in all inbred lines. In addition we combined available and new nucleotide polymorphism data at 29,090 sites to characterize the genetic structure of our panel. We found significant variation for both LER_max and GS among groups defined by our genetic structuring. Tropicals displayed larger GS than Flints while Dents exhibited intermediate values. LER_max followed the opposite trend with greater growth rate in Flints than in Tropicals. In other words, LER_max and GS exhibited a significantly negative correlation (r = − 0.27). However, this correlation was driven by among-group variation rather than within-group variation—it was no longer significant after controlling for structure and kinship among inbreds. Our results indicate that selection on GS may have accompanied ancient maize diffusion from its center of origin, with large DNA content excluded from temperate areas. Whether GS has been targeted by more intense selection during modern breeding within groups remains an open question.

Identification of genetic modifiers of age-at-onset for familial Parkinson's disease

Parkinson's disease (PD) is the most common cause of neurodegenerative movement disorder and the second most common cause of dementia. Genes are thought to have a stronger effect on age-at-onset of PD than on risk, yet there has been a phenomenal success in identifying risk loci but not age-at-onset modifiers. We conducted a genome-wide study for age-at-onset. We analysed familial and non-familial PD separately, per prior evidence for strong genetic effect on age-at-onset in familial PD. GWAS was conducted in 431 unrelated PD individuals with at least one affected relative (familial PD) and 1544 non-familial PD from the NeuroGenetics Research Consortium (NGRC); an additional 737 familial PD and 2363 non-familial PD were used for replication. In familial PD, two signals were detected and replicated robustly: one mapped to LHFPL2 on 5q14.1 (PNGRC = 3E-8, PReplication = 2E-5, PNGRC + Replication = 1E-11), the second mapped to TPM1 on 15q22.2 (PNGRC = 8E-9, PReplication = 2E-4, PNGRC + Replication = 9E-11). The variants that were associated with accelerated onset had low frequencies (<0.02). The LHFPL2 variant was associated with earlier onset by 12.33 [95% CI: 6.2; 18.45] years in NGRC, 8.03 [2.95; 13.11] years in replication, and 9.79 [5.88; 13.70] years in the combined data. The TPM1 variant was associated with earlier onset by 15.30 [8.10; 22.49] years in NGRC, 9.29 [1.79; 16.79] years in replication, and 12.42 [7.23; 17.61] years in the combined data. Neither LHFPL2 nor TPM1 was associated with age-at-onset in non-familial PD. LHFPL2 (function unknown) is overexpressed in brain tumours. TPM1 encodes a highly conserved protein that regulates muscle contraction, and is a tumour-suppressor gene.

An intergenic risk locus containing an enhancer deletion in 2q35 modulates breast cancer risk by deregulating IGFBP5 expression

Breast cancer is the most diagnosed malignancy and the second leading cause of cancer mortality in females. Previous association studies have identified variants on 2q35 associated with the risk of breast cancer. To identify functional susceptibility loci for breast cancer, we interrogated the 2q35 gene desert for chromatin architecture and functional variation correlated with gene expression. We report a novel intergenic breast cancer risk locus containing an enhancer copy number variation (enCNV; deletion) located approximately 400Kb upstream to IGFBP5, which overlaps an intergenic ERα-bound enhancer that loops to the IGFBP5 promoter. The enCNV is correlated with modified ERα binding and monoallelic-repression of IGFBP5 following oestrogen treatment. We investigated the association of enCNV genotype with breast cancer in 1,182 cases and 1,362 controls, and replicate our findings in an independent set of 62,533 cases and 60,966 controls from 41 case control studies and 11 GWAS. We report a dose-dependent inverse association of 2q35 enCNV genotype (percopy OR = 0.68 95%CI 0.55-0.83, P = 0.0002; replication OR = 0.77 95% CI 0.73-0.82, P = 2.1 × 10-19) and identify 13 additional linked variants (r2 > 0.8) in the 20Kb linkage block containing the enCNV (P = 3.2 × 10-15 - 5.6 × 10-17). These associations were independent of previously reported 2q35 variants, rs13387042/rs4442975 and rs16857609, and were stronger for ER-positive than ER-negative disease. Together, these results suggest that 2q35 breast cancer risk loci may be mediating their effect through IGFBP5.

Four Susceptibility Loci for Gallstone Disease Identified in a Meta-analysis of Genome-Wide Association Studies

BACKGROUND & AIMS

A genome-wide association study (GWAS) of 280 cases identified the hepatic cholesterol transporter ABCG8 as a locus associated with risk for gallstone disease, but findings have not been reported from any other GWAS of this phenotype. We performed a large-scale, meta-analysis of GWASs of individuals of European ancestry with available prior genotype data, to identify additional genetic risk factors for gallstone disease.

METHODS

We obtained per-allele odds ratio (OR) and standard error estimates using age- and sex-adjusted logistic regression models within each of the 10 discovery studies (8720 cases and 55,152 controls). We performed an inverse variance weighted, fixed-effects meta-analysis of study-specific estimates to identify single-nucleotide polymorphisms that were associated independently with gallstone disease. Associations were replicated in 6489 cases and 62,797 controls.

RESULTS

We observed independent associations for 2 single-nucleotide polymorphisms at the ABCG8 locus: rs11887534 (OR, 1.69; 95% confidence interval [CI], 1.54-1.86; P = 2.44 × 10(-60)) and rs4245791 (OR, 1.27; P = 1.90 × 10(-34)). We also identified and/or replicated associations for rs9843304 in TM4SF4 (OR, 1.12; 95% CI, 1.08-1.16; P = 6.09 × 10(-11)), rs2547231 in SULT2A1 (encodes a sulfoconjugation enzyme that acts on hydroxysteroids and cholesterol-derived sterol bile acids) (OR, 1.17; 95% CI, 1.12-1.21; P = 2.24 × 10(-10)), rs1260326 in glucokinase regulatory protein (OR, 1.12; 95% CI, 1.07-1.17; P = 2.55 × 10(-10)), and rs6471717 near CYP7A1 (encodes an enzyme that catalyzes conversion of cholesterol to primary bile acids) (OR, 1.11; 95% CI, 1.08-1.15; P = 8.84 × 10(-9)). Among individuals of African American and Hispanic American ancestry, rs11887534 and rs4245791 were associated positively with gallstone disease risk, whereas the association for the rs1260326 variant was inverse.

CONCLUSIONS

In this large-scale GWAS of gallstone disease, we identified 4 loci in genes that have putative functions in cholesterol metabolism and transport, and sulfonylation of bile acids or hydroxysteroids.

A thrifty variant in CREBRF strongly influences body mass index in Samoans

Samoans are a unique founder population with a high prevalence of obesity, making them well suited for identifying new genetic contributors to obesity. We conducted a genome-wide association study (GWAS) in 3,072 Samoans, discovered a variant, rs12513649, strongly associated with body mass index (BMI) (P = 5.3 × 10(-14)), and replicated the association in 2,102 additional Samoans (P = 1.2 × 10(-9)). Targeted sequencing identified a strongly associated missense variant, rs373863828 (p.Arg457Gln), in CREBRF (meta P = 1.4 × 10(-20)). Although this variant is extremely rare in other populations, it is common in Samoans (frequency of 0.259), with an effect size much larger than that of any other known common BMI risk variant (1.36-1.45 kg/m(2) per copy of the risk-associated allele). In comparison to wild-type CREBRF, the Arg457Gln variant when overexpressed selectively decreased energy use and increased fat storage in an adipocyte cell model. These data, in combination with evidence of positive selection of the allele encoding p.Arg457Gln, support a 'thrifty' variant hypothesis as a factor in human obesity.

Common Genetic Variation in CYP17A1 and Response to Abiraterone Acetate in Patients with Metastatic Castration-Resistant Prostate Cancer

Treatment with abiraterone acetate and prednisone (AA/P) prolongs survival in metastatic castration-resistant prostate cancer (mCRPC) patients. We evaluated the genetic variation in CYP17A1 as predictive of response to AA/P. A prospective collection of germline DNA prior to AA/P initiation and follow-up of a mCRPC cohort was performed. Five common single-nucleotide polymorphisms (SNPs) in CYP17A1 identified using a haplotype-based tagging algorithm were genotyped. Clinical outcomes included biochemical response and time to biochemical progression on AA/P. Logistic regression was used to assess the association between tag SNPs and biochemical response. Proportional hazards regression was used to assess the association between tag SNPs and time to biochemical progression. Odds or hazard ratio per minor allele were estimated and p-values below 0.05 were considered statistically significant. Germline DNA was successfully genotyped for four tag SNPs in 87 patients. The median age was 73 years (54-90); the median prostate-specific antigen was 66 ng/dL (0.1-99.9). A single SNP, rs2486758, was associated with lower odds of experiencing a biochemical response (Odds ratio 0.22, 95% confidence interval 0.07-0.63, p = 0.005) and a shorter time to biochemical progression (Hazard ratio 2.23, 95% confidence interval 1.39-3.56, p < 0.001). This tag SNP located in the promoter region of CYP17A1 will need further validation as a predictive biomarker for AA/P therapy.

No evidence of large genetic effects on steroid response in asthma patients

BACKGROUND

Inhaled corticosteroids (ICSs) are considered the most effective anti-inflammatory therapy for asthma control and management; however, there is substantial treatment response variability.

OBJECTIVE

We sought to identify genetic markers of ICS response by conducting the largest pharmacogenetic investigation to date in 2672 ICS-treated patients with asthma.

METHODS

Genotyping and imputation was performed in fluticasone furoate (FF) or fluticasone propionate-treated patients with asthma from 3 phase IIB and 4 phase IIIA randomized, double-blind, placebo-controlled, parallel group, multicenter studies. The primary end point analyzed was change in trough FEV₁ (ΔFEV₁) from baseline to 8 to 12 weeks of treatment.

RESULTS

More than 9.8 million common genetic variants (minor allele frequency ≥ 1%) were analyzed to test for association with ΔFEV₁. No genetic variant met the prespecified threshold for statistical significance.

CONCLUSIONS

This study provides no evidence to confirm previously reported associations between candidate genetic variants and ICS response (ΔFEV₁) in patients with asthma. In addition, no variant satisfied the criterion for genome-wide significance in our study. Common genetic variants are therefore unlikely to prove useful as predictive biomarkers of ICS response in patients with asthma.

Association of a Locus in the CAMTA1 Gene With Survival in Patients With Sporadic Amyotrophic Lateral Sclerosis

IMPORTANCE

Amyotrophic lateral sclerosis (ALS) is a devastating adult-onset neurodegenerative disorder with a poor prognosis and a median survival of 3 years. However, a significant proportion of patients survive more than 10 years from symptom onset.

OBJECTIVE

To identify gene variants influencing survival in ALS.

DESIGN, SETTING, AND PARTICIPANTS

This genome-wide association study (GWAS) analyzed survival in data sets from several European countries and the United States that were collected by the Italian Consortium for the Genetics of ALS and the International Consortium on Amyotrophic Lateral Sclerosis Genetics. The study population included 4256 patients with ALS (3125 [73.4%] deceased) with genotype data extended to 7 174 392 variants by imputation analysis. Samples of DNA were collected from January 1, 1993, to December 31, 2009, and analyzed from March 1, 2014, to February 28, 2015.

MAIN OUTCOMES AND MEASURES

Cox proportional hazards regression under an additive model with adjustment for age at onset, sex, and the first 4 principal components of ancestry, followed by meta-analysis, were used to analyze data. Survival distributions for the most associated genetic variants were assessed by Kaplan-Meier analysis.

RESULTS

Among the 4256 patients included in the analysis (2589 male [60.8%] and 1667 female [39.2%]; mean [SD] age at onset, 59 [12] years), the following 2 novel loci were significantly associated with ALS survival: at 10q23 (rs139550538; P = 1.87 × 10-9) and in the CAMTA1 gene at 1p36 (rs2412208, P = 3.53 × 10-8). At locus 10q23, the adjusted hazard ratio for patients with the rs139550538 AA or AT genotype was 1.61 (95% CI, 1.38-1.89; P = 1.87 × 10-9), corresponding to an 8-month reduction in survival compared with TT carriers. For rs2412208 CAMTA1, the adjusted hazard ratio for patients with the GG or GT genotype was 1.17 (95% CI, 1.11-1.24; P = 3.53 × 10-8), corresponding to a 4-month reduction in survival compared with TT carriers.

CONCLUSIONS AND RELEVANCE

This GWAS robustly identified 2 loci at genome-wide levels of significance that influence survival in patients with ALS. Because ALS is a rare disease and prevention is not feasible, treatment that modifies survival is the most realistic strategy. Therefore, identification of modifier genes that might influence ALS survival could improve the understanding of the biology of the disease and suggest biological targets for pharmaceutical intervention. In addition, genetic risk scores for survival could be used as an adjunct to clinical trials to account for the genetic contribution to survival.

GWAS meta-analysis of 16 852 women identifies new susceptibility locus for endometrial cancer

Endometrial cancer is the most common gynecological malignancy in the developed world. Although there is evidence of genetic predisposition to the disease, most of the genetic risk remains unexplained. We present the meta-analysis results of four genome-wide association studies (4907 cases and 11 945 controls total) in women of European ancestry. We describe one new locus reaching genome-wide significance (P < 5 × 10 -8) at 6p22.3 (rs1740828; P = 2.29 × 10 -8, OR = 1.20), providing evidence of an additional region of interest for genetic susceptibility to endometrial cancer.

An Empirical Comparison of Joint and Stratified Frameworks for Studying G × E Interactions: Systolic Blood Pressure and Smoking in the CHARGE Gene-Lifestyle Interactions Working Group

Studying gene-environment (G × E) interactions is important, as they extend our knowledge of the genetic architecture of complex traits and may help to identify novel variants not detected via analysis of main effects alone. The main statistical framework for studying G × E interactions uses a single regression model that includes both the genetic main and G × E interaction effects (the "joint" framework). The alternative "stratified" framework combines results from genetic main-effect analyses carried out separately within the exposed and unexposed groups. Although there have been several investigations using theory and simulation, an empirical comparison of the two frameworks is lacking. Here, we compare the two frameworks using results from genome-wide association studies of systolic blood pressure for 3.2 million low frequency and 6.5 million common variants across 20 cohorts of European ancestry, comprising 79,731 individuals. Our cohorts have sample sizes ranging from 456 to 22,983 and include both family-based and population-based samples. In cohort-specific analyses, the two frameworks provided similar inference for population-based cohorts. The agreement was reduced for family-based cohorts. In meta-analyses, agreement between the two frameworks was less than that observed in cohort-specific analyses, despite the increased sample size. In meta-analyses, agreement depended on (1) the minor allele frequency, (2) inclusion of family-based cohorts in meta-analysis, and (3) filtering scheme. The stratified framework appears to approximate the joint framework well only for common variants in population-based cohorts. We conclude that the joint framework is the preferred approach and should be used to control false positives when dealing with low-frequency variants and/or family-based cohorts.

The GenABEL Project for statistical genomics

Development of free/libre open source software is usually done by a community of people with an interest in the tool. For scientific software, however, this is less often the case. Most scientific software is written by only a few authors, often a student working on a thesis. Once the paper describing the tool has been published, the tool is no longer developed further and is left to its own device. Here we describe the broad, multidisciplinary community we formed around a set of tools for statistical genomics. The GenABEL project for statistical omics actively promotes open interdisciplinary development of statistical methodology and its implementation in efficient and user-friendly software under an open source licence. The software tools developed withing the project collectively make up the GenABEL suite, which currently consists of eleven tools. The open framework of the project actively encourages involvement of the community in all stages, from formulation of methodological ideas to application of software to specific data sets. A web forum is used to channel user questions and discussions, further promoting the use of the GenABEL suite. Developer discussions take place on a dedicated mailing list, and development is further supported by robust development practices including use of public version control, code review and continuous integration. Use of this open science model attracts contributions from users and developers outside the “core team”, facilitating agile statistical omics methodology development and fast dissemination.

Genome-wide association study of biologically informed periodontal complex traits offers novel insights into the genetic basis of periodontal disease

Genome-wide association studies (GWAS) of chronic periodontitis (CP) defined by clinical criteria alone have had modest success to-date. Here, we refine the CP phenotype by supplementing clinical data with biological intermediates of microbial burden (levels of eight periodontal pathogens) and local inflammatory response (gingival crevicular fluid IL-1β) and derive periodontal complex traits (PCTs) via principal component analysis. PCTs were carried forward to GWAS (∼2.5 million markers) to identify PCT-associated loci among 975 European American adult participants of the Dental ARIC study. We sought to validate these findings for CP in the larger ARIC cohort (n = 821 participants with severe CP, 2031-moderate CP, 1914-healthy/mild disease) and an independent German sample including 717 aggressive periodontitis cases and 4210 controls. We identified six PCTs with distinct microbial community/IL-1β structures, although with overlapping clinical presentations. PCT1 was characterized by a uniformly high pathogen load, whereas PCT3 and PCT5 were dominated by Aggregatibacter actinomycetemcomitans and Porphyromonas gingivalis, respectively. We detected genome-wide significant signals for PCT1 (CLEC19A, TRA, GGTA2P, TM9SF2, IFI16, RBMS3), PCT4 (HPVC1) and PCT5 (SLC15A4, PKP2, SNRPN). Overall, the highlighted loci included genes associated with immune response and epithelial barrier function. With the exception of associations of BEGAIN with severe and UBE3D with moderate CP, no other loci were associated with CP in ARIC or aggressive periodontitis in the German sample. Although not associated with current clinically determined periodontal disease taxonomies, upon replication and mechanistic validation these candidate loci may highlight dysbiotic microbial community structures and altered inflammatory/immune responses underlying biological sub-types of CP.

Antiphospholipid antibodies in a large population-based cohort: genome-wide associations and effects on monocyte gene expression

The antiphospholipid syndrome (APS) is characterised by venous and/or arterial thrombosis and pregnancy morbidity in women combined with the persistent presence of antiphospholipid antibodies (aPL). We aimed to identify genetic factors associated with the presence of aPL in a population based cohort. Furthermore, we wanted to clarify if the presence of aPL affects gene expression in circulating monocytes. Titres of IgG and IgM against cardiolipin, β2glycoprotein 1 (anti-β2GPI), and IgG against domain 1 of β2GPI (anti-domain 1) were determined in approx. 5,000 individuals from the Gutenberg Health Study (GHS) a population based cohort of German descent. Genotyping was conducted on Affymetrix Genome-Wide Human SNP 6.0 arrays. Monocyte gene expression was determined in a subgroup of 1,279 individuals by using the Illumina HT-12 v3 BeadChip. Gene expression data were confirmed in vitro and ex vivo by qRT-PCR. Genome wide analysis revealed significant associations of anti-β2GPI IgG and APOH on chromosome 17, which had been previously identified by candidate gene approaches, and of anti-domain1 and MACROD2 on chromosome 20 which has been listed in a previous GWAS as a suggestive locus associated with the occurrence of anti-β2GPI antibodies. Expression analysis confirmed increased expression of TNFα in monocytes and identified and confirmed neuron navigator 3 (NAV3) as a novel gene induced by aPL. In conclusion, MACROD2 represents a novel genetic locus associated with aPL. Furthermore, we show that aPL induce the expression of NAV3 in monocytes and endothelial cells. This will stimulate further research into the role of these genes in the APS.

Tracking Resilience to Infections by Mapping Disease Space

Infected hosts differ in their responses to pathogens; some hosts are resilient and recover their original health, whereas others follow a divergent path and die. To quantitate these differences, we propose mapping the routes infected individuals take through “disease space.” We find that when plotting physiological parameters against each other, many pairs have hysteretic relationships that identify the current location of the host and predict the future route of the infection. These maps can readily be constructed from experimental longitudinal data, and we provide two methods to generate the maps from the cross-sectional data that is commonly gathered in field trials. We hypothesize that resilient hosts tend to take small loops through disease space, whereas nonresilient individuals take large loops. We support this hypothesis with experimental data in mice infected with Plasmodium chabaudi, finding that dying mice trace a large arc in red blood cells (RBCs) by reticulocyte space as compared to surviving mice. We find that human malaria patients who are heterozygous for sickle cell hemoglobin occupy a small area of RBCs by reticulocyte space, suggesting this approach can be used to distinguish resilience in human populations. This technique should be broadly useful in describing the in-host dynamics of infections in both model hosts and patients at both population and individual levels.

This study shows that infections cause sick hosts to loop through disease space on their return back to health; resilient individuals take tiny loops through this space, minimizing the impact of the infection.

When we get sick, we long for recovery; thus, a major goal of medicine is to promote resilience—the ability of a host to return to its original health following an infection. While in the laboratory we can study the response to infection with precise knowledge of inoculation time and dose, sick patients in the clinic do not have this information. This creates a problem because we can’t easily differentiate between patients who are early in the stages of infection that will develop severe disease from more disease-tolerant patients who present later in the infection. The distinction between these two types of patients is important, as the less disease-tolerant patient would require a more aggressive treatment regime. To determine where patients lie along the infection timeline, we charted “disease maps” that trace a patient’s route through “disease space.” We select symptoms that produce looping graphs as patients grow sick and recover. Using a mouse–malaria model, we demonstrate that less resilient individuals take wider loops through this space, representing a longer infection time with more severe symptoms. We find this looping behavior also applies to humans and suggest that people carrying the sickle cell trait are more resilient to malaria infections.

Genome-wide association and Mendelian randomization study of NT-proBNP in patients with acute coronary syndrome

N-terminal pro-B-type natriuretic peptide (NT-proBNP) is a strong predictor of mortality in coronary artery disease and is widely employed as a prognostic biomarker. However, a causal relationship between NT-proBNP and clinical endpoints has not been established. We have performed a genome-wide association and Mendelian randomization study of NT-proBNP. We used a discovery set of 3740 patients from the PLATelet inhibition and patient Outcomes (PLATO) trial, which enrolled 18 624 patients with acute coronary syndrome (ACS). A further set of 5492 patients, from the same trial, was used for replication. Genetic variants at two novel loci (SLC39A8 and POC1B/GALNT4) were associated with NT-proBNP levels and replicated together with the previously known NPPB locus. The most significant SNP (rs198389, pooled P = 1.07 × 10(-15)) in NPPB interrupts an E-box consensus motif in the gene promoter. The association in SLC39A8 is driven by a deleterious variant (rs13107325, pooled P = 5.99 × 10(-10)), whereas the most significant SNP in POC1B/GALNT4 (rs11105306, pooled P = 1.02 × 10(-16)) is intronic. The SLC39A8 SNP was associated with higher risk of cardiovascular (CV) death (HR = 1.39, 95% CI: 1.08-1.79, P = 0.0095), but the other loci were not associated with clinical endpoints. We have identified two novel loci to be associated with NT-proBNP in patients with ACS. Only the SLC39A8 variant, but not the NPPB variant, was associated with a clinical endpoint. Due to pleotropic effects of SLC39A8, these results do not suggest that NT-proBNP levels have a direct effect on mortality in ACS patients. PLATO Clinical Trial Registration: www.clinicaltrials.gov; NCT00391872.

GENOME-WIDE ASSOCIATION STUDY (GWAS) AND GENOME-WIDE BY ENVIRONMENT INTERACTION STUDY (GWEIS) OF DEPRESSIVE SYMPTOMS IN AFRICAN AMERICAN AND HISPANIC/LATINA WOMEN

BACKGROUND

Genome-wide association studies (GWAS) have made little progress in identifying variants linked to depression. We hypothesized that examining depressive symptoms and considering gene-environment interaction (GxE) might improve efficiency for gene discovery. We therefore conducted a GWAS and genome-wide by environment interaction study (GWEIS) of depressive symptoms.

METHODS

Using data from the SHARe cohort of the Women's Health Initiative, comprising African Americans (n = 7,179) and Hispanics/Latinas (n = 3,138), we examined genetic main effects and GxE with stressful life events and social support. We also conducted a heritability analysis using genome-wide complex trait analysis (GCTA). Replication was attempted in four independent cohorts.

RESULTS

No SNPs achieved genome-wide significance for main effects in either discovery sample. The top signals in African Americans were rs73531535 (located 20 kb from GPR139, P = 5.75 × 10(-8) ) and rs75407252 (intronic to CACNA2D3, P = 6.99 × 10(-7) ). In Hispanics/Latinas, the top signals were rs2532087 (located 27 kb from CD38, P = 2.44 × 10(-7) ) and rs4542757 (intronic to DCC, P = 7.31 × 10(-7) ). In the GEWIS with stressful life events, one interaction signal was genome-wide significant in African Americans (rs4652467; P = 4.10 × 10(-10) ; located 14 kb from CEP350). This interaction was not observed in a smaller replication cohort. Although heritability estimates for depressive symptoms and stressful life events were each less than 10%, they were strongly genetically correlated (rG = 0.95), suggesting that common variation underlying self-reported depressive symptoms and stressful life event exposure, though modest on their own, were highly overlapping in this sample.

CONCLUSIONS

Our results underscore the need for larger samples, more GEWIS, and greater investigation into genetic and environmental determinants of depressive symptoms in minorities.

Genome-wide significant results identified for plasma apolipoprotein H levels in middle-aged and older adults

Apolipoprotein H (ApoH) is a multi-functional plasma glycoprotein that has been associated with negative health outcomes. ApoH levels have high heritability. We undertook a genome-wide association study of ApoH levels using the largest sample to date and replicated the results in an independent cohort (total N = 1,255). In the discovery phase, a meta-analysis of two cohorts, the Sydney Memory and Ageing Study (Sydney MAS) and the Older Australian Twins Study (OATS) (n = 942) revealed genome-wide significant results in or near the APOH gene on chromosome 17 (top SNP, rs7211380, p = 1 × 10⁻¹¹). The results were replicated in an independent cohort, the Hunter Community Study (p < 0.002) (n = 313). Conditional and joint analysis (COJO) confirmed the association of the chromosomal 17 region with ApoH levels. The set of independent SNPs identified by COJO explained 23% of the variance. The relationships between the top SNPs and cardiovascular/lipid/cognition measures and diabetes were assessed in Sydney MAS, with suggestive results observed for diabetes and cognitive performance. However, replication of these results in the smaller OATS cohort was not found. This work provides impetus for future research to better understand the contribution of genetics to ApoH levels and its possible impacts on health.

CYP19A1 fine-mapping and Mendelian randomization: estradiol is causal for endometrial cancer

Candidate gene studies have reported CYP19A1 variants to be associated with endometrial cancer and with estradiol (E2) concentrations. We analyzed 2937 single nucleotide polymorphisms (SNPs) in 6608 endometrial cancer cases and 37 925 controls and report the first genome wide-significant association between endometrial cancer and a CYP19A1 SNP (rs727479 in intron 2, P=4.8×10(-11)). SNP rs727479 was also among those most strongly associated with circulating E2 concentrations in 2767 post-menopausal controls (P=7.4×10(-8)). The observed endometrial cancer odds ratio per rs727479 A-allele (1.15, CI=1.11-1.21) is compatible with that predicted by the observed effect on E2 concentrations (1.09, CI=1.03-1.21), consistent with the hypothesis that endometrial cancer risk is driven by E2. From 28 candidate-causal SNPs, 12 co-located with three putative gene-regulatory elements and their risk alleles associated with higher CYP19A1 expression in bioinformatical analyses. For both phenotypes, the associations with rs727479 were stronger among women with a higher BMI (Pinteraction=0.034 and 0.066 respectively), suggesting a biologically plausible gene-environment interaction.

Genome-wide association analysis identifies TXNRD2, ATXN2 and FOXC1 as susceptibility loci for primary open-angle glaucoma

Primary open angle glaucoma (POAG) is a leading cause of blindness world-wide. To identify new susceptibility loci, we meta-analyzed GWAS results from 8 independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significant SNPs in two Australian studies (1,252 cases and 2,592 controls), 3 European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of top SNPs identified three novel loci: rs35934224[T] within TXNRD2 (odds ratio (OR) = 0.78, P = 4.05×10⁻¹¹ encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] within ATXN2 (OR = 1.17, P = 8.73×10⁻¹⁰), and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76×10⁻¹⁰). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest novel targets for preventative therapies.

Independent Replication of Published Germline Polymorphisms Associated with Urinary Bladder Cancer Prognosis and Treatment Response

BACKGROUND

Many studies investigated the prognostic or predictive relevance of single nucleotide polymorphisms (SNPs) in biologically plausible genes in urinary bladder cancer (UBC) patients. Most published SNP associations have never been replicated in independent patient series.

OBJECTIVE

To independently replicate all previously reported associations between germline SNPs and disease prognosis or treatment response in UBC.

METHODS

A Pubmed search was performed to identify studies published by July 1, 2014 reporting on germline SNP associations with UBC prognosis or treatment response. For the replication series, consisting of 1,284 non-muscle-invasive bladder cancer (NMIBC) and 275 muscle-invasive or metastatic bladder cancer (MIBC) patients recruited through the Netherlands Cancer Registry, detailed clinical data were retrieved from medical charts. Patients were genotyped using a genome-wide SNP array. SNP association with recurrence-free, progression-free, and overall survival (OS) within specific patient and treatment strata was tested using Cox regression analyses.

RESULTS

For only six of the 114 evaluated SNPs, the association with either UBC prognosis or treatment response was replicated at the p <  0.05 level: rs1799793 (ERCC2) and rs187238 (IL18) for BCG recurrence; rs6678136 (RGS4) and rs11585883 (RGS5) for NMIBC progression; rs12035879 (RGS5) and rs2075786 (TERT) for MIBC OS.

CONCLUSIONS

Non-replicated genetic associations in the literature require cautious interpretation. This single replication does not provide definitive proof of association for the six SNPs, and non-replication of other SNPs may result from population-specific effects or the retrospective patient enrollment.

Genome-wide meta-analyses identify novel loci associated with n-3 and n-6 polyunsaturated fatty acid levels in Chinese and European-ancestry populations

Epidemiological studies suggest that levels of n-3 and n-6 long-chain polyunsaturated fatty acids are associated with risk of cardio-metabolic outcomes across different ethnic groups. Recent genome-wide association studies in populations of European ancestry have identified several loci associated with plasma and/or erythrocyte polyunsaturated fatty acids. To identify additional novel loci, we carried out a genome-wide association study in two population-based cohorts consisting of 3521 Chinese participants, followed by a trans-ethnic meta-analysis with meta-analysis results from 8962 participants of European ancestry. Four novel loci (MYB, AGPAT4, DGAT2 and PPT2) reached genome-wide significance in the trans-ethnic meta-analysis (log10(Bayes Factor) ≥ 6). Of them, associations of MYB and AGPAT4 with docosatetraenoic acid (log10(Bayes Factor) = 11.5 and 8.69, respectively) also reached genome-wide significance in the Chinese-specific genome-wide association analyses (P = 4.15 × 10(-14) and 4.30 × 10(-12), respectively), while associations of DGAT2 with gamma-linolenic acid (log10(Bayes Factor) = 6.16) and of PPT2 with docosapentaenoic acid (log10(Bayes Factor) = 6.24) were nominally significant in both Chinese- and European-specific genome-wide association analyses (P ≤ 0.003). We also confirmed previously reported loci including FADS1, NTAN1, NRBF2, ELOVL2 and GCKR. Different effect sizes in FADS1 and independent association signals in ELOVL2 were observed. These results provide novel insight into the genetic background of polyunsaturated fatty acids and their differences between Chinese and European populations.

A Novel Risk Locus at 6p21.3 for Epstein-Barr Virus-Positive Hodgkin Lymphoma

BACKGROUND

A proportion of the genetic variants involved in susceptibility to Hodgkin lymphoma differ by the tumor's Epstein-Barr virus (EBV) status, particularly within the MHC region.

METHODS

We have conducted an SNP imputation study of the MHC region, considering tumor EBV status in 1,200 classical Hodgkin lymphoma (cHL) cases and 5,726 control subjects of European origin. Notable findings were genotyped in an independent study population of 468 cHL cases and 551 controls.

RESULTS

We identified and subsequently replicated a novel association between a common genetic variant rs6457715 and cHL. Although strongly associated with EBV-positive cHL [OR, 2.33; 95% confidence interval (CI), 1.83-2.97; P = 7 × 10(-12)], there was little evidence for association between rs6457715 and the EBV-negative subgroup of cHL (OR, 1.06; 95% CI, 0.92-1.21), indicating that this association was specific to the EBV-positive subgroup (Phet < P = 10(-8)). Furthermore, the association was limited to EBV-positive cHL subgroups within mixed cell (MCHL) and nodular sclerosis subtypes (NSHL), suggesting that the association is independent of histologic subtype of cHL.

CONCLUSIONS

rs6457715, located near the HLA-DPB1 gene, is associated with EBV-positive cHL and suggests this region as a novel susceptibility locus for cHL.

IMPACT

This expands the number of genetic variants that are associated with cHL and provides additional evidence for a critical and specific role of EBV in the etiology of this disease.

Calcium intake, polymorphisms of the calcium-sensing receptor, and recurrent/aggressive prostate cancer

PURPOSE

To assess whether calcium intake and common genetic variants of the calcium-sensing receptor (CASR) are associated with either aggressive prostate cancer (PCa) or disease recurrence after prostatectomy.

METHODS

Calcium intake at diagnosis was assessed, and 65 common single-nucleotide polymorphisms (SNPs) in CASR were genotyped in 886 prostatectomy patients. We investigated the association between calcium intake and CASR variants with both PCa recurrence and aggressiveness (defined as Gleason score ≥4 + 3, stage ≥pT3, or nodal-positive disease).

RESULTS

A total of 285 men had aggressive disease and 91 experienced recurrence. A U-shaped relationship between calcium intake and both disease recurrence and aggressiveness was observed. Compared to the middle quintile, the HR for disease recurrence was 3.07 (95% CI 1.41-6.69) for the lowest quintile and 3.21 (95% CI 1.47-7.00) and 2.97 (95% CI 1.37-6.45) for the two upper quintiles, respectively. Compared to the middle quintile, the OR for aggressive disease was 1.80 (95% CI 1.11-2.91) for the lowest quintile and 1.75 (95% CI 1.08-2.85) for the highest quintile of calcium intake. The main effects of CASR variants were not associated with PCa recurrence or aggressiveness. In the subgroup of patients with moderate calcium intake, 31 SNPs in four distinct blocks of high linkage disequilibrium were associated with PCa recurrence.

CONCLUSIONS

We observed a protective effect of moderate calcium intake for PCa aggressiveness and recurrence. While CASR variants were not associated with these outcomes in the entire cohort, they may be associated with disease recurrence in men with moderate calcium intakes.

Genome-wide Meta-analysis on the Sense of Smell Among US Older Adults

Olfactory dysfunction is common among older adults and affects their safety, nutrition, quality of life, and mortality. More importantly, the decreased sense of smell is an early symptom of neurodegenerative diseases such as Parkinson disease (PD) and Alzheimer disease. However, the genetic determinants for the sense of smell have been poorly investigated. We here performed the first genome-wide meta-analysis on the sense of smell among 6252 US older adults of European descent from the Atherosclerosis Risk in Communities (ARIC) study, the Health, Aging, and Body Composition (Health ABC) study, and the Religious Orders Study and the Rush Memory and Aging Project (ROS/MAP). Genome-wide association study analysis was performed first by individual cohorts and then meta-analyzed using fixed-effect models with inverse variance weights. Although no SNPs reached genome-wide statistical significance, we identified 13 loci with suggestive evidence for an association with the sense of smell (Pmeta < 1 × 10). Of these, 2 SNPs at chromosome 17q21.31 (rs199443 in NSF, P = 3.02 × 10; and rs2732614 in KIAA1267-LRRC37A, P = 6.65 × 10) exhibited cis effects on the expression of microtubule-associated protein tau (MAPT, 17q21.31) in 447 frontal-cortex samples obtained postmortem and profiled by RNA-seq (P < 1 × 10). Gene-based and pathway-enrichment analyses further implicated MAPT in regulating the sense of smell in older adults. Similar results were obtained after excluding participants who reported a physician-diagnosed PD or use of PD medications. In conclusion, we provide preliminary evidence that the MAPT locus may play a role in regulating the sense of smell in older adults and therefore offer a potential genetic link between poor sense of smell and major neurodegenerative diseases.

Quality control, imputation and analysis of genome-wide genotyping data from the Illumina HumanCoreExome microarray

The decreasing cost of performing genome-wide association studies has made genomics widely accessible. However, there is a paucity of guidance for best practice in conducting such analyses. For the results of a study to be valid and replicable, multiple biases must be addressed in the course of data preparation and analysis. In addition, standardizing methods across small, independent studies would increase comparability and the potential for effective meta-analysis. This article provides a discussion of important aspects of quality control, imputation and analysis of genome-wide data from a low-coverage microarray, as well as a straight-forward guide to performing a genome-wide association study. A detailed protocol is provided online, with example scripts available at https://github.com/JoniColeman/gwas_scripts.

Genome-wide meta-analysis reveals common splice site acceptor variant in CHRNA4 associated with nicotine dependence

We conducted a 1000 Genomes-imputed genome-wide association study (GWAS) meta-analysis for nicotine dependence, defined by the Fagerström Test for Nicotine Dependence in 17 074 ever smokers from five European-ancestry samples. We followed up novel variants in 7469 ever smokers from five independent European-ancestry samples. We identified genome-wide significant association in the alpha-4 nicotinic receptor subunit (CHRNA4) gene on chromosome 20q13: lowest P=8.0 × 10(-9) across all the samples for rs2273500-C (frequency=0.15; odds ratio=1.12 and 95% confidence interval=1.08-1.17 for severe vs mild dependence). rs2273500-C, a splice site acceptor variant resulting in an alternate CHRNA4 transcript predicted to be targeted for nonsense-mediated decay, was associated with decreased CHRNA4 expression in physiologically normal human brains (lowest P=7.3 × 10(-4)). Importantly, rs2273500-C was associated with increased lung cancer risk (N=28 998, odds ratio=1.06 and 95% confidence interval=1.00-1.12), likely through its effect on smoking, as rs2273500-C was no longer associated with lung cancer after adjustment for smoking. Using criteria for smoking behavior that encompass more than the single 'cigarettes per day' item, we identified a common CHRNA4 variant with important regulatory properties that contributes to nicotine dependence and smoking-related consequences.

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, WHRadjBMI), 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 WHRadjBMI, 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 WHRadjBMI 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 WHRadjBMI, 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 WHRadjBMI. 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 WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

International genome-wide meta-analysis identifies new primary biliary cirrhosis risk loci and targetable pathogenic pathways

Primary biliary cirrhosis (PBC) is a classical autoimmune liver disease for which effective immunomodulatory therapy is lacking. Here we perform meta-analyses of discovery data sets from genome-wide association studies of European subjects (n=2,764 cases and 10,475 controls) followed by validation genotyping in an independent cohort (n=3,716 cases and 4,261 controls). We discover and validate six previously unknown risk loci for PBC (Pcombined<5 × 10(-8)) and used pathway analysis to identify JAK-STAT/IL12/IL27 signalling and cytokine-cytokine pathways, for which relevant therapies exist.

The role of rare variants in systolic blood pressure: analysis of ExomeChip data in HyperGEN African Americans

Cardiovascular diseases are among the most significant health problems in the United States today, with their major risk factor, hypertension, disproportionately affecting African Americans (AAs). Although GWAS have identified dozens of common variants associated with blood pressure (BP) and hypertension in European Americans, these variants collectively explain <2.5% of BP variance, and most of the genetic variants remain yet to be identified. Here, we report the results from rare-variant analysis of systolic BP using 94,595 rare and low-frequency variants (minor allele frequency, MAF, <5%) from the Illumina exome array genotyped in 2,045 HyperGEN AAs. In addition to single-variant analysis, 4 gene-level association tests were used for analysis: burden and family-based SKAT tests using MAF cutoffs of 1 and 5%. The gene-based methods often provided lower p values than the single-variant approach. Some consistency was observed across these 4 gene-based analysis options. While neither the gene-based analyses nor the single-variant analysis produced genome-wide significant results, the top signals, which had supporting evidence from multiple gene-based methods, were of borderline significance. Though additional molecular validations are required, 6 of the 16 most promising genes are biologically plausible with physiological connections to BP regulation.

A guide to genome-wide association analysis and post-analytic interrogation

This tutorial is a learning resource that outlines the basic process and provides specific software tools for implementing a complete genome‐wide association analysis. Approaches to post‐analytic visualization and interrogation of potentially novel findings are also presented. Applications are illustrated using the free and open‐source R statistical computing and graphics software environment, Bioconductor software for bioinformatics and the UCSC Genome Browser. Complete genome‐wide association data on 1401 individuals across 861,473 typed single nucleotide polymorphisms from the PennCATH study of coronary artery disease are used for illustration. All data and code, as well as additional instructional resources, are publicly available through the Open Resources in Statistical Genomics project: http://www.stat-gen.org. © 2015 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.

Genome-wide association of polycystic ovary syndrome implicates alterations in gonadotropin secretion in European ancestry populations

Polycystic ovary syndrome (PCOS) is a common, highly heritable complex disorder of unknown aetiology characterized by hyperandrogenism, chronic anovulation and defects in glucose homeostasis. Increased luteinizing hormone relative to follicle-stimulating hormone secretion, insulin resistance and developmental exposure to androgens are hypothesized to play a causal role in PCOS. Here we map common genetic susceptibility loci in European ancestry women for the National Institutes of Health PCOS phenotype, which confers the highest risk for metabolic morbidities, as well as reproductive hormone levels. Three loci reach genome-wide significance in the case–control meta-analysis, two novel loci mapping to chr 8p32.1 and chr 11p14.1, and a chr 9q22.32 locus previously found in Chinese PCOS. The same chr 11p14.1 SNP, rs11031006, in the region of the follicle-stimulating hormone B polypeptide (FSHB) gene strongly associates with PCOS diagnosis and luteinizing hormone levels. These findings implicate neuroendocrine changes in disease pathogenesis.

Polycystic Ovary Sydrome is a highly heritable, complex reproductive disorder with unknown underlying genetic factors. Here Hayes and Urbanek et al. identify three loci in European women strongly associated with neuroendocrine changes and disease susceptibility.

A multiancestry study identifies novel genetic associations with CHRNA5 methylation in human brain and risk of nicotine dependence

Nicotine dependence is influenced by chromosome 15q25.1 single nucleotide polymorphisms (SNPs), including the missense SNP rs16969968 that alters function of the α5 nicotinic acetylcholine receptor (CHRNA5) and noncoding SNPs that regulate CHRNA5 mRNA expression. We tested for cis-methylation quantitative trait loci (cis-meQTLs) using SNP genotypes and DNA methylation levels measured across the IREB2-HYKK-PSMA4-CHRNA5-CHRNA3-CHRNB4 genes on chromosome 15q25.1 in the BrainCloud and Brain QTL cohorts [total N = 175 European-Americans and 65 African-Americans (AAs)]. We identified eight SNPs that were significantly associated with CHRNA5 methylation in prefrontal cortex: P ranging from 6.0 × 10(-10) to 5.6 × 10(-5). These SNP-methylation associations were also significant in frontal cortex, temporal cortex and pons: P ranging from 4.8 × 10(-12) to 3.4 × 10(-3). Of the eight cis-meQTL SNPs, only the intronic CHRNB4 SNP rs11636753 was associated with CHRNA5 methylation independently of the known SNP effects in prefrontal cortex, and it was the most significantly associated SNP with nicotine dependence across five independent cohorts (total N = 7858 European ancestry and 3238 AA participants): P = 6.7 × 10(-4), odds ratio (OR) [95% confidence interval (CI)] = 1.11 (1.05-1.18). The rs11636753 major allele (G) was associated with lower CHRNA5 DNA methylation, lower CHRNA5 mRNA expression and increased nicotine dependence risk. Haplotype analyses showed that rs11636753-G and the functional rs16969968-A alleles together increased risk of nicotine dependence more than each variant alone: P = 3.1 × 10(-12), OR (95% CI) = 1.32 (1.22-1.43). Our findings identify a novel regulatory SNP association with nicotine dependence and connect, for the first time, previously observed differences in CHRNA5 mRNA expression and nicotine dependence risk to underlying DNA methylation differences.

Layered genetic control of DNA methylation and gene expression: a locus of multiple sclerosis in healthy individuals

DNA methylation may contribute to the etiology of complex genetic disorders through its impact on genome integrity and gene expression; it is modulated by DNA-sequence variants, named methylation quantitative trait loci (meQTLs). Most meQTLs influence methylation of a few CpG dinucleotides within short genomic regions (<3 kb). Here, we identified a layered genetic control of DNA methylation at numerous CpGs across a long 300 kb genomic region. This control involved a single long-range meQTL and multiple local meQTLs. The long-range meQTL explained up to 75% of variance in methylation of CpGs located over extended areas of the 300 kb region. The meQTL was identified in four samples (P = 2.8 × 10(-17), 3.1 × 10(-31), 4.0 × 10(-71) and 5.2 × 10(-199)), comprising a total of 2796 individuals. The long-range meQTL was strongly associated not only with DNA methylation but also with mRNA expression of several genes within the 300 kb region (P = 7.1 × 10(-18)-1.0 × 10(-123)). The associations of the meQTL with gene expression became attenuated when adjusted for DNA methylation (causal inference test: P = 2.4 × 10(-13)-7.1 × 10(-20)), indicating coordinated regulation of DNA methylation and gene expression. Further, the long-range meQTL was found to be in linkage disequilibrium with the most replicated locus of multiple sclerosis, a disease affecting primarily the brain white matter. In middle-aged adults free of the disease, we observed that the risk allele was associated with subtle structural properties of the brain white matter found in multiple sclerosis (P = 0.02). In summary, we identified a long-range meQTL that controls methylation and expression of several genes and may be involved in increasing brain vulnerability to multiple sclerosis.

A Joint Location-Scale Test Improves Power to Detect Associated SNPs, Gene Sets, and Pathways

Gene-based, pathway, and other multivariate association methods are motivated by the possibility of GxG and GxE interactions; however, accounting for such interactions is limited by the challenges associated with adequate modeling information. Here we propose an easy-to-implement joint location-scale (JLS) association testing framework for single-variant and multivariate analysis that accounts for interactions without explicitly modeling them. We apply the JLS method to a gene-set analysis of cystic fibrosis (CF) lung disease, which is influenced by multiple environmental and genetic factors. We identify and replicate an association between the constituents of the apical plasma membrane and CF lung disease (p = 0.0099 and p = 0.0180, respectively) and highlight a role for the SLC9A3-SLC9A3R1/2-EZR complex in contributing to CF lung disease. Many association studies could benefit from re-analysis with the JLS method that leverages complex genetic architecture for SNP, gene, and pathway identification. Analytical verification, simulation, and additional proof-of-principle applications support our approach.

Genetic variants in the ADAMTS13 and SUPT3H genes are associated with ADAMTS13 activity

We identify rs41314453 as the strongest genetic predictor of ADAMTS13 activity, associated with a decrease of >20%. We present evidence of further independent associations with a common variant in SUPT3H, as well as 5 variants at the ADAMTS13 locus.

Genome-Wide Meta-Analysis of Longitudinal Alcohol Consumption Across Youth and Early Adulthood

The public health burden of alcohol is unevenly distributed across the life course, with levels of use, abuse, and dependence increasing across adolescence and peaking in early adulthood. Here, we leverage this temporal patterning to search for common genetic variants predicting developmental trajectories of alcohol consumption. Comparable psychiatric evaluations measuring alcohol consumption were collected in three longitudinal community samples (N=2,126, obs=12,166). Consumption-repeated measurements spanning adolescence and early adulthood were analyzed using linear mixed models, estimating individual consumption trajectories, which were then tested for association with Illumina 660W-Quad genotype data (866,099 SNPs after imputation and QC). Association results were combined across samples using standard meta-analysis methods. Four meta-analysis associations satisfied our pre-determined genome-wide significance criterion (FDR<0.1) and six others met our 'suggestive' criterion (FDR<0.2). Genome-wide significant associations were highly biological plausible, including associations within GABA transporter 1, SLC6A1 (solute carrier family 6, member 1), and exonic hits in LOC100129340 (mitofusin-1-like). Pathway analyses elaborated single marker results, indicating significant enriched associations to intuitive biological mechanisms, including neurotransmission, xenobiotic pharmacodynamics, and nuclear hormone receptors (NHR). These findings underscore the value of combining longitudinal behavioral data and genome-wide genotype information in order to study developmental patterns and improve statistical power in genomic studies.

WNT10A exonic variant increases the risk of keratoconus by decreasing corneal thickness

Keratoconus is a degenerative eye condition which results from thinning of the cornea and causes vision distortion. Treatments such as ultraviolet (UV) cross-linking have proved effective for management of keratoconus when performed in early stages of the disease. The central corneal thickness (CCT) is a highly heritable endophenotype of keratoconus, and it is estimated that up to 95% of its phenotypic variance is due to genetics. Genome-wide association efforts of CCT have identified common variants (i.e. minor allele frequency (MAF) >5%). However, these studies typically ignore the large set of exonic variants whose MAF is usually low. In this study, we performed a CCT exome-wide association analysis in a sample of 1029 individuals from a population-based study in Western Australia. We identified a genome-wide significant exonic variant rs121908120 (P = 6.63 × 10(-10)) in WNT10A. This gene is 437 kb from a gene previously associated with CCT (USP37). We showed in a conditional analysis that the WNT10A variant completely accounts for the signal previously seen at USP37. We replicated our finding in independent samples from the Brisbane Adolescent Twin Study, Twin Eye Study in Tasmania and the Rotterdam Study. Further, we genotyped rs121908120 in 621 keratoconus cases and compared the frequency to a sample of 1680 unscreened controls from the Queensland Twin Registry. We found that rs121908120 increases the risk of keratoconus two times (odds ratio 2.03, P = 5.41 × 10(-5)).

Genetic contributions to urgency urinary incontinence in women

PURPOSE

We identify genetic variants associated with urgency urinary incontinence in postmenopausal women.

MATERIALS AND METHODS

A 2-stage genome-wide association analysis was conducted to identify variants associated with urgency urinary incontinence. The WHI GARNET substudy with 4,894 genotyped post-reproductive white women was randomly split into independent discovery and replication cohorts. Genome-wide imputation was performed using IMPUTE2 with the 1000 Genomes ALL Phase I integrated variant set as a reference. Controls reported no urgency urinary incontinence at enrollment or followup. Cases reported monthly or greater urgency urinary incontinence and leaked sufficiently to wet/soak underpants/clothes. Logistic regression models were used to predict urgency urinary incontinence case vs control status based on genotype, assuming additive inheritance. Age, obesity, diabetes and depression were included in the models as covariates.

RESULTS

Following quality control, 975,508 single nucleotide polymorphisms in 2,241 cases (discovery 1,102; replication 1,133) and 776 controls (discovery 405, replication 371) remained. Genotype imputation resulted in 9,077,347 single nucleotide polymorphisms and insertions/deletions with minor allele frequency greater than 0.01 available for analysis. Meta-analysis of the discovery and replication samples identified 6 loci on chromosomes 5, 10, 11, 12 and 18 associated with urgency urinary incontinence at p <10(-6). Of the loci 3 were within genes, the zinc finger protein 521 (ZFP521) gene on chromosome 18q11, the ADAMTS16 gene on chromosome 5p15 and the CIT gene on chromosome 12q24. The other 3 loci were intergenic.

CONCLUSIONS

Although environmental factors also likely contribute, this first exploratory genome-wide association study for urgency urinary incontinence suggests that genetic variants in the ZFP521, CIT and ADAMTS16 genes might account for some of the observed heritability of the condition.

NLRC4 Inflammasome Is an Important Regulator of Interleukin-18 Levels in Patients With Acute Coronary Syndromes

Background—

Interleukin 18 (IL-18) promotes atherosclerotic plaque formation and is increased in patients with acute coronary syndromes. However the relative contribution of genetic variants to the IL-18 levels has not been fully determined.

Methods and Results—

Baseline plasma IL-18 levels were measured in 16 633 patients with acute coronary syndrome, of whom 9340 had genetic data that passed genotype quality control. A 2-stage genome-wide association study was performed, followed by combined analyses using >10 million genotyped or imputed genetic markers. Single nucleotide polymorphisms at 3 loci ( IL18, NLRC4 , and MROH6 ) were identified ( P <3.15×10 −8 ) in the discovery cohort (n=3777) and replicated in the remaining patients (n=5563). In the pooled data (discovery+replication cohort), 7 independent associations, in 5 chromosomal regions, were associated with IL-18 levels (minimum P =6.99×10 –72 ). Six single nucleotide polymorphisms are located in predicted promoter regions of which one disrupts a transcription factor binding site. One single nucleotide polymorphism in NLRC4 is a rare missense variant, predicted to be deleterious to the protein. Altogether, the identified genetic variants explained 8% of the total variation in IL-18 levels in the cohort.

Conclusions—

Our results show that genetic variants play an important role in determining IL-18 levels in patients with acute coronary syndrome and we have identified genetic variants located in the IL-18 gene ( IL18 ) or close to genes that are involved in procaspase-1 activation ( NLRC4 and CARD16 , CARD17 , and CARD18 ). These associations also highlight the importance of the NLRC4 inflammasome for IL-18 production in acute coronary syndrome patients.