Fine mapping of a linkage peak with integration of lipid traits identifies novel coronary artery disease genes on chromosome 5

Protein glycosylation in cardiovascular health and disease

Protein glycosylation, which involves the attachment of carbohydrates to proteins, is one of the most abundant protein co-translational and post-translational modifications. Advances in technology have substantially increased our knowledge of the biosynthetic pathways involved in protein glycosylation, as well as how changes in glycosylation can affect cell function. In addition, our understanding of the role of protein glycosylation in disease processes is growing, particularly in the context of immune system function, infectious diseases, neurodegeneration and cancer. Several decades ago, cell surface glycoproteins were found to have an important role in regulating ion transport across the cardiac sarcolemma. However, with very few exceptions, our understanding of how changes in protein glycosylation influence cardiovascular (patho)physiology remains remarkably limited. Therefore, in this Review, we aim to provide an overview of N-linked and O-linked protein glycosylation, including intracellular O-linked N-acetylglucosamine protein modification. We discuss our current understanding of how all forms of protein glycosylation contribute to normal cardiovascular function and their roles in cardiovascular disease. Finally, we highlight potential gaps in our knowledge about the effects of protein glycosylation on the heart and vascular system, highlighting areas for future research.

Genetic Factors for Coronary Heart Disease and Their Mechanisms: A Meta-Analysis and Comprehensive Review of Common Variants from Genome-Wide Association Studies

Genome-wide association studies (GWAS) have discovered 163 loci related to coronary heart disease (CHD). Most GWAS have emphasized pathways related to single-nucleotide polymorphisms (SNPs) that reached genome-wide significance in their reports, while identification of CHD pathways based on the combination of all published GWAS involving various ethnicities has yet to be performed. We conducted a systematic search for articles with comprehensive GWAS data in the GWAS Catalog and PubMed, followed by a meta-analysis of the top recurring SNPs from ≥2 different articles using random or fixed-effect models according to Cochran Q and I² statistics, and pathway enrichment analysis. Meta-analyses showed significance for 265 of 309 recurring SNPs. Enrichment analysis returned 107 significant pathways, including lipoprotein and lipid metabolisms (rs7412, rs6511720, rs11591147, rs1412444, rs11172113, rs11057830, rs4299376), atherogenesis (rs7500448, rs6504218, rs3918226, rs7623687), shared cardiovascular pathways (rs72689147, rs1800449, rs7568458), diabetes-related pathways (rs200787930, rs12146487, rs6129767), hepatitis C virus infection/hepatocellular carcinoma (rs73045269/rs8108632, rs56062135, rs188378669, rs4845625, rs11838776), and miR-29b-3p pathways (rs116843064, rs11617955, rs146092501, rs11838776, rs73045269/rs8108632). In this meta-analysis, the identification of various genetic factors and their associated pathways associated with CHD denotes the complexity of the disease. This provides an opportunity for the future development of novel CHD genetic risk scores relevant to personalized and precision medicine.

Clinical significance of costimulatory molecules CD40/CD40L and CD134/CD134L in coronary heart disease: A case-control study

The aim of the study was to evaluate the potential role of CD40/CD40 ligand (CD40L) and CD134/CD134 ligand (CD134L) in the development of coronary heart disease (CHD) via the performance of a case-control study.The research objects were 234 cases of CHD patients and 120 cases of well-matched normal controls. Following the separation of peripheral blood mononuclear cells (PBMCs), real-time quantitative PCR (qRT-PCR), Western blot, immunohistochemistry, and flow cytometry were applied for the detection of mRNA levels and expression levels of CD40/CD40L and CD134/CD134L; meanwhile, intercellular adhesion molecule-1 (ICAM-1) and Fas protein mRNA levels were detected using qRT-PCR.There was no statistical difference in the comparison of baseline characteristics between groups, indicating comparability between groups. qRT-PCR and Western blot analysis indicated that CD40/CD40L and CD134/CD134L mRNA and protein expression levels were all increased in the CHD group than those in the control group. Flow cytometry further confirmed the similar tendency. Meanwhile, ICAM-1 and Fas protein mRNA levels were elevated in the CHD group and positively correlated with the above parameters. Furthermore, CD40/CD40L expression rates were negatively correlated with gender and different types of CHD. Meanwhile, CD134/CD134L expressions were also higher in male patients, in patients with family history, previous history of hypertension, diabetes, and cerebrovascular diseases.CD40/CD40L and CD134/CD134L are increased and may have potential correlation with clinical pathological features of patients with CHD. Further in-depth exploration of costimulatory molecules for CHD guidance as well as intrinsic mechanisms are needed combined with in vivo and in vitro experiments.

