Versatile Gene-based Association Study Research Articles

Pathway analysis identified a significant association between cell-cell adherens junctions-related genes and non-syndromic cleft lip/palate in 895 Asian case-parent trios

ObjectivesTo study the association between the seven cell-cell adherens junctions (AJs) -related genes (CDH1, CTNND1, CTNNA1, ESRP1, ESRP2, PLEKHA5, and PLEKHA7) and non-syndromic cleft lip with or without cleft palate (NSCL/P) in the Asian population. DesignThe study included 895 NSCL/P case-parent trios of Asian ethnicity drawn from an international consortium established for a genome-wide association study. We performed association analysis by applying the genotypic transmission disequilibrium test for each of the 144 single nucleotide polymorphisms (SNPs), the versatile gene-based association study for the combined effects of all SNPs, all in or around the seven target genes, and the versatile pathway-based approach for the combined effects of these seven target genes, consecutively. ResultsIn our analysis, we found neither a significant SNP-based nor gene-based association with NSCL/P for any of the 144 relevant SNPs or the seven target genes. However, novel evidence of a significant association (P = 6.00 × 10−6) with NSCL/P was observed in the combined effects of these seven target genes in the versatile pathway-based analysis. ConclusionsOur results were consistent with the findings of previously published human sequencing studies and reinforced the role of these cell-cell AJs-related genes in the pathogenesis of NSCL/P that may be replicated using data from other genome-wide association studies. This pathway approach better reflects our current understanding of the biological mechanisms underlying the aetiology of NSCL/P.

Association of genetic variants in autophagy-lysosome pathway genes with susceptibility and survival to prostate cancer.

Previous studies have indicated the connections between autophagy-lysosome pathway genes dysfunction and prostate cancer, but few studies have investigated whether single nucleotide polymorphisms (SNPs) in autophagy-lysosome pathway genes are implicated in prostate cancer risk and survival. Logistic regression analysis and stepwise Cox regression analysis were conducted in 4,662 cases and 3,114 controls from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial. The false positive rate probability (FPRP) method was applied to correct for multiple comparisons. Gene-based analysis was calculated by versatile gene-based association study approach. We found that SLC11A1 rs7573065 significantly increased the risk of prostate cancer [adjusted odds ratio (OR)=1.24, 95% confidence interval (CI)=1.06-1.46, P=7.02×10-3, FPRP=0.082]. Furthermore, rs7573065 was confirmed as the independent predicator of overall survival (OS) for prostate cancer patients [Hazard ratio (HR)=1.30, 95% CI=1.01-1.66, P=0.041]. The significant association between SLC11A1 and prostate cancer risk was calculated by gene-based analysis (P=0.030). We also observed that the mRNA of SLC11A1 in prostate tumor tissues was significantly over-expressed than that in normal tissues. This study suggested that rs7573065 in SLC11A1 was associated with an increased risk and poor OS of prostate cancer. Our findings may provide evidence for genetic variants in autophagy-lysosome pathway as the risk and prognostic biomarkers for prostate cancer.

Identification of Potential Pleiotropic Genes for Immune and Skeletal Diseases Using Multivariate MetaCCA Analysis.

BackgroundImmune and skeletal systems physiologically and pathologically interact with each other. Immune and skeletal diseases may share potential pleiotropic genetics factors, but the shared specific genes are largely unknown.ObjectiveThis study aimed to investigate the overlapping genetic factors between multiple diseases (including rheumatoid arthritis (RA), psoriasis, osteoporosis, osteoarthritis, sarcopenia, and fracture).MethodsThe canonical correlation analysis (metaCCA) approach was used to identify the shared genes for six diseases by integrating genome-wide association study (GWAS)-derived summary statistics. The versatile Gene-based Association Study (VEGAS2) method was further applied to refine and validate the putative pleiotropic genes identified by metaCCA.ResultsAbout 157 (p<8.19E-6), 319 (p<3.90E-6), and 77 (p<9.72E-6) potential pleiotropic genes were identified shared by two immune diseases, four skeletal diseases, and all of the six diseases, respectively. The top three significant putative pleiotropic genes shared by both immune and skeletal diseases, including HLA-B, TSBP1, and TSBP1-AS1 (p<E-300), were located in the major histocompatibility complex (MHC) region. Nineteen of 77 putative pleiotropic genes identified by metaCCA analysis were associated with at least one disease in the VEGAS2 analysis. Specifically, the majority (18) of these 19 putative validated pleiotropic genes were associated with RA.ConclusionThe metaCCA method identified some pleiotropic genes shared by the immune and skeletal diseases. These findings help to improve our understanding of the shared genetic mechanisms and signaling pathways underlying immune and skeletal diseases.

