New variants including ARG1 polymorphisms associated with C-reactive protein levels identified by genome-wide association and pathway analysis

The Potential of Single Nucleotide Polymorphisms (SNPs) as Biomarkers and Their Association with the Increased Risk of Coronary Heart Disease: A Systematic Review

Human genetic analyses and epidemiological studies showed a potential association between several types of gene polymorphism and the development of coronary heart disease (CHD). Many studies on this pertinent topic need to be investigated further to reach an evidence-based conclusion. Therefore, in this current review, we describe several types of gene polymorphisms that are potentially linked to CHD. A systematic review using the databases EBSCO, PubMed, and ScienceDirect databases was searched until October of 2022 to find relevant studies on the topic of gene polymorphisms on risk factors for CHD, especially for the factors associated with single nucleotide polymorphisms (SNPs). The risk of bias and quality assessment was evaluated by Joanna Briggs Institute (JBI) guidelines. From keyword search results, a total of 6243 articles were identified, which were subsequently narrowed to 14 articles using prespecified inclusion criteria. The results suggested that there were 33 single nucleotide polymorphisms (SNPs) that can potentially increase the risk factors and clinical symptoms of CHD. This study also indicated that gene polymorphisms had a potential role in increasing CHD risk factors that were causally associated with atherosclerosis, increased homocysteine, immune/inflammatory response, Low-Density Lipoprotein (LDL), arterial lesions, and reduction of therapeutic effectiveness. In conclusion, the findings of this study indicate that SNPs may increase risk factors for CHD and SNPs show different effects between individuals. This demonstrates that knowledge of SNPs on CHD risk factors can be used to develop biomarkers for diagnostics and therapeutic response prediction to decide successful therapy and become the basis for defining personalized medicine in future.

Protein prediction for trait mapping in diverse populations

Genetically regulated gene expression has helped elucidate the biological mechanisms underlying complex traits. Improved high-throughput technology allows similar interrogation of the genetically regulated proteome for understanding complex trait mechanisms. Here, we used the Trans-omics for Precision Medicine (TOPMed) Multi-omics pilot study, which comprises data from Multi-Ethnic Study of Atherosclerosis (MESA), to optimize genetic predictors of the plasma proteome for genetically regulated proteome-wide association studies (PWAS) in diverse populations. We built predictive models for protein abundances using data collected in TOPMed MESA, for which we have measured 1,305 proteins by a SOMAscan assay. We compared predictive models built via elastic net regression to models integrating posterior inclusion probabilities estimated by fine-mapping SNPs prior to elastic net. In order to investigate the transferability of predictive models across ancestries, we built protein prediction models in all four of the TOPMed MESA populations, African American (n = 183), Chinese (n = 71), European (n = 416), and Hispanic/Latino (n = 301), as well as in all populations combined. As expected, fine-mapping produced more significant protein prediction models, especially in African ancestries populations, potentially increasing opportunity for discovery. When we tested our TOPMed MESA models in the independent European INTERVAL study, fine-mapping improved cross-ancestries prediction for some proteins. Using GWAS summary statistics from the Population Architecture using Genomics and Epidemiology (PAGE) study, which comprises ∼50,000 Hispanic/Latinos, African Americans, Asians, Native Hawaiians, and Native Americans, we applied S-PrediXcan to perform PWAS for 28 complex traits. The most protein-trait associations were discovered, colocalized, and replicated in large independent GWAS using proteome prediction model training populations with similar ancestries to PAGE. At current training population sample sizes, performance between baseline and fine-mapped protein prediction models in PWAS was similar, highlighting the utility of elastic net. Our predictive models in diverse populations are publicly available for use in proteome mapping methods at https://doi.org/10.5281/zenodo.4837327.

Assessing the relationship between high-sensitivity C-reactive protein and kidney function employing mendelian randomization in a Japanese community based J-MICC Study

Background

Inflammation is thought to be a risk factor for kidney disease. However, whether inflammatory status is either a cause or an outcome of chronic kidney disease remains controversial. We aimed to investigate the causal relationship between high-sensitivity C-reactive protein (hs-CRP) and estimated glomerular filtration rate (eGFR) using Mendelian randomization (MR) approaches.

Methods

