Skip to main content

Risk score prediction model based on single nucleotide polymorphism for predicting malaria: a machine learning approach

Abstract

Background

The malaria risk prediction is currently limited to using advanced statistical methods, such as time series and cluster analysis on epidemiological data. Nevertheless, machine learning models have been explored to study the complexity of malaria through blood smear images and environmental data. However, to the best of our knowledge, no study analyses the contribution of Single Nucleotide Polymorphisms (SNPs) to malaria using a machine learning model. More specifically, this study aims to quantify an individual's susceptibility to the development of malaria by using risk scores obtained from the cumulative effects of SNPs, known as weighted genetic risk scores (wGRS).

Results

We proposed an SNP-based feature extraction algorithm that incorporates the susceptibility information of an individual to malaria to generate the feature set. However, it can become computationally expensive for a machine learning model to learn from many SNPs. Therefore, we reduced the feature set by employing the Logistic Regression and Recursive Feature Elimination (LR-RFE) method to select SNPs that improve the efficacy of our model. Next, we calculated the wGRS of the selected feature set, which is used as the model's target variables. Moreover, to compare the performance of the wGRS-only model, we calculated and evaluated the combination of wGRS with genotype frequency (wGRS + GF). Finally, Light Gradient Boosting Machine (LightGBM), eXtreme Gradient Boosting (XGBoost), and Ridge regression algorithms are utilized to establish the machine learning models for malaria risk prediction.

Conclusions

Our proposed approach identified SNP rs334 as the most contributing feature with an importance score of 6.224 compared to the baseline, with an importance score of 1.1314. This is an important result as prior studies have proven that rs334 is a major genetic risk factor for malaria. The analysis and comparison of the three machine learning models demonstrated that LightGBM achieves the highest model performance with a Mean Absolute Error (MAE) score of 0.0373. Furthermore, based on wGRS + GF, all models performed significantly better than wGRS alone, in which LightGBM obtained the best performance (0.0033 MAE score).

Peer Review reports

Background

Malaria is a mosquito-borne infectious disease that can progressively infect and annihilate the red blood cells. The World Health Organization’s report entitled “World Malaria Report” [1] estimated the annual death of 400,000 from malaria, of which two-thirds are children below the age of five. It is also more commonly diagnosed in sub-Saharan Africa, where a significantly large number of cases throughout the endemic region are due to increased mosquitoes' reproduction resulting in more malaria transmission. Furthermore, ineffective prevention methods are still widespread [2]. Thus, much research has been devoted to developing tools for malaria surveillance. It is believed that prevention is better than cure [3], and hence, the implementations of preventative solutions are crucial in determining individuals’ early detection even before the disease symptoms strike.

The development of malaria cases in malaria-endemic regions is not the only known risk factor. According to other observations in the literature [4,5,6], specific inherited genetic traits are also possible risks of malaria contraction among different populations. Thus, genetic susceptibility (an inherited increase in the risk of developing the disease) [7] and genetic resistance (the ability to prevent the risk of developing the disease or a marked reduction in the severity of the symptoms) to malaria have been characterized by genetic mutations with erythrocytes. This includes hemoglobin variants or related sickle cell disease, Duffy blood antigens glycophorins, and blood type groupings. Substantial evidence in the past decade has established that different populations have different susceptibility to malaria due to diverse genetic adaptations and gene selection pressures of malaria in the human genome [8].

A genome-wide association study (GWAS) is an approach used in genetics research that scans genomes of multiple individuals to find specific genetic variations with a particular disease. These genetic variations are markers or contributing risk factors of the disease [9]. Multiple prior studies [10, 11] have established that the Hemoglobin Subunit Beta (HBB) gene to be a significant genetic risk factor for malaria because of its risk alleles. However, HBB variants alone cannot precisely predict the disease as malaria is a complex disease. Many other genes, such as ABO, ATP2B4, G6PD, CD40LG, FY, GYPA, GYPB, GYPC, HBA, HP, SCL4A1, have been associated with malaria susceptibility or resistance across different populations via GWAS [12]. To further combat the disease's spread, it is crucial to understand the contributing genetic risk factors involved in characterizing the disease's complexity [10].

Machine learning is an Artificial Intelligence strategy that uses training algorithms to learn from large datasets, and is the basis of numerous developments, from speech recognition to autonomous self-driving cars [13]. Inspired by the recent breakthroughs of machine learning in various domains, this paper applies machine learning algorithms to study malaria’s complexity by leveraging well-planned malaria research on genetic variants. The Malaria Genomic Epidemiology Network (MalariaGEN) http://www.malariagen.net/ is a community of researchers from more than 20 countries to understand how natural genetic variation in humans and malaria parasites affects the biology and epidemiology of malaria. This network, formally established in 2005, has developed several GWAS [12, 14,15,16,17,18,19,20,21,22] with their collaborating partners to fulfil their goals. Therefore, MalariaGEN produces the essential data required for machine learning algorithms to mine for valuable patterns.

In recent endeavours, machine learning algorithms were used to explore the complexity of malaria, particularly malaria parasites and development stages, through blood smear images [23, 24]. In separate studies, environmental data were collected and trained by machine learning algorithms to link climate change to malaria transmission [25, 26]. Since the aforementioned prediction strategy was deemed successful, we leveraged the power of machine learning to quantify the risk score of an individual’s susceptibility to malaria based on genetic variation, i.e., the contributing risk factors, instead. Currently, the prediction is limited to using advanced statistical methods, such as time series and cluster analysis [27]. Hence, machine learning prediction models based on genetic variation are required to fully explore the disease's potential genetic markers. The machine learning algorithms can detect valuable patterns in complex datasets by applying various optimization, statistical and probabilistic methods to identify the most suitable set of variables to train the model, for example, Single Nucleotide Polymorphisms (SNPs). SNP is the most common type of genetic variation among individuals. We are interested in finding out whether it is possible to quantify an individual’s risk of malaria based on SNP genotype data for facilitating personalized prevention and treatment. Thus, instead of classifying “Is this individual infected by malaria?”, this study predicts a risk score to answer the question “What is the individual risk score towards malaria susceptibility?”. In answering this question, we propose a feature extraction algorithm that aids in selecting the main SNP features representing the genetic risk factors as well as exploring machine learning models’ performance based on SNP genotype data. Weighted genetic risk scores (wGRS) and the combination of wGRS with genotype frequency (wGRS + GF) were calculated as the target variable. The contributions of this paper are summarized below:

  • Proposes a novel feature extraction algorithm using genotype patterns that can aid Logistic Regression and Recursive Feature Elimination (LR-RFE) in selecting the significant features. We also provide the ranked features list obtained from our approach;

  • Proposes a novel machine learning-based model to predict the risk score of an individual developing malaria, and;

  • Proposes a risk score that combines wGRS and genotype frequency (wGRS + GF). We also provide a comprehensive analysis of the experimental results, comparing wGRS with wGRS + GF.