Exploring the Pleiotropic Genes and Therapeutic Targets Associated with Heart Failure and Chronic Kidney Disease by Integrating metaCCA and SGLT2 Inhibitors' Target Prediction.

Background Previous studies have shown that heart failure (HF) and chronic kidney disease (CKD) have common genetic mechanisms, overlapping pathophysiological pathways, and therapeutic drug—sodium-glucose cotransporter 2 (SGLT2) inhibitors. Methods The genetic pleiotropy metaCCA method was applied on summary statistics data from two independent meta-analyses of GWAS comprising more than 1 million people to identify shared variants and pleiotropic effects between HF and CKD. Targets of SGLT2 inhibitors were predicted by SwissTargetPrediction and DrugBank databases. To refine all genes, we performed using versatile gene-based association study 2 (VEGAS2) and transcriptome-wide association studies (TWAS) for HF and CKD, respectively. Gene enrichment and KEGG pathway analyses were used to explore the potential functional significance of the identified genes and targets. Results After metaCCA analysis, 4,624 SNPs and 1,745 genes were identified to be potentially pleiotropic in the univariate and multivariate SNP-multivariate phenotype analyses, respectively. 21 common genes were detected in both metaCCA and SGLT2 inhibitors' target prediction. In addition, 169 putative pleiotropic genes were identified, which met the significance threshold both in metaCCA analysis and in the VEGAS2 or TWAS analysis for at least one disease. Conclusion We identified novel variants associated with HF and CKD using effectively incorporating information from different GWAS datasets. Our analysis may provide new insights into HF and CKD therapeutic approaches based on the pleiotropic genes, common targets, and mechanisms by integrating the metaCCA method, TWAS and VEGAS2 analyses, and target prediction of SGLT2 inhibitors.

Mapping gene and gene pathways associated with coronary artery disease: a CARDIoGRAM exome and multi-ancestry UK biobank analysis

Coronary artery disease (CAD) genome-wide association studies typically focus on single nucleotide variants (SNVs), and many potentially associated SNVs fail to reach the GWAS significance threshold. We performed gene and pathway-based association (GBA) tests on publicly available Coronary ARtery DIsease Genome wide Replication and Meta-analysis consortium Exome (n = 120,575) and multi ancestry pan UK Biobank study (n = 442,574) summary data using versatile gene-based association study (VEGAS2) and Multi-marker analysis of genomic annotation (MAGMA) to identify novel genes and pathways associated with CAD. We included only exonic SNVs and excluded regulatory regions. VEGAS2 and MAGMA ranked genes and pathways based on aggregated SNV test statistics. We used Bonferroni corrected gene and pathway significance threshold at 3.0 × 10–6 and 1.0 × 10–5, respectively. We also report the top one percent of ranked genes and pathways. We identified 17 top enriched genes with four genes (PCSK9, FAM177, LPL, ARGEF26), reaching statistical significance (p ≤ 3.0 × 10–6) using both GBA tests in two GWAS studies. In addition, our analyses identified ten genes (DUSP13, KCNJ11, CD300LF/RAB37, SLCO1B1, LRRFIP1, QSER1, UBR2, MOB3C, MST1R, and ABCC8) with previously unreported associations with CAD, although none of the single SNV associations within the genes were genome-wide significant. Among the top 1% non-lipid pathways, we detected pathways regulating coagulation, inflammation, neuronal aging, and wound healing.

Gene-based analysis of ADHD using PASCAL: a biological insight into the novel associated genes

