Weang Kee Ho

Genetic risk factors for ischaemic stroke and its subtypes (the METASTROKE collaboration): a meta-analysis of genome-wide association studies.

Summary Background Various genome-wide association studies (GWAS) have been done in ischaemic stroke, identifying a few loci associated with the disease, but sample sizes have been 3500 cases or less. We established the METASTROKE collaboration with the aim of validating associations from previous GWAS and identifying novel genetic associations through meta-analysis of GWAS datasets for ischaemic stroke and its subtypes. Methods We meta-analysed data from 15 ischaemic stroke cohorts with a total of 12 389 individuals with ischaemic stroke and 62 004 controls, all of European ancestry. For the associations reaching genome-wide significance in METASTROKE, we did a further analysis, conditioning on the lead single nucleotide polymorphism in every associated region. Replication of novel suggestive signals was done in 13 347 cases and 29 083 controls. Findings We verified previous associations for cardioembolic stroke near PITX2 (p=2·8×10−16) and ZFHX3 (p=2·28×10−8), and for large-vessel stroke at a 9p21 locus (p=3·32×10−5) and HDAC9 (p=2·03×10−12). Additionally, we verified that all associations were subtype specific. Conditional analysis in the three regions for which the associations reached genome-wide significance (PITX2, ZFHX3, and HDAC9) indicated that all the signal in each region could be attributed to one risk haplotype. We also identified 12 potentially novel loci at p<5×10−6. However, we were unable to replicate any of these novel associations in the replication cohort. Interpretation Our results show that, although genetic variants can be detected in patients with ischaemic stroke when compared with controls, all associations we were able to confirm are specific to a stroke subtype. This finding has two implications. First, to maximise success of genetic studies in ischaemic stroke, detailed stroke subtyping is required. Second, different genetic pathophysiological mechanisms seem to be associated with different stroke subtypes. Funding Wellcome Trust, UK Medical Research Council (MRC), Australian National and Medical Health Research Council, National Institutes of Health (NIH) including National Heart, Lung and Blood Institute (NHLBI), the National Institute on Aging (NIA), the National Human Genome Research Institute (NHGRI), and the National Institute of Neurological Disorders and Stroke (NINDS).

Genome-Wide Association Study Identifies a Novel Locus Contributing to Type 2 Diabetes Susceptibility in Sikhs of Punjabi Origin From India

We performed a genome-wide association study (GWAS) and a multistage meta-analysis of type 2 diabetes (T2D) in Punjabi Sikhs from India. Our discovery GWAS in 1,616 individuals (842 case subjects) was followed by in silico replication of the top 513 independent single nucleotide polymorphisms (SNPs) (P < 10−3) in Punjabi Sikhs (n = 2,819; 801 case subjects). We further replicated 66 SNPs (P < 10−4) through genotyping in a Punjabi Sikh sample (n = 2,894; 1,711 case subjects). On combined meta-analysis in Sikh populations (n = 7,329; 3,354 case subjects), we identified a novel locus in association with T2D at 13q12 represented by a directly genotyped intronic SNP (rs9552911, P = 1.82 × 10−8) in the SGCG gene. Next, we undertook in silico replication (stage 2b) of the top 513 signals (P < 10−3) in 29,157 non-Sikh South Asians (10,971 case subjects) and de novo genotyping of up to 31 top signals (P < 10−4) in 10,817 South Asians (5,157 case subjects) (stage 3b). In combined South Asian meta-analysis, we observed six suggestive associations (P < 10−5 to < 10−7), including SNPs at HMG1L1/CTCFL, PLXNA4, SCAP, and chr5p11. Further evaluation of 31 top SNPs in 33,707 East Asians (16,746 case subjects) (stage 3c) and 47,117 Europeans (8,130 case subjects) (stage 3d), and joint meta-analysis of 128,127 individuals (44,358 case subjects) from 27 multiethnic studies, did not reveal any additional loci nor was there any evidence of replication for the new variant. Our findings provide new evidence on the presence of a population-specific signal in relation to T2D, which may provide additional insights into T2D pathogenesis.

Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms

Coronary artery disease (CAD) is a leading cause of morbidity and mortality worldwide. Although 58 genomic regions have been associated with CAD thus far, most of the heritability is unexplained, indicating that additional susceptibility loci await identification. An efficient discovery strategy may be larger-scale evaluation of promising associations suggested by genome-wide association studies (GWAS). Hence, we genotyped 56,309 participants using a targeted gene array derived from earlier GWAS results and performed meta-analysis of results with 194,427 participants previously genotyped, totaling 88,192 CAD cases and 162,544 controls. We identified 25 new SNP–CAD associations (P < 5 × 10−8, in fixed-effects meta-analysis) from 15 genomic regions, including SNPs in or near genes involved in cellular adhesion, leukocyte migration and atherosclerosis (PECAM1, rs1867624), coagulation and inflammation (PROCR, rs867186 (p.Ser219Gly)) and vascular smooth muscle cell differentiation (LMOD1, rs2820315). Correlation of these regions with cell-type-specific gene expression and plasma protein levels sheds light on potential disease mechanisms.

