A genomic variant of ALPK2 is associated with increased liver fibrosis risk in HIV/HCV coinfected women

About the Authors:

Alec T. McIntosh

Roles Formal analysis, Methodology, Project administration, Writing – original draft, Writing – review & editing

Affiliation: Department of Oncology, Georgetown University, Washington, DC, United States of America

ORCID logo https://orcid.org/0000-0002-7764-3608

Renhuizi Wei

Roles Data curation, Formal analysis, Investigation, Software, Visualization

Affiliation: Department of Oncology, Georgetown University, Washington, DC, United States of America

Jaeil Ahn

Roles Formal analysis, Investigation, Project administration, Supervision, Visualization, Writing – review & editing

Affiliation: Department of Biostatistics, Bioinformatics & Biomathematics, Georgetown University Medical Center, Washington, DC, United States of America

ORCID logo https://orcid.org/0000-0001-7998-4759

Brad E. Aouizerat

Roles Data curation, Investigation, Methodology, Resources, Supervision, Writing – review & editing

Affiliation: Bluestone Center for Clinical Research, College of Dentistry, New York University, New York, New York, United States of America

Seble G. Kassaye

Roles Investigation, Resources, Writing – review & editing

Affiliation: Department of Infectious Diseases, Georgetown University Medical Center, Washington, DC, United States of America

Michael H. Augenbraun

Roles Investigation, Resources, Writing – review & editing

Affiliation: Division of Infectious Diseases, Department of Medicine, State University of New York, Downstate Medical Center, Brooklyn, New York, United States of America

Jennifer C. Price

Roles Investigation, Resources

Affiliation: Division of Liver Diseases, Department of Medicine, University of California San Francisco, San Francisco, California, United States of America

Audrey L. French

Roles Investigation, Resources, Writing – review & editing

Affiliation: Division of Infectious Disease, Department of Internal Medicine, Stroger Hospital of Cook County, Chicago, Illinois, United States of America

Stephen J. Gange

Roles Investigation, Resources, Writing – review & editing

Affiliation: Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland, United States of America

ORCID logo https://orcid.org/0000-0001-7842-512X

Kathryn M. Anastos

Roles Investigation, Resources, Writing – review & editing

Affiliation: Departments of Medicine and Epidemiology & Population Health, Albert Einstein College of Medicine, Bronx, New York, United States of America

Radoslav Goldman

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Oncology, Georgetown University, Washington, DC, United States of America, Department of Biochemistry and Molecular & Cell Biology, Georgetown University Medical Center, Washington, DC, United States of America, Clinical Translational Glycoscience Research Center, Georgetown University Medical Center, Washington, DC, United States of America

ORCID logo https://orcid.org/0000-0003-1304-0522

Introduction

The emergence of highly active antiretroviral therapy has transitioned the once acutely fatal human immunodeficiency virus (HIV) infection to a chronic disease. However, the longer survival of persons living with HIV infection presents a new set of morbidities and increased risk for mortality [1–3]. Due to common modes of transmission, those infected with HIV are at higher risk of contracting hepatitis C virus (HCV). Accelerated progression of liver disease in HIV/HCV coinfected patients compared to those with HCV monoinfection is well documented and liver disease has become a leading cause of non-AIDS related death in coinfected individuals [4,5]. The pathogenesis of liver disease in coinfected individuals is multifactorial [6]. Despite substantial progress in identifying risk factors for liver disease progression in coinfected persons, our understanding of risk factors remains incomplete, with genomic factors among those that remain ill-defined. Even with the emergence of direct-acting HCV antivirals, the ability of these agents to regress fibrosis upon HCV clearance remains unclear [7] and the cost of treatment is often inaccessible to at-risk populations [8]. Taken together, these considerations document that liver fibrosis remains a challenge in coinfected and HCV monoinfected patients. In this study, we examined the impact of non-synonymous single nucleotide variants (nsSNV) affecting N-glycosylation both directly at the NxS/T sequons of proteins and indirectly through the enzymes and lectins of glycosylation-related pathways.

Glycosylation is one of the most common and structurally diverse protein modifications and affects protein synthesis, structure, and function [9]. Glycosylation enzymes or the glycoproteins they produce are involved in immune surveillance and host-pathogen interactions [10,11] as well as in the progression of viral liver disease [12]. To elucidate the impacts of glycosylation on the pathogenesis of fibrosis in HIV/HCV coinfected patients, we focused on nsSNV affecting the NxS/T sequons required for the attachment of N-glycans to proteins and thus the potential to change the number of glycans on the protein surface [13,14]. The Women’s Interagency HIV Study (WIHS) is a longitudinal natural history study of HIV infection that features a sufficient number of HIV/HCV coinfected participants and biomarkers of liver disease to evaluate the impact of risk factors for liver disease progression [15], including genetic risk factors [16–18]. Given its longitudinal cohort design, the Fibrosis-4 Index (FIB-4) [19] and the AST to Platelet Ratio Index (APRI) [20] were collected in the WIHS cohort as noninvasive surrogate measures of hepatic fibrosis. These measures were used in lieu of serial tissue biopsy as the means with which to evaluate genetic risk factors in HIV/HCV coinfection [21]. A set of 887 nsSNV were extracted from a genome-wide association study performed in the WIHS cohort. Of these, 278 nsSNV produce an amino acid substitution that alters the potential glycosylation of target proteins by altering the number of N-glycans decorating a protein, which we analyzed in relation to surrogate biomarkers of liver fibrosis. Given that previous work demonstrated that genetic variation in glycosylation-related enzymes and lectins can alter kinetics and binding affinities respectively, the remaining 609 of 887 nsSNV were glycosylation-related proteins evaluated in relation to surrogate biomarkers of liver fibrosis in coinfected participants [22,23].