BackgroundAttention-Deficit Hyperactivity Disorder (ADHD) is a complex neurodevelopmental disorder (NDD) which may significantly impact on the affected individual’s life. ADHD is acknowledged to have a high heritability component (70–80%). Recently, a meta-analysis of GWAS (Genome Wide Association Studies) has demonstrated the association of several independent loci. Our main aim here, is to apply PASCAL (pathway scoring algorithm), a new gene-based analysis (GBA) method, to the summary statistics obtained in this meta-analysis. PASCAL will take into account the linkage disequilibrium (LD) across genomic regions in a different way than the most commonly employed GBA methods (MAGMA or VEGAS (Versatile Gene-based Association Study)). In addition to PASCAL analysis a gene network and an enrichment analysis for KEGG and GO terms were carried out. Moreover, GENE2FUNC tool was employed to create gene expression heatmaps and to carry out a (DEG) (Differentially Expressed Gene) analysis using GTEX v7 and BrainSpan data.ResultsPASCAL results have revealed the association of new loci with ADHD and it has also highlighted other genes previously reported by MAGMA analysis. PASCAL was able to discover new associations at a gene level for ADHD: FEZF1 (p-value: 2.2 × 10− 7) and FEZF1-AS1 (p-value: 4.58 × 10− 7). In addition, PASCAL has been able to highlight association of other genes that share the same LD block with some previously reported ADHD susceptibility genes. Gene network analysis has revealed several interactors with the associated ADHD genes and different GO and KEGG terms have been associated. In addition, GENE2FUNC has demonstrated the existence of several up and down regulated expression clusters when the associated genes and their interactors were considered.ConclusionsPASCAL has been revealed as an efficient tool to extract additional information from previous GWAS using their summary statistics. This study has identified novel ADHD associated genes that were not previously reported when other GBA methods were employed. Moreover, a biological insight into the biological function of the ADHD associated genes across brain regions and neurodevelopmental stages is provided.

Identification of Shared Genes Between Ischemic Stroke and Parkinson's Disease Using Genome-Wide Association Studies.

Ischemic stroke (IS) and Parkinson's disease (PD) are two neurological diseases that often strike individuals of advanced age. Although thought of as a disease of old age, PD can occur in younger patients. In many of these cases, genetic mutations underlie the disease. As with PD, stroke can also have a genetic component. Although many of the risk factors for IS are considered to be modifiable, a significant portion is not, suggesting that some of stroke risk factors may have a genetic origin. Large-scale genome-wide association studies (GWAS) have identified several IS and PD gene variants recently. Converging epidemiologic and pathological evidence suggests that IS and PD may be linked. However, it is still unclear whether these two conditions share a common mechanism. Here, we sought to determine the genetic mechanism underlying the possible association between IS and PD. We conducted a multi-step systemic analysis comprising (1) identification of IS and PD variants validated by known GWAS, (2) two separate gene-based tests using Versatile Gene-based Association Study 2 (VEGAS2) and PLINK, (3) a transcriptome-wide association study (TWAS), and (4) analyses of gene expression using an online tool in Gene Expression Omnibus. Our investigation revealed that IS and PD have in common five shared genes: GPX7, LBH, ZCCHC10, DENND2A, and NUDT14, which pass gene-based tests. Functionally, these genes are expressed differentially in IS and PD patients compared to neurologically healthy control subjects. This genetic overlap may provide clues on how IS and PD are linked mechanistically. This new genetic insight into these two diseases may be very valuable for narrowing the focus of future studies on the genetic basis of IS and PD and for developing novel therapies.

Associations of NADPH oxidase-related genes with blood pressure changes and incident hypertension: The GenSalt Study.

Previous studies have indicated that reactive oxygen species produced by NADPH oxidase (Nox) are important risk factors of hypertension. The current study aims to examine the associations of Nox related genes with longitudinal blood pressure (BP) changes and the risk of incident hypertension in the Genetic Epidemiology Network of Salt Sensitivity (GenSalt) follow-up study. A total of 1,768 participants from 633 families were included in our analysis. Nine BP measurements were obtained in the morning at baseline and during two follow-up visits. The mixed-effect models were used to investigate the associations of 52 tagged single nucleotide polymorphisms in 11 Nox related genes with BP changes and incident hypertension. Gene-based analyses were performed by truncated product method (TPM) and Versatile Gene-based Association Study (VEGAS). Over the 7.2 years of follow-up, systolic BP (SBP) and diastolic BP (DBP) increased, and 32.1% (512) of participants developed hypertension. SNPs rs12094228, rs16861188 and rs12066019 in NCF2 were significantly associated with longitudinal change in SBP (Pinteraction = 1.1 × 10−3, 2.8 × 10−3 and 1.2 × 10−3, respectively). Gene-based analyses revealed that NCF2 was significantly associated with SBP (PTPM = 1.00 × 10−6, PVEGAS = 1.26 × 10−4) and DBP changes (PTPM = 5.84 × 10−4, PVEGAS = 1.04 × 10−3). These findings suggested that NCF2 may play an important role in BP changes over time in the Han Chinese population.

