Summary

Aging is a major risk factor for chronic disease in the human population, but there is little human data on gene expression alterations that accompany the process. We examined human peripheral blood leucocyte in-vivo RNA in a large-scale transcriptomic microarray study (subjects aged 30 to 104 years). We tested associations between probe expression intensity and advancing age (adjusting for confounding factors), initially in a discovery set (n = 458), following-up findings in a replication set (n=240). We confirmed expression of key results by real-time PCR. Of 16,571 expressed probes, only 295 (2%) were robustly associated with age. Just six probes were required for a highly efficient model for distinguishing between young and old (Area Under the Curve in replication set; 95%). The focussed nature of age-related gene expression may therefore provide potential biomarkers of aging. Similarly, only 7 of 1065 biological or metabolic pathways were age-associated, in Gene Set Enrichment Analysis (GSEA), notably including the processing of messenger RNAs (mRNAs); (p<0.002, FDR q<0.05). This is supported by our observation of age-associated disruption to the balance of alternatively-expressed isoforms for selected genes, suggesting that modification of mRNA processing may be a feature of human aging.

Altered expression of oxidative metabolism genes has been described in the skeletal muscle of individuals with type 2 diabetes. Pancreatic beta cells contain low levels of antioxidant enzymes and are particularly susceptible to oxidative stress. In this study, we explored the effect of hyperglycemia-induced oxidative stress on a panel of oxidative metabolism genes in a rodent beta cell line. We exposed INS-1 rodent beta cells to low (5.6 mmol/L), ambient (11 mmol/L), and high (28 mmol/L) glucose conditions for 48 hours. Increases in oxidative stress were measured using the fluorescent probe dihydrorhodamine 123. We then measured the expression levels of a panel of 90 oxidative metabolism genes by real-time PCR. Elevated reactive oxygen species (ROS) production was evident in INS-1 cells after 48 hours (P < 0.05). TLDA analysis revealed a significant (P < 0.05) upregulation of 16 of the 90 genes under hyperglycemic conditions, although these expression differences did not reflect differences in ROS. We conclude that although altered glycemia may influence the expression of some oxidative metabolism genes, this effect is probably not mediated by increased ROS production. The alterations to the expression of oxidative metabolism genes previously observed in human diabetic skeletal muscle do not appear to be mirrored in rodent pancreatic beta cells.

A recent genome-wide association (GWA) study of late-onset Alzheimer's disease (LOAD) identified 15 novel single nucleotide polymorphisms (SNPs) independent of ApoE. We hypothesized that variants associated with LOAD are also associated with poor cognitive function in elderly populations. We measured additive associations between the five most strongly associated LOAD SNPs and grouped Mini Mental State Examination (MMSE) scores. Variants were genotyped in respondents (mean age 79yrs) from the Oxford Healthy Aging project (OHAP) and other sites of the MRC Cognitive Function and Aging Study (MRC-CFAS). In adjusted ordinal logistic models, two variants were associated with poorer cognitive function: rs11622883 (OR=1.14, 95%CI: 1.01 to 1.28, p=0.040) and rs505058 (OR=1.29, 95% CI: 1.02 to 1.64, p=0.036). These SNPs are close to a SERPINA gene cluster and within LMNA respectively. The mechanisms underlying the associations with cognitive impairment and LOAD require further elucidation, but both genes are interesting candidates for involvement in age-related cognitive impairment.

