Research Paper Volume 11, Issue 4 pp 1189—1203

Age and poverty status alter the coding and noncoding transcriptome

Nicole Noren Hooten1, , Michele K. Evans1, ,

  • 1 Laboratory of Epidemiology and Population Science, National Institute on Aging, National Institutes of Health, Baltimore, MD 21224, USA

Received: November 27, 2018       Accepted: February 8, 2019       Published: February 17, 2019      

https://doi.org/10.18632/aging.101823
How to Cite

Copyright: © 2019 Noren Hooten et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Emerging evidence indicates that noncoding RNAs play regulatory roles in aging and disease. The functional roles of long noncoding RNAs (lncRNAs) in physiology and disease are not completely understood. Little is known about lncRNAs in the context of human aging and socio-environmental conditions. Microarray profiling of lncRNAs and mRNAs from peripheral blood mononuclear cells from young and old white (n=16) and African American (AA) males (n=16) living above or below poverty from the Healthy Aging in Neighborhoods of Diversity across the Life Span study revealed changes in both lncRNAs and mRNAs with age and poverty status in white males, but not in AA males. We validated lncRNA changes in an expanded cohort (n=40); CTD-3247F14.2, GAS5, H19, TERC and MEG3 changed significantly with age, whereas AK022914, GAS5, KB-1047C11.2, MEG3 and XLOC_003262 changed with poverty. Mitochondrial function and response to DNA damage and stress were pathways enriched in younger individuals. Response to stress, viral infection, and immune signals were pathways enriched in individuals living above poverty. These data show that both human age and a marker of social adversity influence lncRNA expression, which may provide insight about molecular pathways underlying aging and social factors that affect disparities in aging and disease.

Introduction

Recent attention has focused on discovering new biomarkers that may serve as indicators of both health span and life span. DNA methylation, telomere erosion, DNA damage and repair, mitochondrial copy number and function are considered markers of biological aging [1]. However, it has been difficult to unequivocally identify a single universal biomarker to monitor aging differentiated from disease because there is considerable variation in the signs and symptoms of aging. Studies have also identified noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), as relevant biomarkers of aging. We and others have shown that miRNA levels change during human aging [26]. Interestingly, these changes in miRNA levels can be modulated by interventions such as calorie restriction or metformin [7,8]. miRNAs are also key regulators of cellular processes important for aging such as cellular senescence [2]. Data from model systems also support a role for miRNAs in the aging process. For example, the Caenorhabditis elegans miRNA lin-4 regulates lifespan through modulating expression of its target gene, lin-14 [9]. Since this initial discovery several other miRNAs have been associated with regulating longevity in C. elegans [10].

Most recently, increasing interest in another class of ncRNAs, long ncRNAs (lncRNAs), has gained momentum as these molecules can also regulate gene expression at the transcriptional, post-transcriptional and translational levels [11,12]. Furthermore, lncRNAs have been linked to processes important for aging and age-related disease [1315]. However, little is known about the global changes in lncRNA expression that occur with aging. Aging is a multi-factorial process marked by a fundamental decline in physiological responses and maintenance of tissue homeostasis and integrity. This decline of cell, tissue and organ function leads to an increase in a myriad of age-related diseases including cancer, type 2 diabetes mellitus, autoimmunity, infections, cardiovascular disease and ultimately mortality. Recent findings from model systems suggest that longevity and health span can be modulated by specific changes in gene expression patterns [10,16]. These patterns could be useful markers for identifying individuals in at-risk populations for the premature development of disease or monitoring adverse outcomes that may lead to early mortality.

Social determinants of health significantly influence the trajectory of aging, health status and outcomes among populations at risk for health disparities. It is well known that poverty and low socioeconomic status remain major risk factors for cardiovascular disease, chronic kidney disease and early mortality [1719]. The weathering hypothesis suggests that socioeconomic disadvantage results in accelerated aging or premature declines in health status among African Americans (AAs) [20]. In the United States, AA men are particularly vulnerable to early mortality [2124]. Therefore, it is important to identify the biological mechanisms and biologic risk factors through which social determinants of health accelerate aging phenotypes and trigger premature mortality. Adverse social conditions, including poverty, poor neighborhood conditions, discrimination, and crime are robust risk factors for negative health outcomes. Social adversity has been linked to gene expression changes in both children and adults [2527]. Specifically, challenging social-environmental conditions initiate changes in gene expression among specific gene sets. The observed increase in mRNAs encoding inflammatory proteins and the decrease in mRNAs encoding immune-response proteins, termed the Conserved Transcriptional Response to Adversity (CTRA) [2527], has been documented in several human and non-human models (for review [28]). These data hint that differences in gene expression may underlie racial and socioeconomic disparities in health.

Here we profiled both lncRNAs and mRNAs in the context of human aging, race and poverty. These lncRNAs may serve as potential biomarkers of susceptibility to age-related diseases. In addition, lncRNAs may be important promoters of the change in gene expression that results from exposure to various social-environmental conditions.

Results

lncRNA changes with age and poverty

Although lncRNAs have been studied in the context of senescence and various other hallmarks of aging, little is known about whether lncRNA expression is altered with human age. Furthermore, since poverty is well known to influence accelerated development of age-associated chronic disease and premature mortality [29,30], we chose a sub-cohort of young and old AAs and whites from the Healthy Aging in Neighborhoods of Diversity across the Life Span (HANDLS) study that were both below and above poverty (Table 1). Poverty was designated as below poverty if the self-reported household income was below 125% of the 2004 Health and Human Services Poverty Guidelines at baseline recruitment. For our initial assessment of lncRNAs, we used the Arraystar human lncRNA microarray V3.0, which profiles about 30,586 lncRNAs and 26,109 coding transcripts. We selected young (30.6 yrs) and old males (63.3 years) that were either white or AA both below or above poverty (n=8/group Table 1).

Table 1. Demographic information for human cohorts.