Neurotransmitter Pathway Genes in Cognitive Decline During Aging: Evidence for GNG4 and KCNQ2 Genes.

Experimental studies support the role of neurotransmitter genes in dementia risk, but human studies utilizing single variants in candidate genes have had limited success. We used the gene-based testing program Versatile Gene-based Association Study to assess whether aggregate variation across 6 neurotransmitter pathways influences risk of cognitive decline in 8159 cognitively normal elderly (≥65 years old) adults from 3 community-based cohorts. Common genetic variation in GNG4 and KCNQ2 was associated with cognitive decline. In human brain tissue data sets, both GNG4 and KCNQ2 show higher expression in hippocampus relative to other brain regions; GNG4 expression decreases with advancing age. Both GNG4 and KCNQ2 show highest expression in fetal astrocytes. Genetic variation analyses and gene expression data suggest that GNG4 and KCNQ2 may be associated with cognitive decline in normal aging. Gene-based testing of neurotransmitter pathways may confirm and reveal novel risk genes in future studies of healthy cognitive aging.

PR Interval Associated Genes, Atrial Remodeling and Rhythm Outcome of Catheter Ablation of Atrial Fibrillation—A Gene-Based Analysis of GWAS Data

Background: PR interval prolongation has recently been shown to associate with advanced left atrial remodeling and atrial fibrillation (AF) recurrence after catheter ablation. While different genome-wide association studies (GWAS) have implicated 13 loci to associate with the PR interval as an AF endophenotype their subsequent associations with AF remodeling and response to catheter ablation are unknown. Here, we perform a gene-based analysis of GWAS data to test the hypothesis that PR interval candidate genes also associate with left atrial remodeling and arrhythmia recurrence following AF catheter ablation.Methods and Results: Samples from 660 patients with paroxysmal (n = 370) or persistent AF (n = 290) undergoing AF catheter ablation were genotyped for ~1,000,000 SNPs. Gene-based association was investigated using VEGAS (versatile gene-based association study). Among the 13 candidate genes, SLC8A1, MEIS1, ITGA9, SCN5A, and SOX5 associated with the PR interval. Of those, ITGA9 and SOX5 were significantly associated with left atrial low voltage areas and left atrial diameter and subsequently with AF recurrence after radiofrequency catheter ablation.Conclusion: This study suggests contributions of ITGA9 and SOX5 to AF remodeling expressed as PR interval prolongation, low voltage areas and left atrial dilatation and subsequently to response to catheter ablation. Future and larger studies are necessary to replicate and apply these findings with the aim of designing AF pathophysiology-based multi-locus risk scores.

FAM13A is a modifier gene of cystic fibrosis lung phenotype regulating rhoa activity, actin cytoskeleton dynamics and epithelial-mesenchymal transition

BackgroundCystic fibrosis (CF) lung disease severity is highly variable and dependent on several factors including genetic modifiers. Family with sequence similarity 13 member A (FAM13A) has been previously associated with lung function in the general population as well as in several chronic lung diseases, such as chronic obstructive pulmonary disease (COPD), we examined whether FAM13A is a modifier gene of CF lung phenotype. We also studied how FAM13A may contribute to the physiopathological mechanisms associated with CF. MethodsWe investigated the association of FAM13A with lung function in CF French patients (n=1222) by SNP-wise analysis and Versatile Gene Based Association Study. We also analyzed the consequences of FAM13A knockdown in A549 cells and primary bronchial epithelial cells from CF patients. ResultsWe found that FAM13A is associated with lung function in CF patients. Utilizing lung epithelial A549 cells and primary human bronchial epithelial cells from CF patients we observed that IL-1β and TGFβ reduced FAM13A expression. Knockdown of FAM13A was associated with increased RhoA activity, induction of F-actin stress fibers and regulation of epithelial-mesenchymal transition markers such as E-cadherin, α-smooth muscle actin and vimentin. ConclusionOur data show that FAM13A is a modifier gene of CF lung phenotype regulating RhoA activity, actin cytoskeleton dynamics and epithelial-mesenchymal transition.