A total of 10,521 participants of the Japan Multi-institutional Collaborative Cohort Study was analyzed in this study. We used two-sample MR approaches (the inverse-variance weighted (IVW), the weighted median (WM), and the MR-Egger method) to estimate the effect of genetically determined hs-CRP on kidney function. We selected four and three hs-CRP associated single nucleotide polymorphisms (SNPs) as two instrumental variables (IV): IV_CRP and IV_Asian, based on SNPs previously identified in European and Asian populations. IV_CRP and IV_Asian explained 3.4% and 3.9% of the variation in hs-CRP, respectively.

Results

Using the IV_CRP, genetically determined hs-CRP was not significantly associated with eGFR in the IVW and the WM methods (estimate per 1 unit increase in ln(hs-CRP), 0.000; 95% confidence interval [CI], −0.019 to 0.020 and −0.003; 95% CI, −0.019 to 0.014, respectively). For IV_Asian, we found similar results using the IVW and the WM methods (estimate, 0.005; 95% CI, −0.020 to 0.010 and −0.004; 95% CI, −0.020 to 0.012, respectively). The MR-Egger method also showed no causal relationships between hs-CRP and eGFR (IV_CRP: −0.008; 95% CI, −0.058 to 0.042; IV_Asian: 0.001; 95% CI, −0.036 to 0.036).

Conclusion

Our two-sample MR analyses with different IVs did not support a causal effect of hs-CRP on eGFR.

Common Variants in the ARG1 Gene Contribute to the Risk of Dilated Cardiomyopathy in the Han Chinese Population

Background: Arginase I, encoded by the ARG1 gene, is an enzyme that catalyzes the conversion of arginine to ornithine in the urea cycle; mutations in this gene has recently been reported to be associated with dilated cardiomyopathy (DCM) in Pakistan. The present study aimed to investigate the relationship between ARG1 gene mutations and DCM in the Han Chinese population. Methods: A total of 488 DCM cases and 924 matched-healthy controls were recruited. All subjects were genotyped for 12 tag single nucleotide polymorphisms (SNPs) within the ARG1 gene. Genetic association studies, including SNP and haplotype analyses, were performed. Further analyses were conducted to examine the correlations between the associated SNPs and specific clinical characteristics. Results: Only the rs2781666 and rs2781667 loci in the ARG1 gene were found to be significantly associated with DCM compared to the healthy controls. The risk of DCM at both of these loci for T allele carriers was ∼1.42-fold higher than that for carriers of the alternative alleles. There were significant differences in end-diastolic interventricular septal diameter, end-diastolic left ventricular posterior wall diameter, left ventricular end-diastolic diameter, left ventricular end-systolic diameter, and left ventricular ejection fraction among the genotype distributions of both SNPs. Furthermore, we found that the T alleles at the rs2781666 and rs2781667 loci were significantly associated with DCM in gender subgroups and the subgroup of patients <58 years of age. The haplotype T-T (rs2781666-rs2781667) also showed a significant association with DCM. Conclusion: Our results support the hypothesis that alleles and haplotypes of the ARG1 gene are significantly involved in the etiology of DCM in the Han Chinese population, but further research is necessary to elucidate the mechanism governing this association.

Association analysis between ARG1 gene polymorphisms and idiopathic dilated cardiomyopathy

The current study aimed at investigate the potential association of ARG1 polymorphisms in subjects affected by idiopathic dilated cardiomyopathy (IDCM).We have investigated 352 subjects affected by IDCM and 352 population-matched healthy controls by exploiting case-control study. The serum lipids were quantified using spectrophotometric assay, serum arginase activity was done by enzyme colorimetric assay and 2 polymorphisms (rs2781666 and rs2781667) in ARG1 were typed by polymerase chain reaction (PCR) followed by restriction fragment length polymorphism (RFLP) to find out disease associate allele/haplotype segregating in subjects affected by IDCM.Significantly high arginase activity was found to be associated with IDCM subjects when compared with population-matched healthy controls (P < .0001). The higher arginase level in IDCM subjects is negatively correlated with nitrite and nitrate (r = -0.4687, and r = -0.6435, respectively) in our study. There was a significant difference in the distribution of rs2781666 and rs2781667 genotypes of ARG1 polymorphism in patients and controls (P < .0001). Similarly, variant allele T at both loci showed a significant association with the disease phenotypes (P < .0001). Haplotype TT at rs2781666G/T and rs2781667C/T also showed a significantly association (P < .0001).To our knowledge, this is the first report to show a significant involvement of ARG1 polymorphisms to produce IDCM symptoms in subjects originating in Pakistan.