This study hypothesizes that we can predict the risk score of an individual’s susceptibility to malaria from genetic variants through the proposed feature extraction algorithm with LR-RFE and the analysis of multiple machine learning models. The potential findings will add to our understanding of this topic by exploring genetic variation for different populations to predict an individual’s risk score of developing malaria, i.e., individual malaria risk score, in various regions. This will further create a comprehensive understanding of malaria susceptibility and reaching out to a wider community to combat the spread of malaria.

Methods

Figure 1 depicts the pipeline for identifying the malaria risk score in subjects, starting from the raw data until the evaluation stage of the prediction model. The pipeline consists of four discrete stages of operation: (1) Data Mining and Modeling, (2) Feature Extraction and Selection, (3) Model Development and (4) Model Evaluation. Figure 2 shows the methodology flow chart in detail.

Fig. 1
figure 1

Machine learning pipeline for individual malaria risk score prediction

Fig. 2
figure 2

Methodology flow chart

Dataset

Our study uses the human GWAS data produced from the MalariaGEN Consortial Project 1 entitled: “Genome-wide study of resistance to severe malaria in eleven populations.” The consortial project structure is described in [28], and information on each collaborating partner's studies and field sites is acknowledged on the MalariaGEN website.

We used genotype data of 20,854 individuals (10,791 malaria-affected individuals and 10,063 controls) from 11 worldwide populations (Table 1).

Table 1 Analysed populations and samples

Data mining and modeling

Data preprocessing

Initially, the first stage of the pipeline is converting raw MalariaGEN genotype data into a format that machine learning models can use. We downloaded the data from the European Genome-phenome Archive http://ega-archive.org/. We then identify 122 SNPs relevant to malaria through literature reviews, where 31 academic articles [8, 12, 15,16,17,18,19,20,21,22, 29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49] are reviewed and analysed. Out of 122 SNPs, 18 SNPs are removed due to unreported effect size and unavailability in some populations. A total of 104 SNPs is retained (Additional file 1), with genotypes comprising major allele A and minor allele a. This stage also converts all unparseable values in the data, such as data types and standard format errors, into null representations. We also map 32 kgpIDs to rsIDs and remove 37 samples without detailed information about malaria subtypes. This preprocessing procedure yielded 104 SNP variables from a total of 20,817 samples.

Generally, genotype imputation softwares such as IMPUTE2 [50] and Beagle [51] are deployed for estimating missing genotypes. These software programs impute missing genotypes based on publicly available reference datasets such as 1000 Genomes Project or HapMap 3. However, in our case, imputation needs to be more specific as we are developing a prediction model of individual’s susceptibilities to malaria. Thus, we developed a python program that imputes any missing genotypes based on the population group and malaria subtype from the human GWAS data used in this study. In order to do so, the program first groups individuals based on their countries and then by their malaria subtypes. Finally, a comparison of a total of six SNPs for each missing genotype, i.e., three SNPs before and after the missing loci, is performed before imputing the missing genotype with the most common genotype data.

Moreover, to perform feature extraction and selection separately from model development, we split the preprocessed dataset into two independent datasets, of which 25% were used for feature extraction and selection (5,204 samples; 2,667 malaria-affected individuals and 2,537 controls), and 75% were used for model development (15,613 samples; 8,087 malaria-affected individuals and 7,526 controls).

Dataset for feature extraction and selection

As a preparation for the feature extraction and selection stage, we represent the genotype data as feature and class label data frame. Each SNP itself is a feature, where the columns represent 104 SNPs containing genotype data. The genotype data is formed by two alleles, A and a, generally expressed as AA, Aa, and aa. The last column representing the class label contains the categorical value representing the binary classification of the individuals: 0–case (malaria-affected) and 1–control (healthy). Note that binary classification is only used for feature extraction and selection; for model development, the target variable is the risk score described in the next section.

Dataset for model development

The model is designed to provide individuals with continuous risk scores rather than binary classification. Therefore, we convert the categorical value in the class label by calculating wGRS as the target variable. To further evaluate the model performance of wGRS, we also computed wGRS + GF as the target variable for an identical dataset, which will be compared with wGRS-only model to determine the better performing model.

Genotype patterns and frequency

Several observations in the literature [52, 53] indicate that genotype patterns contribute to disease association. The exploration of genotype patterns is significant for malaria prediction as substantial evidence has established that sickle cell anemia traits can partially prevent malaria [10, 54,55,56]. The genetic trait of sickle cell anemia is found on the recessive allele of the hemoglobin gene. This means that an individual needs to have two copies of the recessive alleles—one from the mother and one from the father—to have this condition. If the alleles are heterozygous, the individual tends to be resistant to the development of malaria. In contrast, if the alleles are homozygous, the individual is susceptible to the development of malaria. Inspired by these evidences, our study will include genotype frequency in the formulation of wGRS + GF to characterize the individual’s malaria risk score; and propose a novel feature extraction algorithm based on genotype patterns and frequencies.

wGRS

The wGRS is calculated for each genotype by multiplying the number of risk alleles (0, 1, 2) by the estimated effect size reported for each variant. The reported effect size was estimated using the logistic regression association tests method and can be found in the association test summary statistics on the MalariaGEN website https://www.malariagen.net/sppl25/. The calculation of wGRS is summarized as: risk allele * effect size = wGRS [57].

wGRS + GF

The wGRS + GF combines the wGRS mentioned above and genotype frequency. Similarly, the wGRS + GF is calculated for each genotype. The Hardy–Weinberg equation is utilized to calculate genotype frequency from genotype data as this equation calculates an individual’s genetic variation at equilibrium. For example, the frequency of allele A is represented by p, and the frequency of allele a is represented by q. Thus, the frequency of genotype AA is calculated as p2, genotype aa is calculated as q2, and genotype Aa is calculated as 2pq. The wGRS + GF is then produced by multiplying the genotype frequency by the wGRS. The calculation of wGRS + GF is summarized as: (risk allele * effect size) * genotype frequency = wGRS + GF. Of note, this approach is a novelty, because researchers to date only used wGRS [58,59,60]. However, the wGRS only considers the risk alleles and the variant effect size, which is insufficient in these aspects for malaria prediction as genotype patterns have been associated with malaria. Thus, it is essential to include genotype frequency to characterize an individual’s malaria risk score further.