Microarray cohorts
White malesYoungOld
n88
Age (mean(SD))30.6 (0.5)63.6 (0.5)
PovStat = Below (%)3 (37.5)4 (50)
AA malesYoungOld
n88
Age (mean(SD))30.6 (0.5)63.5 (0.5)
PovStat = Below (%)4 (50)4 (50)
Validation cohort
White malesYoungOld
n2020
Age (mean(SD))32.1 (1.7)62.6 (1.2)
PovStat = Below (%)9 (45)9 (45)

We found that a substantial number of lncRNAs were significantly altered with human age in white males (1938 lncRNAs; Figure 1A; Supplementary Table 1, Supplementary Table 2). In addition, poverty status also significantly changed the levels of lncRNAs in white males (1807 lncRNAs; Figure 1B; Supplementary Table 3, Supplementary Table 4). The top changes in lncRNA abundance with age and poverty are listed in Tables 2 and 3, respectively and shown in the heat maps in Figure 1C. A list of all significantly altered lncRNAs with age and poverty are found in Supplementary Table 1, Supplementary Table 2, Supplementary Table 3, Supplementary Table 4. Intriguingly, lncRNA changes with poverty were not observed in young individuals, but were robust among old individuals, consistent with the idea that poverty over the life span may greatly influence gene expression (Figure 1D). However, we did not assess childhood poverty or socioeconomic status. Comparison of the lncRNAs changed in abundance between poverty and age indicate that 1,029 lncRNAs that were significantly altered overlap with poverty and age (Figure 1E). We observed a different pattern for lncRNA differences with poverty and age in AA males (Supplementary Figure 1), where far fewer lncRNAs were significantly changed in abundance with age (only 26) and with poverty (245).

Changes in lncRNA abundance with age and poverty. The lncRNA expression levels in white young and old males (A) living above or below poverty (B) were assessed using microarrays (Table 1 for demographic details). Volcano plots show log2 fold change and P value for each lncRNA. Red indicates lncRNAs that were >2-fold change and P0.05. (A,B,D) lncRNA comparisons between groups are indicated. (C) Heat maps indicate the top fold changed lncRNAs with age and poverty. These lncRNAs are also listed in Tables 2 and 3. (E) Venn diagram of the total number of significantly increased or decreased lncRNAs for each comparison. Y, young; O, old, AP, above poverty; BP, below poverty.

Figure 1. Changes in lncRNA abundance with age and poverty. The lncRNA expression levels in white young and old males (A) living above or below poverty (B) were assessed using microarrays (Table 1 for demographic details). Volcano plots show log2 fold change and P value for each lncRNA. Red indicates lncRNAs that were >2-fold change and P<0.05. (A,B,D) lncRNA comparisons between groups are indicated. (C) Heat maps indicate the top fold changed lncRNAs with age and poverty. These lncRNAs are also listed in Tables 2 and 3. (E) Venn diagram of the total number of significantly increased or decreased lncRNAs for each comparison. Y, young; O, old, AP, above poverty; BP, below poverty.

Table 2. Top changes in lncRNA abundance with age.

Increased in YoungDecreased in Young
Gene SymbolFoldP valueSeqnameSourceGene SymbolFoldP valueSeqnameSource
CTD-3247F14.218.500.008ENST00000566457GENCODEKB-1047C11.237.340.000ENST00000517655GENCODE
AK02714510.520.001uc002ywy.3UCSC_knowngeneXLOC_00273031.530.016TCONS_00006917LincRNAs; Cabili et al
RP5-826L7.110.140.002ENST00000446476GENCODERP11-573D15.927.850.000ENST00000577781GENCODE
SDHAP19.570.000ENST00000440850pseudogeneRP11-479J7.227.240.002ENST00000425271GENCODE
RP11-773D16.18.940.000ENST00000496359GENCODEXLOC_00837526.510.002TCONS_00018115LincRNAs; Cabili et al
DIP2A-IT18.670.002NR_046400RefSeqLOC40045622.550.001NR_034095RefSeq
CTD-3074O7.58.470.001ENST00000533502GENCODEXLOC_00805822.390.001TCONS_00017457LincRNAs; Cabili et al
RP11-773D16.17.600.003ENST00000488805GENCODEKCNIP4-IT121.320.002NR_002813RefSeq
RP4-803J11.27.460.005ENST00000418348GENCODEAC018730.421.060.001ENST00000454183GENCODE
LOC1005072177.390.011NR_037600RefSeqXLOC_00406920.960.001TCONS_00008575LincRNAs; Cabili et al
BC0798327.160.000uc003nhj.3UCSC_knowngenechr10:7584675-759482518.280.002chr10:7584675-7594825-LincRNAs; Cabili et al
LINC000857.100.003ENST00000573896GENCODEXLOC_00273018.030.008TCONS_00006916LincRNAs; Cabili et al
RP11-510J16.36.890.013ENST00000564138GENCODEXLOC_01173017.920.000TCONS_00024847LincRNAs; Cabili et al
RP13-582O9.56.830.001ENST00000521207GENCODELOC10013211117.580.002NR_024237RefSeq
RP11-1007O24.36.640.005ENST00000565181GENCODEXLOC_00593917.490.002TCONS_00011248LincRNAs; Cabili et al
XLOC_0032626.590.023TCONS_00006666LincRNAs; Cabili et alMBL1P17.350.002ENST00000480805GENCODE
RP11-473M20.146.550.000ENST00000575139GENCODEHMlincRNA15216.720.002HMlincRNA152-LincRNAs; Cabili et al
XLOC_0083826.530.000TCONS_00018121LincRNAs;Cabili et alXLOC_00078116.360.000TCONS_00001451LincRNAs; Cabili et al
WBSCR226.520.008NR_037776RefSeqRP11-281A20.116.340.002ENST00000449730GENCODE
FJ0090346.410.001uc010wyg.2UCSC_knowngeneRP11-525K10.316.310.001ENST00000568776GENCODE
Top lncRNAs differentially expressed with age and P<0.05, Fold change >2 and FDR <0.1.

