Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79x10^-13), rs74506613 (JMJD1C, P = 1.17x10^-19), rs4782371 (ZFPM1, P = 1.59x10^-9) and rs2639990 (ZADH2, P = 1.72x10^-8), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52x10^-18; rs7043199, VLDLR-AS1, P = 5.12x10^-14) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39x10^-1467; rs1740073, C6orf223, P = 2.34x10^-17; rs6993770, ZFPM2, P = 2.44x10^-60; rs2375981, KCNV2, P = 1.48x10^-100). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.     

New Sequence Variants in HLA Class II/III Region Associated with Susceptibility to Knee Osteoarthritis Identified by Genome-Wide Association Study

Osteoarthritis (OA) is a common disease that has a definite genetic component. Only a few OA susceptibility genes that have definite functional evidence and replication of association have been reported, however. Through a genome-wide association study and a replication using a total of ∼4,800 Japanese subjects, we identified two single nucleotide polymorphisms (SNPs) (rs7775228 and rs10947262) associated with susceptibility to knee OA. The two SNPs were in a region containing HLA class II/III genes and their association reached genome-wide significance (combined P = 2.43×10⁻⁸ for rs7775228 and 6.73×10⁻⁸ for rs10947262). Our results suggest that immunologic mechanism is implicated in the etiology of OA.     

Genome-Wide Association Scan for Variants Associated with Early-Onset Prostate Cancer

Prostate cancer is the most common non-skin cancer and the second leading cause of cancer related mortality for men in the United States. There is strong empirical and epidemiological evidence supporting a stronger role of genetics in early-onset prostate cancer. We performed a genome-wide association scan for early-onset prostate cancer. Novel aspects of this study include the focus on early-onset disease (defined as men with prostate cancer diagnosed before age 56 years) and use of publically available control genotype data from previous genome-wide association studies. We found genome-wide significant (p<5×10⁻⁸) evidence for variants at 8q24 and 11p15 and strong supportive evidence for a number of previously reported loci. We found little evidence for individual or systematic inflated association findings resulting from using public controls, demonstrating the utility of using public control data in large-scale genetic association studies of common variants. Taken together, these results demonstrate the importance of established common genetic variants for early-onset prostate cancer and the power of including early-onset prostate cancer cases in genetic association studies.     

Genome-Wide Association Study Identifies Variants in PMS1 Associated with Serum Ferritin in a Chinese Population

Only a small proportion of genetic variation in serum ferritin has been explained by variant genetic studies, and genome-wide association study (GWAS) for serum ferritin has not been investigated widely in Chinese population. We aimed at exploring the novel genetic susceptibility to serum ferritin, and performed this two stage GWAS in a healthy Chinese population of 3,495 men aged 20–69 y, including 1,999 unrelated subjects in the first stage and 1,496 independent individuals in the second stage. Serum ferritin was measured with electrochemiluminescence immunoassay, and DNA samples were collected for genotyping. A total of 1,940,243 SNPs were tested by using multivariate linear regression analysis. After adjusting for population stratification, age and BMI, the rs5742933 located in the 5′UTR region of PMS1 gene on chromosome 2 was the most significantly associated with ferritin concentrations (P-combined = 2.329×10⁻¹⁰) (β = −0.11, 95% CI: −0.14, −0.07). Moreover, this marker was about 200kb away from the candidate gene SLC40A1 which is responsible for iron export. PMS1 gene was the novel genetic susceptibility to serum ferritin in Chinese males and its relation to SLC40A1 needs further study.     

Genome-Wide Association Study Reveals Two Nucleotide Variants Associated with Educational Attainment in Koreans

Russian Journal of Genetics - Genome-wide association studies for educational attainment have been limited to Europeans, but genetic factors might be largely related to cultural environments and...     

Genetic Variants Associated with Increased Risk of Malignant Pleural Mesothelioma: A Genome-Wide Association Study

