Discovery properties of genome-wide association signals from cumulatively combined data sets

Meta-regression of genome-wide association studies to estimate age-varying genetic effects

Fixed-effect meta-analysis has been used to summarize genetic effects on a phenotype across multiple Genome-Wide Association Studies (GWAS) assuming a common underlying genetic effect. Genetic effects may vary with age (or other characteristics), and not allowing for this in a GWAS might lead to bias. Meta-regression models between study heterogeneity and allows effect modification of the genetic effects to be explored. The aim of this study was to explore the use of meta-analysis and meta-regression for estimating age-varying genetic effects on phenotypes. With simulations we compared the performance of meta-regression to fixed-effect and random -effects meta-analyses in estimating (i) main genetic effects and (ii) age-varying genetic effects (SNP by age interactions) from multiple GWAS studies under a range of scenarios. We applied meta-regression on publicly available summary data to estimate the main and age-varying genetic effects of the FTO SNP rs9939609 on Body Mass Index (BMI). Fixed-effect and random-effects meta-analyses accurately estimated genetic effects when these did not change with age. Meta-regression accurately estimated both main genetic effects and age-varying genetic effects. When the number of studies or the age-diversity between studies was low, meta-regression had limited power. In the applied example, each additional minor allele (A) of rs9939609 was inversely associated with BMI at ages 0 to 3, and positively associated at ages 5.5 to 13. Our findings challenge the assumption that genetic effects are consistent across all ages and provide a method for exploring this. GWAS consortia should be encouraged to use meta-regression to explore age-varying genetic effects.

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at ∼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.

Discerning asthma endotypes through comorbidity mapping

Asthma is a heterogeneous, complex syndrome, and identifying asthma endotypes has been challenging. We hypothesize that distinct endotypes of asthma arise in disparate genetic variation and life-time environmental exposure backgrounds, and that disease comorbidity patterns serve as a surrogate for such genetic and exposure variations. Here, we computationally discover 22 distinct comorbid disease patterns among individuals with asthma (asthma comorbidity subgroups) using diagnosis records for >151 M US residents, and re-identify 11 of the 22 subgroups in the much smaller UK Biobank. GWASs to discern asthma risk loci for individuals within each subgroup and in all subgroups combined reveal 109 independent risk loci, of which 52 are replicated in multi-ancestry meta-analysis across different ethnicity subsamples in UK Biobank, US BioVU, and BioBank Japan. Fourteen loci confer asthma risk in multiple subgroups and in all subgroups combined. Importantly, another six loci confer asthma risk in only one subgroup. The strength of association between asthma and each of 44 health-related phenotypes also varies dramatically across subgroups. This work reveals subpopulations of asthma patients distinguished by comorbidity patterns, asthma risk loci, gene expression, and health-related phenotypes, and so reveals different asthma endotypes.

Data Integration, Imputation, and Meta-analysis for Genome-Wide Association Studies

Growing genomic and phenotypic datasets require different groups around the world to collaborate and integrate these valuable resources to maximize their benefit and increase reference population sizes for genomic prediction and genome-wide association studies (GWAS). However, different studies use different genotyping techniques which requires a synchronizing step for the genotyped variants called "imputation" before combining them. Optimally, different GWAS datasets can be analysed within a meta-analysis, which recruits summary statistics instead of actual data. This chapter describes the general principles for genotypic imputation and meta-GWAS analysis with a description of study designs and command lines required for such analyses.

Meta-analysis of genome-wide association studies for litter size in sheep

Litter size and ovulation rate are important reproduction traits in sheep and have important impacts on the profitability of farm animals. To investigate the genetic architecture of litter size, we report the first meta-analysis of genome-wide association studies (GWAS) using 522 ewes and 564,377 SNPs from six sheep breeds. We identified 29 significant associations for litter size which 27 of which have not been reported in individual GWAS for each population. However, we could confirm the role of BMPR1B in prolificacy. Our gene set analysis discovered biological pathways related to cell signaling, communication, and adhesion. Functional clustering and enrichment using protein databases identified epidermal growth factor-like domain affecting litter size. Through analyzing protein-protein interaction data, we could identify hub genes like CASK, PLCB4, RPTOR, GRIA2, and PLCB1 that were enriched in most of the significant pathways. These genes have a role in cell proliferation, cell adhesion, cell growth and survival, and autophagy. Notably, identified SNPs were scattered on several different chromosomes implying different genetic mechanisms underlying variation of prolificacy in each breed. Given the different layers that make up the follicles and the need for communication and transfer of hormones and nutrients through these layers to the oocyte, the significance of pathways related to cell signaling and communication seems logical. Our results provide genetic insights into the litter size variation in different sheep breeds.

Reviewing and assessing existing meta-analysis models and tools

Over the past few years, meta-analysis has become popular among biomedical researchers for detecting biomarkers across multiple cohort studies with increased predictive power. Combining datasets from different sources increases sample size, thus overcoming the issue related to limited sample size from each individual study and boosting the predictive power. This leads to an increased likelihood of more accurately predicting differentially expressed genes/proteins or significant biomarkers underlying the biological condition of interest. Currently, several meta-analysis methods and tools exist, each having its own strengths and limitations. In this paper, we survey existing meta-analysis methods, and assess the performance of different methods based on results from different datasets as well as assessment from prior knowledge of each method. This provides a reference summary of meta-analysis models and tools, which helps to guide end-users on the choice of appropriate models or tools for given types of datasets and enables developers to consider current advances when planning the development of new meta-analysis models and more practical integrative tools.

Meta-Analysis of SNP-Environment Interaction With Overlapping Data

Meta-analysis, which combines the results of multiple studies, is an important analytical method in genome-wide association studies. In genome-wide association studies practice, studies employing meta-analysis may have overlapping data, which could yield false positive results. Recent studies have proposed models to handle the issue of overlapping data when testing the genetic main effect of single nucleotide polymorphism. However, there is still no meta-analysis method for testing gene-environment interaction when overlapping data exist. Inspired by the methods of testing the main effect of gene with overlapping data, we proposed an overlapping meta-regulation method to address the issue in testing the gene-environment interaction. We generalized the covariance matrices of the regular meta-regression model by employing Lin’s and Han’s correlation structures to incorporate the correlations introduced by the overlapping data. Based on our proposed models, we further provided statistical significance tests of the gene-environment interaction as well as joint effects of the gene main effect and the interaction. Through simulations, we examined type I errors and statistical powers of our proposed methods at different levels of data overlap among studies. We demonstrated that our method well controls the type I error and simultaneously achieves statistical power comparable with the method that removes overlapping samples a priori before the meta-analysis, i.e., the splitting method. On the other hand, ignoring overlapping data will inflate the type I error. Unlike the splitting method that requires individual-level genotype and phenotype data, our proposed method for testing gene-environment interaction handles the issue of overlapping data effectively and statistically efficiently at the meta-analysis level.