Table 3. Top changes in lncRNA abundance with poverty.

Increased in Above PovertyDecreased in Above Poverty
Gene SymbolFoldP valueSeqnameSourceGene SymbolFoldP valueSeqnameSource
D6378515.930.006uc002hlt.1UCSC_knowngeneAK0229149.820.003uc001vvr.1UCSC_knowngene
CTD-3247F14.212.740.028ENST00000566457GENCODELOC4004569.410.039NR_034095RefSeq
XLOC_0032629.130.006TCONS_00006666LincRNAs; Cabili et alchr10:7584675-75948258.880.036chr10:7584675-7594825-LincRNAs; Cabili et al
RP11-143M1.29.020.022ENST00000442069GENCODELA16c-390E6.58.810.006ENST00000566287GENCODE
RP1-283E3.48.860.006ENST00000577672GENCODEAC108142.18.390.048ENST00000507869GENCODE
AC092620.38.230.004ENST00000453636GENCODEBC0079408.220.038uc002xzq.1UCSC_knowngene
RP11-510J16.38.150.006ENST00000564138GENCODELINC004268.170.001ENST00000447147GENCODE
RP11-622K12.18.050.003ENST00000423617GENCODERP11-19P22.58.140.017ENST00000583179GENCODE
RPL187.910.003NR_073022RefSeqAC093391.28.010.016ENST00000446492GENCODE
CTD-2576F9.27.550.003ENST00000568835GENCODEXLOC_0077038.010.032TCONS_00016773LincRNAs; Cabili et al
CTC-338M12.47.040.006ENST00000506340GENCODEOSTCP17.700.005ENST00000522287GENCODE
LOC5506436.540.005NR_015367RefSeqSMCR27.560.025ENST00000456090GENCODE
RP5-826L7.16.390.024ENST00000446476GENCODERP11-10N23.27.380.008ENST00000520562GENCODE
RNASEK6.350.001NR_037715RefSeqAC004562.17.270.024ENST00000420427GENCODE
RAC1P26.260.003ENST00000511164pseudogeneRP11-175K6.17.260.027ENST00000523301GENCODE
WBSCR226.010.013NR_037776RefSeqAC098973.27.210.019ENST00000414382GENCODE
LOC1005072176.000.028NR_037600RefSeqAK1300767.170.039uc001kfc.1UCSC_knowngene
RNASEK-C17ORF495.810.002NR_037717RefSeqRP11-378J18.56.960.011ENST00000433391GENCODE
RPL6P275.660.007ENST00000583065GENCODERP11-558A11.36.910.040ENST00000568458GENCODE
LSM15.650.002NR_045493RefSeqXLOC_0055906.8980.035TCONS_00012060LincRNAs; Cabili et al
Top lncRNAs differentially expressed with poverty and P<0.05, Fold change >2 and FDR <0.1.

To determine the effects of race on lncRNA expression, we compared lncRNAs that were altered by race. We found that 291 lncRNAs were altered between white and AA males. Furthermore, if we categorize based on race, poverty and age, we identified changes in lncRNA levels that occurred with both races comparing above and below poverty status. There were 12 lncRNAs that were altered with poverty in both races and these are indicated in gray boxes in Figure 2. Only one lncRNA, RNF157-AS1, changed significantly with age in both whites and AAs, and one lncRNA, C6orf3, changed significantly with both poverty and age in AAs (Figure 2A).

Comparison of lncRNA level changes with age, poverty and senescence. (A) Significantly changed lncRNA abundance was compared between African American (AA) males living below poverty (BP) and above poverty (AB) to white males (W) living below and above poverty. (B) Changes in lncRNAs levels in whites with poverty and age were also compared to significantly changed senescence-associated lncRNAs (SAL-RNAs) identified previously [31]. The comparisons in (A) were done using the Seqname and (B) using Gene symbol, hence the slight differences in numbers. Arrows indicate the direction of change (higher or lower) in the comparisons. The color of the arrow indicates the comparison examined.

Figure 2. Comparison of lncRNA level changes with age, poverty and senescence. (A) Significantly changed lncRNA abundance was compared between African American (AA) males living below poverty (BP) and above poverty (AB) to white males (W) living below and above poverty. (B) Changes in lncRNAs levels in whites with poverty and age were also compared to significantly changed senescence-associated lncRNAs (SAL-RNAs) identified previously [31]. The comparisons in (A) were done using the Seqname and (B) using Gene symbol, hence the slight differences in numbers. Arrows indicate the direction of change (higher or lower) in the comparisons. The color of the arrow indicates the comparison examined.

Comparison to age-associated lncRNAs

Aging and poverty associated changes in lncRNA levels in whites were also compared to lncRNAs that have previously been described in the literature to be age-associated. Since senescent cells accumulate with age, and have been linked to a number of age-related pathologies, we first compared our dataset to lncRNAs that previously were found to be significantly changed in abundance during cellular senescence, termed senescence-associated lncRNAs (SAL-RNAs) [31]. Several poverty-associated lncRNAs (14) and aging-related lncRNAs (13) overlapped with the various SAL-RNAs (Figure 2B). There were 10 SAL-RNAs that overlapped with both comparisons (Figure 2B).

We also further compared the lncRNAs we found to be significantly different to previously reported age-associated lncRNAs. These lncRNAs were compiled based on the regulation of various hallmarks of aging [13] or age-associated functions, including differences with aging in the achilles tendon [32], specific lncRNAs associated with regulating senescence [33,34], mouse premature aging and muscle wasting [35]. This list consisted of 118 “Age-associated lncRNAs” (Supplementary Table 5). Three lncRNAs were overlapping between these Age-associated lncRNAs and our poverty and aging lncRNAs and three were overlapping just with poverty and four just with age (Figure 3A).

