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.     

The Contribution of GWAS Loci in Familial Dyslipidemias

Familial combined hyperlipidemia (FCH) is a complex and common familial dyslipidemia characterized by elevated total cholesterol and/or triglyceride levels with over five-fold risk of coronary heart disease. The genetic architecture and contribution of rare Mendelian and common variants to FCH susceptibility is unknown. In 53 Finnish FCH families, we genotyped and imputed nine million variants in 715 family members with DNA available. We studied the enrichment of variants previously implicated with monogenic dyslipidemias and/or lipid levels in the general population by comparing allele frequencies between the FCH families and population samples. We also constructed weighted polygenic scores using 212 lipid-associated SNPs and estimated the relative contributions of Mendelian variants and polygenic scores to the risk of FCH in the families. We identified, across the whole allele frequency spectrum, an enrichment of variants known to elevate, and a deficiency of variants known to lower LDL-C and/or TG levels among both probands and affected FCH individuals. The score based on TG associated SNPs was particularly high among affected individuals compared to non-affected family members. Out of 234 affected FCH individuals across the families, seven (3%) carried Mendelian variants and 83 (35%) showed high accumulation of either known LDL-C or TG elevating variants by having either polygenic score over the 90^th percentile in the population. The positive predictive value of high score was much higher for affected FCH individuals than for similar sporadic cases in the population. FCH is highly polygenic, supporting the hypothesis that variants across the whole allele frequency spectrum contribute to this complex familial trait. Polygenic SNP panels improve identification of individuals affected with FCH, but their clinical utility remains to be defined.     

通过二代测序技术开发25个小熊猫的微卫星位点

采用Roche 454基因组测序平台在小熊猫Ailurus fulgens的基因组中进行测序拼装,得到了大量含有微卫星片段的序列。本文选择了其中467条序列设计引物,筛选出了115个微卫星位点,最终选出25个微卫星位点,重复单元包含TAT、GATA、AAACA、AAAACA等,并应用到42只小熊猫个体上进行了分型。观察到的等位基因数为3～10;观察杂合度和期望杂合度分别为0. 143～0. 929和0. 136～0. 869;多态信息含量为0. 128～0. 843,经检测所有位点均未偏离哈迪-温伯格平衡。这些微卫星位点可应用于小熊猫圈养种群的遗传管理,在未来还可进一步作为小熊猫的种群遗传学等研究的有力工具。     

Replication and Meta-Analysis of GWAS Identified Susceptibility Loci in Kawasaki Disease Confirm the Importance of B Lymphoid Tyrosine Kinase (BLK) in Disease Susceptibility

The BLK and CD40 loci have been associated with Kawasaki disease (KD) in two genome-wide association studies (GWAS) conducted in a Taiwanese population of Han Chinese ancestry (Taiwanese) and in Japanese cohorts. Here we build on these findings with replication studies of the BLK and CD40 loci in populations of Korean and European descent. The BLK region was significantly associated with KD susceptibility in both populations. Within the BLK gene the rs2736340-located linkage disequilibrium (LD ) comprising the promoter and first intron was strongly associated with KD, with the combined results of Asian studies including Taiwanese, Japanese, and Korean populations (2,539 KD patients and 7,021 controls) providing very compelling evidence of association (rs2736340, OR = 1.498, 1.354–1.657; P = 4.74×10⁻³¹). We determined the percentage of B cells present in the peripheral blood mononuclear cell (PBMC) population and the expression of BLK in the peripheral blood leukocytes (leukocytes) of KD patients during the acute and convalescent stages. The percentage of B cells in the PBMC population and the expression of BLK in leukocytes were induced in patients in the acute stage of KD. In B cell lines derived from KD patients, and in purified B cells from KD patients obtained during the acute stage, those with the risk allele of rs2736340 expressed significantly lower levels of BLK. These results suggest that peripheral B cells play a pathogenic role during the acute stage of KD. Decreased BLK expression in peripheral blood B cells may alter B cell function and predispose individuals to KD. These associative data suggest a role for B cells during acute KD. Understanding the functional implications may facilitate the development of B cell-mediated therapy for KD.     