Statistical Methods and Software for Substance Use and Dependence Genetic Research

BACKGROUND

Substantial substance use disorders and related health conditions emerged dur-ing the mid-20th century and continue to represent a remarkable 21st century global burden of disease. This burden is largely driven by the substance-dependence process, which is a complex process and is influenced by both genetic and environmental factors. During the past few decades, a great deal of pro-gress has been made in identifying genetic variants associated with Substance Use and Dependence (SUD) through linkage, candidate gene association, genome-wide association and sequencing studies.

METHODS

Various statistical methods and software have been employed in different types of SUD ge-netic studies, facilitating the identification of new SUD-related variants.

CONCLUSION

In this article, we review statistical methods and software that are currently available for SUD genetic studies, and discuss their strengths and limitations.

Genetic determinants of risk in pulmonary arterial hypertension: international genome-wide association studies and meta-analysis

BACKGROUND

Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes.

METHODS

We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11 744 individuals with European ancestry (including 2085 patients). One GWAS used genotypes from 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals. Cross-validation of loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration of survival. All-cause mortality was the primary endpoint in survival analyses.

FINDINGS

A locus near SOX17 (rs10103692, odds ratio 1·80 [95% CI 1·55-2·08], p=5·13 × 10^-15) and a second locus in HLA-DPA1 and HLA-DPB1 (collectively referred to as HLA-DPA1/DPB1 here; rs2856830, 1·56 [1·42-1·71], p=7·65 × 10^-20) within the class II MHC region were associated with pulmonary arterial hypertension. The SOX17 locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25-1·48], p=1·69 × 10^-12; and rs10103692). Functional and epigenomic data indicate that the risk variants near SOX17 alter gene regulation via an enhancer active in endothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reduced SOX17 expression. The HLA-DPA1/DPB1 rs2856830 genotype was strongly associated with survival. Median survival from diagnosis in patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02-8·05]), despite similar baseline disease severity.

INTERPRETATION

This is the first study to report that common genetic variation at loci in an enhancer near SOX17 and in HLA-DPA1/DPB1 is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by rare mutations in SOX17. Further studies are needed to confirm the association between HLA typing or rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improves risk stratification in clinical practice or trials.

FUNDING

UK NIHR, BHF, UK MRC, Dinosaur Trust, NIH/NHLBI, ERS, EMBO, Wellcome Trust, EU, AHA, ACClinPharm, Netherlands CVRI, Dutch Heart Foundation, Dutch Federation of UMC, Netherlands OHRD and RNAS, German DFG, German BMBF, APH Paris, INSERM, Université Paris-Sud, and French ANR.

Comparison of Two Meta-Analysis Methods: Inverse-Variance-Weighted Average and Weighted Sum of Z-Scores

The meta-analysis has become a widely used tool for many applications in bioinformatics, including genome-wide association studies. A commonly used approach for meta-analysis is the fixed effects model approach, for which there are two popular methods: the inverse variance-weighted average method and weighted sum of z-scores method. Although previous studies have shown that the two methods perform similarly, their characteristics and their relationship have not been thoroughly investigated. In this paper, we investigate the optimal characteristics of the two methods and show the connection between the two methods. We demonstrate that the each method is optimized for a unique goal, which gives us insight into the optimal weights for the weighted sum of z-scores method. We examine the connection between the two methods both analytically and empirically and show that their resulting statistics become equivalent under certain assumptions. Finally, we apply both methods to the Wellcome Trust Case Control Consortium data and demonstrate that the two methods can give distinct results in certain study designs.

Power considerations for λ inflation factor in meta-analyses of genome-wide association studies

The genomic control (GC) approach is extensively used to effectively control false positive signals due to population stratification in genome-wide association studies (GWAS). However, GC affects the statistical power of GWAS. The loss of power depends on the magnitude of the inflation factor (λ) that is used for GC. We simulated meta-analyses of different GWAS. Minor allele frequency (MAF) ranged from 0·001 to 0·5 and λ was sampled from two scenarios: (i) random scenario (empirically-derived distribution of real λ values) and (ii) selected scenario from simulation parameter modification. Adjustment for λ was considered under single correction (within study corrected standard errors) and double correction (additional λ corrected summary estimate). MAF was a pivotal determinant of observed power. In random λ scenario, double correction induced a symmetric power reduction in comparison to single correction. For MAF 1·2 and MAF >5%. Our results provide a quick but detailed index for power considerations of future meta-analyses of GWAS that enables a more flexible design from early steps based on the number of studies accumulated in different groups and the λ values observed in the single studies.

Recurrent Coding Sequence Variation Explains Only A Small Fraction of the Genetic Architecture of Colorectal Cancer

Whilst common genetic variation in many non-coding genomic regulatory regions are known to impart risk of colorectal cancer (CRC), much of the heritability of CRC remains unexplained. To examine the role of recurrent coding sequence variation in CRC aetiology, we genotyped 12,638 CRCs cases and 29,045 controls from six European populations. Single-variant analysis identified a coding variant (rs3184504) in SH2B3 (12q24) associated with CRC risk (OR = 1.08, P = 3.9 × 10(-7)), and novel damaging coding variants in 3 genes previously tagged by GWAS efforts; rs16888728 (8q24) in UTP23 (OR = 1.15, P = 1.4 × 10(-7)); rs6580742 and rs12303082 (12q13) in FAM186A (OR = 1.11, P = 1.2 × 10(-7) and OR = 1.09, P = 7.4 × 10(-8)); rs1129406 (12q13) in ATF1 (OR = 1.11, P = 8.3 × 10(-9)), all reaching exome-wide significance levels. Gene based tests identified associations between CRC and PCDHGA genes (P < 2.90 × 10(-6)). We found an excess of rare, damaging variants in base-excision (P = 2.4 × 10(-4)) and DNA mismatch repair genes (P = 6.1 × 10(-4)) consistent with a recessive mode of inheritance. This study comprehensively explores the contribution of coding sequence variation to CRC risk, identifying associations with coding variation in 4 genes and PCDHG gene cluster and several candidate recessive alleles. However, these findings suggest that recurrent, low-frequency coding variants account for a minority of the unexplained heritability of CRC.

