key: cord-325155-lqzgz6p3
authors: Gallo, Juan E.; Ochoa, Juan E.; Warren, Helen R.; Misas, Elizabeth; Correa, Monica M.; Gallo-Villegas, Jaime A.; Bedoya, Gabriel; Aristizábal, Dagnóvar; McEwen, Juan G.; Caulfield, Mark J.; Parati, Gianfranco; Clay, Oliver K.
title: Hypertension and the roles of the 9p21.3 risk locus: classic findings and new association data
date: 2020-09-15
journal: nan
DOI: 10.1016/j.ijchy.2020.100050
sha: 
doc_id: 325155
cord_uid: lqzgz6p3

Background The band 9p21.3 contains an established genomic risk zone for cardiovascular disease (CVD). Since the initial 2007 Wellcome Trust Case Control Consortium study (WTCCC), the increased CVD risk associated with 9p21.3 has been confirmed by multiple studies in different continents. However, many years later there was still no confirmed report of a corresponding association of 9p21.3 with hypertension, a major CV risk factor, nor with blood pressure (BP). Theory In this contribution, we review the bipartite haplotype structure of the 9p21.3 risk locus: one block is devoid of protein-coding genes but contains the lead CVD risk SNPs, while the other block contains the first exon and regulatory DNA of the gene for the cell cycle inhibitor p15. We consider how findings from molecular biology offer possibilities of an involvement of p15 in hypertension etiology, with expression of the p15 gene modulated by genetic variation from within the 9p21.3 risk locus. Results We present original results from a Colombian study revealing moderate but persistent association signals for BP and hypertension within the classic 9p21.3 CVD risk locus. These SNPs are mostly confined to a ‘hypertension island’ that spans less than 60 kb and coincides with the p15 haplotype block. We find confirmation in data originating from much larger, recent European BP studies, albeit with opposite effect directions. Conclusion Although more work will be needed to elucidate possible mechanisms, previous findings and new data prompt reconsidering the question of how variation in 9p21.3 might influence hypertension components of cardiovascular risk.

Legend to Graphical Abstract: Schematic depiction of the main observation presented and discussed in this article. Two adjacent haplotype blocks characterize the 9p21.3 cardiovascular risk locus: left, the block or island containing the first part of the p15 gene and its wellcharacterized promoter, in which we observed clearly elevated associations (red) with blood pressure (DBP, SBP) and/or hypertension in a Colombian and a European study sample, and right, the block hypertension and BP association 'Hypertension island' * (haplotype block < 60 kb) Lead CVD risk SNPs (haplotype block < 60 kb)

Significance of association strongest CVD association current threshold for genome-wide reporting ( p = 5 x 10 -8 )

fair CVD association contains start and promoter of p15 gene contains no protein-coding genes

It had long been a widely held belief that common genetic variation in the established cardiovascular risk locus of cytogenetic band 9p21.3, as discovered and delimited via corresponding genetic associations in 2007 [1, 2, 3] , is not reproducibly associated with high blood pressure or hypertension, a prime risk factor for cardiovascular disease [4, 5] .

The pleiotropic nature of the ≈ 120-kb classic risk locus of 9p21.3, the strongest known single contributor in the human genome to genetic risk of cardiovascular disease [6, 7, 8] , viewed together with the complexity of the genetic and molecular regulation of hypertension, prompted us to reopen the question if this region might, after all, be associated with a contribution of hypertension to the increased cardiovascular risk that characterizes the locus.

In the Theory section, we briefly depict the classic two-block structure of the region in the light of current knowledge, and review some findings and hypotheses that could admit a step of a hypertension etiology being modulated at this locus.

In the Results we first present results from an original, small Colombian association study focusing on 9p21.3 variation. We then present largely corroborating results from much larger, recently published European studies, which we identified in a second step. Taken together, the evidence accumulating so far suggests that common genetic variation in a well-described block (< 60 kb) of the 9p21.3 risk region -and, more specifically, in the regulatory DNA of the p15 gene it harbors -may play a role in promoting hypertension, for example via vascular modifications in resistance arteries (arterioles). 