Reporting Correct p Values in VEGAS Analyses

VEGAS (versatile gene-based association study) is a popular methodological framework to perform gene-based tests based on summary statistics from single-variant analyses. The approach incorporates linkage disequilibrium information from reference panels to account for the correlation of test statistics. The gene-based test can utilize three different types of tests. In 2015, the improved framework VEGAS2, using more detailed reference panels, was published. Both versions provide user-friendly web- and offline-based tools for the analysis. However, the implementation of the popular top-percentage test is erroneous in both versions. The p values provided by VEGAS2 are deflated/anti-conservative. Based on real data examples, we demonstrate that this can increase substantially the rate of false-positive findings and can lead to inconsistencies between different test options. We also provide code that allows the user of VEGAS to compute correct p values.

No Evidence That Genetic Variation in the Myeloid-Derived Suppressor Cell Pathway Influences Ovarian Cancer Survival.

Background: The precise mechanism by which the immune system is adversely affected in cancer patients remains poorly understood, but the accumulation of immunosuppressive/protumorigenic myeloid-derived suppressor cells (MDSCs) is thought to be a prominent mechanism contributing to immunologic tolerance of malignant cells in epithelial ovarian cancer (EOC). To this end, we hypothesized genetic variation in MDSC pathway genes would be associated with survival after EOC diagnoses.Methods: We measured the hazard of death due to EOC within 10 years of diagnosis, overall and by invasive subtype, attributable to SNPs in 24 genes relevant in the MDSC pathway in 10,751 women diagnosed with invasive EOC. Versatile Gene-based Association Study and the admixture likelihood method were used to test gene and pathway associations with survival.Results: We did not identify individual SNPs that were significantly associated with survival after correction for multiple testing (P < 3.5 × 10-5), nor did we identify significant associations between the MDSC pathway overall, or the 24 individual genes and EOC survival.Conclusions: In this well-powered analysis, we observed no evidence that inherited variations in MDSC-associated SNPs, individual genes, or the collective genetic pathway contributed to EOC survival outcomes.Impact: Common inherited variation in genes relevant to MDSCs was not associated with survival in women diagnosed with invasive EOC. Cancer Epidemiol Biomarkers Prev; 26(3); 420-4. ©2016 AACR.

Abstract WP204: Genetic Variant in VCAM1 Mediates Acute Infarct Size in Ischemic Stroke Patients

Introduction: Even though neuroinflammation is increasingly recognized as an essential contributor to ischemic brain injury, the exact underlying mechanisms remain unclear. Higher expression of the vascular cell adhesion molecule 1 (VCAM-1) increases leukocyte-brain endothelium interaction and has been associated with larger infarct size following acute ischemic stroke (AIS). The activation of the TGF-β signaling pathways can down-regulate VCAM-1 expression and ameliorate deleterious tissue outcome during neuroinflammation. Hypothesis: We sought to investigate whether genetic variation in the TGF-beta pathway and adhesion molecule genes is associated with acute stroke lesion size. Methods: We completed genome-wide association (GWA) testing and diffusion-weighted imaging lesion volume (DWIv) analysis in a discovery cohort of 532 AIS patients of European ancestry enrolled within 48 hours of symptom onset. An independent European ancestry cohort of 724 AIS patients with GWA data and automated DWIv served as replication cohort. GWA testing per SNP was performed using linear regression modeling of natural log-transformed DWIv adjusted for age, sex and relevant principal components. We selected 42 inflammatory genes in the TGF-beta pathway for a gene-based analysis using VEGAS (Versatile Gene-based -Association Study) software. A pre-specified discovery phase Bonferroni-corrected threshold was set at p&lt;0.001. In the replication phase, 14 SNPs overlapping were tested at the Bonferroni-corrected p-value threshold of p&lt;0.004. Results: Of all genes in the TGF-beta pathway, VCAM1 (p=0.0006) was significantly associated with DWIv in the discovery AIS cohort (age: 66 ± 14.9 years, sex: 63.4% male, DWIv: 2.2cm 3 (IQR: 0.6-11.7cm 3 )). A single SNP within the VCAM1 gene boundaries (rs3176876 (BETA=0.2341, p=0.003)) was significantly associated with DWIv in the replication AIS cohort (age: 65 ± 14.1 years, sex: 68.6% male, DWIv: 4cm 3 (IQR: 1.3-17.2cm 3 )). Conclusion: The genetic variant rs3176876 in the VCAM1 gene is associated with larger infarct lesion size measured on brain MRI of AIS patients. These findings suggest genetic contribution to the pro-inflammatory mechanisms in acute cerebral ischemia and warrant further investigation.

