 Research article
 Open Access
 Open Peer Review
 Published:
A comparison of machine learning techniques for classification of HIV patients with antiretroviral therapyinduced mitochondrial toxicity from those without mitochondrial toxicity
BMC Medical Research Methodology volume 19, Article number: 216 (2019)
Abstract
Background
Antiretroviral therapy (ART) has significantly reduced HIVrelated morbidity and mortality. However, therapeutic benefit of ART is often limited by delayed drugassociated toxicity. Nucleoside reverse transcriptase inhibitors (NRTIs) are the backbone of ART regimens. NRTIs compete with endogenous deoxyribonucleotide triphosphates (dNTPs) in incorporation into elongating DNA chain resulting in their cytotoxic or antiviral effect. Thus, the efficacy of NRTIs could be affected by direct competition with endogenous dNTPs and/or feedback inhibition of their metabolic enzymes. In this paper, we assessed whether the levels of ribonucleotides (RN) and dNTP pool sizes can be used as biomarkers in distinguishing between HIVinfected patients with ARTinduced mitochondrial toxicity and HIVinfected patients without toxicity.
Methods
We used data collected through a casecontrol study from 50 subjects. Cases were defined as HIVinfected individuals with clinical and/or laboratory evidence of mitochondrial toxicity. Each case was age, gender, and race matched with an HIVpositive without evidence of toxicity. We used a range of machine learning procedures to distinguish between patients with and without toxicity. Using resampling methods like Monte Carlo kfold cross validation, we compared the accuracy of several machine learning algorithms applied to our data. We used the algorithm with highest classification accuracy rate in evaluating the diagnostic performance of 12 RN and 14 dNTP pool sizes as biomarkers of mitochondrial toxicity.
Results
We used eight classification algorithms to assess the diagnostic performance of RN and dNTP pool sizes distinguishing HIV patients with and without NRTIassociated mitochondrial toxicity. The algorithms resulted in crossvalidated classification rates of 0.65–0.76 for dNTP and 0.72–0.83 for RN, following reduction of the dimensionality of the input data. The reduction of input variables improved the classification performance of the algorithms, with the most pronounced improvement for RN. Complex treebased methods worked the best for both the deoxyribose dataset (Random Forest) and the ribose dataset (Classification Tree and AdaBoost), but it is worth noting that simple methods such as Linear Discriminant Analysis and Logistic Regression were very competitive in terms of classification performance.
Conclusions
Our finding of changes in RN and dNTP pools in participants with mitochondrial toxicity validates the importance of dNTP pools in mitochondrial function. Hence, levels of RN and dNTP pools can be used as biomarkers of ARTinduced mitochondrial toxicity.
Background
Although significant progress has been made to contain the AIDS epidemic, the Joint United Nations Program on HIV/AIDS (UNAIDS) estimates that there are currently over 36.9 [34.3–41.4] million people globally living with HIV [UNAIDS, 2015]. According to UNAIDS, 1.20 [0.98–1.60] million people died from AIDSrelated causes worldwide in 2014. Antiretroviral therapy (ART) has significantly reduced HIVrelated morbidity and mortality worldwide, however, its therapeutic benefit is compromised by delayed drugassociated toxicity [1]. A recent update from the Italian of Cohort Naïve Antiretrovirals (ICONA) Foundation study showed the 1year probability of discontinuation of ART due to toxicity was 19% for patients who initiated ART between 2008 and 2014 [2]. Nucleoside reverse transcriptase inhibitors (NRTIs), the first class of ART used in the treatment of individuals with HIV, have been associated with toxicities that mirror symptoms in individuals with mitochondrial disorders. The presentation of NRTIinduced mitochondrial toxicity includes lactic acidosis, lipodystrophy, peripheral neuropathies, cardiomyopathies, and pancytopenia [1, 3,4,5,6]. These manifestations have been ascribed to the inhibitory effect of NRTIs on polymerase gamma (Polɣ), the enzyme that replicates mitochondrial DNA (mtDNA). Inhibition of Polɣ leads to depletion of mtDNA content and subsequent mitochondrial dysfunction. However, there are emerging reports that other classes of ART such as nonnucleoside reverse transcriptase inhibitors (NNRTIs) and protease inhibitors (PIs) can cause mitochondrial dysfunction through Polɣindependent mechanisms [7, 8].
We recently reported that HIVinfected patients on NRTIbased ART with mitochondrial toxicity tended to have decreased concentrations of both ribonucleotide (RN) and deoxyribonucleotide (dNTP) pool sizes (precursors of DNA synthesis) in their peripheral blood mononuclear cells (PBMCs) [9]. Interestingly, patients with mitochondrial toxicity had statistically significant higher mRNA expression of Polɣ compared to patients without mitochondrial toxicity. In this paper, we aimed to determine the effectiveness of RN and dNTP pool sizes as biomarkers in distinguishing between HIV patients with and without mitochondrial toxicity. More specifically, the objective of this paper is twofold: First, we assessed whether the levels of RN and dNTP pool sizes can be used as biomarkers in distinguishing between HIVinfected patients with ARTinduced mitochondrial toxicity and HIVinfected patients without toxicity. Second, we assessed whether a range of machine learning procedures could distinguish between patients with and without toxicity. In this study, where we have both the limited sample size and dataset with varying correlation structure, it will be particularly prudent to consider many different machine learning algorithms to assess their performances in a meaningful way. We used previously acquired data on ARTinduced mitochondrial toxicity to evaluate the relative performance of eight different machine learning procedures. Using resampling methods like Monte Carlo kfold cross validation, we compared the accuracy of the machine learning algorithms applied to our data. We ran a large set of simulation studies to evaluate the performance of the machine learning algorithms for varying sample sizes and correlation between RN and dNTP pool sizes. We used the algorithm with highest classification accuracy rate in evaluating the diagnostic performance of RN and dNTP pool sizes as biomarkers of mitochondrial toxicity.
Methods
Study participants
The study included 50 HIVinfected individuals on NRTIbased ART regimens enrolled through a casecontrol study from April 2011 to March 2013. Twentyfive of the HIVinfected individuals were diagnosed by their providers as having mitochondrial toxicity based on at least one clinical or laboratory findings associated with NRTI toxicity (i.e., lactic acidemia, pancreatitis, peripheral neuropathy, lipodystrophy, creatinine, anemia, renal dysfunction, elevated liver enzymes, hyperlipidemia, amylase, and lipase) [1, 3,4,5,6]. The study participants were recruited from the Adults AIDS Care Programs at YaleNew Haven Hospital. At enrollment participants answered a brief survey comprised of questions regarding past medical history and demographic characteristics. Medical records of participants were reviewed for date of HIV diagnosis, medication history, date of diagnosis of toxicity, complete blood count, serum chemistries, liver function test, lipid profile, urinalysis, HIV RNA copy number, and CD4+ Tcell count at study entry. Blood was collected from all study participants at the time of enrollment for extraction and quantification of ribonucleotide (RN) and deoxyribonucleotide (dNTP) pool sizes. The study protocol was approved by the Institutional Review Board of the Yale University School of Medicine. Informed and written consents were obtained from study participations prior to enrollment in the study. The rationale, organization, and recruitment of the subjects, biological procedures used for extraction and quantification RN and dNTP pools have been described previously [9].
Statistical methods
Data classification methods
There are many classification methods, some of which can be very sophisticated and stateoftheart [10]. In this paper, we focused on the most common methods which are readily available and can easily be implemented in most statistical packages. We briefly described the classification methods used in the paper below:
A. Linear and quadratic discriminant analysis
The linear discriminant analysis (LDA) and the quadratic discriminant analysis (QDA) are classical statistical classification methods [11]. In classifying two groups, both methods incorporate log likelihood ratio based on normal distribution, with the main difference that LDA assumes the equal covariance while QDA does not.
B. Knearest neighbor
The knearest neighbor (KNN) is a simple classification method, using the k nearest points of the input to predict the response [12]. Here, we first predict the value Y based on x, by
where N_{k}(x) is the neighborhood of x defined by the k closest input points x_{i} to x (using Euclidean distance in most cases), and k values of y_{i} corresponding to x_{i}, where y_{i} = 0 or y_{i} = 1, at each neighborhood N_{k}(x). Hence, \( \hat{Y}(x) \) is the proportion of time that y_{i} = 1 in the neighborhood N_{k}(x) for the point x, and we make the final decision as 1 if \( \hat{Y}(x) \) > 0.5 and as 0 otherwise, at each point x. We may use different values of k for the kneighbor. As k decreases, the model can be more accurate on the training data, however it runs the risk of overfitting.
C. Logistic regression
Logistic regression is a simple method that can be used to predict the outcome of the input variables [13]. We denote x = (x_{1}, ..., x_{p}) as the input variables and y as the response (y can be any two values). If we let y = 0 denote failure and y = 1 denote success, we have that
or
which we interpret as the probability of success (y = 1) given the data (x). The betas β = (β_{1}, …, β_{p}) and β_{0} are the coefficient values estimated from the model. Given input data, many statistical packages fit a logistic regression model and return the coefficient values. Then, given new data values x, we can use the above equation to estimate Pr(y = 1 x), which will be between zero and one, inclusive. Given the input x of an individual, we classify the individual as success if Pr(y = 1 x) > c where c is a cutoff, usually set at c = 0.5. Note that this method is linear in terms of parameters and will have difficulty fitting any data with a large number of input variables (p), even if p is smaller than the number of observations (n). In addition, if the input data follow exactly the pattern of the outcome variable, we have the socalled “complete separation” problem, which is again pronounced if we have large p. However, some of these shortcomings can easily be overcome with Firth’s bias adjustment [14].
D. Support vector machine
A support vector machine (SVM) is a machine learning method for binary classification. It uses a linear separating hyperplane f(x) = x^{T}β + β_{0} to split the pdimensional sample space into two groups. The SVM transforms nonlinear classification into a simpler linear classification problem, using a kernel function K(x, x_{i}), with the separating hyperplane
and the classification criteria, sign[f(x)]. The optimization criterion is that we maximize the margin (support vector) of the separating hyperplane to obtain the optimal separation, where margin is defined as \( M=\frac{1}{\mid \left\beta \right\mid } \), where \( {\left\beta \right}_2=\sqrt{\beta_1^2+\cdots +{\beta}_p^2} \). The typical choice for the kernel is a radial basis (Gaussian) kernel, K(x, x_{i}) = exp(γx − x_{i}^{2}) which is the default for most SVM software. Other kernels such as the polynomial kernel K(x, x_{i}) = (1 +〈x, x_{i}〉)^{d}, or the neural network (hyperbolic tangent) kernel K(x, x_{i}) = tanh(κ_{1}〈x, x_{i}〉 + κ_{2}) may be used, but the Gaussian kernel is the most popular because it derives many desirable properties from its association with the Gaussian distribution. As the SVM method involves nonlinear kernel and optimization, it can be very computationally intensive as compared to logistic regression, but the SVM is designed for a large number of input variables (p > n), which logistic regression cannot handle.
E. Classification trees
The classification tree method looks for the best classification of data by splitting each variable recursively and finding the optimal combination [15]. In other words, if we are given the data y and x = (x_{1}, ..., x_{p}) the classification tree looks for the best split points (t_{1}, ..., t_{p}) that gives us the decision rule. For example, if we have three input variables x = (x_{1}, x_{2}, x_{3}) x_{1}, x_{2}, x_{3} each taking values between 0 and 10, then the classification tree algorithm may provide the split points t_{1} = 5, t_{2}, = 8, t_{3} = 4.5 such that we declare an input (x_{1}, x_{2}, x_{3}) as success if x_{1} ≥ 5, x_{2} ≤ 8, x_{3} ≥ 4.5, and failure otherwise. The determination of split points largely depends on the algorithm used, for which we have many choices. Nevertheless, we see that it is easy to understand conceptually and is a popular method for classification. On the other hand, since any treebased algorithm will involve recursively partitioning all variables to find optimal splits, generating the classification tree will also be time consuming.
F. AdaBoost
The classification trees, although simple to understand and implement, may be too naïve, especially if the input variables have complicated structure. The AdaBoost algorithm is a refinement of the classification tree, where the classification trees are fit recursively to determine the final classifier via a majority vote [16]. The majority voting is defined as taking the results of many subclassifiers (here, the recursively fit trees) and making the final decision 0 or 1 from the majority of times that the subclassifiers vote 0 or 1 [17]. The AdaBoost uses the weighted majority voting in determining the final decision, where the 0 or 1 votes are weighted based on the accuracy of subclassifiers. Here, the weights are defined as log((err/(1err)) at each iteration, where err is defined as weighted error rate [10]. Hence, as an iterative procedure that updates the weight at each iteration, the AdaBoost will almost always improve over classification tree. However, the AdaBoost will be even more computationally intensive than the classification trees, requiring more time to run.
G. Random Forest
The random forest is another method based on the classification tree [18]. Similar to AdaBoost, the random forest also involves taking a majority vote of a sample of trees to determine the final classifier, but the building of the trees involves a statistical technique known as the bootstrap. The bootstrap method repeatedly resampling the data with replacement to obtain an ensemble of trees from the original sample. To obtain the final classifier, we take a majority vote of the bootstrapped trees. Hence, the random forest is yet another recursive partitioning method that involves resampling of data. As with the AdaBoost, its performance will improve over classification tree results; however, the computation burden is greater.
II. Data dimension reduction methods
Classifiers, such as logistic regression, work best when the number of input variables (p) is smaller than the number of observations (n). However, if p > n or if p is moderately large (as is the case for our problem), we cannot correctly fit some of the classifier with the full input. Therefore, reducing the input dimension sometimes increases the performance of the classifier, and is a good practice in general for a statistical analysis. We describe here some of the dimension reduction techniques.
A. Principal component analysis
The most common method of reducing the dimension of input variable is the principal component analysis (PCA) [19]. In the PCA, one transforms the original variables into a set of uncorrelated orthogonal basis vectors (components), which explain monotonically decreasing amounts of variance in the original data. However, the interpretation of these bases vectors becomes difficult because the biological meaning of these principle components is not always clear, so in this paper, we adopted a technique that simply selects the most important input variables, rather than transforming them.
B. Shrinkage methods (ridge, LASSO, and elastic net)
Recently, the variable selection methods such as the least absolute shrinkage and selection operator (LASSO) and elastic net have gained popularity in highdimensional statistical problems [20, 21]. We used the glmnet function from glmnet package in R, which implements the elastic net where the LASSO is a special case [22]. In the glmnet function, it is flexible to implement the ridge regression (α = 0) and LASSO (α = 1), and anything in between (elastic net, 0 ≤ α ≤ 1). In other words, the elastic net solves
Hence, the α term can be flexible to control how variables are selected based on the coefficients β.
For our paper, we focused on LASSO (α = 1). In LASSO, the procedure simply selects variables based on the L1 hard threshold of the coefficients, ‖β‖_{1} (as opposed to ridge regression which only smooths out, or shrinks, the coefficients [10]). Any other α values smaller than one will also have the effect of shrinkage, which we do not want here. Since the (α = 1) is fixed, it remains to calculate the other parameter, λ, to complete the process of selecting the variables via LASSO. The glmnet function in R recommends that the users view the entire solution path consisting of results from all possible λ values, but practically this is unfeasible so the program selects two plausible λ values: lambda.min, the value of lambda that gives minimum mean crossvalidated error, or lambda.1se, the largest value of lambda such that the error is within one standarderror of the minimum—the so called “onestandarderror” rule [22]. Hence, we consider two sets of selected variables based on lambda.min and lambda.1se.
Normally, one assigns test and validation sets within sample to perform classifications. However, since we have a limited sample size, we employ the kfold cross validation (CV) [10]. For the kfold CV, the sample is randomly divided into k roughly equal sized sets, and one of the K sets is left out. A classifier is fit with the K − 1 sets (training set) and validated using the remaining (left out) set. We repeat this for each of the K sets; the kfold CV classification rate is obtained by averaging the k individual classification rates. Now, we repeat the process M times (Monte Carlo simulation of the kfold CV), because the partition of the kfold is random so that each time it gives us a different result. Hence, we obtain a more reliable estimate of the kfold CV by repeating the procedure M times and also obtain the distribution of the kfold CV classification rates, including the mean and the standard deviation (standard error). We set M = 1,000 and report the 1000 Monte Carlo mean of the kfold CV classification rates. The computation times for the eight classifiers range 2–3 s for all algorithms except for AdaBoost (177 s), random forest (9 s) and KNN (< 1 s).
Results
In this section, we present results from our simulation studies and analyses of the mitochondrial toxicity data. To obtain the classification rates, we fit the models using each of the machine learning algorithms described in the methods section and predict the outcome using simulated as well as the original mitochondrial toxicity data.
Simulation studies
Simulation procedures
We ran a large set of simulation studies to assess the impact of sample size and correlation on the performance of the machine learning algorithms. We considered various sample sizes (n_{1} = n_{2} = 25, n_{1} = n_{2} = 50, n_{1} = n_{2} = 100). In each simulation, the outcome variable D ∈ {0, 1} (without toxicity or with toxicity) was created (for example, n_{1} = 25 were assigned 0 while n_{2} = 25 were assigned 1). For subjects with and without toxicity, we generated 12 (p = 12) correlated predictors from a multivariate distribution with mean vectors \( \overrightarrow{1} \) and \( \overrightarrow{0} \), respectively. Three distributions were considered: (1) multivariate normal, (2) multivariatet with 3 degrees of freedom, and (3) multivariate Cauchy (t with 1 degree of freedom). We considered other distributions besides the normal distribution as some classifiers (such as LDA and SVM) perform overwhelmingly well under normality because of their algorithm structures. Due to the autoregressive nature of the mitochondrial toxicity data, we considered First Order Autoregressive AR(1) in addition to compound symmetry when specifying the covariance structure in our simulation. That is, for subjects with and without toxicity the predictors were generated from the aforementioned multivariate distributions with the following correlation matrices: (1) compound symmetry with offdiagonal correlation (ρ) values ranging from 0 to 0.9; and (2) AR(1) with ρ = 0.4 and ρ = 0.8. We have considered the following classifiers: Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), KNearest Neighbor (KNN), Support Vector Machine (SVM), Classification Tree (CART), AdaBoost (ADA), Random Forest (RF), and Logistic Regression (LOGIT). Depending on the study, one may need to finetune the parameter values of the classifiers to obtain desirable results. We have tried many ways of carefully tuning the parameter values, but the tuning led to either overfitting or inaccuracy; the default values yielded a good balance between the two tradeoffs and performed well in our setting. Thus, we used default parameter values with the exception of KNN (we used k = 3 instead of the default k = 1).
Simulation results
The average classification rates of our simulation studies are presented in Fig. 1 and the corresponding standard errors were all less than 0.05. The classification rates were obtained using 5fold CV and 1000 Monte Carlo runs. As shown in the figure, the main findings of our simulation results were: (1) all the classification algorithms work better when data were least correlated, and the sample sizes were large, (2) classifiers perform better when data were normally distributed, and (3) although no single algorithm clearly out performed the others, LDA, SVM and logistic regression methods worked well under normality. KNN and treebased methods, particularly random forest, also worked well under nonnormal distributions.
We also performed simulations to assess the impact of sample size and correlation structure on data dimension reduction algorithms. More specifically, we assessed how sample size and correlation structure impacted the number of components chosen out of the 12 (p = 12) simulated variables. Tables 1 presents the number of components chosen (out of p = 12) by LASSO method. As expected, as the correlation between the 12 variables increases, the algorithm selects fewer number of variables. The AR(1) type correlation structure results in more selected variables compared to the compound symmetry correlation structure. The results from t and Cauchy distributions are unstable for small sample sizes thereby yielding unexpected results. Increasing in sample size resulted in an increased number of components selected.
Data analyses
The study included 25 HIVinfected patients on NRTIbased ARTs with mitochondrial toxicity (cases), and 25 HIVinfected patients on NRTIbased ART without mitochondrial toxicity (positive controls). The median age of participants was 53 years (interquartile range (IQR), 50–57) with 60% of them being males. The race distribution among participants was 60% African American, 28% White, and 12% Hispanics. The median duration of mitochondrial toxicity was 2.2 years (IQR, 1.3–7.0). About 52% of the cases had only one manifestation of toxicity, while 48% had multiple toxicities.
A total of 12 ribose and 14 deoxyribose variables were considered in the analyses. Table 2 presents the 5fold Monte Carlo average classification rates with corresponding standard deviations in parentheses. As shown in Table 2, the algorithms resulted in crossvalidated classification rates of 0.54–0.76 for dNTP and 0.59–0.69 for RN. We have also reduced the dimension of the data using the LASSO method and applied classifiers on the reduced data. The LASSO method selects 10 deoxyribose and 2–4 ribose variables depending on whether lambda.min and lambda.1se was used as a selection criterion. For the deoxyribose data, both lambda.min and lambda.1se selects ten variables (dATP, dADP, dAMP, dTTP, dTDP, dTMP, dGTP, dGMP, dCDP, dCMP). However, for the Ribose data, lambda.min selects four variables (ATP, CDP, UTP, UMP) and lambda.1se selects two variables (ATP, UTP). Because the correlations between the deoxyribose input variables are lower than those of the ribose input variables, it is no surprise that LASSO selects most of deoxyribose input variables. Table 3 presents the 5fold Monte Carlo average classification rates with corresponding standard deviations in parentheses for the reduced data. In general, the reduction of input variables indeed seems to improve the classification performance, particularly for the ribose dataset. The crossvalidate classification rates for the ten selected dRN variables ranges from 0.65–0.76. However, crossvalidate classification rates using the selected two (i.e., ATP, UTP) or four (ATP, CDP, UTP, UMP) RN variables ranged from 0.72–0.83.
As for which machine learning method to recommend in a classification, the treebased methods work the best for both the deoxyribose dataset (RF) and the ribose dataset (CART and AdaBoost). But it is worthwhile noting that the simple methods such as LDA and logistic regression were very competitive in terms of classification performance. Therefore, we recommend the use of treebased methods for this application, but if one is not comfortable with technicalities of such methods, then the simple methods such as LDA and logistic regression work as well.
Discussion
Our study aimed to evaluate whether ARTinduced mitochondrial dysfunction assessed by RN and dNTP pool sizes can be used as biomarkers in distinguishing HIV patients with mitochondrial toxicity from those without toxicity. We used eight classification algorithms to assess the diagnostic performance of RN and dNTP pool sizes distinguishing HIV patients with and without NRTIassociated mitochondrial toxicity. The algorithms resulted in crossvalidated classification rates of 0.54–0.76 for dNTP and 0.59–0.69 for RN. dNTPs can be synthesized via two pathways: salvage and de novo pathways [23]. In the de novo pathway, RN is reduced to dNTP. Since there are two main sources of dNTPs, there may not be onetoone relationship between RN and dNTP pools. This could explain the low classification rates observed in our analyses.
The reduction of input variables improved the classification performance with most of the classification algorithms. The improvement was more pronounced for RN. Due to the high correlation among the RN variables, the data reduction technique results in fewer RN variables as compared to the 10 selected dNTP variables. The crossvalidate classification rates for the 10 selected dNTP variables ranges from 0.65–0.76. However, crossvalidate classification rates using the selected two (i.e., ATP, UTP) or four (ATP, CDP, UTP, UMP) RN variables ranged from 0.72–0.83. The concentration of the dNTPs available during DNA replication is critical for the fidelity of DNA replication [24, 25]. dNTPs are essential precursors for DNA synthesis and perturbations in the absolute and relative concentrations of the 4 dNTPs (dATP, dTTP, dGTP and dCTP) increase mutation rates by reducing the fidelity of DNA synthesis [23]. Furthermore, imbalance in dNTP pools has been associated with mitochondrial DNA (mtDNA) mutagenesis in cell culture and animal models [26, 27]. Maintaining dNTP pools in the cell is critical not only for DNA replication but also for cell cycle control, protooncogene function, mitochondrial function, defense against viral infections, DNA mismatch repair (MMR), telomere length, mitochondrial function [24, 28, 29]. Thus, our finding of changes in RN and dNTP pools in participants with mitochondrial toxicity validates the importance of dNTP pools in mitochondrial function. Hence, RN and dNTP pools can be used as biomarkers of ARTinduced mitochondrial toxicity.
In cells, the concentration of RN is several folds higher than the concentration of dNTP [30,31,32]. Because of the abundance of RN in the cells, its quantification is much easier than dNTP. Thus, ATP can be measured easily as biomarker of ARTinduced mitochondrial dysfunction. ATP production in the cell is predominantly from mitochondria through oxidative phosphorylation (OXPHOS). Thus, ATP concentration could serve as a sensitive marker of mitochondrial function. Further studies are needed to validate our finding that ATP concentration can serve as a biomarker of ARTinduced mitochondrial dysfunction. This could lead to the development of pointofcare assay for the diagnosis and monitoring of ARTinduced mitochondrial toxicity.
Currently, there is no goldstandard assay for diagnosis of ARTinduced mitochondrial toxicity. ARTinduced mitochondrial toxicity is diagnosed by a combination of clinical symptoms, laboratory testing, and imaging studies, and sometime a tissue biopsy to demonstrate mitochondrial damage [8, 9]. Tissue biopsy is deemed to be the most accurate method. However, because it is invasive and cost prohibitive, it is seldom used in clinical practice. Moreover, not all affected organs or tissues are easily accessible for biopsy. Therefore, most providers resort to stopping the perceived offending medication to see if the clinical manifestations resolve. This practice is not optimum as it can results in the emergence of drugresistant strains of HIV [33] and could lead to inappropriate use of secondline medications. There is a need for a noninvasive, costeffective biomarker for ARTinduced toxicity to prevent unnecessary interruptions in ART and to guide use of secondline regimens. Our finding that intracellular concentration of ATP determined by Liquid Chromatography with tandem mass spectrometry (LC/MS/MS) could be a biomarker of mitochondrial toxicity is promising. Mitochondria are responsible for ATP production through oxidative phosphorylation. Therefore, ARTinduced mitochondrial dysfunction is likely to compromise ATP synthesis capacity of mitochondria. Arguments against routine use of LC/MS/MS for measuring ATP as marker of ARTinduced mitochondrial toxicity are cost and laborintensive nature of the procedure. This can be circumvented by: (1) using available and easy to use fluorometric enzymelinked assay kits for quantifying intracellular ATP levels; or (2) using ATP concentration measured LC/MS/MS as gold standard and compare to available in vitro assays of mitochondrial function to determine which in vitro biomarkers best correlate with the level ATP measured by LC/MS/MS and, therefore, serve as the “goto” assay(s) for diagnosing and monitoring ARTinduced mitochondrial toxicity.
Conclusions
We used a range of machine learning procedures to distinguish between HIV patients with and without toxicity. Using resampling methods like Monte Carlo kfold cross validation, we compared the accuracy of several machine learning algorithms applied to our data. We used the algorithm with highest classification accuracy rate in evaluating the diagnostic performance of RN and dNTP pool sizes as biomarkers of mitochondrial toxicity. The algorithms resulted in crossvalidated classification rates of 0.65–0.76 for dNTP and 0.72–0.83 for RN. Our finding of changes in RN and dNTP pools in participants with mitochondrial toxicity validates the importance of dNTP pools in mitochondrial function. Hence, levels of RN and dNTP pools can be used as biomarkers of ARTinduced mitochondrial toxicity. There are none invasive and costeffective assays to measure intracellular ATP concentration that could be used to monitor or diagnose ARTinduced mitochondrial toxicity.
Availability of data and materials
The dataset used in the manuscript is available from the corresponding author on reasonable request.
Abbreviations
 ADA:

AdaBoost
 AIDS:

Acquired immunodeficiency syndrome
 AR (1):

First order autoregressive
 ART:

Antiretroviral therapy
 CART:

Classification and regression tree
 CV:

Cross validation
 DNA:

Deoxyribonucleic acid
 dNTPs:

deoxyribonucleotide triphosphates
 dNTTP:

deoxyribonucleotide
 HIV:

Human immunodeficiency virus
 IQR:

Interquartile range
 KNN:

Knearest neighbor
 LASSO:

Least absolute shrinkage and selection operator
 LDA:

Linear discriminant analysis
 LOGIT:

Logistic regression
 mtDNA:

mitochondrial DNA
 NRTIs:

Nucleoside reverse transcriptase inhibitors
 PCA:

Principal component analysis
 PIs:

protease inhibitors
 QDA:

Quadratic discriminant analysis
 RF:

Random forest
 RN:

Ribonucleotides
 SVM:

Support vector machine
References
 1.
Gardner K, Hall PA, Chinnery PF, Payne BAI. HIV treatment and associated mitochondrial pathology: review of 25 years of in vitro, animal, and human studies. Toxicol Pathol. 2014;42:811–22.
 2.
Di Biagio A, CozziLepri A, Prinapori R, Angarano G, Gori A, Quirino T, et al. Discontinuation of initial antiretroviral therapy in clinical practice: moving toward individualized therapy. J Acquir Immune Defic Syndr. 2016;71:263–71.
 3.