Methods

Ethics statement

The parent WIHS study and this sub-study conformed to the procedures for informed written consent approved by institutional review boards (IRB) at all sponsoring organizations and to human-experimentation guidelines set forth by the United States Department of Health and Human Services, and finally reviewed and approved by the local IRBs in the Baltimore/DC region. The IRBs at the Chicago site were Cook County CORE Center IRB, Cook County Health & Hospitals System IRB, Northwestern University IRB, Northwestern Memorial Hospital IRB, Rush University IRB, Rush-Presbyterian-St. Lukes IRB, University of Illinois at Chicago IRB. The IRBs at the Los Angeles/Southern California site were University of California at Los Angeles Medical Center IRB, Western IRB, University of Hawaii/Kapi’olani Health Research Institute IRB, University of Southern California IRB. The IRB at the New York Brooklyn site was SUNY Downstate Medical Center Institutional Review Board. The IRBs at the New York Bronx site were Albert Einstein College of Medicine IRB, Beth Israel IRB, Mt. Sinai Medical Center IRB. The IRBs at the San Francisco Bay Area site were Alameda Health System IRB, Sutter Health IRB, University of California, San Francisco IRB, California Committee for the Protection of Human Subjects (CPHS). The IRBs at the Washington District of Columbia site were Georgetown University Medical Center IRB, Howard University IRB, Inova Health Systems IRB, Montgomery County Health Department IRB, Whitman Walker Clinic IRB.

Study population

The WIHS is an active, multicenter prospective study of the natural history of HIV infection among women with or at risk for HIV-infection in the United States. Established in 1994, a total of 4,982 women (3,677 HIV-seropositive) have been enrolled. At semi-annual visits, participants completed socio-demographic and medical questionnaires, laboratory testing, and a limited physical examination. Data included in this analysis include age, race and ethnicity, continuous clinical measures of liver fibrosis (APRI and FIB-4), plasma HIV RNA viral titer, HCV infection status (at all visits), and HCV viral titer (at baseline only) [24]. HCV status was defined as “positive” when participants had both positive HCV antibody and detectable HCV RNA at their baseline visit. Those that were defined as “negative” had no detectable HCV antibody upon testing. Self-reported race and ethnicity was used to define four groups: “White” (non-Hispanic), “African American” (non-Hispanic), “Hispanic”, or “Other”. In addition to self-reported race and ethnicity, genomic estimates of ancestry were derived using principle component analysis of ancestry informative genetic markers [25].

Genotyping

Genotype data in WIHS was generated from genomic DNA from peripheral blood mononuclear cells using the Infinium Omni2.5 BeadChip (Illumina, San Diego, CA, USA) [24]. Of the 10,141 genes (602 glycogene + 9,539 NxS/T-containing genes) known to exist in the human genome, data on 1,029 nsSNVs (698 glycogene + 331 NxS/T) spanning 660 genes (349 glycogene + 311 NxS/T) from 2,120 WIHS participants were available for analysis [13]. Of these 2,120 women, 262 participants were either HCV antibody positive in the absence of detectable HCV RNA (N = 156) or HCV serostatus was never assessed at baseline (N = 106) and therefore excluded (S1 Fig). Of the 331 NxS/T nsSNVs, 44 were mono-allelic (non-polymorphic) and therefore excluded. Among the 287 nsSNVs that remained, missing genotypes across 1,858 patients were imputed with the most common genotype for each nsSNVs (mean participants imputed per nsSNV±SD = 5±6.04). Of the 698 glycogene nsSNVs, 75 were monomorphic and excluded. Among the 623 nsSNVs that remained, missing genotypes across 1,858 patients were imputed with the most common genotype for each nsSNV (mean participants imputed per nsSNV±SD = 5±6.21). A minor allele frequency (MAF) threshold ≥0.001 (≥0.1%) was applied using the MAFs of the coinfected population (N = 661) to further isolate SNVs sufficiently frequent to allow for statistical analysis. Using this method, an additional 9 NxS/T and 14 glycogenes were eliminated from analysis. After applying all exclusion criteria, clinical and genotype (887 glycogene [609 nsSNV, 278 NxS/T nsSNV] spanning 564 genes) data were available for 1,858 women (S1 Fig). The nsSNV reference sequence identifiers (rsID), major (or common) and minor (or rare) alleles, and MAF across serotypes are provided as supplementary tables for the NxS/T (S1 Table) and glycogene (S2 Table) nsSNVs utilized in our analysis.

Serologic markers of fibrosis

Fibrosis-4 Index (FIB-4) and the AST to Platelet Ratio Index (APRI) were used as measures of hepatic fibrosis as described in previous literature [21]. For the FIB-4 index, scores <1.45 and >3.25 indicate a high negative and a high positive predictive value for advanced fibrosis respectively [19]. For APRI, scores <0.5 have a high negative predictive value for liver disease while scores >0.7 and >2 indicate a high positive predictive value for moderate and severe hepatic fibrosis respectively [20]. Genetic polymorphisms were independently analyzed against each continuous surrogate index (i.e., APRI, FIB-4) at baseline to identify variants with stronger associations that would manifest across both current clinical diagnostic resources.

Statistical analysis