Validation of age and poverty associated level changes in lncRNAs abundance. (A) lncRNAs significantly changed in abundance with age and poverty in white males were compared to lncRNAs identified in the literature to be age-associated. Overlapping lncRNAs are shown. (B-C) RNA was isolated from PBMCs from young (mean=32.1 yrs) and old (mean=62.6 yrs) white males above and below poverty (n=40; see Table 1 for demographics). RNA was reverse transcribed and lncRNA levels were quantified by RT-qPCR.

Figure 3. Validation of age and poverty associated level changes in lncRNAs abundance. (A) lncRNAs significantly changed in abundance with age and poverty in white males were compared to lncRNAs identified in the literature to be age-associated. Overlapping lncRNAs are shown. (B-C) RNA was isolated from PBMCs from young (mean=32.1 yrs) and old (mean=62.6 yrs) white males above and below poverty (n=40; see Table 1 for demographics). RNA was reverse transcribed and lncRNA levels were quantified by RT-qPCR.

Validation of lncRNA changes with age and poverty

Several of these overlapping lncRNAs were chosen to further validate in an expanded cohort of individuals (n=40; Table 1). In addition, we also validated the top lncRNAs that were significantly changed in abundance with either age or poverty (Tables 2, 3). All significantly changed lncRNAs are listed in Supplementary Table 1, Supplementary Table 2, Supplementary Table 3, Supplementary Table 4. We designed primers for the top lncRNAs showing the most robust changes in levels and for various “Age-associated lncRNAs” including H19, MEG3, TERC and GAS5. We designed multiple primer sets for XLOC-002730 and LOC400456, but we were unable to amplify these lncRNAs. We chose to further validate in the expanded cohort a total of 13 lncRNAs. RNF157-AS1 was included since levels of this lncRNA were changed differentially with age in AAs and whites, but it was not significantly altered in this expanded cohort. Importantly, CTD-3247F14.2 (LNCipedia: lnc-EGR3-1), GAS5, H19, TERC and MEG3 levels were all significantly changed with age (Figure 3B). AK022914 (DUXAP9),GAS5, KB-1047C11.2 (LNCipedia: C8orf37-AS1:1), MEG3 and XLOC_003262 (LNCipedia: lnc-KY-1) levels were significantly different with poverty (Figure 3C). Interestingly, GAS5 binds and acts as a cytoplasmic decoy for the glucocorticoid receptor (GR) [36]. This observation is intriguing given that transcriptional changes in GR signaling pathways are part of the Conserved Transcriptional Response to Adversity (CTRA) in mice and humans [28].

mRNA changes with age and poverty

In addition to examining lncRNA levels, we also wanted to further understand differences in the levels of expressed mRNAs with poverty and age. To do this, we examined mRNA expression levels by microarray and found that, similar to the lncRNAs, there were many significantly different mRNAs with both age and poverty in white males (Figure 4A,B; Supplementary Table 6, Supplementary Table 7, Supplementary Table 8, Supplementary Table 9). The top changes in mRNA abundance with age and poverty are listed in Tables 4 and 5, respectively and shown in the heat maps in Figure 4C. A list of all significantly altered mRNAs with age and poverty are found in Supplementary Table 6, Supplementary Table 7, Supplementary Table 8, Supplementary Table 9. Changes in expressed mRNAs with poverty were also more robust in the older individuals compared to the younger individuals, similar to our results with lncRNAs (Figure 4D). Many of the differentially abundant mRNAs were overlapping between poverty and age (Figure 4E). Interestingly, comparison of the top mRNAs that were more abundant in both young individuals and those living above poverty showed that out of the top 20 hits 6 overlapped including IL1B, EIF1AY, CCL3L3, FAM96A, HLA-DRB1, and ATRAID (Figure 4C; Tables 4, 5). Out of the top 20 mRNAs that were less abundant with age and poverty status ETS1, C20orf152, NANOG, and SLC18A3 were all overlapping (Figure 4C; Tables 4, 5).

Changes in mRNA expression levels with age and poverty. The levels of expressed mRNAs in white young and old males (A) living above or below poverty (B) were assessed using microarrays. Volcano plots show log2 fold change and P value for each mRNA. Red indicates mRNAs that changed were >2-fold, with P0.05. (A,B,D) Comparisons of significantly different mRNAs between groups are indicated. (C) Heat maps indicate the top fold changed mRNAs with age and poverty. These mRNAs are also listed in Tables 4 and 5. (E) Venn diagram of the total number of significantly increased or decreased mRNAs for each comparison. Y, young; O, old, AP, above poverty; BP, below poverty

Figure 4. Changes in mRNA expression levels with age and poverty. The levels of expressed mRNAs in white young and old males (A) living above or below poverty (B) were assessed using microarrays. Volcano plots show log2 fold change and P value for each mRNA. Red indicates mRNAs that changed were >2-fold, with P<0.05. (A,B,D) Comparisons of significantly different mRNAs between groups are indicated. (C) Heat maps indicate the top fold changed mRNAs with age and poverty. These mRNAs are also listed in Tables 4 and 5. (E) Venn diagram of the total number of significantly increased or decreased mRNAs for each comparison. Y, young; O, old, AP, above poverty; BP, below poverty

Table 4. Top changes in mRNA abundance with age.