Association in a Chinese population of a genetic variation in the early B-cell factor 1 gene with coronary artery disease

Background

Early B-cell factor 1 (EBF1) is a transcription factor expressed primarily during early B cell development. Previous studies have shown EBF1 regulates blood glucose and lipid metabolism in mice with diabetes and central adiposity. Recently, a genetic variation (rs36071027) located in an EBF1 gene intron was associated with carotid artery intima-media thickness. However, whether this polymorphism is actually linked with coronary artery disease (CAD) and its severity remains unclear.

Methods

This study includes 293 CAD cases and 262 controls without CAD. All participants were devided into two groups based on their coronary angiography results. A polymerase chain reaction-ligase detection reaction was used to identify genotypes at rs36071027, and CAD patients were further divided into subgroups with one-, two-, or three-vessel stenosis reflective of CAD severity.

Results

The frequency of the rs36071027 TT genotype was significantly higher in CAD cases versus controls (4.8% vs. 1.5%, 95% CI: 1.13-10.81 P = 0.029). Subjects with a variant genotype T allele had an increased risk of CAD compared to C allele carriers (additive model: 95% CI: 1.13-2.23, P = 0.008). After adjustment for cardiovascular risk factors, analysis of the additive and dominant models involving rs36071027 also revealed that T allele carriers had a significantly higher risk for CAD than C allele carriers (additive model: OR 1.56, 95% CI 1.10–2.22, P = 0.013; dominant model: OR 1.60, 95% CI 1.07–2.41, P = 0.023). Furthermore, both diabetes and the CT + TT rs36071027 genotype were significantly associated with three-vessel stenosis.

Conclusion

Our results in a Chinese population suggest that the TT genotype and T alleles in rs36071027 in the EBF1 gene are associated with an increased risk of CAD and its severity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12872-017-0489-2) contains supplementary material, which is available to authorized users.

Replication of Early B-cell Factor 1 (EBF1) Gene-by-psychosocial Stress Interaction Effects on Central Adiposity in a Korean Population

Objectives

Central obesity plays a major role in the development of many chronic diseases, including cardiovascular disease and cancer. Chronic stress may be involved in the pathophysiology of central obesity. Although several large-scale genome-wide association studies have reported susceptibility genes for central adiposity, the effects of interactions between genes and psychosocial stress on central adiposity have rarely been examined. A recent study focusing on Caucasians discovered the novel gene early B-cell factor 1 (EBF1), which was associated with central obesity-related traits via interactions with stress levels. We aimed to evaluate EBF1 gene-by-stress interaction effects on central adiposity traits, including visceral adipose tissue (VAT), in Korean adults.

Methods

A total of 1467 Korean adults were included in this study. We selected 22 single-nucleotide polymorphisms (SNPs) in the EBF1 gene and analyzed their interactions with stress on central adiposity using additive, dominant, and recessive genetic modeling.

Results

The four SNPs that had strong linkage disequilibrium relationships (rs10061900, rs10070743, rs4704967, and rs10056564) demonstrated significant interactions with the waist-hip ratio in the dominant model (p_int<0.007). In addition, two other SNPs (rs6556377 and rs13180086) were associated with VAT by interactions with stress levels, especially in the recessive genetic model (p_int<0.007). As stress levels increased, the mean values of central adiposity traits according to SNP genotypes exhibited gradual but significant changes (p<0.05).

Conclusions

These results suggest that the common genetic variants for EBF1 are associated with central adiposity through interactions with stress levels, emphasizing the importance of managing stress in the prevention of central obesity.

An Efficient Stepwise Statistical Test to Identify Multiple Linked Human Genetic Variants Associated with Specific Phenotypic Traits