Replication of ZNF804A gene variant associations with risk of heroin addiction

Heroin addiction is heritable, but few specific genetic variants have been reproducibly associated with this disease. The zinc finger protein 804A (ZNF804A) gene is a biologically plausible susceptibility gene for heroin addiction, given its function as a transcription factor in human brain. Novel associations of two common ZNF804A single nucleotide polymorphisms (SNPs), rs7597593 and rs1344706, with heroin addiction have been reported in Han Chinese. Both SNPs have also been implicated for regulating ZNF804A expression in human brain, including the addiction-relevant dorsolateral prefrontal cortex. In this independent replication study, we tested the rs7597593 and rs1344706 SNP genotypes and their corresponding haplotypes for association with heroin addiction using cases drawn from the Urban Health Study and population controls: total N = 10 757 [7095 European Americans (EAs) and 3662 African Americans (AAs)]. We independently replicated both ZNF804A SNP associations in EAs: the rs7597593-T (P = 0.016) and rs1344706-A (P = 0.029) alleles both being associated with increased risk of heroin addiction, consistent with the prior report. Neither SNP was associated in AAs alone, but meta-analysis across both ancestry groups resulted in significant associations for rs1344706-A [P = 0.016, odds ratio (95% confidence interval) = 1.13 (1.02-1.25)] and its haplotype with rs7597593-T [P = 0.0067, odds ratio (95% confidence interval) = 1.16 (1.04-1.29)]. By showing consistent associations across independent studies and diverse ancestry groups, our study provides evidence that these two ZNF804A SNPs and their risk haplotype are among the few replicable genetic associations with heroin addiction.

Statistical and Computational Methods for Genetic Diseases: An Overview

The identification of causes of genetic diseases has been carried out by several approaches with increasing complexity. Innovation of genetic methodologies leads to the production of large amounts of data that needs the support of statistical and computational methods to be correctly processed. The aim of the paper is to provide an overview of statistical and computational methods paying attention to methods for the sequence analysis and complex diseases.

Multiethnic genome-wide association study of cerebral white matter hyperintensities on MRI

BACKGROUND

The burden of cerebral white matter hyperintensities (WMH) is associated with an increased risk of stroke, dementia, and death. WMH are highly heritable, but their genetic underpinnings are incompletely characterized. To identify novel genetic variants influencing WMH burden, we conducted a meta-analysis of multiethnic genome-wide association studies.

METHODS AND RESULTS

We included 21 079 middle-aged to elderly individuals from 29 population-based cohorts, who were free of dementia and stroke and were of European (n=17 936), African (n=1943), Hispanic (n=795), and Asian (n=405) descent. WMH burden was quantified on MRI either by a validated automated segmentation method or a validated visual grading scale. Genotype data in each study were imputed to the 1000 Genomes reference. Within each ethnic group, we investigated the relationship between each single-nucleotide polymorphism and WMH burden using a linear regression model adjusted for age, sex, intracranial volume, and principal components of ancestry. A meta-analysis was conducted for each ethnicity separately and for the combined sample. In the European descent samples, we confirmed a previously known locus on chr17q25 (P=2.7×10(-19)) and identified novel loci on chr10q24 (P=1.6×10(-9)) and chr2p21 (P=4.4×10(-8)). In the multiethnic meta-analysis, we identified 2 additional loci, on chr1q22 (P=2.0×10(-8)) and chr2p16 (P=1.5×10(-8)). The novel loci contained genes that have been implicated in Alzheimer disease (chr2p21 and chr10q24), intracerebral hemorrhage (chr1q22), neuroinflammatory diseases (chr2p21), and glioma (chr10q24 and chr2p16).

CONCLUSIONS

We identified 4 novel genetic loci that implicate inflammatory and glial proliferative pathways in the development of WMH in addition to previously proposed ischemic mechanisms.

An empirical comparison of meta-analysis and mega-analysis of individual participant data for identifying gene-environment interactions

For analysis of the main effects of SNPs, meta-analysis of summary results from individual studies has been shown to provide comparable results as "mega-analysis" that jointly analyzes the pooled participant data from the available studies. This fact revolutionized the genetic analysis of complex traits through large GWAS consortia. Investigations of gene-environment (G×E) interactions are on the rise since they can potentially explain a part of the missing heritability and identify individuals at high risk for disease. However, for analysis of gene-environment interactions, it is not known whether these methods yield comparable results. In this empirical study, we report that the results from both methods were largely consistent for all four tests; the standard 1 degree of freedom (df) test of main effect only, the 1 df test of the main effect (in the presence of interaction effect), the 1 df test of the interaction effect, and the joint 2 df test of main and interaction effects. They provided similar effect size and standard error estimates, leading to comparable P-values. The genomic inflation factors and the number of SNPs with various thresholds were also comparable between the two approaches. Mega-analysis is not always feasible especially in very large and diverse consortia since pooling of raw data may be limited by the terms of the informed consent. Our study illustrates that meta-analysis can be an effective approach also for identifying interactions. To our knowledge, this is the first report investigating meta-versus mega-analyses for interactions.

On individual genome-wide association studies and their meta-analysis

Individual genome-wide association (GWA) studies and their meta-analyses represent two approaches for identifying genetic loci associated with complex diseases/traits. Inconsistent findings and non-replicability between individual GWA studies and meta-analyses are commonly observed, hence posing the critical question as to how to interpret their respective results properly. In this study, we performed a series of simulation studies to investigate and compare the statistical properties of the two approaches. Our results show that (1) as expected, meta-analysis of larger sample size is more powerful than individual GWA studies under the ideal setting of population homogeneity among individual studies; (2) under the realistic setting of heterogeneity among individual studies, detection of heterogeneity is usually difficult and meta-analysis (even with the random-effects model) may introduce elevated false positive and/or negative rates; (3) despite relatively small sample size, well-designed individual GWA study has the capacity to identify novel loci for complex traits; (4) replicability between meta-analysis and independent individual studies or between independent meta-analyses is limited, and thus inconsistent findings are not unexpected.

Multivariate meta-analysis of the association of G-protein beta 3 gene (GNB3) haplotypes with cardiovascular phenotypes