Feature extraction and selection

Feature extraction

We proposed an algorithm to extract genotype, i.e., SNP features, into a valuable set of information that aids feature selection. As mentioned earlier, for all subjects, genotype data is usually represented as AA, Aa, and aa, by two alleles A (major allele) and a (minor allele). Such representation leads to repeated data that may affect analysis performance, and thus feature extraction and normalization are required. The proposed algorithm focuses on genotype patterns where we find pattern frequencies based on populations due to its strong linkage between genotype patterns and malaria. This approach is a novelty because the current state-of-the-art methods substitute genotype data with numerical values {0, 1} for the major and minor allele [61]; or with numerical values {0, 1, 2} for homozygous major, heterozygous and homozygous minor [62,63,64,65]. Simultaneously, some other studies rely on odd ratios calculated from genotype data [60], and some are combined with clinical features [66].

For each population, the proposed feature extraction algorithm makes a left-to-right pass over the preprocessed dataset to compute the frequency counts of the genotype patterns. The genotype patterns and frequency counts are stored in a dictionary. Next, we compute the total of all the genotype patterns by summing up all the counts in the dictionary. The total is used for data normalization. Another left-to-right pass is then made over each population's preprocessed dataset, but this time, we retrieve the genotype pattern counts from the dictionary in O(1) time to divide it by the total, which we define the resulting value as the pattern frequency. Each pattern frequency is stored based on the SNP feature list. The pseudocode of the proposed feature extraction algorithm is presented in Fig. 3. The most common feature extraction approach that substitutes homozygous major, heterozygous and homozygous minor with numerical values {0, 1, 2} is used for benchmark comparisons.

Fig. 3
figure 3

Pseudocode of the proposed feature extraction algorithm

Feature selection

To create an accurate model that relies on the most relevant features (namely, features that contribute the most for malaria prediction), we evaluated the feature importance, i.e., feature dependence of the model, through Logistic Regression and Recursive Feature Elimination (LR-RFE). The analysis is based on the feature extraction algorithm described above, followed by the Recursive Feature Elimination (RFE) algorithm, in which the Logistic Regression (LR) coefficient is utilized in the core of the model to perform feature selection.

The LR model’s coefficients have been widely utilized for feature importance estimation [67]. Each coefficient represents a score, known as the feature importance score, which describes the significance level between the feature and the target variable. The higher the coefficient, the more relevant the feature is to the target variable. In other words, coefficients can be utilized to determine the important and unimportant features to avoid overfitting [67] and are thus useful for prediction [68]. The RFE model ranks the 104 features based on their importance scores obtained from the LR model into a list, in which the first position represents the most significant feature, while the least important feature is ranked on the last position. The features are ranked via their importance scores, where the least important features are iteratively eliminated through remodeling until the required number of features is retained. This iterative process results in a list of ranked features, i.e., from the most to the least important features.

In this regard, a recent study [69] used LR-RFE to rank feature importance and selection to find the optimal feature set in breast cancer prediction. The study focused on cytological characteristics obtained from the breast fine needle aspiration test. Several models were developed to predict breast cancer using different number of features. It was concluded that LR-RFE contributed to better classification performance and improved model accuracy. Therefore, we chose to use LR-RFE on genetic variation to observe the effectiveness of improving the accuracy of risk score prediction.

Subsequently, to explore different sets of features for predicting malaria risk score, the results obtained are then utilized to generate 104 feature sets based on feature importance ranking. For example, the first feature set consists of the top one feature, whereas the second feature set consists of the top two features, and so on as shown in Fig. 4. Finally, the last feature set consisting of all 104 features is used to compare the efficacy of other feature sets. For each feature set, the cumulative effects of features are used as the target variable, whereby the wGRS and wGRS + GF are summed, respectively. Figure 5 shows the high-level pseudocode up to this stage.

Fig. 4
figure 4

Overview of genotype-pattern-frequency-based features

Fig. 5
figure 5

High-level pseudocode of the feature extraction and selection stage

Model development

Preprocessed datasets from the previous stage are split into two independent datasets, namely Test-1 (15%) and Test-2 (85%). The study in [70] recommended the split percentage.

Test-1 and Test-2 are split into seven equally-sized random groups by using sevenfold cross-validation to prevent overfitting. Three machine learning models are trained, namely, Light Gradient Boosting Machine (LightGBM), eXtreme Gradient Boosting (XGBoost), and Ridge regression. We choose these machine learning algorithms because they are rarely used to analyse SNP genotype data and have been proven to have higher efficiency and faster training speed in other machine learning domains. For each machine learning model, the Test-1 dataset is utilized to obtain the best parameters through a grid search. These best parameters are then used on the Test-2 dataset for the machine learning prediction models.

Model evaluation

Performance metrics

The final stage of the pipeline in Fig. 1 evaluates the performance of each model in predicting individual malaria risk scores. The outcome is a continuous risk score, so the Mean Absolute Error (MAE) and Root Mean Squared Error (RMSE) metrics are utilized [71]. The MAE measures the average residual error between the target and predicted values. On the other hand, the RMSE measures the square root of the average squared residual error between the target and predicted values. Unlike other metrics, MAE and RMSE are negatively-oriented scores, where a smaller value indicates better model performance. In other words, the lower the MAE or RMSE value, the higher the prediction accuracy.

All code was developed using the Python programming language, and simulations were performed on a machine with Apple M1 Max processor and 32 GB of memory. All methods were carried out in accordance with relevant guidelines and regulations.

Results

Optimization of feature importance ranking

Figures 6 and 7 highlight the feature importance ranking and importance scores for all 104 SNPs, starting from the most important to the least important in the context of malaria risk score prediction. For readability purposes, we only show the absolute value of the scores. The ranks in Fig. 6 are computed using the proposed feature extraction algorithm with LR-RFE as detailed in the Methods section. (We also provide a list of these features ranked in ascending order in the Additional file 2.) By contrast, Fig. 7 is based on the feature extraction algorithm that substitutes homozygous major, heterozygous, and homozygous minor genotypes with numerical values {0, 1, 2} with LR-RFE, as benchmark comparisons.