Recent advances in genotyping methodologies have allowed genome-wide association studies (GWAS) to accurately identify genetic variants that associate with common or pathological complex traits. Although most GWAS have focused on associations with single genetic variants, joint identification of multiple genetic variants, and how they interact, is essential for understanding the genetic architecture of complex phenotypic traits. Here, we propose an efficient stepwise method based on the Cochran-Mantel-Haenszel test (for stratified categorical data) to identify causal joint multiple genetic variants in GWAS. This method combines the CMH statistic with a stepwise procedure to detect multiple genetic variants associated with specific categorical traits, using a series of associated I × J contingency tables and a null hypothesis of no phenotype association. Through a new stratification scheme based on the sum of minor allele count criteria, we make the method more feasible for GWAS data having sample sizes of several thousands. We also examine the properties of the proposed stepwise method via simulation studies, and show that the stepwise CMH test performs better than other existing methods (e.g., logistic regression and detection of associations by Markov blanket) for identifying multiple genetic variants. Finally, we apply the proposed approach to two genomic sequencing datasets to detect linked genetic variants associated with bipolar disorder and obesity, respectively.

A Guide for a Cardiovascular Genomics Biorepository: the CATHGEN Experience

The CATHeterization GENetics (CATHGEN) biorepository was assembled in four phases. First, project start-up began in 2000. Second, between 2001 and 2010, we collected clinical data and biological samples from 9334 individuals undergoing cardiac catheterization. Samples were matched at the individual level to clinical data collected at the time of catheterization and stored in the Duke Databank for Cardiovascular Diseases (DDCD). Clinical data included the following: subject demographics (birth date, race, gender, etc.); cardiometabolic history including symptoms; coronary anatomy and cardiac function at catheterization; and fasting chemistry data. Third, as part of the DDCD regular follow-up protocol, yearly evaluations included interim information: vital status (verified via National Death Index search and supplemented by Social Security Death Index search), myocardial infarction (MI), stroke, rehospitalization, coronary revascularization procedures, medication use, and lifestyle habits including smoking. Fourth, samples were used to generate molecular data. CATHGEN offers the opportunity to discover biomarkers and explore mechanisms of cardiovascular disease.

Gene by stress genome-wide interaction analysis and path analysis identify EBF1 as a cardiovascular and metabolic risk gene

We performed gene-environment interaction genome-wide association analysis (G × E GWAS) to identify SNPs whose effects on metabolic traits are modified by chronic psychosocial stress in the Multi-Ethnic Study of Atherosclerosis (MESA). In Whites, the G × E GWAS for hip circumference identified five SNPs within the Early B-cell Factor 1 (EBF1) gene, all of which were in strong linkage disequilibrium. The gene-by-stress interaction (SNP × STRESS) term P-values were genome-wide significant (Ps = 7.14E-09 to 2.33E-08, uncorrected; Ps = 1.99E-07 to 5.18E-07, corrected for genomic control). The SNP-only (without interaction) model P-values (Ps = 0.011-0.022) were not significant at the conventional genome-wide significance level. Further analysis of related phenotypes identified gene-by-stress interaction effects for waist circumference, body mass index (BMI), fasting glucose, type II diabetes status, and common carotid intimal-medial thickness (CCIMT), supporting a proposed model of gene-by-stress interaction that connects cardiovascular disease (CVD) risk factor endophenotypes such as central obesity and increased blood glucose or diabetes to CVD itself. Structural equation path analysis suggested that the path from chronic psychosocial stress to CCIMT via hip circumference and fasting glucose was larger (estimate = 0.26, P = 0.033, 95% CI = 0.02-0.49) in the EBF1 rs4704963 CT/CC genotypes group than the same path in the TT group (estimate = 0.004, P = 0.34, 95% CI = -0.004-0.012). We replicated the association of the EBF1 SNPs and hip circumference in the Framingham Offspring Cohort (gene-by-stress term P-values = 0.007-0.012) as well as identified similar path relationships. This observed and replicated interaction between psychosocial stress and variation in the EBF1 gene may provide a biological hypothesis for the complex relationship between psychosocial stress, central obesity, diabetes, and cardiovascular disease.

Integrated microRNA-mRNA analysis of coronary artery disease

Although patients with coronary artery disease (CAD) have a high mortality rate, the pathogenesis of CAD is still poorly understood. The purpose of this study was to explore the underlying molecular mechanisms and potential target molecules for CAD. The platelet miRNA (GSE28858) and blood mRNA (GSE42148) expression profiles of patients with CAD and healthy controls were downloaded from Gene Expression Omnibus. Differentially expressed miRNAs and genes (DEGs) were identified by significant analysis of microarray algorithm after data preprocessing. Furthermore, the miRNA-target gene regulatory network was constructed based on miRecords database. The spearman correlation coefficients (ρ) between miRNAs and their target genes were calculated. Six up- (miR-340, miR-545, miR-451, miR454-5p, miR-624 and miR-585) and four down-regulated (miR-199a, miR-17-3p, miR-154 and miR-339) miRNAs were screened. Total 295 target genes of miR-545, miR-451, miR-585 and miR-154 were predicted. Among these 295 target genes, 7 genes were DEGs. Further analysis showed miR-545-TFEC and miR-585-SPOCK1 were highly positively correlated (ρ = 0.808091264; ρ = 0.874680776) in CAD samples. Therefore, differentially expressed miRNAs might participate in the pathogenesis of CAD by regulating their target genes.

