Yanxia Wu

Regulation of CARD8 Expression by ANRIL and Association of CARD8 Single Nucleotide Polymorphism rs2043211 (p.C10X) With Ischemic Stroke

Background and Purpose— ANRIL has long been considered as the strongest candidate gene at the 9p21 locus, robustly associated with stroke and coronary artery disease. However, the underlying molecular mechanism remains unknown. The present study works to elucidate such a mechanism. Methods— Using expression quantitative loci analysis, we identified potential genes whose expression may be influenced by genetic variation in ANRIL. To verify the identified gene(s), knockdown and overexpression of ANRIL were evaluated in human umbilical vein endothelial cells and HepG2 cells. Ischemic stroke and coronary artery disease risk were then evaluated in the gene(s) demonstrated to be mediated by ANRIL in 3 populations of Chinese Han ancestry: 2 ischemic stroke populations consisting of the Central China cohort (903 cases and 873 controls) and the Northern China cohort (816 cases and 879 controls) and 1 coronary artery disease cohort consisting of 772 patients and 873 controls. Results— Expression quantitative loci analysis identified CARD8 among others, with knockdown of ANRIL expression decreasing CARD8 expression and overexpression of ANRIL increasing CARD8 expression. The minor T allele of a previously identified CARD8 variant (rs2043211) was found to be significantly associated with a protective effect of ischemic stroke under the recessive model in 2 independent stroke cohorts. No significant association was found between rs2043211 and coronary artery disease. Conclusions— CARD8 is a downstream target gene regulated by ANRIL. Single nucleotide polymorphism rs2043211 in CARD8 is significantly associated with ischemic stroke. ANRIL may increase the risk of ischemic stroke through regulation of the CARD8 pathway.

The Same Chromosome 9p21.3 Locus Is Associated With Type 2 Diabetes and Coronary Artery Disease in a Chinese Han Population

OBJECTIVE Recent genome-wide association studies (GWAS) revealed that a 9p21.3 locus was associated with type 2 diabetes. In this study, we carried out a large-scale case-control study in the GeneID Chinese Han population to 1) further replicate the association of 9p21.3 type 2 diabetes GWAS single nucleotide polymorphisms (SNPs) and 2) assess the association of these SNPs with coronary artery disease. RESEARCH DESIGN AND METHODS Three SNPs (rs2383208, rs10811661, and rs10757283) were genotyped in two GeneID cohorts of 3,167 Chinese Han individuals. Case-control association design was used to determine the association of the SNPs with type 2 diabetes and coronary artery disease. Gensini scores were calculated in the coronary artery disease subjects and were tested for association with the variants. Multivariate logistic regressions were performed on association studies. RESULTS The association between two of the three SNPs and type 2 diabetes was replicated in the GeneID population (rs2383208, P = 0.936; rs10811661-T, P = 0.02, odds ratio [OR] = 1.23; rs10757283-C, P = 0.003, OR = 1.30). The same two SNPs also contributed to the risk of coronary artery disease (CAD) (rs10811661-T, P = 0.002, OR = 1.19; rs10757283-C, P = 0.003, OR = 1.18). In addition, rs10757283 was associated with severity of coronary atherosclerosis estimated by the Gensini scoring system (risk allele C, quantitative-trait regression adjusted P = 0.002). CONCLUSIONS For the first time to our knowledge, our results indicated that the same 9p21.3 locus, represented by SNPs rs10811661 and rs10757283, contributed to the risk of type 2 diabetes and coronary artery disease in our GeneID Chinese Han population.