Asbestos exposure is the main risk factor for malignant pleural mesothelioma (MPM), a rare aggressive tumor. Nevertheless, only 5–17% of those exposed to asbestos develop MPM, suggesting the involvement of other environmental and genetic risk factors.To identify the genetic risk factors that may contribute to the development of MPM, we conducted a genome-wide association study (GWAS; 370,000 genotyped SNPs, 5 million imputed SNPs) in Italy, among 407 MPM cases and 389 controls with a complete history of asbestos exposure. A replication study was also undertaken and included 428 MPM cases and 1269 controls from Australia.Although no single marker reached the genome-wide significance threshold, several associations were supported by haplotype-, chromosomal region-, gene- and gene-ontology process-based analyses. Most of these SNPs were located in regions reported to harbor aberrant alterations in mesothelioma (SLC7A14, THRB, CEBP350, ADAMTS2, ETV1, PVT1 and MMP14 genes), causing at most a 2–3-fold increase in MPM risk. The Australian replication study showed significant associations in five of these chromosomal regions (3q26.2, 4q32.1, 7p22.2, 14q11.2, 15q14).Multivariate analysis suggested an independent contribution of 10 genetic variants, with an Area Under the ROC Curve (AUC) of 0.76 when only exposure and covariates were included in the model, and of 0.86 when the genetic component was also included, with a substantial increase of asbestos exposure risk estimation (odds ratio, OR: 45.28, 95% confidence interval, CI: 21.52–95.28).These results showed that genetic risk factors may play an additional role in the development of MPM, and that these should be taken into account to better estimate individual MPM risk in individuals who have been exposed to asbestos.     

Gene-Specific DNA Methylation Association with Serum Levels of C-Reactive Protein in African Americans

A more thorough understanding of the differences in DNA methylation (DNAm) profiles in populations may hold promise for identifying molecular mechanisms through which genetic and environmental factors jointly contribute to human diseases. Inflammation is a key molecular mechanism underlying several chronic diseases including cardiovascular disease, and it affects DNAm profile on both global and locus-specific levels. To understand the impact of inflammation on the DNAm of the human genome, we investigated DNAm profiles of peripheral blood leukocytes from 966 African American participants in the Genetic Epidemiology Network of Arteriopathy (GENOA) study. By testing the association of DNAm sites on CpG islands of over 14,000 genes with C-reactive protein (CRP), an inflammatory biomarker of cardiovascular disease, we identified 257 DNAm sites in 240 genes significantly associated with serum levels of CRP adjusted for age, sex, body mass index and smoking status, and corrected for multiple testing. Of the significantly associated DNAm sites, 80.5% were hypomethylated with higher CRP levels. The most significant Gene Ontology terms enriched in the genes associated with the CRP levels were immune system process, immune response, defense response, response to stimulus, and response to stress, which are all linked to the functions of leukocytes. While the CRP-associated DNAm may be cell-type specific, understanding the DNAm association with CRP in peripheral blood leukocytes of multi-ethnic populations can assist in unveiling the molecular mechanism of how the process of inflammation affects the risks of developing common disease through epigenetic modifications.     

Non-Homologous End-Joining Pathway Associated with Occurrence of Myocardial Infarction: Gene Set Analysis of Genome-Wide Association Study Data

Purpose

DNA repair deficiencies have been postulated to play a role in the development and progression of cardiovascular disease (CVD). The hypothesis is that DNA damage accumulating with age may induce cell death, which promotes formation of unstable plaques. Defects in DNA repair mechanisms may therefore increase the risk of CVD events. We examined whether the joints effect of common genetic variants in 5 DNA repair pathways may influence the risk of CVD events.

Methods

The PLINK set-based test was used to examine the association to myocardial infarction (MI) of the DNA repair pathway in GWAS data of 866 subjects of the GENetic DEterminants of Restenosis (GENDER) study and 5,244 subjects of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) study. We included the main DNA repair pathways (base excision repair, nucleotide excision repair, mismatch repair, homologous recombination and non-homologous end-joining (NHEJ)) in the analysis.

Results

The NHEJ pathway was associated with the occurrence of MI in both GENDER (P = 0.0083) and PROSPER (P = 0.014). This association was mainly driven by genetic variation in the MRE11A gene (P_GENDER = 0.0001 and P_PROSPER = 0.002). The homologous recombination pathway was associated with MI in GENDER only (P = 0.011), for the other pathways no associations were observed.

Conclusion

This is the first study analyzing the joint effect of common genetic variation in DNA repair pathways and the risk of CVD events, demonstrating an association between the NHEJ pathway and MI in 2 different cohorts.     

Gene-Based Association Analysis Identified Novel Genes Associated with Bone Mineral Density

