Genome-Wide Association Studies Suggest Shared Polymorphisms Are Associated with Severity of Sickle Cell Anemia and Exceptional Longevity.

Blood ◽  
2008 ◽  
Vol 112 (11) ◽  
pp. 1446-1446
Author(s):  
Paola Sebastiani ◽  
Nadia Timofeev ◽  
Steven H. Hartley ◽  
Daniel Dworkis ◽  
Lindsay Farrer ◽  
...  

Abstract Genome-wide association studies (GWAS) allow an assessment of associations between single nucleotide polymorphisms (SNPs) and phenotypes or traits of interest in a non-hypothesis driven manner. Previously, based on limited candidate gene association analysis, we showed that survival in sickle cell anemia and exceptional longevity (EL) in the general population share common genetic modifiers (Blood, 52a, 2007). This preliminary result suggested that aging mechanisms and associated genes might play a role in the variability of sickle cell anemia. Using GWAS, we now report strong evidence supporting this conjecture. We conducted a GWAS using an Illumina platform that permits genotyping up to 1 million haplotype-tagging SNPs spread across the genome, as well as other types of genetic variation, in large populations. We used the Illumina 610K SNP array to discover SNPs that are associated with different degrees of severity of sickle cell anemia in 684 patients. Patients were assigned to either a severe or mild disease category based on an integrated measure of sickle cell anemia severity that was determined by a network model that assigns a score predicting the risk of death (Blood110: 272, 2007). In parallel, we used the Illumina 370K SNP and the Illumina 1M SNP arrays to discover SNPs associated with EL in 877 centenarians enrolled in the New England Centenarian Study and 1,850 younger controls. In both studies, each SNP was tested for association with the traits of severe or less severe sickle cell anemia and EL using Bayesian tests of general, dominant and recessive associations (BMC Genet.9, 2008). We then identified those SNPs satisfying these 3 criteria: at least one model of association was 10 times more likely than no association in the GWAS of EL; the same model of association was at least 3 times more likely (because of the smaller sample size) than no association in the GWAS of sickle cell anemia severity, the same allele was more frequent in centenarians and in sickle cell anemia patients with milder disease. This analysis identified 140 SNPs in more than 50 genes and some intergenic regions that showed robust and consistent associations. This number is more than twice the number that would be expected by chance. Among the most ‘significant’ genes with associated SNPs were ARFGEF2, ADAMTS12, DOK5, DPP10, FGF21, KCNQ1, IRF4, MYO3B NAIF1, TNNI3K; more than one SNP was found in ARFGEF2, NAIF1, DPP10, SORCS3, TNNI3K. KCNQ1 has a putative role in blood circulation and regulation of heart contraction. The frequency of the common genotype for SNP rs108961 increases by almost 60% in sickle cell anemia patients with severe disease (27% versus 43%). The same common genotype in random Caucasian controls has frequency 34% that decreases to 29% in centenarians. Mutations in this gene are associated with long and short QT syndrome, with familial atrial fibrillation, heart disease and sudden death. SNPs in 2 of the genes (HAO2, a peroxisome protein involved in fatty acid oxidation, and MAP2K1, a MAP kinase involved in multiple biochemical signals) that were significantly associated with both sickle cell disease severity and EL in our earlier candidate gene studies, were also associated in the GWAS. GWAS also revealed significant association with CDKN2A, a cyclin-dependent kinase that has been associated with Type 2 diabetes, risk of myocardial infarction and triglyceride levels in several GWAS, and with FGF21, the fibroblast growth factor 21 precursor that has been shown to regulate glucose metabolism. CDKN2A has been associated with disease free survival in other studies. Common metabolic pathways are likely to influence the chance of developing complications of Mendelian and multigenic diseases and the likelihood of achieving EL. This might explain the commonality of genes whose SNPs are associated with the vascular complications of sickle cell anemia, arteriosclerosis and diabetes. A new paradigm suggests that hitherto unexpected genetic differences modulate a limited number of pathways that form a common route toward determining good health and disease.

Blood ◽  
2010 ◽  
Vol 115 (9) ◽  
pp. 1815-1822 ◽  
Author(s):  
Nadia Solovieff ◽  
Jacqueline N. Milton ◽  
Stephen W. Hartley ◽  
Richard Sherva ◽  
Paola Sebastiani ◽  
...  

