Meta-analysis identifies robust association between SNP rs17465637 in MIA3 on chromosome 1q41 and coronary artery disease

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

Association study between polymorphisms in MIA3, SELE, SMAD3 and CETP genes and coronary artery disease in an Iranian population

Background

Coronary artery disease (CAD) is the most common heart disease. Several studies have shown association between some polymorphism in different genes with CAD. Finding this association can be used in order to early diagnosis and prevention of CAD.

Method

101 CAD patients with ≥ 50% luminal stenosis of any coronary vessel as case group and 111 healthy individuals as control group were selected. the polymorphisms were evaluated by ARMS-PCR and RFLP-PCR methods.

Result

The results of this study show that there is no significant association between rs17228212, rs17465637, and rs708272 and risk of CAD. But there is significant association between risk of CAD and rs5355 (p-value = 0.022) and rs3917406 (p-value = 0.006) in total cases, and rs5882 (p-value = 0.001) in male cases.

Conclusions

Our findings revealed a significant interaction between CETP SNPs and CETP activity for affecting HDL-C levels. The SELE gene is a known cell adhesion molecule with a significant role in inflammation. Studies about possible linkage between SELE gene polymorphisms and the development of CAD are conflicting. We have found a significant association between polymorphisms of SELE gene and risk of CAD.

TANGO1 marshals the early secretory pathway for cargo export

TANGO1 protein facilitates the endoplasmic reticulum (ER) export of large cargoes that cannot be accommodated in 60 nm transport vesicles. It assembles into a ring in the plane of the ER membrane to create a distinct domain. Its lumenal portion collects and sorts folded cargoes while its cytoplasmic domains collar COPII coats, recruit retrograde COPI-coated membranes that fuse within the TANGO1 ring, thus opening a tunnel for cargo transfer from the ER into a growing export conduit. This mode of cargo transfer bypasses the need for vesicular intermediates and is used to export the most abundant and bulky cargoes. The evolution of TANGO1 and its activities defines the difference between yeast and animal early secretory pathways.

GATA2 regulates the CAD susceptibility gene ADTRP rs6903956 through preferential interaction with the G allele

Myocardial infarction (MI) is a frequent outcome of coronary artery disease (CAD) and the key factor contributing to worldwide disability and death. Genetic factors contribute to the pathogenesis of CAD/MI, and SNP rs6903956 in the ADTRP gene was first found associated with CAD/MI in the Chinese Han population, which was successfully replicated in other cohorts. However, whether rs6903956 is a functional SNP and its risk mechanism to CAD/MI remains unknown. The ADTRP gene-encoded androgen-dependent TFPI regulating protein regulates vascular endothelial cell function, endothelial-monocyte adhesion, and thrombosis. The allele A of rs6903956, in particular, is associated with lower ADTRP mRNA levels in lymphocytes. In the current study, we found that SNP rs6903956 exhibits allelic differences in transcriptional activity by interacting with GATA2. Also, the A allele conferred a greater risk of CAD and MI, lowered transcriptional activity, and GATA2 binding ability as compared to the G allele. Our findings provide details on how rs6903956 regulates the expression of ADTRP and may provide novel insights into CAD pathology and susceptibility.

Feedback regulation of coronary artery disease susceptibility gene ADTRP and LDL receptors LDLR/CD36/LOX-1 in endothelia cell functions involved in atherosclerosis