Genetic analysis of hsCRP in American Indians: The Strong Heart Family Study

Background

Increased serum levels of C-reactive protein (CRP), an important component of the innate immune response, are associated with increased risk of cardiovascular disease (CVD). Multiple single nucleotide polymorphisms (SNP) have been identified which are associated with CRP levels, and Mendelian randomization studies have shown a positive association between SNPs increasing CRP expression and risk of colon cancer (but thus far not CVD). The effects of individual genetic variants often interact with the genetic background of a population and hence we sought to resolve the genetic determinants of serum CRP in a number of American Indian populations.

Methods

The Strong Heart Family Study (SHFS) has serum CRP measurements from 2428 tribal members, recruited as large families from three regions of the United States. Microsatellite markers and MetaboChip defined SNP genotypes were incorporated into variance components, decomposition-based linkage and association analyses.

Results

CRP levels exhibited significant heritability (h² = 0.33 ± 0.05, p<1.3 X 10⁻²⁰). A locus on chromosome (chr) 6, near marker D6S281 (approximately at 169.6 Mb, GRCh38/hg38) showed suggestive linkage (LOD = 1.9) to CRP levels. No individual SNPs were found associated with CRP levels after Bonferroni adjustment for multiple testing (threshold <7.77 x 10⁻⁷), however, we found nominal associations, many of which replicate previous findings at the CRP, HNF1A and 7 other loci. In addition, we report association of 46 SNPs located at 7 novel loci on chromosomes 2, 5, 6(2 loci), 9, 10 and 17, with an average of 15.3 Kb between SNPs and all with p-values less than 7.2 X 10⁻⁴.

Conclusion

In agreement with evidence from other populations, these data show CRP serum levels are under considerable genetic influence; and include loci, such as near CRP and other genes, that replicate results from other ethnic groups. These findings also suggest possible novel loci on chr 6 and other chromosomes that warrant further investigation.

Association between CPR-related genetic variants and risk of ischemic stroke: a nested case-control study

Background: Serum C-reactive protein (CRP) has been reported to be associated with risk of ischemic vascular disease including ischemic stroke. Genome-wide association studies have revealed several gene variants related to CRP concentration. Methods: We investigated genetic variants in CRP-related genes associated with ischemic stroke in a nested case-control study with 138 ischemic stroke cases and 276 controls. We sequenced the whole coding region of six CPR-related genes and selected eligible SNPs. Three genetic models (additive, dominant and recessive) were calculated by a multivariable conditional logistic regression to estimate the association between SNPs and risk of ischemic stroke. We also calculated gene-environment interactions by using a crossover analysis. Results: Three out of 10 eligible SNPs were shown to be associated with risk of ischemic stroke. rs1800947 in CRP gene (additive model: OR = 2.08, 95% CI: 1.00-4.23) and rs1169288 in HNF1A gene (additive model: OR = 1.45, 95% CI: 1.03-2.06) were associated with an increased risk of ischemic stroke. rs440446 in APOE gene (additive model: OR = 0.63, 95%CI: 0.44-0.88) was associated with a decreased risk of ischemic stroke. Genetic risk scores models including SC-GRS and OR-GRS both showed a significant association with risk of ischemic stroke. These three SNPs interacted with smoking and red meat intake. Conclusions: Our study showed genetic variants of CRP-related genes were associated with risk of ischemic stroke. Our findings could provide useful data for the etiology of ischemic stroke.

ARG1 single nucleotide polymorphisms rs2781666 and rs2781665 confer risk of Type 2 diabetes mellitus