Replication of Candidate SNP Survival Analyses and Gene-Based Tests of Association with Survival Outcomes after an Unrelated Donor Blood or Marrow Transplant: Results from the Discovery-BMT Study

▪Over the past decade, there have been numerous published candidate gene association studies of non-HLA variants, genes and/or pathways in relation to BMT survival outcomes. Most of these studies investigated a small number of candidate single nucleotide polymorphisms (SNPs), chosen based on known gene function, in a limited number of cases without replication. Here we test the association of previously reported SNPs with survival outcomes (overall survival (OS), progression-free survival (PFS), death due to disease (DD) or transplant-related mortality (TRM)) after URD BMT using data from the DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to one-Year mortality after BMT) study which included 2,052 European descent recipient-unrelated donor pairs (Cohort 1) and 763 recipient-unrelated donor pairs (Cohort 1) reported to the Center for International Blood and Marrow Transplant Research (CIBMTR) between 2000-2009 and 2009-2011, respectively. DISCOVeRY-BMT is a Genome Wide Association Study (GWAS) that includes over 9 million typed or imputed SNPs. This represents the largest study to date and is well powered to test both SNP and gene-based association hypotheses.We conducted an extensive literature review of PubMed using MESH terms related to transplant and hematologic disease which yielded 76 candidate gene association studies that investigated a total of 46 genes with any survival outcome. Of the 46 genes, 11 genes (CTLA4, CYP2B6, CYP2C19, GSTM1, GSTT1, IL23R, IL6, MTHFR, NOD2/CARD15, TNFSRF1B, and VDR) were studied in at least 2 publications. A total of 16 SNPs in 8 of these genes reported a significant (P< 0.05) association with at least one survival outcome (OS, PFS, DD or TRM). We tested these 16 SNPs for an association with the survival outcome originally reported in each article (Table 1). All association models included age at BMT, diagnosis (AML, ALL, MDS), disease status at BMT, cell source (blood, bone marrow) and year of BMT. P-values from Cohorts 1 and 2 were combined using a fixed-effects model in METAL software. Since the candidate gene studies have a gene-based hypothesis, we also performed gene-level testing using the VErsatile Gene-Based Association Study (VEGAS2) software to measure the aggregate effect of all available typed and imputed SNPs from our GWAS data in these 8 genes.The candidate gene association studies included heterogeneous patient populations of autologous and/or allogeneic (related or unrelated donors) BMT for a variety of hematologic diseases. Most candidate gene studies performed genotyping on both donors and recipients, and SNPs were tested either independently for donor or recipient genomes, or considering the total number of variant alleles the recipient-donor pair had at a given locus (0-4) (Table 1).The 16 SNPs from the previous candidate gene association studies were not significantly associated at a nominal P<0.05 to either the survival outcome reported in the original published paper, or to any survival outcome (Table 1). Likewise, the VEGAS2 gene-based analysis did not reveal any statistically significant associations for any of the 8 genes at P<0.05. Our study also corroborated prior reports of SNPs and genes that were not associated with survival outcomes.Candidate gene association studies use a priori biological knowledge to select genes that are likely associated with outcome and subsequently SNPs within the gene thought likely to affect function based on knowledge at the time of the investigation. Our findings suggest that neither the genetic variation tested in these candidate genes, nor other individual SNPs in the gene independently or collectively are not associated with survival outcomes in recipients of URD BMT. Our failure to replicate the prior candidate gene studies may be due to these SNPs being dependent on the underlying disease (ie. AML, ALL, etc), the donor relation (related, unrelated, autologous), or that there is no true association of these SNPs with the outcomes studied.#These authors contributed equally [Display omitted] DisclosuresMcCarthy:Gamida Cell: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Karyopharm: Honoraria, Membership on an entity's Board of Directors or advisory committees; Onyx: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bristol Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; The Binding Site: Honoraria, Membership on an entity's Board of Directors or advisory committees; Sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees. Hahn:Novartis: Equity Ownership; NIH: Research Funding. Sucheston-Campbell:NIH/NCI: Research Funding.