Brinkman K, ter Hofstede HJ, Burger DM, Smeitink JA, Koopmans PP. Adverse effects of reverse transcriptase inhibitors: mitochondrial toxicity as common pathway. AIDS. 1998;12:1735–44.
 4.
López S, Miró O, Martínez E, Pedrol E, RodríguezSantiago B, Milinkovic A, et al. Mitochondrial effects of antiretroviral therapies in asymptomatic patients. Antivir Ther (Lond). 2004;9:47–55.
 5.
Montaner JSG, Côté HCF, Harris M, Hogg RS, Yip B, Harrigan PR, et al. Nucleosiderelated mitochondrial toxicity among HIVinfected patients receiving antiretroviral therapy: insights from the evaluation of venous lactic acid and peripheral blood mitochondrial DNA. Clin Infect Dis. 2004;38(Suppl 2):S73–9.
 6.
Moyle G. Clinical manifestations and management of antiretroviral nucleoside analogrelated mitochondrial toxicity. Clin Ther. 2000;22:911–36 discussion 898.
 7.
Deng W, Baki L, Yin J, Zhou H, Baumgarten CM. HIV protease inhibitors elicit volumesensitive cl− current in cardiac myocytes via mitochondrial ROS. J Mol Cell Cardiol. 2010;49:746–52.
 8.