Genetic polymorphisms mapped in the ARG1 locus (chr6:131894344-131905472) and their functional effects on type 2 diabetes mellitus (T2DM) have not been thoroughly elucidated to date. The present study aimed to investigate an association between variant alleles at ARG1 locus and T2DM in patients. Two ARG1 single nucleotide polymorphisms (SNPs) were characterized in a representative sample of 500 patients with T2DM and 500 healthy volunteers. Serum lipid profile was studied by spectrophotometric analysis, while serum arginase-1 concentrations were determined by an enzyme-linked immunosorbent assay. The regions, encompassing target SNPs (rs2781665 and rs2781666), were amplified by polymerase chain reaction and genotypes were assigned by restriction digestions. A statistically significant increase was observed in the serum hs-CRP and arginase-1 levels in the subjects with T2DM than in controls (P <0.0001; for each). The variant genotypes of rs2781666 and rs2781665 were significantly associated with T2DM when compared with controls (P< 0.0001). Moreover, type 2 diabetic patients showed higher frequencies of T allele at rs2781666 and rs2781665 compared to the controls (OR = 1.7; 95 % CI=1.31-2.13; P <0.0001, and OR = 1.9; 95 % CI=1.45-2.38; P <0.0001, respectively). Haplotype T-T (chr6: 131893247-131893559) mapped at rs2781665-A/T and rs2781666-G/T displays higher frequency in the subjects when compared to the healthy ethnically-matched control samples (P <0.0001). We wish to propose, the first ever observation to our knowledge that concluding high levels of arginase-1 and the ARG1 polymorphisms are possible causes to confer/augment the risk of T2DM in subjects originates in Pakistan.

GWAS-Supported CRP Gene Polymorphisms and Functional Outcome of Large Artery Atherosclerotic Stroke in Han Chinese

Elevated C-reactive protein (CRP) levels increase the risk of poor functional disability in patients with ischemic stroke (IS). This study aimed to investigate the association between CRP gene polymorphisms and 3-month functional disability of large artery atherosclerotic (LAA) stroke in Han Chinese. Patients with first-ever LAA IS were prospectively enrolled in Nanjing Stroke Registry Program between August 2013 and October 2015. Five single-nucleotide polymorphisms (SNPs) (rs876537, rs2794520, rs3093059, rs7553007 and rs11265260) in CRP gene related to CRP levels in Asian by genome-wide association study were genotyped. The functional outcome at 3 months after the index stroke was assessed by the modified Rankin scale. Associations between genotypes and functional outcome of LAA IS were analyzed with logistic regression model. A total of 690 eligible patients (507 males) were evaluated. SNPs rs11265260 (multivariate-adjusted, p = 0.022), rs2794520 (multivariate-adjusted, p = 0.036) and rs3093059 (multivariate-adjusted, p = 0.027) were significantly associated with elevated CRP in acute IS. Two SNPs, rs3093059 (dominant model: adjusted OR 2.49; 95% CI 1.55-4.00; recessive model: adjusted OR 3.67; 95% CI 1.22-11.03) and rs11265260 (dominant model: adjusted OR 2.51; 95% CI 1.56-4.02; recessive model: adjusted OR 4.70; 95% CI 1.63-13.56) independently predicted 3-month poor outcome of first-ever LAA IS, after adjusting for covariates. In addition, haplotype analysis indicated that haplotype GCTGC (adjusted OR 1.76; 95% CI 1.05-2.95; p = 0.031) increased the poor outcome risk. SNPs rs3093059 and rs11265260 in CRP gene may influence the 3-month functional outcome of first-ever LAA IS in Han Chinese.

Hepatocyte nuclear factors as possible C-reactive protein transcriptional inducer in the liver and white adipose tissue of rats with experimental chronic renal failure

Inflammation related to chronic kidney disease (CKD) is an important clinical problem. We recently determined that hepatocyte nuclear factor 1α (HNF1α) was upregulated in the livers of chronic renal failure (CRF) rats—experimental model of CKD. Considering that the promoter region of gene encoding C-reactive protein (CRP) contains binding sites for HNF1α and that the loss-of-function mutation in the Hnfs1α leads to significant reduction in circulating CRP levels, we hypothesized that HNF1α can activate the Crp in CRF rats. Here, we found coordinated upregulation of genes encoding CRP, interleukin-6 (IL-6), HNF1α, and HNF4α in the livers and white adipose tissue (WAT) of CRF rats, as compared to the pair-fed and control animals. This was accompanied by elevated serum levels of CRP and IL-6. CRP and HNFs’ mRNA levels correlated positively with CRP and HNFs’ protein levels in the liver and WAT. Similar upregulation of the Crp, Il-6, and Hnfs in the liver and WAT and increased serum CRP and IL-6 concentrations were found in lipopolysaccharide (LPS)-induced systemic inflammation in rats. Moreover, silencing HNF1α in HepG2 cells by small interfering RNA led to decrease in CRP mRNA levels. Our results suggests that (a) HNFs act in concert with IL-6 in the upregulation of CRP production by the liver and WAT, leading to an increase in circulating CRP concentration in CRF rats and (b) CRF-related inflammation plays an important role in the upregulation of genes that encode HNFs and CRP in the liver and WAT of CRF rats.