A high level of low-density lipoprotein cholesterol (LDL) is one of the most important risk factors for coronary artery disease (CAD), the leading cause of death worldwide. However, a low concentration of LDL may be protective. Genome-wide association studies revealed that variation in ADTRP gene increased the risk of CAD. In this study, we found that a low concentration of oxidized-LDL induced the expression of ADTRP. Further analyses showed that knockdown of the expression of LDL receptor genes LDLR, CD36, or LOX-1 significantly downregulated ADTRP expression, whereas overexpression of LDLR/CD36/LOX-1 markedly increased ADTRP expression through the NF-κB pathway. Like ADTRP, LDLR, CD36 and LOX-1 were all involved in endothelial cell (EC) functions relevant to the initiation of atherosclerosis. Downregulation of LDLR/CD36/LOX-1 promoted monocyte adhesion to ECs and transendothelial migration of monocytes by increasing expression of ICAM-1, VCAM-1, E-selectin and P-selectin, decreased EC proliferation and migration, and increased EC apoptosis, thereby promoting the initiation of atherosclerosis. Opposite effects were observed with the overexpression of ADTRP and LDLR/CD36/LOX-1 in ECs. Interestingly, through the NF-κB and AKT pathways, overexpression of ADTRP significantly upregulated the expression of LDLR, CD36, and LOX-1, and knockdown of ADTRP expression significantly downregulated the expression of LDLR, CD36, and LOX-1. These data suggest that ADTRP and LDL receptors LDLR/CD36/LOX-1 positively regulate each other, and form a positive regulatory loop that regulates endothelial cell functions, thereby providing a potential protective mechanism against atherosclerosis. Our findings provide a new molecular mechanism by which deregulation of ADTRP and LDLR/CD36/LOX-1 promote the development of atherosclerosis and CAD.

MIA SH3 Domain ER Export Factor 3 Deficiency Prevents Neointimal Formation by Restoring BAT-Like PVAT and Decreasing VSMC Proliferation and Migration

Abnormal proliferation and migration of vascular smooth muscle cells (VSMCs) and excessive accumulation of dysfunctional PVAT are hallmarks of pathogenesis after angioplasty. Recent genome-wide association studies reveal that single-nucleotide polymorphism (SNP) in MIA3 is associated with atherosclerosis-relevant VSMC phenotypes. However, the role of MIA3 in the vascular remodeling response to injury remains unknown. Here, we found that expression of MIA3 is increased in proliferative VSMCs and knockdown of MIA3 reduces VSMCs proliferation, migration, and inflammation, whereas MIA3 overexpression promoted VSMC migration and proliferation. Moreover, knockdown of MIA3 ameliorates femoral artery wire injury-induced neointimal hyperplasia and increases brown-like perivascular adipocytes. Collectively, the data suggest that MIA3 deficiency prevents neointimal formation by decreasing VSMC proliferation, migration, and inflammation and maintaining BAT-like perivascular adipocytes in PVAT during injury-induced vascular remodeling, which provide a potential therapeutic target for preventing neointimal hyperplasia in proliferative vascular diseases.

Statistical and Functional Studies Identify Epistasis of Cardiovascular Risk Genomic Variants From Genome-Wide Association Studies

Background

Epistasis describes how gene‐gene interactions affect phenotypes, and could have a profound impact on human diseases such as coronary artery disease (CAD). The goal of this study was to identify gene‐gene interactions in CAD using an easily generalizable multi‐stage approach.

Methods and Results

Our forward genetic approach consists of multiple steps that combine statistical and functional approaches, and analyze information from global gene expression profiling, functional interactions, and genetic interactions to robustly identify gene‐gene interactions. Global gene expression profiling shows that knockdown of ANRIL (DQ485454) at 9p21.3 GWAS (genome‐wide association studies) CAD locus upregulates TMEM100 and TMEM106B. Functional studies indicate that the increased monocyte adhesion to endothelial cells and transendothelial migration of monocytes, 2 critical processes in the initiation of CAD, by ANRIL knockdown are reversed by knockdown of TMEM106B, but not of TMEM100. Furthermore, the decreased monocyte adhesion to endothelial cells and transendothelial migration of monocytes induced by ANRIL overexpression was reversed by overexpressing TMEM106B. TMEM106B expression was upregulated by >2‐fold in CAD coronary arteries. A significant association was found between variants in TMEM106B (but not in TMEM100) and CAD (P=1.9×10⁻⁸). Significant gene‐gene interaction was detected between ANRIL variant rs2383207 and TMEM106B variant rs3807865 (P=0.009). A similar approach also identifies significant interaction between rs6903956 in ADTRP and rs17465637 in MIA3 (P=0.005).