Genetic factors contribute to the variation of bone mineral density (BMD), which is a major risk factor of osteoporosis. The aim of this study was to identify more “novel” genes for BMD. Based on the publicly available SNP-based P values, we performed an initial gene-based analysis in a total of 32,961 individuals. Furthermore, we performed differential expression, pathway and protein-protein interaction analyses to find supplementary evidence to support the significance of the identified genes. About 21,695 genes for femoral neck (FN)-BMD and 21,683 genes for lumbar spine (LS)-BMD were analyzed using gene-based association analysis. A total of 35 FN-BMD associated genes and 53 LS-BMD associated genes were identified (P < 2.3×10^-6) after Bonferroni correction. Among them, 64 genes have not been reported in previous SNP-based genome-wide association studies. Differential expression analysis further supported the significant associations of 14 genes with FN-BMD and 19 genes with LS-BMD. Especially, WNT3 and WNT9B in the Wnt signaling pathway for FN-BMD were further supported by pathway analysis and protein-protein interaction analysis. The present study took the advantage of gene-based association method to perform a supplementary analysis of the GWAS dataset and found some BMD-associated genes. The evidence taken together supported the importance of Wnt signaling pathway genes in determining osteoporosis. Our findings provided more insights into the genetic basis of osteoporosis.     

Validation of Type 2 Diabetes Risk Variants Identified by Genome-Wide Association Studies in Han Chinese Population: A Replication Study and Meta-Analysis

Background

Several genome-wide association studies (GWAS) involving European populations have successfully identified risk genetic variants associated with type 2 diabetes mellitus (T2DM). However, the effects conferred by these variants in Han Chinese population have not yet been fully elucidated.

Methods

We analyzed the effects of 24 risk genetic variants with reported associations from European GWAS in 3,040 Han Chinese subjects in Taiwan (including 1,520 T2DM cases and 1,520 controls). The discriminative power of the prediction models with and without genotype scores was compared. We further meta-analyzed the association of these variants with T2DM by pooling all candidate-gene association studies conducted in Han Chinese.

Results

Five risk variants in IGF2BP2 (rs4402960, rs1470579), CDKAL1 (rs10946398), SLC30A8 (rs13266634), and HHEX (rs1111875) genes were nominally associated with T2DM in our samples. The odds ratio was 2.22 (95% confidence interval, 1.81-2.73, P<0.0001) for subjects with the highest genetic score quartile (score>34) as compared with subjects with the lowest quartile (score<29). The incoporation of genotype score into the predictive model increased the C-statistics from 0.627 to 0.657 (P<0.0001). These estimates are very close to those observed in European populations. Gene-environment interaction analysis showed a significant interaction between rs13266634 in SLC30A8 gene and age on T2DM risk (P<0.0001). Further meta-analysis pooling 20 studies in Han Chinese confirmed the association of 10 genetic variants in IGF2BP2, CDKAL1, JAZF1, SCL30A8, HHEX, TCF7L2, EXT2, and FTO genes with T2DM. The effect sizes conferred by these risk variants in Han Chinese were similar to those observed in Europeans but the allele frequencies differ substantially between two populations.

Conclusion

We confirmed the association of 10 variants identified by European GWAS with T2DM in Han Chinese population. The incorporation of genotype scores into the prediction model led to a small but significant improvement in T2DM prediction.     

Two Variants of the C-Reactive Protein Gene Are Associated with Risk of Pre-Eclampsia in an American Indian Population

Background

The etiology of pre-eclampsia (PE) is unknown; but it is accepted that normal pregnancy represents a distinctive challenge to the maternal immune system. C-reactive protein is a prominent component of the innate immune system; and we previously reported an association between PE and the CRP polymorphism, rs1205. Our aim was to explore the effects of additional CRP variants. The IBC (Cardiochip) genotyping microarray focuses on candidate genes and pathways related to the pathophysiology of cardiovascular disease.

Methods

This study recruited 140 cases of PE and 270 matched controls, of which 95 cases met criteria as severe PE, from an American Indian community. IBC array genotypes from 10 suitable CRP SNPs were analyzed. A replication sample of 178 cases and 427 controls of European ancestry was also genotyped.

Results

A nominally significant difference (p value <0.05) was seen in the distribution of discordant matched pairs for rs3093068; and Bonferroni corrected differences (P<0.005) were seen for rs876538, rs2794521, and rs3091244. Univariate conditional logistic regression odds ratios (OR) were nominally significant for rs3093068 and rs876538 models only. Multivariate logistic models with adjustment for mother''s age, nulliparity and BMI attenuated the effect (OR 1.58, P = 0.066, 95% CI 0.97–2.58) for rs876538 and (OR 2.59, P = 0.050, 95% CI 1.00–6.68) for rs3093068. An additive risk score of the above two risk genotypes shows a multivariate adjusted OR of 2.04 (P = 0.013, 95% CI 1.16–3.56). The replication sample also demonstrated significant association between PE and the rs876538 allele (OR = 1.55, P = 0.01, 95% CI 2.16–1.10). We also show putative functionality for the rs876538 and rs3093068 CRP variants.