Demographic and clinical information, including selected cardiovascular risk factors, was collected for all participants in a study conducted in Medellín, Colombia (see Supplementary Data). Systolic (S) and diastolic (D) blood pressure (BP) levels were defined by the average of two conventional auscultatory BP measurements through a mercury manometer according to an approach in line with current European guidelines [9] . Presence of hypertension was identified based on physician's diagnosis, prescription of antihypertensive drugs, average office SBP ≥ 140 and/or DBP ≥ 90 mmHg, or any combination of these possibilities.

A total of 357 participants were genotyped. In view of this small sample size and in line with traditional case-control methodology [10] , [11, Chapter 9] , the highest and lowest DBP tertiles of the population were slightly overrepresented in the genotyped individuals.

DNA was extracted from white blood cells following standard salting-out procedures and genotyped by LGC Group, UK with KASP TM technology (https://www. biosearchtech.com). All SNPs genotyped were biallelic in our study sample. SNP rs10738605 is triallelic in some other 1KG populations (not CLM), its rare third allele frequency being 0.0014 in AMR, 0.006 in EUR and 0.0008 in AFR.

Statistical analyses were performed considering BP either as a continuous or a categorical variable (above or below the current [9] threshold for hypertension diagnosis). Effect sizes were defined as differences of mean values for continuous variables (e.g., SBP, DBP), and as odds ratios for dichotomous variables (e.g., presence or absence of hypertension). Variables considered as possible covariates in the statistical analyses included gender, age, smoking and body mass index (BMI). To calculate association p-values and effect sizes between SNP states and phenotypic variables we used the R package snpassoc [12, 13] . Models fitted included codominant, dominant, recessive, overdominant and log-additive models.

Participants signed an informed consent formulated for this genotype-phenotype association study. The study was approved by the Ethical Committee of the Corporación para Investigaciones Biológicas, Medellín, Colombia. All procedures were performed in agreement with the Declaration of Helsinki.

Additional Material and Methods are given in Supplementary Material S1.

Methods of European studies were as described in the publications presenting those studies [14, 15, 16] . The originally described risk locus of 9p21.3 has a bipartite or two-segment structure, in which each of the two mutually adjacent segments is dominated by a distinct haplotype block. Although the haplotypes (i.e., states) of the two blocks are correlated, via a linkage disequilibrium originally estimated as around 0.6 (for average |D |) [2] , their frequencies are different. The minor allele frequency (MAF) plots shown in Figure 1 depict individual SNP-and haplotype-level variation in the locus. Two of the superpopulations that were sequenced and genotyped by the 1000 Genomes Project (1KG; [17]) are represented: AMR (populations from the Americas having a Native American ancestry component) and EUR (populations of European ancestry). The 'left-hand' block of the risk region, colored green in Figure 1 , is composed of essentially three haplotype classes, giving one low and two higher, similar minor allele frequencies (separated by a green horizontal bar). The 'right-hand' block is dominated by two haplotype classes and has fairly constant allele frequencies close to 50%; the lead CVD risk SNPs of the region are at its rightmost fringe, denoted by ochre triangles.

The left-hand block includes the upstream regulatory DNA and start of the protein-coding gene encoding p15 (transcribed leftward). The Figure shows a big difference between AMR and EUR in this block. The Figure also shows a matrix, in 0/1 haplotype (schema) coding, for 16 SNPs that are representative of variation within this block, as obtained by 1000 Genomes' sampling of 94 individuals (i.e., 188 chromosomes) from CLM (Medellín, Colombia, an AMR population).

Where three or more major haplotype classes dominate a region, individual bialellic SNPs can provide only aggregated resolution of what is happening at the haplotype level [18] . Here, SNP classes s1 (blue identifiers), s2 (red) or s3 provide information on, respectively, the contrast of haplotype class c2, c1 or c3 versus the remaining two haplotype classes of the haplotype block.

The CVD risk locus in 9p21.3 is a confirmed risk locus not only for cardiovascular disease but also for other diseases, including cancers. The pleiotropies appear to be partly antagonistic, in that the risk allele for one disease can be the protective allele for another disease. For example, the NHGRI-EBI GWAS Catalog (www.ebi.ac.uk/gwas) reports SNPs of the left block where the protective allele for an age-related condition (e.g., CVD, glaucoma) would be the risk allele for a cancer-related condition (e.g., breast cancer, glioma, pediatric brain tumor, endometriosis; see Supplementary Material S2.1). Antagonistic pleiotropies have been noted or proposed before for this locus ( [19] ; see also [20] ), with data suggesting a corresponding positive selection in some populations having Native American ancestry [21] .

Recent literature on the risk locus of 9p21.3 and its associations has emphasized the enigmatic long noncoding RNA (lncRNA) gene ANRIL, which has regulatory roles that are still only partly understood, but that can affect expression of proteincoding genes in the risk locus [22] .

Earlier molecular biology studies had focused on two protein-coding genes for cyclin-dependent kinase inhibitors, which when activated can arrest or inhibit the cell cycle at the G1 stage: the gene CDKN2A encodes both p14 ARF (p14) and p16 INK4a (p16), and CDKN2B encodes p15 INK4b (p15).

A third protein-coding gene, MTAP, which is located outside the risk locus, is also relevant for understanding the region. MTAP encodes a key enzyme in the biosynthetic pathway of polyamines (which plays an important role in the stabilization of atherosclerotic plaques); there is also evidence that MTAP may be partly regulated from within the 9p21.3 CVD risk locus despite the intervening distance [23, 24, 16] . 3.4. The p15 gene Some common genetic variation in the p15 haplotype block is in regulatory DNA, where it could in principle modulate transcription of the p15 gene. Thus, the common SNP rs2069418, in an established promoter of the p15 gene, borders a tetranucleotide that was critical for C/EBPβ binding and transcription of the p15 gene in experiments in vitro ( Figure 2 ) [25, 26, 27] ; its common variation might also modulate a predicted stemloop (see Supplementary Material S2.2). The functionally strategic position of SNP rs2069418 was also noted before, but in the apparently antagonistic context of glioma risk [28] .

The well-studied promoter of the p15 gene [29] is a 'battlefield' in which the proto-oncogene c-Myc and TGFβ/Smad struggle to control the p15 gene. Indeed, control of this gene can be important: where expression of the cell cycle inhibitor and tumor suppressor p15 is suppressed, cancer risk can rise [27] . TGFβ is a main 7 J o u r n a l P r e -p r o o f regulator of blood vessel development and maintenance, but it acts through several alternative 'arms', or pathways, e.g., via crosstalk. We ask when and where the arm of TGFβ's activities that utilizes p15, as part of the TGFβ/Smad pathway, might be employed under physiological conditions in a context that would be relevant for the genesis of hypertension. Figure 3 shows a rough sketch of elements of three pathways in which TGFβ can play a role and/or modulate hypertension etiologies: the TGFβ/Smad pathway with target gene p15, in the central lane, and two other pathways that are trans-activatable by 'crosstalk'' with TGFβ, respectively mediated by the type 1 angiotensin-II receptor (AT1R / classic RAS pathway, left) and the epidermal growth factor receptor (EGFR, right).

The central and rightmost of these three pathways and their crosstalk, for which roles in hypertension etiologies remain to be elucidated in more detail, are examined in this context in Supplementary Material S2.3; refs. [27, 30, 31, 32] describe experimental results that may be relevant.

In addition to such cues from experiments, an independent line of evidence supporting a likely role for the central pathway comes from large-scale wholegenome association studies and metastudies that can now identify whole networks showing evidence of collective association with blood pressure [15, 16] . The results indicate that the TGFβ/Smad pathway plays an important role in BP regulation, as it is enriched for BP-associated genes [16] . It then seems a plausible hypothesis that p15, as an important TGFβ/Smad-responsive target, could be involved in the modulating of blood pressure.

Although the complexity and possible sensitivity of such a network can render difficult the prediction of how a genetic change at a given locus will affect BP or the risk of developing hypertension, already Figure 3 hints that the p15 context puts us in a locus and a scenario where key players, such as cell cycle inhibitors, are only a step away from etiologically familiar routes to pathogenesis of hypertension or its prevention. 

We conducted a local, small study in the city of Medellín, Colombia (n=357 individuals without overt CVD) to explore possible associations between hypertension and 9p21.3 (67 common SNPs spanning > 300 kb). One of the motivations for the study was that we knew of no 9p21.3 studies of a population living in Latin America. It was immediately clear from the sample size that even for minor allele frequencies (MAF) above 20% we would find at best nominally significant evidence for associations from this study alone (0.05 > p ≥ 5 × 10 −8 ; [33] ). Still, we wished to see what could be achieved by collecting and analyzing local data of modest size and integrating them into the global knowledge base [34, p. 6] , [35] .

The screening of 9p21.3 for associations in the Colombian study sample revealed an 'island' spanning 16 consecutive SNPs (58.2 kb), in which SBP and DBP levels as well as the presence of hypertension gave consistent association signals, and the island coincided almost exactly with the 'left' haplotype block of the 9p21.3 risk locus (Figures 1 and 2 ; see Theory). Effect sizes for the 16 SNPs in the island were approximately allele-additive. SNPs in the immediately flanking regions gave only sporadic or no association signals (Supplementary Material S3.1). For instance, despite the small sample size, SNP rs564398 in the island achieved p values of 2 × 10 −6 for hypertension, and 2 × 10 −3 / 5 × 10 −4 for SBP/DBP increase, with a per-allele Odds Ratio of 2.44 for hypertension and a mean difference (effect size) of 2.94 / 3.14 mmHg for SBP/DBP levels, after correcting for sex and age. Comparable results were obtained when including smoking and/or BMI as covariates (Supplementary Material S1.4 and Supplementary Data).

We compared the inferred risk and protector alleles for hypertension or blood pressure at our genotyped SNPs with the risk and protector alleles of disease phenotypes that had been previously reported at those same SNPs. Risk alleles were defined as those for which mean difference or beta was positive (continuous variables), for which effect size was greater than one (dichotomous variables), or that were explicitly labeled as 'risk allele'. Previously published associations listed in the NHGRI-EBI GWAS database suggested two mutually 'antagonistic' groups (see section 3.2, Supplementary Material S2.1, or [19, Supplementary Information]), the risk alleles for the CVD-or aging-related conditions (e.g., coronary heart disease, or intracranial aneurism) being the protector alleles for cancerrelated conditions (e.g., glioma). In our Colombian study and in the 9p21.3 island where we observed the best association signals, the inferred risk alleles for hypertension or BP (e.g., A, the major allele, for rs7865618) were also the previ-9 J o u r n a l P r e -p r o o f ously reported risk alleles for CVD-or aging-related conditions at all SNPs that we genotyped.

Although our study focused on blood pressure variables and hypertension, we did consider, and subjected to the same basic phenotype-genotype analysis, several dozen (not all independent) alternative candidates for the phenotypic/outcome variable of primary interest, selected on the basis of their hemodynamic or other relevance to the development of hypertension (listed in Supplementary Data). However, the BP and hypertension variables consistently gave by far the best support for an association at the SNPs we genotyped (see Supplementary Data for the results obtained for rs7044859 and rs2069418, representing the two equivalence classes of SNPs s1 and s2, which are characterized by high linkage disequilibrium).

A review of the data from recently published large European studies [14, 15, 16] with detailed look-ups of results in this region of interest, again showed higher DBP association signals in the same island that we had delimited in the Colombian study than in its flanking DNA, with p-values that almost reached genome-wide significance. Thus, in the studies by Warren et al. [15] and by Evangelou et al. [16] (n > 750,000), the two lowest p-values for DBP, 1.23 × 10 −7 and 9.22 × 10 −7 , were found, respectively, at the 'left' fringe of the island's haplotype block (rs3217992) and at a SNP in the interior of the island that we had not genotyped in the Colombian study (rs7874604).

Furthermore, in the European blood pressure studies [15, 16] genome-wide significance of DBP associations was attained, outside of the classic 9p21.3 CVD risk locus and its flanking regions, in the next gene MTAP (see Figure 2 and Theory), with a lowest p-value of 1.3 × 10 −10 for the sentinel SNP rs4364717 (red asterisk and red horizontal bar at left in Figure 2 ; see also the LocusZoom plot in Supplementary Material S3.2) . The BP association with the MTAP gene had not been noted in previous studies, and indeed the region had not been sufficiently covered in earlier imputation panels.

In the earlier study by Ehret et al ( [14] , n > 180,000), the best p-values observed were within the island, and again for DBP, although they did not reach values below 0.002 (see Supplementary Material, subsection S3.2). We were even able to detect concordant results as far back as 2011, in a study screening close to 2.5 million SNPs for BP associations [36] . The SNPs in the island again had much lower p-values than those in the lead CVD risk zone and regions flanking the island. Thus, the SNP that had the lowest p-value for hypertension in our Colombian study sample, rs564398, had similarly the lowest p-value for DBP in the island in ref. [36] , namely 0.0161 as its genome-wide meta-analysis p-value corrected for genomic control.

A notable difference with respect to the Colombian results was the effect direction. Indeed, directions of the BP effects in the European studies [15, 16] were consistently opposite to those observed in the Colombian study at the same SNPs. In other words, the BP effects were also opposite to the established CVD effects for those SNPs in the literature. The now established CVD effect directions within the island are the same in European populations, deduced via direct genotyping of island SNPs, as they are in populations with a Latin American or Native American ancestry component (e.g., [37, 38] ), where they can be inferred from genotyped SNPs in linkage disequilibrium located close to the island. Thus, in the European study samples, but not in the Colombian study sample, the inferred protector alleles for hypertension (i.e., earlier in disease etiology) are the known, constant risk alleles for cardiovascular disease or events (i.e., later in the etiology).

Three other European reports we found, in ref. [39] and in the PhenoScanner and Roslin GeneATLAS association databases, had varying significance and independence, but the reported effect directions appear compatible with our findings (Supplementary Material S3.3).

J o u r n a l P r e -p r o o f

The results from the trans-ethnic investigations we describe here are compatible with the persisting presence, in a South American and in a European population, of associations between BP levels/hypertension and genetic variation specifically in a previously identified haplotype block or 'island' nested within the cardiovascular risk locus of band 9p21.3. In the European population, dense coverage of the association screening also allowed the recognition of an association peak for BP in a neighboring gene, MTAP, which is located at about 140 kb from the risk locus.

According to the data collected so far, the hypertension or BP effect within the CVD risk locus 9p21.3 would have the opposite direction in the South American compared to the European population, whereas the CVD effect appears to be the same in both populations, according to published reports [37, 38] . To our knowledge, however, our study is the first screen of the 9p21.3 'island' for hypertension or BP associations in a South American population, and we also know of no dedicated CVD association study with dense coverage in this part of the 9p21.3 risk zone. It is therefore too early to generalize or extrapolate from this one small study, and we must await confirmation from independent genotypephenotype screening in future, if possible also in South American populations having other ancestry proportions, such as those of Peru or Bolivia (see, for example, [40] ). The larger networks in which the p15 gene is embedded (e.g., the pathways of Figure 3 ) are likely to be complex and resilient, and not just genetic background but also environment (e.g., diet) and local epidemiology (e.g., disease prevalences) could influence routing or rebalancing of the networks and thus affect associations.

According to theoretical arguments, genotype-phenotype associations that change direction of effect between geographic, population or other environments (e.g., in the form of crossing reaction norms [41, 42] ) could actually be favored, because a variant that is inferior in all environments would be rendered unstable and would eventually be eliminated from the gene pool [43, 41] , [44, section 5] . In line with this view, trans-ethnically significant BP associations can exhibit sizable directional inconsistency rates. For example, among 562 SNPs genome-wide that were associated with BP traits (SBP, DBP, or pulse pressure) at a high significance of p < 10 −10 in a European sample (n = 757,601) in ref. [16, Supplementary Table 18] , 20 of the 32 SNPs having p < 0.01 also in a South Asian population, and 6 of the 22 SNPs having p < 0.01 in an African population, had effect directions that were opposite to those in the European population.

As a final note, a wider view of physical or regulatory gene interactions may 12 help understand potential roles of the CVD risk locus in modulating blood pressure or hypertension. One line of investigation that could be relevant is represented by Hi-C data that capture interactions between different chromosomal regions. Thus, Hi-C data reported in ref. [16, Supplementary Table 8] show, at least in mesenchymal stem cells, an interaction cluster that includes the classic risk region and other regions from almost all of cytogenetic band 9p21.3 (4.7/5.7 Mb; for details see Supplementary Material S3.4). The clustering regions include (or immediately flank) the 9p21.3 interferon gene group, relevant in inflammation, and two genes (GADD45G and FOCAD) that have been considered of interest in cardiovascular as well as cancer contexts. Of particular note, the cluster also includes the MTAP and DMRTA1 genes, which have been recently reported to be associated with DBP [15] and SBP [16] , respectively, and are located at a distance on either side of the classic risk region. Such findings suggest that a full understanding of the potential role of p15 and the classic CVD risk locus in modulating hypertension and CVD might require a wider view that takes into account the regulatory interactions among genes distributed across, and possibly beyond, the 5.7-Mb cytogenetic band 9p21.3.

Additional (extended) Materials & Methods, Theory, and Results are included in Supplementary Material sections S1, S2, and S3, respectively. De-identified genotype data and detailed association results of the Medellín study are given in a Supplementary Data file (Excel workbook) and described in Supplementary Material section S4.

J o u r n a l P r e -p r o o f pomethylation that is disrupted in some cancers [45] ) illustrate the functional importance of the region around SNP rs2069418. Bottom: functional p15 promoter binding sites (boxes), critical subsequences (beige) and attempted experimental replacements [25, 26] that compromised normal transcription of p15 (bottom subsequences).