Conclusions

We demonstrate 2 pairs of epistatic interactions between GWAS loci for CAD and offer important insights into the genetic architecture and molecular mechanisms for the pathogenesis of CAD. Our strategy has broad applicability to the identification of epistasis in other human diseases.

The Mediterranean diet reduces the genetic risk of chromosome 9p21 for myocardial infarction in an Asian population community cohort

The interaction of genetic susceptibility and dietary habits in cardiovascular disease (CVD) remains undetermined. The purpose of this study was to investigate whether a Mediterranean dietary style modified the genetic risk of developing CVD in a Chinese cohort. A total of 2098 subjects with dietary information from a Chinese community cohort (CVDFACTS) were enrolled. Candidate genes, including SNP markers rs1333049 (CDKN2B, 9p21.3), rs17465637 (MIA3, 1q41) and rs501120 (CXCL12, 10q11.21), were genotyped to analyze the association with future CVD. The impact of dietary pattern was also analyzed according to adherence to the diet using the Mediterranean Diet Score (MDS). After an average follow-up of 7.8 years, only the C risk allele of rs1333049 at chromosome 9p21.3 was associated with a higher risk of MI with either an additive [HR = 1.78, 95% CI:1.23-2.5] or a recessive model [HR = 2.40, 95% CI: 1.42-4.04], and the CC genotype had a higher risk of developing MI (p = 0.009, log-rank test). There was no significant difference in the association of the lipid profile with future CV outcomes among the MDS tertiles. However, the high MI risk of the CC genotype in individuals consuming a less healthy diet (MDS1) (HR: 6.39, 95% CI: 1.74-23.43) significantly decreased to 2.38 (95% CI: 0.57-10.04) in individuals consuming a healthier diet (MDS3), indicating that a healthier dietary pattern (higher MDS) modified the risk of developing MI in carriers of variants in CDKN2B. In conclusion, genetic variants of CDKN2B at 9p21 were significantly associated with future MI risk in a Chinese cohort, and the genetic risk of MI could be modified by a healthier diet.

Identifying genetic markers associated with susceptibility to cardiovascular diseases

The development of cardiovascular diseases (CVDs) is due to a complex interaction between the genome and the environment. Understanding how genetic differences in individuals contribute to their susceptibility to CVDs can help guide practitioners to give the best advice to achieve a favorable outcome for the patient. As genome technologies evolve, genotyping of individuals could be available to all patients using a simple saliva test. Large-scale genome-wide association studies and meta analyses have provided powerful insights into polymorphisms that may be predictive of disease and an individual's response to certain nutrients, but moving forward it is imperative that these insights can be applied in the medical setting to reduce the incidence and mortality of CVDs.

Cardiovascular diseases (CVDs) are the leading cause of death worldwide, and while most CVDs can be prevented by adopting a healthy lifestyle, this is only half the story. Evidence suggests changes in an individual's genes or DNA can cause some form of CVDs, highlighting a complex relationship between genes and the environment. Genotyping, a process used to determine genetic differences within an individual's DNA, can provide doctors with relevant information to identify individuals who are at high risk of developing CVDs. This would allow treatment to begin early and encourage individuals to adopt a healthy lifestyle to reduce their risk.

GWAS implicated risk variants in different genes contribute additively to increase the risk of coronary artery disease (CAD) in the Pakistani subjects

Background

Coronary artery disease (CAD) remains the single most important cause of mortality worldwide. Many candidate and GWAS genetic variants have been identified in the recent years. In the current study, we selected six SNPs from various genes that have originally been identified in GWAS studies and examined the association of SNPs individually and as a genetic risk score (GRS) with CAD and blood lipid levels in the Pakistani subjects.

Methods