Conclusion

The CRP variants, rs876538 and rs3093068, previously associated with other cardiovascular disease phenotypes, show suggestive association with PE in this American Indian population, further supporting a possible role for CRP in PE.     

Genetic Polymorphisms of Stromal Interaction Molecule 1 Associated with the Erythrocyte Sedimentation Rate and C-Reactive Protein in HLA-B27 Positive Ankylosing Spondylitis Patients

Ankylosing spondylitis (AS) is a chronic inflammation of the sacroiliac joints, spine and peripheral joints. The development of ankylosing spondylitis is still unclear. Genetics factors such as human leukocyte antigen HLA-B27 and ERAP1 have been widely reported to associate to AS susceptibility. In this study, we enrolled 361 AS patients and selected four tagging single nucleotides polymorphisms (tSNPs) at STIM1 gene. The correlation between STIM1 genetic polymorphisms and AS activity index (BASDAI, BASFI, BAS-G) as well as laboratory parameters of inflammation (erythrocyte sedimentation rate (ESR) and C-reactive protein (CRP)) were tested. Our results indicated that HLA-B27 positive AS patients who are carrying the minor allele homozygous G/G genotype of SNP rs3750996 significantly associated with a higher level of ESR in serum. Furthermore, rs3750996/rs3750994 pairwise allele analysis indicated that G-C haplotypes also significantly correlated with higher level of ESR as well as CRP. These findings provide a better understanding of STIM1 genetic contribution to the pathogenesis of AS.     

XRCC1 Polymorphisms and Urinary 8-Hydroxydeoxyguanine Levels Are Associated with Urothelial Carcinoma

The aim of this study was to examine the associations between the combined effects of urinary 8-Hydroxydeoxyguanine (8-OHdG) level and polymorphisms of XRCC1 Arg194Trp and XRCC1 Arg399Gln on the risk of urothelial carcinoma (UC). We conducted a hospital-based case-control study that included 168 cases of UC and 336 age- and gender-matched healthy controls. We used polymerase chain reaction and restriction fragment length polymorphism analyses to examine the genotypes of XRCC1 Arg194Trp and XRCC1 Arg399Gln. We used a competitive in vitro enzyme-linked immunosorbent assay to determine urinary 8-OHdG levels. The XRCC1 399 Gln/Gln genotype and the XRCC1 194 Arg/Arg genotype were positively correlated to UC (OR [95%CI] = 2.27 [1.20–4.27] and 1.59 [1.06–2.36], respectively). Urinary 8-OHdG levels were associated with UC in a dose-dependent manner. Participants with the XRCC1 (Arg399Gln) Gln/Gln genotype or the G-C/A-C haplotype of XRCC1 and a high urinary 8-OHdG level had a significantly higher risk of UC than those with the Arg/Arg + Arg/Gln genotype or the G-T haplotype and a low urinary 8-OHdG level. This is the first study to investigate the combined effect of urinary 8-OHdG level and XRCC1 polymorphisms on UC risk. The findings are especially meaningful for participants with XRCC1 399Gln or XRCC1 Arg194 genotypes and a high urinary 8-OHdG level, since these variables are associated with an increased risk of UC.     

Levels of C-Reactive Protein Associated with High and Very High Cardiovascular Risk Are Prevalent in Patients with Rheumatoid Arthritis

Objective

C-reactive protein (CRP) levels>3 mg/L and>10 mg/L are associated with high and very high cardiovascular risk, respectively, in the general population. Because rheumatoid arthritis (RA) confers excess cardiovascular mortality, we determined the prevalence of these CRP levels among RA patients stratified on the basis of their RA disease activity.

Methods

We evaluated physician and patient global assessments of disease activity, tender and swollen 28 joint counts, erythrocyte sedimentation rate (ESR), and CRP measured in a single clinic visit for 151 RA patients. Disease activity was calculated using the Clinical Disease Activity Index (CDAI) and the Disease Activity Score 28 Joints (DAS28-ESR and DAS28-CRP).

Results