Legend to Figure 3 : Schematic diagram sketching three well-characterized pathways (3 vertical lanes) in which TGFβ plays a role via signaling and/or transactivation/crosstalk and which may act to promote or prevent hypertension. At the top of each lane, preparation steps needed for a master product that is essential for the pathway's activation are summarized. Asterisks and daggers indicate gene products in which common SNPs in or in the gene's vicinity have been reported to be associated with blood pressure or hypertension in ref. [15, 16] (*; (*) for pathway) or have been reviewed as being associated with BP or hypertension in ref. [46] ( †). Dashed horizontal arrows indicate experimentally observed or hypothesized crosstalk/transactivation between pathways (not just receptors). Not shown, for simplicity, is another potentially relevant system, descending alongside the classic RAS pathway at left, namely the 'nonclassical' RAS, composed primarily of the AngII/Ang III-AT2R pathway and the ACE2-Ang-(1-7)-AT7R axis [47] , which generally counteracts the effects of a stimulated classic Ang II-AT1R axis as shown in the Figure rs1063192  rs3217986  rs2069418  rs495490  rs7044859  rs496892  rs615552  rs10965215  rs564398  rs662463  rs7865618  rs1360590  rs17694493  rs10738605  rs11790231  rs2184061 1KG-CLM,% 1KG-GBR,% Class  s2  s2  s1  s1  s2  s1  s2  s3  s2  s1  s3  s1  s2  188 