Six hundred twenty-four (404 cases and 219 controls) subjects were genotyped for variants rs10757274 in CDKN2A gene, rs17465637 in MIA3 gene, rs7025486 in DAB2IP gene, rs17228212 in SMAD3 gene, rs981887 in MRAS gene and rs1746048 in CXCL12 gene, by TaqMan and KASPar allele discrimination techniques. Serum lipid parameters were measured using commercially available kits. Statistical analyses were done using SPSS version 22.

Results

Individually, the single SNPs were not associated with CAD (p < 0.05). However, the combined GRS of 6 SNPs was significantly higher in cases than controls (4.89 ± 0.11 vs 4.58 ± 0.08, p = 0.024). Among blood lipids, GRS showed significant positive association with serum triglycerides levels (p = 0.022).

Conclusion

The GRS was quantitatively associated with CAD risk and showed association with serum triglycerides levels, suggesting that the mechanism of these variants is likely to be in part at least through creating an atherogenic lipid profile in subjects carrying high numbers of risk alleles.

Electronic supplementary material

The online version of this article (10.1186/s12944-018-0736-2) contains supplementary material, which is available to authorized users.

The importance of melanoma inhibitory activity gene family in the tumor progression of oral cancer

Oral squamous cell carcinoma has a high potential for locoregional invasion and nodal metastasis. Consequently, early detection of such malignancies is of immense importance. The melanoma inhibitory activity (MIA) gene family comprises MIA, MIA2, transport and Golgi organization protein 1 (TANGO), and otoraplin (OTOR). These members of the MIA gene family have a highly conserved Src homology 3 (SH3)-like structure. Although the molecules of this family share 34-45% amino acid homology and 47-59% cDNA sequence homology, those members, excluding OTOR, play different tumor-associated functions. MIA has a pivotal role in the progression and metastasis of melanoma; MIA2 and TANGO have been suggested to possess tumor-suppressive functions; and OTOR is uniquely expressed in cochlea of the inner ear. Therefore, the definite functions of the MIA gene family in cancer cells remain unclear. Since the members of the MIA gene family are secreted proteins, these molecules might be useful tumor markers that can be detected in the body fluids, including serum and saliva. In this review, we described the molecular biological functions of the MIA gene family in oral cancer.

Significant genetic association of a functional TFPI variant with circulating fibrinogen levels and coronary artery disease

The tissue factor pathway inhibitor (TFPI) gene encodes a protease inhibitor with a critical role in regulation of blood coagulation. Some genomic variants in TFPI were previously associated with plasma TFPI levels, however, it remains to be further determined whether TFPI variants are associated with other coagulation factors. In this study, we carried out a large population-based study with 2313 study subjects for blood coagulation data, including fibrinogen levels, prothrombin time (PT), activated partial thromboplastin time (APTT), and thrombin time (TT). We identified significant association of TFPI variant rs10931292 (a functional promoter variant with reduced transactivation) with increased plasma fibrinogen levels (P = 0.017 under a recessive model), but not with PT, APTT or TT (P > 0.05). Using a large case-control association study population with 4479 CAD patients and 3628 controls, we identified significant association between rs10931292 and CAD under a recessive model (OR 1.23, P = 0.005). For the first time, we show that a TFPI variant is significantly associated with fibrinogen levels and risk of CAD. Our finding contributes significantly to the elucidation of the genetic basis and biological pathways responsible for fibrinogen levels and development of CAD.

Androgen inhibits key atherosclerotic processes by directly activating ADTRP transcription