For descriptive summaries, continuous variables were summarized using means and standard deviations while categorical variables were summarized using frequency counts and percentages. HIV and HCV RNA status were categorized into 4 serostatus groups: both HIV/HCV noninfected, HIV monoinfected, HCV monoinfected, and HIV/HCV coinfected. Distributions of APRI, FIB-4, and HIV RNA viral load among the four serostatus groups were summarized using descriptive statistics and were compared using chi-square analyses or one-way Analysis of Variance (ANOVA). APRI and FIB-4 scores used for analysis were obtained at baseline visits (pre-2003) at a time prior to the broad use of HCV therapy. For each variant from the 278 NxS/T and 609 glycogene nsSNV that met inclusion criteria, a separate multiple variable linear regression model was constructed for the HIV/HCV coinfected patient population for continuous APRI and FIB-4 outcomes, respectively. Explanatory variables included each variant, age, race and ethnicity, HIV viral load, and HCV viral load where HIV and HCV viral load were normalized using log2(x+1) transformation. Model results were near indistinguishable when adjusted for either HCV viral load or HCV status. Genomic estimates of racial and ethnicity were estimated using principle component analysis of 185 ancestry informative markers (SNV) selected to differentiate major racial and ethnic groups in the cohort (i.e., European, African, Hispanic) [25] The first three principle components (PC1, PC2, PC3), which explain >90% variation, were utilized to adjust for genomic estimates of race and ethnicity in the aforementioned linear regression models. HIV and HCV viral load were normalized using log2(x+1) transformation. To reassure our findings were not sensitive to confounding factors of liver fibrosis, we performed a sensitivity analysis on the linear regression models with adjustment for additional factors including Hepatitis B (HBV) infection status (categorical variable; positive or negative for chronic infection by serology), participant reported alcohol usage (categorical variable; abstainers [0 drinks/week], mild users [1–7 drinks/week], or heavy users [>7 drinks/week]), and Body Mass Index (BMI, a surrogate marker for Non-Alcoholic Fatty Liver Disease (NAFLD)). The false discovery rate (FDR) was used to account for multiple testing (i.e., 887 linear regression each for FIB-4 and APRI). All reported p-values are two-sided, and SNV associations with FDR<0.05 were considered statistically significant. To investigate whether the impact of SNV varies across the four serogroups, we examined SNV*serostatus interactions and computed 95% confidence intervals for each SNV genotype*serostatus configuration. Statistical analyses were performed using SAS version 9.4 (SAS Institute, Cary, NC, USA) and R version 3.3 (R Core Team, Vienna, Austria).

Results

Of the 75% of women who were HIV seropositive (1386 of 1858), 52% (n = 725) were HIV monoinfected and 48% (n = 661) were HIV/HCV coinfected (Table 1).

[Figure omitted. See PDF.]

Table 1. Baseline population characteristics of patients by serostatus.

https://doi.org/10.1371/journal.pone.0247277.t001