The identification of multiple signals at individual loci could explain additional phenotypic variance (‘missing heritability’) of common traits, and help identify causal genes. We examined gene expression levels as a model trait because of the large number of strong genetic effects acting in cis. Using expression profiles from 613 individuals, we performed genome-wide single nucleotide polymorphism (SNP) analyses to identify cis-expression quantitative trait loci (eQTLs), and conditional analysis to identify second signals. We examined patterns of association when accounting for multiple SNPs at a locus and when including additional SNPs from the 1000 Genomes Project. We identified 1298 cis-eQTLs at an approximate false discovery rate 0.01, of which 118 (9%) showed evidence of a second independent signal. For this subset of 118 traits, accounting for two signals resulted in an average 31% increase in phenotypic variance explained (Wilcoxon P< 0.0001). The association of SNPs with cis gene expression could increase, stay similar or decrease in significance when accounting for linkage disequilibrium with second signals at the same locus. Pairs of SNPs increasing in significance tended to have gene expression increasing alleles on opposite haplotypes, whereas pairs of SNPs decreasing in significance tended to have gene expression increasing alleles on the same haplotypes. Adding data from the 1000 Genomes Project showed that apparently independent signals could be potentially explained by a single association signal. Our results show that accounting for multiple variants at a locus will increase the variance explained in a substantial fraction of loci, but that allelic heterogeneity will be difficult to define without resequencing loci and functional work.

BACKGROUND

It is recognized that FMR1 premutation expansions are associated with premature ovarian failure (POF), but the role of smaller repeats at the boundary of premutation and normal is less clear.

METHODS

We have therefore investigated the incidence of these intermediate sized FMR1 CGG repeats (35–58 repeats) in a series of 366 women ascertained because of menopause before the age of 40.

RESULTS

We found no significant difference in the incidence of intermediates in cases compared with controls. Thus, we were unable to replicate previous studies showing a positive association, despite a significantly larger sample size.

CONCLUSIONS

We therefore conclude that intermediate sized FMR1 CGG repeat alleles should not be considered a high-risk factor for POF based on current evidence.

Women become infertile approximately 10 years before menopause, and as more women delay childbirth into their 30s, the number of women who experience infertility is likely to increase. Tests that predict the timing of menopause would allow women to make informed reproductive decisions. Current predictors are only effective just prior to menopause, and there are no long-range indicators. Age at menopause and early menopause (EM) are highly heritable, suggesting a genetic aetiology. Recent genome-wide scans have identified four loci associated with variation in the age of normal menopause (40–60 years). We aimed to determine whether theses loci are also risk factors for EM. We tested the four menopause-associated genetic variants in a cohort of approximately 2000 women with menopause ≤45 years from the Breakthrough Generations Study (BGS). All four variants significantly increased the odds of having EM. Comparing the 4.5% of individuals with the lowest number of risk alleles (two or three) with the 3.0% with the highest number (eight risk alleles), the odds ratio was 4.1 (95% CI 2.4–7.1, P = 4.0 × 10−7). In combination, the four variants discriminated EM cases with a receiver operator characteristic area under the curve of 0.6. Four common genetic variants identified by genome-wide association studies, had a significant impact on the odds of having EM in an independent cohort from the BGS. The discriminative power is still limited, but as more variants are discovered they may be useful for predicting reproductive lifespan.

Background The human paraoxonase (PON1) protein detoxifies certain organophosphates, and the PON1 Q192R polymorphism (rs662) affects PON1 activity. Groups with higher dose exposure to organophosphate sheep dips or first Gulf War nerve toxins reported poorer health if they had 192R, and these associations have been used to exemplify Mendelian randomization analysis. However, a reported association of 192R with depression in a population-based study of older women recently cast doubt on the specificity of the higher dose findings. We aimed to examine associations between the PON1 Q192R polymorphism and self-reported poor health and depression in two independent population-based studies.

Methods We used logistic regression models to examine the associations in men and women aged 60–79 years from the English Longitudinal Study of Ageing (ELSA, n = 3158) and InCHIANTI (n = 761) population studies. Outcomes included the Center for Epidemiologic Studies Depression (CES-D) scale, self-rated general health status and (in ELSA only) diagnoses of depression.

Results The PON1 Q192R polymorphism was not associated with self-reported poor health {meta-analysis: odds ratio (OR) = 1.01 [confidence interval (CI) 0.91–1.13], P = 0.80} or depressive symptoms in either study or in meta-analyses [CES-D: OR = 1.01 (CI 0.87–1.17), P = 0.90]. There was also no association with histories of diagnosed depression in ELSA [OR = 1.03 (CI 0.82–1.30), P = 0.80].