The objective of the present study was to review previous investigations on the association of haplotypes in the G-protein β3 subunit (GNB3) gene with representative cardiovascular risk factors/phenotypes: hypertension, overweight, and variation in the systolic and diastolic blood pressures (SBP and DBP, respectively) and as well as body mass index (BMI). A comprehensive literature search was undertaken in Pubmed, Web of Science, EMBASE, Biological Abstracts, LILACS and Google Scholar to identify potentially relevant articles published up to April 2011. Six genetic association studies encompassing 16,068 participants were identified. Individual participant data were obtained for all studies. The three most investigated GNB3 polymorphisms (G-350A, C825T and C1429T) were considered. Expectation-maximization and generalized linear models were employed to estimate haplotypic effects from data with uncertain phase while adjusting for covariates. Study-specific results were combined through a random-effects multivariate meta-analysis. After carefully adjustments for relevant confounding factors, our analysis failed to support a role for GNB3 haplotypes in any of the investigated phenotypes. Sensitivity analyses excluding studies violating Hardy-Weinberg expectations, considering gender-specific effects or more extreme phenotypes (e.g. obesity only) as well as a fixed-effects "pooled" analysis also did not disclose a significant influence of GNB3 haplotypes on cardiovascular phenotypes. We conclude that the previous cumulative evidence does not support the proposal that haplotypes formed by common GNB3 polymorphisms might contribute either to the development of hypertension and obesity, or to the variation in the SBP, DBP and BMI.

Clinical implications of family history of prostate cancer and genetic risk single nucleotide polymorphism (SNP) profiles in an active surveillance cohort

OBJECTIVES

To explore the potential prognostic role of family history (FH) of prostate cancer and prostate cancer risk single nucleotide polymorphisms (SNPs) in patients undergoing active surveillance (AS) for prostate cancer. This is the first study to date, which has investigated the potential prognostic role of SNP profiles in an AS cohort

PATIENTS AND METHODS

FH data were collected from patients in the Royal Marsden Hospital AS study. In all, 39 prostate cancer-risk SNPs identified from published genome wide association studies (GWAS) were genotyped using the Sequenom Platform and TaqMan™ assays from available DNA. The cumulative genetic-risk scores for each patient were then calculated using the weighted effect estimated from previous GWAS (log-additive model). FH status and the genetic-risk scores were assessed against adverse outcomes in AS, time to treatment and adverse histology on repeat biopsy, using univariable and multivariable Cox regression models to address time to treatment; and binary logistic regression to address biopsy upgrade.

RESULTS

Of 471 patients, 55 (13.6%) had adverse histology on repeat biopsies and 145 (30.8%) had deferred treatment. On univariate analysis, there was no significant relationship between FH of prostate cancer in any degree of relation, and adverse histology or time to treatment. For risk score analyses, 386 patients' DNA was studied; and there was also no relationship found between the calculated genetic risk scores and adverse histology or time to treatment (P = 0.573 and P = 0.965, respectively). The retrospective study design and the few events were the main limitation of the study.

CONCLUSIONS

There is currently insufficient data to support the use of FH status or prostate cancer SNP profile risk scores as prognostic factors in AS and these should not be used to influence management decisions. As more genetic variants are discovered this may change and should be reassessed in multicentre AS cohorts.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm

BACKGROUND

Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown.

AIM

To identify genetic variants associated with forearm BMD and forearm fractures.

METHODS

BMD at distal radius, measured by dual-energy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a meta-analysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls.

RESULTS

We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10(-8)) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10(-9)). The gene-based association test suggested an association between MEF2C and forearm BMD (p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2.

CONCLUSIONS

These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

The power of meta-analysis in genome-wide association studies

Meta-analysis of multiple genome-wide association (GWA) studies has become common practice over the past few years. The main advantage of this technique is the maximization of power to detect subtle genetic effects for common traits. Moreover, one can use meta-analysis to probe and identify heterogeneity in the effect sizes across the combined studies. In this review, we systematically appraise and evaluate the characteristics of GWA meta-analyses with 10,000 or more subjects published up to June 2012. We provide an overview of the current landscape of variants discovered by GWA meta-analyses, and we discuss and assess with extrapolations from empirical data the value of larger meta-analyses for the discovery of additional genetic associations and new biology in the future. Finally, we discuss some emerging logistical and practical issues related to the conduct of meta-analysis of GWA studies.

Meta-analysis methods for genome-wide association studies and beyond

Meta-analysis of genome-wide association studies (GWASs) has become a popular method for discovering genetic risk variants. Here, we overview both widely applied and newer statistical methods for GWAS meta-analysis, including issues of interpretation and assessment of sources of heterogeneity. We also discuss extensions of these meta-analysis methods to complex data. Where possible, we provide guidelines for researchers who are planning to use these methods. Furthermore, we address special issues that may arise for meta-analysis of sequencing data and rare variants. Finally, we discuss challenges and solutions surrounding the goals of making meta-analysis data publicly available and building powerful consortia.

Genome-wide association study meta-analysis of European and Asian-ancestry samples identifies three novel loci associated with bipolar disorder

Meta-analyses of bipolar disorder (BD) genome-wide association studies (GWAS) have identified several genome-wide significant signals in European-ancestry samples, but so far account for little of the inherited risk. We performed a meta-analysis of ∼750,000 high-quality genetic markers on a combined sample of ∼14,000 subjects of European and Asian-ancestry (phase I). The most significant findings were further tested in an extended sample of ∼17,700 cases and controls (phase II). The results suggest novel association findings near the genes TRANK1 (LBA1), LMAN2L and PTGFR. In phase I, the most significant single nucleotide polymorphism (SNP), rs9834970 near TRANK1, was significant at the P=2.4 × 10(-11) level, with no heterogeneity. Supportive evidence for prior association findings near ANK3 and a locus on chromosome 3p21.1 was also observed. The phase II results were similar, although the heterogeneity test became significant for several SNPs. On the basis of these results and other established risk loci, we used the method developed by Park et al. to estimate the number, and the effect size distribution, of BD risk loci that could still be found by GWAS methods. We estimate that >63,000 case-control samples would be needed to identify the ∼105 BD risk loci discoverable by GWAS, and that these will together explain <6% of the inherited risk. These results support previous GWAS findings and identify three new candidate genes for BD. Further studies are needed to replicate these findings and may potentially lead to identification of functional variants. Sample size will remain a limiting factor in the discovery of common alleles associated with BD.

Multilocus family-based association analysis of seven candidate polymorphisms with essential hypertension in an african-derived semi-isolated brazilian population

Background. It has been widely suggested that analyses considering multilocus effects would be crucial to characterize the relationship between gene variability and essential hypertension (EH). Objective. To test for the presence of multilocus effects between/among seven polymorphisms (six genes) on blood pressure-related traits in African-derived semi-isolated Brazilian populations (quilombos). Methods. Analyses were carried out using a family-based design in a sample of 652 participants (97 families). Seven variants were investigated: ACE (rs1799752), AGT (rs669), ADD2 (rs3755351), NOS3 (rs1799983), GNB3 (rs5441 and rs5443), and GRK4 (rs1801058). Sensitivity analyses were further performed under a case-control design with unrelated participants only. Results. None of the investigated variants were associated individually with both systolic and diastolic BP levels (SBP and DBP, respectively) or EH (as a binary outcome). Multifactor dimensionality reduction-based techniques revealed a marginal association of the combined effect of both GNB3 variants on DBP levels in a family-based design (P = 0.040), whereas a putative NOS3-GRK4 interaction also in relation to DBP levels was observed in the case-control design only (P = 0.004). Conclusion. Our results provide limited support for the hypothesis of multilocus effects between/among the studied variants on blood pressure in quilombos. Further larger studies are needed to validate our findings.

Pourcain B, Xiao Y, Timpson NJ, Smith GD, Reid IR, Ring SM, Sambrook PN, Karlsson M, Dennison EM, Kemp JP, Danoy P, Sayers A, Wilson SG, Nethander M, McCloskey E, Vandenput L, Eastell R, Liu J, Spector T, Mitchell BD, Streeten EA, Brommage R, Pettersson-Kymmer U, Brown MA, Ohlsson C, Richards JB, Lorentzon M. WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)<P<5.9 × 10(-4)) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture.