CRP-level-associated polymorphism rs1205 within the CRP gene is associated with 2-hour glucose level: The SAPPHIRe study

C-reactive protein (CRP) encoded by CRP gene is a reflection of systemic inflammation. Many studies associated CRP level with diabetes and glucose levels, but the association of CRP gene with these traits is unclear. We conducted a cross-sectional study consisting of 945 siblings from 330 families collected by the Stanford Asian Pacific Program in Hypertension and Insulin Resistance (SAPPHIRe) to investigate associations between CRP polymorphisms, circulating CRP, diabetes, and glucose levels. Five single-nucleotide polymorphisms were analyzed: rs3093059, rs2794521, rs1417938, rs1800947, and rs1205. The generalized estimating equation approach was used to deal with correlated data within families. CRP level was positively correlated with diabetes prevalence and levels of fasting and 2-hour glucose (each P < 0.008). Alleles C at rs3093059 and G at rs1205 were associated with elevated CRP level (each P < 1.2 × 10⁻⁶). Allele C at rs3093059 was associated with fasting glucose (β = 0.20, P = 0.045) and G at rs1205 was associated with 2-hour glucose (β = 0.46, P = 0.00090) post oral glucose tolerance test, but only the latter passed Bonferroni correction. No polymorphism was associated with diabetes. Since 2-hour glucose is an indicator of glucose tolerance, this study indicated CRP gene is associated with glucose intolerance.

Systematic analysis of coronary artery disease datasets revealed the potential biomarker and treatment target

Coronary artery disease caused about 1 of every 7 deaths in the United States and early prevention was potential to decrease the incidence and mortality. We aimed to figure the genes involving in the coronary artery disease using meta-anlaysis. Five datasets of coronary heart disease from GEO series were retrieved and data preprocessing and quality control were carried out. Moderated t-test was used to decide the differentially expressed genes for a single dataset. And the combined p-value using systematic-analysis methods were conducted using MetaDE. The pathway enrichment was carried out using Reactome database. Protein-protein interactions of the identified differentially expressed genes were also analyzed using STRING v10.0 online tool. After removing unidentified or intermediate samples and a total of 238 cases and 189 matched or partially matched control from five microarray datasets were retrieved from GEO. Six different quality control measures were calculated and PCA biplots were plotted in order to visualize the quantitative measure. The first two PCs captured 91% of the variance and we decided to include all of the datasets for systematic analysis. Using the FDR cut-off as 0.1, nine genes, including LFNG, ID3, PLA2G7, FOLR3, PADI4, ARG1, IL1R2, NFIL3 and MGAM, were differentially expressed according to maxP. Their protein-protein interactions showed that they were closely connected and 24 Reactome pathways were related to coronary artery disease. We concluded that pathways related to immune responses, especially neutrophil degranulation, were associated with coronary heart disease.

Serum gamma-glutamyl transferase and risk of type 2 diabetes in the general Korean population: a Mendelian randomization study