Conclusions We found no evidence of an association between the PON1 Q192R polymorphism and poor general or mental health in two independent population-based studies. Neither the claimed Q192R association with depression in the general population nor its theoretical implications were supported.

We conducted a meta-analysis of genome-wide association data to detect genes influencing age at menarche in 17,510 women. The strongest signal was at 9q31.2 (P = 1.7 × 10−9), where the nearest genes include TMEM38B, FKTN, FSD1L, TAL2 and ZNF462. The next best signal was near the LIN28B gene (rs7759938; P = 7.0 × 10−9), which also influences adult height. We provide the first evidence for common genetic variants influencing female sexual maturation.

Background

A common variant at chromosome 9p21 (tagged by the rs1333049 or rs10757278 SNP) is strongly associated with Myocardial Infarction (MI) and major arterial aneurysms. An association with Peripheral Arterial Disease (PAD) was also reported in a sample aged <75 years, but this disappeared on removal of respondents with a MI history, resulting in an odds ratio for PAD of 1.09 (p=0.075). We aimed to estimate the association of this variant with Ankle Brachial Index (ABI) and PAD in three older populations.

Methods and Results

We used data from the InCHIANTI, Baltimore Longitudinal Study of Aging and Health, Aging and Body Composition studies. In 2,630 Caucasian individuals (mean age 76.4 years) the C allele at rs1333049 was associated with lower mean ABI measures and with increased prevalence of PAD. These associations remained after removal of baseline and incident MI cases over a 6 year follow-up for both ABI (−0.017 ABI units, 95% CI: −0.03- −0.01, p=1.3×10−4) and PAD (per allele OR: 1.29, 95% CI: 1.06–1.56, p=0.012). These associations also remained after adjustment for known atherosclerosis risk factors including Diabetes Mellitus, smoking, hypercholesterolemia and hypertension.

Conclusions

The C allele at rs1333049 is associated with an increased prevalence of Peripheral Arterial Disease and lower mean Ankle Brachial Index. This association was independent of the presence of diagnosed MI and atherosclerotic risk factors in 3 older Caucasian populations.

Aims

There are a large number of common genetic variants that have been robustly associated with low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, or triglyceride concentrations. The majority of these have been identified or confirmed in recent genome-wide association studies, but few studies have assessed the combined effect of known lipid variants. We hypothesized that these variants would influence both the need for interventions and myocardial infarction (MI) outcomes. We aimed to estimate combined effects of proven SNPs on LDL, HDL, and triglyceride concentrations and MI history in a representative older population.

Methods and results

In the InCHIANTI Study of Aging (age ≥65 years), we calculated individual dyslipidaemia risk allele counts for increased LDL (range 4–14, n = 594), reduced HDL (5–16, n = 635), and increased triglycerides (7–16, n = 611). Lipid levels were compared with ATPIII National Cholesterol Education Panel (NCEP) intervention guidelines. Individual variants and the APOE haplotype explained <2.1% of the variance in their respective lipid concentrations, with the exception of the CETP SNP rs1800775 and HDL levels (4.76%). Combined risk allele counts outperformed the largest single-SNP effects for LDL (explaining 7.1% of variance) and triglycerides (4.8%), but not HDL (3.4%). Risk alleles were divided as near as possible into quartiles. The 31% of respondents with 10 or more LDL increasing alleles were more likely to have LDL levels above the intervention threshold (OR 3.00, 95% CI 1.67–5.39, P = 2.5 × 10−4), compared with the 21% with 7 or less risk alleles. Similarly, the 35% with 13 or more triglyceride risk alleles were more likely to exceed NCEP intervention thresholds (OR 2.98, 95% CI 1.43–6.22, P = 0.004) compared with the 24% with 10 or less alleles. The number of individuals reporting an MI event was small (n = 67), but an event was more common in the 36% of respondents who had the highest combined risk allele score for all three lipids (OR 3.68, 95% CI 1.21–11.2, P = 0.021) compared with the lowest risk 22%.

Conclusion

In a representative older population, the cumulative effects of proven LDL- and triglyceride-altering genetic variants increased the odds of crossing the lipid-level threshold for therapeutic intervention by approximately three-fold.