Abstract In a genome-wide association study of 848 blacks with sickle cell anemia, we identified single nucleotide polymorphisms (SNPs) associated with fetal hemoglobin concentration. The most significant SNPs in a discovery sample were tested in a replication set of 305 blacks with sickle cell anemia and in subjects with hemoglobin E or β thalassemia trait from Thailand and Hong Kong. A novel region on chromosome 11 containing olfactory receptor genes OR51B5 and OR51B6 was identified by 6 SNPs (lowest P = 4.7E−08) and validated in the replication set. An additional olfactory receptor gene, OR51B2, was identified by a novel SNP set enrichment analysis. Genome-wide association studies also validated a previously identified SNP (rs766432) in BCL11A, a gene known to affect fetal hemoglobin levels (P = 2.6E−21) and in Thailand and Hong Kong subjects. Elements within the olfactory receptor gene cluster might play a regulatory role in γ-globin gene expression.


Blood ◽  
2008 ◽  
Vol 112 (11) ◽  
pp. 2471-2471 ◽  
Author(s):  
Nadia Timofeev ◽  
Paola Sebastiani ◽  
Steven H. Hartley ◽  
Clinton T. Baldwin ◽  
Martin H. Steinberg

Abstract Fetal hemoglobin (HbF) is the major genetic modulator of sickle cell anemia. Candidate gene-based and genome-wide association studies (GWAS) have provided strong evidence that single nucleotide polymorphisms (SNPs) linked to genes on chromosome 6q (HBS1L-MYB) and 2p (BCL11A), along with elements in cis to HBG help determine HbF concentration in untreated patients with sickle cell anemia and β thalassemia, and in normal individuals. The HbF response to hydroxyurea (HU) varies considerably among treated patients, even when compliance with treatment is good and patients are treated under controlled conditions. This suggests that genetic factors might affect the response to treatment with this agent. In the Multicenter Study of Hydroxyurea, 299 patients were randomized to receive either HU titrated to maximum tolerated doses, or a placebo, and HbF levels were measured before and at the completion of the randomized phase of the study. In 123 HU-treated patients, we completed GWAS using Illumina 370K chips that include approximately 350,000 haplotype tagging SNPs, and studied the association of SNPs with the change in HbF from baseline levels to levels measured at the end of the active treatment portion of the study. We conducted a GWAS using the analytical program PLINK, of approximately 273K SNPs with minor allele frequency >0.05, using linear regression and an additive model of inheritance. We selected for further investigation those SNPs with association that reached 0.05 significance, after we adjusted for sex. Because of the limited sample size that results in relatively large p-values, no single SNP reached so-called genome-wide significance after correcting for multiple comparisons using a Bonferroni correction (p-value <10-7) or 5% false discovery rate. Two SNPs had an association with p value <10–6 and 27 SNPs reached at least 10–5 significance. Noticeably, the SNP rs6899351 in FABP7 in 6q22.31 was associated with the largest increment in HbF after treatment with HU (6.9% change per copy of allele G, p-value 4 ×10–5). We also identified 2 SNPs in PDE7B (6q23.3) that were significantly associated with positive changes of HbF and 3 SNPs in MAP7 (6q23.3) that were significantly associated with a reduction of HbF after treatment. Using candidate gene association studies, we had previously shown that PDE7B and MAP7 were significantly associated with differential expression of HbF in sickle cell anemia. These new GWAS results suggest a regulatory role for these genes, or this region of chromosome 6q, in the HbF response to HU in sickle cell anemia. Analysis of the distribution of significant SNPs per chromosome also showed that chromosome 20 had a larger number of significant SNPs than expected at random, especially in CST9, one of a family of protease inhibitors. CST9 is tagged by 3 SNPs in the 370K array (rs2983639, rs2983640, rs10485646), 2 of which were associated with significant positive changes in HbF after treatment with HU and one with significant negative changes of HbF. Specifically, the average increase in HbF was 1.5% for each copy of allele G for SNP rs2983639 (p = 0.025), and 1.6% for each copy of allele A for SNP rs2983640 (p = 0.047), while the level of HbF decreased by approximately 1.5% for each copy of allele A for SNP rs10485646 (p = 0.035). The SNP rs2983640 is an exon variant that produces the amino acid change F-L. Although these SNPs do not individually reach genome-wide significance, cumulatively they provide strong evidence of association, as the probability that they are all simultaneously associated by chance is 10-4. Furthermore, we identified significant variants in other genes that belong to the same family of type 2 cysteine protease inhibitors, specifically 2 SNPs in CTS3 and 1 SNP in CTS5. Although the small sample size and the large number of SNPs tested suggest caution until these results are replicated in independent patient treatment groups, these preliminary findings suggest that type 2 cystatin genes and pseudogenes are associated with the HbF response to HU. If confirmed, it might be possible to use results like these to build a prognostic model of the HbF response to HU in sickle cell anemia.