Median CRP level was 5.3 mg/L. 68% of patients had CRP>3 mg/L, and 25% had CRP>10 mg/L. Of those with 0–1 swollen joints (n = 56), or 0–1 tender joints (n = 81), 64% and 67%, respectively, had CRP>3 mg/L, and 23% and 20%, respectively, had CRP>10 mg/L. Of those with remission or mildly active disease by CDAI (n = 58), DAS28-ESR (n = 39), or DAS28-CRP (n = 70), 49–66% had CRP>3 mg/L, and 10–14% had CRP>10 mg/L. Of patients with moderate disease activity by CDAI (n = 51), DAS28-ESR (n = 78), or DAS28-CRP (n = 66), 67–73% had CRP>3 mg/L, and 25–33% had CRP>10 mg/L.

Conclusion

Even among RA patients whose disease is judged to be controlled by joint counts or standardized disease scores, a substantial proportion have CRP levels that are associated high or very high risk for future cardiovascular events in the general population.     

Identification of Genome-Wide Variants and Discovery of Variants Associated with Brassica rapa Clubroot Resistance Gene Rcr1 through Bulked Segregant RNA Sequencing

Clubroot, caused by Plasmodiophora brassicae, is an important disease on Brassica species worldwide. A clubroot resistance gene, Rcr1, with efficacy against pathotype 3 of P. brassicae, was previously mapped to chromosome A03 of B. rapa in pak choy cultivar “Flower Nabana”. In the current study, resistance to pathotypes 2, 5 and 6 was shown to be associated with Rcr1 region on chromosome A03. Bulked segregant RNA sequencing was performed and short read sequences were assembled into 10 chromosomes of the B. rapa reference genome v1.5. For the resistant (R) bulks, a total of 351.8 million (M) sequences, 30,836.5 million bases (Mb) in length, produced 120-fold coverage of the reference genome. For the susceptible (S) bulks, 322.9 M sequences, 28,216.6 Mb in length, produced 109-fold coverage. In total, 776.2 K single nucleotide polymorphisms (SNPs) and 122.2 K insertion / deletion (InDels) in R bulks and 762.8 K SNPs and 118.7 K InDels in S bulks were identified; each chromosome had about 87% SNPs and 13% InDels, with 78% monomorphic and 22% polymorphic variants between the R and S bulks. Polymorphic variants on each chromosome were usually below 23%, but made up 34% of the variants on chromosome A03. There were 35 genes annotated in the Rcr1 target region and variants were identified in 21 genes. The numbers of poly variants differed significantly among the genes. Four out of them encode Toll-Interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat proteins; Bra019409 and Bra019410 harbored the higher numbers of polymorphic variants, which indicates that they are more likely candidates of Rcr1. Fourteen SNP markers in the target region were genotyped using the Kompetitive Allele Specific PCR method and were confirmed to associate with Rcr1. Selected SNP markers were analyzed with 26 recombinants obtained from a segregating population consisting of 1587 plants, indicating that they were completely linked to Rcr1. Nine SNP markers were used for marker-assisted introgression of Rcr1 into B. napus canola from B. rapa, with 100% accuracy in this study.     

A Genome Wide Association Study Identifies Common Variants Associated with Lipid Levels in the Chinese Population

Plasma lipid levels are important risk factors for cardiovascular disease and are influenced by genetic and environmental factors. Recent genome wide association studies (GWAS) have identified several lipid-associated loci, but these loci have been identified primarily in European populations. In order to identify genetic markers for lipid levels in a Chinese population and analyze the heterogeneity between Europeans and Asians, especially Chinese, we performed a meta-analysis of two genome wide association studies on four common lipid traits including total cholesterol (TC), triglycerides (TG), low-density lipoprotein cholesterol (LDL) and high-density lipoprotein cholesterol (HDL) in a Han Chinese population totaling 3,451 healthy subjects. Replication was performed in an additional 8,830 subjects of Han Chinese ethnicity. We replicated eight loci associated with lipid levels previously reported in a European population. The loci genome wide significantly associated with TC were near DOCK7, HMGCR and ABO; those genome wide significantly associated with TG were near APOA1/C3/A4/A5 and LPL; those genome wide significantly associated with LDL were near HMGCR, ABO and TOMM40; and those genome wide significantly associated with HDL were near LPL, LIPC and CETP. In addition, an additive genotype score of eight SNPs representing the eight loci that were found to be associated with lipid levels was associated with higher TC, TG and LDL levels (P = 5.52×10^-16, 1.38×10^-6 and 5.59×10^-9, respectively). These findings suggest the cumulative effects of multiple genetic loci on plasma lipid levels. Comparisons with previous GWAS of lipids highlight heterogeneity in allele frequency and in effect size for some loci between Chinese and European populations. The results from our GWAS provided comprehensive and convincing evidence of the genetic determinants of plasma lipid levels in a Chinese population.