Increased in YoungDecreased in Young
Gene SymbolFoldP valueSeqnameGene SymbolFoldP valueSeqname
IL1B13.230.027NM_000576CSMD230.480.002ENST00000241312
EIF1AY12.580.027ENST00000361365PLCD325.300.000NM_133373
CCL3L311.020.009NM_001001437KRTAP9-924.110.003ENST00000394008
FAM216A10.800.002NM_013300KRTAP9-419.410.002NM_033191
FGFBP210.320.002NM_031950C20orf7817.730.003NM_001242671
FAM96A9.420.003NM_001014812C4orf2614.590.000NM_178497
HIGD2A9.330.003NM_138820LRR114.360.000NM_203467
TPI18.500.000NM_000365PARP1614.220.001NM_017851
CLEC4D8.460.004NM_080387ZNF44214.170.002NM_030824
SPTLC27.950.001NM_004863CABLES113.310.005NM_138375
HLA-DRB17.930.014NM_002124IFIT311.770.002NM_001549
TBXAS17.890.018NM_001166254ETS111.350.018NM_005238
COMMD37.840.000NM_012071CR111.210.004NM_000651
TXN7.630.000NM_003329C20orf15210.510.002ENST00000349339
ATRAID7.630.020NM_016085NANOG9.890.005NM_024865
XRCC57.450.016NM_021141SLC18A39.550.003NM_003055
AP1M17.410.000ENST00000291439ADCY69.170.004NM_015270
PDXK7.210.005NM_003681HOXC48.890.000ENST00000430889
ERH7.150.000NM_004450KRTAP2-48.600.002NM_033184
GSTO17.030.004NM_001191003DKK18.560.002NM_012242
Top mRNAs differentially expressed with age and P<0.05, Fold change >2 and FDR <0.1.

Table 5. Top changes in mRNA abundance with poverty.

Increased in Above PovertyDecreased in Above Poverty
Gene SymbolFoldP valueSeqnameGene SymbolFoldP valueSeqname
EIF1AY16.030.014ENST00000361365ETS111.590.018NM_005238
IL1B14.940.020NM_000576AMMECR1L8.720.006NM_031445
CST711.060.004ENST00000480798CDC14B8.020.011NM_003671
HLA-DQB110.060.003NM_002123DDR17.640.005NM_001954
ATRAID9.660.008NM_016085NANOG7.340.021NM_024865
TREM19.460.005NM_018643SOX126.430.016NM_006943
FLJ446359.260.005NM_207422MBD3L16.390.004ENST00000305625
HLA-DRB19.050.008NM_002124GTPBP35.660.022NM_133644
SAT18.750.004NM_002970SLC18A35.630.035NM_003055
NFKBIA8.730.005NM_020529ZDHHC85.530.026NM_013373
CCL3L38.670.023NM_001001437C20orf1525.500.043ENST00000349339
CLEC4E8.660.003NM_014358DCN5.470.003NM_133504
CUTC8.640.004NM_015960SYT45.310.029NM_020783
OAZ18.620.004NM_004152PTF1A5.280.041NM_178161
RPS88.440.001NM_001012FAM153B5.270.006ENST00000253490
JTB8.310.004NM_006694MAGEA65.240.001NM_175868
ISCU8.160.002NM_213595CPSF15.230.010NM_013291
NPM18.070.007NM_199185DCTN15.140.015NM_001135041
FAM96A8.010.003NM_032231CT45A15.070.005NM_001017417
SMYD27.970.002NM_020197CPN14.980.005NM_001308
Top mRNAs differentially expressed with poverty and P<0.05, Fold change >2 and FDR <0.1.

Differences in mRNA abundance were also examined in AA males. There were very few mRNAs showing significant changes in abundance with age in AA males (Supplementary Figure 1). More mRNAs were significantly different comparing individuals above and below poverty (Supplementary Figure 1).

Age and poverty pathway analysis

To gain additional information about the biological pathways that may be regulated by the mRNAs changing in abundance with poverty and age, we performed gene ontology (GO) analysis. In younger individuals, we found enrichment of pathways associated with mitochondrial function and response to DNA damage and stress (Figure 5A). In older individuals, pathways related to chemokine production and development were enriched. In the context of poverty, pathways related to development and differentiation were enriched in individuals living below poverty (Figure 5B). We found that the top pathways enriched among those living above poverty were those implicated in the response to stress, immune stimuli, and viral infection (Figure 5B). This finding is not unexpected because mRNAs increased in abundance in response to social adversity, are part of the inflammatory response cascade, and mRNAs decreased in abundance with social adversity are noted to be related to the innate interferon antiviral response as well as synthesis of antibodies [27].

Pathway analysis for age and poverty. Differentially expressed mRNAs with age and poverty were used for gene ontology (GO) analysis using categories derived from Gene Ontology. Top biological pathways are shown for age (A) and poverty (B).

Figure 5. Pathway analysis for age and poverty. Differentially expressed mRNAs with age and poverty were used for gene ontology (GO) analysis using categories derived from Gene Ontology. Top biological pathways are shown for age (A) and poverty (B).

Discussion

Here, we have examined differences in expressed lncRNAs and mRNAs with age and poverty status in both white and AA males. In our initial genome-wide profiling, we found that white males had larger significant differences in expression levels of lncRNAs and mRNAs with poverty and age, as compared to AA males, where these changes were far more modest. This indicates that the influence of age and poverty on gene expression patterns differs between racial groups. It is not clear what drives these differences, but they may correlate with the wide gap in health status and health outcomes between AAs and whites particularly between AA men and white men. AA males living below poverty in the HANDLS cohort, and in other studies, are particularly vulnerable to early mortality compared to their white counterparts and even when compared to AA men who live above the poverty line despite adjustments for various lifestyle cofactors [21,22]. Given this disparity, it is important to examine whether and how gene expression or other genomic factors may act as a biologic transduction pathway for the social determinants of health.

We found altered expression levels in a number of lncRNAs with age and poverty in white males. Although lncRNAs have been studied in the context of aging, in general most studies have focused on a specific lncRNA or on those that relate to various hallmarks of aging. Here, we have performed a comprehensive study of lncRNAs that change in abundance with human age in men whose life expectancy at birth and life span is shorter than women. We also compared these lncRNAs to a growing list of age-associated lncRNAs from the literature. Several of these overlapping lncRNAs were further validated in an expanded cohort of individuals; among these, GAS5, H19, TERC and MEG3 levels were all significantly changed with age (Figure 3B). GAS5 and MEG3 levels were significantly different with poverty. These data further confirm that these are true age-associated lncRNAs, particularly since these lncRNAs were identified as being aging-associated in other tissues or model systems. Several of these lncRNAs regulate aging processes including telomere function (TERC), autophagy (MEG3), protein trafficking (GAS5), and epigenetic change and proliferation (H19) [13]. Investigators have found that many of these lncRNAs participate in regulating age-related disease [15].