Low androgen levels are associated with an increased risk of coronary artery disease (CAD), thrombosis and myocardial infarction (MI), suggesting that androgen has a protective role. However, little is known about the underlying molecular mechanism. Our genome-wide association study identified the ADTRP gene encoding the androgen-dependent TFPI regulating protein as a susceptibility gene for CAD and MI. The expression level of ADTRP was regulated by androgen, but the molecular mechanism is unknown. In this study, we identified the molecular mechanism by which androgen regulates ADTRP expression and tested the hypothesis that androgen plays a protective role in cardiovascular disease by activating ADTRP expression. Luciferase assays with an ADTRP promoter luciferase reporter revealed that androgen regulated ADTRP transcription in a dose- and time-dependent manner, and the effect was abolished by three different androgen inhibitors, including pyrvinium pamoate, bicalutamide, and cyproterone acetate. Chromatin-immunoprecipitation showed that the androgen receptor bound to a half androgen response element (ARE, TGTTCT) located at +324bp from the ADTRP transcription start site. The ARE is required for concentration-dependent transcriptional activation of ADTRP. HL-60 monocyte adhesion to EAhy926 endothelial cells (ECs) and transmigration across the EC layer, the two processes critical to development of CAD and MI, were inhibited by androgen, but the effect was rescued by ADTRP siRNA and exacerbated by overexpression of ADTRP and its downstream genes PIK3R3 and MIA3. These data suggest that one molecular mechanism by which androgen confers protection against CAD is stimulation of ADTRP expression.

Identification of a molecular signaling gene-gene regulatory network between GWAS susceptibility genes ADTRP and MIA3/TANGO1 for coronary artery disease

Coronary artery disease (CAD) is the leading cause of death worldwide. GWAS have identified >50 genomic loci for CAD, including ADTRP and MIA3/TANGO1. However, it is important to determine whether the GWAS genes form a molecular network. In this study, we have uncovered a novel molecular network between ADTRP and MIA3/TANGO1 for the pathogenesis of CAD. We showed that knockdown of ADTRP expression markedly down-regulated expression of MIA3/TANGO1. Mechanistically, ADTRP positively regulates expression of PIK3R3 encoding the regulatory subunit 3 of PI3K, which leads to activation of AKT, resulting in up-regulation of MIA3/TANGO1. Both ADTRP and MIA3/TANGO1 are involved in endothelial cell (EC) functions relevant to atherosclerosis. Knockdown of ADTRP expression by siRNA promoted oxidized-LDL-mediated monocyte adhesion to ECs and transendothelial migration of monocytes, inhibited EC proliferation and migration, and increased apoptosis, which was reversed by expression of constitutively active AKT1 and MIA3/TANGO1 overexpression, while the over-expression of ADTRP in ECs blunted these processes. Knockdown of MIA3/TANGO1 expression also promoted monocyte adhesion to ECs and transendothelial migration of monocytes, and vice versa for overexpression of MIA3/TANGO1. We found that ADTRP negatively regulates the levels of collagen VII and ApoB in HepG2 and endothelial cells, which are downstream regulatory targets of MIA3/TANGOI. In conclusion, we have uncovered a novel molecular signaling pathway for the pathogenesis of CAD, which involves a novel gene-gene regulatory network. We show that ADTRP positively regulates PIK3R3 expression, which leads to activation of AKT and up-regulation of MIA3/TANGO1, thereby regulating endothelial cell functions directly relevant to atherosclerosis.

Coronary artery disease susceptibility gene ADTRP regulates cell cycle progression, proliferation, and apoptosis by global gene expression regulation

The ADTRP gene encodes the androgen-dependent TFPI-regulating protein and is a susceptibility gene for contrary artery disease (CAD). We performed global gene expression profiling for ADTRP knock-down using microarrays in human HepG2 cells. Follow-up real-time RT-PCR analysis demonstrated that ADTRP knock-down regulates a diverse set of genes, including upregulation of seven histone genes, downregulation of multiple cell cycle genes (CCND1, CDK4, and CDKN1A), and upregulation of apoptosis genes (CASP7 and PDCD2) in HepG2 cells and endothelial cells. Consistently, ADTRP increases the number of S phase cells during cell cycle, promotes cell proliferation, and inhibits apoptosis. Our study provides novel insights into the function of ADTRP and biological pathways involving ADTRP, which may be involved in the pathogenesis of CAD.

TANGO1 and Mia2/cTAGE5 (TALI) cooperate to export bulky pre-chylomicrons/VLDLs from the endoplasmic reticulum

Santos et al. show that TANGO1 and a TANGO1-like protein, TALI, bind each other and function together as receptors to export bulky ApoB-containing lipid particles from the endoplasmic reticulum. However, TANGO1-mediated export of bulky collagens by the same cells is TALI independent.