Genome-wide association study using family-based cohorts identifies the WLS and CCDC170/ESR1 loci as associated with bone mineral density.

BackgroundOsteoporosis is a common and debilitating bone disease that is characterised by a low bone mineral density (BMD), a highly heritable trait. Genome-wide association studies (GWAS) have proven to be very successful in identifying common genetic variants associated with BMD adjusted for age, gender and weight, however a large portion of the genetic variance for this trait remains unexplained. There is evidence to suggest significant genetic correlation between body size traits and BMD. It has also recently been suggested that unintended bias can be introduced as a result of adjusting a phenotype for a correlated trait. We performed a GWAS meta-analysis in two populations (total n = 6,696) using BMD data adjusted for only age and gender, in an attempt to identify genetic variants associated with BMD including those that may have potential pleiotropic effects on BMD and body size traits.ResultsWe observed a single variant, rs2566752, associated with spine BMD at the genome-wide significance level in the meta-analysis (P = 3.36 × 10−09). Logistic regression analysis also revealed an association between rs2566752 and fracture rate in one of our study cohorts (P = 0.017, n = 5,654). This is an intronic variant located in the wntless Wnt ligand secretion mediator (WLS) gene (1p31.3), a known BMD locus which encodes an integral component of the Wnt ligand secretion pathway. Bioinformatics analyses of variants in moderate LD with rs2566752 produced strong evidence for a regulatory role for the variants rs72670452, rs17130567 and rs1430738. Expression quantitative trait locus (eQTL) analysis suggested that the variants rs12568456 and rs17130567 are associated with expression of the WLS gene in whole blood, cerebellum and temporal cortex brain tissue (P = 0.034–1.19 × 10−23). Gene-wide association testing using the VErsatile Gene-based Association Study 2 (VEGAS2) software revealed associations between the coiled-coil domain containing 170 (CCDC170) gene, located adjacent to the oestrogen receptor 1 (ESR1) gene, and BMD at the spine, femoral neck and total hip sites (P = 1.0 × 10−06, 2.0 × 10−06 and 2.0 × 10−06 respectively).ConclusionsGenetic variation at the WLS and CCDC170/ESR1 loci were found to be significantly associated with BMD adjusted for only age and gender at the genome-wide level in this meta-analysis.

Polymorphisms in genes in the androgen pathway and risk of Barrett's esophagus and esophageal adenocarcinoma.

The strong male predominance in Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) remains inadequately explained, but sex hormones might be involved. We hypothesized that single nucleotide polymorphisms (SNPs) in the androgen pathway influence risk of developing BE and EAC. This genetic-epidemiological analysis included 14 studies from Australia, Europe and North America. Polymorphisms in 16 genes coding for the androgen pathway were analyzed using a gene-based approach: versatile gene-based test association study. This method evaluates associations between a trait and all SNPs within a specific gene rather than each SNP marker individually as in a conventional GWAS. The data were stratified for sex, body-mass index, waist-to-hip ratio, tobacco smoking and gastroesophageal reflux status. Included were data from 1,508 EAC patients, 2,383 BE patients and 2,170 control participants. SNPs within the gene CYP17A1 were associated with risk of BE in the sexes combined (p = 0.002) and in males (p = 0.003), but not in females separately (p = 0.3). This association was found in tobacco smokers (p = 0.003) and in BE patients without reflux (p = 0.004), but not in nonsmokers (p = 0.2) or those with reflux (p = 0.036). SNPs within JMJD1C were associated with risk of EAC in females (p = 0.001). However, none of these associations replicated in a subsequent sample. Fourteen other genes studied did not reach statistically significant levels of association with BE, EAC or the combination of BE and EAC, after correcting for the number of genes included in the analysis. In conclusion, genetic variants in the androgen-related genes CYP17A1 and JMJD1C might be associated with risk of BE and EAC, respectively, but replication data with larger sample sizes are needed.