Fig. 6
figure 6

Feature importance ranking of all 104 SNPs, computed using the proposed feature extraction algorithm with LR-RFE

Fig. 7
figure 7

Feature importance ranking of all 104 SNPs, computed using the benchmark feature extraction algorithm with LR-RFE

The results of two different feature extraction algorithms with LR-RFE are summarized in Fig. 8, which shows the comparative feature importance of all 104 SNPs. These results indicate that LR-RFE ranks features and calculates their importance scores differently depending on the feature extraction algorithm in use. For example, rs334 was ranked as the 1st feature with an importance score of 6.224 in Fig. 6. However, when comparing with the benchmark in Fig. 7, it was ranked as the 1st feature with an importance score of 1.1314 instead. The higher the score, the more prominent the feature is in predicting malaria risk score. Of note, previous findings from MalariaGEN [12, 15, 18,19,20,21,22] and prior studies [10, 11] have proven that the HBB gene is the major genetic risk factor for malaria, and rs334 is an SNP from HBB. Thus, the proposed algorithm with LR-RFE can be considered a promising method, where it can identify rs334 with a higher importance score compared to the benchmark algorithm.

Fig. 8
figure 8

Comparison of feature importance scores using different feature extraction algorithms with LR-RFE: (1) proposed algorithm and (2) benchmark algorithm

Analysis of model prediction results

We train and evaluate three machine learning regression models to predict individual risk scores towards malaria susceptibility. Results across different feature sets based on wGRS and wGRS + GF as target variables are reported in Additional files 3 and 4, respectively. These results show that reducing the feature set sizes lowers the RMSE scores in contrast to increasing the sizes, which lowers the MAE scores. This is as expected because, in RMSE, the errors are squared before being averaged, which gives higher errors more weight resulting in the metric being more sensitive to isolated outliers than MAE [71]. In other words, increasing the feature set size increases the number of outliers, which in turn increases the number of errors. Hence, MAE is a better indicator of the model's performance than the RMSE and is used in this study as the evaluation metric.

When based on wGRS alone, the best performing model is LightGBM, which achieves an MAE score of 0.0373 when trained on a single feature, i.e., rs334, compared to the MAE score of 1.1104 when trained on all the 104 features. In contrast, based on wGRS + GF as the target variable, the MAE scores indicate that all models performed significantly better than wGRS alone. LightGBM obtains the best performance, yielding an MAE score of 0.0033 when only a single feature, i.e., rs334, is utilized. Compared to wGRS, the MAE scores obtained across different feature sets based on wGRS + GF as the target variable is much lower, indicating higher accuracy.

Regardless of whether wGRS or wGRS + GF are utilized, all models achieve the lowest MAE scores when trained solely on a feature, i.e., rs334. On the other hand, they achieve the highest MAE scores when trained with all the 104 features. Furthermore, there is a marginal difference between the scores obtained by using the default parameters and the best parameters. We also provide a graph (Fig. 9) containing the MAE scores against the number of features based on wGRS and wGRS + GF.

Fig. 9
figure 9

Performance analysis of the wGRS-based and wGRS + GF-based models with respect to MAE scores and feature sets

Discussion

The genetic-based method of predicting the malaria risk is a powerful and feasible method that can strengthen the existing prevention strategies for malaria elimination. Traditional GWAS has achieved fruitful results in exploring the genetic risk associated with malaria through SNPs identification in specific populations [72]. GWAS is conducted to understand the disease, genes, and pathways, but it is not designed to predict whether an individual may develop the disease. On the other hand, machine learning methods can identify interacting genetic variants [61, 73] and, thus, are used in predicting complex genetic diseases, such as breast cancer [64], asthma [74], and Crohn’s disease [75].

This study puts forward a hypothesis that machine learning methods can quantify individual risk scores of susceptibility to malaria based on genetic variants. We propose a genotype-pattern-frequency-based feature extraction algorithm with LR-RFE for feature selection, where the importance of genotype patterns in malaria prediction was essentially highlighted in the Methods section. To evaluate the performance of our proposed method, we implemented the most common feature extraction algorithm with LR-RFE for feature selection as a benchmark comparison. We observe that LR-RFE ranks features and calculates their importance scores differently depending on different feature extraction algorithms. Thus, a suitable feature extraction algorithm may be an essential step to discover the most important markers in building an accurate prediction model. In particular, when using the proposed feature extraction algorithm, LR-RFE ranked rs334 as the most contributing feature with an importance score of 6.224, while the baseline has a score of 1.1314. The importance score is crucial as the scores highlight the significant level of features contributing to malaria risk prediction that will affect the final prediction result. Therefore, the higher the score, the more prominent the feature is in predicting malaria risk.

This is an important result because previous findings from MalariaGEN [12, 15, 18,19,20,21,22] and prior studies [10, 11] have proven that the HBB gene is the major genetic risk factor for malaria, and rs334 is an SNP from HBB. Hence, the proposed feature extraction algorithm with LR-RFE appears to be a promising method to extract the significant contributing risk factors to malaria.

We calculate wGRS and wGRS + GF as target variables to learn and model the relationship between SNPs. These scores enable us to understand the degree to which a genetic marker, i.e., an SNP, is associated with malaria via machine learning. This includes training and analysing three machine learning regression models: LightGBM, XGBoost, and Ridge regression, to predict the risk score of an individual’s susceptibility to malaria through genetic variants information only. Machine learning models are expected to simulate the relationship between SNP genotype data and target variables, explaining more genetic components of complex diseases such as malaria. Compared with GWAS research alone, building a model to predict malaria risk scores among different populations will be an essential advancement for disease prediction [76].

The wGRS based results indicate that LightGBM is the best performing model obtaining an MAE score of 0.0373 when trained solely on feature rs334 in contrast to obtaining an MAE score of 1.1104 despite being trained on all 104 features. We further compare the performance of wGRS with wGRS + GF to identify the better performing model. Surprisingly, the MAE scores obtained across different feature sets for all models is much lower than solely utilizing wGRS, indicating significantly better performance. The best performing model is LightGBM, which achieves an MAE score of 0.0033 when trained solely on feature rs334. The wGRS + GF approach is a novelty, as to date, researchers have only used wGRS [58,59,60], which is insufficient for malaria risk score prediction as genotype patterns are essential in malaria prediction. It is also worth noting that when the number of features utilized to train the model increases, the data correlation becomes complex. Thus, the performance of the model decreases. Furthermore, as the model's MAE score is very low, these findings support our hypothesis that genetic variants are efficient markers of the disease and, therefore, may be used for future machine learning predictions.