Genomics of post-prandial lipidomic phenotypes in the Genetics of Lipid lowering Drugs and Diet Network (GOLDN) study

BACKGROUND

Increased postprandial lipid (PPL) response to dietary fat intake is a heritable risk factor for cardiovascular disease (CVD). Variability in postprandial lipids results from the complex interplay of dietary and genetic factors. We hypothesized that detailed lipid profiles (eg, sterols and fatty acids) may help elucidate specific genetic and dietary pathways contributing to the PPL response.

METHODS AND RESULTS

We used gas chromatography mass spectrometry to quantify the change in plasma concentration of 35 fatty acids and 11 sterols between fasting and 3.5 hours after the consumption of a high-fat meal (PPL challenge) among 40 participants from the GOLDN study. Correlations between sterols, fatty acids and clinical measures were calculated. Mixed linear regression was used to evaluate associations between lipidomic profiles and genomic markers including single nucleotide polymorphisms (SNPs) and methylation markers derived from the Affymetrix 6.0 array and the Illumina Methyl450 array, respectively. After the PPL challenge, fatty acids increased as well as sterols associated with cholesterol absorption, while sterols associated with cholesterol synthesis decreased. PPL saturated fatty acids strongly correlated with triglycerides, very low-density lipoprotein, and chylomicrons. Two SNPs (rs12247017 and rs12240292) in the sorbin and SH3 domain containing 1 (SORBS1) gene were associated with b-Sitosterol after correction for multiple testing (P≤4.5*10(-10)). SORBS1 has been linked to obesity and insulin signaling. No other markers reached the genome-wide significance threshold, yet several other biologically relevant loci are highlighted (eg, PRIC285, a co-activator of PPARa).

CONCLUSIONS

Integration of lipidomic and genomic data has the potential to identify new biomarkers of CVD risk.

A draft de novo genome assembly for the northern bobwhite (Colinus virginianus) reveals evidence for a rapid decline in effective population size beginning in the Late Pleistocene

Wild populations of northern bobwhites (Colinus virginianus; hereafter bobwhite) have declined across nearly all of their U.S. range, and despite their importance as an experimental wildlife model for ecotoxicology studies, no bobwhite draft genome assembly currently exists. Herein, we present a bobwhite draft de novo genome assembly with annotation, comparative analyses including genome-wide analyses of divergence with the chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genomes, and coalescent modeling to reconstruct the demographic history of the bobwhite for comparison to other birds currently in decline (i.e., scarlet macaw; Ara macao). More than 90% of the assembled bobwhite genome was captured within <40,000 final scaffolds (N50 = 45.4 Kb) despite evidence for approximately 3.22 heterozygous polymorphisms per Kb, and three annotation analyses produced evidence for >14,000 unique genes and proteins. Bobwhite analyses of divergence with the chicken and zebra finch genomes revealed many extremely conserved gene sequences, and evidence for lineage-specific divergence of noncoding regions. Coalescent models for reconstructing the demographic history of the bobwhite and the scarlet macaw provided evidence for population bottlenecks which were temporally coincident with human colonization of the New World, the late Pleistocene collapse of the megafauna, and the last glacial maximum. Demographic trends predicted for the bobwhite and the scarlet macaw also were concordant with how opposing natural selection strategies (i.e., skewness in the r-/K-selection continuum) would be expected to shape genome diversity and the effective population sizes in these species, which is directly relevant to future conservation efforts.

Genetic discoveries in AD using CSF amyloid and tau