In addition to confirming these age-associated lncRNAs, we also identified many novel age- and poverty-associated lncRNAs (Tables 2, 3). The lncRNAs displaying the largest increases and decreases in abundance were further validated in our expanded cohort; CTD-3247F14.2 (LNCipedia: lnc-EGR3-1) levels were significantly changed with age, while AK022914 (DUXAP9),KB-1047C11.2 (LNCipedia: C8orf37-AS1:1), and XLOC_003262 (LNCipedia: lnc-KY-1) levels were significantly changed with poverty. These data shed new light on various lncRNAs that may be important biological transducers of stressors related to poverty and age, as well as to the loss of homeostasis that may result from both. Future work lies in determining whether the expression of these lncRNAs is a consequence or contributor to the aging process. Interestingly, we found that the levels of one lncRNA, RNF157-AS1, were altered differentially with age in AAs and in whites. It will be interesting in the future to further investigate how this lncRNA expression is altered with age and race.

In summary, we have identified lncRNAs that are altered in abundance with both poverty and age in whites and AAs. As our knowledge of the functional roles of lncRNAs expands, it will be important to determine how these ncRNAs contribute to the aging process and to the biological transduction of social adversity and social determinants that result in health disparities.

Materials and Methods

Clinical participants

Participants were chosen from the Healthy Aging in Neighborhoods of Diversity across the Life Span (HANDLS) study of the National Institute on Aging Intramural Research Program, National Institutes of Health (NIH). The study has been approved by the Institutional Review Board of the National Institute of Environmental Health Sciences, NIH. All participants provided written informed consent. HANDLS is a Baltimore, Maryland based longitudinal epidemiologic study focused on examining the interaction of race and socioeconomic status on aging and age-related health disparities [37]. A sub-cohort of young (~30 yrs) and old (~64 yrs) white or African American (AA) males above or below poverty were chosen for the microarray (n=8/group) (Table 1). An expanded cohort of white males above and below poverty were chosen for RT-qPCR validation (n=20/group; Table 1). Below poverty was designated if the self-reported household income was below 125% of the 2004 Health and Human Services Poverty Guidelines at baseline recruitment.

PBMCs

Fasting blood samples are obtained from participants in the morning and collected in 8-ml Vacutainer® heparinized vials (BD, Franklin Lakes, NJ). Peripheral blood mononuclear cells were isolated within 3 hours of phlebotomy as detailed previously [4]. After isolation, PBMCs are aliquoted and stored at -80oC.

RNA isolation and microarray

Total RNA from PBMCs was isolated using TRIzol™ according to manufacturer’s instructions. RNA quality and quantity were analyzed using a NanoDrop ND-1000 spectrophotometer and standard denaturing agarose gel electrophoresis.

lncRNA and mRNA microarrays were performed by Arraystar Inc. The Arraystar Human LncRNA Microarray V3.0 contains lncRNAs identified through public transcriptome databases (Refseq, UCSC knowngenes, Gencode, etc), as well as from the literature. In brief, sample preparation and microarray hybridization were performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology) with minor modifications. mRNA was purified from total RNA after removal of rRNA (mRNA-ONLY™ Eukaryotic mRNA Isolation Kit, Epicentre). Subsequently, each sample was amplified and transcribed into fluorescent cRNA along the entire length of the transcripts without 3’ bias utilizing a random priming method (Arraystar Flash RNA Labeling Kit, Arraystar). The labeled cRNAs were hybridized onto the Human LncRNA Array v3.0 (8 x 60K, Arraystar). After washing and fixing of the slides, the arrays were scanned by the Agilent Scanner G2505C.