The remaining 25% of women (N = 472) were HCV/HIV noninfected and 22% (N = 104) were HCV monoinfected. The proportion of non-Hispanic African-Americans differed significantly between the four serogroups (P<0.001) (Table 1). Post-hoc contrasts indicate that, in comparison to non-Hispanic African Americans, Hispanics were less likely (OR = 0.66 [95% CI, 0.499 to 0.881], PFDR = 0.03) and non-Hispanic Whites were more likely (OR = 1.82 [95% CI, 1.272 to 2.609], PFDR = 0.006) to be in the coinfected serogroup as compared to the HIV monoinfected group. HCV-infected women were older on average (P<0.001) than those not infected with the virus (Table 1). Post-hoc contrasts indicate that in comparison to the noninfected serogroup, older subjects were more likely to be monoinfected with HIV (OR = 1.08 [95% CI, 1.059 to 1.100], PFDR <0.001), monoinfected with HCV (OR = 1.21 [95% CI, 1.172, to 1.252], PFDR <0.001) or coinfected (OR = 1.22 [95% CI, 1.197 to 1.252], PFDR <0.001). Distribution of the APRI and FIB-4 fibrotic indexes, stratified by their clinically relevant cutoffs, also differed between the groups (P<0.001) (Table 1). Six pairwise comparisons of mean liver scores were performed between designated serogroups for each fibrosis metric (Fig 1). On average, coinfected participants had higher APRI and FIB-4 measures when compared to the HCV monoinfected population (APRI, 0.36, 95% CI, 0.064 to 0.651; FIB-4, 0.78, 95% CI, 0.375 to 1.184) and the HIV monoinfected population (APRI, 0.64, 95% CI, 0.487 to 0.786; FIB-4, 1.19, 95% CI, 0.984 to 1.397) (Fig 1). All other comparisons between serogroups were found to be statistically significant (95% confidence interval) except the mean APRI scores for noninfected vs. HIV monoinfected (-0.15, 95% CI, -0.325 to 0.031) and HIV monoinfected vs. HCV monoinfected (-0.28, 95% CI, -0.571 to 0.012) which follow the expected trend but do not achieve statistical significance (Fig 1). With regard to possible confounding factors, BMI, HBV serological status, and alcohol usage are also shown stratified between the four HIV/HCV serostatuses (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Mean liver fibrosis scores for each viral serotype.

The mean APRI and FIB-4 scores at baseline visit + SD are shown for 1,858 women stratified by definitive viral serostatus. The number of patients (n) for each viral serostatus are indicated in Table 1. Pairwise comparisons of liver outcomes were performed for all viral-serotype group-comparisons within each liver metric. Apart from those comparisons labeled not significant (ns), all other comparisons within fibrosis index aside from those labeled as significant were also below the 0.05 p-value threshold.

https://doi.org/10.1371/journal.pone.0247277.g001

Utilizing the larger cohort for race and age adjustment helped to ensure that we had the population size to account for the impact of these variables on liver fibrosis. Whereas the first principle component of the PCA of genomic markers of ancestry (PC1) appeared to adequately differentiate African-Americans (non-Hispanic), Hispanics, and Whites (non-Hispanic) (Fig 2A), we took a conservative approach and included the first three principle components (PC1, PC2, PC3). Although adjusting for race/ethnicity using self-report identified the same genetic associations, we opted to employ the 3 PCs to better account for confounding that can occur with self-reported race and ethnicity. For example, many participants self-reporting as “Hispanic” or “Other” in the coinfected subgroup present with a genetic background more consistent with the African-American (non-Hispanic) and White (non-Hispanic) clusters (Fig 2B).

[Figure omitted. See PDF.]

Fig 2. Genomic estimates of race/ethnicity using principle component analysis.

The 1,858 WIHS participants utilized in this study were stratified using 185 ancestral SNV markers from across the genome. PC1 and PC2 are depicted above. Genomic estimates of race and ethnicity were compared to reported race and ethnicity (A) and the distribution of HIV/HCV coinfected participants of focus are displayed (B).

https://doi.org/10.1371/journal.pone.0247277.g002

After adjustment for multiple testing, seven nsSNV mapping to separate genes met the a priori criterion of PFDR <0.05 for at least one of the biomarkers of liver fibrosis (Table 2). Of these, only rs52828316 (MAN2A2) was significant by PFDR for both indices. For all nsSNV, two copies of the minor allele (i.e., minor allele homozygotes) was associated with increases in APRI and/or FIB-4 (Table 2). Upon sensitivity analysis of additional adjustable factors of HBV Status, alcohol usage, and BMI (N = 616), all statistically significant nsSNVs were maintained, and APRI and FIB-4 estimates for all but rs1800472 (TGFB1) were within 10% of that of the initial model (S3 Table). Two additional nsSNV were found to be significant when additionally adjusting for HBV status, alcohol, and BMI, namely rs3745925 (MADCAM1) and rs2307145 (IL12RB2) (S3 Table). ALPK2 rs3809973 was evaluated further in relation to hepatic fibrosis among the four viral serogroups because it was the only variant sufficiently frequent (n>100) to result in all three genotypic groups (i.e., major allele homozygotes, heterozygotes, minor allele homozygotes), and the ALPK2 variant was not sensitive to adjustable factors on the linear regression model.

[Figure omitted. See PDF.]

Table 2. Baseline screen of multiple linear regression models.

https://doi.org/10.1371/journal.pone.0247277.t002

In order to see if the ALPK2 variant’s tentative association with increased liver fibrosis in coinfected women was preferentially impacted by the effects of either virus, we analyzed each genotype of the ALPK2 variant across all viral serostatuses utilized in the cohort. For the ALPK2 rs3809973, HIV/HCV coinfected participants who were homozygous for the minor allele had significantly higher mean APRI and FIB-4 scores relative to the coinfected participants homozygous for the major allele (APRI, 0.61, 95% CI, 0.334 to 0.875; FIB-4, 0.74, 95% CI, 0.336 to 1.144) and compared to those HCV monoinfected participants homozygous for the minor allele (APRI, 0.79, 95% CI, 0.370 to 1.200; FIB-4, 1.24, 95% CI, 0.820 to 1.650) (Fig 3A and 3B). Evaluation of the association of ALPK2 rs3809973 with APRI and FIB-4 by racial and ethnic group (i.e., non-Hispanic African-American, non-Hispanic Caucasian, Hispanic) revealed similar patterns of association with one exception in non-Hispanic African-American where no difference by genotypic group was observed with FIB-4 (data not shown). The heterozygotes of coinfected individuals displayed no significant increases in APRI or FIB-4 relative to the coinfected major allele homozygotes (APRI, 0.10, 95% CI, -0.128 to 0.332; FIB-4, 0.15, 95% CI, -0.195 to 0.492). These coinfected population findings for the ALPK2 nsSNV (rs3809973) were observed against a background of similar CD4+ T-cell percentage, measured relative to other leukocyte types, and detectable HIV viral loads at baseline (Fig 3C and 3D). HIV monoinfected participants homozygous for the minor allele displayed significantly increased FIB-4 (0.45, 95% CI, 0.028 to 0.879) scores relative to noninfected participants homozygous for the minor allele, but the finding was not significant for the same comparison using APRI (0.00, 95% CI, -0.422 to 0.429) (Fig 3A and 3B). For the HIV monoinfected genotype, a slight significant difference in the CD4+ T-cell percentage (P = 0.007) accompanied this FIB-4 finding. Impacts of variant ALPK2 allele burden between genotypes in the noninfected and HCV monoinfected serogroups were unremarkable (Fig 3A and 3B). A parallel set of analyses were conducted with adjustment for the additional factors of HBV status, alcohol use, and BMI. Overall, little differences from the present analysis was observed (S2 Fig). However, when comparing the heterozygotes of coinfected individuals to the coinfected major allele homozygotes, FIB-4 was found to be significantly increased (0.19, 95% CI, 0.006 to 0.382) with the APRI measure of the same comparison borderline insignificant (0.17, 95% CI, -0.022 to 0.353) when adjusting for these additional factors (S2 Fig).

[Figure omitted. See PDF.]

Fig 3. rs3809973 (ALPK2) impact on APRI and FIB-4 outcomes.

Mean baseline APRI (A) and FIB-4 (B) score with 95% confidence interval shown for the allele pairs of each serotype (Total N = 1,858). The number of samples representing each genotype for the respective serogroups is displayed as follows: homozygous for the major allele (green), heterozygous (orange), and homozygous for the minor allele (blue). Comparisons were shown, for reference, relative to the major allele homozygote of the noninfected serogroup. ANOVA was used to compare relative CD4+ T-cell percentages (C) or HIV/HCV viral loads (D) between the genotypes of each serogroup. Asterisks indicate two-sided P-values below 0.05 (*) and 0.01 (**) respectively.

https://doi.org/10.1371/journal.pone.0247277.g003

Discussion

This study identified several novel genetic associations among glycogenes and biomarkers of liver fibrosis. Average liver fibrosis scores in the coinfected serogroup analyzed at baseline were significantly higher than that of either HIV or HCV monoinfected participants as shown in Fig 1, recapitulating the fibrotic trends described in the literature for coinfected populations [6,26]. Of the 887 nsSNV assessed that either directly (nsSNV affecting NxS/T sequons) or indirectly (nsSNV in glycosylation or lectin genes) alter glycosylation pathways or products, we found seven nsSNV that were associated with an increased risk of hepatic fibrosis among the HIV/HCV coinfected population. We observed higher APRI and FIB-4 values (indicative of greater fibrosis) in coinfected participants homozygous for the ALPK2 rs3809973 minor allele when compared with participants homozygous for the ALPK2 rs3809973 major allele irrespective of HIV and HCV viral load, BMI, HBV Status, or alcohol usage (Fig 3 and S2 Fig). As elevations in viral HIV titers have been correlated with increased liver injury in the coinfected [27], our findings suggest that the genetic risk of hepatic fibrosis among coinfected individuals carrying the ALPK2 rs3809973 risk allele is not confounded by these titers. The lack of association in coinfected participants between the ALPK2 variant heterozygote and increased liver fibrosis suggested the need for variant homozygosity for the detrimental impacts of the variant to manifest in the coinfected and conforms with a recessive mode of inheritance. As shown in Fig 3, among HIV monoinfected participants, the relative CD4+ T-cell percentage among leukocytes was inversely associated with FIB-4. Furthermore, coinfected participants that were homozygous for the minor allele were significantly increased in terms of fibrosis risk when compared to HCV monoinfected participants of the same genotype; suggesting a need for general HIV-associated immunosuppression to produce the detrimental fibrosis risk increase. Our comparisons of ALPK2 rs3809973 among serogroups at baseline therefore leads us to speculate that the risk allele’s impact on liver fibrosis is reliant on both HCV-mediated liver damage and its perpetuation by HIV-mediated immunosuppression.

The ALPK2 gene and rs3809973 have not been linked to a pathological liver phenotype to date. Although the association of ALPK2 rs3809973 with fibrosis in the setting of HIV/HCV coinfection warrants replication, some circumstantial evidence suggests a potential role for ALPK2 in risk for viral hepatic fibrosis. Hepatic Stellate cells (HSC) lining the perisinusoidal space of the liver are the fibrogenic cell type responsible for extracellular matrix deposition and subsequent fibrosis in the setting of hepatic inflammation. Although this fibrogenic response is designed to protect against liver injury and typically reverses after the hepatic insult subsides, progressive inflammation and the chronic activation of HSCs can lead to cirrhosis and an increased risk of developing hepatocellular carcinoma [28,29]. Signaling molecules and pathways involved in HSC activation have been described, but most notable is the expression of toll-like receptors on HSCs enables their activation upon exposure to structurally conserved microbial-derived products such as lipopolysaccharide [30]. Although incompletely characterized, the ALPK2 protein was first identified in cardiomyocytes as an integral regulator of heart development [31]. However, most studies of ALPK2 have focused on the pathological processes of the gastrointestinal tract. ALPK2 is associated with luminal apoptosis in colorectal cancer cell lines [32]. Luminal shedding is a phenomenon of the enteric innate immune system designed to maintain the function of the gut barrier by preventing the invasion of virulent microbes into systemic circulation. The increased apoptosis of enterocytes in the absence of their concomitant replacement can yield an increasingly permeable intestinal membrane that leaves the host vulnerable to augmented translocation of microbial products [33]. Since increased microbial translocation is a hallmark of HIV infection and lipopolysaccharide activates HSCs [34], any augmentation of this translocation process would exacerbate liver fibrosis in the context of existing hepatitis. In addition, a genome-wide association study of inflammatory bowel disease, a condition in which pathological epithelial shedding is observed, found the mRNA of ALPK2 to be up-regulated in inflamed mucosa compared to control samples [35]. Therefore, we speculate that ALPK2 rs3809973 could enhance liver fibrosis by interfering with enteric immunity. Even though this variant generates a novel NxS/T sequon, augmentation of glycan attachment onto the amino acid substitution encoded by the minor allele has not been demonstrated. Structural studies indicate that the variant does not interrupt the kinase domain nor any residues known to be post-translationally modified, glycosylation or otherwise [31,36]. Therefore, the impact of the K829N substitution encoded by rs3809973 on ALPK2 structure and function warrant investigation.

Although the low minor allele frequency of the remaining 6 nsSNV associated with increased liver fibrosis (i.e., CCR7 rs2228015, GLT8D2 rs17035120, MAN2A2 rs52828316, MBL2 rs1800450, MCOLN2 rs77452813, and TGFB1 rs1800472) precluded detailed modeling, their previous associations with liver disease, metabolism, and viral immunity in other literature provide impetus for further investigation. TGFB1 has long been linked to enhanced hepatocyte destruction and HSC activation [37]. GLT8D2 and MAN2A2 have been implicated in non-alcoholic fatty liver disease [38] and coinfected liver disease [39], respectively. Finally, CCR7, MBL2, and MCOLN2 have all been connected to the modulation of host cellular and immune responses in viral infection [40–42]. Taken together, all seven nsSNV and their cognate genes have plausible roles in liver disease pathology and warrant replication in an independent sample(s).

We note that there are limitations to this study. As this is an exclusively female cohort, future studies should aim to validate these findings in men. Although we identified that there are potentially more than 42,000 nsSNV that either directly altered N-linked glycosylation sites (NxS/T) or that may indirectly alter the function of enzymes or lectins interacting with glycoproteins, only a subset (2.5%) were available for analysis in the cohort due to a combination of low nsSNV allele frequency in the population and/or lack of inclusion of the majority of nsSNV on the commercial array used. A number of factors hindered our attempts at a longitudinal assessment of the genetic risk factors on the coinfected population. For one, as the standard of care for HIV treatment changed for coinfected WIHS participants longitudinally (S3A and S3B Fig), adjusting for shifting antiviral treatment regimens at different time points to meet this standard of care was a challenge. Secondly, as this cohort was initially designed to investigate HIV progression, the HCV virologic status was only determined at baseline and HCV clearance was not assessed further, thereby making longitudinal analyses of HCV-infection status susceptible to misclassification bias. In addition, we had no way of determining the duration of HCV infection at the time of serologic assessment. Along similar lines, survival bias and the availability of more effective and less toxic HIV treatments for new enrollees over time in WIHS may have attenuated the relationship between susceptibility alleles and fibrosis [27]. The rate of loss to follow-up was also problematic in that nearly one third of the coinfected participants at baseline were no longer in the cohort after four visits (2 years) (S3C and S3D Fig). With regard to our sensitivity analysis, only 32 patients (2%) were chronically infected with HBV and this may not have been sufficient power to account for the influence of HBV on liver fibrosis in the linear regression model. Although BMI is a practical means of estimating obesity, it was used in the sensitivity analysis as an imperfect measure of NAFLD and may not fully account for the impact of the condition of liver fibrosis. Along similar lines, alcohol usage was self-reported and may be subject to under reporting in the sensitivity analysis, as the largest fraction of subjects from each of the virally infected serogroups in Table 1 reported to abstain from all alcohol consumption. Finally, APRI and FIB-4 scores are correlated with liver fibrosis risk, but both are imperfect surrogate measures. The utility of these surrogate fibrosis metrics has been validated against biopsy and other measures (i.e. transient elastography) in previous studies [43], but were not validated by such methods at present.

Given the goal of improving the prognosis of liver disease in the HIV/HCV coinfected population, we have identified candidate genes that may participate in hepatic fibrosis. Although these genetic associations require replication, demonstration of the impact of each of the seven nsSNVs on N-linked glycosylation of the site and its cognate protein warrants investigation. Longitudinal HCV testing and assessment of viremia along with additional studies in other HIV/HCV coinfected cohorts may permit better stratification of serogroups and to perform longitudinal analysis of these candidate gene nsSNV.

Supporting information

S1 Fig. Genotyping and data processing workflow.

https://doi.org/10.1371/journal.pone.0247277.s001

(TIF)

S2 Fig. rs3809973 (ALPK2) impact on APRI and FIB-4 outcomes with additional HBV, alcohol, and BMI adjustment.

Mean baseline APRI (A) and FIB-4 (B) score with 95% confidence interval shown for the allele pairs of each serotype (Total N = 1,790). The number of samples representing each genotype for the respective serogroups is displayed as follows: homozygous for the major allele (green), heterozygous (orange), and homozygous for the minor allele (blue). Comparisons were shown, for reference, relative to the major allele homozygote of the noninfected serogroup. Different from Fig 3, Outcome was additionally adjusted for factors influencing liver fibrosis including HBV infection status, Alcohol usage, and BMI. ANOVA was used to compare relative CD4+ T-cell percentages (C) or HIV/HCV viral loads (D) between the genotypes of each serogroup. Asterisks indicate two-sided P-values below 0.05 (*) and 0.01 (**) respectively.

https://doi.org/10.1371/journal.pone.0247277.s002

(TIF)

S3 Fig. Longitudinal trends characterizing the HIV/HCV coinfected population.

The values are summarized annually at every second visit. (A) Composition of the anti-HIV drug regimens changing in line with treatment guidelines; no therapy (red), monotherapy (blue), combination therapy (green), or highly active antiretroviral therapy (HAART) (purple). (B) Mean CD4+ percentages (black symbols) and mean HIV viral load (red symbols). (C) The median APRI and (D) FIB-4 indexes.

https://doi.org/10.1371/journal.pone.0247277.s003

(TIF)

S1 Table. NxS/T nsSNV utilized in analysis.

https://doi.org/10.1371/journal.pone.0247277.s004

(XLSX)

S2 Table. Glycogene nsSNV utilized in analysis.

https://doi.org/10.1371/journal.pone.0247277.s005

(XLSX)

S3 Table. Baseline screen of multiple linear regression models with additional adjustment for Alcohol, HBV, and BMI (NAFLD).

https://doi.org/10.1371/journal.pone.0247277.s006

(DOCX)

Acknowledgments

Data in this manuscript were collected by the Women’s Interagency HIV Study, now the MACS/WIHS Combined Cohort Study (MWCCS). The contents of this publication are solely the responsibility of the authors and do not represent the official views of the National Institutes of Health (NIH). MWCCS (Principal Investigators): Atlanta CRS (Ighovwerha Ofotokun, Anandi Sheth, and Gina Wingood), U01-HL146241; Bronx CRS (Kathryn Anastos and Anjali Sharma), U01-HL146204; Brooklyn CRS (Deborah Gustafson and Tracey Wilson), U01-HL146202; Data Analysis and Coordination Center (Gypsyamber D’Souza, Stephen Gange and Elizabeth Golub), U01-HL146193; Chicago-Cook County CRS (Mardge Cohen and Audrey French), U01-HL146245; Connie Wofsy Women’s HIV Study, Northern California CRS (Bradley Aouizerat and Phyllis Tien), U01-HL146242; Metropolitan Washington CRS (Seble Kassaye and Daniel Merenstein), U01-HL146205; Miami CRS (Maria Alcaide, Margaret Fischl, and Deborah Jones), U01-HL146203; UAB-MS CRS (Mirjam-Colette Kempf and Deborah Konkle-Parker), U01-HL146192; UNC CRS (Adaora Adimora), U01-HL146194. Additionally, the authors would like to thank Jing Wu and Yiwen Wang for their data assembly and analysis efforts.

Citation: McIntosh AT, Wei R, Ahn J, Aouizerat BE, Kassaye SG, Augenbraun MH, et al. (2021) A genomic variant of ALPK2 is associated with increased liver fibrosis risk in HIV/HCV coinfected women. PLoS ONE 16(3): e0247277. https://doi.org/10.1371/journal.pone.0247277

References

1. Sherman KE, Peters MG, Thomas D. Human immunodeficiency virus and liver disease: A comprehensive update. Hepatol Commun. 2017;1(10):987–1001. pmid:30838978

2. Cahill S, Valadéz R. Growing older with HIV/AIDS: new public health challenges. Am J Public Health. 2013;103(3):e7–e15. pmid:23327276

3. Palella FJ, Baker RK, Moorman AC, et al. Mortality in the highly active antiretroviral therapy era: changing causes of death and disease in the HIV outpatient study. J Acquir Immune Defic Syndr. 2006;43(1):27–34. pmid:16878047

4. Salmon-Ceron D, Lewden C, Morlat P, et al. Liver disease as a major cause of death among HIV infected patients: role of hepatitis C and B viruses and alcohol. J Hepatol. 2005;42(6):799–805. pmid:15973779

5. Avihingsanon A, Jitmitraparp S, Tangkijvanich P, et al. Advanced liver fibrosis by transient elastography, fibrosis 4, and alanine aminotransferase/platelet ratio index among Asian hepatitis C with and without human immunodeficiency virus infection: role of vitamin D levels. J Gastroenterol Hepatol. 2014;29(9):1706–1714. pmid:24730732

6. Chrysanthidis T, Loli G, Metallidis S, Germanidis G. Mechanisms of accelerated liver fibrosis in HIV-HCV coinfection. AIDS Rev. 2017;19(3):148–155. pmid:28926561

7. Lledó GM, Carrasco I, Benítez-Gutiérrez LM, et al. Regression of liver fibrosis after curing chronic hepatitis C with oral antivirals in patients with and without HIV coinfection. AIDS. 2018;32(16):2347–2352. pmid:30096074

8. Henry B. Drug Pricing & Challenges to Hepatitis C Treatment Access. J Health Biomed Law. 2018;14:265–283. pmid:30258323

9. Schwarz F, Aebi M. Mechanisms and principles of N-linked protein glycosylation. Curr Opin Struct Biol. 2011;21(5):576–582. pmid:21978957

10. Rudd PM. Glycosylation and the Immune System. Science. 2001;291(5512):2370–2376. pmid:11269318

11. Crocker PR, Paulson JC, Varki A. Siglecs and their roles in the immune system. Nat Rev Immunol. 2007;7(4):255–266. pmid:17380156

12. Van Breedam W, Pöhlmann S, Favoreel HW, de Groot RJ, Nauwynck HJ. Bitter-sweet symphony: glycan-lectin interactions in virus biology. FEMS Microbiol Rev. 2014;38(4):598–632. pmid:24188132

13. Mazumder R, Morampudi KS, Motwani M, Vasudevan S, Goldman R. Proteome-Wide Analysis of Single-Nucleotide Variations in the N-Glycosylation Sequon of Human Genes. Dalby AR, ed. PLoS ONE. 2012;7(5):e36212. pmid:22586465

14. Fan Y, Hu Y, Yan C, et al. Loss and gain of N-linked glycosylation sequons due to single-nucleotide variation in cancer. Sci Rep. 2018;8(1):4322. pmid:29531238

15. Bacon MC, von Wyl V, Alden C, et al. The Women’s Interagency HIV Study: an observational cohort brings clinical sciences to the bench. Clin Diagn Lab Immunol. 2005;12(9):1013–1019. pmid:16148165

16. Sarkar M, Aouzierat B, Bacchetti P, et al. Association of IFNL3 and IFNL4 polymorphisms with liver-related mortality in a multiracial cohort of HIV/HCV-coinfected women. J Viral Hepat. 2015;22(12):1055–1060. pmid:26115445

17. French AL, Martin JW, Evans CT, et al. Macrophage Activation and the Tumor Necrosis Factor Cascade in Hepatitis C Disease Progression Among HIV-Infected Women Participating in the Women’s Interagency HIV Study. J Acquir Immune Defic Syndr. 2017;76(4):438–444. pmid:29077674

18. Urabe Y, Ochi H, Kato N, et al. A genome-wide association study of HCV-induced liver cirrhosis in the Japanese population identifies novel susceptibility loci at the MHC region. J Hepatol. 2013;58(5):875–882. pmid:23321320

19. Sterling RK, Lissen E, Clumeck N, et al. Development of a simple noninvasive index to predict significant fibrosis in patients with HIV/HCV coinfection. Hepatology. 2006;43(6):1317–1325. pmid:16729309

20. Lin Z-H, Xin Y-N, Dong Q-J, et al. Performance of the aspartate aminotransferase-to-platelet ratio index for the staging of hepatitis C-related fibrosis: an updated meta-analysis. Hepatology. 2011;53(3):726–736. pmid:21319189

21. Kuniholm MH, Gao X, Xue X, et al. The relation of HLA genotype to hepatitis C viral load and markers of liver fibrosis in HIV-infected and HIV-uninfected women. J Infect Dis. 2011;203(12):1807–1814. pmid:21606539

22. Hansen L, Lind-Thomsen A, Joshi HJ, et al. A glycogene mutation map for discovery of diseases of glycosylation. Glycobiology. 2015;25(2):211–224. pmid:25267602

23. Bouwman LH, Roos A, Terpstra OT, et al. Mannose binding lectin gene polymorphisms confer a major risk for severe infections after liver transplantation. Gastroenterology. 2005;129(2):408–414. pmid:16083697

24. Shendre A, Wiener HW, Irvin MR, et al. Genome-wide admixture and association study of subclinical atherosclerosis in the Women’s Interagency HIV Study (WIHS). PLoS ONE. 2017;12(12):e0188725. pmid:29206233

25. Frasco MA, Mack WJ, Van Den Berg D, et al. Underlying genetic structure impacts the association between CYP2B6 polymorphisms and response to efavirenz and nevirapine. AIDS. 2012;26(16):2097–2106. pmid:22951632

26. Mohsen AH, Easterbrook PJ, Taylor C, et al. Impact of human immunodeficiency virus (HIV) infection on the progression of liver fibrosis in hepatitis C virus infected patients. Gut. 2003;52(7):1035–1040. pmid:12801963

27. Bräu N, Salvatore M, Ríos-Bedoya CF, et al. Slower fibrosis progression in HIV/HCV-coinfected patients with successful HIV suppression using antiretroviral therapy. J Hepatol. 2006;44(1):47–55. pmid:16182404

28. Lee YA, Wallace MC, Friedman SL. Pathobiology of liver fibrosis: a translational success story. Gut. 2015;64(5):830–841. pmid:25681399

29. O’Rourke JM, Sagar VM, Shah T, Shetty S. Carcinogenesis on the background of liver fibrosis: Implications for the management of hepatocellular cancer. World J Gastroenterol. 2018;24(39):4436–4447. pmid:30357021

30. Tsuchida T, Friedman SL. Mechanisms of hepatic stellate cell activation. Nat Rev Gastroenterol Hepatol. 2017;14(7):397–411. pmid:28487545

31. Ryazanov AG, Pavur KS, Dorovkov MV. Alpha-kinases: a new class of protein kinases with a novel catalytic domain. Curr Biol. 1999;9(2):R43–45. pmid:10021370

32. Yoshida Y, Tsunoda T, Doi K, et al. ALPK2 is crucial for luminal apoptosis and DNA repair-related gene expression in a three-dimensional colonic-crypt model. Anticancer Res. 2012;32(6):2301–2308. pmid:22641666

33. Yu LC-H, Wang J-T, Wei S-C, Ni Y-H. Host-microbial interactions and regulation of intestinal epithelial barrier function: From physiology to pathology. World J Gastrointest Pathophysiol. 2012;3(1):27–43. pmid:22368784

34. Brenchley JM, Price DA, Schacker TW, et al. Microbial translocation is a cause of systemic immune activation in chronic HIV infection. Nat Med. 2006;12(12):1365–1371. pmid:17115046

35. Palmieri O, Creanza TM, Bossa F, et al. Genome-wide Pathway Analysis Using Gene Expression Data of Colonic Mucosa in Patients with Inflammatory Bowel Disease. Inflamm Bowel Dis. 2015;21(6):1260–1268. pmid:25901971

36. Hornbeck PV, Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 2015;43(Database issue):D512–520. pmid:25514926

37. Dooley S, ten Dijke P. TGF-β in progression of liver disease. Cell Tissue Res. 2012;347(1):245–256. pmid:22006249

38. Zhan Y, Zhao F, Xie P, et al. Mechanism of the effect of glycosyltransferase GLT8D2 on fatty liver. Lipids Health Dis. 2015;14:43. pmid:25952508

39. Gupta P, Liu B, Wu JQ, et al. Genome-wide mRNA and miRNA analysis of peripheral blood mononuclear cells (PBMC) reveals different miRNAs regulating HIV/HCV co-infection. Virology. 2014;450–451:336–349. pmid:24503097

40. Hayasaka H, Kobayashi D, Yoshimura H, Nakayama EE, Shioda T, Miyasaka M. The HIV-1 Gp120/CXCR4 axis promotes CCR7 ligand-dependent CD4 T cell migration: CCR7 homo- and CCR7/CXCR4 hetero-oligomer formation as a possible mechanism for up-regulation of functional CCR7. PLoS ONE. 2015;10(2):e0117454. pmid:25688986

41. Eisen S, Dzwonek A, Klein NJ. Mannose-binding lectin in HIV infection. Future Virol. 2008;3(3):225–233. pmid:21218140

42. Rinkenberger N, Schoggins JW. Mucolipin-2 Cation Channel Increases Trafficking Efficiency of Endocytosed Viruses. mBio. 2018;9(1). pmid:29382735

43. Papadopoulos N, Vasileiadi S, Papavdi M, et al. Liver fibrosis staging with combination of APRI and FIB-4 scoring systems in chronic hepatitis C as an alternative to transient elastography. Ann Gastroenterol. 2019;32(5):498–503. pmid:31474797