Challenges and opportunities in genome-wide environmental interaction (GWEI) studies

The interest in performing gene-environment interaction studies has seen a significant increase with the increase of advanced molecular genetics techniques. Practically, it became possible to investigate the role of environmental factors in disease risk and hence to investigate their role as genetic effect modifiers. The understanding that genetics is important in the uptake and metabolism of toxic substances is an example of how genetic profiles can modify important environmental risk factors to disease. Several rationales exist to set up gene-environment interaction studies and the technical challenges related to these studies-when the number of environmental or genetic risk factors is relatively small-has been described before. In the post-genomic era, it is now possible to study thousands of genes and their interaction with the environment. This brings along a whole range of new challenges and opportunities. Despite a continuing effort in developing efficient methods and optimal bioinformatics infrastructures to deal with the available wealth of data, the challenge remains how to best present and analyze genome-wide environmental interaction (GWEI) studies involving multiple genetic and environmental factors. Since GWEIs are performed at the intersection of statistical genetics, bioinformatics and epidemiology, usually similar problems need to be dealt with as for genome-wide association gene-gene interaction studies. However, additional complexities need to be considered which are typical for large-scale epidemiological studies, but are also related to "joining" two heterogeneous types of data in explaining complex disease trait variation or for prediction purposes.

Meta-analysis of genetic association studies under heterogeneity

In multi-cohort genetic association studies or meta-analysis, associations of genetic variants with complex traits across cohorts may be heterogeneous because of genuine genetic diversity or differential biases or errors. To detect the associations of genes with heterogeneous associations across cohorts, new global fixed-effect (FE) and random-effects (RE) meta-analytic methods have been recently proposed. These global methods had improved power over both traditional FE and RE methods under heterogeneity in limited simulation scenarios and data application, but their usefulness in a wide range of practical situations is not clear. We assessed the performance of these methods for both binary and quantitative traits in extensive simulations and applied them to a multi-cohort association study. We found that these new approaches have higher power to detect mostly the very small to small associations of common genetic variants when associations are highly heterogeneous across cohorts. They worked well when both the underlying and assumed genetic models are either multiplicative or dominant. But, they offered no clear advantage for less common variants unless heterogeneity was substantial. In conclusion, these new meta-analytic methods can be used to detect the association of genetic variants with high heterogeneity, which can then be subjected to further exploration, in multi-cohort association studies and meta-analyses.

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10(-8)). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10(-4), Bonferroni corrected), of which six reached P < 5 × 10(-8), including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

Methods for meta-analyses of genome-wide association studies: critical assessment of empirical evidence

There has been a steep increase in the number of meta-analyses of genome-wide association (GWA) studies aimed at identifying genetic variants with increasingly smaller effects, but pressure to publish findings of new genetic associations has limited the time available for careful consideration of all of their methodological aspects. The authors surveyed the literature (2007-2010) to provide empirical evidence on the methods used in GWA meta-analyses, including their organization, requirements about the uniformity of methods used in primary studies, methods for data pooling, investigation of between-study heterogeneity, and quality of reporting. This review showed that a great variety of methods are being used, but the rationale for their choice is often unclear. It also highlights how important methodological aspects have received insufficient attention, potentially leading to missed opportunities for improving gene discovery and characterization. Evaluation of power to replicate findings was inadequate, and the number of variants selected for replication was not associated with replication sample size. A low proportion of GWA meta-analyses investigated the presence and magnitude of heterogeneity, even when there was little uniformity in methods used in primary studies. More methodological work is required before clear guidance can be offered as to optimal methods or tradeoffs between alternative methods. However, there is a clear need for guidelines for reporting the results of GWA meta-analyses.

Interpreting meta-analyses of genome-wide association studies

Meta-analysis is an increasingly popular tool for combining multiple genome-wide association studies in a single analysis to identify associations with small effect sizes. The effect sizes between studies in a meta-analysis may differ and these differences, or heterogeneity, can be caused by many factors. If heterogeneity is observed in the results of a meta-analysis, interpreting the cause of heterogeneity is important because the correct interpretation can lead to a better understanding of the disease and a more effective design of a replication study. However, interpreting heterogeneous results is difficult. The standard approach of examining the association p-values of the studies does not effectively predict if the effect exists in each study. In this paper, we propose a framework facilitating the interpretation of the results of a meta-analysis. Our framework is based on a new statistic representing the posterior probability that the effect exists in each study, which is estimated utilizing cross-study information. Simulations and application to the real data show that our framework can effectively segregate the studies predicted to have an effect, the studies predicted to not have an effect, and the ambiguous studies that are underpowered. In addition to helping interpretation, the new framework also allows us to develop a new association testing procedure taking into account the existence of effect.