J o u r n a l P r e -p r o o f rs7044859

Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls

Genomewide association analysis of coronary artery disease

A common allele on chromosome 9 associated with coronary heart disease

Chromosome 9p21.3 locus for coronary artery disease: how little we know

Alu elements in ANRIL non-coding RNA at chromosome 9p21 modulate atherogenic cell functions through trans-regulation of gene networks

Recent studies of the human chromosome 9p21 locus, which is associated with atherosclerosis in human populations. Arterioscler

From genotype to phenotype in human atherosclerosis -recent findings

Genomic risk prediction of coronary artery disease in 480,000 adults: implications for primary prevention

ESC/ESH Guidelines for the management of arterial hypertension

Haplotypes of the beta-2 adrenergic receptor associate with high diastolic blood pressure in the Caerphilly prospective study

Principles of Statistical Inference

snpassoc: an R package to perform whole genome association studies

Basic statistical analysis in genetic case-control studies

The genetics of blood pressure regulation and its target organs from association studies in 342,415 individuals

Genomewide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk

Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits

Genomes Project Consortium. A global reference for human genetic variation

Toward multiple SNP motif analyses of loci associated with phenotypic traits

Antagonistic pleiotropy and mutation accumulation influence human senescence and disease

Chromosome 9p21 SNPs associated with multiple disease phenotypes correlate with ANRIL expression