BlasGarcía A, Apostolova N, Ballesteros D, Monleón D, Morales JM, Rocha M, et al. Inhibition of mitochondrial function by efavirenz increases lipid content in hepatic cells. Hepatology. 2010;52:115–25.
 9.
Selvaraj S, Ghebremichael M, Li M, Foli Y, LangsBarlow A, Ogbuagu A, et al. Antiretroviral therapyinduced mitochondrial toxicity: potential mechanisms beyond polymeraseγ inhibition. Clin Pharmacol Ther. 2014;96:110–20.
 10.
Hastie T, Tibshirani R, Friedman J. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. 2nd edition. New York: SpringerVerlag; 2009. //www.springer.com/us/book/9780387848570. .
 11.
Fisher RA. The use of multiple measurements in taxonomic problems. Ann Eugenics. 1936;7:179–88.
 12.
Cover T, Hart P. Nearest neighbor pattern classification. IEEE Trans Inf Theory. 1967;13:21–7.
 13.
Hosmer D, Lemeshow S. Introduction to the Logistic Regression Model. In: Applied Logistic Regression. John Wiley & Sons, Ltd; 2005. p. 1–30. doi:https://doi.org/10.1002/0471722146.ch1.
 14.
Firth D. Bias reduction of maximum likelihood estimates. Biometrika. 1993;1:27–38.
 15.