Genome-wide association studies are an effective means of identifying genetic variants that are associated with diseases. Although many associated loci have been identified, those loci account for only a small fraction of the genetic contribution to the disease. The remaining contribution may be accounted by loci with very small effect sizes, so small that tens of thousands of samples are needed to identify them. Since it is costly to conduct a study collecting such a large sample, a practical alternative is to combine multiple independent studies in a single analysis called meta-analysis. However, many factors, such as genetic or environmental factors, can differ between the studies combined in a meta-analysis. These factors can cause the effect size of the causal variant to differ between the studies, a phenomenon called heterogeneity. If heterogeneity exists in the data of a meta-analysis, interpreting the meta-analysis results is an important but difficult task. In this paper, we propose a method that helps such interpretation, in addition to a new association testing procedure that is powerful when heterogeneity exists.

Consistency of genome-wide associations across major ancestral groups

It is not well known whether genetic markers identified through genome-wide association studies (GWAS) confer similar or different risks across people of different ancestry. We screened a regularly updated catalog of all published GWAS curated at the NHGRI website for GWAS-identified associations that had reached genome-wide significance (p ≤ 5 × 10(-8)) in at least one major ancestry group (European, Asian, African) and for which replication data were available for comparison in at least two different major ancestry groups. These groups were compared for the correlation between and differences in risk allele frequencies and genetic effects' estimates. Data on 108 eligible GWAS-identified associations with a total of 900 datasets (European, n = 624; Asian, n = 217; African, n = 60) were analyzed. Risk-allele frequencies were modestly correlated between ancestry groups, with >10% absolute differences in 75-89% of the three pairwise comparisons of ancestry groups. Genetic effect (odds ratio) point estimates between ancestry groups correlated modestly (pairwise comparisons' correlation coefficients: 0.20-0.33) and point estimates of risks were opposite in direction or differed more than twofold in 57%, 79%, and 89% of the European versus Asian, European versus African, and Asian versus African comparisons, respectively. The modest correlations, differing risk estimates, and considerable between-association heterogeneity suggest that differential ancestral effects can be anticipated and genomic risk markers may need separate further evaluation in different ancestry groups.

What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations

BACKGROUND

Robust replication is a sine qua non for the rigorous documentation of proposed associations in the genome-wide association (GWA) setting. Currently, associations of common variants reaching P ≤ 5 × 10(-8) are considered replicated. However, there is some ambiguity about the most suitable threshold for claiming genome-wide significance.

METHODS

We defined as 'borderline' associations those with P > 5 × 10(-8) and P ≤ 1 × 10(-7). The eligible associations were retrieved using the 'Catalog of Published Genome-Wide Association Studies'. For each association we assessed whether it reached P ≤ 5 × 10(-8) with inclusion of additional data from subsequent GWA studies.

RESULTS

Thirty-four eligible genotype-phenotype associations were evaluated with data and clarifications contributed from diverse investigators. Replication data from subsequent GWA studies could be obtained for 26 of them. Of those, 19 associations (73%) reached P ≤ 5 × 10(-8) for the same or a related trait implicating either the exact same allele or one in very high linkage disequilibrium and 17 reached P < 10(-8). If the seven associations that did not reach P ≤ 5 × 10(-8) when additional data were considered are assumed to have been false-positives, the false-discovery rate for borderline associations is estimated to be 27% [95% confidence interval (CI) 12-48%]. For five associations, the current P-value is > 10(-6) [corresponding false-discovery rate 19% (95% CI 7-39%)].

CONCLUSION

A substantial proportion, but not all, of the associations with borderline genome-wide significance represent replicable, possibly genuine associations. Our empirical evaluation suggests a possible relaxation in the current GWS threshold.

Genome-wide meta-analysis identifies novel multiple sclerosis susceptibility loci

OBJECTIVE

To perform a 1-stage meta-analysis of genome-wide association studies (GWAS) of multiple sclerosis (MS) susceptibility and to explore functional consequences of new susceptibility loci.

METHODS

We synthesized 7 MS GWAS. Each data set was imputed using HapMap phase II, and a per single nucleotide polymorphism (SNP) meta-analysis was performed across the 7 data sets. We explored RNA expression data using a quantitative trait analysis in peripheral blood mononuclear cells (PBMCs) of 228 subjects with demyelinating disease.

RESULTS

We meta-analyzed 2,529,394 unique SNPs in 5,545 cases and 12,153 controls. We identified 3 novel susceptibility alleles: rs170934(T) at 3p24.1 (odds ratio [OR], 1.17; p = 1.6 × 10(-8)) near EOMES, rs2150702(G) in the second intron of MLANA on chromosome 9p24.1 (OR, 1.16; p = 3.3 × 10(-8)), and rs6718520(A) in an intergenic region on chromosome 2p21, with THADA as the nearest flanking gene (OR, 1.17; p = 3.4 × 10(-8)). The 3 new loci do not have a strong cis effect on RNA expression in PBMCs. Ten other susceptibility loci had a suggestive p < 1 × 10(-6) , some of these loci have evidence of association in other inflammatory diseases (ie, IL12B, TAGAP, PLEK, and ZMIZ1).

INTERPRETATION

We have performed a meta-analysis of GWAS in MS that more than doubles the size of previous gene discovery efforts and highlights 3 novel MS susceptibility loci. These and additional loci with suggestive evidence of association are excellent candidates for further investigations to refine and validate their role in the genetic architecture of MS.

Defining the power limits of genome-wide association scan meta-analyses

Large-scale meta-analyses of genome-wide association scans (GWAS) have been successful in discovering common risk variants with modest and small effects. The detection of lower frequency signals will undoubtedly require concerted efforts of at least similar scale. We investigate the sample size-dictated power limits of GWAS meta-analyses, in the presence and absence of modest levels of heterogeneity and across a range of different allelic architectures. We find that data combination through large-scale collaboration is vital in the quest for complex trait susceptibility loci, but that effect size heterogeneity across meta-analyzed studies drawn from similar populations does not appear to have a profound effect on sample size requirements.

A note on the use of the generalized odds ratio in meta-analysis of association studies involving bi- and tri-allelic polymorphisms

Background

The generalized odds ratio (GOR) was recently suggested as a genetic model-free measure for association studies. However, its properties were not extensively investigated. We used Monte Carlo simulations to investigate type-I error rates, power and bias in both effect size and between-study variance estimates of meta-analyses using the GOR as a summary effect, and compared these results to those obtained by usual approaches of model specification. We further applied the GOR in a real meta-analysis of three genome-wide association studies in Alzheimer's disease.