Polymorphisms in Genes of Relevance for Oestrogen and Oxytocin Pathways and Risk of Barrett's Oesophagus and Oesophageal Adenocarcinoma: A Pooled Analysis from the BEACON Consortium.

BackgroundThe strong male predominance in oesophageal adenocarcinoma (OAC) and Barrett’s oesophagus (BO) continues to puzzle. Hormonal influence, e.g. oestrogen or oxytocin, might contribute.MethodsThis genetic-epidemiological study pooled 14 studies from three continents, Australia, Europe, and North America. Polymorphisms in 3 key genes coding for the oestrogen pathway (receptor alpha (ESR1), receptor beta (ESR2), and aromatase (CYP19A1)), and 3 key genes of the oxytocin pathway (the oxytocin receptor (OXTR), oxytocin protein (OXT), and cyclic ADP ribose hydrolase glycoprotein (CD38)), were analysed using a gene-based approach, versatile gene-based test association study (VEGAS).ResultsAmong 1508 OAC patients, 2383 BO patients, and 2170 controls, genetic variants within ESR1 were associated with BO in males (p = 0.0058) and an increased risk of OAC and BO combined in males (p = 0.0023). Genetic variants within OXTR were associated with an increased risk of BO in both sexes combined (p = 0.0035) and in males (p = 0.0012). We followed up these suggestive findings in a further smaller data set, but found no replication. There were no significant associations between the other 4 genes studied and risk of OAC, BO, separately on in combination, in males and females combined or in males only.ConclusionGenetic variants in the oestrogen receptor alpha and the oxytocin receptor may be associated with an increased risk of BO or OAC, but replication in other large samples are needed.

Gene-based meta-analysis of genome-wide association studies implicates new loci involved in obesity.

To date, genome-wide association studies (GWASs) have identified >100 loci with single variants associated with body mass index (BMI). This approach may miss loci with high allelic heterogeneity; therefore, the aim of the present study was to use gene-based meta-analysis to identify regions with high allelic heterogeneity to discover additional obesity susceptibility loci. We included GWAS data from 123 865 individuals of European descent from 46 cohorts in Stage 1 and Metabochip data from additional 103 046 individuals from 43 cohorts in Stage 2, all within the Genetic Investigation of ANthropometric Traits (GIANT) consortium. Each cohort was tested for association between ∼2.4 million (Stage 1) or ∼200 000 (Stage 2) imputed or genotyped single variants and BMI, and summary statistics were subsequently meta-analyzed in 17 941 genes. We used the 'VErsatile Gene-based Association Study' (VEGAS) approach to assign variants to genes and to calculate gene-based P-values based on simulations. The VEGAS method was applied to each cohort separately before a gene-based meta-analysis was performed. In Stage 1, two known (FTO and TMEM18) and six novel (PEX2, MTFR2, SSFA2, IARS2, CEP295 and TXNDC12) loci were associated with BMI (P < 2.8 × 10(-6) for 17 941 gene tests). We confirmed all loci, and six of them were gene-wide significant in Stage 2 alone. We provide biological support for the loci by pathway, expression and methylation analyses. Our results indicate that gene-based meta-analysis of GWAS provides a useful strategy to find loci of interest that were not identified in standard single-marker analyses due to high allelic heterogeneity.

Beyond the usual suspects: a multidimensional genetic exploration of infant attachment disorganization and security

Although the environmental influences on infant attachment disorganization and security are well-studied, little is known about their heritability. Candidate gene studies have shown small, often non-replicable effects. In this study, we gathered the largest sample (N = 657) of ethnically homogenous, 14-month-old children with both observed attachment and genome-wide data. First, we used a Genome-Wide Association Study (GWAS) approach to identify single nucleotide polymorphisms (SNPs) associated with attachment disorganization and security. Second, we annotated them into genes (Versatile Gene-based Association Study) and functional pathways. Our analyses provide evidence of novel genes (HDAC1, ZNF675, BSCD1) and pathways (synaptic transmission, cation transport) associated with attachment disorganization. Similar analyses identified a novel gene (BECN1) but no distinct pathways associated with attachment security. The results of this first extensive, exploratory study on the molecular-genetic basis of infant attachment await replication in large, independent samples.