Genetic loci associated with coronary artery disease harbor evidence of selection and antagonistic pleiotropy

Long noncoding RNA ANRIL: Lnc-ing genetic variation at the chromosome 9p21 locus to molecular mechanisms of atherosclerosis

Expression of Chr9p21 genes CDKN2B (p15 INK4b ), CDKN2A (p16 INK4a , p14 ARF ) and MT AP in human atherosclerotic plaque

Cytoplasmic MTAP expression loss detected by immunohistochemistry correlates with 9p21 homozygous deletion detected by FISH in pleural effusion cytology of mesothelioma

TGFβ influences Myc, Miz-1 and Smad to control the CDK inhibitor p15 INK4b

C/EBPβ at the core of the TGFβ cytostatic response and its evasion in metastatic breast cancer cells

TGF-β family signaling in tumor suppression and cancer progression

Variants in the CDKN2B and RTEL1 regions are associated with high-grade glioma susceptibility

Battle at the p15 promoter

TGF-beta regulation by Emilin1: new links in the etiology of hypertension

Loss of EMILIN-1 enhances arteriolar myogenic tone through TGF-β (transforming growth factor-β)-dependent transactivation of EGFR (epidermal growth factor receptor) and is relevant for hypertension in mice and humans

Dissecting the role of epidermal growth factor receptor catalytic activity during liver regeneration and hepatocarcinogenesis