We also note that both wGRS and wGRS + GF approaches achieve the best performing models when solely trained on feature rs334. This result is as expected, as multiple prior studies [10,11,12, 15, 18,19,20,21,22] have identified SNP rs334 from the HBB gene as malaria's main genetic risk factor. However, it is insufficient to train on a single feature as recent studies have shown that each SNP is associated with disease development only to an extent, and complex interactions between features may improve the predictive ability of the model [10, 77]. Moreover, many genes are involved in the genetic basis of malaria's susceptibility or resistance, making prediction complicated on several levels [10]. For instance, even though the rs334 feature may provide the most potent predictive power, the remaining features may improve the prediction further.

In addition, we compare the prediction performance of different feature combinations determined using the feature extraction algorithms with LR-RFE. To represent small and large feature sets respectively, we used sets 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100. The results are summarized in Table 2. We note that the feature set generated using the proposed feature extraction algorithm with LR-RFE achieved slightly better MAE scores when the feature sizes were between 50 and 70. In contrast, the feature set generated using the baseline with LR-RFE obtained slightly better MAE scores when the feature sizes were between 20 and 40. However, there is only a marginal difference in the MAE scores obtained for feature sizes 10, 80, 90, and 100.

Table 2 Comparison of prediction performance between different feature combinations determined using the feature extraction algorithms with LR-RFE

Our findings further reaffirm that prediction performance is affected by the complex interactions between the susceptibility or resistance features (i.e., SNP) in predicting malaria risk. This is because each SNP is not equally important due to its susceptibility or resistance levels. Therefore, when combined, the impact of the SNPs interactions differs from one feature set to another. Hence, we hardly justify which feature combination is the most important for malaria prediction as each feature has a particular contribution to malaria development.

Extremely statistically significant differences in MAE scores, i.e., p-value < 0.0001, are obtained using wGRS and wGRS + GF via one-way Analysis of Variance (ANOVA) test. Of note, we performed a normality test for the data distribution to ensure no violation of the one-way ANOVA test's normality assumption. We used the skewness and kurtosis index to identify the normality of the data, obtaining the range of [-0.4204, + 0.4005] for skewness and [-0.8025, -0.2566] for kurtosis. The data is considered normal if the skewness falls in the range of [‐2, + 2] [78, 79] and the kurtosis falls in the range of [‐7, + 7] [78, 79]. Hence, we can conclude the MAE scores of the models are normally distributed. To further compare the two methods, we also provide Table 3 summarizing the p-values obtained with the best parameters and Table 4 summarizing the p-values obtained with the default parameters.

Table 3 P-values obtained with best parameters
Table 4 P-values obtained with default parameters

This study has obtained promising results in predicting the risk score of individual susceptibility to malaria. We believe that our findings hold great promise for individual malaria risk score prediction and contribute in bridging the implementation gap between healthcare practitioners and computer scientists. The pipeline of the model’s development emphasized in this study can fully be reproduced and thus, be used as a base towards model retraining with new data.

Our study, based on 20,817 individuals from eleven populations, provides a basis for further exploration and improvement of the machine learning models, where individuals from different continents can be included. In addition, to explore more comprehensive solutions and strengthen our proposed method, alternate hypotheses may be considered, namely, developing a machine learning model that integrates blood group and SNP genotype data of malaria-affected individuals to characterize the genetic component and complexity of the disease. There is substantial evidence showing that blood group A is very susceptible to malaria, whereas blood group O can prevent malaria [80, 81].

Conclusions

We developed machine learning-based prediction models which utilized SNPs genotype data and calculated wGRS and wGRS + GF as target variables to quantify the risk score of an individual’s susceptibility to malaria. More precisely, we employed an approach that consists of the proposed genotype-pattern-frequency-based feature extraction algorithm with LR-RFE to identify the SNPs’ significant level implication to malaria. Results show that this approach identified and ranked SNP rs334 (a major genetic risk factor for malaria proved by previous studies) as the most contributing feature with an importance score of 6.224 compared to the baseline which only yields an importance score of 1.1314. (This suggests that our approach has the potential to discover significant genetic markers for other diseases as well.) Furthermore, LightGBM, a tree-based model, is the best-performing model in this study. It is also found that compared with wGRS alone, the model trained based on SNP genotype data and wGRS + GF obtains a lower MAE scores, and it is also a novelty as far as the literature on risk scores is concerned. To a larger extent, we have shown a promising method that demonstrates how machine learning can augment the insights derived from GWAS to quantify an individual’s risk score for a particular disease.

Availability of data and materials