Procollagens, pre-chylomicrons, and pre–very low-density lipoproteins (pre-VLDLs) are too big to fit into conventional COPII-coated vesicles, so how are these bulky cargoes exported from the endoplasmic reticulum (ER)? We have shown that TANGO1 located at the ER exit site is necessary for procollagen export. We report a role for TANGO1 and TANGO1-like (TALI), a chimeric protein resulting from fusion of MIA2 and cTAGE5 gene products, in the export of pre-chylomicrons and pre-VLDLs from the ER. TANGO1 binds TALI, and both interact with apolipoprotein B (ApoB) and are necessary for the recruitment of ApoB-containing lipid particles to ER exit sites for their subsequent export. Although export of ApoB requires the function of both TANGO1 and TALI, the export of procollagen XII by the same cells requires only TANGO1. These findings reveal a general role for TANGO1 in the export of bulky cargoes from the ER and identify a specific requirement for TALI in assisting TANGO1 to export bulky lipid particles.

Genomic variant in CAV1 increases susceptibility to coronary artery disease and myocardial infarction

BACKGROUND

The CAV1 gene encodes caveolin-1 expressed in cell types relevant to atherosclerosis. Cav-1-null mice showed a protective effect on atherosclerosis under the ApoE(-/-) background. However, it is unknown whether CAV1 is linked to CAD and MI in humans. In this study we analyzed a tagSNP for CAV1 in intron 2, rs3807989, for potential association with CAD.

METHODS AND RESULTS

We performed case-control association studies in three independent Chinese Han populations from GeneID, including 1249 CAD cases and 841 controls in Population I, 1260 cases and 833 controls in Population II and 790 cases and 1212 controls in Population III (a total of 3299 cases and 2886 controls). We identified significant association between rs3807989 and CAD in three independent populations and in the combined population (Padj = 2.18 × 10(-5), OR = 1.19 for minor allele A). We also detected significant association between rs3807989 and MI (Padj = 5.43 × 10(-5), OR = 1.23 for allele A). Allele A of SNP rs3807989 was also associated with a decreased level of LDL cholesterol. Although rs3807989 is a tagSNP for both CAV1 and nearby CAV2, allele A of SNP rs3807989 was associated with an increased expression level of CAV1 (both mRNA and protein), but not CAV2.

CONCLUSIONS

The data in this study demonstrated that rs3807989 at the CAV1/CAV2 locus was associated with significant risk of CAD and MI by increasing expression of CAV1 (but not CAV2). Thus, CAV1 becomes a strong candidate susceptibility gene for CAD/MI in humans.

Potential Signals of Natural Selection in the Top Risk Loci for Coronary Artery Disease: 9p21 and 10q11

BACKGROUND

Coronary artery disease (CAD) is a complex disease and the leading cause of death in the world. Populations of different ancestry do not always share the same risk markers. Natural selective processes may be the cause of some of the population differences detected for specific risk mutations.

OBJECTIVE

In this study, 384 single nucleotide polymorphisms (SNPs) located in four genomic regions associated with CAD (1p13, 1q41, 9p21 and 10q11) are analysed in a set of 19 populations from Europe, Middle East and North Africa and also in Asian and African samples from the 1000 Genomes Project. The aim of this survey is to explore for the first time whether the genetic variability in these genomic regions is better explained by demography or by natural selection.

RESULTS

The results indicate significant differences in the structure of genetic variation and in the LD patterns among populations that probably explain the population disparities found in markers of susceptibility to CAD.

CONCLUSIONS

The results are consistent with potential signature of positive selection in the 9p21 region and of balancing selection in the 9p21 and 10q11. Specifically, in Europe three CAD risk markers in the 9p21 region (rs9632884, rs1537371 and rs1333042) show consistent signals of positive selection. The results of this study are consistent with a potential selective role of CAD in the configuration of genetic diversity in current human populations.

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.