Findings

For bi-allelic polymorphisms, the GOR performs virtually identical to a standard multiplicative model of analysis (e.g. per-allele odds ratio) for variants acting multiplicatively, but augments slightly the power to detect variants with a dominant mode of action, while reducing the probability to detect recessive variants. Although there were differences among the GOR and usual approaches in terms of bias and type-I error rates, both simulation- and real data-based results provided little indication that these differences will be substantial in practice for meta-analyses involving bi-allelic polymorphisms. However, the use of the GOR may be slightly more powerful for the synthesis of data from tri-allelic variants, particularly when susceptibility alleles are less common in the populations (≤10%). This gain in power may depend on knowledge of the direction of the effects.

Conclusions

For the synthesis of data from bi-allelic variants, the GOR may be regarded as a multiplicative-like model of analysis. The use of the GOR may be slightly more powerful in the tri-allelic case, particularly when susceptibility alleles are less common in the populations.

Multilocus analyses of seven candidate genes suggest interacting pathways for obesity-related traits in Brazilian populations

We investigated whether variants in major candidate genes for food intake and body weight regulation contribute to obesity-related traits under a multilocus perspective. We studied 375 Brazilian subjects from partially isolated African-derived populations (quilombos). Seven variants displaying conflicting results in previous reports and supposedly implicated in the susceptibility of obesity-related phenotypes were investigated: β2-adrenergic receptor (ADRB2) (Arg16Gly), insulin induced gene 2 (INSIG2) (rs7566605), leptin (LEP) (A19G), LEP receptor (LEPR) (Gln223Arg), perilipin (PLIN) (6209T > C), peroxisome proliferator-activated receptor-γ (PPARG) (Pro12Ala), and resistin (RETN) (-420 C > G). Regression models as well as generalized multifactor dimensionality reduction (GMDR) were employed to test the contribution of individual effects and higher-order interactions to BMI and waist-hip ratio (WHR) variation and risk of overweight/obesity. The best multilocus association signal identified in the quilombos was further examined in an independent sample of 334 Brazilian subjects of European ancestry. In quilombos, only the PPARG polymorphism displayed significant individual effects (WHR variation, P = 0.028). No association was observed either with the risk of overweight/obesity (BMI ≥ 25 kg/m2), risk of obesity alone (BMI ≥ 30 kg/m2) or BMI variation. However, GMDR analyses revealed an interaction between the LEPR and ADRB2 polymorphisms (P = 0.009) as well as a third-order effect involving the latter two variants plus INSIG2 (P = 0.034) with overweight/obesity. Assessment of the LEPR-ADRB2 interaction in the second sample indicated a marginally significant association (P = 0.0724), which was further verified to be limited to men (P = 0.0118). Together, our findings suggest evidence for a two-locus interaction between the LEPR Gln223Arg and ADRB2 Arg16Gly variants in the risk of overweight/obesity, and highlight further the importance of multilocus effects in the genetic component of obesity.

Genome-wide meta-analysis identifies regions on 7p21 (AHR) and 15q24 (CYP1A2) as determinants of habitual caffeine consumption

We report the first genome-wide association study of habitual caffeine intake. We included 47,341 individuals of European descent based on five population-based studies within the United States. In a meta-analysis adjusted for age, sex, smoking, and eigenvectors of population variation, two loci achieved genome-wide significance: 7p21 (P = 2.4 × 10(-19)), near AHR, and 15q24 (P = 5.2 × 10(-14)), between CYP1A1 and CYP1A2. Both the AHR and CYP1A2 genes are biologically plausible candidates as CYP1A2 metabolizes caffeine and AHR regulates CYP1A2.

Statistically significant meta-analyses of clinical trials have modest credibility and inflated effects

OBJECTIVE

To assess whether nominally statistically significant effects in meta-analyses of clinical trials are true and whether their magnitude is inflated.

STUDY DESIGN AND SETTING

Data from the Cochrane Database of Systematic Reviews 2005 (issue 4) and 2010 (issue 1) were used. We considered meta-analyses with binary outcomes and four or more trials in 2005 with P<0.05 for the random-effects odds ratio (OR). We examined whether any of these meta-analyses had updated counterparts in 2010. We estimated the credibility (true-positive probability) under different prior assumptions and inflation in OR estimates in 2005.

RESULTS

Four hundred sixty-one meta-analyses in 2005 were eligible, and 80 had additional trials included by 2010. The effect sizes (ORs) were smaller in the updating data (2005-2010) than in the respective meta-analyses in 2005 (median 0.85-fold, interquartile range [IQR]: 0.66-1.06), even more prominently for meta-analyses with less than 300 events in 2005 (median 0.67-fold, IQR: 0.54-0.96). Mean credibility of the 461 meta-analyses in 2005 was 63-84% depending on the assumptions made. Credibility estimates changed >20% in 19-31 (24-39%) of the 80 updated meta-analyses.

CONCLUSIONS

Most meta-analyses with nominally significant results pertain to truly nonnull effects, but exceptions are not uncommon. The magnitude of observed effects, especially in meta-analyses with limited evidence, is often inflated.

Replication of genome-wide discovered breast cancer risk loci in the Cypriot population

Genome-wide association studies (GWAS) have identified associations with robust statistical support for influencing breast cancer susceptibility. Most GWAS and replications have been conducted in Northern European populations and to a lesser extent in Asians, and Ashkenazi Jews. It is important to evaluate whether these variants confer risk across different populations and also to assess the magnitude of risk conferred. The aim of this study was to evaluate previously GWAS-identified breast cancer risk variants in the Cypriot population. Eleven GWAS-discovered single nucleotide polymorphisms (SNPs) were analyzed for association with breast cancer in 1,109 Cypriot female breast cancer patients and 1,177 healthy female controls. Four of the 11 SNPs evaluated were found to be nominally significantly associated (P < 0.05) with breast cancer risk in the Cypriot population. Based on estimated power, five associations would be expected to be nominally significant. The correlation coefficient of effect sizes (per-allele odds ratio) between the Cypriot population and the original GWAS populations where these SNPs had been discovered was 0.58 (P = 0.064), while allele frequencies were very similar (r = 0.88, P < 0.001). Overall, we show modest concordance for breast cancer GWAS-discovered alleles and their effect sizes in the Cypriot population. The effects sizes of GWAS-discovered SNPs need to be verified separately in different populations.

Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