Acquired array images were analyzed using the Agilent Feature Extraction software (version 11.0.1.1). Quantile normalization and subsequent data processing were performed using the GeneSpring GX v12.0 software package (Agilent Technologies). After quantile normalization of the raw data, ComBat was used to adjust batch effects. After normalization of the raw data, lncRNAs and mRNAs that at least 8 out of 16 samples have flags in present or marginal were chosen for further data analysis. Differentially expressed lncRNAs and mRNAs with statistical significance with Fold Change >= 2.0, P value <= 0.05 were identified through Volcano Plot filtering between two groups. All significantly changed lncRNAs are listed in Supplementary Table 1, Supplementary Table 2, Supplementary Table 3, Supplementary Table 4 and mRNAs in Supplementary Table 6, Supplementary Table 7, Supplementary Table 8, Supplementary Table 9. Differentially expressed mRNAs were used for gene ontology (GO) analysis using categories derived from Gene Ontology (http://www.geneontology.org). Biological pathways are shown. Fisher’s exact test was used to calculate significance and P value denotes the significance of GO term enrichment in the differentially expressed genes. Microarray data can be accessed at GEO (Accession Number: GSE123500).

RT-qPCR

Cryopreserved PBMCs were thawed quickly and washed with PBS. Total RNA was isolated using TRIzol™ according to manufacturer’s instructions with an inclusion of a DNase treatment step. RNA quality and quantity were analyzed using a NanoDrop 2000c. RNA was reverse transcribed using random hexamers and reverse transcriptase (Invitrogen). RT-qPCR was performed using SYBR green master mix and lncRNA specific primers. LncRNA primer sequences were designed using NCBI reference sequences or sequences from LNCipedia (www.lncipedia.org). Primer sequences and detailed information are listed in Supplementary Table 10. In some cases, the LNCipedia nomenclature differs from the nomenclature used here and these alternative names are also listed in the text and in Supplementary Table 10. NEAT1 and MEG3 primers were described previously [38,39]. GAS5 has 29 different transcripts and primers were designed using sequences from microarray. Similar data was observed for transcript variants GAS5_019, GAS5_24 and GAS5_026 and data only from GAS5_019 is shown. For other lncRNAs with multiple variants, primers were designed against exons that overlapping between the variants. Reactions were run on an Applied Biosystems 7500 Real-Time PCR System using default settings. lncRNA expression was normalized to an average of HPRT and UBC. Previously we have found that these mRNAs are least variable for aging studies in PBMCs [40]. lncRNA levels were examined for Gaussian distribution by measuring kurtosis and skewness and outliers for each lncRNA were excluded from the analysis using Grubb’s test with an alpha of 0.05.

Acknowledgements

We would like to thank Alan B. Zonderman for help picking the HANDLS sub-cohort used in this study, the HANDLS participants and medical staff and Althaf Lohani for technical assistance. We also would like to thank Elin Lehrmann for help with the GEO submission.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This study was supported by the Intramural Research Program of the National Institute on Aging, National Institutes of Health.

References

  • 1. 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]
  • 2. Abdelmohsen K, Gorospe M. Noncoding RNA control of cellular senescence. Wiley Interdiscip Rev RNA. 2015; 6:615–29. https://doi.org/10.1002/wrna.1297 [PubMed]
  • 3. Jung HJ, Suh Y. MicroRNA in Aging: From Discovery to Biology. Curr Genomics. 2012; 13:548–57. https://doi.org/10.2174/138920212803251436 [PubMed]
  • 4. 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]
  • 5. Noren Hooten N, Fitzpatrick M, Wood WH3rd, 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]
  • 6. Smith-Vikos T, Liu Z, Parsons C, Gorospe M, Ferrucci L, Gill TM, Slack FJ. A serum miRNA profile of human longevity: findings from the Baltimore Longitudinal Study of Aging (BLSA). Aging (Albany NY). 2016; 8:2971–87. https://doi.org/10.18632/aging.101106 [PubMed]
  • 7. Noren Hooten N, Martin-Montalvo A, Dluzen DF, Zhang Y, Bernier M, Zonderman AB, Becker KG, Gorospe M, de Cabo R, Evans MK. Metformin-mediated increase in DICER1 regulates microRNA expression and cellular senescence. Aging Cell. 2016; 15:572–81. https://doi.org/10.1111/acel.12469 [PubMed]
  • 8. Smith-Vikos T, de Lencastre A, Inukai S, Shlomchik M, Holtrup B, Slack FJ. MicroRNAs mediate dietary-restriction-induced longevity through PHA-4/FOXA and SKN-1/Nrf transcription factors. Curr Biol. 2014; 24:2238–46. https://doi.org/10.1016/j.cub.2014.08.013 [PubMed]
  • 9. Boehm M, Slack F. A developmental timing microRNA and its target regulate life span in C. elegans. Science. 2005; 310:1954–57. https://doi.org/10.1126/science.1115596 [PubMed]
  • 10. Smith-Vikos T, Slack FJ. MicroRNAs and their roles in aging. J Cell Sci. 2012; 125:7–17. https://doi.org/10.1242/jcs.099200 [PubMed]
  • 11. Kopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell. 2018; 172:393–407. https://doi.org/10.1016/j.cell.2018.01.011 [PubMed]
  • 12. Noh JH, Kim KM, McClusky WG, Abdelmohsen K, Gorospe M. Cytoplasmic functions of long noncoding RNAs. Wiley Interdiscip Rev RNA. 2018; 9:e1471. https://doi.org/10.1002/wrna.1471 [PubMed]
  • 13. Grammatikakis I, Panda AC, Abdelmohsen K, Gorospe M. Long noncoding RNAs(lncRNAs) and the molecular hallmarks of aging. Aging (Albany NY). 2014; 6:992–1009. https://doi.org/10.18632/aging.100710 [PubMed]
  • 14. Kim C, Kang D, Lee EK, Lee JS. Long Noncoding RNAs and RNA-Binding Proteins in Oxidative Stress, Cellular Senescence, and Age-Related Diseases. Oxid Med Cell Longev. 2017; 2017:2062384. https://doi.org/10.1155/2017/2062384 [PubMed]
  • 15. Kim J, Kim KM, Noh JH, Yoon JH, Abdelmohsen K, Gorospe M. Long noncoding RNAs in diseases of aging. Biochim Biophys Acta. 2016; 1859:209–21. https://doi.org/10.1016/j.bbagrm.2015.06.013 [PubMed]
  • 16. Stegeman R, Weake VM. Transcriptional Signatures of Aging. J Mol Biol. 2017; 429:2427–37. https://doi.org/10.1016/j.jmb.2017.06.019 [PubMed]
  • 17. Crews DC, Charles RF, Evans MK, Zonderman AB, Powe NR. Poverty, race, and CKD in a racially and socioeconomically diverse urban population. Am J Kidney Dis. 2010; 55:992–1000. https://doi.org/10.1053/j.ajkd.2009.12.032 [PubMed]
  • 18. Odutayo A, Gill P, Shepherd S, Akingbade A, Hopewell S, Tennankore K, Hunn BH, Emdin CA. Income Disparities in Absolute Cardiovascular Risk and Cardiovascular Risk Factors in the United States, 1999-2014. JAMA Cardiol. 2017; 2:782–90. https://doi.org/10.1001/jamacardio.2017.1658 [PubMed]
  • 19. Stringhini S, Carmeli C, Jokela M, Avendaño M, Muennig P, Guida F, Ricceri F, d’Errico A, Barros H, Bochud M, Chadeau-Hyam M, Clavel-Chapelon F, Costa G, et al, and LIFEPATH consortium. Socioeconomic status and the 25 × 25 risk factors as determinants of premature mortality: a multicohort study and meta-analysis of 1·7 million men and women. Lancet. 2017; 389:1229–37. https://doi.org/10.1016/S0140-6736(16)32380-7 [PubMed]
  • 20. Geronimus AT. The weathering hypothesis and the health of African-American women and infants: evidence and speculations. Ethn Dis. 1992; 2:207–21. [PubMed]
  • 21. Sloan FA, Ayyagari P, Salm M, Grossman D. The longevity gap between Black and White men in the United States at the beginning and end of the 20th century. Am J Public Health. 2010; 100:357–63. https://doi.org/10.2105/AJPH.2008.158188 [PubMed]
  • 22. Zonderman AB, Mode NA, Ejiogu N, Evans MK. Race and poverty status as a risk for overall mortality in community-dwelling middle-aged adults. JAMA Intern Med. 2016; 176:1394–95. https://doi.org/10.1001/jamainternmed.2016.3649 [PubMed]
  • 23. Kochanek K, Murphy S, Xu J, Arias E. Mortality in the United States, 2016. NCHS Data Brief, no 293. Hyattsville, MD: National Center for Health Statistics. 2017.
  • 24. Bilal U, Diez-Roux AV. Troubling Trends in Health Disparities. N Engl J Med. 2018; 378:1557–58. https://doi.org/10.1056/NEJMc1800328 [PubMed]
  • 25. Gaye A, Gibbons GH, Barry C, Quarells R, Davis SK. Influence of socioeconomic status on the whole blood transcriptome in African Americans. PLoS One. 2017; 12:e0187290. https://doi.org/10.1371/journal.pone.0187290 [PubMed]
  • 26. Powell ND, Sloan EK, Bailey MT, Arevalo JM, Miller GE, Chen E, Kobor MS, Reader BF, Sheridan JF, Cole SW. Social stress up-regulates inflammatory gene expression in the leukocyte transcriptome via β-adrenergic induction of myelopoiesis. Proc Natl Acad Sci USA. 2013; 110:16574–79. https://doi.org/10.1073/pnas.1310655110 [PubMed]
  • 27. Cole SW. Human social genomics. PLoS Genet. 2014; 10:e1004601. https://doi.org/10.1371/journal.pgen.1004601 [PubMed]
  • 28. Slavich GM, Cole SW. The Emerging Field of Human Social Genomics. Clin Psychol Sci. 2013; 1:331–48. https://doi.org/10.1177/2167702613478594 [PubMed]
  • 29. Hero JO, Zaslavsky AM, Blendon RJ. The United States Leads Other Nations In Differences By Income In Perceptions Of Health And Health Care. Health Aff (Millwood). 2017; 36:1032–40. https://doi.org/10.1377/hlthaff.2017.0006 [PubMed]
  • 30. Seeman TE, Crimmins E. Social environment effects on health and aging: integrating epidemiologic and demographic approaches and perspectives. Ann N Y Acad Sci. 2001; 954:88–117. https://doi.org/10.1111/j.1749-6632.2001.tb02749.x [PubMed]
  • 31. Abdelmohsen K, Panda A, Kang MJ, Xu J, Selimyan R, Yoon JH, Martindale JL, De S, Wood WH3rd, Becker KG, Gorospe M. Senescence-associated lncRNAs: senescence-associated long noncoding RNAs. Aging Cell. 2013; 12:890–900. https://doi.org/10.1111/acel.12115 [PubMed]
  • 32. Peffers MJ, Fang Y, Cheung K, Wei TK, Clegg PD, Birch HL. Transcriptome analysis of ageing in uninjured human Achilles tendon. Arthritis Res Ther. 2015; 17:33. https://doi.org/10.1186/s13075-015-0544-2 [PubMed]
  • 33. Puvvula PK, Desetty RD, Pineau P, Marchio A, Moon A, Dejean A, Bischof O. Long noncoding RNA PANDA and scaffold-attachment-factor SAFA control senescence entry and exit. Nat Commun. 2014; 5:5323. https://doi.org/10.1038/ncomms6323 [PubMed]
  • 34. Bianchessi V, Badi I, Bertolotti M, Nigro P, D’Alessandra Y, Capogrossi MC, Zanobini M, Pompilio G, Raucci A, Lauri A. The mitochondrial lncRNA ASncmtRNA-2 is induced in aging and replicative senescence in Endothelial Cells. J Mol Cell Cardiol. 2015; 81:62–70. https://doi.org/10.1016/j.yjmcc.2015.01.012 [PubMed]
  • 35. Lai KM, Gong G, Atanasio A, Rojas J, Quispe J, Posca J, White D, Huang M, Fedorova D, Grant C, Miloscio L, Droguett G, Poueymirou WT, et al. Diverse Phenotypes and Specific Transcription Patterns in Twenty Mouse Lines with Ablated LincRNAs. PLoS One. 2015; 10:e0125522. https://doi.org/10.1371/journal.pone.0125522 [PubMed]
  • 36. Mourtada-Maarabouni M, Pickard MR, Hedge VL, Farzaneh F, Williams GT. GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer. Oncogene. 2009; 28:195–208. https://doi.org/10.1038/onc.2008.373 [PubMed]
  • 37. 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]
  • 38. Mondal T, Subhash S, Vaid R, Enroth S, Uday S, Reinius B, Mitra S, Mohammed A, James AR, Hoberg E, Moustakas A, Gyllensten U, Jones SJ, et al. MEG3 long noncoding RNA regulates the TGF-β pathway genes through formation of RNA-DNA triplex structures. Nat Commun. 2015; 6:7743. https://doi.org/10.1038/ncomms8743 [PubMed]
  • 39. Yoon JH, De S, Srikantan S, Abdelmohsen K, Grammatikakis I, Kim J, Kim KM, Noh JH, White EJ, Martindale JL, Yang X, Kang MJ, Wood WH3rd, et al. PAR-CLIP analysis uncovers AUF1 impact on target RNA fate and genome integrity. Nat Commun. 2014; 5:5248. https://doi.org/10.1038/ncomms6248 [PubMed]
  • 40. Cabili MN, Trapnell C, Goff L, Koziol M, Tazon-Vega B, Regev A, Rinn JL. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 2011; 25:1915–27. https://doi.org/10.1101/gad.17446611 [PubMed]