Background

The proinflammatory cytokine interleukin-18 (IL-18) is associated with major disabling conditions, although whether as byproduct or driver is unclear. The role of common variation in the IL-18 gene on serum concentrations and functioning in old age is unknown.

Methods

We used 1671 participants aged 65–80 years from two studies: the InCHIANTI study and wave 6 of the Iowa-Established Populations for Epidemiological Study of the Elderly (EPESE). We tested three common polymorphisms against IL-18 concentration and measures of functioning.

Results

In the InCHIANTI study, a 1 standard deviation increase in serum IL-18 concentrations was associated with an increased chance of being in the 20% of slowest walkers (odds ratio 1.45; 95% confidence interval, 1.17–1.80; p = .0007) and 20% of those with poorest function based on the Short Physical Performance Battery Score (odds ratio 1.52; 95% confidence interval, 1.22–1.89; p =.00016) in age sex adjusted logistic regression models. There was no association with Activities of Daily Living (p = .26) or Mini-Mental State Examination score (p = .66). The C allele of the IL-18 polymorphism rs5744256 reduced serum concentrations of IL-18 by 39 pmol/mL per allele (p = .00001). The rs5744256 single nucleotide polymorphism was also associated with shorter walk times in InCHIANTI (n =662, p =.016) and Iowa-EPESE (n =995, p =.026). In pooled ranked models rs5744256 was also associated with higher SPPB scores (n =1671, p = .019). Instead of adjusting for confounders in the IL-18 walk time association, we used rs5744256 in a Mendelian randomization analysis: The association remained in instrumental variable models (p = .021).

Conclusion

IL-18 concentrations are associated with physical function in 65- to 80-year-olds. A polymorphism in the IL-18 gene alters IL-18 concentrations and is associated with an improvement in walk speed. IL-18 may play an active role in age-related functional impairment, but these findings need independent replication.

There is considerable evidence that human genetic variation influences gene expression. Genome-wide studies have revealed that mRNA levels are associated with genetic variation in or close to the gene coding for those mRNA transcripts – cis effects, and elsewhere in the genome – trans effects. The role of genetic variation in determining protein levels has not been systematically assessed. Using a genome-wide association approach we show that common genetic variation influences levels of clinically relevant proteins in human serum and plasma. We evaluated the role of 496,032 polymorphisms on levels of 42 proteins measured in 1200 fasting individuals from the population based InCHIANTI study. Proteins included insulin, several interleukins, adipokines, chemokines, and liver function markers that are implicated in many common diseases including metabolic, inflammatory, and infectious conditions. We identified eight Cis effects, including variants in or near the IL6R (p = 1.8×10−57), CCL4L1 (p = 3.9×10−21), IL18 (p = 6.8×10−13), LPA (p = 4.4×10−10), GGT1 (p = 1.5×10−7), SHBG (p = 3.1×10−7), CRP (p = 6.4×10−6) and IL1RN (p = 7.3×10−6) genes, all associated with their respective protein products with effect sizes ranging from 0.19 to 0.69 standard deviations per allele. Mechanisms implicated include altered rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA), variation in gene copy number (CCL4L1) and altered transcription (GGT1). We identified one novel trans effect that was an association between ABO blood group and tumour necrosis factor alpha (TNF-alpha) levels (p = 6.8×10−40), but this finding was not present when TNF-alpha was measured using a different assay , or in a second study, suggesting an assay-specific association. Our results show that protein levels share some of the features of the genetics of gene expression. These include the presence of strong genetic effects in cis locations. The identification of protein quantitative trait loci (pQTLs) may be a powerful complementary method of improving our understanding of disease pathways.

Author Summary