Significant Association Between CAV1 Variant rs3807989 on 7p31 and Atrial Fibrillation in a Chinese Han Population

Background

Recent genome-wide association studies (GWAS) in European ancestry populations revealed several genomic loci for atrial fibrillation (AF). We previously replicated the 4q25 locus (PITX2) and 16q22 locus (ZFHX3) in the Chinese population, but not the KCNN3 locus on 1q21. With single-nucleotide polymorphism rs3807989 in CAV1 encoding caveolin-1, however, controversial results were reported in 2 Chinese replication studies.

Methods and Results

Six remaining AF genetic loci from GWAS, including rs3807989/CAV1, rs593479/PRRX1, rs6479562/C9orf3, rs10824026/SYNPO2L, rs1152591/SYNE2, and rs7164883/HCN4, were analyzed in a Chinese Han population with 941 cases and 562 controls. Only rs3807989 showed significant association with AF (P_adj=4.77×10⁻⁵), and the finding was replicated in 2 other independent populations with 709 cases and 2175 controls, 463 cases and 644 controls, and the combined population with a total of 2113 cases and 3381 controls (P_adj=2.20×10⁻⁹; odds ratio [OR]=1.34 for major allele G). Meta-analysis, together with data from previous reports in Chinese and Japanese populations, also showed a significant association between rs3807989 and AF (P=3.40×10⁻⁴; OR=1.24 for allele G). We also found that rs3807989 showed a significant association with lone AF in 3 independent populations and in the combined population (P_adj=3.85×10⁻⁸; OR=1.43 for major allele G).

Conclusions

The data in this study revealed a significant association between rs3807989 and AF in the Chinese Han population. Together with the findings that caveolin-1 interacts with potassium channels Kir2.1, KCNH2, and HCN4 and sodium channels Nav1.5 and Nav1.8, CAV1 becomes a strong candidate susceptibility gene for AF across different ethnic populations. This study is the first to show a significant association between rs3807989 and lone AF.

Time series analyses of hand, foot and mouth disease integrating weather variables

Background

The past decade witnessed an increment in the incidence of hand foot mouth disease (HFMD) in the Pacific Asian region; specifically, in Guangzhou China. This emphasized the requirement of an early warning system designed to allow the medical community to better prepare for outbreaks and thus minimize the number of fatalities.

Methods

Samples from 1,556 inpatients (hospitalized) and 11,004 outpatients (non-admitted) diagnosed with HFMD were collected in this study from January 2009 to October 2013. Seasonal Autoregressive Integrated Moving Average (SARIMA) model was applied to establish high predictive model for inpatients and outpatient as well as three viral serotypes (EV71, Pan-EV and CA16). To integrate climate variables in the data analyses, data from eight climate variables were simultaneously obtained during this period. Significant climate variable identified by correlation analyses was executed to improve time series modeling as external repressors.

Results

Among inpatients with HFMD, 248 (15.9%) were affected by EV71, 137 (8.8%) were affected by Pan-EV+, and 436 (28.0%) were affected by CA16. Optimal Univariate SARIMA model was identified: (2,0,3)(1,0,0)₅₂ for inpatients, (0,1,0)(0,0,2)₅₂ for outpatients as well as three serotypes (EV71, (1,0,1)(0,0,1)₅₂; CA16, (1,0,1)(0,0,0)₅₂; Pan-EV, (1,0,1)(0,0,0)₅₂). Using climate as our independent variable, precipitation (PP) was first identified to be associated with inpatients (r = 0.211, P = 0.001), CA16-serotype (r = 0.171, P = 0.007) and outpatients (r = 0.214, P = 0.01) in partial correlation analyses, and was then shown a significant lag in cross-autocorrelation analyses. However, inclusion of PP [lag -3 week] as external repressor showed a moderate impact on the predictive performance of the SARIMA model described here-in.

Conclusion

Climate patterns and HFMD incidences have been shown to be strongly correlated. The SARIMA model developed here can be a helpful tool in developing an early warning system for HFMD.