Research Paper Volume 11, Issue 15 pp 5518—5534
Frailty in middle age is associated with frailty status and race-specific changes to the transcriptome
- 1 Laboratory of Epidemiology and Population Sciences, National Institute on Aging, National Institutes of Health, Baltimore, MD 21224, USA
- 2 Laboratory of Genetics and Genomics, National Institute on Aging, National Institutes of Health, Baltimore, MD 21224, USA
Received: March 5, 2019 Accepted: July 27, 2019 Published: August 8, 2019https://doi.org/10.18632/aging.102135
How to Cite
Copyright © 2019 Prince et al. This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY) 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Frailty is an aging-associated syndrome resulting from diminished capacity to respond to stressors and is a significant risk factor for disability and mortality. Although frailty is usually studied in old age, it is present in mid-life. Given the increases in mortality statistics among middle-aged Americans, understanding molecular drivers of frailty in a younger, diverse cohort may facilitate identifying pathways for early intervention. We analyzed frailty-associated, genome-wide transcriptional changes in middle-aged blacks and whites. Next generation RNA sequencing was completed using total RNA from peripheral blood mononuclear cells (n = 16). We analyzed differential gene expression patterns and completed a parametric analysis of gene set enrichment (PAGE). Differential gene expression was validated using RT-qPCR (n = 52). We identified 5,082 genes differentially expressed with frailty. Frailty altered gene expression patterns and biological pathways differently in blacks and whites, including pathways related to inflammation and immunity. The validation study showed a significant two-way interaction between frailty, race, and expression of the cytokine IL1B and the transcription factor EGR1. The glucose transporter, SLC2A6, the neutrophil receptor, FCGR3B, and the accessory protein, C17orf56, were decreased with frailty. These results suggest that there may be demographic dependent, divergent biological pathways underlying frailty in middle-aged adults.
Frailty is a clinical syndrome that has been defined by several different assessments as a loss of physiological reserve, thus hindering recovery from endogenous and exogenous stressors [1–3]. Frail individuals have increased risk for several adverse outcomes including: falls, disability, complications following surgery, and premature mortality [4–7]. The most recent data finds that the overall prevalence of frailty in adults 65 years old and older is approximately 15% . In the United States, the prevalence of frailty in the elderly is highest in women, racial/ethnic minorities, and individuals at lower incomes [8–10].
While frailty has been primarily examined as a clinical entity of the aged, frailty can occur in adults 45 – 64 years old. We have previously shown in a longitudinal cohort study of blacks and whites in Baltimore that frailty prevalence among these individuals is approximately 7% in adults 35-44 years old and over 13% in whites between the ages of 45-54 years old . Other studies in young to middle-aged diverse cohorts, show that frailty, grip strength, and low appendicular skeletal muscle mass are associated with poor health outcomes [10–12]. Frailty, because of its link to risk for mortality, is an increasingly important entity to study in younger cohorts in light of the increase mid-life mortality among Americans across racial and ethnic groups . With an expected increase in the US 45 – 64 year old population by 15% by 2060, it is imperative to identify functional biological markers of frailty to identify potential interventions that could aid in ameliorating the phenotype.
Identifying a single operational definition for frailty in a population has, thus far, resulted in close to 30 different assessments . The Cardiovascular Health Study (CHS)/Fried Index  is a popular and well-validated assessment that characterizes frailty by phenotype based on the presence of: shrinking (unintentional weight loss of > 10 pounds by direct measurement of weight), weakness, poor endurance and energy (self reported), slowness (time to walk 15 feet), and low physical activity (weighted score of kilocalories expended per week). Although seminal to the field of frailty research, diagnosing frailty using the CHS frail phenotype requires an in person physical with techniques that are not typically available to primary care physicians. The International Academy on Nutrition and Aging FRAIL scale is an easily administered questionnaire  that has been validated in community based, racially diverse populations consisting of middle-aged adults [9,10,15]. Most clinical and molecular studies assessing frailty utilizing the CHS index and have focused on older, largely white cohorts.
Complementary to diagnosing frailty based on phenotypic changes is to identify biomarkers for frailty. Potential dysregulated pathways associated with frailty include: inflammation, apoptosis, calcium homeostasis, neuronal signaling (central nervous system and peripheral nervous system), hormone regulation, and gene expression/transcriptional regulation . Although characterized as a chronic state of inflammation , circulating cytokine concentrations vary amongst frail populations which is a reflection of the complexity of the syndrome that could be attributed to genetic variability and participant demographics . Thus, examining gene expression profiles can serve as a promising molecular pathway to investigate frailty. The Vitality 90+ study , a longitudinal study that focuses on longevity in the oldest old, is the largest study to date to examine frailty-associated changes to the transcriptome of nonagenarians. However, the focus of this study was to examine transcriptional changes in the context of survival and mortality and was not necessarily focused on frailty-associated changes .
Previous studies do not examine whether there are differences in the underlying biological processes associated with the frail phenotype in the context of race and younger age groups. In this study, we have chosen to examine frailty-associated, genome-wide transcriptional changes in frail and non-frail middle-aged blacks and whites. This will permit us to either confirm that racially diverse middle-aged individuals with frailty have similar gene expression profiles as older European ancestry adults from previously published work  or to identify novel gene expression profiles that may be specific for younger or non-European ancestry individuals. We hypothesize that gene expression is influenced by frailty status and demographic characteristics including race.
Frailty-associated global changes to gene expression
The overall study design and RNA sequencing and analysis workflow are shown in Figure 1. We classified frailty in blacks and whites in the RNA sequencing and validation cohorts (Figure 1A) utilizing the FRAIL scale. To identify genome-wide transcriptional changes associated with frailty in middle-aged adults, we isolated total RNA from peripheral blood mononuclear cells (PBMCs) from 16 non-frail, pre-frail/frail blacks and whites between the ages of 45-49 years old. (Table 1A). We performed next generation RNA sequencing and used TopHat to align the sequences to human genome NCBI build Version 38 (Figure 1B). The percentage of mapped reads ranged between 66-92% and the total mapped reads ranged between 14 and 44 million. RNA-seq data was analyzed as shown in Figure 1B.
Figure 1. Workflow and schematic representation of RNA sequencing and RT-qPCR analysis. Overview of study design for RNA sequencing (n=16) and RT-qPCR in validation cohort (n=52) (A), and RNA sequencing analysis for alignment and identification of frailty-associated genes (B).
Table 1. Demographic information for frailty cohorts.
|A||RNA sequencing cohort|
|Age (mean (sd))||48.09 (1.21)||47.85 (1.84)|
|Sex = Women (%)||4 (50.0)||4 (50.0)|
|Race = White (%)||4 (50.0)||4 (50.0)|
|PovStat = Below (%)||2 (25.0)||4 (50.0)|
|Age (mean (sd))||47.51 (1.59)||47.23 (1.26)||47.04 (1.62)||47.58 (1.80)|
|Sex = Women (%)||7 (53.8)||6 (46.2)||7 (53.8)||6 (46.2)|
|PovStat = Below (%)||6 (46.2)||3 (23.1)||8 (61.5)||5 (38.5)|
Frailty was associated with differential expression of 5,028 genes (Figure 2A; Supplementary Table 1). When stratifying by race, as displayed in the Venn diagrams, 124 genes in blacks, 396 in whites, and 109 overlapping genes between blacks and whites were significantly and differentially expressed in the context of frailty (Figure 2B). Genes significantly different between non-frail and frail blacks and whites are listed in Supplementary Table 2 and Supplementary Table 3, respectively. When comparing the number of genes that were higher in abundance in frail individuals, we identified 139 genes in blacks, 53 genes in whites, and 18 overlapping genes between blacks and whites (Figure 2C). When comparing the number of genes with lower abundance in frail individuals, we identified 51 genes in blacks, 409 genes in whites, and 25 overlapping in blacks and whites (Figure 2D).
Figure 2. Global gene expression changes by frailty and race. Total RNA was isolated from PBMCs from non-frail and frail blacks and whites (n=16). Differential gene expression was assessed using RNA sequencing. The graph shows log2 fold change of genes significantly altered with frailty (5,082 genes) (A). A list of all these genes can be found in Supplementary Table 1. Venn diagrams of the total number of significant differentially expressed genes in blacks and whites with frailty (B). Significantly increased (C; up arrow) and decreased (D; down arrow) genes in blacks and whites with frailty.
To identify pathways associated with frailty in blacks and whites, we used the quantile normalized Fragments Per Kilobase of transcript per Million mapped reads (FPKM) to calculate Z-scores [19,20] (Figure 1B). The Z-scores were transformed (natural logarithm) and used to calculate Z-ratios that were imputed into PAGE analysis to identify relevant biological processes. We identified enriched gene ontology (GO) sets associated with biological processes comparing non-frail and frail blacks (Figure 3) and whites (Figure 4). Importantly, there were fewer pathways identified in blacks when compared to whites. The most enriched pathways in frail blacks were fever, immune response, inflammatory response, chemotaxis, and negative regulation of hormone secretion. The most downregulated pathways in frail blacks were oxygen transport, phosphate transport, sensory perception of taste, cell matrix adhesion, and induction of an organ. Conversely, in frail whites, the most enriched pathways were sensory perception of taste, tRNA processing, detection of chemical stimulus involved in sensory perception of bitter taste, response to acidic pH, and mRNA transport. The most downregulated pathways in frail whites were inflammatory response, chemotaxis, immune response, cell-cell signaling, and positive regulation of ossification. We highlighted GO terms in red that were associated with inflammation, as the pathway Z-scores for frailty varied greatly between blacks and whites (Table 2). This suggests that inflammatory pathways could be a driver for racial differences in frailty. There were 19 overlapping pathways comparing blacks to whites, however most of the pathway Z-scores were in opposite directions, only 2 were in the same direction: chemotaxis and cell matrix adhesion (Table 2). The complete list of biological processes, stratified by race, can be found in Supplementary Table 4 and Supplementary Table 5.
Figure 3. Frailty-associated biological pathways in blacks. Differentially expressed genes between non-frail and frail blacks were imputed into Parametric Analysis of Gene Set Enrichment (PAGE) analysis. Significantly changed gene sets associated with biological processes, organized by Z-score, are shown here. GO terms associated with inflammation and immune response are highlighted in red.
Figure 4. Frailty-associated biological pathways in whites. Differentially expressed genes between non-frail and frail whites were imputed into Parametric Analysis of Gene Set Enrichment (PAGE) analysis. Significantly changed gene sets, associated with biological processes, organized by Z-score, are shown here. GO terms associated with inflammation and immune response are highlighted in red.
Table 2. Frailty-associated gene ontology terms with overlapping relationships in blacks and whites.
|blacks||Gene Ontology Term||whites|
|6.29||INTERLEUKIN 1 RECEPTOR BINDING||-4.98|
|4.42||NEGATIVE REGULATION OF HORMONE SECRETION||-4.18|
|4.34||CELL CELL SIGNALING||-5.41|
|3.69||HUMORAL IMMUNE RESPONSE||-2.87|
|3.16||REGULATION OF PROGRESSION THROUGH CELL CYCLE||-2.98|
|2.54||POSITIVE REGULATION OF CELL PROLIFERATION||-3.46|
|2.17||T HELPER 1 TYPE IMMUNE RESPONSE||-3.21|
|-1.84||CELL MATRIX ADHESION||-2.38|
|-2.42||SENSORY PERCEPTION OF TASTE||2.62|
|Gene ontology Z-scores for non-frail versus frail blacks or whites are indicated.|
Interaction network analysis
To visualize functional connections between frailty-associated genes, we used protein-protein interaction analysis as a complement to the PAGE analysis . The significant, differentially expressed, frailty-associated genes referenced in Figure 2B were imputed into the STRING database. In blacks, we identified 2 functional clusters, inflammation and signaling: G protein coupled receptor (GPCR) and transcription factor (TF) signaling (Figure 5A). In whites, we identified 3 functional clusters, inflammation, chemokines, and signaling: G protein coupled receptor (GPCR) and receptor tyrosine kinase (RTK) signaling (Figure 5B). The analysis yielded different, frailty-associated protein-protein interaction networks for blacks and whites.
Figure 5. Interaction network analysis of frailty-associated differential gene expression. The interaction network for significant, differentially expressed protein coding genes with frailty in blacks (A) and whites (B) is shown. Network nodes represent proteins and lines (edges) represent protein-protein interactions. The solid lines represent direct interactions between proteins, dashed lines represent indirect interactions between proteins, grey lines represent putative protein interactions. The line colors represent the molecular action type: green (activation), dark blue (binding), black (reaction), red (inhibition), purple (catalysis), yellow (transcriptional regulation). The action effect is represented by the shape at the end of the line: arrow head (positive), perpendicular line (negative), dot (unspecified). The prominent functional clusters in blacks and whites are indicated.
Gene expression changes associated with frailty and race
We designed a validation cohort to increase our sample size and to confirm our results using RT-qPCR. In our validation cohort (Table 1B), we increased our sample size to 52, including 14 participants from the RNA sequencing cohort. We selected genes based on the following criteria: pathways of interest, a significant log2 fold change of + 3, and the presence within functional clusters from network analysis (Figure 5). In addition, we chose genes from the Vitality 90+ study, which is a longitudinal study of successful aging in nonagenarians in Finland . Previously, gene expression was examined in this cohort in the context of survival and mortality. Among the individuals examined were frail individuals. We reanalyzed this data set to identify frailty- associated genes and found 3 genes associated with frailty (FAM116B, SLC2A6, and C17orf56). Therefore, we included these genes in our data set.
Gene specific primers were designed and are listed in Supplementary Table 6. Of the 36 genes (TNF, SLC2A6, RARA, OSM, LRCH3, IL15RA, IL6, IL1B, HBB, FCGR3B, FAM116B, CXCR1, CXCL1, CALCA, C17orf56, PF4, CXCL3, FPR3, CXCL2, PTGS2, ICAM1, GBP1, IRF4, TRIM25, OAS3, EGR1, IFIT1, IFIT3, RSAD2, IFI44, IFI44L, NOTCH1, CDKN1A, CCL5, MAPKAPK2, FPR2) we observed statistically significant changes in the expression of 8 genes. Frailty was associated with significantly lower levels of solute carrier family 2, facilitated glucose transporter member 6 (SLC2A6), Fc fragment of IgG receptor IIIb (FCGR3B), and TEPSIN, adaptor related protein complex 4 accessory protein (C17orf56) (Figure 6A). Significantly higher levels of oncostatin M (OSM), C-X-C motif chemokine ligand 1 (CXCL1), and interleukin 6 (IL6) (Figure 6B) were identified in whites when compared to blacks, independent of frailty status which was not significantly associated with these genes. We also identified a significant two-way interaction between frailty and race for the expression of interleukin 1 beta (IL1B) and early growth response 1 (EGR1) (Figure 7). Sex was not a significant factor in any of the models.
Figure 6. Frailty and race-associated changes in gene expression in the validation cohort. Total RNA was isolated from PBMCs from non-frail and frail blacks and whites in the validation cohort (Table 1B; n=52). Gene expression was analyzed using RT-qPCR with gene specific primers (Refer to Supplementary Table 6). The scatter plots show the relative expression (log2 transformed) in non-frail vs frail (A) and blacks (B) vs whites (W) in this same cohort (B). The open bars represent the mean and error bars show standard error of the mean. Significance was determined using linear regression models on the log2 transformed values.
Figure 7. Frailty-associated changes in gene expression with race. Total RNA was isolated from PBMCs from non-frail and frail blacks and whites in the validation cohort (n=52). Gene expression was analyzed using RT-qPCR with gene specific primers. The scatter plots show the relative expression (log2 transformed) in non-frail and frail blacks (B) and whites (W). The open bars represent the mean and error bars show standard error of the mean. There is a significant two-way interaction between frailty status and race for IL1B (p=0.041) and EGR1 (p=0.019). Significance was determined using linear regression models on the log2 transformed values.
Correlation of IL6 and IL1B gene expression with serum levels
To examine the relationship between IL6 and IL1B gene expression and circulating cytokine levels, we measured serum IL6 and IL1β by ELISA. We did not identify a significant correlation between circulating cytokine levels and gene expression of IL6 (r = -0.05; p = 0.76) or IL1B (r = -0.25; p = 0.08).
In this present study, we have examined the transcriptome of non-frail versus frail middle-aged, community dwelling blacks and whites. Using RNA-seq, we identified over 5,000 differentially expressed genes associated with frailty. There were a larger number of genes differentially expressed in the context of frailty among whites when compared to blacks with a small number of genes that overlapped between the groups. However, there were more genes whose expression was increased among frail blacks when compared to frail whites. In contrast, there were more genes whose expression was decreased among frail whites when compared to frail blacks. Pathway analysis of these differentially expressed genes identified enriched gene ontology sets associated with biological processes; in blacks the most upregulated frailty associated biologic pathways were those associated with chemotaxis, inflammatory and immune response. In contrast, those same pathways were the most down-regulated among frail whites. While there were 19 overlapping pathways between blacks and whites, the directional pattern of these pathways was different. Using RT-qPCR, we validated expression patterns in genes that had a significant log2 fold change in expression, genes that were in functional clusters from our network analysis, genes that were in pathways of interest, or genes that we identified as frail associated and differentially expressed in the literature. Among the 36 genes studied, we identified 8 genes whose expression was influenced by frailty, race or the interaction of frailty and race. These included: solute carrier family 2, facilitated glucose transporter member 6 (SLC2A6), Fc fragment of IgG receptor IIIb (FCGR3B), and TEPSIN, adaptor related protein complex 4 accessory protein (C17orf56) oncostatin M (OSM), C-X-C motif chemokine ligand 1 (CXCL1), interleukin 6 (IL6), interleukin 1 beta (IL1B) and early growth response 1 (EGR1).
Previous work has been done to examine gene expression in frailty. Given that our cohort is middle-aged and diverse, it is noteworthy that we have identified such a large number of genes that are differentially expressed with frailty. This is also important in the context of other studies examining the transcriptome of frailty. Although there are few studies examining the transcriptome of frailty, as part of the Vitality study, gene expression was analyzed in the context of developing a mortality-predicting signature. Frailty status was included as a variable in the model, but gene expression was not analyzed based on frailty status. We reanalyzed these data to identify frailty-associated genes and included these genes in our validation gene set (SLC2A6, FAM116B, and C17orf56). Not only were these genes significantly associated with frailty in the European ancestry, elderly, frail in the Vitality study, they were also significantly associated with frailty in our cohort of middle-aged blacks and whites, suggesting that these genes may have important roles in frailty regardless of age or race. El Assar et al.  examined gene expression in 350 individuals of European ancestry, non-frail and frail, who were 65 years old and older within the Toledo Study of Healthy Aging, a part of the Frailomic Initiative . In this study, the expression of 21 genes associated with aging or the response to oxidative stress were examined in the context of frailty. We also examined expression of these genes in our cohort and identified significant differential expression within our RNA sequencing cohort of four genes of the 21 genes reported by the Toledo Study investigators, TXNRD1, HMOX2, PTGS2, and EGLN3. However, we did not further examine gene expression of TXNRD1, HMOX2, and EGLN3 in the validation cohort because the log2 fold change was less than 1, thus not meeting our selection criteria. However, PTGS2 was one of the putative frailty-associated genes selected in the Toledo Study and met our criteria for validation. However, there were no differences in PTGS2 levels across frailty and race groups in our cohort. The failure to replicate the Toledo findings may be related to differences in the characteristics between the 2 cohorts. The Toledo Study Cohort is substantially older with a mean age of 76.5 years for the overall cohort and 83.0 years for frail participants. The Toledo cohort is a non-diverse European ancestry group and importantly only 22% of the 350 participants are frail. The HANDLS cohort is approximately 36 years younger, consists of blacks and whites, and importantly is balanced for percentage of frail individuals. Another group examined gene expression in robust, pre-frail and frail individuals in the Singapore Longitudinal Study of Aging in response to influenza vaccination . Microarray of peripheral blood mononuclear cells prior to vaccination showed no gene expression differences between the robust, pre-frail or frail participants. Influenza vaccination altered gene expression, however, there were no significant differences between groups that exceeded two-fold differential expression.
Race influenced the transcriptional profiles in our cohort. In our RNA sequencing analysis and cohort, we observed significant gene expression differences between non-frail and frail blacks and whites. There were a larger number of differentially expressed genes associated with frailty in whites when compared to blacks (Figure 2B). The most robust differences were in genes that were more abundant in frail, white individuals compared to non-frail individuals. Many of these genes are significantly more abundant only in whites and not in blacks. The parametric analysis of gene set enrichment (PAGE) and network analyses in the context of frailty showed that the enriched biological processes in pre-frail/frail, middle-aged blacks are associated with immune response, inflammatory response, and chemotaxis (Figure 3), which is consistent with the idea that frailty is a chronic state of inflammation that adversely impacts several organ systems . Interestingly, the opposite pattern was observed in whites in that biological pathways associated with immune response, inflammatory response, and chemotaxis were diminished in pre-frail/frail middle-aged whites (Figure 4) which is consistent with racial differences associated with inflammation-related processes . Enriched biological processes in pre-frail/frail whites are sensory perception of taste, tRNA processing, and sensory detection of bitter taste (Figure 4). Although decreased smell and taste perception is associated with frailty in an elderly Japanese population , our data suggest that these processes are also impacted in middle-aged, frail whites. Collectively, these data suggest that there are differing molecular signatures for frailty between blacks and whites. This may aid in identifying the underlying differences that may explain the frailty prevalence differences between blacks and whites under the age of 50 .
Our RT-qPCR validation study in an expanded cohort of non-frail and frail individuals focused on the 36 genes that we found were significantly changed with frailty status, race or race by frailty in our cohort, and that we found by reanalysis of data from the Vitality 90+ study . The three genes we validated from the Vitality study were SLC2A6, FAM116B, and C17orf56. SLC2A6 encodes for a recently identified solute carrier family 2 facilitated glucose transporter member 6 and is highly expressed in lymphocytes; however, the physiological substrate has yet to be identified . SLC2A6 is a member of class III glucose transporters and possesses dileucine motif that facilitates its translocation to organelle membranes. This suggests that SLC2A6 could facilitate glucose transport across intracellular organelles and warrants further investigation in the context of frailty. Increased SLC2A6 levels were found in frail nonagenarians , but in our validation cohort, we found lower levels of SLC2A6 in frail, middle-aged individuals. These differences could be attributed to racial homogeneity and advanced age of the Vitality cohort in contrast to ours.
Decreased levels of FCGR3B were associated with frail individuals in our validation cohort. FCGR3B encodes for a membrane bound receptor that is important for calcium mobilization and neutrophil degranulation. Dysregulation of calcium homeostasis is a hallmark of aging and aging-related diseases  while dysregulated neutrophil degranulation is a feature of several inflammatory diseases . Therefore, it is interesting to hypothesize that lower levels of FCGR3B may contribute to frailty through dysregulation of calcium homeostasis and through inflammation, both well-established contributors to the frailty phenotype.
C17orf56 encodes for an accessory protein that is a part of a protein complex that plays an important role in endocytosis and clathrin and non-clathrin vesicle formation [28,29]. Here, we found that C17orf56 levels were lower in frail individuals which is contrary to what was observed with frailty in the Vitality 90+ nonagenarian cohort. A hallmark of aging is loss of proteostasis . Maintaining protein homeostasis is mediated, in part by endocytosis , making it interesting to speculate that dysregulation of proteostsis may occur in younger frail individuals. Future work lies in experimentally testing this idea.
In addition to validating gene expression patterns of frailty-associated genes identified in the Vitality 90+ and Toledo studies, we selected genes (OSM, CXCL1, and IL6) that encode cytokines that were found in biological pathways we identified in our PAGE analysis (Table 2). Oncostatin M (OSM) is a diverse cytokine that has recently been identified as a member of the interleukin 6 (IL6) subfamily and also stimulates IL6 production . The chemokine (C-X-C motif) ligand 1 (CXCL1) has a role in inflammation and neutrophil recruitment . Elevated levels of IL6, a proinflammatory cytokine , have been detected in old, frail cohorts and have been considered a potential frailty biomarker . Interestingly, increased expression of OSM, CXCL1, IL6 were not associated with frailty in our cohort. However, race was the driver for increased expression of all three genes in whites when compared to blacks. These data suggest that previously identified markers of frailty may not be appropriate markers for frailty in middle-aged adults or among blacks cohorts.
We also observed interactions between race, frail status, and expression of interleukin 1 beta (IL1B) and the early growth response 1 (EGR1) genes. We observed that frail whites have lower levels of both IL1B and EGR1 compared to non-frail whites. Interestingly, this observation was only in whites, as expression of IL1B and EGR1 are increased with frailty in blacks. IL1B is an important mediator of the inflammatory response and fever induction . EGR1, however, has been investigated in the central nervous system as it is a transcription factor that is required for differentiation, associated with neuroplasticity , and synaptic exocytosis . Although neither gene has been previously studied in frailty, and the functional role of EGR1 has not been extensively investigated in peripheral blood cells, the results from this study indicate that these genes may be novel biomarkers of frailty in younger cohorts.
This current study builds upon our understanding of frailty at the molecular level as we examined the transcriptome of middle-aged blacks and whites. Interestingly, the results from the RNA sequencing, PAGE analysis, and interaction network analysis indicate that race is a driver for the changes that we see in the transcriptional profile of frail, middle-aged adults. This is important considering that we found previously in HANDLS that frailty prevalence in adults under 50 is greater in whites than blacks .
Our study is limited by the sample size, however, similar sample sizes have previously been reported [39–42]. The strengths of the study are the use of next generation sequencing combined with pathway analysis and RT-qPCR validation in a diverse, community-based cohort.
Frailty is a risk factor for disability and mortality. Previous clinical and epidemiologic work in frailty, ours and others, highlights the need to clinically identify frailty earlier in the lifespan. The current work suggests that it is imperative to study the molecular aspects of frailty prior to old age to enable the identification of early molecular drivers of frailty that may aid in intervention to prevent adverse outcomes.
Study participants and frailty assessment
Participants are from the Healthy Aging in Neighborhoods of Diversity across the Life Span (HANDLS) study of the National Institute on Aging Intramural Research Program (NIA IRP), National Institutes of Health. HANDLS is an epidemiologic, longitudinal study in which age-related health disparities are examined in the context of behavioral, psychosocial, and environmental influences . The HANDLS cohort includes 3720 white and black residents of Baltimore, Maryland between the ages of 30-64 at baseline. HANDLS has been approved by the Institutional Review Board of the National Institute on Environmental Health Sciences and all participants provided written, informed consent. Participants within the RNA sequencing (n = 16) and validation cohorts (n = 52; 14 overlap with the RNA sequencing cohort) were stratified across sex, race, frailty status, and were between 45-49 years old (Table 1 A and B). FRAIL scores for the different domains were the same for blacks and whites (Supplementary Table 7). Additionally, when stratifying based on the types of illnesses within the illness domain, there were no differences between blacks and whites (Table 3).
Table 3. Comorbid illnesses in frail blacks and whites in the validation cohort.
|Diabetes = Yes (%)||7 (53.8)||2 (15.4)||0.10|
|Hypertension = Yes (%)||10 (76.9)||7 (53.8)||0.41|
|Myocardial Infarction = Yes (%)*||0 (100.0)||0 (100.0)||1|
|Congestive Heart Failure = Yes (%)*||0 (100.0)||0 (100.0)||1|
|Angina = Yes (%)*||2 (16.7)||2 (16.7)||1|
|Asthma = Yes (%)*||2 (16.7)||4 (33.3)||0.64|
|Stroke = Yes (%)*||0 (0.0)||0 (0.0)||1|
|Kidney Disease = Yes (%)*||0 (0.0)||2 (16.7)||0.46|
|Cancer = Yes (%)*||0 (0.0)||2 (16.7)||0.46|
|Chronic Lung Disease = Yes (%)*||2 (16.7)||4 (33.3)||0.64|
|Arthritis= Yes (%)*||5 (41.7)||4 (33.3)||1|
|Chronic Diseases (mean (sd))*||2.17 (1.75)||2.25 (1.60)||0.90|
|*Information for these illnesses were missing from one black and one white individual. Fisher’s exact test was used to analyze differences for categorical variables and Student’s t-test was used for the continuous variable.|
Frailty was assessed previously  using the FRAIL scale  with adaptions for the loss of weight domain [10,44]. We assessed loss of weight from responses to item two of the Center for Epidemiologic Studies Depression scale (CES-D) . Weight loss was considered present when participants responded occasionally (3–4 days a week) or mostly (5–7 days a week) to the following question: Over the past week did you not feel like eating or have a poor appetite?. The International Academy on Nutrition and Aging FRAIL scale consists of five domains: fatigue, resistance, ambulation, illnesses, and loss of weight. The frail score ranges from 0 – 5, in which one point can be received for each component. A score of 0 is non-frail, 1-2 is pre-frail, and 3-5 is frail. In the RNA sequencing cohort, 8 participants were non-frail, 2 were pre-frail, and 6 were frail.
Next generation sequencing
Peripheral blood mononuclear cells (PBMCs) were isolated from fasting blood samples and stored at -80oC . Total RNA was isolated using TRIzol® (Life Technologies) according to the manufacturer’s instructions with the addition of a DNase treatment and overnight precipitation. Samples were resuspended in 20 µl RNase-free water and frozen at -80°C until further use. Next generation RNA-seq was performed at the Johns Hopkins University Deep Sequencing and Microarray Core Facility. Libraries were prepared with the TruSeq stranded total RNA with Ribo-Zero Gold kit (RS-122-2301). The libraries were multiplexed into 4 plex per pool and sequenced on Illumina HiSeq 2500 and the associated Illumina RTA software (Version number 1.18.66) was used to process the images and generate base-calls to perform the primary analysis. The FASTQ file was then used for sequence input for alignment to the human genome NCBI build v38 using TopHat version 2.21 (Figure 1B). Transcripts were assembled using Cufflinks version 2.2.1, Cuffmerge version 2.2.1 and differential gene expression was determined using Cuffdiff version 2.2.1. Significance was determined using the q-value (FDR adjusted p-value using Benjamini-Hochberg correction for multiple-testing) less than 0.05 and also log2 fold change > 0.58 or < -0.58. RNA-seq data can be accessed at GEO (Accession number: GSE129534). All genes significantly associated with frailty are listed in Supplementary Table 1, Supplementary Table 2 and Supplementary Table 3 list genes significantly altered by frailty in blacks and whites, respectively.
Quantile normalized FPKM values were used to calculate Z-scores (log transformed). Z-ratios were then obtained and imputed into PAGE analysis to analyze gene ontology (GO) terms (Figure 1B) . Significance was based on a false discovery rate of < 0.30, p < 0.05, and at least three genes for that significant gene set. In order to visualize functional interactions between significant, differentially expressed genes, we used the Search Tool for the Retrieval of Interacting Genes (STRING) [21,48]. The interactions were analyzed according to the highest confidence score (0.900) clustered using k-means clustering and disconnected nodes were removed.
Reverse transcription and real-time, quantitative PCR
Total RNA was isolated from PBMCs as described above. 500 µg of total RNA was reverse transcribed into cDNA using random hexamers and Superscript II Reverse Transcriptase (Invitrogen). Reactions were performed with 2X SYBR Green Master Mix and gene-specific primers (Supplement Table 6). Reverse transcription and quantitative real- time PCR (RT-qPCR) reactions were performed on a 7900HT Fast Real-Time PCR System according to the manufacturer’s protocol. mRNA levels were normalized to the average of GAPDH and ACTB. Gene expression was calculated using the 2−ΔΔCt method .
IL6 and IL1β enzyme-linked immunosorbent assays
IL6 and IL1β levels were measured in serum following the manufacturer’s guidelines (Quantikine ELISA Human IL1B/IL-1F2 and Human IL-6; R&D systems). Optical density was measured using the SpectraMax M2 system (Molecular Devices) along with SoftMax Pro (Version 7). IL6 and IL1β concentrations were calculated according to the respective standards. Log transformed values were used to examine the correlation between serum and PBMC mRNA levels using Pearson’s correlation.
The balanced design across sex, race and frailty status for the validation cohort provided 84% power to detect differences in relative expression that accounted for 15% of the total variance. Linear regression models were used to examine the relationship of mRNA levels with frailty status and race. Models were built using forward and backward selection to assess possible interactions among the variables. Log transformations (base 2) were applied for non-normally distributed variables. The association between frailty domains and race were determined using Fisher’s exact tests. Analyses were performed using R . A p-value less than 0.05 was considered to be statistically significant.
The authors would like to thank the HANDLS medical staff, especially physical performance tester, Mrs. Claire Jefferson, and the HANDLS participants. We would also like to thank Dr. Salman Tajuddin for his expertise with R Studio.
Conflicts of Interest
The authors declare that they do not have conflicts of interest.
This work is supported by the National Institute on Aging Intramural Research Program, NIH, Project AG 000513. CSP was also funded through a postdoctoral fellowship through the National Institute on Minority Health and Health Disparities (NIMHD) Intramural Research Program.
- 1. Chen X, Mao G, Leng SX. Frailty syndrome: an overview. Clin Interv Aging. 2014; 9:433–41. https://doi.org/10.2147/CIA.S45300 [PubMed]
- 2. Morley JE. Frailty and sarcopenia in elderly. Wien Klin Wochenschr. 2016 (Suppl 7); 128:439–45. https://doi.org/10.1007/s00508-016-1087-5 [PubMed]
- 3. Fried LP, Tangen CM, Walston J, Newman AB, Hirsch C, Gottdiener J, Seeman T, Tracy R, Kop WJ, Burke G, McBurnie MA, and Cardiovascular Health Study Collaborative Research Group. Frailty in older adults: evidence for a phenotype. J Gerontol A Biol Sci Med Sci. 2001; 56:M146–56. https://doi.org/10.1093/gerona/56.3.M146 [PubMed]
- 4. Xue QL, Walston JD, Fried LP, Beamer BA. Prediction of risk of falling, physical disability, and frailty by rate of decline in grip strength: the women’s health and aging study. Arch Intern Med. 2011; 171:1119–21. https://doi.org/10.1001/archinternmed.2011.252 [PubMed]
- 5. Ensrud KE, Ewing SK, Taylor BC, Fink HA, Stone KL, Cauley JA, Tracy JK, Hochberg MC, Rodondi N, Cawthon PM, and Study of Osteoporotic Fractures Research Group. Frailty and risk of falls, fracture, and mortality in older women: the study of osteoporotic fractures. J Gerontol A Biol Sci Med Sci. 2007; 62:744–51. https://doi.org/10.1093/gerona/62.7.744 [PubMed]
- 6. Makary MA, Segev DL, Pronovost PJ, Syin D, Bandeen-Roche K, Patel P, Takenaga R, Devgan L, Holzmueller CG, Tian J, Fried LP. Frailty as a predictor of surgical outcomes in older patients. J Am Coll Surg. 2010; 210:901–08. https://doi.org/10.1016/j.jamcollsurg.2010.01.028 [PubMed]
- 7. McAdams-DeMarco MA, Suresh S, Law A, Salter ML, Gimenez LF, Jaar BG, Walston JD, Segev DL. Frailty and falls among adult patients undergoing chronic hemodialysis: a prospective cohort study. BMC Nephrol. 2013; 14:224. https://doi.org/10.1186/1471-2369-14-224 [PubMed]
- 8. Bandeen-Roche K, Seplaki CL, Huang J, Buta B, Kalyani RR, Varadhan R, Xue QL, Walston JD, Kasper JD. Frailty in Older Adults: A Nationally Representative Profile in the United States. J Gerontol A Biol Sci Med Sci. 2015; 70:1427–34. https://doi.org/10.1093/gerona/glv133 [PubMed]
- 9. Morley JE, Malmstrom TK, Miller DK. A simple frailty questionnaire (FRAIL) predicts outcomes in middle aged African Americans. J Nutr Health Aging. 2012; 16:601–08. https://doi.org/10.1007/s12603-012-0084-2 [PubMed]
- 10. Griffin FR, Mode NA, Ejiogu N, Zonderman AB, Evans MK. Frailty in a racially and socioeconomically diverse sample of middle-aged Americans in Baltimore. PLoS One. 2018; 13:e0195637. https://doi.org/10.1371/journal.pone.0195637 [PubMed]
- 11. Malmstrom TK, Miller DK, Herning MM, Morley JE. Low appendicular skeletal muscle mass (ASM) with limited mobility and poor health outcomes in middle-aged African Americans. J Cachexia Sarcopenia Muscle. 2013; 4:179–86. https://doi.org/10.1007/s13539-013-0106-x [PubMed]
- 12. Thorpe RJ, Simonsick E, Zonderman A, Evans MK. Association between Race, Household Income and Grip Strength in Middle- and Older-Aged Adults. Ethn Dis. 2016; 26:493–500. https://doi.org/10.18865/ed.26.4.493 [PubMed]
- 13. Woolf SH, Chapman DA, Buchanich JM, Bobby KJ, Zimmerman EB, Blackburn SM. Changes in midlife death rates across racial and ethnic groups in the United States: systematic analysis of vital statistics. BMJ. 2018; 362:k3096. https://doi.org/10.1136/bmj.k3096 [PubMed]
- 14. Dent E, Kowal P, Hoogendijk EO. Frailty measurement in research and clinical practice: A review. Eur J Intern Med. 2016; 31:3–10. https://doi.org/10.1016/j.ejim.2016.03.007 [PubMed]
- 15. Woo J, Yu R, Wong M, Yeung F, Wong M, Lum C. Frailty Screening in the Community Using the FRAIL Scale. J Am Med Dir Assoc. 2015; 16:412–19. https://doi.org/10.1016/j.jamda.2015.01.087 [PubMed]
- 16. Cardoso AL, Fernandes A, Aguilar-Pimentel JA, de Angelis MH, Guedes JR, Brito MA, Ortolano S, Pani G, Athanasopoulou S, Gonos ES, Schosserer M, Grillari J, Peterson P, et al. Towards frailty biomarkers: candidates from genes and pathways regulated in aging and age-related diseases. Ageing Res Rev. 2018; 47:214–77. https://doi.org/10.1016/j.arr.2018.07.004 [PubMed]
- 17. Ferrucci L, Fabbri E. Inflammageing: chronic inflammation in ageing, cardiovascular disease, and frailty. Nat Rev Cardiol. 2018; 15:505–22. https://doi.org/10.1038/s41569-018-0064-2 [PubMed]
- 18. Jylhävä J, Raitanen J, Marttila S, Hervonen A, Jylhä M, Hurme M. Identification of a prognostic signature for old-age mortality by integrating genome-wide transcriptomic data with the conventional predictors: the Vitality 90+ Study. BMC Med Genomics. 2014; 7:54. https://doi.org/10.1186/1755-8794-7-54 [PubMed]
- 19. Cheadle C, Cho-Chung YS, Becker KG, Vawter MP. Application of z-score transformation to Affymetrix data. Appl Bioinformatics. 2003; 2:209–17. [PubMed]
- 20. Cheadle C, Vawter MP, Freed WJ, Becker KG. Analysis of microarray data using Z score transformation. J Mol Diagn. 2003; 5:73–81. https://doi.org/10.1016/S1525-1578(10)60455-2 [PubMed]
- 21. Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, Jensen LJ, von Mering C. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011; 39:D561–68. https://doi.org/10.1093/nar/gkq973 [PubMed]
- 22. El Assar M, Angulo J, Carnicero JA, Walter S, García-García FJ, López-Hernández E, Sánchez-Puelles JM, Rodríguez-Mañas L. Frailty Is Associated With Lower Expression of Genes Involved in Cellular Response to Stress: Results From the Toledo Study for Healthy Aging. J Am Med Dir Assoc. 2017; 18:734.e1–7. https://doi.org/10.1016/j.jamda.2017.04.019 [PubMed]
- 23. Narang V, Lu Y, Tan C, Camous XF, Nyunt SZ, Carre C, Mok EW, Wong G, Maurer-Stroh S, Abel B, Burdin N, Poidinger M, Tambyah PA, et al. Influenza Vaccine-Induced Antibody Responses Are Not Impaired by Frailty in the Community-Dwelling Elderly With Natural Influenza Exposure. Front Immunol. 2018; 9:2465. https://doi.org/10.3389/fimmu.2018.02465 [PubMed]
- 24. Miller MA, Cappuccio FP. Ethnicity and inflammatory pathways - implications for vascular disease, vascular risk and therapeutic intervention. Curr Med Chem. 2007; 14:1409–25. https://doi.org/10.2174/092986707780831131 [PubMed]
- 25. Somekawa S, Mine T, Ono K, Hayashi N, Obuchi S, Yoshida H, Kawai H, Fujiwara Y, Hirano H, Kojima M, Ihara K, Kim H. Relationship between Sensory Perception and Frailty in a Community-Dwelling Elderly Population. J Nutr Health Aging. 2017; 21:710–14. https://doi.org/10.1007/s12603-016-0836-5 [PubMed]
- 26. Vrhovac I, Breljak D, Sabolic I. Glucose transporters in the mammalian blood cells. Period Biol. 2014; 116:131–38.
- 27. Lacy P. Mechanisms of degranulation in neutrophils. Allergy Asthma Clin Immunol. 2006; 2:98–108. https://doi.org/10.1186/1710-1492-2-3-98 [PubMed]
- 28. Takatsu H, Futatsumori M, Yoshino K, Yoshida Y, Shin HW, Nakayama K. Similar subunit interactions contribute to assembly of clathrin adaptor complexes and COPI complex: analysis using yeast three-hybrid system. Biochem Biophys Res Commun. 2001; 284:1083–89. https://doi.org/10.1006/bbrc.2001.5081 [PubMed]
- 29. Hirst J, Bright NA, Rous B, Robinson MS. Characterization of a fourth adaptor-related protein complex. Mol Biol Cell. 1999; 10:2787–802. https://doi.org/10.1091/mbc.10.8.2787 [PubMed]
- 30. López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013; 153:1194–217. https://doi.org/10.1016/j.cell.2013.05.039 [PubMed]
- 31. Jackson MP, Hewitt EW. Cellular proteostasis: degradation of misfolded proteins by lysosomes. Essays Biochem. 2016; 60:173–80. https://doi.org/10.1042/EBC20160005 [PubMed]
- 32. M. T. A. M. (2003). Oncostatin M, a multifunctional cytokine. Reviews of Physiology, Biochemistry, and Pharmacology: Springer, Berlin, Heidelberg), pp. 39-52.
- 33. Kobayashi Y. The role of chemokines in neutrophil biology. Front Biosci. 2008; 13:2400–07. https://doi.org/10.2741/2853 [PubMed]
- 34. Scheller J, Chalaris A, Schmidt-Arras D, Rose-John S. The pro- and anti-inflammatory properties of the cytokine interleukin-6. Biochim Biophys Acta. 2011; 1813:878–88. https://doi.org/10.1016/j.bbamcr.2011.01.034 [PubMed]
- 35. Li H, Manwani B, Leng SX. Frailty, inflammation, and immunity. Aging Dis. 2011; 2:466–73. [PubMed]
- 36. Lopez-Castejon G, Brough D. Understanding the mechanism of IL-1β secretion. Cytokine Growth Factor Rev. 2011; 22:189–95. https://doi.org/10.1016/j.cytogfr.2011.10.001 [PubMed]
- 37. Knapska E, Kaczmarek L. A gene for neuronal plasticity in the mammalian brain: Zif268/Egr-1/NGFI-A/Krox-24/TIS8/ZENK? Prog Neurobiol. 2004; 74:183–211. https://doi.org/10.1016/j.pneurobio.2004.05.007 [PubMed]
- 38. Petersohn D, Thiel G. Role of zinc-finger proteins Sp1 and zif268/egr-1 in transcriptional regulation of the human synaptobrevin II gene. Eur J Biochem. 1996; 239:827–34. https://doi.org/10.1111/j.1432-1033.1996.0827u.x [PubMed]
Dluzen DF, Noren Hooten N, De S, Wood WH
3rd, Zhang Y, Becker KG, Zonderman AB, Tanaka T, Ferrucci L, Evans MK. Extracellular RNA profiles with human age. Aging Cell. 2018; 17:e12785. https://doi.org/10.1111/acel.12785 [PubMed]
- 40. Heyer EE, Deveson IW, Wooi D, Selinger CI, Lyons RJ, Hayes VM, O’Toole SA, Ballinger ML, Gill D, Thomas DM, Mercer TR, Blackburn J. Diagnosis of fusion genes using targeted RNA sequencing. Nat Commun. 2019; 10:1388. https://doi.org/10.1038/s41467-019-09374-9 [PubMed]
- 41. Soreq L, Salomonis N, Guffanti A, Bergman H, Israel Z, Soreq H. Whole transcriptome RNA sequencing data from blood leukocytes derived from Parkinson’s disease patients prior to and following deep brain stimulation treatment. Genom Data. 2014; 3:57–60. https://doi.org/10.1016/j.gdata.2014.11.009 [PubMed]
- 42. Gogliotti RG, Fisher NM, Stansley BJ, Jones CK, Lindsley CW, Conn PJ, Niswender CM. Total RNA Sequencing of Rett Syndrome Autopsy Samples Identifies the M4 Muscarinic Receptor as a Novel Therapeutic Target. J Pharmacol Exp Ther. 2018; 365:291–300. https://doi.org/10.1124/jpet.117.246991 [PubMed]
- 43. Evans MK, Lepkowski JM, Powe NR, LaVeist T, Kuczmarski MF, Zonderman AB. Healthy aging in neighborhoods of diversity across the life span (HANDLS): overcoming barriers to implementing a longitudinal, epidemiologic, urban study of health, race, and socioeconomic status. Ethn Dis. 2010; 20:267–75. [PubMed]
- 44. Theou O, Brothers TD, Mitnitski A, Rockwood K. Operationalization of frailty using eight commonly used scales and comparison of their ability to predict all-cause mortality. J Am Geriatr Soc. 2013; 61:1537–51. https://doi.org/10.1111/jgs.12420 [PubMed]
- 45. Radloff LS. The CES-D Scale:A Self-Report Depression Scale for Research in the General Population. Appl Psychol Meas. 1977; 1:385–401. https://doi.org/10.1177/014662167700100306
- 46. Noren Hooten N, Abdelmohsen K, Gorospe M, Ejiogu N, Zonderman AB, Evans MK. microRNA expression patterns reveal differential expression of target genes with age. PLoS One. 2010; 5:e10724. https://doi.org/10.1371/journal.pone.0010724 [PubMed]
Noren Hooten N, Fitzpatrick M, Wood WH
3rd, De S, Ejiogu N, Zhang Y, Mattison JA, Becker KG, Zonderman AB, Evans MK. Age-related changes in microRNA levels in serum. Aging (Albany NY). 2013; 5:725–40. https://doi.org/10.18632/aging.100603 [PubMed]
- 48. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, Jensen LJ, von Mering C. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017; 45:D362–68. https://doi.org/10.1093/nar/gkw937 [PubMed]
- 49. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25:402–08. https://doi.org/10.1006/meth.2001.1262 [PubMed]
- 50. Team RC. 2018. R: A Language and Environment for Statistical Computing. Vienna, Austria.