Identification of single nucleotide polymorphisms in the p21 (CDKN1A) gene and correlations with longevity in the Italian population

Introduction

Longevity in humans is believed to be a multifactorial condition to which both genetic and environmental factors are likely to contribute. Centenarians are often spared from age-related diseases, such as cardiovascular disease, Alzheimer's disease, diabetes mellitus, and cancer. The rate of aging and maximum lifespan varies among species, and therefore it has been postulated to be at least in part under genetic control [1-2]. Epidemiological data indicate the presence of a strong familiar component of longevity that is largely determined by genetics. Thus, progeroid syndromes of accelerated aging have known genetic causes [3-4]. A number of possible associations between longevity and allelic variants of genes have been described. Estimates of the heritability of human lifespan vary from 10-50% with the most common finding being that about a third of human lifespan may be heritable. The rest is due to environmental exposure, accidents and injuries, lifestyle and chance. Very long life, to beyond the age of 90 years, appears to have an even stronger genetic basis [5], which explains why centenarians and near-centenarians tend to cluster in families.

Theories on aging postulate that aging is a remodeling process, where the body of survivors progressively adapts to internal and external damaging agents, to which they are exposed over several decades. Thus, stress response and adaptation mechanisms play a fundamental role in the aging process and have an impact on individual lifespan. Centenarians' capability to live such extraordinarily long lives is in large part due to genetic variations that either affect the rate of aging or decrease the susceptibility to age-associated diseases.

Some of the most promising candidate genes appear to be those involved in stress response. An interesting possible candidate is p21 (CDKN1A) which has been shown to be involved both in stress response mechanisms and in the expression of genes implicated in age-related diseases. p21 is best known as a stress-inducible cyclin-dependent kinase inhibitor, which triggers cell growth arrest associated with senescence and damage response. Some evidence suggests that the effects of p21 inductionon gene expression in senescent cells may contribute to the pathogenesisof cancer and age-relateddiseases. In particular, p21 expression was found to upregulatemultiple genes that have not only been associated with senescence but also implicatedin age-related diseases, including Alzheimer'sdisease, atherosclerosis, amyloidosis, arthritis and cancer, thus suggesting that p21 induction by stress may play a causal role in these diseases [6]. The role of p21 in cell senescence and its possible implication in the risk ofage-related diseases suggests that allelic variations in this gene may have an impact on the lifespan. The goal of this study was to identify p21 polymorphisms and to determine whether they may be associated with longevity.

Results

Discussion

In the present study, we have investigated possible associations between longevity and SNPs in the p21 (CDKN1A) gene, which plays a role in stress response and cell senescence, and increased expression of which was shown to upregulate genes implicated in several age-related diseases [6]. By sequencing the three exons and 5 kb of the promoter region of p21, we have identified many previously unknown or uncharacterized SNPs in the p21 promoter. We have tested the activity of the promoter derived from an LLI and containing a cluster of minor alleles in the region between -4547 bp and -3489 bp, where a novel p53 binding site has been recently identified [7] and found no changes in the basal activity or p53 responsiveness of this promoter. It should be noted, however, that there are many p53 independent physiological signals that induce p21, where the response of the two alleles may potentially be different. Interestingly the p21 promoter is induced by some signals involved in stress response and inflammation (such as TGFβ, INFγ, IL-6) that, as discussed elsewhere, contribute to the pathogenesis of many age-related diseases [9-11].

In a large case-control design study, we have compared the frequency of 41 SNPs spanning the p21 gene between large populations of LLI and younger Italian individuals. Only two of 41 SNPs showed a statistically significant difference between the two populations, and remarkably, these were the only two exon-derived SNPs. A 6-SNP haplotype comprising the common alleles of these two SNPs was strongly overrepresented among the centenarians relative to the control population. One of these exonic SNPs (rs1801270) changes the amino acid sequence of p21 from Ser to Arg at codon 31, and the other (rs1059234) leads to a C->T transition 20 nucleotides downstream of the stop codon in the 3' untranslated region. Remarkably, several studies sug-gested that the rare alleles of these two SNPs may play a role in different types of cancer [12-17]. In particular, Li et al. [12] have shown that, in non-Hispanic whites, the rare versions of the rs1801270 and rs1059234 alleles are associated to an increased risk susceptibility to squamous cell carcinoma, individually and in combination. In addition, Mousses at al. [18] observed that the rare alleles of these two SNPs were under-represented in breast cancer and sarcoma patients whose tumors possessed somatic p53 mutations, as compared to tumors without p53 mutations, suggesting that these alleles could influence p21 functions in a p53-independent manner.