Molecular Basis of Gene-Gene Interaction: Cyclic Cross-Regulation of Gene Expression and Post-GWAS Gene-Gene Interaction Involved in Atrial Fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia at the clinic. Recent GWAS identified several variants associated with AF, but they account for <10% of heritability. Gene-gene interaction is assumed to account for a significant portion of missing heritability. Among GWAS loci for AF, only three were replicated in the Chinese Han population, including SNP rs2106261 (G/A substitution) in ZFHX3, rs2200733 (C/T substitution) near PITX2c, and rs3807989 (A/G substitution) in CAV1. Thus, we analyzed the interaction among these three AF loci. We demonstrated significant interaction between rs2106261 and rs2200733 in three independent populations and combined population with 2,020 cases/5,315 controls. Compared to non-risk genotype GGCC, two-locus risk genotype AATT showed the highest odds ratio in three independent populations and the combined population (OR=5.36 (95% CI 3.87-7.43), P=8.00×10-24). The OR of 5.36 for AATT was significantly higher than the combined OR of 3.31 for both GGTT and AACC, suggesting a synergistic interaction between rs2106261 and rs2200733. Relative excess risk due to interaction (RERI) analysis also revealed significant interaction between rs2106261 and rs2200733 when exposed two copies of risk alleles (RERI=2.87, P<1.00×10-4) or exposed to one additional copy of risk allele (RERI=1.29, P<1.00×10-4). The INTERSNP program identified significant genotypic interaction between rs2106261 and rs2200733 under an additive by additive model (OR=0.85, 95% CI: 0.74-0.97, P=0.02). Mechanistically, PITX2c negatively regulates expression of miR-1, which negatively regulates expression of ZFHX3, resulting in a positive regulation of ZFHX3 by PITX2c; ZFHX3 positively regulates expression of PITX2C, resulting in a cyclic loop of cross-regulation between ZFHX3 and PITX2c. Both ZFHX3 and PITX2c regulate expression of NPPA, TBX5 and NKX2.5. These results suggest that cyclic cross-regulation of gene expression is a molecular basis for gene-gene interactions involved in genetics of complex disease traits.

Candidate Pathway-Based Genome-Wide Association Studies Identify Novel Associations of Genomic Variants in the Complement System Associated With Coronary Artery Disease

Background—Genomic variants identified by genome-wide association studies (GWAS) explain <20% of heritability of coronary artery disease (CAD), thus many risk variants remain missing for CAD. Identification of new variants may unravel new biological pathways and genetic mechanisms for CAD. To identify new variants associated with CAD, we developed a candidate pathway-based GWAS by integrating expression quantitative loci analysis and mining of GWAS data with variants in a candidate pathway. Methods and Results—Mining of GWAS data was performed to analyze variants in 32 complement system genes for positive association with CAD. Functional variants in genes showing positive association were then identified by searching existing expression quantitative loci databases and validated by real-time reverse transcription polymerase chain reaction. A follow-up case–control design was then used to determine whether the functional variants are associated with CAD in 2 independent GeneID Chinese populations. Candidate pathway-based GWAS identified positive association between variants in C3AR1 and C6 and CAD. Two functional variants, rs7842 in C3AR1 and rs4400166 in C6, were found to be associated with expression levels of C3AR1 and C6, respectively. Significant association was identified between rs7842 and CAD (P=3.99×10−6; odds ratio, 1.47) and between rs4400166 and CAD (P=9.30×10−3; odds ratio, 1.24) in the validation cohort. The significant findings were confirmed in the replication cohort (P=1.53×10−5; odds ratio, 1.37 for rs7842; P=8.41×10−3; odds ratio, 1.21 for rs4400166). Conclusions—Integration of GWAS with biological pathways and expression quantitative loci is effective in identifying new risk variants for CAD. Functional variants increasing C3AR1 and C6 expression were shown to confer significant risk of CAD for the first time.

Association of SNP Rs9943582 in APLNR with Left Ventricle Systolic Dysfunction in Patients with Coronary Artery Disease in a Chinese Han GeneID Population.