BACKGROUND

Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naïve strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication.

METHODS

Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (τ²).

RESULTS

Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to τ² whenever the susceptibility allele is common (MAF ≥ 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model.

CONCLUSION

Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.

A compendium of genome-wide associations for cancer: critical synopsis and reappraisal

Since 2007, genome-wide association (GWA) studies have identified numerous well-supported, novel genetic risk loci for common cancers; however, there are concerns that this technology is reaching its limits. We provide an overview of GWA-identified genetic associations with solid tumors. We simulated the distribution of population risk alleles for colorectal, prostate, testicular, and thyroid cancers based on genetic variants identified in GWA studies. We also evaluated whether statistical power to detect typical genetic effects could be improved with studies performing GWA analyses of all available samples rather than multistage designs. Fifty-six eligible articles yielded 92 eligible associations between cancer phenotypes and genetic variants with a median per-allele odds ratio (OR) of 1.22 (interquartile range = 1.15-1.36). Half of the associations pertained to prostate, colorectal, or breast cancer. Individuals at the upper quartile of simulated risk had only 2.1- to 4.2-fold higher relative risk than those in the lower quartile. Comprehensive evaluation of currently available samples with GWA platforms would yield few additional variants with per-allele OR = 1.4, but many more variants with OR = 1.2 could be detected; statistical power to detect weak associations (OR = 1.07) would still be negligible. The GWA approach is effective in identifying common genetic variants with moderate effect; however, identifying loci with very small effects and rare variants will require major new efforts. At present, the utility of GWA-identified risk loci in risk stratification for cancer is limited.

The Gene, Environment Association Studies consortium (GENEVA): maximizing the knowledge obtained from GWAS by collaboration across studies of multiple conditions

Genome-wide association studies (GWAS) have emerged as powerful means for identifying genetic loci related to complex diseases. However, the role of environment and its potential to interact with key loci has not been adequately addressed in most GWAS. Networks of collaborative studies involving different study populations and multiple phenotypes provide a powerful approach for addressing the challenges in analysis and interpretation shared across studies. The Gene, Environment Association Studies (GENEVA) consortium was initiated to: identify genetic variants related to complex diseases; identify variations in gene-trait associations related to environmental exposures; and ensure rapid sharing of data through the database of Genotypes and Phenotypes. GENEVA consists of several academic institutions, including a coordinating center, two genotyping centers and 14 independently designed studies of various phenotypes, as well as several Institutes and Centers of the National Institutes of Health led by the National Human Genome Research Institute. Minimum detectable effect sizes include relative risks ranging from 1.24 to 1.57 and proportions of variance explained ranging from 0.0097 to 0.02. Given the large number of research participants (N>80,000), an important feature of GENEVA is harmonization of common variables, which allow analyses of additional traits. Environmental exposure information available from most studies also enables testing of gene-environment interactions. Facilitated by its sizeable infrastructure for promoting collaboration, GENEVA has established a unified framework for genotyping, data quality control, analysis and interpretation. By maximizing knowledge obtained through collaborative GWAS incorporating environmental exposure information, GENEVA aims to enhance our understanding of disease etiology, potentially identifying opportunities for intervention.

Genome-wide meta-analyses identify multiple loci associated with smoking behavior

Consistent but indirect evidence has implicated genetic factors in smoking behavior. We report meta-analyses of several smoking phenotypes within cohorts of the Tobacco and Genetics Consortium (n = 74,053). We also partnered with the European Network of Genetic and Genomic Epidemiology (ENGAGE) and Oxford-GlaxoSmithKline (Ox-GSK) consortia to follow up the 15 most significant regions (n > 140,000). We identified three loci associated with number of cigarettes smoked per day. The strongest association was a synonymous 15q25 SNP in the nicotinic receptor gene CHRNA3 (rs1051730[A], beta = 1.03, standard error (s.e.) = 0.053, P = 2.8 x 10(-73)). Two 10q25 SNPs (rs1329650[G], beta = 0.367, s.e. = 0.059, P = 5.7 x 10(-10); and rs1028936[A], beta = 0.446, s.e. = 0.074, P = 1.3 x 10(-9)) and one 9q13 SNP in EGLN2 (rs3733829[G], beta = 0.333, s.e. = 0.058, P = 1.0 x 10(-8)) also exceeded genome-wide significance for cigarettes per day. For smoking initiation, eight SNPs exceeded genome-wide significance, with the strongest association at a nonsynonymous SNP in BDNF on chromosome 11 (rs6265[C], odds ratio (OR) = 1.06, 95% confidence interval (Cl) 1.04-1.08, P = 1.8 x 10(-8)). One SNP located near DBH on chromosome 9 (rs3025343[G], OR = 1.12, 95% Cl 1.08-1.18, P = 3.6 x 10(-8)) was significantly associated with smoking cessation.

Genome-wide association studies--data generation, storage, interpretation, and bioinformatics

Genome-wide association studies (GWAS) have had great success in identifying common genetic determinants of disease. One of the challenges posed by GWAS is the analysis of the large amount of data generated. This review aims to provide the non-geneticists with an overview of the different steps entailed in analysis of GWAS data, with an emphasis on popular bioinformatics tools available. GWAS data generation, analysis, and interpretation will be covered.

Advances in genomic analysis of stroke: what have we learned and where are we headed?

As a result of technological advances, the genomic analysis of stroke has shifted from candidate gene association studies to genome-wide association studies (GWAS). Agnostic GWAS evaluate up to 90% of common genetic variation in a single experiment, creating an improved framework for identifying novel genetic leads for biochemical and cellular mechanisms underlying stroke. Given the ubiquity of the GWAS approach, it has become essential for stroke researchers and clinicians to be able to interpret GWAS results. Thus, we review the basic elements of design, methods, presentation, and interpretation of GWAS in the context of stroke research. In 8 recent stroke GWAS reports, no single locus has been identified in 2 GWAS at a genome-wide level of significance. Additionally, no significant association signal between stroke and a locus with previous evidence from candidate gene studies of stroke has been identified yet. Some caveats of the approach and future directions for stroke genomics are discussed, including the use of intermediate phenotypes, Mendelian randomization, phenomics, and deep resequencing. Intelligent, appropriately powered, multidisciplinary studies incorporating knowledge from clinical medicine, epidemiology, genetics, and molecular biology will be required to fully characterize the genomic contributors to stroke.