Elevated gamma-glutamyl transferase (GGT) levels are associated with higher risk of type 2 diabetes in observational studies, but the underlying causal relationship is still unclear. Here, we tested a hypothesis that GGT levels have a causal effect on type 2 diabetes risk using Mendelian randomization. Data were collected from 7640 participants in a South Korean population. In a single instrumental variable (IV) analysis using two stage least squares regression with the rs4820599 in the GGT1 gene region as an instrument, one unit of GGT levels (IU/L) was associated with 11% higher risk of type 2 diabetes (odds ratio (OR) = 1.11, 95% confidence interval (CI): 1.04 to 1.19). In a multiple IV analysis using seven genetic variants that have previously been demonstrated to be associated with GGT at a genome-wide level of significance, the corresponding estimate suggested a 2.6% increase in risk (OR = 1.026, 95% CI: 1.001 to 1.052). In a two-sample Mendelian randomization analysis using genetic associations with type 2 diabetes taken from a trans-ethnic GWAS study of 110 452 independent samples, the single IV analysis confirmed an association between the rs4820599 and type 2 diabetes risk (P-value = 0.04); however, the estimate from the multiple IV analysis was compatible with the null (OR = 1.007, 95% CI: 0.993 to 1.022) with considerable heterogeneity between the causal effects estimated using different genetic variants. Overall, there is weak genetic evidence that GGT levels may have a causal role in the development of type 2 diabetes.

Pauci- and Multibacillary Leprosy: Two Distinct, Genetically Neglected Diseases

After sustained exposure to Mycobacterium leprae, only a subset of exposed individuals develops clinical leprosy. Moreover, leprosy patients show a wide spectrum of clinical manifestations that extend from the paucibacillary (PB) to the multibacillary (MB) form of the disease. This "polarization" of leprosy has long been a major focus of investigation for immunologists because of the different immune response in these two forms. But while leprosy per se has been shown to be under tight human genetic control, few epidemiological or genetic studies have focused on leprosy subtypes. Using PubMed, we collected available data in English on the epidemiology of leprosy polarization and the possible role of human genetics in its pathophysiology until September 2015. At the genetic level, we assembled a list of 28 genes from the literature that are associated with leprosy subtypes or implicated in the polarization process. Our bibliographical search revealed that improved study designs are needed to identify genes associated with leprosy polarization. Future investigations should not be restricted to a subanalysis of leprosy per se studies but should instead contrast MB to PB individuals. We show the latter approach to be the most powerful design for the identification of genetic polarization determinants. Finally, we bring to light the important resource represented by the nine-banded armadillo model, a unique animal model for leprosy.

A new susceptibility locus for myocardial infarction, hypertension, type 2 diabetes mellitus, and dyslipidemia on chromosome 12q24

We examined the role of hepatic nuclear factor-1 alpha (HNF1a) gene polymorphism on coronary artery disease (CAD) traits in 4631 Saudi angiographed individuals (2419 CAD versus 2212 controls) using TaqMan assay on ABI Prism 7900HT sequence detection system. Following adjustment for confounders, the rs2259820_CC (1.19 (1.01–1.42); P = 0.041), rs2464196_TT (1.19 (1.00–1.40); P = 0.045), and rs2259816_T (1.13 (1.01–1.26); P = 0.031) were associated with MI. The rs2259820_T (1.14 (1.03–1.26); P = 0.011) and rs2464196_C (1.12 (1.02–1.24); P = 0.024) were associated with type 2 diabetes mellitus (T2DM), while the rs2393791_T (1.14 (1.01–1.28); P = 0.032), rs7310409_G (1.16 (1.03–1.30); P = 0.013), and rs2464196_AG+GG (1.25 (1.05–1.49); P = 0.012) were implicated in hypertension. Hypertriglyceridemia was linked to the rs2393791_T (1.14 (1.02–1.27); P = 0.018), rs7310409_G (1.12 (1.01–1.25); P = 0.031), rs1169310_G (1.15 (1.04–1.28); P = 0.010), and rs1169313_CT+TT (1.24 (1.06–1.45); P = 0.008) and high low density lipoprotein-cholesterol levels were associated with rs2259820_T (1.23 (1.07–1.41); P = 0.004), rs2464196_T (1.22 (1.06–1.39); P = 0.004), and rs2259816_T (1.18 (1.02–1.36); P = 0.023). A 7-mer haplotype CATATAC (χ² = 7.50; P = 0.0062), constructed from the studied SNPs, was associated with MI, and CATATA implicated in T2DM (χ² = 3.94; P = 0.047). Hypertriglyceridemia was linked to TGCGGG (χ² = 4.26; P = 0.039), and obesity to ACGGGT (χ² = 5.04; P = 0.025). Our results suggest that the HNF1a is a common susceptibility gene for MI, T2DM, hypertension, and dyslipidemia.