Evaluating the Performance of Fine-Mapping Strategies at Common Variant GWAS Loci

The growing availability of high-quality genomic annotation has increased the potential for mechanistic insights when the specific variants driving common genome-wide association signals are accurately localized. A range of fine-mapping strategies have been advocated, and specific successes reported, but the overall performance of such approaches, in the face of the extensive linkage disequilibrium that characterizes the human genome, is not well understood. Using simulations based on sequence data from the 1000 Genomes Project, we quantify the extent to which fine-mapping, here conducted using an approximate Bayesian approach, can be expected to lead to useful improvements in causal variant localization. We show that resolution is highly variable between loci, and that performance is severely degraded as the statistical power to detect association is reduced. We confirm that, where causal variants are shared between ancestry groups, further improvements in performance can be obtained in a trans-ethnic fine-mapping design. Finally, using empirical data from a recently published genome-wide association study for ankylosing spondylitis, we provide empirical confirmation of the behaviour of the approximate Bayesian approach and demonstrate that seven of twenty-six loci can be fine-mapped to fewer than ten variants.     

Variants Identified by Hepatocellular Carcinoma and Chronic Hepatitis B Virus Infection Susceptibility GWAS Associated with Survival in HBV-Related Hepatocellular Carcinoma

Recent genome-wide association studies (GWAS) have identified several common susceptibility loci associated with the risk of hepatocellular carcinoma (HCC) or chronic hepatitis B infection (CHB). However, the relationship between these genetic variants and survival of patients with hepatitis B virus (HBV)-related HCC is still unknown. In this study, 22 single nucleotide polymorphisms (SNPs) were genotyped among 330 HBV-related HCC patients using the MassARRAY system from Sequenom. Cox proportional hazards regression was used to examine the effects of genotype on survival time under an additive model with age, sex, smoking status and clinical stage as covariates. We identified four SNPs on 6p21 (rs1419881 T>C, rs7453920 G>A,rs3997872 G>A and rs7768538 T>C), and two SNPs on 8p12 (rs2275959 C>T and rs7821974 C>T) significantly associated with survival time of HBV-related HCC patients. Our results suggest that HCC or CHB susceptibility loci might also affect the prognosis of patients with HBV-related HCC.     

Evaluation of European Schizophrenia GWAS Loci in Asian Populations via Comprehensive Meta-Analyses

Schizophrenia is a severe and highly heritable neuropsychiatric disorder. Recent genetic analyses including genome-wide association studies (GWAS) have implicated multiple genome-wide significant variants for schizophrenia among European populations. However, many of these risk variants were not largely validated in other populations of different ancestry such as Asians. To validate whether these European GWAS significant loci are associated with schizophrenia in Asian populations, we conducted a systematic literature search and meta-analyses on 19 single nucleotide polymorphisms (SNPs) in Asian populations by combining all available case-control and family-based samples, including up to 30,000 individuals. We employed classical fixed (or random) effects inverse variance weighted methods to calculate summary odds ratios (ORs) and 95 % confidence intervals (CIs). Among the 19 GWAS loci, we replicated the risk associations of nine markers (e.g., SNPs at VRK2, ITIH3/4, NDST3, NOTCH4) surpassing significance level (two-tailed P?<?0.05), and three additional SNPs in MIR137 and ZNF804A also showed trend associations (one-tailed P?<?0.05). These risk associations are in the same directions of allelic effects between Asian replication samples and initial European GWAS findings, and the successful replications of these GWAS loci in a different ethnic group provide stronger evidence for their clinical associations with schizophrenia. Further studies, focusing on the molecular mechanisms of these GWAS significant loci, will become increasingly important for understanding of the pathogenesis to schizophrenia.     

Optimizing Selection on Multiple Identified Quantitative Trait Loci in Population with Overlapping Generations