p <5 × 10 −8 has emerged as a standard of statistical significance for genome-wide association studies

Promoting Cardiovascular Health in the Developing World: A Critical Challenge to Achieve Global Health

The missing diversity in human genetic studies

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk

Susceptibility locus for clinical and subclinical coronary artery disease at chromosome 9p21 in the multi-ethnic ADVANCE study

The effect of chromosome 9p21 variants on cardiovascular disease may be modified by dietary intake: evidence from a case/control and a prospective study

Association of a chromosome locus 9p21.3 CDKN2B-AS1 variant rs4977574 with hypertension: the TAM-RISK study

Natural selection on genes related to cardiovascular health in high-altitude adapted Andeans

The Evolution of Life Histories

Introduction to Genetic Analysis

Trade-offs in life-history evolution

The Causes of Evolution

Selective variegated methylation of the p15 CpG island in acute myeloid leukemia

Elastin microfibril interface-located protein 1, transforming growth factor beta, and implications on cardiovascular complications

Nonclassical renin-angiotensin system and renal function

Abstract: Schematic depiction of the main observation from association data presented or discussed in this article. Two adjacent haplotype blocks characterize the 9p21.3 cardiovascular risk locus: left, the block or island containing the first part of the p15 gene and its well-studied promoter, characterized by elevated associations (red) with blood pressure (DBP, SBP) and/or hypertension in Colombian and European studies, and right, the block containing lead cardiovascular risk