GUESS-ing polygenic associations with multiple phenotypes using a GPU-based evolutionary stochastic search algorithm

Genome-wide association studies (GWAS) yielded significant advances in defining the genetic architecture of complex traits and disease. Still, a major hurdle of GWAS is narrowing down multiple genetic associations to a few causal variants for functional studies. This becomes critical in multi-phenotype GWAS where detection and interpretability of complex SNP(s)-trait(s) associations are complicated by complex Linkage Disequilibrium patterns between SNPs and correlation between traits. Here we propose a computationally efficient algorithm (GUESS) to explore complex genetic-association models and maximize genetic variant detection. We integrated our algorithm with a new Bayesian strategy for multi-phenotype analysis to identify the specific contribution of each SNP to different trait combinations and study genetic regulation of lipid metabolism in the Gutenberg Health Study (GHS). Despite the relatively small size of GHS (n = 3,175), when compared with the largest published meta-GWAS (n>100,000), GUESS recovered most of the major associations and was better at refining multi-trait associations than alternative methods. Amongst the new findings provided by GUESS, we revealed a strong association of SORT1 with TG-APOB and LIPC with TG-HDL phenotypic groups, which were overlooked in the larger meta-GWAS and not revealed by competing approaches, associations that we replicated in two independent cohorts. Moreover, we demonstrated the increased power of GUESS over alternative multi-phenotype approaches, both Bayesian and non-Bayesian, in a simulation study that mimics real-case scenarios. We showed that our parallel implementation based on Graphics Processing Units outperforms alternative multi-phenotype methods. Beyond multivariate modelling of multi-phenotypes, our Bayesian model employs a flexible hierarchical prior structure for genetic effects that adapts to any correlation structure of the predictors and increases the power to identify associated variants. This provides a powerful tool for the analysis of diverse genomic features, for instance including gene expression and exome sequencing data, where complex dependencies are present in the predictor space.

Comparative review of methods for handling drop‐out in longitudinal studies

Longitudinal data analysis is frequently complicated by drop-out. In this paper we consider several methods for dealing with drop-out afflicted data. Along with a general comparison, particular attention is paid to the consequences of model misspecification. The purpose of our approach is two-fold. We first deliberate the form of the drop-out model and compare two alternatives. Furthermore, the extent to which each method is dependent on its core assumptions is assessed through scenarios where one or more such assumptions are compromised. Second, the extent to which we can identify adequacy of model fit is investigated via recently developed diagnostics. These twin targets are pursued via simulation scenarios and application to a schizophrenia trial of over 500 patients with near 50 per cent drop-out.

Candidate Gene Association Study of Coronary Artery Calcification in Chronic Kidney Disease: Findings From the CRIC Study (Chronic Renal Insufficiency Cohort)

OBJECTIVES This study sought to identify loci for coronary artery calcification (CAC) in patients with chronic kidney disease (CKD). BACKGROUND CKD is associated with increased CAC and subsequent coronary heart disease (CHD), but the mechanisms remain poorly defined. Genetic studies of CAC in CKD may provide a useful strategy for identifying novel pathways in CHD. METHODS We performed a candidate gene study (∼2,100 genes; ∼50,000 single nucleotide polymorphisms [SNPs]) of CAC within the CRIC (Chronic Renal Insufficiency Cohort) study (N = 1,509; 57% European, 43% African ancestry). SNPs with preliminary evidence of association with CAC in CRIC were examined for association with CAC in the PennCAC (Penn Coronary Artery Calcification) (N = 2,560) and AFCS (Amish Family Calcification Study) (N = 784) samples. SNPs with suggestive replication were further analyzed for association with myocardial infarction (MI) in the PROMIS (Pakistan Risk of Myocardial Infarction Study) (N = 14,885). RESULTS Of 268 SNPs reaching p < 5 × 10(-4) for CAC in CRIC, 28 SNPs in 23 loci had nominal support (p < 0.05 and in same direction) for CAC in PennCAC or AFCS. Besides chr9p21 and COL4A1, known loci for CHD, these included SNPs having reported genome-wide association study association with hypertension (e.g., ATP2B1). In PROMIS, 4 of the 23 suggestive CAC loci (chr9p21, COL4A1, ATP2B1, and ABCA4) had significant associations with MI, consistent with their direction of effect on CAC. CONCLUSIONS We identified several loci associated with CAC in CKD that also relate to MI in a general population sample. CKD imparts a high risk of CHD and may provide a useful setting for discovery of novel CHD genes and pathways.