A method was developed to model and optimize selection on multiple identified quantitative trait loci (QTLs) and polygenic estimated breeding value, in order to maximize a weighted sum of cumulative response to selection over multiple years in a population with overlapping generations. The model allows for a population with multiple sex-age classes, different number of age class between sires and dams, and varied genetic contribution of the age class. The optimization problem was formulated as a multiple-stage optimal control problem and solved by a forward and backward iteration loop. The practical utility of this method was illustrated in an example of pig breeding population with overlapping generations. The selection response of this method was compared with standard QTL selection and conventional best linear unbiased prediction (BLUP) selection. Simulation results show that optimal selection achieved greater selection response than either standard QTL or conventional BLUP selections. The influence of population structure on optimal selection was significant. Optimal QTL selection and standard QTL selection were more favorable in a population with overlapping generations than discrete generations, and obtained more benefits relative to conventional BLUP selection in a population with overlapping generations. Optimal QTL selection relative to conventional BLUP selection is also more favorable following increase of genetic contribution of two-year-old boars and sows in a population with overlapping generations.     

Re-sequencing Expands Our Understanding of the Phenotypic Impact of Variants at GWAS Loci

Genome-wide association studies (GWAS) have identified >500 common variants associated with quantitative metabolic traits, but in aggregate such variants explain at most 20–30% of the heritable component of population variation in these traits. To further investigate the impact of genotypic variation on metabolic traits, we conducted re-sequencing studies in >6,000 members of a Finnish population cohort (The Northern Finland Birth Cohort of 1966 [NFBC]) and a type 2 diabetes case-control sample (The Finland-United States Investigation of NIDDM Genetics [FUSION] study). By sequencing the coding sequence and 5′ and 3′ untranslated regions of 78 genes at 17 GWAS loci associated with one or more of six metabolic traits (serum levels of fasting HDL-C, LDL-C, total cholesterol, triglycerides, plasma glucose, and insulin), and conducting both single-variant and gene-level association tests, we obtained a more complete understanding of phenotype-genotype associations at eight of these loci. At all eight of these loci, the identification of new associations provides significant evidence for multiple genetic signals to one or more phenotypes, and at two loci, in the genes ABCA1 and CETP, we found significant gene-level evidence of association to non-synonymous variants with MAF<1%. Additionally, two potentially deleterious variants that demonstrated significant associations (rs138726309, a missense variant in G6PC2, and rs28933094, a missense variant in LIPC) were considerably more common in these Finnish samples than in European reference populations, supporting our prior hypothesis that deleterious variants could attain high frequencies in this isolated population, likely due to the effects of population bottlenecks. Our results highlight the value of large, well-phenotyped samples for rare-variant association analysis, and the challenge of evaluating the phenotypic impact of such variants.     

世代重叠群体多个QTL最优化选择

一种扩展的方法能够在一个世代重叠的群体内对多个数量性状位点选择进行最优化,目的是为了在整个计划期内获得最大的累积反应加权和。该模型允许群体有多个性别年龄组、公母畜间有不同的年龄组数、各年龄组有不同的遗传贡献。整个最优化问题被描述成一个多阶段系统优化控制问题,通过一个向前和向后的迭代循环解决。用一个世代重叠的实际育种猪群的参数来评价该方法的选择效果,并和标准QTL选择和常规BLUP选择进行比较。模拟结果表明,优化选择要优于标准QTL选择和常规BLUP选择。群体结构对优化选择的影响比较明显。优化QTL选择和标准QTL选择在世代重叠的群体内比在世代离散的群体内的选择优势更明显,相对于常规BLUP选择,能够获得更大的选择优势。在世代重叠群体内随着2岁公畜遗传贡献的增大,优化选择相对于常规BLUP选择的优势越明显。     

ezQTL: A Web Platform for Interactive Visualization and Colocalization of QTLs and GWAS Loci