1KG-EUR and 1KG-AMR, representing mainly European origin and mainly or partly Native American origin, respectively. Color-shaded matrix shows 16-SNP haplotype motifs (rows of red 0's/major allele and 1's/minor allele: master motifs, yellow shading: one-SNP mutants) of the 94 1KG-AMR individuals (188 haplotypes) from Medellín, 1KG-CLM, which correspond in ancestry and admixture to the 357 individuals M357-CLM studied here from the same city (SNP allele frequencies are shown below matrix). The 3 major haplotype classes are shaded light green (c1), light blue (c2) and sand (c3); by inclusion/exclusion they define the 3 SNP classes s1(blue rs identifiers), s2 (red rs identifiers) and s3. Relative haplotype frequencies are shown at right for 1KG-CLM, corresponding to our association results from Medellín, and for the British population

This figure illustrates a nested-plot overview of genotyped 9p21.3 SNPs, the association signals for hypertension in the Medellín study, and elements of a causal hypothesis involving a SNP in the p15 promoter. Colored bars and asterisks (top track) show extent and sentinel/lead SNP of regions of strongest BP (red) or cardiovascular risk (blue) association in the studies by Warren

We thank Professors I. King Jordan and Greg Gibson (Georgia Institute of Technology) for suggestions and references early in this study, Drs. Edwin García and Juan Ramón González for advice on association analysis, and Richard Jacobs and Jane McDougall (LGC Genomics UK) for genotyping assistance. We also thank two anonymous referees of a earlier version of this paper for suggestions. This 

 rs2069418   rs2069418   rs7023954  rs15735  rs1134871  rs3731257  rs3088440  rs11515  rs3731249  rs3731239  rs3814960  rs36228834  rs7036656*  rs2811711  rs2518723  rs3218022  rs3218020  rs2811712  rs3218018  rs3218012  rs3218005  rs3217992  rs1011970 rs2383207  rs10811650 rs1537375  rs16905599 rs1333045  rs10116277 rs10738610  rs10965227 rs10757278  rs1547705 rs1333049  rs10965228 rs1333050  rs1333040 rs10965243  rs1537370 rs2891169  rs7857345 rs2383208  rs10757274 rs7045889  rs4977574 rs10811659  rs2891168 rs10757282  rs10125231 rs10811661  rs2383206