The use of cerebrospinal fluid levels of Aβ42 and pTau181 as endophenotypes for genetic studies of Alzheimer's disease (AD) has led to successful identification of both rare and common AD risk variants. In addition, this approach has provided meaningful hypotheses for the biological mechanisms by which known AD risk variants modulate the disease process. In this article we discuss these successes and outline challenges to effective and continued applications of this approach. We contrast the statistical power of this approach with traditional case-control designs and discuss solutions to address challenges in quality control and data analysis for these phenotypes. Finally, we discuss the potential for the use of this approach with larger samples as well as the incorporation of next generation sequencing and for future work with other endophenotypes for AD.

The PSEN1, p.E318G variant increases the risk of Alzheimer's disease in APOE-ε4 carriers

The primary constituents of plaques (Aβ42/Aβ40) and neurofibrillary tangles (tau and phosphorylated forms of tau [ptau]) are the current leading diagnostic and prognostic cerebrospinal fluid (CSF) biomarkers for AD. In this study, we performed deep sequencing of APP, PSEN1, PSEN2, GRN, APOE and MAPT genes in individuals with extreme CSF Aβ42, tau, or ptau levels. One known pathogenic mutation (PSEN1 p.A426P), four high-risk variants for AD (APOE p.L46P, MAPT p.A152T, PSEN2 p.R62H and p.R71W) and nine novel variants were identified. Surprisingly, a coding variant in PSEN1, p.E318G (rs17125721-G) exhibited a significant association with high CSF tau (p = 9.2 × 10(-4)) and ptau (p = 1.8 × 10(-3)) levels. The association of the p.E318G variant with Aβ deposition was observed in APOE-ε4 allele carriers. Furthermore, we found that in a large case-control series (n = 5,161) individuals who are APOE-ε4 carriers and carry the p.E318G variant are at a risk of developing AD (OR = 10.7, 95% CI = 4.7-24.6) that is similar to APOE-ε4 homozygous (OR = 9.9, 95% CI = 7.2.9-13.6), and double the risk for APOE-ε4 carriers that do not carry p.E318G (OR = 3.9, 95% CI = 3.4-4.4). The p.E318G variant is present in 5.3% (n = 30) of the families from a large clinical series of LOAD families (n = 565) and exhibited a higher frequency in familial LOAD (MAF = 2.5%) than in sporadic LOAD (MAF = 1.6%) (p = 0.02). Additionally, we found that in the presence of at least one APOE-ε4 allele, p.E318G is associated with more Aβ plaques and faster cognitive decline. We demonstrate that the effect of PSEN1, p.E318G on AD susceptibility is largely dependent on an interaction with APOE-ε4 and mediated by an increased burden of Aβ deposition.

Meta-analyses of KIF6 Trp719Arg in coronary heart disease and statin therapeutic effect

Aims

The goal of our study is to assess the contribution of KIF6 Trp719Arg to both the risk of CHD and the efficacy of statin therapy in CHD patients.

Methods and Results

Meta-analysis of 8 prospective studies among 77,400 Caucasians provides evidence that 719Arg increases the risk of CHD (P<0.001, HR = 1.27, 95% CI = 1.15–1.41). However, another meta-analysis of 7 case-control studies among 65,200 individuals fails to find a significant relationship between Trp719Arg and the risk of CHD (P = 0.642, OR = 1.02, 95% CI = 0.95–1.08). This suggests that the contribution of Trp719Arg to CHD varies in different ethnic groups. Additional meta-analysis also shows that statin therapy only benefit the vascular patients carry 719Arg allele (P<0.001, relative ratio (RR) = 0.60, 95% CI = 0.54–0.67). To examine the role of this genetic variant in CHD risk in Han Chinese, we have conducted a case-control study with 289 CHD cases, 193 non-CHD controls, and 329 unrelated healthy volunteers as healthy controls. On post hoc analysis, significant allele frequency difference of 719Arg is observed between female CHD cases and female total controls under the dominant model (P = 0.04, χ² = 4.228, df = 1, odd ratio (OR) = 1.979, 95% confidence interval (CI) = 1.023–3.828). Similar trends are observed for post hoc analysis between female CHD cases and female healthy controls (dominant model: P = 0.04, χ² = 4.231, df = 1, OR = 2.015, 95% CI = 1.024–3.964). Non-genetic CHD risk factors are not controlled in these analyses.

Conclusions

Our meta-analysis demonstrates the role of Trp719Arg of KIF6 gene in the risk of CHD in Caucasians. The meta-analysis also suggests the role of this variant in statin therapeutic response in vascular diseases. Our case-control study suggests that Trp719Arg of KIF6 gene is associated with CHD in female Han Chinese through a post hoc analysis.