Breiman L. Classification and regression trees. Routledge. 2017. https://doi.org/10.1201/9781315139470.
 16.
Freund Y, Schapire RE. A decisiontheoretic generalization of online learning and an application to boosting. J Comput Syst Sci. 1997;55:119–39.
 17.
James G, Witten D, Hastie T, Tibshirani R. An introduction to statistical learning with applications in R. New York: Springer; 2013.
 18.
Breiman L. Random Forests Machine Learning. 2001;45:5–32.
 19.
Rencher AC, Christensen WF. Methods of multivariate analysis. Third ed. Hoboken, New Jersey: Wiley; 2012.
 20.
Tibshirani R. Regression shrinkage and selection via the Lasso. J R Stat Soc Ser B Methodol. 1996;58:267–88.
 21.
Zou H, Hastie T. Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2005;67:301–20.
 22.
Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33:1–22.
 23.
Van Rompay AR, Johansson M, Karlsson A. Phosphorylation of nucleosides and nucleoside analogs by mammalian nucleoside monophosphate kinases. Pharmacol Ther. 2000;87:189–98.
 24.
Mathews CK, Song S. Maintaining precursor pools for mitochondrial DNA replication. FASEB J. 2007;21:2294–303.
 25.
Reichard P. Interactions between deoxyribonucleotide and DNA synthesis. Annu Rev Biochem. 1988;57:349–74.
 26.