2021 ◽  
Vol 12 ◽  
Author(s):  
Marlon Caicedo ◽  
Eduardo D. Munaiz ◽  
Rosa A. Malvar ◽  
José C. Jiménez ◽  
Bernardo Ordas

Senescence is an important trait in maize (Zea mais L.), a key crop that provides nutrition values and a renewable source of bioenergy worldwide. Genome-wide association studies (GWAS) can be used to identify causative genetic variants that influence the major physiological measures of senescence, which is used by plants as a defense mechanism against abiotic and biotic stresses affecting its performance. We measured four physiological and two agronomic traits that affect senescence. Six hundred seventy-two recombinant inbred lines (RILs) were evaluated in two consecutive years. Thirty-six candidate genes were identified by genome-wide association study (GWAS), and 11 of them were supported by additional evidence for involvement in senescence-related processes including proteolysis, sugar transport, and sink activity. We identified a candidate gene, Zm00001d043586, significantly associated with chlorophyll, and independently studied its transcription expression in an independent panel. Our results showed that Zm00001d043586 affects chlorophyl rate degradation, a key determinant of senescence, at late plant development stages. These results contribute to better understand the genetic relationship of the important trait senescence with physiology related parameters in maize and provide new putative molecular markers that can be used in marker assisted selection for line development.


Blood ◽  
2007 ◽  
Vol 110 (11) ◽  
pp. 150-150
Author(s):  
Paola Sebastiani ◽  
Ling Wang ◽  
Thomas Perls ◽  
Dellara F. Terry ◽  
Monty Montano ◽  
...  

Abstract Phenotypic heterogeneity is a well known characteristic of sickle cell anemia. Patients have different rates of hemolysis-related complications, like pulmonary hypertension, priapism and leg ulceration, and viscosity/vasoocclusion-related complications, like painful episodes, acute chest syndrome and osteonecrosis; they also have variation in levels of HbF and hematocrit. To integrate individual disease variables into a global measure of severity, we developed a Bayesian network model that described the complex associations of 25 clinical and laboratory variables, deriving a score that we used to define disease severity (0, least severe to 1, most severe) as the risk of death within 5 years (Sebastiani et al, Blood 2007). This initial network, validated in 2 unrelated patient populations, did not incorporate the genetic heterogeneity that is likely to modulate its components. Accordingly, we studied the association of single nucleotide polymorphisms (964 SNPs) in candidate genes (315 genes) using a Bayesian beta regression model of the severity score in 741 HBB glu6val homozygotes, aged more than 18 years. Forty-three SNPs in about 25 genes were associated with disease severity. Some associated SNPs tag genes that affect nitric oxide and oxidative biology and the endothelium, such as NOS1, ASS, KL, HMOX1, ECE1, KDR, FLT1. Homozygosity for an intronic SNP in ECE1 is associated with a increase of severity (OR=3.5). As expected, some associations were consistent with our previous findings. For example, the same SNP in ECE1 and TGFBR3, that was highly predictive of severity, was also strongly associated with sickle cell stroke (Sebastiani et al, Nature Genet 2005). Also, the association with severity of genes in the TGF-beta signaling pathway, including BMP6 and TGFBR3, were also associated with individual disease complications. Other associated genes play a less obvious role in the pathobiology of disease, e.g., HAO2, but are very strongly associated with the phenotype of severity (probability of a chance association, for HAO2, 10−6). Several of the genes associated with severity, including KL, PRKCA, FLT1 and MET have been related to aging, as suggested by gene expression profiling and studies in model organisms for aging. In genome-wide studies of the genetic basis of exceptional longevity, we found associations with some of the same genes that were associated with severity in sickle cell anemia. Perhaps increased oxidative stress, and the relentless progression of vasculopathy in sickle cell anemia, cause accelerated tissue damage that is modulated by a set of genes similar to those involved in the normal aging process. We suggest that the disease severity score can be used as a phenotype integrating many features of the disease, for genetic association studies. As we add the results of unbiased genome-wide association studies to capture polymorphisms not included in candidate gene studies, we can develop a predictive network with even greater reliability than one using only clinical and laboratory variables. Such networks might also identify pathways that could be targeted to alter the course of disease.


