Association of HIV-1 Infection and Antiretroviral Therapy With Type 2 Diabetes in the Hispanic Population of the Rio Grande Valley, Texas, USA
Juan Carlos Lopez-Alvarenga1† 
Dora A. Martinez2† 
Alvaro Diaz-Badillo1 
Liza D. Morales1 
Rector Arya1 
Christopher P. Jenkinson1 
Joanne E. Curran1 
Donna M. Lehman3 
John Blangero1 
Ravindranath Duggirala1 
Srinivas Mummidi1*‡ 
Ruben D. Martinez2*‡- 1Department of Human Genetics, South Texas Diabetes and Obesity Institute, School of Medicine, University of Texas Rio Grande Valley, Edinburg, TX, United States
- 2Valley AIDS Council, Harlingen, TX, United States
- 3Department of Medicine, University of Texas Health San Antonio, San Antonio, TX, United States
The Rio Grande Valley (RGV) in South Texas has one of the highest prevalence of obesity and type 2 diabetes (T2D) in the United States (US). We report for the first time the T2D prevalence in persons with HIV (PWH) in the RGV and the interrelationship between T2D, cardiometabolic risk factors, HIV-related indices, and antiretroviral therapies (ART). The PWH in this study received medical care at Valley AIDS Council (VAC) clinic sites located in Harlingen and McAllen, Texas. Henceforth, this cohort will be referred to as Valley AIDS Council Cohort (VACC). Cross-sectional analyses were conducted using retrospective data obtained from 1,827 registries. It included demographic and anthropometric variables, cardiometabolic traits, and HIV-related virological and immunological indices. For descriptive statistics, we used mean values of the quantitative variables from unbalanced visits across 20 months. Robust regression methods were used to determine the associations. For comparisons, we used cardiometabolic trait data obtained from HIV-uninfected San Antonio Mexican American Family Studies (SAMAFS; N = 2,498), and the Mexican American population in the National Health and Nutrition Examination Survey (HHANES; N = 5,989). The prevalence of T2D in VACC was 51% compared to 27% in SAMAFS and 19% in HHANES, respectively. The PWH with T2D in VACC were younger (4.7 years) and had lower BMI (BMI 2.43 units less) when compared to SAMAFS individuals. In contrast, VACC individuals had increased blood pressure and dyslipidemia. The increased T2D prevalence in VACC was independent of BMI. Within the VACC, ART was associated with viral load and CD4+ T cell counts but not with metabolic dysfunction. Notably, we found that individuals with any INSTI combination had higher T2D risk: OR 2.08 (95%CI 1.67, 2.6; p < 0.001). In summary, our results suggest that VACC individuals may develop T2D at younger ages independent of obesity. The high burden of T2D in these individuals necessitates rigorously designed longitudinal studies to draw potential causal inferences and develop better treatment regimens.
Introduction
Current antiretroviral therapies (ART) have dramatically improved the quality of life and increased the life span of persons with HIV (PWH). However, when compared to the general population, PWH have an increased incidence of common age-associated degenerative diseases, including metabolic disorders such as type 2 diabetes (T2D) and cardiovascular disease (CVD). The Hispanic/Latino community comprises ~23% of the PWH and 27% of the ~36,400 newly diagnosed HIV cases in the US in 2018. A recent study estimated that the median Latino disparity in the HIV prevalence to be 2.4 (1). Hispanics also have a higher prevalence of cardiometabolic risk factors that predispose them to T2D when compared to non-Hispanic whites (NHW). For example, Lorenzo et al. compared the non-HIV infected Mexican American population with non-Hispanic whites, using Matsuda and insulinogenic indexes to assess insulin resistance. They showed that Mexican Americans have an odds ratio of 2.3 (95%CI 1.5, 3.34) for insulin resistance (2). Thus, HIV-1 infected Hispanics are at an even greater risk of developing T2D and T2D-related complications. However, few studies have systematically examined the incidence of these metabolic disorders in HIV-infected Hispanics, especially in the U.S.-Mexico border region who encounter unique socio-cultural, economic, migration-related challenges.
Prior to the introduction of the new generation of ART, the incidence of diabetes was found to be 4-fold higher in HIV-positive men when compared to HIV-negative men in the Multicenter AIDS Cohort Study (MACS) (3). While the recent ART regimens have less dramatic effects, there appears to be growing evidence of increasing T2D prevalence in HIV-positive individuals in the US. In a recent large cross-sectional study (N = 17,946,580) that included 90,900 persons with HIV (60.8% male), Birabaharan et al. found that the overall prevalence of T2D in the HIV group was 22.1% (N = 20,080) which is much higher than the T2D prevalence in the US population which is 14.9%. Notably, the prevalence of T2D was higher in the group classified as “others,” which includes the Hispanic population (27.6 vs. 21.3% in the white population). Nansseu et al. performed a systematic review and metanalysis of 44 studies conducted across the globe and projected that in a “cohort of 1,000 PWH starting ART, about 50 may have T2D after 3.7 years of follow-up while around 457 may develop prediabetes within the same period” (4). Other metabolic abnormalities that predispose HIV-infected individuals to T2D have also been reported with currently used ART regimens. A recent large-scale study (N = 14,084) showed that 22% of lean PWH individuals became overweight just after 3 years on ART (5). A similar trend was observed in overweight people who became obese. In another study, Sax et al. performed a pooled analysis of weight gain in eight randomized controlled clinical trials of treatment-naïve PWH (N = ~5,000) and found that demographic factors (e.g., female, black), HIV related factors at the start of the treatment (e.g., low CD4, high viral load), and ART regimen used may play a role in weight gain (6). Thus, the contributions to weight gain in ART-treated PWH could be multifactorial, and both HIV/ART and traditional risk factors could lead to increased weight gain and metabolic risk in PWH. This is of great concern as a majority of MAs are overweight or obese with increased genetic susceptibility to T2D and other metabolic diseases and, when infected with HIV, may have a poorer response to T2D medications compared to non-Hispanic whites (7). The increased risk to T2D could be attributed to a family history of diabetes, impaired glucose tolerance, obesity, insulin resistance, and genetic factors.
The Rio Grande Valley (RGV) of South Texas has one of the highest prevalence rates of obesity and T2D in the US (8, 9). Thus, the intersection of HIV-1 infection and metabolic disorders such as T2D impose enormous health burden in the PWH in the RGV. However, there is a lack of focused studies in the PWH in the RGV especially with regard to complications arising due to ART and their management. This is important as the Hispanics in the RGV, a majority of whom are Mexican Americans (10), face unique socio-economic challenges and cultural practices and have a genetic predisposition to develop metabolic diseases. Here, we report for the first time the prevalence of T2D and its associated risk factors in PWH in the RGV, who are all enrolled at a single clinic, and contrast them with the HIV-negative Mexican Americans in the US and South Texas. We also report the interrelationship between cardiometabolic risk factors, T2D risk, HIV viral load, various immunological parameters, and treatment regimens in this cohort.
Methods
The Valley AIDS Council Cohort
All the individuals received medical care under the supervision of two physicians at the two Valley AIDS Council (VAC) clinic sites located in Harlingen and McAllen, TX. The RGV mainly comprises of four counties (Cameron, Hidalgo, Starr, and Willacy) that are in the Southeastern region of Texas bordering Mexico and the overall percentage of Hispanics is ~90% (10). Henceforth, this cohort of PWH in the RGV will be referred to as the Valley AIDS Council Cohort (VACC). The VAC is the only Ryan White funded agency providing medical care and supportive services for people living with HIV/AIDS in regions of Texas that are south and east of San Antonio, TX, including the cities of Corpus Christi, McAllen, Harlingen and Brownsville. In 1994, a Ryan White Title III grant enabled the agency to provide primary health care for HIV/AIDS patients in this region, as a primary healthcare outpatient clinic. The VAC has provided continuous services, including HIV testing and prevention efforts, since its start in 1987 and currently serves a large patient population (~2,000).
Data Collection
Valley AIDS Council Cohort
All analyses were conducted using deidentified data extracted from a retrospective database using protocols approved by the Institutional Review Board (IRB) of the University of Texas Rio Grande Valley (UTRGV). The VAC houses four servers at the Westbrook Clinic in Harlingen that work in tandem to run eClinicalWorks (Electronic Medical Records). The VAC also uses ARIES (AIDS Regional Information and Evaluation System) and CMBHS (Clinical Management for Behavioral Health Services), both of which are State databases where patient information is entered as part of the clinical routine. A total of 1,827 records with unique IDs were abstracted from medical charts which were documented between November 2017–October 2019. The following variables were extracted from the database: demographic (age, sex, and marital status), anthropometric [height, weight, and body mass index (BMI)] and cardiometabolic traits (blood pressure, fasting glucose, HbA1c, and lipids), and HIV-related virological and immunological indices (viral load, absolute and relative CD4+ and CD8+ T cell counts, and CD4/CD8 ratio). For analysis, the WHO obesity classification system based on BMI (computed with kg/m2 and categorized as normal weight <25, overweight between 25 and 30, obesity class I between 30 and 35, class II between 35 and 40 and class III >40) was used (11). We used the ADA definition for defining T2D, which includes one of the following conditions: a clinical history of T2D, a fasting glucose ≥ 126 mg/dL, or HbA1C > 6.4% (12).
San Antonio Mexican American Family Studies and NHANES 2018 Data
For comparison with the VACC, we used data from known HIV-1 negative cohorts with Mexican American ancestry, namely the San Antonio Mexican American Family Studies (SAMAFS) and a subset of National Health and Nutrition Examination Survey (NHANES) individuals with reported Mexican American ancestry, henceforth referred to as the Hispanic Health and Nutrition Examination Survey (HHANES). The SAMAFS refers to the combination of two San Antonio-based family studies called the San Antonio Family Heart Study (SAFHS) (13, 14) and the San Antonio Family Diabetes and Gallbladder Study (SAFDGS) (15, 16). The SAFHS included 1,431 individuals in 42 large pedigrees at baseline from low-income Mexican Americans selected at random without regard to the presence or absence of disease. The SAFDGS enrolled 579 individuals distributed across 32 large pedigrees at baseline. These families were ascertained by probands with T2D who were low-income Mexican American men and women. Both study-specific protocols were approved by the IRBs of the University of Texas Health San Antonio (UTHSA) and the UTRGV. The major objective of SAMAFS was to investigate the genetics of complex diseases such as T2D, CVD, obesity and their related traits such as BMI, glucose, insulin, lipids, and blood pressure measures using standard protocols (14, 17–20). The T2D status was diagnosed according to the criteria of the World Health Organization: fasting glucose levels ≥ 126 mg/dL and/or 2-h glucose levels ≥ 200 mg/dL (21) and/or the American Diabetes Association criteria: fasting glucose levels ≥ 126 mg/dL (22). Participants who did not meet these criteria but who reported that they were under treatment with either oral antidiabetic agents or insulin and who gave a history of diabetes were also considered to have T2D. The SAMAFS participants have been followed in a mixed longitudinal fashion, up to a maximum of five visits. The original families have been expanded through additional recruitment and eight new families were added (i.e., SAFDGS) over the years. The data used for this study correspond to those measured at the last clinic examination, and individuals with T2D status information and associated cardiometabolic data were considered for the analysis (N = 2,498). The NHANES data were derived from the surveys conducted during 2017–2018. The HHANES data included measured glucose, HbA1c, blood pressure, triglycerides, and HDL-C for comparison with VACC. The T2D status for the individuals enrolled in NHANES was obtained from the publicly available databases (23).
Statistical Analyses
We used 1,827 entries from VACC for cross-sectional data analyses. We used mean values of the quantitative variables from unbalanced visits across 20 months. Continuous variables were aggregated to obtain mean values and if necessary, were transformed to obtain normal distributions. An illustrative example of square root transformation to normalize CD4 counts is shown in Supplementary Figure 1. Descriptive statistics were used to evaluate the demographic and biochemical variables related to T2D and cardiometabolic risk factors in the VACC. Student t-tests, adjusted for variance, were performed to contrast PWH with and without T2D. Robust regression methods with robust Huber-White sandwich method HC3 (an approximation to a jackknife estimator using squared residuals divided by the square of 1-h, the letter “h” represents the leverage values in the “hat” matrix) were used to determine the associations between various immunological and HIV-related indices (24). The association between the treatment regimen and T2D, and various cardiometabolic risk factors, i.e.,: LDL-C higher than 110 mg/dL, systolic blood pressure > 140 mmHg, diastolic blood pressure > 90 mmHg, and HDL-C adjusted by sex (females <50 mg/dL and males <40 mg/dL) were performed using logistic regression to obtain odds ratios (ORs) (95% CI). Each model was adjusted by sex, age, WHO BMI categories, viral loads (<50 copies/mL, between 50 and 1,000 copies/mL; and >1,000 copies/mL), and absolute CD4+ counts (<500 cells/mm3, between 500 and 900 cells/mm3, and >900 cells/mm3). A backward stepwise regression with p-values <0.20 was performed. The variables collected from the VACC were compared with the cardiometabolic trait data from SAMAFS (N = 2,498) after adjusting for pedigree structure and the HHANES data (N = 5,989) using robust regressions. Descriptive statistics such as means, standard deviations, and frequencies were used to compare VACC, SAMAFS, and HHANES data. ANOVA with robust HC3 method was used for comparing the continuous variables between the three cohorts after adjusting for pedigree structure. All analyses were performed with STATA/SE 16.0 (Stata Corp. LLC. College Station TX. USA).
Results
Descriptive Statistics of the VACC
We extracted a total of 2,959 entries from the patient registries in the VAC database, out of which 1,827 entries were suitable for further analysis. All entries were from patient registries in the VAC database entered between November 2017 to October 2019. However, not all variables were available for each registry, and the sample size varied from trait to trait based on the available information. Hence, maximal information was provided for a given trait by using the available data for description. Please note that the use of a value was conditioned by the presence of other variables (for example, the descriptive tables for age require the presence of sex and age from a specific sample; meanwhile, the simple description of sex in the same sample was not conditioned by the presence of age) giving slightly different mean values from descriptive tables and the text. The median number of visits was 8 (Q1 5.3, Q3: 9.2) for each individual. The following data were available: sex (1,362/1,707; 80% male) and T2D status (917/1807; 51%). The mean (SD) age was 42.1 (12.7) years old, BMI 29.4 (6.1), HDL-C 44.5 (12.7) mg/dL, LDL-C 96 (25) mg/dL, triglycerides 176 (121.6) mg/dL. The most common treatment was integrase stand transfer inhibitors (INSTIs) in any combination was 64% (918/1,425) of the PWH in the VACC. The most frequent combination was INSTI+Nucleoside/Nucleotide reverse transcriptase inhibitors (NRTI) in 22% (317/1,425) of the individuals. Other treatment regimens used were combinations containing Protease inhibitors (PIs) and NRTI + Non-nucleoside reverse transcriptase inhibitors (NNRTI). The absolute CD4+ counts were 589 (315) cells/mm3, absolute CD8+ counts 909 (427) cells/mm3, and the mean CD4/CD8 ratio was 0.76 (0.49). The frequency of BMI categories as defined by the WHO (11) were as follows: normal 375 (23%), overweight 626 (38%), obesity class I were 393 (24%), class II were 147 (9%), and class III were 97 (6%). Males had higher blood pressure and lower BMI, HDL-C, CD4/CD8 ratio compared with females. Table 1 summarizes the descriptive conditional statistics of the demographic and clinical data and comparative analysis between the sexes in the VACC.
Table 1. General description of the Valley AIDS Council Cohort (VACC).
HIV-1 Related Indices and Their Relation to VACC Metabolic Traits
We then evaluated several HIV-1 related indices in the VACC and their relationship to various metabolic traits. Eighty percent of the VACC individuals had viral loads lower than 50 copies/mL (N = 1,182), six percent showed viral loads between 50 and 1,000 copies/mL (N = 92), and fourteen percent showed >1,000 copies/mL (N = 212). The virus load was an independent predictor of CD4+ counts, CD8+ counts, % CD4, % CD8, and CD4/CD8 ratio adjusted by age, sex, presence of T2D, lipid concentration, and duration of enrollment in the clinic (Table 2). There was no difference in viral loads between men and women (p = 0.34). Figure 1 shows a violin plot depicting the relationship between CD4/CD8 ratio and viral load, and as expected, higher viral loads were negatively associated with CD4/CD8 ratios. We also found that increasing BMI was positively associated with increasing CD4+ cell counts (Figure 2). T2D was found in 51% of the VACC individuals (917/1,807), with no differences by sex (females OR = 0.89, 95%CI: 0.68, 1.13; p = 0.306). However, for example, participants with T2D had higher CD8+ cell % (Cohen-d 25.6%, p < 0.001) and less LDL-C (Cohen-d 0.24, p < 0.001) compared to those without T2D (Table 3).
Table 2. Factors associated with laboratory immunological parameters in the VACC¶.
Figure 1. Violin plot of the relationship between CD4/CD8 ratio and viral load. Square root transformed CD4/CD8 ratios are shown on the y-axis which decrease as a function of increasing viral loads (<50 copies, 50–1,000 copies and >1,000 copies) shown on x-axis [b = −0.13 (se 0.01), p < 0.001].
Figure 2. A histogram depicting an increase of absolute CD4 lymphocytes with an increase of body weight. The regression with CD4+ T cell counts (square root transformed) was b = 0.48 (robust HC3 se = 0.16), p = 0.004. NW, normal weight; OW, overweight; OB-C1, obesity-category 1; OB-C2, obesity-category 2; and OB-C3, obesity category 3.
Table 3. Descriptive statistics of the VACC individuals stratified by T2D status¶.
Association Between ART Regimens and T2D in VACC
ART was associated with decreased viral load and increased CD4+ T cell counts but not with metabolic dysfunction. However, individuals on different ART regimens showed differences in the association with T2D (Table 4). Notably, we found that individuals on any combination regimen with INSTI were associated with increased risk to T2D [2.08 (95%CI 1.67, 2.6; p < 0.001)]. Other treatment regimens like NRTI + NNRTI and any combination regimens with PIs did not show association with T2D (Table 4). When stepwise logistic regression was used, with adjustment for sex, age, BMI, virus load, and treatments to explain the presence of cardiometabolic risk factors (dichotomous variables described in the method section), T2D was associated with viral loads >1,000 [OR: 1.56 (95%CI: 1.09, 2.22), p = 0.014], and CD4 lower than 500 [OR: 1.36 (95%CI: 0.96, 1.93), p = 0.088]. Any treatment with INSTI-containing combination was associated with LDL-C [OR 0.62 (95%CI: 0.48, 0.80), p < 0.001], systolic blood pressure [OR 0.71 (95% CI: 0.54, 0.95), p = 0.021] and, triglyceride concentration > 150 mg/dL [OR: 0.85 (95%CI 0.67, 1.07)].
Table 4. Treatments used and association with T2D in the cross-sectional analysis of VACC¶.
Comparative Analysis of Metabolic Traits in the VACC and Other MA Populations
We have limited our comparisons between cohorts to the currently available data from the VACC. Both HHANES and SAMAFS have a higher proportion of female participants when compared to VACC. The mean (SD) differences between the VACC, HHANES and SAMAFS individuals in the traits shown below are: BMI, 29.4 (6), 30.9 (6), and 31.1 (7) (p < 0.001); triglycerides, 176 (121), 137 (109), and 151 (140) mg/dL (p < 0.001); and HDL-C, 44 (13), 53 (15), 48 (14) mg/dL (p < 0.001), respectively. This post hoc analysis was done with Tukey-Kramer for pairwise comparisons, all adjusted for sex, age, and T2D. Thus, the VACC individuals had lower BMI but significantly higher triglycerides and lower HDL-C. Descriptive statistics organized by sex from SAMAFS and HHANES individuals are shown in Supplementary Tables 1, 2. Both sexes showed Cohen-d >0.3 for diastolic blood pressure and HDL-C, and about 0.2 for BMI and triglycerides, with men showing a lower effect than females. In SAMAFS, T2D showed effects greater than Cohen-d 0.5 for age, systolic blood pressure, and triglycerides; Cohen-d 0.37 for BMI and 0.24 for HDL-C. Tables 5, 6 show the descriptive statistics for these traits in both SAMAFS and HHANES individuals stratified by their T2D status.
Table 5. General description and stratification by T2D status of selected variables from the SAMAFS¶.
Table 6. General description and stratification by T2D status of selected variables from the HHANES¶.
The prevalence of T2D in VACC was 51% (i.e., 917/1,807) compared to 27% (i.e., 672/2,498) in SAMAFS and 23% (i.e., 262/1,157) in HHANES (p < 0.001). The males in the VACC had a greater risk of developing T2D, which is similar to that seen in the SAMAFS. However, the PWH with T2D in the VACC were younger (4.7 years or 0.3 standard deviations) and had lower BMI (BMI 2.43 units less or 3.7 standard deviations) when compared to SAMAFS. When comparing the T2D groups and adjusting for sex and pedigree structure, individuals in the VACC were younger and had lower BMI, p-value <0.0001 (Table 7 and Figure 3) but showed the highest DBP levels and lowest HDL-C concentrations (Figure 3). The relationship between the triglycerides and BMI showed the classic inverted U-shape relationship in all three cohorts. However, the HDL-C levels decreased by the BMI strata and were found to be lowest for the VACC individuals (Table 7 and Figure 4).
Table 7. General linear models for age, BMI, blood pressure, triglycerides and HDL-C.
Figure 3. Violin plots depicting the relationship between T2D status and various demographic and cardiometabolic traits in the VACC and HIV-1 negative HHANES and SAMAFS individuals. (A) Violin plot showing relationship between T2D status and age in the three cohorts. Participants with T2D were older than non T2D [b = 13.1 (se 1.3), p < 0.001], but the VACC individuals showed interaction with T2D [b = −9.5 (se 1.9), p < 0.001]. (B) Violin plots showing relationship between age-, sex-, and cohort-adjusted BMI and T2D status in three different cohorts. Higher BMI was associated with T2D status [b = 2.8 (se 0.45) p < 0.001] in HHANES and SAMAFS, but interaction between VACC and T2D shows lower BMI when compared with other cohorts [b = −3.5 (se 0.55) p < 0.001]. (C,D) Violin plots showing relationship between SBP and DBP adjusted by age, sex, cohort, BMI and T2D status. For VACC, there was an interaction seen for DBP [b = 3.8 (se 1.0) p < 0.001], but not for SBP. (E,F) Violin plots showing relationship between triglycerides and HDL-C adjusted by age, sex, cohort, BMI and T2D status. T2D was associated with higher triglycerides concentrations [b = 28.2 (se 12.3), p = 0.02] and lower HDL-C [b = −3.5 (se 0.9) p < 0.001]. The HDL-C showed interaction for VACC [b = 2.2 (se 1.1) p = 0.042]. HHANES: The Hispanic component of NHANES 2018 cohort; SAMAFS, San Antonio Mexican American Family Studies; VACC, Valley AIDS Council Cohort. The continuous horizontal line represents the mean value without adjustment.
Figure 4. Violin plots showing relationship between triglycerides (A) and HDL-C (B) with BMI categories after adjusting for age, sex, cohort, and T2D status. The presence of diabetes was associated with higher triglyceride concentration [b = 28.2 (se 12.3), p = 0.02] and lower HDL-C [b = −3.5 (se 0.9) p < 0.001]. The HDL-C showed interaction for VACC [b = 2.2 (se 1.1) p = 0.042]. HHANES: The Hispanic component of NHANES 2018 cohort; SAMAFS, San Antonio Mexican American Family Studies; VACC, Valley AIDS Council Cohort. NW, normal weight; OW, overweight; OB-C1, obesity-category 1; OB-C2, obesity-category 2; and OB-C3, obesity category 3. The continuous horizontal line represents the mean value without adjustment.
Discussion
Our results suggest that the PWH in the RGV developed T2D at a younger age when compared to HIV-negative individuals in the SAMAFS and HHANES. There were striking differences between the cohorts: the PWH in VACC were younger and had lower BMI but had a higher prevalence of metabolic problems when compared to individuals in HHANES and SAMAFS. Furthermore, the increased T2D risk in PWH was independent of BMI. Patients with better virological control, irrespective of the treatment regimen, had a more favorable metabolic profile.
Previous studies have shown that insulin resistance and T2D prevalence are increased in individuals who have HIV-1 infection (3). However, there are only limited studies that have examined the prevalence of HIV and T2D in the RGV region. A recent analysis by De La Garza et al. showed that the average HIV incidence rates increased by 115% per 100,000 within the RGV from 2007 to 2015, which is higher than other Texas counties with large populations (26). RGV region also has a high prevalence of adult T2D and has been estimated to be as high as 30.7% compared to a nationwide prevalence of 12.3% (8, 9). Remarkably, Restrepo et al. found that the prevalence of T2D in individuals affected with tuberculosis in the Cameron and Hidalgo counties to be 39% (27). The presence of preexisting metabolic diseases and the de novo metabolic complications due to HIV-1 infection and ART could contribute to unusually high T2D prevalence seen in the VACC.
While recently introduced ART regimens have greatly improved metabolic outcomes, there have been reports of persistent metabolic dysfunction in PWH (28). For example, a longitudinal study conducted by Gomes et al. suggests that ART may increase T2D risk in Hispanics (25). They examined the development of new-onset impaired fasting glucose (IFG) following ART initiation in PWH in a Dominican HIV Hispanic cohort who were initiating ART (N = 153; ≤ 180 days prior to starting ART). At baseline, they found that 6% had T2D and 16% had IFG. Strikingly, in this cohort, 46 developed IFG after 329/1,000 person-years follow up. Most of the individuals developed metabolic complications within 12 months following the ART start, although only 13% were on PIs. Recent data also suggests that ART regimens that are initiated with INSTIs or PIs may confer a greater T2D risk, and this may be mediated through weight gain (29). They also reported that this effect is more pronounced in raltegravir-based initiators when compared to those based on NNRTIs. Our findings in VACC suggest that any regimen containing INSTI was associated with an increased prevalence of T2D. However, in the absence of longitudinal data, we cannot attribute this increased prevalence directly to a specific regimen. Recent studies have shown that INSTI-based regimens can increase T2D risk through multiple mechanisms, including an increased deposition of central adipose tissue and adverse effects on adipose tissue that may result in insulin resistance (30, 31). While some PI-based regimens have been associated with deleterious effects on insulin secretion (32), the impact of the integrase inhibitors on insulin secretion is not well-understood. It has been speculated that INSTIs may disrupt insulin secretion by interacting with magnesium ions (33) that are critical for electrical activity and insulin secretion in pancreatic beta-cells (34). Thus, it is imperative to conduct well-controlled longitudinal studies to understand the extent of metabolic dysfunction associated with various treatment regimens in the VACC to implement individualized preventative measures to reduce future cardiometabolic risk.
Obesity is thought to be associated with a chronic inflammatory state in the adipose tissue which is mediated through a number of inflammatory signals/cells [e.g., M1 macrophages, inflammatory cytokines (TNF-alpha and IL-6), and activated CD4+ and CD8+ T cells] (35), HIV/ART further exacerbate this condition through further modulation of inflammatory cytokines and mitochondrial dysfunction (36). Prior to the introduction of the newer ART regimens, lipodystrophy was a common feature in HIV-infected persons (37, 38). In contrast, recent studies have shown the increasing incidence of obesity and overweight in this group, which can be as high as 70% (39–43). In contrast to these findings, VACC individuals had a lower mean BMI when compared to individuals enrolled in SAMAFS, HHANES and RGV-based cohorts or their subsets (9, 44). However, in contrast to some RGV-based HIV-negative cohorts (9), we did not find a difference in the obesity prevalence in the diabetic and non-diabetic individuals in the VACC. Hernandez-Romieu et al. also found similar trends in the NHANES individuals. Using 2009–2010 data, they found that diabetes prevalence is much higher in HIV-infected adults when compared to the general U.S. adult population in the absence of obesity (45). Also, they also found that HIV-infected individuals had a higher prevalence of diabetes at younger ages. This higher T2D prevalence in PWH may be related to a chronic inflammatory state induced by HIV-1 infection and ART-induced insulin resistance.
Recent studies have suggested that the major classes of ART have differential effects on weight gain (31, 46–50) and abnormalities in fat and lipid storage (51). For example, Menard et al., reported that there was a mean weight gain of ~3 kg in a cohort (N = 462) receiving an integrase inhibitor-based ART regimen (dolutegravir) after a follow-up of an average of 276 days (46). Using a large study population (ACTG A5257 clinical trial; N = 1,809) in which the treatment naïve HIV were randomized into three different treatments: raltegravir (integrase inhibitors) or atazanavir/ritinovir (PIs) or darunavir/ritonavir (PIs), all in combination with tenofovir disoproxil fumerate/emtricitabine, Bhagwat et al., reported that blacks on integrase inhibitors had larger increases in waist circumference when compared to non-blacks (31). In the same study, they found that women had a greater increase in waist circumference when compared to men. In the VACC, we could not find any relationship between obesity prevalence and the broad ART regimens, and future longitudinal studies are needed to understand this phenomenon. Remarkably, we found that increasing BMI had a positive association with absolute CD4+ cell counts. This finding is in agreement with a previous study that also showed a positive relationship between higher BMI and CD4+ cell recovery in HIV-1 infected people on ART (52). Notably, such a positive relationship between higher BMI levels and CD4+ counts is also described in the absence of HIV-1 infection (53). The relevance of such an increase in CD4+ counts with increasing BMI in PWH needs additional studies.
Our study has some important limitations that are commonly encountered when extracting data from clinical records, including missing data, unavailability of records that are relevant for metabolic diseases such as waist circumference, physical activity, dietary habits, acculturation factors, co-infection with HCV, mode of disease transmission, prediabetes, among others. An important limitation is that we have relied on BMI as a major determinant to assess metabolic health in the VACC individuals in absence of other more appropriate obesity measures. While body composition measures such as waist circumference and fat mass would be more suitable to assess metabolic dysfunction in non-obese PWH due to HIV and ART-induced fat alterations, BMI categories have been successfully used to evaluate metabolic health in absence of such measures. For example, Lake et al. showed that in the Multicenter AIDS Cohort Study (MACS) normal weight and overweight HIV-positive men are much less likely to be metabolically healthy when compared to non-obese HIV-negative men (54). Given that the recent reports suggest that integrase inhibitor-based regimens could result in increased accumulation of both subcutaneous and visceral adipose tissue (55), we plan to conduct additional assessments for body composition in the VACC participants for future studies. Another limitation of our study is that the data were available for only 20 months due to a recent switch from a different electronic medical records database. Future studies will incorporate these additional factors and longitudinal designs to gain further insights into T2D pathogenesis in the VACC. Nevertheless, our analyses recapitulated several key aspects of HIV disease pathogenesis, including an inverse relationship between viral load and CD4+ cell count and CD4/CD8 ratio, thus suggesting our cross-sectional data are highly informative.
In summary, the prevalence of T2D in the PWH in the Rio Grande Valley is higher than the national and South Texas Hispanic populations. Remarkably, this increased prevalence is not associated with a higher BMI, suggesting that factors independent of obesity may be in part responsible for disease pathogenesis. Despite their lower BMI, these individuals showed increased dyslipidemia and hypertension, reflecting poor metabolic health. Further cross-sectional and longitudinal studies are necessary to develop personalized interventions to improve metabolic health and reduce future cardiovascular-related morbidity and mortality.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://wwwn.cdc.gov/nchs/nhanes/continuousnhanes/default.aspx?BeginYear=2017 and dbGaP Study Accession: phs000847.v1.p1.
Ethics Statement
The studies involving human participants were reviewed and approved by Institutional Review Board, University of Texas Rio Grande Valley. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
JL-A, DAM, RD, SM, and RDM contributed to the study conception and study design. DAM and RDM were responsible for the clinical care of the PWH and supervised the collection of clinical data from VACC. SAMAFS data was collected under the supervision of DL, RD, JEC, and JB. JL-A performed the statistical analysis with contributions from AD-B, SM, and RD. AD-B, LM, RA, and CPJ contributed to the data analysis and interpretation of the data. JL-A wrote the initial draft of the manuscript, which was critically revised by DAM, RD, SM, and RDM. All authors contributed to the article and approved the submitted version.
Funding
The work was supported in part by the National Institutes of Health grants, R01 AI119131, P01 HL045522, R01 DK047482, and R01 DK053889, and the McAllen Medical Center Foundation Endowed Professor in Diabetes Research Award, School of Medicine, UTRGV, Edinburg, TX.
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.
Acknowledgments
We dedicate this manuscript to the memory of Dr. Virginia L. Scofield (1948–2017), School of Medicine, UTRGV, whose untiring efforts have led to the initiation of this study. We acknowledge the help of the staff of the VAC and Ms. Jennifer Schneider in the extraction of data from the VAC database. We thank the participants of the San Antonio Family Heart Study and the San Antonio Family Diabetes/Gallbladder Study for their continued cooperation and participation in our research programs.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.676979/full#supplementary-material
Supplementary Figure 1. Histograms showing the distribution of raw variables (A) and after transformation (B). The frequency of CD4 absolute counts is in the upper panel. The lower panel shows the frequency of the squared root variable. The transformation resulted in an approximately normal distribution.
Supplementary Table 1. General description and stratification by sex of selected variables from the SAMAFS¶.
Supplementary Table 2. General description and contrast by sex of selected variables from HHANES¶ 2018 data.
References
1. Benbow ND, Aaby DA, Rosenberg ES, Brown CH. County-level factors affecting Latino HIV disparities in the United States. PLoS ONE. (2020) 15:e0237269. doi: 10.1371/journal.pone.0237269
2. Lorenzo C, Hazuda HP, Haffner SM. Insulin resistance and excess risk of diabetes in Mexican-Americans: the San Antonio Heart Study. J Clin Endocrinol Metab. (2012) 97:793–9. doi: 10.1210/jc.2011-2272
3. Brown TT, Cole SR, Li X, Kingsley LA, Palella FJ, Riddler SA, et al. Antiretroviral therapy and the prevalence and incidence of diabetes mellitus in the multicenter AIDS cohort study. Arch Intern Med. (2005) 165:1179–84. doi: 10.1001/archinte.165.10.1179
4. Nansseu JR, Bigna JJ, Kaze AD, Noubiap JJ. Incidence and risk factors for prediabetes and diabetes mellitus among HIV-infected adults on antiretroviral therapy: a systematic review and meta-analysis. Epidemiology. (2018) 29:431–41. doi: 10.1097/EDE.0000000000000815
5. Koethe JR, Jenkins CA, Lau B, Shepherd BE, Justice AC, Tate JP, et al. North American and design: rising obesity prevalence and weight gain among adults starting antiretroviral therapy in the United States and Canada. AIDS Res Hum Retroviruses. (2016) 32:50–8. doi: 10.1089/aid.2015.0147
6. Sax PE, Erlandson KM, Lake JE, McComsey GA, Orkin C, Esser S, et al. Weight gain following initiation of antiretroviral therapy: risk factors in randomized comparative clinical trials. Clin Infect Dis. (2020) 71:1379–89. doi: 10.1093/cid/ciz999
7. Han JH, Gordon K, Womack JA, Gibert CL, Leaf DA, Rimland D, et al. Comparative effectiveness of diabetic oral medications among HIV-infected and HIV-uninfected veterans. Diabetes Care. (2017) 40:218–25. doi: 10.2337/dc16-0718
8. Millard AV, Graham MA, Mier N, Moralez J, Perez-Patron M, Wickwire B, et al. Diabetes screening and prevention in a high-risk, medically isolated border community. Front Public Health. (2017) 5:135. doi: 10.3389/fpubh.2017.00135
9. Watt GP, Fisher-Hoch SP, Rahbar MH, McCormick JB, Lee M, Choh AC, et al. Mexican American and South Asian population-based cohorts reveal high prevalence of type 2 diabetes and crucial differences in metabolic phenotypes. BMJ Open Diabetes Res Care. (2018) 6:e000436. doi: 10.1136/bmjdrc-2017-000436
10. Ryabov I, Merino S. Recent demographic change in the Rio Grande Valley of Texas: the importance of domestic migration. J Borderlands Stud. (2017) 32:211–31. doi: 10.1080/08865655.2016.1195704
11. Weir CB, Jan A. BMI classification percentile and cut off points. In: StatPearls. Treasure Island, FL: StatPearls Publishing (2020).
12. American Diabetes Association. 2. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes-2020. Diabetes Care. 43(Suppl 1):S14–31. doi: 10.2337/dc20-S002
13. Olvera RL, Bearden CE, Velligan DI, Almasy L, Carless MA, Curran JE, et al. Common genetic influences on depression, alcohol, and substance use disorders in Mexican-American families. Am J Med Genet B Neuropsychiatr Genet. (2011) 156B:561–8. doi: 10.1002/ajmg.b.31196
14. Mitchell BD, Kammerer CM, Blangero J, Mahaney MC, Rainwater DL, Dyke B, et al. Genetic and environmental contributions to cardiovascular risk factors in Mexican Americans. The San Antonio Family Heart Study. Circulation. (1996) 94:2159–70. doi: 10.1161/01.CIR.94.9.2159
15. Hunt KJ, Lehman DM, Arya R, Fowler S, Leach RJ, Goring HH, et al. Genome-wide linkage analyses of type 2 diabetes in Mexican Americans: the San Antonio Family Diabetes/Gallbladder Study. Diabetes. (2005) 54:2655–62. doi: 10.2337/diabetes.54.9.2655
16. Puppala S, Dodd GD, Fowler S, Arya R, Schneider J, Farook VS, et al. A genomewide search finds major susceptibility loci for gallbladder disease on chromosome 1 in Mexican Americans. Am J Hum Genet. (2006) 78:377–92. doi: 10.1086/500274
17. Arya R, Blangero J, Williams K, Almasy L, Dyer TD, Leach RJ, et al. Factors of insulin resistance syndrome–related phenotypes are linked to genetic locations on chromosomes 6 and 7 in nondiabetic mexican-americans. Diabetes. (2002) 51:841–7. doi: 10.2337/diabetes.51.3.841
18. Comuzzie AG, Hixson JE, Almasy L, Mitchell BD, Mahaney MC, Dyer TD, et al. A major quantitative trait locus determining serum leptin levels and fat mass is located on human chromosome 2. Nat Genet. (1997) 15:273–6. doi: 10.1038/ng0397-273
19. Duggirala R, Blangero J, Almasy L, Dyer TD, Williams KL, Leach RJ, et al. Stern: linkage of type 2 diabetes mellitus and of age at onset to a genetic location on chromosome 10q in Mexican Americans. Am J Hum Genet. (1999) 64:1127–40. doi: 10.1086/302316
20. Curran JE, McKay DR, Winkler AM, Olvera RL, Carless MA, Dyer TD, et al. Identification of pleiotropic genetic effects on obesity and brain anatomy. Hum Hered. (2013) 75:136–43. doi: 10.1159/000353953
21. World Health Organization. Definition, Diagnosis and Classification of Diabetes Mellitus and Its Complications: Report of a WHO Consultation. Part 1, Diagnosis and Classification of Diabetes Mellitus. Geneva: World Health Organization (1999).
22. The Expert Committee on Diagnosis and Classification of Diabetes Mellitus. Report of the The Expert Committee on Diagnosis and Classification of Diabetes Mellitus. Diabetes Care. (1998) 21:S5–19. doi: 10.2337/diacare.21.1.S5
23. Centers for Disease Control and Prevention (CDC). National Center for Health Statistics (NCHS). Hyattsville, MD: National Health and Nutrition Examination Survey Data (2020).
24. Huber P. The behavior of maximum likelihood estimates under nonstandard condition. In: LeCam NM, Neyman, J, editors. Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics Probability. Oakland, CA: University of California Press (1967).
25. Gomes A, Reyes EV, Garduno LS, Rojas R, Mir Mesejo G, Del Rosario E, et al. Incidence of diabetes mellitus and obesity and the overlap of comorbidities in HIV+ hispanics initiating antiretroviral therapy. PLoS ONE. (2016) 11:e0160797. doi: 10.1371/journal.pone.0160797
26. De La Garza R, Rodrigo H, Fernandez F, Roy U. The increase of HIV-1 infection, neurocognitive impairment, and type 2 diabetes in The Rio Grande Valley. Curr HIV Res. (2019) 17:377–87. doi: 10.2174/1570162X17666191029162235
27. Restrepo BI, Camerlin AJ, Rahbar MH, Wang W, Restrepo MA, Zarate I, et al. Cross-sectional assessment reveals high diabetes prevalence among newly-diagnosed tuberculosis cases. Bull World Health Organ. (2011) 89:352–9. doi: 10.2471/BLT.10.085738
28. Willig AL, Overton ET. Metabolic complications and glucose metabolism in HIV infection: a review of the evidence. Curr HIV/AIDS Rep. (2016) 13:289–96. doi: 10.1007/s11904-016-0330-z
29. Rebeiro PF, Jenkins CA, Bian A, Lake JE, Bourgi K, Moore RD, et al. Risk of incident diabetes mellitus, weight gain, and their relationships with integrase inhibitor-based initial antiretroviral therapy among persons with HIV in the US and Canada. Clin Infect Dis. (2020). doi: 10.1093/cid/ciaa1403. [Epub ahead of print].
30. Gorwood J, Bourgeois C, Pourcher V, Pourcher G, Charlotte F, Mantecon M, et al. The integrase inhibitors dolutegravir and raltegravir exert pro-adipogenic and profibrotic effects and induce insulin resistance in human/simian adipose tissue and human adipocytes. Clin Infect Dis. (2020) 71:e549–60. doi: 10.1093/cid/ciaa259
31. Bhagwat P, Ofotokun I, McComsey GA, Brown TT, Moser C, Sugar CA, et al. Changes in waist circumference in HIV-infected individuals initiating a raltegravir or protease inhibitor regimen: effects of sex and race. Open Forum Infect Dis. (2018) 5:ofy201. doi: 10.1093/ofid/ofy201
32. Lee GA, Rao MN, Grunfeld C. The effects of HIV protease inhibitors on carbohydrate and lipid metabolism. Curr Infect Dis Rep. (2004) 6:471–82. doi: 10.1007/s11908-004-0067-5
33. Fong PS, Flynn DM, Evans CD, Korthuis PT. Integrase strand transfer inhibitor-associated diabetes mellitus: a case report. Int J STD AIDS. (2017) 28:626–8. doi: 10.1177/0956462416675107
34. Kostov K. Effects of magnesium deficiency on mechanisms of insulin resistance in type 2 diabetes: focusing on the processes of insulin secretion and signaling. Int J Mol Sci. (2019) 20:1351. doi: 10.3390/ijms20061351
35. Gonzalez-Muniesa P, Martinez-Gonzalez MA, Hu FB, Despres JP, Matsuzawa Y, Loos RJF, et al. Obesity. Nat Rev Dis Primers. (2017) 3:17034. doi: 10.1038/nrdp.2017.34
36. Lake JE. The fat of the matter: obesity and visceral adiposity in treated HIV infection. Curr HIV/AIDS Rep. (2017) 14:211–9. doi: 10.1007/s11904-017-0368-6
37. Carr A, Samaras K, Burton S, Law M, Freund J, Chisholm DJ, et al. A syndrome of peripheral lipodystrophy, hyperlipidaemia and insulin resistance in patients receiving HIV protease inhibitors. AIDS. (1998) 12: F51–8. doi: 10.1097/00002030-199807000-00003
38. Carr A, Samaras K, Thorisdottir A, Kaufmann GR, Chisholm DJ, A D. Cooper: diagnosis, prediction, and natural course of HIV-1 protease-inhibitor-associated lipodystrophy, hyperlipidaemia, and diabetes mellitus: a cohort study. Lancet. (1999) 353:2093–9. doi: 10.1016/S0140-6736(98)08468-2
39. Achhra AC, Mocroft A, Reiss P, Sabin C, Ryom L, de Wit S, et al. Short-term weight gain after antiretroviral therapy initiation and subsequent risk of cardiovascular disease and diabetes: the D:A:D study. HIV Med. (2016) 17:255–68. doi: 10.1111/hiv.12294
40. Buchacz K, Baker RK, Palella FJ Jr, Shaw L, Patel P, Lichtenstein KA, et al. Investigators VOS, Disparities in prevalence of key chronic diseases by gender and race/ethnicity among antiretroviral-treated HIV-infected adults in the US. Antivir Ther. (2013) 18:65–75. doi: 10.3851/IMP2450
41. Erlandson KM, Taejaroenkul S, Smeaton L, Gupta A, Singini IL, Lama JR, et al. A randomized comparison of anthropomorphic changes with preferred and alternative efavirenz-based antiretroviral regimens in diverse multinational settings. Open Forum Infect Dis. (2015) 2:ofv095. doi: 10.1093/ofid/ofv095
42. Hernandez D, Kalichman S, Cherry C, Kalichman M, Washington C, Grebler T. Dietary intake and overweight and obesity among persons living with HIV in Atlanta Georgia. AIDS Care. (2017) 29:767–71. doi: 10.1080/09540121.2016.1238441
43. Levy ME, Greenberg AE, Hart R, Powers Happ L, Hadigan C, Castel A, et al. Committee: High burden of metabolic comorbidities in a citywide cohort of HIV outpatients: evolving health care needs of people aging with HIV in Washington, DC. HIV Med. (2017) 18:724–35. doi: 10.1111/hiv.12516
44. Manusov EG, Diego VP, Smith J, Garza JR 2nd, Lowdermilk J, Blangero J, et al. UniMovil: a mobile health clinic providing primary care to the colonias of the Rio Grande Valley, South Texas. Front Public Health. (2019) 7:215. doi: 10.3389/fpubh.2019.00215
45. Hernandez-Romieu AC, Garg S, Rosenberg ES, Thompson-Paul AM, Skarbinski J. Is diabetes prevalence higher among HIV-infected individuals compared with the general population? Evidence from MMP and NHANES 2009-2010. BMJ Open Diabetes Res Care. (2017) 5:e000304. doi: 10.1136/bmjdrc-2016-000304
46. Menard A, Meddeb L, Tissot-Dupont H, Ravaux I, Dhiver C, Mokhtari S, et al. Dolutegravir and weight gain: an unexpected bothering side effect? AIDS. (2017) 31:1499–500. doi: 10.1097/QAD.0000000000001495
47. Norwood J, Turner M, Bofill C, Rebeiro P, Shepherd B, Bebawy S, et al. Brief report: weight gain in persons with HIV Switched from efavirenz-based to integrase strand transfer inhibitor-based regimens. J Acquir Immune Defic Syndr. (2017) 76:527–31. doi: 10.1097/QAI.0000000000001525
48. Bourgi K, Rebeiro PF, Turner M, Castilho JL, Hulgan T, Raffanti SP, et al. Greater weight gain in treatment-naive persons starting dolutegravir-based antiretroviral therapy. Clin Infect Dis. (2020) 70:1267–74. doi: 10.1093/cid/ciz407
49. Kerchberger AM, Sheth AN, Angert CD, Mehta CC, Summers NA, Ofotokun I, et al. Weight gain associated with integrase stand transfer inhibitor use in women. Clin Infect Dis. (2020) 71:593–600. doi: 10.1093/cid/ciz853
50. Kline MDA, Daniels C, Xu X, Sunil T, Ganesan A, Agan BK, et al. Antiretroviral therapy anchor-based trends in body mass index following treatment initiation among military personnel with HIV. Mil Med. (2020) 186:279–85. doi: 10.1093/milmed/usaa416
51. Koethe JR, Lagathu C, Lake JE, Domingo P, Calmy A, Falutz J, et al. HIV and antiretroviral therapy-related fat alterations. Nat Rev Dis Primers. (2020) 6:48. doi: 10.1038/s41572-020-0181-1
52. Palermo B, Bosch RJ, Bennett K, Jacobson JM. Body mass index and CD4+ T-lymphocyte recovery in HIV-infected men with viral suppression on antiretroviral therapy. HIV Clin Trials. (2011) 12:222–7. doi: 10.1310/hct1204-222
53. Womack J, Tien PC, Feldman J, Shin JH, Fennie K, Anastos K, et al. Obesity and immune cell counts in women. Metabolism. (2007) 56:998–1004. doi: 10.1016/j.metabol.2007.03.008
54. Lake JE, Li X, Palella FJ Jr, Erlandson KM, Wiley D, Kingsley L, et al. Metabolic health across the BMI spectrum in HIV-infected and HIV-uninfected men. AIDS. (2018) 32:49–57. doi: 10.1097/QAD.0000000000001651
Keywords: South Texas, AIDS, type 2 diabetes, HIV, Mexican Americans, antiretroviral treatment, Rio Grande Valley, obesity
Citation: Lopez-Alvarenga JC, Martinez DA, Diaz-Badillo A, Morales LD, Arya R, Jenkinson CP, Curran JE, Lehman DM, Blangero J, Duggirala R, Mummidi S and Martinez RD (2021) Association of HIV-1 Infection and Antiretroviral Therapy With Type 2 Diabetes in the Hispanic Population of the Rio Grande Valley, Texas, USA. Front. Med. 8:676979. doi: 10.3389/fmed.2021.676979
Received: 06 March 2021; Accepted: 26 May 2021;
Published: 05 July 2021.
Edited by:
Alfredo Ulloa-Aguirre, National Autonomous University of Mexico, MexicoReviewed by:
Liliana Munoz, National Institute of Medical Science and Nutrition Salvador Zubiran, MexicoAlvaro López Iñiguez, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán (INCMNSZ), Mexico
Copyright © 2021 Lopez-Alvarenga, Martinez, Diaz-Badillo, Morales, Arya, Jenkinson, Curran, Lehman, Blangero, Duggirala, Mummidi and Martinez. 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: Srinivas Mummidi, srinivas.mummidi@utrgv.edu; Ruben D. Martinez, rmartinez@westbrookclinic.org
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share senior authorship