Genome-Wide Association Analysis of Young-Onset Stroke Identifies a Locus on Chromosome 10q25 Near HABP2

Background and Purpose— Although a genetic contribution to ischemic stroke is well recognized, only a handful of stroke loci have been identified by large-scale genetic association studies to date. Hypothesizing that genetic effects might be stronger for early- versus late-onset stroke, we conducted a 2-stage meta-analysis of genome-wide association studies, focusing on stroke cases with an age of onset <60 years. Methods— The discovery stage of our genome-wide association studies included 4505 cases and 21 968 controls of European, South-Asian, and African ancestry, drawn from 6 studies. In Stage 2, we selected the lead genetic variants at loci with association P<5×10−6 and performed in silico association analyses in an independent sample of ⩽1003 cases and 7745 controls. Results— One stroke susceptibility locus at 10q25 reached genome-wide significance in the combined analysis of all samples from the discovery and follow-up stages (rs11196288; odds ratio =1.41; P=9.5×10−9). The associated locus is in an intergenic region between TCF7L2 and HABP2. In a further analysis in an independent sample, we found that 2 single nucleotide polymorphisms in high linkage disequilibrium with rs11196288 were significantly associated with total plasma factor VII–activating protease levels, a product of HABP2. Conclusions— HABP2, which encodes an extracellular serine protease involved in coagulation, fibrinolysis, and inflammatory pathways, may be a genetic susceptibility locus for early-onset stroke.

Loss of Cardio-Protective Effects at the ADAMTS7 Locus Due to Gene-Smoking Interactions.

Background: Common diseases such as coronary heart disease (CHD) are complex in etiology. The interaction of genetic susceptibility with lifestyle factors may play a prominent role. However, gene-lifestyle interactions for CHD have been difficult to identify. Here, we investigate interaction of smoking behavior, a potent lifestyle factor, with genotypes that have been shown to associate with CHD risk. Methods: We analyzed data on 60 919 CHD cases and 80 243 controls from 29 studies for gene-smoking interactions for genetic variants at 45 loci previously reported to be associated with CHD risk. We also studied 5 loci associated with smoking behavior. Study-specific gene-smoking interaction effects were calculated and pooled using fixed-effects meta-analyses. Interaction analyses were declared to be significant at a P value of <1.0×10–3 (Bonferroni correction for 50 tests). Results: We identified novel gene-smoking interaction for a variant upstream of the ADAMTS7 gene. Every T allele of rs7178051 was associated with lower CHD risk by 12% in never-smokers (P=1.3×10–16) in comparison with 5% in ever-smokers (P=2.5×10–4), translating to a 60% loss of CHD protection conferred by this allelic variation in people who smoked tobacco (interaction P value=8.7×10–5). The protective T allele at rs7178051 was also associated with reduced ADAMTS7 expression in human aortic endothelial cells and lymphoblastoid cell lines. Exposure of human coronary artery smooth muscle cells to cigarette smoke extract led to induction of ADAMTS7. Conclusions: Allelic variation at rs7178051 that associates with reduced ADAMTS7 expression confers stronger CHD protection in never-smokers than in ever-smokers. Increased vascular ADAMTS7 expression may contribute to the loss of CHD protection in smokers.

Dropouts in the AB/BA crossover design

Missing data arise in crossover trials, as they do in any form of clinical trial. Several papers have addressed the problems that missing data create, although almost all of these assume that the probability that a planned observation is missing does not depend on the value that would have been observed; that is, the data are missing at random (MAR). In many applications, this assumption is likely to be untenable; in which case, the data are missing not at random (MNAR). We investigate the effect on estimates of the treatment effect that assume data are MAR when data are actually MNAR. We also propose using the assumption of no carryover treatment effect, which is usually required for this design, to permit the estimation of a treatment effect when data are MNAR. The results are applied to a trial comparing two treatments for neuropathic pain and show that the estimate of treatment effect is sensitive to the assumption of MAR.

Dropout in crossover and longitudinal studies: Is complete case so bad?

We discuss inference for longitudinal clinical trials subject to possibly informative dropout. A selection of available methods is reviewed for the simple case of trials with two timepoints. Using data from two such clinical trials, each with two treatments, we demonstrate that different analysis methods can at times lead to quite different conclusions from the same data. We investigate properties of complete-case estimators for the type of trials considered, with emphasis on interpretation and meaning of parameters. We contrast longitudinal and crossover designs and argue that for crossover studies there are often good reasons to prefer a complete case analysis. More generally, we suggest that there is merit in an approach in which no untestable assumptions are made. Such an approach would combine a dropout analysis, an analysis of complete-case data only, and a careful statement of justified conclusions.