According to the data in the NCBI database (Entrez SNP), the allele frequency distribution of the SNPs rs1801270 and rs1059234 is highly variable and ethnic-specific. The frequency of the rs1801270 minor allele varies from 0.021 (Europeans) to 0.47 (Asians) with African-Americans, Sub-Saharan Africans, and Hispanics having intermediate values. Similar values are found for the rs1059234 rare allele, whose frequencies vary from 0.021 (Europeans) to 0.45 (Asians), with African-Americans, Sub-Saharan Africans, and Hispanics also having intermediate values. This high ethnicity-related variability can greatly complicate the interpretation of disease-associated studies, especially those conducted in multi-ethnic societies. In our large-scale study, we have used a relatively homogenous Italian population. Our results suggest that the presence of the rs1801270 and rs1059234 rare alleles combined may be detrimental to longevity and therefore negatively selected in Italian LLI. Further large-scale studies could be useful to compensate for genetic heterogeneity within the Italian population and to clarify the potential role of these SNPs in limiting the lifespan.

The mechanisms underlying the potential detrimental effect of the rare alleles of rs1801270 and rs1059234 are presently unknown. The obvious hypotheses are that the amino acid change at codon 31, which was proposed to abolish p21 phosphorylation at Ser 31 [19], could modulate its abilities to arrest the cell cycle or to induce transcription of genes implicated in age-related diseases, and that a nucleotide change in the 3' UTR could affect p21 mRNA stability or translational efficiency. These possibilities remain to be tested in future studies.

Longevity in humans can be defined as a multifactorial condition to which both genetic and environmental factors are likely to contribute. Twin studies have shown that genetic differences account for about a quarter of the variance in adult human lifespan. Despite the challenges of studying complex traits such as lifespan, studies have been reporting alleles that were significantly associated with human longevity. One of the best examples is APOE whose association has been reproduced consistently [20-22]. The compression of morbidity hypothesis proposed by James Fries in 1980 [23] postulates that as the limit of human lifespan is approached, the onset and duration of lethal impairment compresses toward the end of life. This ‘compression' is observed in the majority of centenarians who are often spared from age-related diseases, specifically cardiovascular disease, Alzheimer disease, diabetes mellitus, and cancer. To achieve their extreme age, centenarians likely lack numerous gene variants that are associated with age-related diseases and they may be more likely to carry protective variants as well. Our finding that the frequency of specific minor alleles of p21 is decreased among Italian centenarians lends additional support to this concept.

Materials and Methods

Subjects. In the pilot study we used ten DNA samples from LLI over 90 y.o. (white Americans of European descent), including five females and five males. The DNA samples, obtained from the NIA Aging Cell Repository DNA panel, were obtained from CORIELL bank (Camden, NJ). The second group of samples comprised DNA from 92 non-centenarian subjects, belonging to Utah/CEPH population, provided by CORIELL bank. These samples were initially collected from Utah residents with ancestry from northern and western Europe. In the large case-control study, 184 Italians with exceptional longevity (mean age, 100.88±1.77 years) and a control group (38.97±12.21) were recruited by the Bologna group in Central Italy, after checking for ethnicity and ancestor origins. The sex ratio in the Italian samples was 7 female to 1 male in the centenarian group and 2:1 in the control group. A full socioeconomic, quality of life and health status assessment was performed. Primary criteria for inclusion in the study were good health (for centenarians, subjects categorized A or B according to Franceschi et al. [24] were included), physical activity and absence of major diseases.

p21 genotyping. We sequenced the three exons (68 bp, 450 bp, 1600 bp) and a 5 kb promoter region of the p21 (CDKN1A) gene in the DNA of ten LLI. The exons and the promoter were amplified by PCR in overlapping fragments of ~ 400 bp each. The PCR primer pairs are listed in Table 4, where p21pro.1 is the furthest from the transcription start site. Primers were designed using the primer3 software available at http://frodo.wi.mit.edu. The primers were sized between 22-24 bases with a Tm of 69-71^oC and a GC content of 40-60 %. The primers were checked for loops, hairpins and 3' complementarity. The selected primers were synthesized by idtDNA (Coralville, IA). The genomic DNA templates were added to a master mix containing 2 μl of each primer (10 μM), 5 μl Buffer (New England Biolabs, Ipswich, MA), 2 μl of Taq polymerase (Invitrogen), and water to the volume of 50 μl. The following PCR profile was used: preincubation for 2 min at 96°C, 40 cycles of 30 sec at 95°C, 30 sec at the primer-specific annealing temperature (Table 4) and 15 sec at 72°C, followed by final incubation for 5 min at 72°C. Once amplified, the fragments were purified with Millipore columns (Millipore Billerica, MA) and then sequenced. Sequencing was performed with ABI 3730 DNA analyzer, using a Big Dye protocol with Zymo column-purified products. Complete sequences were aligned, assembled and compared using the Clone Manager program. For verification, visual inspection of sequence profiles for each candidate SNP was carried out. At least two overlapping DNA templates amplified with different primers were used for identification of each candidate SNP. In addition to the SNP candidate approach, the latest data available on the HapMap were analyzed, in order to select an appropriate number of tagSNPs to cover at least 90% of the genetic information in the locus (r2=0.9 and MAF>0.05). The tagSNP and candidate SNP sets were used for genotyping.

