- 1Division of Infectious Disease, Emory School of Medicine, Atlanta, GA, United States
- 2Department of Epidemiology, Rollins School of Public Health, Emory University, Atlanta, GA, United States
- 3Division of Infectious Diseases Department of Medicine and Center for Clinical Epidemiology and Biostatistics Perelman School of Medicine University of Pennsylvania, Philadelphia, PA, United States
- 4Boston University Medical School, Boston, MA, United States
- 5Connecticut Veteran Health System, West Haven, CT, United States
- 6Yale University School of Medicine, New Haven, CT, United States
- 7Cardiovascular Medicine Division and Tennessee Valley Healthcare System, Vanderbilt University Medical Center, Nashville, TN, United States
- 8Centers for Disease Control and Prevention, Atlanta, GA, United States
- 9Atlanta Veterans Affairs Health Care System, Decatur, GA, United States
- 10Hubert Department of Global Health, Rollins School of Public Health, Atlanta, GA, United States
- 11Emory Vaccine Center, Yerkes National Primate Research Center, Emory University, Atlanta, GA, United States
Background: Liver disease (LD) is an important cause of morbidity and mortality for people with HIV (PWH). The molecular factors linked with LD in PWH are varied and incompletely characterized. We performed an epigenome-wide association study (EWAS) to identify associations between DNA methylation (DNAm) and biomarkers of liver function—aspartate transaminase, alanine transaminase, albumin, total bilirubin, platelet count, FIB-4 score, and APRI score—in male United States veterans with HIV.
Methods: Blood samples and clinical data were obtained from 960 HIV-infected male PWH from the Veterans Aging Cohort Study. DNAm was assessed using the Illumina 450K or the EPIC 850K array in two mutually exclusive subsets. We performed a meta-analysis for each DNAm site measured by either platform. We also examined the associations between four measures of DNAm age acceleration (AA) and liver biomarkers.
Results: Nine DNAm sites were positively associated with serum albumin in the meta-analysis of the EPIC and 450K EWAS after correcting for multiple testing. Four DNAm sites (cg16936953, cg18942579, cg01409343, and cg12054453), annotated within the TMEM49 and four of the remaining five sites (cg18181703, cg03546163, cg20995564, and cg23966214) annotated to SOCS3, FKBP5, ZEB2, and SAMD14 genes, respectively. The DNAm site, cg12992827, was not annotated to any known coding sequence. No significant associations were detected for the other six liver biomarkers. Higher PhenoAA was significantly associated with lower level of serum albumin (β = -0.007, p-value = 8.6 × 10–4, CI: -0.011116, -0.002884).
Conclusion: We identified epigenetic associations of both individual DNAm sites and DNAm AA with liver function through serum albumin in men with HIV. Further replication analyses in independent cohorts are warranted to confirm the epigenetic mechanisms underlying liver function and LD in PWH.
Introduction
The widespread availability and use of effective combination antiretroviral therapy to treat HIV has transformed the outlook for people with HIV (PWH), many of whom now live longer and are facing the challenges of aging with HIV. The burden of non-AIDS related chronic disease (NACD) is increasingly recognized in this patient population as an urgent concern that needs to be addressed to improve care for PWH.
Liver disease (LD) accounts for an estimated 13%–18% all-cause mortality in PWH, and PWH have a high prevalence of traditional causes of LD (Palella et al., 2006; Smith et al., 2010; Smith et al., 2014; Kaspar and Sterling, 2017). Co-infection with hepatitis B and C viruses is estimated between 5%–30% based on cohort studies (Konopnicki et al., 2005; Shepard et al., 2005), and up to 30%–40% of PWH have signs of non-alcoholic fatty liver disease (NAFLD) (Lemoine et al., 2012). Poorly controlled HIV-infection itself is an independent risk factor for liver fibrosis (Kim et al., 2016a), which is the most common outcome of chronic liver injury. Even in individuals with well-controlled HIV (undetectable HIV RNA) and immunologic recovery (CD4+ cell counts > 500 cells/mm3), oxidative stress, mitochondrial injury, toxic metabolite accumulation, gut microbial translocation, systemic inflammation, senescence, and nodular regenerative hyperplasia [reviewed in (Kim et al., 2016a)] contribute to liver injury and fibrosis. Given all of the factors contributing to LD among PWH, there is a need to improve prevention, risk prediction, and treatment of LD in PWH.
Progress in elucidating the biological mechanisms underlying LD has been challenging due to the complexity of the multiple roles played by the liver in homeostasis and metabolic functions. Epigenome-wide association study (EWAS) have the potential to identify novel mechanisms of LD and generate new hypotheses that may provide insights into prevention and treatment of LD. Previous studies on the role of DNA methylation (DNAm) in LD have focused on non-alcoholic fatty LD (NAFLD), since determinants of NAFLD are not as easily identifiable as LD resulting from e.g., alcoholism or viral hepatitis [reviewed in (Zhang et al., 2021)]. Studies on the association between alcoholic liver disease and DNAm have identified genes which may be associated with the incidence and progression of hepatocellular carcinoma (HCC) (French, 2013; Varela-Rey et al., 2013; Zakhari, 2013; Rosen et al., 2018). Other EWAS have uncovered DNAm sites and genes associated with the development of HCC in patients with chronic hepatitis B infections (Su et al., 2007; Nishida et al., 2008; Zhao et al., 2014). PWH are significantly underrepresented in these studies. Moreover, most studies have focused on LD as an overall clinical diagnosis and have not assessed potential associations between DNAm and specific metabolic functions of the liver. Biomarkers are useful surrogates of the liver’s metabolic functions and may reflect its capacity to perform specific functions. The absence of data on the relationship between DNAm and specific markers of liver function, particularly among PWH, represents an important gap in the literature which we hope the present study will begin to fill.
DNAm age acceleration (AA) is a novel biomarker of biological aging, and predicts age-related disease outcomes and mortality (Bell et al., 2019). Several studies have also shown that PWH exhibit elevated DNAm AA above the average level exhibited by non-infected people of equivalent chronological age. (Horvath and Levine, 2015; Rickabaugh et al., 2015). These observations suggest that DNAm AA may partially explain the early onset of chronic comorbid conditions among PWH (Levine et al., 2015a). Relative acceleration of DNAm age from blood-based assays is often measured by four metrics: intrinsic epigenetic age acceleration (IEAA) (Levine et al., 2015a; Horvath and Levine, 2015), extrinsic epigenetic age acceleration (EEAA) (Hannum et al., 2013; Levine et al., 2015b), phenotypic age acceleration (PhenoAA) (Levine et al., 2018), and Grim age acceleration (GrimAA) (Lu et al., 2019).
The present study included available markers of liver function from the Veterans Aging Cohort Study (VACS): aspartate transaminase (AST), alanine transaminase (ALT), serum albumin, total bilirubin, platelet count, fibrosis-4 (FIB-4) score, and the AST-to-platelet-ratio index (APRI) score. We conducted epigenetic association analyses between individual DNAm sites, measures of DNAm AA (IEAA, EEAA, PhenoAA and GrimAA) and these seven biomarkers of liver function in a cohort of male veterans with HIV infection in the United States.
Methods
Study population
A total of 960 male veterans with HIV at the time of blood collection were included from the VACS, having both phenotypic and epigenetic data available from previous studies (Justice et al., 2006; Chen et al., 2020). Participants with cancer at the time of blood collection were excluded because the cancer treatment and disease can influence liver function and DNAm which is beyond the scope of the present study. The VACS consists of electronic medical record data from patients with HIV receiving care at Veterans Health Administration (VA) medical facilities across the United States Each PWH is matched on age, sex, race/ethnicity, and site to two persons without HIV. Data include hospital and outpatient diagnoses (recorded using International Classification of Diseases, Ninth Revision [ICD-9] and Tenth [ICD-10] codes), procedures (recorded using Current Procedural Terminology [CPT] codes), laboratory results, and pharmacy data (Justice et al., 2006). Information on age, race, smoking status, body mass index (BMI), diabetes status (glucose level ≥ 200 mg/dl on 2 separate occasions or glucose level ≥ 200 mg/dl on 1 occasion plus treatment with an oral hypoglycemic or insulin for ≥ 30 days) (Mathur et al., 2019), ever infection with HBV (defined as HBV surface antigen positive, acute resolved HBV or ICD-9 code for HBV diagnosis), ever infection with HCV (defined as HCV antibody positive regardless of HCV RNA, ICD-9 code for HCV diagnosis, or HCV genotype), and antiretroviral (ART) use at time of blood draw were obtained from VA’s clinical data warehouse. Alcohol use status was obtained from the VACS survey (Justice et al., 2006). At the visit of blood draw, AST, ALT, serum albumin, total bilirubin, platelet count, CD4+ T-cell count, and plasma HIV-1 RNA viral load were measured for each patient using approved clinical assays from certified clinical laboratories. Participants were categorized as HIV suppressed if their HIV RNA was < 200 copies/mm3 or as unsuppressed VL otherwise. A fibrosis-4 (FIB-4) score and an AST-to-platelet ratio index (APRI) score were calculated for each participant following the standard equations (Wai et al., 2003; Sterling et al., 2006). All participants provided written consent for the use of their data.
DNAm data generation, processing, and quality control
The genome-wide DNAm profiles were measured using genome DNA samples extract from whole blood. Genomic DNA extraction, epityping, data processing, and quality control were performed as previously described. (Chen et al., 2020). Genome-wide DNAm levels for 473 of the participants were assessed using the Infinium HumanMethylation450 (450K) array platform (Illumina), while those for the remaining 487 participants were assessed using the Infinium HumanMethylationEPIC (EPIC) array platform (Illumina) at the Yale Center for Genomic Analysis. Quality control procedures were followed as previously described in published studies (Zhang et al., 2017; Solomon et al., 2018; Chen et al., 2020). We compared the characteristics between the 450K and EPIC subsets using t-tests for continuous variables and chi-square tests for categorical variables.
DNAm sites for 412,583 and 846,604 DNAm sites were included in the analyses for the dataset acquired with the 450K and EPIC arrays, respectively. Differences in the proportions of the six main leukocyte cell types present in whole blood (CD4+ T cells, CD8+ T cells, monocytes, B cells, granulocytes, and natural killer cells) across samples are well described confounders of associations between DNAm in the blood and many phenotypes. To account for this, the proportions of these six cell types for each participant were determined based on the top cell-type-specific DNAm sites in a reference panel of known proportions following the standard algorithm through the minfi package in R (Houseman et al., 2012; Moran et al., 2016). The estimated cell type proportions were then controlled for in all EWAS analyses.
The association between each covariate included in the EWAS model—race, smoking status, BMI, diabetes status, hazardous alcohol use (Freiberg et al., 1999), ever HCV infection, ever HBV infection, ART use, CD4 count, viral suppression, and leukocyte cell-type proportions, top ten principal components—and each liver marker was assessed using a linear model controlling for chronological age as in the following model.
Individual liver marker ∼ covariate + age
Associations between each liver marker and chronological age were also assessed using a univariate model.
Individual liver marker ∼ age
Regardless of statistical significance, all covariates listed above were included in the final EWAS model to account for the reasonable possibility that they might confound the relationship between DNAm and the liver markers as in the following model.
Individual liver marker ∼ DNAm + age + race +current smoking + BMI + diabetes + hazardous alcohol use + ever infection with HCV + ever infection with HBV + ART use + CD4+ T-cell count + VL + leukocyte cell type proportions + top ten principal components.
Correlations among liver biomarkers
To identify correlations among the biomarkers of liver function included in this study, a Spearman correlation coefficient and corresponding p-value were calculated for every pair of selected biomarkers using the R statistical software package. To normalize the distributions of strongly right-skewed biomarkers for this and all subsequent analyses, a natural log transformation was performed on the measured ALT, AST, and calculated FIB-4 and APRI scores.
Associations of DNAm age acceleration with liver biomarkers
DNAm AA for each participant was measured using the IEAA, EEAA, PhenoAA, and GrimAA metrics as specified by original reporting articles (Chen et al., 2016; Levine et al., 2018; Lu et al., 2019). IEAA is calculated as the residual of a person’s DNAm level in blood cells across a set of DNAm sites proposed by Horvath after regressing DNAm level on chronological age and controlling for proportions of different leukocytes within the sample (Levine et al., 2015a; Horvath and Levine, 2015). EEAA is calculated similarly, except that a different set of DNAm sites proposed by Hanum are used while controlling for proportions of a different set of cell types (Hannum et al., 2013; Levine et al., 2018). PhenoAA, in contrast, is calculated by first estimating phenotypic age (in years) with a linear regression model that uses clinical variables as inputs, then estimating DNAm age (in years) based on methylation levels at a set of DNAm sites proposed by Levine et al. (2018) and finally calculating the residual that results from regressing calculated DNAm age on calculated phenotypic age. GrimAA is calculated in a similar way as PhenoAA but utilizes levels of plasma biomarkers indicative of physiological stress to estimate phenotypic age and a different set of DNAm sites to assess DNAm age (Lu et al., 2019). Each of these measures of DNAm AA perform differently depending on the outcome being evaluated and can complement each other to assess the biological aging process.
The association between DNAm AA, as measured by each of these four metrics, and each of the selected biomarkers of liver function was then assessed by linear regression while controlling for all covariates included in the final EWAS model.
Individual liver marker ∼ DNAm AA + age + race +current smoking + BMI + diabetes + hazardous alcohol use + ever infection with HCV + ever infection with HBV + ART use + CD4+ T-cell count + VL + leukocyte cell type proportions + top ten principal components.
Meta-analysis of EWAS for liver biomarkers
To investigate the association between each biomarker of liver function and individual DNAm sites across the autosomal chromosomes, data from the 450K dataset were analyzed in parallel with data from the EPIC dataset using the same model, and a fixed effect model meta-analysis was conducted using inverse variance weighted effect size method. The final model was adjusted for age, race, current smoking, BMI, diabetes, hazardous alcohol use, ever infection with HCV, ever infection with HBV, ART use, CD4+ T-cell count, VL, leukocyte cell type proportions, and the top ten principal components within the group being analyzed. Where necessary for any given combination of DNAm site and biomarker, participants missing values required for the model were excluded from the analysis for that combination.
A total of 385,062 DNAm sites measured by both the 450K, and EPIC arrays were included in the EWAS meta-analysis. For DNAm sites not covered by both platforms, results were obtained only from the cohort profiled by the covering platform and no meta-analysis was conducted. Separate false-discovery rate adjustments (q < 0.05) were conducted for each liver biomarker. All analyses were performed in R (v4.0.3).
Results
Participant characteristics
After data processing and quality control, the analysis dataset for those profiled using the 450K and EPIC platforms included observations from 473 to 487 individuals, respectively. Characteristics of the two EWAS subsets are summarized in Table 1. All participants were male veterans with HIV who were never diagnosed with cancer and had an average age of 51.2 ± 7.5 years. The chronological age between the 450K and EPIC sub-cohorts was moderately different (Δ = 1.1 years, p-value = 0.031). The distribution of two sub-cohorts did not differ significantly (p-value > 0.05) in race (Black vs. non-Black) or in prevalence of current smoking, diabetes, hazardous alcohol use, ever infection with HBV, ever infection with HCV, or viral load suppression. They also did not differ significantly in their average BMI or CD4+ T-cell count. Average biomarker values from the 450K cohort differed significantly from those from the EPIC cohort for ALT (Δ = 4.7 units/L, p-value = 1.93 × 10–3), AST (Δ = 6.4 units/L, p-value = 7.34 × 10–5), serum albumin (Δ = -0.11 mg/dl, p-value = 2.89 × 10–4), FIB-4 score (Δ = 0.24, p-value = 8.49 × 10–4), and APRI score (Δ = 0.112, p-value = 8.35 × 10–4). Total bilirubin and platelet count did not differ significantly between the two sub-cohorts.
Correlations among liver biomarkers
To illustrate the overall correlation between seven liver function markers, we summarized the pair-wise correlation in Figure 1. Observed values for ALT, AST, FIB-4 score, and APRI score all correlated positively and strongly with each other (Spearman correlation coefficient ρ > 0.6), except for ALT and FIB-4 score with moderate correlation (ρ = 0.36) (Figure 1). Platelet count correlated negatively with AST, ALT, FIB-4 score, and APRI score to varying degrees and weakly with albumin (ρ = 0.13), but not with total bilirubin (ρ = -0.06). Serum albumin correlated negatively with AST, APRI score, and FIB-4 score (−0.25 ≤ ρ ≤ −0.18) but not with ALT (ρ = −0.01). Total bilirubin did not correlate significantly with serum albumin, platelet count, FIB-4 score, or APRI score (p > 0.05) and correlated only weakly with AST and ALT (ρ < 0.14). Distributions of the biomarkers before and after transformation, are presented in Supplementary Figure S1.
FIGURE 1. Correlation matrix based on Spearman coefficients for selected liver markers among all study participants. The strength and direction of each correlation is indicated by the size and color of its representative circle in the matrix. A legend for the colors is included beneath the matrix. Values of the Spearman correlation coefficients are also presented within the circles. Cross-off circles indicate that the observed correlation coefficient is not statistically significant after Bonferroni correction to an overall a = 0.05.
Covariate associations with liver biomarkers
Associations between covariates included in the final EWAS model and the selected liver biomarkers are presented in Table 2. Only “ever infection with HCV” was associated with all seven selected liver biomarkers when controlling for age, while diabetes status was not associated with any of the selected markers. All liver markers associated with at least four of the covariates included in the model excluding cell type proportions. Regardless of the p-value for the observed association between each covariate and each liver biomarker, all covariates were consistently included in the final epigenetic analysis model for each liver biomarker. Participants with complete data of all covariates were analyzed.
TABLE 2. Associations between DNA methylation age acceleration and selected biomarkers of liver function.
Associations of DNAm age acceleration with liver biomarkers
Serum albumin was significantly associated with PhenoAA (beta = −0.014, p-value = 1.6 × 10–13) in an unadjusted model (Figure 2). That association remained significant after adjusting for all covariates in the final regression model (beta = −0.007, p-value = 8.6 × 10–4) (Table 3). No other metric of DNAm AA (i.e., IEAA, EEAA, and GrimAA) was associated with any other selected biomarker of liver function in the adjusted model (Table 3).
FIGURE 2. Scatter plot of serum albumin values across values for phenotypic age acceleration (PhenoAA) among the entire cohort. Red dots indicate participants from 450K and green dots indicate participants from EPIC. An unadjusted best-fit line is included (b = -0.01, p-value = 1.6 × 10–13).
TABLE 3. DNAm sites are significantly associated with selected liver markers after meta-analysis and FDR correction.
EWAS and meta-analysis of liver biomarkers
Meta-analysis of the EWAS results from the EPIC and 450K datasets for each selected liver biomarker and each DNAm site included in both platforms identified that hypermethylation at nine DNAm sites were significantly associated with increased serum albumin among male veterans with HIV after adjusting for all covariates (Table 3; Figure 3A). A regional plot for the section of chromosome 17 that contains five significant DNAm sites is presented in Figure 3B. Quantile-quantile analysis of the expected and observed p-values from the meta-analysis for serum albumin showed very moderate global inflation (inflation factor of 1.02) of the nominal p-values (Figure 3C), so no further corrections were applied. Notably, among the top one hundred DNAm sites most strongly associated with (i.e., with the lowest p-values) serum albumin in the meta-analysis, the beta coefficients obtained from the two EWAS sub-cohorts (the EPIC and 450K) show a strong, positive correlation (beta = 0.89, p-value < 0.001) (Figure 3D).
FIGURE 3. (A) Manhattan plot of unadjusted p-values from meta-analysis of EWAS results for serum albumin. The red line indicates the unadjusted p-value that corresponds to a threshold for FDR significance at Q < 0.05. DNA methylation was not significantly associated with serum albumin at any DNAm site after a more restrictive Bonferroni correction to an overall a = 0.05. (B) Regional Manhattan plot of unadjusted p-values from meta-analysis of EWAS results for CpG sites in chromosome 17 and serum albumin. The red line indicates the unadjusted p-value that corresponds to a threshold for FDR significance (among all CpG sites from across entire genome) at Q < 0.05. DNA methylation was not significantly associated with serum for any CpG site after a more restrictive Bonferroni correction to an overall a = 0.05. (C) Quantile-quantile plot of unadjusted p-values from meta-analysis of EWAS results for serum albumin. The global inflation factor was 1.02. Red lines represent a distribution of null associations with a 95% confidence interval. (D) Correlation of beta coefficients obtained from separate EWAS of serum albumin among the EPIC and 450K cohorts for the 100 DNAm sites with the smallest unadjusted p-values after meta-analysis. Beta coefficients from the separate analyses are positively correlated (r = 0.89, p-value < 0.001), indicating that the signs of observed beta coefficient are consistent between the two datasets. Circles representing p-value pairs that were significant after meta-analysis and FDR correction to Q < 0.05 are highlighted in red.
Four of the nine DNAm sites was positively associated with serum albumin—cg16936953, cg18942579, cg01409343, and cg12054453—annotated to be located within the TMEM49 gene. A 10% increase in methylation at each of those sites corresponded to average increases in serum albumin of 0.12 g/dl (95% CI: 0.08, 0.16 g/dl), 0.14 g/dl (95% CI: 0.09, 0.19 g/dl), 0.17 g/dl (95% CI: 0.11, 0.23 g/dl), and 0.10 g/dl (95% CI: 0.06, 0.14 g/dl), respectively. Four of the remaining five DNAm sites—cg18181703, cg03546163, cg20995564, and cg23966214—annotated to correspond to the SOCS3, FKBP5, ZEB2, and SAMD14 genes, respectively. A 10% increase in methylation at each of those DNAm sites corresponded to average increases in serum albumin of 0.19 g/dl (95% CI: 0.12, 0.26 g/dl), 0.11 g/dl (95% CI: 0.07, 0.15 g/dl), 0.12 g/dl (95% CI: 0.07, 0.17 g/dl), and 0.35 g/dl (95% CI: 0.21, 0.49 g/dl). The remaining DNAm site, cg12992827, was not annotated to correspond to any known coding sequence. A 10% increase in methylation at that site corresponded to an average increase in serum albumin of 0.19 g/dl (95% CI: 0.11, 0.26 g/dl). A quantile-quantile plot and Manhattan plot of the nominal p-values obtained from meta-analysis of the EPIC and 450K sub-cohorts are presented in Supplementary Figures S2, S3 for the other six liver biomarkers. Additionally, statistics and annotations for the ten DNAm sites with the lowest p-values after meta-analysis are presented in Supplementary Table S1 for each tested liver marker. Among the DNAm sites that were included in only one of the two methylation assay platforms used for this study, and that were not included in the meta-analysis, none was significantly associated with any of the liver biomarkers after correcting for multiple testing. Statistics and annotations for the ten DNAm sites with the lowest p-values among those that could not be meta-analyzed are presented for each liver biomarker in Supplementary Table S2.
Discussion
In this study we evaluated the epigenetic associations with seven clinically relevant biomarkers of liver function—aspartate transaminase, alanine transaminase, serum albumin, total bilirubin, platelet count, FIB-4 score, and APRI score—among male United States veterans with HIV. We identified that nine DNAm sites mapped to the genes TMEM49, SOCS3, FKBP5, ZEB2, and SAMD14 were significantly associated with serum albumin levels. We also demonstrate an association between higher PhenoAA (an aging marker for DNAm AA) and low serum levels of albumin.
Serum albumin is a highly abundant protein primarily synthesized by the liver. Synthesis of albumin occurs in polysomes of hepatocytes at a rate of 10 g/day–15 g/day and accounts of 10% of liver protein synthesis (Garcovich et al., 2009). Albumin is not stored in the liver and therefore is not released on demand. Instead in situations of increased need in healthy adults, hepatocytes can increase albumin synthesis by up to 300% (Garcovich et al., 2009). Impaired albumin synthesis and function have been reported in many liver diseases (Carvalho and Machado, 2018), and low serum levels of albumin may be a useful indicator of LD. Several studies have shown that serum albumin is an important prognostic factor and a significant predictor of death in patients with cirrhosis (Carvalho and Machado, 2018) and also a predictor of serious non-AIDS events in PWH (Ronit et al., 2018a; Ronit et al., 2018b).
Four of the nine serum albumin associated DNAm sites are located near the TMEM49 gene (also known as VMP1). DNAm site TMEM49-cg16936953 was also associated with HIV infection, C-reactive protein and inflammatory bowel disease (Gross et al., 2016; Ligthart et al., 2016; Ventham et al., 2016). TMEM49 encodes a transmembrane protein that drives cellular autophagy (Molejon et al., 2013) and is overexpressed in pancreatitis-affected acinar cells (Vaccaro et al., 2003). Cellular autophagy, therefore, may be linked to regulation of serum albumin. Previous studies have shown that serum albumin suppresses autophagy via mTOR activation and that depriving cultured cells of serum albumin induces autophagy which is protective against the accumulation of harmful reactive oxygen species (ROS).
SOCS3 encodes a protein that inhibits cytokine signaling in the STAT pathway in response to elevated cytokine levels (Jo et al., 2005; Carow and Rottenberg, 2014). Suppression of SOCS3 contributes to liver fibrosis by increasing fibrogenic signaling via STAT3-mediated upregulation of tissue growth factor (TGF)-β (Ogata et al., 2006; Jadid et al., 2018; Dees et al., 2020). In addition, obesity has been associated with downregulation of SOCS3, which raises the possibility that an association between DNAm near SOCS3 and liver dysfunction could result from the well described relationship between obesity and liver dysfunction. DNAm site SOCS3-cg18181703 was associated with BMI in several recent EWAS (Ali et al., 2016; Mendelson et al., 2017; Wahl et al., 2017; Xu et al., 2018; Sun et al., 2019). In this cohort of PWH, SOCS3-cg18181703 was associated with BMI (p-value:0.044) further supporting this hypothesis. Epigenetic associations between SOCS3-cg18181703 and C-reactive protein, inflammatory bowel disease, type-2 diabetes and cognitive abilities (Chambers et al., 2015; Ligthart et al., 2016; Marioni et al., 2018; Juvinao-Quintero et al., 2021) also suggest potential pleiotropic effects of SOCS3 gene.
FKBP5 encodes an immunoregulatory protein that contributes to basic protein folding and trafficking (Zannas et al., 2016). FKBP5 has been demonstrated to positively regulate the stress response and drive acquisition of metabolic disorders including obesity, insulin resistance, and diabetes (Sidibeh et al., 2018). Moreover, FKBP5 is thought to contribute to liver dysfunction (Kusumanchi et al., 2021) and in one study, deletion of FKBP5 protected knock-out mice from fatty liver disease despite high fat diets (Stechschulte et al., 2016). ZEB2 encodes a zinc-finger DNA-binding protein expressed in hepatocytes (Cai et al., 2012). Studies have demonstrated that repression of ZEB2 expression by microRNAs induces apoptosis in hepatocytes (Zhao et al., 2018), suggesting that ZEB2 might influence preservation of liver integrity in the face of hepatocellular damage. The remaining DNAm site associated positively with serum albumin levels was located near the SMAD14 gene. The function of the SMAD14 gene is unknown, but it has been linked to gastric cancer and hypothesized to act as a tumor suppressor gene (Xu et al., 2020). No studies have associated SMAD14 or its tumor suppressor properties to either general liver health and function or specifically to serum albumin regulation.
Despite the association between PhenoAA and serum albumin, no other significant associations were noted between the four assessed metrics of DNAm AA—PhenoAA, GrimAA, IEAA, and EEAA—and other biomarkers. Compared to the other metrics of DNAm aging, PhenoAA is unique in that it is highly predictive of morbidity and mortality even after adjusting for chronological age (Liu et al., 2019). This may explain its significant association with albumin, which is also a prognostic marker of survival in LD (Carvalho and Verdelho Machado, 2018). It is plausible that specific DNAm sites may regulate the expression of genes in a manner that influences behavior of the selected biomarkers of liver health, either by impacting the progression of conditions that damage the liver and inhibit function—like fibrosis or hepatocellular carcinoma—or by affecting pathways that are specifically involved in regulation of a specific liver biomarkers.
Low serum level of albumin, though an important prognostic biomarker for LD (Carvalho and Verdelho Machado, 2018), is not specific. Thus, the associations we have uncovered in this study may be mediated by other non-hepatic factors in PWH. For example, hypoalbuminemia is associated with inflammation and synthesis of albumin is suppressed by pro-inflammatory cytokines TNF-α and Interleukin-1 (Dinarello, 1984; Perlmutter et al., 1986), both of which are important mediators of chronic inflammation in PWH (Deeks et al., 2013; De Pablo-Bernal et al., 2014). It is thus presumptuous to infer that the association between DNAm and albumin reflects underlying LD alone and is not combination of multiple interacting factors.
Strengths and limitations
This study combines data from two different microarray platforms to assess consistent epigenetic associations of overlapping DNAm sites. This enabled a robust meta-analysis of epigenetic associations with a large sample size. The significant associations are likely to be true positive findings, and should be replicated in future studies. Secondly, we examine multiple biomarkers which when combined reflect overall liver function and individually inform on the liver’s capacity to perform specific metabolic functions. This allowed us to examine DNAm sites across a wider range of genes that may link to specific aspects of liver function. Thirdly, unlike other EWAS of liver function or LD, our cohort included exclusively people with HIV, a group which has a higher burden of liver dysfunction and related morbidity. Our findings contribute to the identification of genes that may be relevant to liver health and function in the context of HIV infection.
One limitation of our study is that EWAS discovery was sub-optimally powered for a subset of DNAm sites included in only one of the platforms. Meta-analysis results from DNAm sites measured on both platforms were highly consistent, suggesting that combining data from the two platforms successfully augments the power to detect significant associations. Future replication and meta-analysis with extended epigenomic coverage would improve the discovery of liver function-associated DNAm sites. Secondly, our study population was restricted to male veterans with HIV, limiting the generalizability of our findings to women with HIV. Thirdly, we used DNAm levels measured in blood cells to assess associations between DNAm and biomarkers of liver function, which may not represent DNAm patterns in the liver cells that are responsible for regulating these biomarkers. The lack of significant associations between DNAm and AST, ALT, total serum bilirubin, platelet count, FIB-4 score and, APRI score might result from an absence of epigenetic associations with liver functions of large effects, suggesting better powered future studies are needed to identify epigenetic associations with liver functions with smaller effects. Also, the potential for misclassifying liver fibrosis with the use of clinical scores like FIB-4 and APRI as has been shown in hepatitis B infection (Kim et al., 2016b) and may have limited our ability to detect an association with these markers.
Conclusion
We identified epigenetic associations of both individual DNAm sites and DNAm AA with liver function through serum albumin in men with HIV. EWAS may inform on disease pathogenesis and generate new hypotheses for predicting disease progression and treatment of LD in HIV. Identifying specific genes linked with liver function among PWH can inform disease progression to improve the health of this patient population. Further replication analyses in independent cohorts are warranted to confirm and improve the discovery of the epigenetic markers of liver function and LD in PWH.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: Due to US Department of Veterans Affairs (VA) regulations and our ethics agreements, the analytic data sets used for this study are not permitted to leave the VA firewall without a Data Use Agreement. This limitation is consistent with other studies based on VA data. However, VA data are made freely available to researchers with an approved VA study protocol. For more information, please visit https://www.virec.research.va.gov or contact the VA Information Resource Center at VIReC@va.gov. Requests to access these datasets should be directed to yan.v.sun@emory.edu.
Ethics statement
The datasets for this article are not publicly available due to concerns of ethics agreement with the US Department of Veterans Affairs (VA). These data are part of the VACS database. Requests to access these datasets should be directed to the corresponding author, YS, yan.v.sun@emory.edu. The institutional review board approved the original VACS cohort, and all participants provided written informed consent.
Author contributions
BT and ML—Original manuscript draft and data analysis ZW, JC, QH—Data analysis VR, KA, KX, MF, MG, VM, AJ—Manuscript revision YS—Study design, data and, analysis and, manuscript revision.
Funding
This work was supported by the following grants: BT, ML, ZW, JC, QH, VM and YS received support from the National Institute of Health (R01DK125187), VM and YS received support from Emory CFAR (P30AI050409), KX received support from the National Institute on Drug Abuse (R01DA042691, R01DA047820, R01DA047063), KS received support from NIH K01HL134147, VACS - NIH NIAAA U24 AA020794, U01 AA020790, U10 AA013566 completed.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1020871/full#supplementary-material
Abbreviations
AA, Age acceleration; ALT, Alanine amino transferase; APRI, AST to platelet ratio index; ART, Antiretroviral therapy; AST, Aspartate amino transferase; BMI, Body mass index; CpG, 5′—Cytosine—phosphate—Guanine—3′; DNAm, Deoxyribonucleic acid methylation; EEAA, Extrinsic epigenetic age acceleration; EWAS, Epigenome-wide association study; FIB, Fibrosis; FKBP, FK506 binding protein; HBV, Hepatitis B virus; HCC, Hepatocellular carcinoma; HCV, Hepatitis C virus; HIV, Human immunodeficiency virus; IEAA, Intrinsic epigenetic age acceleration; LD, Liver Disease; NACD, Non-AIDS related chronic disease; NAFLD, Non-Alcohol fatty liver disease; PWH, People with HIV; SOCS, suppressor of cytokine signaling; STAT, Signal transducer and activator of transcription; TGF, Tissue growth factor; TMEM, transmembrane; TNFα, tumor necrosis factor- alpha; VACS, Veterans aging cohort; ZEB, Zinc Finger E-Box Binding Homeobox 1.
References
Ali, O., Cerjak, D., Kent, J. W., James, R., Blangero, J., Carless, M. A., et al. (2016). Methylation of SOCS3 is inversely associated with metabolic syndrome in an epigenome-wide association study of obesity. Epigenetics 11, 699–707. doi:10.1080/15592294.2016.1216284
Bell, C. G., Lowe, R., Adams, P. D., Baccarelli, A. A., Beck, S., Bell, J. T., et al. (2019). DNA methylation aging clocks: Challenges and recommendations. Genome Biol. 20, 249. doi:10.1186/s13059-019-1824-y
Cai, M. Y., Luo, R. Z., Chen, J. W., Pei, X. Q., Lu, J. B., Hou, J. H., et al. (2012). Overexpression of ZEB2 in peritumoral liver tissue correlates with favorable survival after curative resection of hepatocellular carcinoma. PLoS One 7, e32838. doi:10.1371/journal.pone.0032838
Carow, B., and Rottenberg, M. E. (2014). SOCS3, a major regulator of infection and inflammation. Front. Immunol. 5, 58. doi:10.3389/fimmu.2014.00058
Carvalho, J. R., and Machado, M. V. (2018). New insights about albumin and liver disease. Ann. Hepatol. 17, 547–560. doi:10.5604/01.3001.0012.0916
Carvalho, J. R., and Verdelho Machado, M. (2018). New insights about albumin and liver disease. Ann. Hepatol. 17, 547–560. doi:10.5604/01.3001.0012.0916
Chambers, J. C., Loh, M., Lehne, B., Drong, A., Kriebel, J., Motta, V., et al. (2015). Epigenome-wide association of DNA methylation markers in peripheral blood from Indian asians and Europeans with incident type 2 diabetes: A nested case-control study. Lancet. Diabetes Endocrinol. 3, 526–534. doi:10.1016/S2213-8587(15)00127-8
Chen, B. H., Marioni, R. E., Colicino, E., Peters, M. J., Ward-Caviness, C. K., Tsai, P. C., et al. (2016). DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY) 8, 1844–1865. doi:10.18632/aging.101020
Chen, J., Huang, Y., Hui, Q., Mathur, R., Gwinn, M., So-Armah, K., et al. (2020). Epigenetic associations with estimated glomerular filtration rate among men with human immunodeficiency virus infection. Clin. Infect. Dis. 70, 667–673. doi:10.1093/cid/ciz240
De Pablo-Bernal, R. S., Ruiz-Mateos, E., Rosado, I., Dominguez-Molina, B., Alvarez-Rios, A. I., Carrillo-Vico, A., et al. (2014). TNF-α levels in HIV-infected patients after long-term suppressive cART persist as high as in elderly, HIV-uninfected subjects. J. Antimicrob. Chemother. 69, 3041–3046. doi:10.1093/jac/dku263
Deeks, S. G., Tracy, R., and Douek, D. C. (2013). Systemic effects of inflammation on health during chronic HIV infection. Immunity 39, 633–645. doi:10.1016/j.immuni.2013.10.001
Dees, C., Pötter, S., Zhang, Y., Bergmann, C., Zhou, X., Luber, M., et al. (2020). TGF-β-induced epigenetic deregulation of SOCS3 facilitates STAT3 signaling to promote fibrosis. J. Clin. Invest. 130, 2347–2363. doi:10.1172/JCI122462
Dinarello, C. A. (1984). Interleukin-1 and the pathogenesis of the acute-phase response. N. Engl. J. Med. 311, 1413–1418. doi:10.1056/NEJM198411293112205
Freiberg, M. S., McGinnis, K. A., Kraemer, K., Samet, J. H., Conigliaro, J., CuRtis Ellison, R., et al. (1999). The association between alcohol consumption and prevalent cardiovascular diseases among HIV-infected and HIV-uninfected men. J. Acquir. Immune Defic. Syndr. 53, 247–253. doi:10.1097/QAI.0b013e3181c6c4b7
French, S. W. (2013). Epigenetic events in liver cancer resulting from alcoholic liver disease. Alcohol Res. 35, 57–67.
Garcovich, M., Zocco, M. A., and Gasbarrini, A. (2009). Clinical use of albumin in hepatology. Blood Transfus. 7, 268–277. doi:10.2450/2008.0080-08
Gross, A. M., Jaeger, P. A., Kreisberg, J. F., Licon, K., Jepsen, K. L., Khosroheidari, M., et al. (2016). Methylome-wide analysis of chronic HIV infection reveals five-year increase in biological age and epigenetic targeting of HLA. Mol. Cell 62, 157–168. doi:10.1016/j.molcel.2016.03.019
Hannum, G., Guinney, J., Zhao, L., Zhang, L., Hughes, G., Sadda, S., et al. (2013). Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol. Cell 49, 359–367. doi:10.1016/j.molcel.2012.10.016
Horvath, S., and Levine, A. J. (2015). HIV-1 infection accelerates age according to the epigenetic clock. J. Infect. Dis. 212, 1563–1573. doi:10.1093/infdis/jiv277
Houseman, E. A., Accomando, W. P., Koestler, D. C., Christensen, B. C., Marsit, C. J., Nelson, H. H., et al. (2012). DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinforma. 13, 86. doi:10.1186/1471-2105-13-86
Jadid, F. Z., Chihab, H., Alj, H. S., Elfihry, R., Zaidane, I., Tazi, S., et al. (2018). Control of progression towards liver fibrosis and hepatocellular carcinoma by SOCS3 polymorphisms in chronic HCV-infected patients. Infect. Genet. Evol. 66, 1–8. doi:10.1016/j.meegid.2018.08.027
Jo, D., Liu, D., Yao, S., Collins, R. D., and Hawiger, J. (2005). Intracellular protein therapy with SOCS3 inhibits inflammation and apoptosis. Nat. Med. 11, 892–898. doi:10.1038/nm1269
Justice, A. C., Dombrowski, E., Conigliaro, J., Fultz, S. L., Gibson, D., Madenwald, T., et al. (2006). Veterans aging cohort study (VACS): Overview and description. Med. Care 44, S13–S24. doi:10.1097/01.mlr.0000223741.02074.66
Juvinao-Quintero, D. L., Marioni, R. E., Ochoa-Rosales, C., Russ, T. C., Deary, I. J., van Meurs, J. B. J., et al. (2021). DNA methylation of blood cells is associated with prevalent type 2 diabetes in a meta-analysis of four European cohorts. Clin. Epigenetics 13, 40. doi:10.1186/s13148-021-01027-3
Kaspar, M. B., and Sterling, R. K. (2017). Mechanisms of liver disease in patients infected with HIV. BMJ Open Gastroenterol. 4, e000166. doi:10.1136/bmjgast-2017-000166
Kim, H. N., Nance, R., Van Rompaey, S., Delaney, J. C., Crane, H. M., Cachay, E. R., et al. (2016). Poorly controlled HIV infection: An independent risk factor for liver fibrosis. J. Acquir. Immune Defic. Syndr. 72, 437–443. doi:10.1097/QAI.0000000000000992
Kim, W. R., Berg, T., Asselah, T., Flisiak, R., Fung, S., Gordon, S. C., et al. (2016). Evaluation of APRI and FIB-4 scoring systems for non-invasive assessment of hepatic fibrosis in chronic Hepatitis B patients. J. Hepatol. 64, 773–780. doi:10.1016/j.jhep.2015.11.012
Konopnicki, D., Mocroft, A., de Wit, S., Antunes, F., Ledergerber, B., Katlama, C., et al. (2005). Hepatitis B and HIV: Prevalence, AIDS progression, response to highly active antiretroviral therapy and increased mortality in the EuroSIDA cohort. AIDS 19, 593–601. doi:10.1097/01.aids.0000163936.99401.fe
Kusumanchi, P., Liang, T., Zhang, T., Ross, R. A., Han, S., Chandler, K., et al. (2021). Stress-responsive gene FK506-binding protein 51 mediates alcohol-induced liver injury through the hippo pathway and chemokine (C-X-C motif) ligand 1 signaling. Hepatology 74, 1234–1250. doi:10.1002/hep.31800
Lemoine, M., Serfaty, L., and Capeau, J. (2012). From nonalcoholic fatty liver to nonalcoholic steatohepatitis and cirrhosis in HIV-infected patients: Diagnosis and management. Curr. Opin. Infect. Dis. 25, 10–16. doi:10.1097/QCO.0b013e32834ef599
Levine, M. E., Hosgood, H. D., Chen, B., Absher, D., Assimes, T., and Horvath, S. (2015). DNA methylation age of blood predicts future onset of lung cancer in the women's health initiative. Aging (Albany NY) 7, 690–700. doi:10.18632/aging.100809
Levine, M. E., Lu, A. T., Bennett, D. A., and Horvath, S. (2015). Epigenetic age of the pre-frontal cortex is associated with neuritic plaques, amyloid load, and Alzheimer's disease related cognitive functioning. Aging (Albany NY) 7, 1198–1211. doi:10.18632/aging.100864
Levine, M. E., Lu, A. T., Quach, A., Chen, B. H., Assimes, T. L., Bandinelli, S., et al. (2018). An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY) 10, 573–591. doi:10.18632/aging.101414
Ligthart, S., Marzi, C., Aslibekyan, S., Mendelson, M. M., Conneely, K. N., Tanaka, T., et al. (2016). DNA methylation signatures of chronic low-grade inflammation are associated with complex diseases. Genome Biol. 17, 255. doi:10.1186/s13059-016-1119-5
Liu, Z., Kuo, P-L., Horvath, S., Crimmins, E., Ferrucci, L., and Levine, M. (2019). A new aging measure captures morbidity and mortality risk across diverse subpopulations from nhanes IV: A cohort study. PLoS Med. 15, e1002718. doi:10.1371/journal.pmed.1002718
Lu, A. T., Quach, A., Wilson, J. G., Reiner, A. P., Aviv, A., Raj, K., et al. (2019). DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany NY) 11, 303–327. doi:10.18632/aging.101684
Marioni, R. E., McRae, A. F., Bressler, J., Colicino, E., Hannon, E., Li, S., et al. (2018). Meta-analysis of epigenome-wide association studies of cognitive abilities. Mol. Psychiatry 23, 2133–2144. doi:10.1038/s41380-017-0008-y
Mathur, R., Hui, Q., Huang, Y., Gwinn, M., So-Armah, K., Freiberg, M. S., et al. (2019). DNA methylation markers of type 2 diabetes mellitus among male veterans with or without human immunodeficiency virus infection. J. Infect. Dis. 219, 1959–1962. doi:10.1093/infdis/jiz023
Mendelson, M. M., Marioni, R. E., Joehanes, R., Liu, C., Hedman, A. K., Aslibekyan, S., et al. (2017). Association of body mass index with DNA methylation and gene expression in blood cells and relations to cardiometabolic disease: A mendelian randomization approach. PLoS Med. 14, e1002215. doi:10.1371/journal.pmed.1002215
Molejon, M. I., Ropolo, A., Re, A. L., Boggio, V., and Vaccaro, M. I. (2013). The VMP1-Beclin 1 interaction regulates autophagy induction. Sci. Rep. 3, 1055. doi:10.1038/srep01055
Moran, S., Arribas, C., and Esteller, M. (2016). Validation of a DNA methylation microarray for 850, 000 CpG sites of the human genome enriched in enhancer sequences. Epigenomics 8, 389–399. doi:10.2217/epi.15.114
Nishida, N., Nagasaka, T., Nishimura, T., Ikai, I., Boland, C. R., and Goel, A. (2008). Aberrant methylation of multiple tumor suppressor genes in aging liver, chronic hepatitis, and hepatocellular carcinoma. Hepatology 47, 908–918. doi:10.1002/hep.22110
Ogata, H., Chinen, T., Yoshida, T., Kinjyo, I., Takaesu, G., ShiraisHi, H., et al. (2006). Loss of SOCS3 in the liver promotes fibrosis by enhancing STAT3-mediated TGF-beta1 production. Oncogene 25, 2520–2530. doi:10.1038/sj.onc.1209281
Palella, F. J., Baker, R. K., Moorman, A. C., Chmiel, J. S., Wood, K. C., Brooks, J. T., et al. (2006). Mortality in the highly active antiretroviral therapy era: Changing causes of death and disease in the HIV outpatient study. J. Acquir. Immune Defic. Syndr. 43, 27–34. doi:10.1097/01.qai.0000233310.90484.16
Perlmutter, D. H., Dinarello, C. A., Punsal, P. I., and Colten, H. R. (1986). Cachectin/tumor necrosis factor regulates hepatic acute-phase gene expression. J. Clin. Invest. 78, 1349–1354. doi:10.1172/JCI112721
Rickabaugh, T. M., Baxter, R. M., Sehl, M., Sinsheimer, J. S., Hultin, P. M., Hultin, L. E., et al. (2015). Acceleration of age-associated methylation patterns in HIV-1-infected adults. PLoS One 10, e0119201. doi:10.1371/journal.pone.0119201
Ronit, A., Hatleberg, C. I., Ryom, L., Bonnet, F., El-Sadr, W., Reiss, P., et al. (2018). Associations between serum albumin and serious non-AIDS events among people living with HIV. AIDS 32, 1837–1848. doi:10.1097/QAD.0000000000001900
Ronit, A., Sharma, S., Baker, J. V., Mngqibisa, R., Delory, T., Caldeira, L., et al. (2018). Serum albumin as a prognostic marker for serious non-AIDS endpoints in the strategic timing of antiretroviral treatment (START) study. J. Infect. Dis. 217, 405–412. doi:10.1093/infdis/jix350
Rosen, A. D., Robertson, K. D., Hlady, R. A., Muench, C., Lee, J., Philibert, R., et al. (2018). DNA methylation age is accelerated in alcohol dependence. Transl. Psychiatry 8, 182. doi:10.1038/s41398-018-0233-4
Shepard, C. W., Finelli, L., and Alter, M. J. (2005). Global epidemiology of hepatitis C virus infection. Lancet. Infect. Dis. 5, 558–567. doi:10.1016/S1473-3099(05)70216-4
Sidibeh, C. O., Pereira, M. J., Abalo, X. M., J Boersma, G., Skrtic, S., Lundkvist, P., et al. (2018). FKBP5 expression in human adipose tissue: Potential role in glucose and lipid metabolism, adipogenesis and type 2 diabetes. Endocrine 62, 116–128. doi:10.1007/s12020-018-1674-5
Smith, C., Sabin, C. A., Lundgren, J. D., Thiebaut, R., Weber, R., et al. (2010). Factors associated with specific causes of death amongst HIV-positive individuals in the D:A:D Study. AIDS 24, 1537–1548. doi:10.1097/QAD.0b013e32833a0918
Smith, C. J., Ryom, L., Weber, R., Morlat, P., Pradier, C., Reiss, P., et al. (2014). Trends in underlying causes of death in people with HIV from 1999 to 2011 (D:A:D): A multicohort collaboration. Lancet 384, 241–248. doi:10.1016/S0140-6736(14)60604-8
Solomon, O., MacIsaac, J., Quach, H., Tindula, G., Kobor, M. S., Huen, K., et al. (2018). Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children. Epigenetics 13, 655–664. doi:10.1080/15592294.2018.1497386
Stechschulte, L. A., Qiu, B., Warrier, M., Hinds, T. D., Zhang, M., Gu, H., et al. (2016). FKBP51 null mice are resistant to diet-induced obesity and the PPARγ agonist rosiglitazone. Endocrinology 157, 3888–3900. doi:10.1210/en.2015-1996
Sterling, R. K., Lissen, E., Clumeck, N., Sola, R., Correa, M. C., Montaner, J., et al. (2006). Development of a simple noninvasive index to predict significant fibrosis in patients with HIV/HCV coinfection. Hepatology 43, 1317–1325. doi:10.1002/hep.21178
Su, P. F., Lee, T. C., Lin, P. J., Lee, P. H., Jeng, Y. M., Chen, C. H., et al. (2007). Differential DNA methylation associated with Hepatitis B virus infection in hepatocellular carcinoma. Int. J. Cancer 121, 1257–1264. doi:10.1002/ijc.22849
Sun, D., Zhang, T., Su, S., Hao, G., Chen, T., Li, Q. Z., et al. (2019). Body mass index drives changes in DNA methylation: A longitudinal study. Circ. Res. 125, 824–833. doi:10.1161/CIRCRESAHA.119.315397
Vaccaro, M. I., Grasso, D., Ropolo, A., Iovanna, J. L., and Cerquetti, M. C. (2003). VMP1 expression correlates with acinar cell cytoplasmic vacuolization in arginine-induced acute pancreatitis. Pancreatology 3, 69–74. doi:10.1159/000069150
Varela-Rey, M., Woodhoo, A., Martinez-Chantar, M. L., Mato, J. M., and Lu, S. C. (2013). Alcohol, DNA methylation, and cancer. Alcohol Res. 35, 25–35.
Ventham, N. T., Kennedy, N. A., Adams, A. T., Kalla, R., Heath, S., O'Leary, K. R., et al. (2016). Integrative epigenome-wide analysis demonstrates that DNA methylation may mediate genetic risk in inflammatory bowel disease. Nat. Commun. 7, 13507. doi:10.1038/ncomms13507
Wahl, S., Drong, A., Lehne, B., Loh, M., Scott, W. R., Kunze, S., et al. (2017). Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity. Nature 541, 81–86. doi:10.1038/nature20784
Wai, C. T., Greenson, J. K., Fontana, R. J., Kalbfleisch, J. D., Marrero, J. A., Conjeevaram, H. S., et al. (2003). A simple noninvasive index can predict both significant fibrosis and cirrhosis in patients with chronic hepatitis C. Hepatology 38, 518–526. doi:10.1053/jhep.2003.50346
Xu, K., Zhang, X., Wang, Z., Hu, Y., and Sinha, R. (2018). Epigenome-wide association analysis revealed that SOCS3 methylation influences the effect of cumulative stress on obesity. Biol. Psychol. 131, 63–71. doi:10.1016/j.biopsycho.2016.11.001
Xu, X., Chang, X., Xu, Y., Deng, P., Wang, J., Zhang, C., et al. (2020). SAMD14 promoter methylation is strongly associated with gene expression and poor prognosis in gastric cancer. Int. J. Clin. Oncol. 25, 1105–1114. doi:10.1007/s10147-020-01647-4
Zannas, A. S., Wiechmann, T., Gassen, N. C., and Binder, E. B. (2016). Gene-stress-epigenetic regulation of FKBP5: Clinical and translational implications. Neuropsychopharmacology 41, 261–274. doi:10.1038/npp.2015.235
Zhang, X., Asllanaj, E., Amiri, M., Portilla-Fernandez, E., Bramer, W. M., Nano, J., et al. (2021). Deciphering the role of epigenetic modifications in fatty liver disease: A systematic review. Eur. J. Clin. Invest. 51, e13479. doi:10.1111/eci.13479
Zhang, X., Hu, Y., Justice, A. C., Li, B., Wang, Z., Zhao, H., et al. (2017). DNA methylation signatures of illicit drug injection and hepatitis C are associated with HIV frailty. Nat. Commun. 8, 2243. doi:10.1038/s41467-017-02326-1
Zhao, Y., Xue, F., Sun, J., Guo, S., Zhang, H., Qiu, B., et al. (2014). Genome-wide methylation profiling of the different stages of Hepatitis B virus-related hepatocellular carcinoma development in plasma cell-free DNA reveals potential biomarkers for early detection and high-risk monitoring of hepatocellular carcinoma. Clin. Epigenetics 6, 30. doi:10.1186/1868-7083-6-30
Keywords: EWAS, liver disease, HIV, biomarkers, albumin (ALB)
Citation: Titanji BK, Lee M, Wang Z, Chen J, Hui Q, Lo Re III V, So-Armah K, Justice AC, Xu K, Freiberg M, Gwinn M, Marconi VC and Sun YV (2022) Epigenome-wide association study of biomarkers of liver function identifies albumin-associated DNA methylation sites among male veterans with HIV. Front. Genet. 13:1020871. doi: 10.3389/fgene.2022.1020871
Received: 16 August 2022; Accepted: 22 September 2022;
Published: 11 October 2022.
Edited by:
Haitao Wang, National Cancer Institute, United StatesReviewed by:
Qin Xu, National Institute of Allergy and Infectious Diseases (NIH), United StatesYaqiang Cao, National Institutes of Health, United States
Dawei Zhou, Wexner Medical Center, The Ohio State University, United States
Copyright © 2022 Titanji, Lee, Wang, Chen, Hui, Lo Re III, So-Armah, Justice, Xu, Freiberg, Gwinn, Marconi and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yan V. Sun, yan.v.sun@emory.edu
†These authors contributed equally to this work and share first authorship