Genomics ◽  
2009 ◽  
Vol 93 (5) ◽  
pp. 415-419 ◽  
Author(s):  
Stefan Wilkening ◽  
Bowang Chen ◽  
Justo Lorenzo Bermejo ◽  
Federico Canzian

2021 ◽  
Author(s):  
Derek W Linskey ◽  
David C Linskey ◽  
Howard L McLeod ◽  
Jasmine A Luzum

The primary research approach in pharmacogenetics has been candidate gene association studies (CGAS), but pharmacogenomic genome-wide association studies (GWAS) are becoming more common. We are now at a critical juncture when the results of those two research approaches, CGAS and GWAS, can be compared in pharmacogenetics. We analyzed publicly available databases of pharmacogenetic CGAS and GWAS (i.e., the Pharmacogenomics Knowledgebase [PharmGKB®] and the NHGRI-EBI GWAS catalog) and the vast majority of variants (98%) and genes (94%) discovered in pharmacogenomic GWAS were novel (i.e., not previously studied CGAS). Therefore, pharmacogenetic researchers are not selecting the right candidate genes in the vast majority of CGAS, highlighting a need to shift pharmacogenetic research efforts from CGAS to GWAS.


Blood ◽  
2006 ◽  
Vol 108 (11) ◽  
pp. 1221-1221
Author(s):  
Paola Sebastiani ◽  
Maria M. Abad-Grau ◽  
Alberto A. Riva ◽  
Vikki G. Nolan ◽  
Efthimia Melista ◽  
...  

Abstract The complex regulatory environment controlling HbF levels provides ample opportunity for modulating expression of this key modifier of sickle cell anemia (HbSS). Identifing genetic modulators by single gene or biochemical pathway-based candidate gene association studies is limited by the a priori hypothesis that they are responsible for modulation of HbF expression and thus, fail to capture the totality of possible modulators. In contrast, genome-wide association (GWA) studies using hundreds of thousands of SNPs across the entire genome can capture an unbiased assessment of gene association or linkage disequilibrium (LD) with a phenotype. A disadvantage of GWA is the cost when analyzing thousands of subjects. We developed a GWA screening strategy using pooled DNA samples followed by analysis using the Sentrix® HumanHap300 genotyping beadchip that includes over 317,000 tag SNPs, many derived from the Phase I of the HapMap project. These SNPs are within 10kb of a gene, or in regions that are evolutionary conserved, and incorporate the majority of the variation in regions of the genome with high LD by capturing approximately 80% of the loci genotyped in HapMap Phase I and II. DNA concentration for pooling was measured by RNAase-P assays and each pool was run in duplicate. Pools of approximately 60 DNA samples each were analyzed; pools included age, sex and β-globin gene haplotype-matched HbSS patients with the 1st and 4th quartiles of HbF levels. Raw data were processed to compute initial estimates of the allele frequencies. Estimates were normalized and SNPs with a minor allele frequency >0.15 were analyzed for allelic association using a Bayesian approach. Our analysis integrated the measure of association based on the allele frequencies, estimated from the pooled DNA samples, with the strength of LD between consecutive SNPs, information on minor allele frequencies derived from the HapMap project and an estimate of the variability between pools based on preliminary experimental data. By leveraging on this exogenous information we could reduce false positive associations without imposing unnecessarily stringent corrections for multiple comparisons that may lead to substantial loss of power. We identified more than 100 novel genes potentially involved with regulation of HbF and with roles in signal transduction, cell adhesion, and regulation of transcription. These genes are distributed unevenly among the whole genome, and chromosome 8 appears to be the most "enriched" for genes associated with regulation of HbF. Consistently with previous findings, our analysis confirms an involvement of genes in 6q22, 8q12 and 11p15.5 in the regulation of HbF. For further validation, we compared allele frequencies based on pooled DNA samples with those from individual genotyping of several SNPs. The high degree of agreement suggests that GWA based on pooled DNA samples can provide a cost-effective approach to the identification of the full array of genes responsible for HbF expression and can identify new candidates that must be further examined in follow-up genetic and biological studies.


Sign in / Sign up

Export Citation Format

Share Document