One of the central dogmas of molecular genetics is that DNA is transcribed to RNA which is translated to protein and alterations to proteins can influence human diseases. Genome-wide association studies have recently revealed many new DNA variants that influence human diseases. To complement these efforts, several genome-wide studies have established that DNA variation influences mRNA expression levels. Loci influencing mRNA levels have been termed “eQTLs”. In this study we have performed the first genome-wide association study of the third piece in this jigsaw – the role of DNA variation in relation to protein levels, or “pQTLs”. We analysed 42 proteins measured in blood fractions from the InCHIANTI study. We identified eight cis effects including common variants in or near the IL6R, CCL4, IL18, LPA, GGT1, SHBG, CRP and IL1RN genes, all associated with blood levels of their respective protein products. Mechanisms implicated included altered transcription (GGT1) but also rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA) and variation in gene copy number (CCL4). Blood levels of many of these proteins are correlated with human diseases and the identification of “pQTLs” may in turn help our understanding of disease.

p16INK4a is active in cell senescence, ageing and tumor suppression. Deletion of the small p16INK4a / ARF / p15INK4b region occurs in many cancers. We screened 25 common polymorphisms across the region and 3 related genes for associations with physical functioning in older people.

In an initial sample of 938 (aged 65 to 80yrs) from the EPIC study (Norfolk, UK) the rs2811712 SNP minor allele (located between the shared p16INK4a / ARF locus and p15INK4b) was associated with reduced physical impairment. This association remained after testing an additional 1319 EPIC-Norfolk samples (p-value=0.013, total n=2257), and on independent replication in the InCHIANTI Study (n=709, p=0.015), and at one sided significance in Iowa-EPESE (n=419, p=0.079). Overall (n=3372) the prevalence of severely limited physical function was 15.0% in common homozygotes and 7.0% in rare homozygotes (per minor allele Odds Ratio=1.48 95%CI: 1.17−1.88, p =0.001, adjusted for age, sex and study). This estimate was similar excluding screening set 1 (OR=1.45: 95%CI 1.09−1.92, p=0.010, n=2434).

These findings require further replication, but provide the first direct evidence that the p16INK4a / ARF / p15INK4b genetic region and the senescence machinery are active in physical ageing in heterogeneous human populations. The mechanism involved may be via greater cellular restorative activity and reduced stem cell senescence.

Background

A polymorphism in the transcription factor 7-like 2 (TCF7L2) gene has been found to be associated with type 2 diabetes in case-control studies. We aimed to estimate associations of the marker rs7903146 (C/T) polymorphism with fasting glucose, lipids, diabetes prevalence and complications in an older general population.

Methods

In total, 944 subjects aged ≥ 65 years from the population representative InCHIANTI study were enrolled in this study. Those with fasting blood glucose of ≥ 7 mmol/l or physician diagnosis were considered diabetic. Cut-off points for impaired fasting glucose (IFG) were ≥ 5.6 mmol/l to < 7 mmol/l.

Results

In the general population sample, minor (T) allele carriers of rs7903146 had higher fasting blood glucose (FBG) (p = 0.028) but lower fasting insulin (p = 0.030) and HOMA2b scores (p = 0.001), suggesting poorer beta-cell function. T allele carriers also had smaller waist circumference (p = 0.009), lower triglyceride levels (p = 0.006), and higher high-density lipoprotein cholesterol (p = 0.008).

The prevalence of diabetes or IFG was 32.4% in TT carriers and 23.3% in CC carriers; adjusted OR = 1.67 (95% confidence interval 1.05 to 2.65, p = 0.031). Within the diabetic and IFG groups, fewer T allele carriers had metabolic syndrome features (p = 0.047) or had experienced a myocardial infarction (p = 0.037). Conversely, T allele carriers with diabetes had poorer renal function (reduced 24-hour creatinine clearance, p = 0.013), and possibly more retinopathy (p = 0.067). Physician-diagnosed dementia was more common in the T carriers (in diabetes p = 0.05, with IFG p = 0.024).

Conclusion

The TCF7L2 rs7903146 polymorphism is associated with lower insulin levels, smaller waist circumference, and lower risk lipid profiles in the general elderly population. Patients with diabetes who are carriers of the minor allele are less likely to have metabolic-syndrome features, but may experience more microvascular complications, although the number of cases was small. If replicated, these findings may have implications for developing treatment approaches tailored by genotype.