|Home | About | Journals | Submit | Contact Us | Français|
CXCL13 and CXCR5 are a chemokine and receptor pair whose interaction is critical for naïve B cell trafficking and activation within germinal centers. We sought to determine whether CXCL13 levels are elevated prior to HIV-associated non-Hodgkin B-cell lymphoma (AIDS-NHL), and whether polymorphisms in CXCL13 or CXCR5 are associated with AIDS-NHL risk and CXCL13 levels in a large cohort of HIV-infected men.
CXCL13 levels were measured in sera from 179 AIDS-NHL cases and 179 controls at three time-points. TagSNPs in CXCL13 (n=16) and CXCR5 (n=11) were genotyped in 183 AIDS-NHL cases and 533 controls. Odds ratios (OR) and 95% confidence intervals (CIs) for the associations between one unit increase in log CXCL13 levels and AIDS-NHL, as well as tagSNP genotypes and AIDS-NHL, were computed using logistic regression. Mixed linear regression was used to estimate mean ratios (MR) for the association between tagSNPs and CXCL13 levels.
CXCL13 levels were elevated >3 years (OR=3.24, 95% CI=1.90–5.54), 1–3 years (OR=3.39, 95% CI=1.94–5.94) and 0–1 year (OR=3.94, 95% CI=1.98–7.81) prior to an AIDS-NHL diagnosis. The minor allele of CXCL13 rs355689 was associated with reduced AIDS-NHL risk (ORTCvsTT=0.65; 95% CI=0.45–0.96) and reduced CXCL13 levels (MRCCvsTT=0.82, 95% CI=0.68–0.99). The minor allele of CXCR5 rs630923 was associated with increased CXCL13 levels (MRAAvsTT=2.40, 95% CI=1.43–4.50).
CXCL13 levels were elevated preceding an AIDS-NHL diagnosis, genetic variation in CXCL13 may contribute to AIDS-NHL risk, and CXCL13 levels may be associated with genetic variation in CXCL13 and CXCR5.
CXCL13 may serve as a biomarker for early AIDS-NHL detection.
Among HIV-infected people, B cell non-Hodgkin lymphoma (AIDS-NHL) is currently the most commonly diagnosed cancer (1). The depletion of CD4+ T cells during chronic HIV infection contributes to the development of some AIDS-NHL, such as primary central nervous system lymphoma (PCNSL), through the loss of immunoregulatory control over Epstein-Barr virus (EBV) infected B cells (2). Some of the more common non-PCNSL (or systemic) AIDS-NHL subtypes, including diffuse large B cell lymphoma (DLBCL), however, tend to occur in individuals with relatively high CD4+ T cell numbers and a smaller proportion of these tumors are EBV-positive (2, 3). These cancers are thought to arise from downstream effects of chronic B cell activation, which is also a well-documented consequence of chronic HIV infection (4–8).
Chronic B cell activation results in an over-expression of the DNA-editing enzyme activation induced cytidine deaminase (AID), the actions of which are normally reserved for the highly regulated immunoglobulin gene (Ig) diversification reactions of class switch recombination and somatic hypermutation in germinal centers (9). High levels of AID expression can result in aberrant mutations to non-Ig genes and may cause the seminal translocations and mutations seen in AIDS-NHL (10, 11). Studies have shown significant associations between elevated levels of B cell activation-associated biomarkers, including many cytokines and soluble receptors, and AIDS-NHL risk (12–18). Although the exact mechanism driving chronic B cell activation is not known, one possibility that remains to be explored is the potential role of chemokines that regulate B cell trafficking.
CXCL13 (BCA-1, BLC) is a chemokine produced by follicular dendritic and T helper (particularly Th17) cells, in secondary lymphoid organs (19, 20). When bound to its receptor, CXCR5 (BLR1), CXCL13 plays a central role in homeostatic trafficking of antigen-naïve B cells into and within follicles of secondary lymphoid organs. This homing of naïve B cells is necessary for antigen exposure and activation within the germinal center reaction, which are essential elements in the development and structure of secondary lymphoid organs and in the differentiation of B cells into antibody-producing plasma cells (20–22). Aberrant CXCL13 and CXCR5 expression occurs in the setting of HIV infection preventing normal B cell migration, development, and differentiation (22–26). There is also growing evidence that CXCL13 and CXCR5 participate in the pathogenesis of multiple subtypes of lymphomas (27–37). Here we test three hypotheses in studies nested within a large prospective cohort of HIV-infected men: 1) elevated CXCL13 serum levels are associated with increased AIDS-NHL risk, 2) variation in the genes coding for CXCL13 and CXCR5 is associated with increased or decreased AIDS-NHL risk, and 3) variation in the genes coding for CXCL13 and CXCR 5 is associated with CXCL13 serum levels.
The Multicenter AIDS Cohort Study (MACS), established in 1983, includes 6,972 men who have sex with men from four metropolitan areas recruited between 1984 and 2003 (38–41). Participants are re-contacted semiannually for an in-person interview, physical exam, and specimen collection. Antiretroviral use at each visit is summarized according to the Department of Health and Human Services/Kaiser Panel to define HAART usage (42). HIV seropositivity and plasma load, and CD4+ T cell counts, are measured at nearly all study visits, and sera and cell pellets are collected and stored in central repositories (43). Protocols and questionnaires have been approved by the Institutional Review Board from each center and all participants provided written informed consent.
AIDS-NHL is ascertained in the MACS through self-report with confirmation by pathology records and state cancer registries, or is identified at autopsy. A detailed description of the design and participant selection for the serum marker association study is provided elsewhere (12). Briefly, cases included all MACS participants with AIDS-NHL diagnosed prior to April, 2003 with at least one available pre-diagnosis serum specimen (n=179). One HIV-infected control was selected per case and matched on duration of HIV infection based on known date of infection (21 cases), or date of entry (± 1 year) into MACS for men seroprevalent at enrollment (158 cases), and sample availability at equivalent time-points. Using prospectively collected and stored serum, up to three serum samples were obtained for each participant: that which were collected >3 years pre-NHL (n=147), 1–3 years pre-NHL (n=148), and 0–1 year pre-NHL (n=98) in cases and follow-up matched time-points in controls (n=145, n=145, and n=103, respectively).
Cases in the genetic association study included all HIV-infected MACS participants with AIDS-NHL diagnosed prior to July, 2010 with available cells for DNA extraction (n=183). Among the 183 selected cases, 172 (94%) overlapped with cases included in the serum study, while 11 and 7 cases were unique to the genetic and serum studies, respectively (Figure 1). These differences in case numbers are due to lack of available specimen (DNA or serum) and identification of new cases between 2003 and 2010. Controls in the genetic study were selected from all HIV-infected MACS participants with available specimens and matched individually to each case in up to a 3:1 ratio on recruitment period, race, HIV positive follow-up time (± 1 year), and CD4+ T cell count (in categories of 0–49, 50–99, 100–199, 200–349, 350–499, and 500+) at last measurement pre-AIDS-NHL diagnosis for cases and follow-up matched time-point in controls. Among the 533 selected controls, 50 (9%) were included as controls in the serum study.
To investigate the association between genetic variation in CXCL13 and CXCR5 and CXCL13 serum levels, we included all the cases that were included in both the genetic study and serum study (n=172). Since the overlap of controls was low (n=50), we genotyped additional DNA specimens from serum study controls. In total, there were 172 cases and 174 controls with both CXCL13 serum levels and CXCL13/CXCR5 genotyping results included in this analysis.
Archival tissue blocks from biopsy, surgery, or autopsy were obtained for 99 AIDS-NHL cases. EBV detection was performed using EBER in situ hybridization or LMP1 immunohistochemistry in 87 cases, as described previously (44). A sample was considered positive if positive for EBER or LMP1.
Serum levels of CXCL13 were measured in sera by ELISA (R&D Systems, Minneapolis, MN) using an automated plate washer and VersaMax microplate reader and software (Molecular Devices, Sunnyvale, CA), according to manufacturer’s protocol. No samples were below the lower detection limit of 7.8 pg/ml. The interassay coefficient of variation was 8%. Data for additional cytokines and immune markers (sCD30, sCD27, IL6, TNFα, IL10, neopterin, and IP10) were included in multivariate regression models. Methods for measurement of these markers and their association with AIDS-NHL risk have been described previously (12, 45).
All SNPs within coding and non-coding regions of CXCL13 and CXCR5 plus 10 kilobases of flanking sequence on each end of each gene were identified using the European descent genotype data from the International HAPMAP Project (46). We considered only SNPs with a minor allele frequency (MAF) ≥5% which included 40 out of 50 SNPs in CXCL13 and 18 out of 20 SNPs in CXCR5. A pairwise r2≥0.80 was used to delineate groups of highly correlated SNPs within each gene, and one SNP (i.e. tagSNP) per group was selected for genotyping (47). Priority was given to potential tagSNPs with high Illumina design scores and those located in exons or other putative functional regions. Genomic DNA was extracted from cell pellets using the QIAmp DNA Blood Mini Kit according to manufacturer’s protocol (Qiagen Inc., Valencia, CA). Genotyping was performed using a customized GoldenGate® assay according to manufacturer’s protocol (Illumina Inc., San Diego, CA). Two to three positive controls (samples known to be heterozygous or homozygous for each allele based on sequencing) and negative controls (wells containing no DNA) were included in each reaction plate. Quality control (QC) replicate DNA aliquots for 5% of study subjects were distributed throughout the reaction plates. Analysis of these 44 replicate pairs revealed a high concordance proportion of 99.97%. Laboratory personnel were blinded to all research information about the samples, including the identities of the QC replicate aliquots.
One CXCR5 tagSNP (rs676925) failed to be genotyped. All other tagSNP genotypes were tested for consistency with Hardy-Weinberg equilibrium (HWE) in the control sample. Two CXCL13 tagSNPs (rs17002760 and rs171388) had HWE p-values <0.001 and were excluded from further analysis. All other tagSNPs had an HWE p-value >0.01, MAF ≥5% in our study population, and call rates greater than 93% (average call rate=98.6%). All pairwise r2 values between tagSNPs were <0.60.
Means and frequencies were calculated for select covariates separately for cases and controls. CD4+ T cell slope was defined as the change in number of CD4+ T cells per year and was calculated using pre-HAART CD4+ T cell counts for men with at least two measurements. Only CD4+ T cell counts starting with the third HIV positive study visit were included in the calculation. For participants who seroconverted during follow-up, HIV RNA at set point refers to the average HIV RNA measurement 12–24.5 months after seroconversion such that the interim slope approximated zero. For seroprevalent men recruited in 1984, HIV at set point was estimated as the measurement obtained at study visit 3 or 4. Unmatched comparisons of mean natural log transformed (loge) CXCL13 levels at each time-point between cases and controls was performed using t-tests, reported as back-transformed geometric means. Spearman’s correlation coefficients (ρ) were calculated for all pairwise associations between biomarkers at each sampling time-point for cases and controls separately.
Odds ratios (ORs) and 95% confidence intervals (95% CIs) for the association between loge CXCL13 (continuous variable) at each time-point and AIDS-NHL risk were computed using random effects multivariate logistic regression models, and included 179 cases and 179 controls. The ORs represent risk of AIDS-NHL associated with one unit increase in loge CXCL13. The matching by design between each case/control pair was incorporated into the models by adding a random effect term allowing for correlation within the sets and assuming independence between sets. Covariates were included in the models if they were strong predictors of AIDS-NHL risk and plausibly related to HIV disease progression or CXCL13 levels. These included age at AIDS-NHL diagnosis in cases or reference date in controls, absolute CD4+ T cell counts from each visit where CXCL13 was measured, loge HIV RNA levels at set point, having an AIDS diagnosis prior to AIDS-NHL diagnosis or reference date, and having been treated with HAART prior to AIDS-NHL diagnosis or reference date. Reference date in controls was determined by duration of HIV-positive follow-up time of the case in the matched set. Multiple imputation was used to estimate missing covariate data (48). Additionally, multivariate models were estimated including loge CXCL13 levels, covariates, and seven other serum biomarkers related to B cell activation, which were previously found to be strongly associated with AIDS-NHL (sCD30, sCD27, IL6, TNFα, IL10, neopterin, and IP10) (12, 45).
To estimate the association between genotypes of CXCL13 and CXCR5 tagSNPs and AIDS-NHL risk, we computed multivariate ORs and 95% CIs using random effects multivariate logistic regression models, which included 183 cases and 533 controls, and were adjusted for age at AIDS-NHL diagnosis or reference date, absolute CD4+ T cell counts from the visit closest and prior to AIDS-NHL diagnosis date or reference date, loge HIV viral RNA levels at set point, having an AIDS diagnosis prior to AIDS-NHL diagnosis or reference date, having exposure to HAART prior to AIDS-NHL diagnosis or reference date, and race/ethnicity.
The association between genotypes of tagSNPs and mean loge CXCL13 levels was modeled using repeated measures (mixed) linear regression to estimate mean ratios (MR) and 95% CI, and included 172 cases and 174 controls and the following covariates as fixed effects: age at AIDS-NHL diagnosis or reference date, absolute CD4+ T cell counts from each visit where CXCL13 was measured, loge HIV viral RNA levels at set point, having an AIDS diagnosis prior to AIDS-NHL diagnosis or reference date, having exposure to HAART prior to AIDS-NHL diagnosis or reference date, and race/ethnicity. Two random effect variables were included, one to account for the non-independence of observations on the same subject across the three visits and a second to account for the non-independence between each case and its matched control. The variance of loge CXCL13 explained by each tagSNP was estimated by subtracting the r2 value for a model that included the parameters for each tagSNP from the r2 for a model that excluded the tagSNP.
Cases and controls were of similar age, and the majority were White, non-Hispanic (Table 1). Due to differences in control matching between the serum marker study and the genetic association study, there were larger differences between cases and controls from the serum study compared to the genetic study for several variables related to HIV disease progression, including median HIV RNA at setpoint, having a prior AIDS illness, CD4+ T cell slope, and CD4+ T cell count. Less than 10% of the study population was treated with HAART. The majority of tumors were systemic (68%), most of which were DLBCL. Fewer than 50% of cases had adequate tumor tissue available for EBV testing. A majority of tested specimens were EBV positive (68%), and a larger proportion of PCNSL tumors (89%) were EBV-positive compared to systemic tumors (58%).
Mean CXCL13 levels were consistent among the three time-points in the controls (74, 80, and 83 pg/mL, Figure 2). Mean CXCL13 levels were significantly higher in cases compared to controls at each time-point (p<0.001) and increased as time approached AIDS-NHL diagnosis date (123, 153, 178 pg/mL). In the adjusted models, a one unit increase in loge CXCL13 levels was significantly associated with AIDS-NHL risk at all three time-points (>3 years: OR=3.24, 95% CI=1.90–5.54; 1–3 years: OR=3.39, 95% CI=1.94–5.94; 0–1 year: OR=3.94, 95% CI=1.98–7.81, Table 2). In stratified analyses, the association between CXCL13 levels and AIDS-NHL appeared to be stronger for systemic lymphoma compared to PCNSL at the >3 years and 1–3 years time-points, although the differences in ORs were not statistically significant. At the time-point furthest from AIDS-NHL diagnosis, the association between CXCL13 and AIDS-NHL was stronger for EBV-negative compared to EBV-positive AIDS-NHLs. As time approached NHL diagnosis date, higher ORs were observed for EBV-positive tumors, although the differences in ORs by EBV status were not statistically significant.
In multivariate models that included CXCL13 levels and seven other serum biomarkers, the highest statistically significant OR was observed for CXCL13 (OR=2.68, 95% CI=1.38–5.26) at the >3 years time-point (Table 3). At the 1–3 year time-point, none of the biomarkers were independently significantly associated with AIDS-NHL, however in strength of association CXCL13 was second highest following neopterin. At the 0–1 year time-point, neopterin was the most strongly associated with AIDS-NHL risk (OR=5.18, 95% CI=1.36–19.8). Many of these biomarkers were positively correlated, however no two exhibited correlation greater than ρ=0.72, with most between ρ=0.20 and ρ=0.50 (Supplementary Table 1).
AIDS-NHL risk was inversely associated with carriership of the minor allele of CXCL13 rs355689; OR=0.66, 95%, CI=0.46–0.97 for heterozygotes and OR=0.87, 95% CI=0.5–1.55 for homozygotes, compared to major allele homozygotes (Table 4). CXCL13 levels were inversely associated with carriership of the minor allele of rs355689; MR=0.82, 95% CI=0.68–0.99 for minor allele homozygotes compared to major allele homozygotes (Figure 3, Supplementary Table 2). Two additional SNPs were associated with CXCL13 levels, CXCL13 rs17002733 and CXCR5 rs630923 (Figure 3), although these risk estimates were based on only a few minor allele homozygote study subjects. The percentage of variation in loge CXCL13 levels explained by rs355689, rs17002733, and rs630923 was 2.3%, 4.2%, and 3.2%, respectively. CXCL13 levels were strongly associated with AIDS-NHL risk at all three time-points within each genotype strata for rs355689, rs1700273, and rs630923 (Supplementary Table 3).
There is accumulating evidence that CXCL13 and CXCR5 are involved in the pathogenesis of lymphoma in general, and in particular, of extranodal lymphomas which comprise the majority of AIDS-NHL (3). Here, we show that CXCL13 serum levels predict B cell non-Hodgkin lymphoma risk in a cohort of HIV-infected men. The consistent, significant elevation of CXCL13 at three separate time-points across three or more years prior to the diagnosis of AIDS-NHL, and the significant association between CXCL13 and AIDS-NHL after adjustment for AIDS-NHL risk factors and other biomarkers, suggests that CXCL13 may be one of the drivers of the chronic B cell activation phenotype observed in AIDS-NHL.
In a healthy system, CXCL13 is constitutively expressed by T helper cells in lymphoid follicles where it regulates B cell homing (20). CXCL13 has been detected in follicular lymphomas and extranodal lymphomas occurring in a variety of organs including the brain, eyes, stomach, and skin (29, 32)(30)(33)(34). The source of CXCL13 in these lymphomas appears to vary by tumor site and includes malignant B cells, neighboring cells in the tumor environment, and cells produced elsewhere in the body. Malignant B cells of extranodal lymphomas express high levels of CXCR5, (27, 28, 30, 32, 35), the CXCL13 receptor which is normally only expressed in lymphoid follicles (49).
There are at least three possible explanations for our observation that CXCL13 is elevated prior to AIDS-NHL. One possibility is that CXCL13 and CXCR5 drive the formation of ectopic germinal centers, permitting chronic B cell activation and the accumulation of mutations and translocations which contribute to the pathogenesis of lymphoma (50). Alternatively, CXCL13 may be responsible for attracting CXCR5-expressing malignant B cells to these extranodal sites leading to the formation of tumors (27, 37). Lastly, it is possible that malignant B cells from the developing tumors, and/or tumor-infiltrating macrophages and/or T cells, are secreting CXCL13, and that CXCL13 is not contributing to the etiology of these tumors, but may be promoting their growth. However, the strong and consistent association between CXCL13 and AIDS-NHL risk for over three years suggests an etiological role of CXCL13 in AIDS-NHL.
Despite the accumulating in vitro evidence for CXCL13 and CXCR5 in the pathogenesis of lymphoma, there is only one prior epidemiologic study on CXCL13 and lymphoma, in HIV-infected or uninfected populations. Widney et. al. reported that CXCL13 serum levels are elevated up to 2.5 years (mean of 8.2 months) prior to cancer diagnosis in a small sample (n=46) of AIDS-NHL cases when compared to similarly sized groups of either AIDS controls without cancer, or HIV-infected controls (27). Our study confirms and extends this observation utilizing a larger sample of AIDS-NHL cases and HIV-infected controls sampled at three time-points prior to cancer diagnosis.
In stratified models, CXCL13 was significantly associated with systemic lymphoma but not PCNSL, and the point estimates were also higher for systemic lymphoma. These findings are consistent with what has been observed for other B cell activation associated biomarkers and AIDS-NHL (12). Among those cases with known EBV tumor status (n=87, 49%), CXCL13 was more strongly associated with EBV-negative tumors compared to EBV positive tumors at the time-point furthest from NHL diagnosis. In our study, 58% of systemic tumors and 89% of PCNSL tumors were EBV-positive. These findings are consistent with the hypothesis that systemic AIDS-NHLs are more likely to arise from the downstream effects of chronic B cell activation occurring over several years while PCNSL tumors arise principally due to the oncogenic properties of uncontrolled EBV reactivation in severely immunocompromised individuals (3, 15).
Mean CXCL13 levels in cases, but not controls, increased with time to diagnosis, while the ORs remained consistent. Previous studies have shown that CXCL13 serum levels are elevated in HIV-infected patients compared to healthy controls (24), and increase with HIV disease progression (25). Thus, our observation regarding increasing CXCL13 levels with time to diagnosis may reflect a higher rate of HIV disease progression in cases, which is accounted for in the CD4+ T cell count- and HIV RNA level–adjusted OR estimates. A previous study showed that B cell migration towards CXCL13 was higher in HIV-infected patients with low CD4+ T cell counts (<350 cells/μL) compared to patients with higher cell counts (>350 cells/μL) or healthy controls, suggesting an increased functional effect of CXCL13 in advanced HIV disease (24). Alternatively, the developing tumor may be responsible for the increasing CXCL13 levels with time to diagnosis in cases. However, the strong and consistent association between CXCL13 and AIDS-NHL risk for over three years argues against this possibility.
The seven other biomarkers (sCD30, sCD27, IL6, TNFα, IL10, neopterin, and IP10) included in the multivariate biomarker model were selected from over 30 biomarkers with strong (adjusted ORs~2.0 or greater), statistically significant and consistent associations with AIDS-NHL risk at multiple time-points (12, 45). Though mostly positively correlated, these biomarkers were not strongly correlated, which likely reflects their diverse biological functions. All of these biomarkers are thought to induce or result from B cell immune activation, yet they are secreted from, expressed on, and bind to receptors on various cell subsets including Th2 cells, Th17 cells, follicular helper T cells, and monocytes/macrophages (Mo/MΦ).
Th17 cells are a recently-identified subset of pro-inflammatory CD4+ T cells, distinct from Th1 or Th2 cells, which have an established pathogenic role in autoimmune disease and an emerging role in cancer initiation, particularly in the setting of infection (51). Th17 cells express the retinoic acid-related orphan receptor, (ROR)γt, which promotes the gene expression pattern that characterizes the Th17 phenotype (52). There is increasing evidence in support of the hypothesis that Th17 cells drive chronic B cell immune activation and lymphomagenesis and our findings add support to this hypothesis (45, 53). CXCL13, IL6, and TNFα are key examples of Th17-associated chemokines/cytokines (54), and these biomarkers exhibited the strongest association with AIDS-NHL in our study at the >3 year pre-NHL time-point in the multivariate biomarker model.
It has also recently been observed that microbial translocation causes systemic immune activation in HIV and is associated with HIV progression (55–57). Neopterin, which is produced exclusively by macrophages, was the biomarker most strongly associated with AIDS-NHL risk at the visit closest to AIDS-NHL diagnosis date, suggesting that Mo/MΦ activation, possibly due to microbial translocation, coincides with the presentation of AIDS-NHL, but may not necessarily be driving the process.
Our study population was composed nearly entirely of men who were unexposed to HAART prior to AIDS-NHL diagnosis or reference date. HAART, which results in HIV RNA suppression and partial restoration of immune competence, may alter the association between CXCL13 and AIDS-NHL. In recent work, serum CXCL13 levels were significantly reduced in HIV-infected men following HAART initiation compared to their pre-HAART levels, yet still remained higher in comparison to HIV-uninfected controls (58). Thus, this incomplete normalization of CXCL13 by HAART suggests that there is an ongoing environment of chronic B cell activation that could predispose to AIDS-NHL for HAART users and non-users alike.
We did not find strong evidence for an association between genetic variation in either CXCL13 or CXCR5 and AIDS-NHL risk. This does not rule out the possibility that rare variants, structural variation, or other SNPs in these genes not captured by our tagSNPs, could be associated with AIDS-NHL. Interestingly, the minor allele of rs355689 was associated with both a reduced risk of AIDS-NHL and also a lower level of CXCL13. Rs355689 is an intronic SNP and is not in LD with any other known common SNPs in CXCL13 or in neighboring genes according to HAPMAP (59), nor is it predicted to be located in a splicing domain (60). This SNP may be marking some unknown functional variation in CXCL13, or these findings may be due to chance. Furthermore, rs355689 does not appear to modify the association between CXCL13 levels and AIDS-NHL risk (Supplementary Table 3).
Polymorphisms in chemokines and their receptors are important in HIV risk and progression, as well as in the development of AIDS-NHL. The delta 32 deletion polymorphism of chemokine receptor CCR5, protective against HIV-1, is associated with decreased risk of AIDS-NHL, while the stromal cell-derived factor 1 chemokine polymorphism (SDF1-3′A), associated with variable progression rates of HIV, is associated with increased risk of AIDS-NHL (61, 62). A recent study examined the association between four CXCR5 SNPs and NHL in an HIV-uninfected Chinese population (63). Three of the SNPs (rs78440425, rs148351692, and rs80202369) have minor allele frequencies less than 1% in Caucasians and rs6421571 was greater than 10kb upstream of the CXCR5 transcription start site, thus none of these SNPs were genotyped or represented by tagSNPs in our study. Additionally, rs6421571 was not in LD with any of the tagSNPs in our study. The authors reported a significant increased risk of NHL associated with the minor allele of rs80202369 and rs6421571 which raises the possibility that there may be functional variation in or near CXCR5 with consequences for NHL development.
Two other SNPs were significantly associated with CXCL13 levels. Minor allele homozygote carriers for CXCL13 rs17002733 had significantly reduced levels of CXCL13 compared to major allele homozygote carriers. Similar to rs355689, rs17002733 is an intronic SNP, which is neither in LD with any other known SNPs nor has any predicted functional consequence (59, 60). Minor allele homozygote carriers for CXCR5 rs630923 had significantly increased levels (greater than 2-fold) of CXCL13 compared to major allele homozygote carriers. Interestingly, rs630923 is located in the 5′ flanking region and is predicted to be involved in transcriptional regulation by affecting transcription factor binding sites (60). Homozygotes for the rs630923 minor allele may have decreased cell surface expression of CXCR5 leading to an excess of CXCL13 serum levels. Disruption of B cell homing into follicles could cause a reduction in B cell activation. The near 4-fold reduction of AIDS-NHL risk observed for rs630923 minor allele homozygote carriers (non-statistically significant) supports this theory and our hypothesis that AIDS-NHL develops from chronic B cell hyperactivation. However, given that the minor allele homozygote carriers represent a small number of cases and controls in our study, and the confidence intervals surrounding the MR and OR for this SNP are very wide, we cannot rule out the possibility of chance. The percent of variation in CXCL13 explained by our SNPs is small, particularly when compared to the large percent variation explained by variables related to HIV disease progression: 4.2% for rs17002733 compared to 31.7% for CD4+ T cell count, for example. This suggests that individual variation in other risk factors may be driving the differential expression in CXCL13 between cases and controls.
The main strength of this study is the inclusion of a large sample of known AIDS-NHL cases with detailed longitudinal covariate data and specimen availability. In addition to CXCL13, this study population has been characterized with respect to a number of important biomarkers related to AIDS-NHL, permitting us to evaluate our data in this context. Although we controlled for HIV progression with HIV RNA levels and CD4+ T cell numbers, our serum marker ORs could be affected by residual confounding by HIV progression. However, in preliminary analyses, we also included CD4+ T cell slope in our regression models and did not observe any appreciable difference in the ORs (64). For our genetic association study, the inclusion of tagSNPs allowed us to capture all known common variation in our candidate genes. However, due to the sample size, these analyses were underpowered to detect weak or moderate genetic effects. Complementing our tagSNP-AIDS-NHL association testing was the evaluation of potential functional consequence of genetic variation on CXCL13 levels.
As HIV-infected individuals are living longer and to older ages, chronic B cell activation will continue to increase risk of AIDS-NHL. We have identified CXCL13 as an important and novel early pre-diagnosis biomarker for AIDS-NHL that potentially may be used in risk assessment, early detection, or prevention of AIDS-NHL. Our study also raises the possibility that genetic variation in CXCL13 and CXCR5 could influence CXCL13 levels and AIDS-NHL risk, although these findings require validation due to the small sample size.
Dr. Jeffrey Gornbein (UCLA) provided consultation on the statistical analyses. Mr. Joe DeYoung at the UCLA Neuroscience Genomics Core genotyped the CXCL13/CXCR5 SNPs. The authors would like to thank the MACS participants, without whom this and many other studies would not be possible.
This work was supported by grants from the National Institutes of Health (K07-CA-140360, R25-CA-87949, an NCI supplement to U01-AI-035040, a UCLA CFAR grant P30-AI-028697, R01-CA-057152, R01-CA-121195, and the UCLA Older Americans Independence Center grant P30-AG-02874). The Multicenter AIDS Cohort Study (MACS) is funded by the National Institute of Allergy and Infectious Diseases, with additional supplemental funding from the National Cancer Institute (UO1-AI-35042, UL1-RR025005 [GCRC], UO1-AI-35043, UO1-AI-35039, UO1-AI-35040, UO1-AI-35041). Samples and data in this manuscript were collected by the MACS, with centers (Principal Investigators) at The Johns Hopkins Bloomberg School of Public Health (Joseph B. Margolick, Lisa P. Jacobson), Howard Brown Health Center, Feinberg School of Medicine, Northwestern University, and Cook County Bureau of Health Services (John P. Phair, Steven M. Wolinsky), University of California, Los Angeles (Roger Detels, Otoniel Martínez-Maza), and University of Pittsburgh (Charles R. Rinaldo). Some of this work was carried out in the facilities of the UCLA AIDS Institute, which were supported, in part, by funds from the James B. Pendleton Charitable Trust and the McCarthy Family Foundation.
Conflicts of Interest: None of the authors have a conflict of interest to report.
AUTHORSHIPContribution: S.K.H. and O.M.M. were responsible for the conception and design, acquisition of data, analysis and interpretation of data, and writing and revision of the manuscript. W.Z., S.C.C., L.M., D.W., and E.V. were responsible for generation of laboratory data, and reviewing/revising of the manuscript. E.C.B. was responsible for the conception and design and reviewing/revising of the manuscript. D.C. and M.S. participated in data analysis and reviewing/revising of the manuscript. L.P.J., J.H.B., S.W., C.R.R., R.F.A., R.D., and Z.F.Z. were responsible for the conception and design, acquisition of data, and reviewing/revising the manuscript.