Heart failure affects 1–2% of the adult population worldwide and coronary artery disease (CAD) is the underlying etiology of heart failure in 70% of the patients. The pathway of apelin and its apelin receptor (APJ) was implicated in the pathogenesis of heart failure in animal models, but a similar role in humans is unknown. We studied a functional variant, rs9943582 (-154G/A), at the 5’-untranslated region, that was associated with decreased expression of the APJ receptor gene (APLNR) in a population consisting of 1,751 CAD cases and 1,022 controls. Variant rs9943582 was not associated with CAD, but among CAD patients, it showed significant association with left ventricular systolic dysfunction (431 CAD patients with left ventricular systolic dysfunction (LV ejection fraction or LVEF< 40%) versus 1,046 CAD patients without LV systolic dysfunction (LVEF>50%) (P-adj = 6.71×10-5, OR = 1.43, 95% CI, 1.20–1.70). Moreover, rs9943582 also showed significant association with quantitative echocardiographic parameters, including left ventricular end-diastolic diameter (effect size: increased 1.67±0.43 mm per risk allele A, P = 1.15×10-4), left atrial size (effect size: increased 2.12±0.61 mm per risk allele A, P = 9.56×10-4) and LVEF (effect size: decreased 2.59±0.32 percent per risk allele A, P = 7.50×10-15). Our findings demonstrate that allele A of rs9943582 was significantly associated with left ventricular systolic dysfunction, left ventricular end-diastolic diameter, the left atrial diameter and LVEF in the CAD population, which suggests an important role of the apelin/APJ system in the pathology of heart failure associated with ischemic heart disease.

Analysis of the genetic association between IL27 variants and coronary artery disease in a Chinese Han population.

Interleukin-27 (IL-27) is an important cytokine in inflammatory diseases, including coronary artery disease (CAD). To explore the precise role of IL-27 in CAD, we investigated the genetic association between IL27 and CAD in the GeneID Chinese Han population. A two-stage case control association analysis was performed for 3075 CAD cases and 2802 controls. Logistic regression analysis was used to adjust the traditional risk factors for CAD. Results showed that a promoter variant, rs153109, tended to be marginally associated with CAD in the discovery population (Padj = 0.028, OR = 1.27, 95%CI: 1.03–1.58). However, this association was not replicated in the validation stage (Padj = 0.559, OR = 1.04, 95%CI: 0.90–1.21). In addition, when we classified the combined population into two subgroups according to the age at disease onset or disease state, we again obtained no significant associations. Finally, we estimated the severity of coronary stenosis using the Gensini Scoring system and determined that the rs153109 genotypes were still not associated with the Gensini scores of the CAD patients. In conclusion, our study failed to find an association between common variants in the functional region of IL27 and CAD in a Chinese Han population, which indicated that IL-27 might only be an inflammatory marker during the development of CAD.

Candidate Pathway-Based GWAS Identifies Novel Associations of Genomic Variants in the Complement System Associated with Coronary Artery Disease

Background—Genomic variants identified by genome-wide association studies (GWAS) explain <20% of heritability of coronary artery disease (CAD), thus many risk variants remain missing for CAD. Identification of new variants may unravel new biological pathways and genetic mechanisms for CAD. To identify new variants associated with CAD, we developed a candidate pathway-based GWAS by integrating expression quantitative loci analysis and mining of GWAS data with variants in a candidate pathway. Methods and Results—Mining of GWAS data was performed to analyze variants in 32 complement system genes for positive association with CAD. Functional variants in genes showing positive association were then identified by searching existing expression quantitative loci databases and validated by real-time reverse transcription polymerase chain reaction. A follow-up case–control design was then used to determine whether the functional variants are associated with CAD in 2 independent GeneID Chinese populations. Candidate pathway-based GWAS identified positive association between variants in C3AR1 and C6 and CAD. Two functional variants, rs7842 in C3AR1 and rs4400166 in C6, were found to be associated with expression levels of C3AR1 and C6, respectively. Significant association was identified between rs7842 and CAD (P=3.99×10−6; odds ratio, 1.47) and between rs4400166 and CAD (P=9.30×10−3; odds ratio, 1.24) in the validation cohort. The significant findings were confirmed in the replication cohort (P=1.53×10−5; odds ratio, 1.37 for rs7842; P=8.41×10−3; odds ratio, 1.21 for rs4400166). Conclusions—Integration of GWAS with biological pathways and expression quantitative loci is effective in identifying new risk variants for CAD. Functional variants increasing C3AR1 and C6 expression were shown to confer significant risk of CAD for the first time.

Genomic Variants in NEURL, GJA1 and CUX2 Significantly Increase Genetic Susceptibility to Atrial Fibrillation.

Atrial fibrillation (AF) is the most common arrhythmia. In 2014, two new meta-GWAS identified 5 AF loci, including the NEURL locus, GJA1 locus, CAND2 locus, and TBX5 locus in the European ancestry populations and the NEURL locus and CUX2 locus in a Japanese population. The TBX5 locus for AF was reported by us in 2013 in the Chinese population. Here we assessed the association between AF and SNPs in the NEURL, GJA1, CAND2 and CUX2 loci in the Chinese Han population. We carried out a large case-control association study with 1,164 AF patients and 1,460 controls. Significant allelic and genotypic associations were identified between NEURL variant rs6584555 and GJA1 variant rs13216675 and AF. Significant genotypic association was found between CUX2 SNP rs6490029 and AF. No association was found between CAND2 variant rs4642101 and AF, which may be due to an insufficient power of the sample size for rs4642101. Together with our previous findings, seven of fifteen AF loci (<50%) identified by GWAS in the European ancestry populations conferred susceptibility to AF in the Chinese population, and explained approximately 14.5% of AF heritability. On the other hand, two AF loci identified in the Japanese population were both replicated in the Chinese population.

Genetic Regulation of the Thymic Stromal Lymphopoietin (TSLP)/TSLP Receptor (TSLPR) Gene Expression and Influence of Epistatic Interactions Between IL-33 and the TSLP/TSLPR Axis on Risk of Coronary Artery Disease

The thymic stromal lymphopoietin (TSLP)/TSLP receptor (TSLPR) axis is involved in multiple inflammatory immune diseases, including coronary artery disease (CAD). To explore the causal relationship between this axis and CAD, we performed a three-stage case-control association analysis with 3,628 CAD cases and 3,776 controls using common variants in the genes TSLP, interleukin 7 receptor (IL7R), and TSLPR. Three common variants in the TSLP/TSLPR axis were significantly associated with CAD in a Chinese Han population [rs3806933T in TSLP, Padj = 4.35 × 10−5, odds ratio (OR) = 1.18; rs6897932T in IL7R, Padj = 1.13 × 10−7, OR = 1.31; g.19646A>GA in TSLPR, Padj = 2.04 × 10−6, OR = 1.20]. Reporter gene analysis demonstrated that rs3806933 and rs6897932 could influence TSLP and IL7R expression, respectively. Furthermore, the “T” allele of rs3806933 might increase plasma TSLP levels (R2 = 0.175, P < 0.01). In a stepwise procedure, the risk for CAD increased by nearly fivefold compared with the maximum effect of any single variant (Padj = 6.99 × 10−4, OR = 4.85). In addition, the epistatic interaction between TSLP and IL33 produced a nearly threefold increase in the risk of CAD in the combined model of rs3806933TT-rs7025417TT (Padj = 3.67 × 10−4, OR = 2.98). Our study illustrates that the TSLP/TSLPR axis might be involved in the pathogenesis of CAD through upregulation of mRNA or protein expression of the referenced genes and might have additive effects on the CAD risk when combined with IL-33 signaling.

Candidate Pathway-Based Genome-Wide Association Studies Identify Novel Associations of Genomic Variants in the Complement System Associated With Coronary Artery DiseaseCLINICAL PERSPECTIVE

Background—Genomic variants identified by genome-wide association studies (GWAS) explain <20% of heritability of coronary artery disease (CAD), thus many risk variants remain missing for CAD. Identification of new variants may unravel new biological pathways and genetic mechanisms for CAD. To identify new variants associated with CAD, we developed a candidate pathway-based GWAS by integrating expression quantitative loci analysis and mining of GWAS data with variants in a candidate pathway. Methods and Results—Mining of GWAS data was performed to analyze variants in 32 complement system genes for positive association with CAD. Functional variants in genes showing positive association were then identified by searching existing expression quantitative loci databases and validated by real-time reverse transcription polymerase chain reaction. A follow-up case–control design was then used to determine whether the functional variants are associated with CAD in 2 independent GeneID Chinese populations. Candidate pathway-based GWAS identified positive association between variants in C3AR1 and C6 and CAD. Two functional variants, rs7842 in C3AR1 and rs4400166 in C6, were found to be associated with expression levels of C3AR1 and C6, respectively. Significant association was identified between rs7842 and CAD (P=3.99×10−6; odds ratio, 1.47) and between rs4400166 and CAD (P=9.30×10−3; odds ratio, 1.24) in the validation cohort. The significant findings were confirmed in the replication cohort (P=1.53×10−5; odds ratio, 1.37 for rs7842; P=8.41×10−3; odds ratio, 1.21 for rs4400166). Conclusions—Integration of GWAS with biological pathways and expression quantitative loci is effective in identifying new risk variants for CAD. Functional variants increasing C3AR1 and C6 expression were shown to confer significant risk of CAD for the first time.