The SNP genotyping was performed using Sequenom's chip-based matrix assisted laser desorption/ionization time-of-flight MS (DNA MASSARRAY). This technology performs allele-specific single-base primer extension reactions, which allow differentiation of homozygous normal, heterozygous mutant and homozygous mutant samples (iPLEX assay).  The MassEXTEND primers anneal up to the polymorphic site and are extended with one single base. The allele product masses depend on the SNP allele base and thus they are easily distinguished with the mass spectrometer. DNA of the Utah/CEPH population was analyzed as a service by Sequenom (San Diego, CA) on PCR-derived extension products from individual DNA samples. Large-scale analysis of Italian samples was performed by the Bologna group. Cases and controls were always analyzed on the same chip to avoid potential artifacts caused by chip-specific miscalls.

Plasmid constructs. The plasmid p21-PGL.4.10-luc containing the 5 kb p21 promoter comprising the common alleles and driving firefly luciferase expression has been previously described [7]. This plasmid was used to replace a cluster of common SNP alleles in the promoter with the minor alleles, contained within a fragment of ~2.1 kb, which was amplified by double-round PCR from genomic DNA of a LLI. PCR was carried out using a proofreading polymerase, Phusion™ Hot Start High-Fidelity DNA Polymerase (New England Biolabs, Ipswich, MA). The primed template was pre-formed in the presence of 5X Phusion GC buffer (New England Biolabs, Ipswich, MA) and 200 μM of each dNTP, 0.5 μM primers, and 1 U of Phusion DNA polymerase. The template being GC-rich, 3% DMSO was added to optimize the product yield. The samples were incubated as follows: preincubation for 30 sec at 98°C, 30 cycles of 10 sec at 98°C, 30 sec at 66°C and 3 min at 72°C and one final incubation for 5 min at 72°C. In the first round, the 5kb PCR product of the p21 promoter was amplified using the following primers: p21-4997F TACAAACATTGGGTGGGGCGAGTC p21-R-44 CTCCGGCTCCACAAGGAACTGACTT In the second round, this PCR product was used as a template to generate a PCR product of ~2.1 kb using the following primers: p21-4497F TACAAACATTGGGTGGGG CGAGTC p21-5R GACAAAATAGCCACCAGCCTCTTCT The latter PCR product was digested with AatII and Sph restriction enzymes (New England Biolabs), and cloned into p21-PGL.4.10-luc plasmid digested with the same enzymes, replacing the corresponding fragment containing the common alleles. The resulting plasmid was sequence-verified and designated p21R-luc.

Promoter analysis by transient transfection. HCT116 colon carcinoma cells, both wild type and p53-/- sublines [8] (a gift of Dr. B. Vogelstein, Johns Hopkins University) were grown in DMEM with Earle's salts supplemented with 10% FCS and 2 mM L-glutamine in a humidified 95% air 5% CO2 incubator. Cells were seeded in 12-well tissue culture plates for 24 h prior to transfection. When 70% confluent, the cells were transfected with 1 μg of the indicated promoter-reporter plasmids, together with pRL-TK Renilla luciferase expressing plasmid (Promega, Madison, WI) to normalize for transfection efficiency, at a ratio of 10:1 test vector:standard vector. Transfections were performed in triplicate, using FuGENE6 (Roche Molecular Biochemicals). A precipitate was formed using 3 μl of FuGENE6/μg of transfected DNA and the transfection mixture was diluted up to 1 ml with serum-free medium. After incubation at 22°C for 10 min, the DNA/FuGENE6 mixture was added to cells. Cells were harvested 48 h after transfection, and firefly and Renilla luciferase activities were measured.

Statistical analysis. We tested departurefrom Hardy-Weinberg equilibrium [25] in the controls by a ² testusing P = 0.01 as threshold. This threshold was chosen basedon anticonservativeness of this test as noted by Wigginton [25]. All SNPs (except for rs4711458, rs4714003, CDKN1A7, rs6920453) were in Hardy-Weinbergequilibrium. Each SNP was tested both with basic association testing based on comparing allele frequencies between cases and controls (asymptotic and empirical p-value to control for multiple testing), and with Conchran-Armitage trend test in a dominant, recessive and general model.

Multi-locus haplotype analysis was performed by using a sliding-window approach implemented in PLINK [26] for multi-loci of 4 or 6 SNPs size. Multimarker haplo- types have been estimated by using the E-M algorithm implemented in the software. An LD map of the region has been produced by using the software Haploview [27] in order to identify LD blocks within the typed markers. The haplotype blocks identified have been tested in the same way, i.e. phased and used for frequency estimation by E-M algorithm and chi-square testing of the frequencies.