Genome-wide association studies(GWAS) have identified thousands of genomic loci associated with complex diseases and traits, including cancer. The vast majority of common traitassociated variants identified via GWAS fall in non-coding regions of the genome, posing a challenge in elucidating the causal variants, genes, and mechanisms involved. Expression quantitative trait locus(e QTL) and other molecular QTL studies have been valuable resources in identifying candidate causal genes from GWAS loc...     

Replication Study in a Japanese Population of Six Susceptibility Loci for Type 2 Diabetes Originally Identified by a Transethnic Meta-Analysis of Genome-Wide Association Studies

AimWe performed a replication study in a Japanese population to evaluate the association between type 2 diabetes and six susceptibility loci (TMEM154, SSR1, FAF1, POU5F1, ARL15, and MPHOSPH9) originally identified by a transethnic meta-analysis of genome-wide association studies (GWAS) in 2014.MethodsWe genotyped 7,620 Japanese participants (5,817 type 2 diabetes patients and 1,803 controls) for each of the single nucleotide polymorphisms (SNPs) using a multiplex polymerase chain reaction invader assay. The association of each SNP locus with the disease was evaluated using logistic regression analysis.ResultsOf the six SNPs examined in this study, four (rs6813195 near TMEM154, rs17106184 in FAF1, rs3130501 in POU5F1 and rs4275659 near MPHOSPH9) had the same direction of effect as in the original reports, but two (rs9505118 in SSR1 and rs702634 in ARL15) had the opposite direction of effect. Among these loci, rs3130501 and rs4275659 were nominally associated with type 2 diabetes (rs3130501; p = 0.017, odds ratio [OR] = 1.113, 95% confidence interval [CI] 1.019–1.215, rs4275659; p = 0.012, OR = 1.127, 95% CI 1.026–1.238, adjusted for sex, age and body mass index), but we did not observe a significant association with type 2 diabetes for any of the six evaluated SNP loci in our Japanese population.ConclusionsOur results indicate that effects of the six SNP loci identified in the transethnic GWAS meta-analysis are not major among the Japanese, although SNPs in POU5F1 and MPHOSPH9 loci may have some effect on susceptibility to type 2 diabetes in this population.     

Pleiotropy among Common Genetic Loci Identified for Cardiometabolic Disorders and C-Reactive Protein

Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10^-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes.     

Variants Identified in a GWAS Meta-Analysis for Blood Lipids Are Associated with the Lipid Response to Fenofibrate

A recent large-scale meta-analysis of genome-wide studies has identified 95 loci, 59 of them novel, as statistically significant predictors of blood lipid traits; we tested whether the same loci explain the observed heterogeneity in response to lipid-lowering therapy with fenofibrate. Using data from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN, n = 861) we fit linear mixed models with the genetic markers as predictors and high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol, total cholesterol, and triglyceride concentrations as outcomes. For all four traits, we analyzed both baseline levels and changes in response to treatment with fenofibrate. For the markers that were significantly associated with fenofibrate response, we fit additional models evaluating potential epistatic interactions. All models were adjusted for age, sex, and study center as fixed effects, and pedigree as a random effect. Statistically significant associations were observed between the rs964184 polymorphism near APOA1 (P-value≤0.0001) and fenofibrate response for HDL and triglycerides. The association was replicated in the Pharmacogenetics of Hypertriglyceridemia in Hispanics study (HyperTG, n = 267). Suggestive associations with fenofibrate response were observed for markers in or near PDE3A, MOSC1, FLJ36070, CETP, the APOE-APOC1-APOC4-APOC2, and CILP2. Finally, we present strong evidence for epistasis (P-value for interaction =  0.0006 in GOLDN, 0.05 in HyperTG) between rs10401969 near CILP2 and rs4420638 in the APOE-APOC1-APOC4-APOC2 cluster with total cholesterol response to fenofibrate. In conclusion, we present evidence linking several novel and biologically relevant genetic polymorphisms to lipid lowering drug response, as well as suggesting novel gene-gene interactions in fenofibrate pharmacogenetics.     

Replication Study in a Japanese Population to Evaluate the Association between 10 SNP Loci,Identified in European Genome-Wide Association Studies,and Type 2 Diabetes