Song S, Pursell ZF, Copeland WC, Longley MJ, Kunkel TA, Mathews CK. DNA precursor asymmetries in mammalian tissue mitochondria and possible contribution to mutagenesis through reduced replication fidelity. Proc Natl Acad Sci U S A. 2005;102:4990–5.
 27.
Song S, Wheeler LJ, Mathews CK. Deoxyribonucleotide pool imbalance stimulates deletions in HeLa cell mitochondrial DNA. J Biol Chem. 2003;278:43893–6.
 28.
Mathews CK. DNA precursor metabolism and genomic stability. FASEB J. 2006;20:1300–14.
 29.
Kunz BA, Kohalmi SE, Kunkel TA, Mathews CK, McIntosh EM, Reidy JA. International Commission for Protection against Environmental Mutagens and Carcinogens. Deoxyribonucleoside triphosphate levels: a critical factor in the maintenance of genetic stability. Mutat Res. 1994;318:1–64.
 30.
SmitskampWilms E, Pinedo HM, Veerman G. Ruiz van Haperen VW. Peters GJ Postconfluent multilayered cell line cultures for selective screening of gemcitabine Eur J Cancer. 1998;34:921–6.
 31.
Traut TW. Physiological concentrations of purines and pyrimidines. Mol Cell Biochem. 1994;140:1–22.
 32.