The datasets analysed during the current study are available in the MalariaGEN Consortial Project 1 (https://www.malariagen.net/data/genome-wide-study-resistance-severe-malaria-eleven-populations-version-2); EGAD00010001799, EGAD00010001738, EGAD00010001818, EGAD00010001737, EGAD00010001739, EGAD00010001734, EGAD00010001736, EGAD00010001740, EGAD00010001743, EGAD00010001741, EGAD00010001742, EGAD00010001733, EGAD00010001735.

Abbreviations

ANOVA:

Analysis of variance

GWAS:

Genome-wide association study

HBB:

Hemoglobin subunit beta

LightGBM:

Light gradient boosting machine

LR:

Logistic regression

LR-RFE:

Logistic regression and recursive feature elimination

MAE:

Mean absolute error

MalariaGEN:

Malaria genomic epidemiology network

RFE:

Recursive feature elimination

RMSE:

Root mean squared error

SNP:

Single nucleotide polymorphism

wGRS:

Weighted genetic risk score

wGRS + GF:

Weighted genetic risk score with genotype frequency

XGBoost:

EXtreme gradient boosting

References

  1. World Health Organization. World malaria report 2020: 20 years of global progress and challenges. World Health Organization; 2020. Available from: https://www.who.int/docs/default-source/malaria/world-malaria-reports/9789240015791-double-page-view.pdf?sfvrsn=2c24349d_5.

  2. Childs LM, Cai FY, Kakani EG, Mitchell SN, Paton D, Gabrieli P, et al. Disrupting mosquito reproduction and parasite development for malaria control. PLoS Pathog. 2016;12(12): e1006060.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Tizifa TA, Kabaghe AN, McCann RS, Van den Berg H, Van Vugt M, Phiri KS. Prevention efforts for malaria. Curr Trop Med Rep. 2018;5(1):41–50.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Fortin A, Stevenson MM, Gros P. Susceptibility to malaria as a complex trait: big pressure from a tiny creature. Hum Mol Genet. 2002;11(20):2469–78.

    Article  CAS  PubMed  Google Scholar 

  5. Mackinnon MJ, Mwangi TW, Snow RW, Marsh K, Williams TN. Heritability of malaria in Africa. PLoS Med. 2005;2(12): e340.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Driss A, Hibbert JM, Wilson NO, Iqbal SA, Adamkiewicz TV, Stiles JK. Genetic polymorphisms linked to susceptibility to malaria. Malar J. 2011;10(1):271.

    Article  PubMed  PubMed Central  Google Scholar 

  7. National Cancer Institute. NCI Dictionary of Cancer Terms. National Cancer Institute; 2021. Available from: https://www.cancer.gov/publications/dictionaries/cancer-terms/def/genetic-susceptibility.

  8. De Mendonça VR, Goncalves MS, Barral-Netto M. The host genetic diversity in malaria infection. J Trop Med. 2012;2012: 940616.

    Article  PubMed  PubMed Central  Google Scholar 

  9. National Human Genome Research Institute. Genome-Wide Association Studies (GWAS). National Human Genome Research Institute; 2021. Available from: https://www.genome.gov/genetics-glossary/Genome-Wide-Association-Studies.

  10. Kwiatkowski DP. How malaria has affected the human genome and what human genetics can teach us about malaria. Am J Hum Genet. 2005;77(2):171–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Gouagna LC, Bancone G, Yao F, Yameogo B, Dabiré KR, Costantini C, et al. Genetic variation in human HBB is associated with Plasmodium falciparum transmission. Nat Genet. 2010;42(4):328–31.

    Article  CAS  PubMed  Google Scholar 

  12. Rockett KA, Clarke GM, Fitzpatrick K, Hubbart C, Jeffreys AE, Malaria Genomic Epidemiology Network, et al. Reappraisal of known malaria resistance loci in a large multicenter study. Nat Genet. 2014;46(11):1197–204.

    Article  PubMed Central  CAS  Google Scholar 

  13. Qiu J, Wu Q, Ding G, Xu Y, Feng S. A survey of machine learning for big data processing. EURASIP J Adv Signal Process. 2016;2016(1):1–6.

    Google Scholar 

  14. Network MGE. A global network for investigating the genomic epidemiology of malaria. Nature. 2008;456(7223):732–7.

    Article  CAS  Google Scholar 

  15. Ndila CM, Uyoga S, Macharia AW, Nyutu G, Peshu N, Ojal J, et al. Human candidate gene polymorphisms and risk of severe malaria in children in Kilifi, Kenya: a case-control association study. Lancet Haematol. 2018;5(8):e333–45.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Busby GB, Band G, Le QS, Jallow M, Bougama E, Mangano VD, et al. Admixture into and within sub-Saharan Africa. Elife. 2016;5: e15266.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Network MGE. A novel locus of resistance to severe malaria in a region of ancient balancing selection. Nature. 2015;526(7572):253–7.

    Article  CAS  Google Scholar 

  18. Shelton JM, Corran P, Risley P, Silva N, Hubbart C, Jeffreys A, et al. Genetic determinants of anti-malarial acquired immunity in a large multi-centre study. Malar J. 2015;14(1):333.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Manjurano A, Sepúlveda N, Nadjm B, Mtove G, Wangai H, Maxwell C, et al. USP38, FREM3, SDC1, DDC, and LOC727982 gene polymorphisms and differential susceptibility to severe malaria in Tanzania. J Infect Dis. 2015;212(7):1129–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Toure O, Konate S, Sissoko S, Niangaly A, Barry A, Sall AH, et al. Candidate polymorphisms and severe malaria in a Malian population. PLoS ONE. 2012;7(9): e43987.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Jallow M, Teo YY, Small KS, Rockett KA, Deloukas P, Clark TG, et al. Genome-wide and fine-resolution association analysis of malaria in West Africa. Nat Genet. 2009;41(6):657–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Band G, Le QS, Jostins L, Pirinen M, Kivinen K, Jallow M, et al. Imputation-based meta-analysis of severe malaria in three African populations. PLoS Genet. 2013;9(5): e1003509.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Poostchi M, Silamut K, Maude RJ, Jaeger S, Thoma G. Image analysis and machine learning for detecting malaria. Transl Res. 2018;194:36–55.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Das DK, Ghosh M, Pal M, Maiti AK, Chakraborty C. Machine learning approach for automated screening of malaria parasite using light microscopic images. Micron. 2013;45:97–106.

    Article  PubMed  Google Scholar 

  25. Sharma V, Kumar A, Panat L, Karajkhede G, Lele A. Malaria outbreak prediction model using machine learning. Int J Adv Res Comput Eng Technol. 2015;4(12):4415–9.

    Google Scholar 

  26. Thakur S, Dharavath R. Artificial neural network based prediction of malaria abundances using big data: a knowledge capturing approach. Clin Epidemiol Glob Health. 2019;7(1):121–6.

    Article  Google Scholar 

  27. Sarkar RR, Chatterjee C. Application of different time series models on epidemiological data-comparison and predictions for malaria prevalence. SM J Biom Biostat. 2017;2(4):1022.

    Google Scholar 

  28. Network MGE. New insights into malaria susceptibility from the genomes of 17,000 individuals from Africa, Asia, and Oceania. Nat Commun. 2019;10:1–19.

    CAS  Google Scholar 

  29. Clark TG, Fry AE, Auburn S, Campino S, Diakite M, Green A, et al. Allelic heterogeneity of G6PD deficiency in West Africa and severe malaria susceptibility. Eur J Hum Genet. 2009;17(8):1080–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wilson JN, Rockett K, Jallow M, Pinder M, Sisay-Joof F, Newport M, et al. Analysis of IL10 haplotypic associations with severe malaria. Genes Immun. 2005;6(6):462–6.

    Article  CAS  PubMed  Google Scholar 

  31. Apinjoh TO, Anchang-Kimbi JK, Njua-Yafi C, Ngwai AN, Mugri RN, Clark TG, et al. Association of candidate gene polymorphisms and TGF-beta/IL-10 levels with malaria in three regions of cameroon: a case–control study. Malar J. 2014;13(1):236–47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Manjurano A, Clark TG, Nadjm B, Mtove G, Wangai H, Sepulveda N, et al. Candidate human genetic polymorphisms and severe malaria in a Tanzanian population. PLoS ONE. 2012;7(10): e47463.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Eid NA, Hussein AA, Elzein AM, Mohamed HS, Rockett KA, Kwiatkowski DP, Ibrahim ME. Candidate malaria susceptibility/protective SNPs in hospital and population-based studies: the effect of sub-structuring. Malar J. 2010;9(1):119.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Ayodo G, Price AL, Keinan A, Ajwang A, Otieno MF, Orago AS, et al. Combining evidence of natural selection with association analysis increases power to detect malaria-resistance variants. Am J Hum Genet. 2007;81(2):234–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Atkinson A, Barbier M, Afridi S, Fumoux F, Rihet P. Evidence for epistasis between hemoglobin C and immune genes in human P. falciparum malaria: a family study in Burkina Faso. Genes Immun. 2011;12(6):481–9.

    Article  CAS  PubMed  Google Scholar 

  36. Dewasurendra RL, Suriyaphol P, Fernando SD, Carter R, Rockett K, Corran P, et al. Genetic polymorphisms associated with anti-malarial antibody levels in a low and unstable malaria transmission area in southern Sri Lanka. Malar J. 2012;11(1):281.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Mombo LE, Ntoumi F, Bisseye C, Ossari S, Lu CY, Nagel RL, et al. Human genetic polymorphisms and asymptomatic Plasmodium falciparum malaria in Gabonese schoolchildren. Am J Trop Med Hyg. 2003;68(2):186–90.

    Article  CAS  PubMed  Google Scholar 

  38. Walley AJ, Aucan C, Kwiatkowski D, Hill AV. Interleukin-1 gene cluster polymorphisms and susceptibility to clinical malaria in a Gambian case–control study. Eur J Hum Genet. 2004;12(2):132–8.

    Article  CAS  PubMed  Google Scholar 

  39. Santos SD, Clark TG, Campino S, Suarez-Mutis MC, Rockett KA, Kwiatkowski DP, et al. Investigation of host candidate malaria-associated risk/protective SNPs in a Brazilian Amazonian population. PLoS ONE. 2012;7(5): e36692.

    Article  CAS  Google Scholar 

  40. Gelabert P, Olalde I, De-Dios T, Civit S, Lalueza-Fox C. Malaria was a weak selective force in ancient Europeans. Sci Rep. 2017;7(1):1377.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Caetano CP, Kraaijenbrink T, Tuladhar NM, Driem GLV, Knijff P, Tyler-Smith C, et al. Nepalese populations show no association between the distribution of malaria and protective alleles. J Mol Genet Med. 2006;2(1):101–6.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Ravenhall M, Campino S, Sepúlveda N, Manjurano A, Nadjm B, Mtove G, et al. Novel genetic polymorphisms associated with severe malaria and under selective pressure in North-eastern Tanzania. PLoS Genet. 2018;14(1): e1007172.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Kariuki SM, Rockett K, Clark TG, Reyburn H, Agbenyega T, Taylor TE, et al. The genetic risk of acute seizures in African children with falciparum malaria. Epilepsia. 2013;54(6):990–1001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Flori L, Delahaye NF, Iraqi FA, Hernandez-Valladares M, Fumoux F, Rihet P. TNF as a malaria candidate gene: polymorphism-screening and family-based association analysis of mild malaria attack and parasitemia in Burkina Faso. Genes Immun. 2005;6(6):472–80.

    Article  CAS  PubMed  Google Scholar 

  45. Clark TG, Diakite M, Auburn S, Campino S, Fry AE, Green A, et al. Tumor necrosis factor and lymphotoxin-α polymorphisms and severe malaria in African populations. J Infect Dis. 2009;199(4):569–75.

    Article  CAS  PubMed  Google Scholar 

  46. Dunstan SJ, Rockett KA, Quyen NT, Teo YY, Thai CQ, Hang NT, et al. Variation in human genes encoding adhesion and proinflammatory molecules are associated with severe malaria in the Vietnamese. Genes Immun. 2012;13(6):503–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Maiga B, Dolo A, Touré O, Dara V, Tapily A, Campino S, et al. Human candidate polymorphisms in sympatric ethnic groups differing in malaria susceptibility in Mali. PLoS ONE. 2013;8(10): e75675.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Diakite M, Achidi EA, Achonduh O, Craik R, Djimde AA, Evehe MS, et al. Host candidate gene polymorphisms and clearance of drug-resistant Plasmodium falciparum parasites. Malar J. 2011;10(1):250.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Sepúlveda N, Manjurano A, Campino SG, Lemnge M, Lusingu J, Olomi R, et al. Malaria host candidate genes validated by association with current, recent, and historical measures of transmission intensity. J Infect Dis. 2017;216(1):45–54.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009;5(6): e1000529.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Browning BL, Zhou Y, Browning SR. A one-penny imputed genome from next-generation reference panels. Am J Hum Genet. 2018;103(3):338–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Long Q, Zhang Q, Ott J. Detecting disease-associated genotype patterns. BMC Bioinform. 2009;10(Suppl 1):S75.

    Article  CAS  Google Scholar 

  53. Nielsen DA, Ji F, Yuferov V, Ho A, Chen A, Levran O, et al. Genotype patterns that contribute to increased risk for or protection from developing heroin addiction. Mol Psychiatry. 2008;13(4):417–28.

    Article  CAS  PubMed  Google Scholar 

  54. Archer NM, Petersen N, Clark MA, Buckee CO, Childs LM, Duraisingh MT. Resistance to Plasmodium falciparum in sickle cell trait erythrocytes is driven by oxygen-dependent growth inhibition. Proc Natl Acad Sci. 2018;115(28):7350–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Williams TN, Mwangi TW, Roberts DJ, Alexander ND, Weatherall DJ, Wambua S, et al. An immune basis for malaria protection by the sickle cell trait. PLoS Med. 2005;2(5): e128.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Luzzatto L. Sickle cell anaemia and malaria. Mediterr J Hematol Infect Dis. 2012;4(1): e2012065.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Hüls A, Krämer U, Carlsten C, Schikowski T, Ickstadt K, Schwender H. Comparison of weighting approaches for genetic risk scores in gene-environment interaction studies. BMC Genet. 2017;18(1):115.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Paré G, Mao S, Deng WQ. A machine-learning heuristic to improve gene score prediction of polygenic traits. Sci Rep. 2017;7(1):12665.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Placek K, Benatar M, Wuu J, Rampersaud E, Hennessy L, Van Deerlin VM, et al. Machine learning suggests polygenic risk for cognitive dysfunction in amyotrophic lateral sclerosis. EMBO Mol Med. 2021;13(1): e12595.

    Article  CAS  PubMed  Google Scholar 

  60. Zhang SJ, Meng P, Zhang J, Jia P, Lin J, Wang X, et al. Machine learning models for genetic risk assessment of infants with non-syndromic orofacial cleft. Genom Proteoms Bioinform. 2018;16(5):354–64.

    Article  Google Scholar 

  61. Romagnoni A, Jégou S, Steen KV, Wainrib G, Hugot JP. Comparative performances of machine learning methods for classifying Crohn Disease patients using genome-wide genotyping data. Sci Rep. 2019;9(1):10351.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Gaudillo J, Rodriguez JJ, Nazareno A, Baltazar LR, Vilela J, Bulalacao R, et al. Machine learning approach to single nucleotide polymorphism-based asthma prediction. PLoS ONE. 2019;14(12): e0225574.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Guzzetta G, Jurman G, Furlanello C. A machine learning pipeline for quantitative phenotype prediction from genotype data. BMC Bioinform. 2010;11(Suppl 8):S3.

    Article  Google Scholar 

  64. Behravan H, Hartikainen JM, Tengström M, Pylkäs K, Winqvist R, Kosma VM, et al. Machine learning identifies interacting genetic variants contributing to breast cancer risk: a case study in Finnish cases and controls. Sci Rep. 2018;8(1):13149.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Huang LC, Hsu SY, Lin E. A comparison of classification methods for predicting Chronic Fatigue Syndrome based on genetic data. J Transl Med. 2009;7(1):81.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Ma X, Wu Y, Zhang L, Yuan W, Yan L, Fan S, et al. Comparison and development of machine learning tools for the prediction of chronic obstructive pulmonary disease in the Chinese population. J Transl Med. 2020;18(1):146.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Kulynych B, Troncoso C. Feature importance scores and lossless feature pruning using Banzhaf power indices; 2017. arXiv preprint arXiv:1711.04992.

  68. Saarela M, Jauhiainen S. Comparison of feature importance measures as explanations for classification models. SN Appl Sci. 2021;3(2):1–2.

    Article  Google Scholar 

  69. Mathew TE. A logistic regression with recursive feature elimination model for breast cancer diagnosis. Int J Emerg Technol. 2019;10(3):55–63.

    Google Scholar 

  70. Joseph VR. Optimal ratio for data splitting. Stat Anal Data Min: The ASA Data Sci Journal. 2022;15(4):531–8.

  71. Chai T, Draxler RR. Root mean square error (RMSE) or mean absolute error (MAE)? Arguments against avoiding RMSE in the literature. Geosci Model Dev. 2014;7(3):1247–50.

    Article  Google Scholar 

  72. Visscher PM, Wray NR, Zhang Q, Sklar P, McCarthy MI, Brown MA, et al. 10 years of GWAS discovery: biology, function, and translation. Am J Hum Genet. 2017;101(1):5–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Okser S, Pahikkala T, Aittokallio T. Genetic variants and their interactions in disease risk prediction–machine learning and network perspectives. BioData Min. 2013;6(1):5.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Finkelstein J, Jeong IC. Machine learning approaches to personalize early prediction of asthma exacerbations. Ann N Y Acad Sci. 2017;1387(1):153–65.

    Article  PubMed  Google Scholar 

  75. Waljee AK, Wallace BI, Cohen-Mekelburg S, Liu Y, Liu B, Sauder K, et al. Development and validation of machine learning models in prediction of remission in patients with moderate to severe Crohn disease. JAMA Netw Open. 2019;2(5): e193721.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Szymczak S, Biernacka JM, Cordell HJ, González-Recio O, König IR, Zhang H, et al. Machine learning in genome-wide association studies. Genet Epidemiol. 2009;33(Suppl 1):S51–7.

    Article  PubMed  Google Scholar 

  77. Altshuler D, Daly MJ, Lander ES. Genetic mapping in human disease. Science. 2008;322(5903):881–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Byrne BM. Structural equation modeling with AMOS: basic concepts, applications, and programming. 2nd ed. New York: Taylor & Francis Group; 2010.

    Google Scholar 

  79. Hair JF, Black WC, Babin BJ, Anderson RE. Multivariate data analysis: a global perspective. 7th ed. New Jersey: Pearson Education Inc; 2010.

    Google Scholar 

  80. Rowe JA, Handel IG, Thera MA, Deans AM, Lyke KE, Koné A, et al. Blood group O protects against severe Plasmodium falciparum malaria through the mechanism of reduced rosetting. Proc Natl Acad Sci. 2007;104(44):17471–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Rowe JA, Opi DH, Williams TN. Blood groups and malaria: fresh insights into pathogenesis and identification of targets for intervention. Curr Opin Hematol. 2009;16(6):480–7.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This study makes use of data generated by MalariaGEN. A full list of the investigators who contributed to the generation of the data is available from www.MalariaGEN.net. Funding for this project was provided by Wellcome Trust (WT077383/Z/05/Z) and the Bill & Melinda Gates Foundation through the Foundation of the National Institutes of Health (566) as part of the Grand Challenges in Global Health Initiative.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

TKY performed analysis, interpreted the results, designed the algorithm and drafted the manuscript. JD supervised the project, interpreted the results, designed the algorithm and revised the manuscript. KSW revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jasbir Dhaliwal.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Table S1. General information on the 104 SNPs used in this study.

Additional file 2.

Table S2. A list of features ranked in ascending order obtained from our approach.

Additional file 3.

Table S3. Results of different feature sets based on SNPs genotype data with wGRS as the target variable. Performance metrics are based on MAE and RMSE scores.

Additional file 4.

Table S4. Results of different feature sets based on SNPs genotype data with wGRS+GF as the target variable. Performance metrics are based on MAE and RMSE scores.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tai, K.Y., Dhaliwal, J. & Wong, K. Risk score prediction model based on single nucleotide polymorphism for predicting malaria: a machine learning approach. BMC Bioinformatics 23, 325 (2022). https://doi.org/10.1186/s12859-022-04870-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12859-022-04870-0

Keywords