AimWe performed a replication study in a Japanese population to evaluate the association between type 2 diabetes and 7 susceptibility loci originally identified by European genome-wide association study (GWAS) in 2012: ZMIZ1, KLHDC5, TLE1, ANKRD55, CILP2, MC4R, and BCAR1. We also examined the association of 3 additional loci: CCND2 and GIPR, identified in sex-differentiated analyses, and LAMA1, which was shown to be associated with non-obese European type 2 diabetes.MethodsWe genotyped 6,972 Japanese participants (4,280 type 2 diabetes patients and 2,692 controls) for each of the 10 single nucleotide polymorphisms (SNPs): rs12571751 in ZMIZ1, rs10842994 near KLHDC5, rs2796441 near TLE1, rs459193 near ANKRD55, rs10401969 in CILP2, rs12970134 near MC4R, rs7202877 near BCAR1, rs11063069 near CCND2, rs8108269 near GIPR, and rs8090011 in LAMA1 using a multiplex polymerase chain reaction invader assay. The association of each SNP locus with the disease was evaluated using a logistic regression analysis.ResultsAll SNPs examined in this study had the same direction of effect (odds ratio > 1.0, p = 9.77 × 10^-4, binomial test), as in the original reports. Among them, rs12571751 in ZMIZ1 was significantly associated with type 2 diabetes [p = 0.0041, odds ratio = 1.123, 95% confidence interval 1.037–1.215, adjusted for sex, age and body mass index (BMI)], but we did not observe significant association of the remaining 9 SNP loci with type 2 diabetes in the present Japanese population (p ≥ 0.005). A genetic risk score, constructed from the sum of risk alleles for the 7 SNP loci identified by un-stratified analyses in the European GWAS meta-analysis were associated with type 2 diabetes in the present Japanese population (p = 2.3 × 10^-4, adjusted for sex, age and BMI).ConclusionsZMIZ1 locus has a significant effect on conferring susceptibility to type 2 diabetes also in the Japanese population.     

Quantitative Trait Loci for Salinity Tolerance Identified under Drained and Waterlogged Conditions and Their Association with Flowering Time in Barley (Hordeum vulgare. L)

Introduction

Salinity is one of the major abiotic stresses affecting crop production via adverse effects of osmotic stress, specific ion toxicity, and stress-related nutritional disorders. Detrimental effects of salinity are also often exacerbated by low oxygen availability when plants are grown under waterlogged conditions. Developing salinity-tolerant varieties is critical to overcome these problems, and molecular marker assisted selection can make breeding programs more effective.

Methods

In this study, a double haploid (DH) population consisting of 175 lines, derived from a cross between a Chinese barley variety Yangsimai 1 (YSM1) and an Australian malting barley variety Gairdner, was used to construct a high density molecular map which contained more than 8,000 Diversity Arrays Technology (DArT) markers and single nucleotide polymorphism (SNP) markers. Salinity tolerance of parental and DH lines was evaluated under drained (SalinityD) and waterlogged (SalinityW) conditions at two different sowing times.

Results

Three quantitative trait loci (QTL) located on chromosome 1H, single QTL located on chromosomes 1H, 2H, 4H, 5H and 7H, were identified to be responsible for salinity tolerance under different environments. Waterlogging stress, daylight length and temperature showed significant effects on barley salinity tolerance. The QTL for salinity tolerance mapped on chromosomes 4H and 7H, QSlwd.YG.4H, QSlwd.YG.7H and QSlww.YG.7H were only identified in winter trials, while the QTL on chromosome 2H QSlsd.YG.2H and QSlsw.YG.2H were only detected in summer trials. Genes associated with flowering time were found to pose significant effects on the salinity QTL mapped on chromosomes 2H and 5H in summer trials. Given the fact that the QTL for salinity tolerance QSlsd.YG.1H and QSlww.YG.1H-1 reported here have never been considered in the literature, this warrants further investigation and evaluation for suitability to be used in breeding programs.