van Moorsel CJ, Bergman AM, Veerman G, Voorn DA. Ruiz van Haperen VW, Kroep JR, et al. differential effects of gemcitabine on ribonucleotide pools of twentyone solid tumour and leukaemia cell lines. Biochim Biophys Acta. 2000;1474:5–12.
 33.
ArnedoValero M, Garcia F, Gil C, Guila T, Fumero E, Castro P, et al. Risk of selecting de novo drugresistance mutations during structured treatment interruptions in patients with chronic HIV infection. Clin Infect Dis. 2005;41(6):883–90.
Acknowledgements
The authors would like to thank the patients who participated in the study and all the providers and nursing staff at the YaleNew Haven Hospital for making the study possible.
Funding
This analysis was supported by grants from the National Institutes of Health NIH/NIAID P30 AI060354. The parent study was supported by the National Institute of Allergy and Infectious Diseases at the National Institutes of Health [K08AI074404]. The funding body had no role in the design of the study; the collection, analyses, and interpretation of data, and in writing manuscript.
Consent for publication
Not applicable.
Author information
Affiliations
Contributions
M.G, J.L and V.G analyzed data and wrote manuscript. E.P. wrote manuscript and provided guidance on the clinical interpretation of the findings. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Musie Ghebremichael.
Ethics declarations
Ethics approval and consent to participate
The study protocol was approved by the Institutional Review Board of the Yale University School of Medicine. Informed and written consents were obtained from study participations prior to enrollment in the study.
Competing interests
The authors do not have a commercial or other association that might pose a conflict of interest, i.e., authors declare that they have no competing interests to disclose.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Lee, J.S., Paintsil, E., Gopalakrishnan, V. et al. A comparison of machine learning techniques for classification of HIV patients with antiretroviral therapyinduced mitochondrial toxicity from those without mitochondrial toxicity. BMC Med Res Methodol 19, 216 (2019) doi:10.1186/s128740190848z
Received:
Accepted:
Published:
Keywords
 HIV/AIDS
 Machine learning
 Classification
 Dimension reduction
 Mitochondrial toxicity
 Antiretroviral therapy