Screening and identification of potential novel biomarker for diagnosis of complicated Plasmodium vivax malaria
Journal of Translational Medicine volume 16, Article number: 272 (2018)
In the recent years Plasmodium vivax has been reported to cause severe infections associated with mortality. Clinical evaluation has limited accuracy for the early identification of the patients progressing towards the fatal condition. Researchers have tried to identify the serum and the plasma-based indicators of the severe malaria. Discovery of MicroRNA (miRNA) has opened up an era of identification of early biomarkers for various infectious and non-infectious diseases. MicroRNAs (miRNA) are the small non-coding RNA molecules of length 19–24 nts and are responsible for the regulation of the majority of human gene expressions at post transcriptional level.
We identified the differentially expressed miRNAs by microarray and validated the selected miRNAs by qRT-PCR. We assessed the diagnostic potential of these up-regulated miRNAs for complicated P. vivax malaria. Futher, the bioinformtic analysis was performed to construct protein–protein and mRNA–miRNA networks to identify highly regulated miRNA.
In the present study, utility of miRNA as potential biomarker of complicated P. vivax malaria was explored. A total of 276 miRNAs were found to be differentially expressed by miRNA microarray and out of which 5 miRNAs (hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378-5p, hsa-miR-194-5p and hsa-miR-3667-5p) were found to be significantly up-regulated in complicated P. vivax malaria patients using qRT-PCR. The diagnostic potential of these 5 miRNAs were found to be significant with sensitivity and specificity of 60–71% and 69–81% respectively and area under curve (AUC) of 0.7 (p < 0.05). Moreover, in silico analysis of the common targets of up-regulated miRNAs revealed UBA52 and hsa-miR-7977 as majorly regulated hubs in the PPI and mRNA–miRNA networks, suggesting their putative role in complicated P. vivax malaria.
miR-7977 might act as a potential biomarker for differentiating complicated P. vivax malaria from uncomplicated type. The elevated levels of miR-7977 may have a role to play in the disease pathology through UBA52 or TGF-beta signalling pathway.
Malaria is a one of the most important life-threatening parasitic infection affecting human beings. Out of the five species of Plasmodium genus responsible for causing infection in humans, Plasmodium falciparum and Plasmodium vivax are the two most important responsible for causing the majority of the malaria infections . Globally 4% of the malaria infections are caused by P. vivax alone and India is one of five countries responsible for causing 85% of P.vivax infections in 2016 . Over the past few years, the belief of P. vivax responsible for causing benign infections in humans is being challenged by the reports of morbidity and mortality caused by the severe P. vivax infections [3,4,5,6]. There have been increasing reports of the presentation of severe P. vivax malaria which includes the symptoms ranging from altered sensorium, seizures, cerebral malaria, jaundice, acute respiratory distress syndrome (ARDS) shock, acute kidney injury (AKI) and severe anaemia [7,8,9,10,11].
In the recent years with the implementation of molecular diagnosis, it has become possible to demonstrate P. vivax as a solo cause for the underlying multi organ dysfunctions and the life-threatening conditions similar to those caused by P. falciparum . However, microscopy is routinely being used for the diagnosis of malaria and remains the gold standard. Also rapid antigen detection test (RDT) and molecular tests (PCR and LAMP) have also been used for the diagnosis of malaria . But we do not have any technique that can be used in order to differentiate the complicated P. vivax malaria from the uncomplicated P. vivax malaria. Therefore, early diagnosis of these complicated P. vivax and accurate treatment is an important means to prevent the progression of uncomplicated P. vivax malaria to a severe form and ultimately death.
An ideal biomarker for a disease should be less-invasive, stable and present in body fluids in enough amounts for easy diagnosis of the disease . miRNA could represent itself as an ideal marker for being stable, single stranded, small, non-coding, RNA molecule which is evolutionary conserved and has an important role in regulating the translation of mRNA. miRNAs are 19–24 nts long, regulating majority of the human gene expressions at the post transcriptional level by targeted RNA degradation and translational arrest . The researchers have established a significant correlation between the miRNA and the cause of the underlying disease. Recently, the potential role of miRNA as biomarker has been widely studied in many infectious and non-infectious diseases [16,17,18,19]. In various viral, bacterial and parasitic infections the dysregulation of these miRNA has been elucidated and well correlated with the underlying disease pathologies . Among parasitic diseases, in schistosome infection, the potential role of circulating miR-223 as new biomarker and the assessment of the response to chemotherapy have been reported . Tiwari et al.  have reported that consistent dysregulation of 85 miRNA in Leishmania infection with the specific involvement of 10 miRNAs in the regulation of macrophage effector functions in infected macrophages. Among the Apicomplexan parasites, in T. gondii-infected mice, plasma miRNAs were found to be differentially expressed suggesting the potential role as an early biomarker of T. gondii infection . The role of miRNAs in host response to Cryptosporidium has also been well elucidated in number of functional studies . However, a few studies have examined the role of miRNA in malaria specifically. In malaria the pathophysiological role of miRNA has not been well established . Although, there are few studies explaining the probable role of these miRNA as a putative biomarker for P. vivax malaria but none explains the potential for their use as a biomarker for severe P. vivax malaria. Here, in the present study, the differential miRNA expression in whole blood samples of complicated P. vivax group of patients were studied using miRNA microarray and the potential of miRNA was then explored to be used as an early predictor of complications in P. vivax malaria.
Sample classification for miRNA microarray and qRT-PCR
In total, 48 whole blood samples have been collected for miRNA studies. The samples have been classified into four groups as shown in Table 1. Group I included healthy controls which are negative for P. vivax by multiplex nested PCR. Group II and III patients were positive for P. vivax by microscopy, confirmed by molecular techniques (multiplex nested PCR, Real time PCR and LAMP) [26, 27]. The samples in these groups were classified as complicated and uncomplicated P. vivax on the basis of the WHO based criteria for severe malaria patients . Group IV consisted of P. falciparum positive patients, which were confirmed by microscopy and multiplex nested PCR. Group I and Group IV subjects were enrolled for miRNA microarray only. The demographic and clinical details from the enrolled patients were obtained at the time of sample collection.
Approximately, 4 ml of blood sample was taken from all subjects by a trained practitioner. For microscopy and molecular detection methods ~ 1.0 ml of the sample was collected in EDTA vacutainer. And 2.5 ml of the sample was collected in the PAXgene blood RNA tubes (PreAnalytics-BD, New Jersey, USA) and kept at − 20 °C for 24 h and later stored at − 80 °C for further use.
On the day of RNA isolation, the samples were allowed to thaw completely at room temperature (18–25 °C) for at least 2–3 h. From the whole blood samples, the total RNA was extracted using the PAXgene blood RNA kit (Qiagen, Germany) as per the manufacturer’s instruction. The concentration and purity of the total RNA was estimated using NanoDrop-1000 Spectrophotometer (Thermo Fisher Scientific Inc., Wilmington, DE) by measuring 260/280 ratio. The RNA with ratio of 1.9–2.2 were further processed for real time PCRs. The integrity of the respective samples was also checked by doing the native gel electrophoresis.
miRNA microarray using Affymetrix miRNA 4.0 array (format 100 midi)
For miRNA microarray, the total RNA was isolated as mentioned above following manufacturer’s instructions. The RNA integrity number was assessed using Agilent’s Bioanalyzer 2100 (Agilent Technologies Inc., USA) for each sample. As per the affymetrix miRNA protocol of miRNA 4.0 array; the steps of fragmentation and hybridization was performed followed by the scanning using Affymetrix scanner 3000 7G. The miRNA 4.0 array provided a comprehensive coverage for miRNAs for human (Affymetrix, Santa Clara, CA, USA). The microarray raw data was extracted from the files (raw intensity file) generated during scanning of slides. These raw data sets were separately analyzed using Expression Console and Trancriptome console software followed by a differential miRNA expression, fold change and cluster analysis.
Endogenous control selection for normalization of qRT-PCR data
As per the previous published studies of the endogenous controls (EC) used for whole blood samples, in total five genes (3 snoRNA and 2 miRNA) namely U6, U44, U48, miR425 and miR16 were selected to test their potential as an endogenous control in our study [16, 29, 30]. In total of 18 samples (10 complicated and 8 uncomplicated) were used to identify the suitable endogenous controls for qRT-PCR.
Candidate EC stability analysis
geNorm, NormFinder, bestkeeper alogrithms and ΔCt ± SD methods were used to assess the stability of the candidate ECs. geNorm was used to rank the endogenous controls as per their stability values (M) representing the variation in expression of each EC in comparison to each other . Normfinder is a Microsoft excel add-in tool that take in account both inter and intra group variation (complicated and uncomplicated). This gives the stability value; the lower the value the more stable will be the EC gene. The converted quantities (2−ΔCt (ΔCt = the corresponding Cq value-minimum Cq)  were used for analysis in both geNorm and Normfinder analysis.
Candidate miRNA selection criteria for validation by quantitative Real time PCR (qRT-PCR)
A profile of differentially expressed small non-coding RNA was obtained from miRNA microarray with 100% miRBase v20 coverage. For validation by qRT-PCR, the selection of mature miRNAs was made by making venn between the complicated vs uncomplicated P. vivax, complicated P. vivax vs healthy and complicated P. vivax vs P. falciparum. The miRNAs which were then found to be common in any of the two or three groups with a fold change of FC ≥ 3.5 and p < 0.05 and miRNAs in complicated vs uncomplicated P. vivax group with a fold change of FC ≥ 5.0 and p < 0.05 were then selected. Out of the selected miRNA, those miRNAs were removed from further analysis which were found to be common in P. falciparum vs healthy and uncomplicated P. vivax vs healthy groups as shown in Fig. 1.
A total of 12 mature miRNAs were selected for validation by qRT-PCR. The primer sequences used for these 12 mature miRNAs are shown in Table 2. Specific miRNA forward and universal reverse primers were used for qRT-PCR.
miRNA validation: cDNA synthesis and quantitative Real time PCR (qRT-PCR)
For cDNA synthesis, 1 µg of total RNA sample for each isolate was reverse transcribed in a final volume of 20 µl using the miRNA 1st-Strand cDNA Synthesis Kit (Agilent Technologies, Santa Clara, California) as per manufacturer’s instruction. Briefly, polyadenylation step was first performed at 37 °C for 30 min and termination of the reaction was carried out at 95 °C for 5 min. After obtaining the polyadenylated product, the 1st strand cDNA synthesis was performed at 55 °C for 5 min, 25 °C for 15 min, 42 °C for 30 min and finally incubated at 95 °C for 5 min to terminate reverse transcription reaction. The aliquots of the obtained cDNAs was diluted to 15-fold and stored at − 20 °C for further use. The selected 12 miRNAs were quantified using Maxima SYBR Green/ROX qPCR Master Mix (2X) (Thermo Scientific, USA) in a 20 µl reaction, consisting of 1× mastermix, 200 nM of primers (miRNA specific forward and universal reverse) and 2 µl of cDNA. The housekeepings were run simultaneously with each sample reaction. The qRT-PCR was performed in 96 well plate using the Applied Biosystems 7500 step one real-time PCR system (Applied Biosystems, Carlsbad, CA, USA) using a three-step protocol of 40 cycles of 95 °C for 10 s, 58 °C for 15 s and data was acquired at 72 °C for 15 s. Each sample was run in triplicates. ΔΔCt method was used to carry out the relative quantitation for obtaining the results for each miRNA. The ΔCts were calculated as difference of Ct values between the target miRNAs and the housekeeping gene.
Diagnostic potential of putative biomarker of complicated P. vivax malaria
The ROC curve was generated for the miRNAs, which were significantly up-regulated by qRT-PCR in the complicated Plasmodium vivax group. The area under curve (AUC) value and 95% confidence intervals were calculated to estimate the sensitivity and specificity of significantly up-regulated miRNAs. All the statistical analysis was carried out using Graph Pad prism. version 5.01.
Target gene selection, functional and pathway enrichment analysis
The qRT-PCR up-regulated miRNAs were then further selected for target gene analysis using TargetScan . The list of potential targeted genes was obtained individually for each miRNA. The DAVID server with default setting was used to analyze the involvement of target genes (of all up-regulated miRNA individually and combined) in biological processes and for functional analysis . The KEGG and reactome pathways were used to exemplify the pathways of the predicted miRNA target genes using functional annotation tool of DAVID using human genes as background set. The cut off of p < 0.01 was used for the identification of potential enriched KEGG pathways.
Network construction and analysis
The target genes common between any of the two miRNAs were then further used for the construction of protein- protein interaction network using STRING database . To study the protein–protein interaction the data obtained from STRING was visualized in Cytoscape Software (Cytoscape Software, Version 2.8.2, Seattle, USA) using the group attributed layout based on the degree of all nodes . The mRNA–miRNA network was constructed using nodes with degree of more than 20 in cytoscape and visualized using degree sorted circular layout.
To determine the differential expression of mature miRNAs between the complicated and uncomplicated P. vivax groups, Mann–Whitney U test was applied. The p < 0.05 was considered as statistical significance. All the statistical analysis was carried out using Graph Pad prism. version 5.01.
The male to female ratio and median age [interquartile range (IQR)] of the patients in complicated, uncomplicated group was found to be 2.7:1, 20 (11–27) years and 0.8:1, 12 (6–29) years respectively. In the healthy group (n = 4), two males and two females with median age of 16 (4–28) years were enrolled in the study and they were not found to have any clinical symptoms of malaria or any other major illness at the time of sample collection. In P. falciparum group, the mean age of the patients was 19 years and they were possessing the acute febrile illness with the presence of severe thrombocytopenia, epistaxis, severe anaemia, jaundice with renal impairment. The 20 uncomplicated P. vivax patients were found to have general symptoms of malaria such as chills and rigors associated fever with general body weakness. The 22 complicated P. vivax patients recruited in the present study had one or the other complications as defined by the WHO criteria for severe malaria. The major complications presented in the 22 complicated group of patients included jaundice 27% (n = 6), seizures with altered sensorium 27% (n = 6), renal impairment, abnormal bleeding, severe anemia 9% (n = 2 each), and ARDS in one patient.
Differentially expressed miRNA in complicated Plasmodium vivax
To interrogate all miRNAs in miRBase Release 20, Affymetrix miRNA 4.0 array is designed. On comparing the profile of miRNA obtained in microarray in the groups, complicated vs uncomplicated P. vivax, complicated P. vivax vs Healthy and complicated P. vivax vs P. falciparum a total of 276 miRNA were found to be differentially expressed (data shown in Additional file 1: Table S1). Out of 276, 155 matute miRNA (32 up-regulated and 123 downregulated), 35 mature miRNAs (12 up-regulated and 23 downregulated), 30 mature miRNAs (22 up-regulated and 8 downregulated) were found to be differentially expressed in complicated vs uncomplicated P. vivax, complicated P. vivax vs Healthy and complicated P. vivax vs P. falciparum. The expression patterns of these miRNAs are displayed in the clustered heat map shown in Fig. 2.
Candidate miRNA selection for quantitative real time PCR validation
In total 12 mature miRNAs (hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378a-5p, hsa-miR-194-5p, hsa-miR-3667-5p, hsa-miR-3651, hsa-miR-192-5p, hsa-miR-19b-3p, hsa-miR-181b-5p, hsa-miR-29a-3p, hsa-miR-200c-3p, hsa-miR-155-5p) were selected for qRT-PCR validation. Six with FC ≥ 3.5 and p < 0.05 were selected on the basis of their presence in any of the two groups (Fig. 3) and 6 miRNAs with FC ≥ 5.0 and p < 0.05 in complicated vs uncomplicated P. vivax group was chosen for validation in a larger number of samples by qRT-PCR.
Selection of candidate housekeeping for quantitative real time PCR (qRT-PCR)
The geNorm, Normfinder, BestKeeper and Ct values ± SD were used to assess the stability of five selected genes (U48, U6, U44, miR-425 and miR-16) (Table 3).
In term of stability (M value), U48 and U6 came out to be most stable gene with M value of 1.3 and 1.4 respectively. In geNorm following the procedure of step-wise exclusion of the candidate gene, miR-16 turned out to worst reference gene with the highest M value of 2.2 (Fig. 4a).
As per the Normfinder analysis, the U48 was the most stable gene with the stability index (ρ) of 0.050 and was found to have the least values of intra-group variability (0.023). The stability of the two genes U6 and U44 was found to be the worst with high ρ values of 0.80 and 0.095 (Fig. 4b).
In terms of variability, U6 displayed the lowest standard deviation value of 0.83 in Bestkeeper analysis. Correlation analysis was performed and BestKeeper index (BKI) was calculated for the candidate genes. On the basis of BKI, U48 and miR-425 showed the highest correlation of r > 0.87 and r > 0.91 and p = 0.001 and U6 had a lower correlation value (r < 0.45 and p = 0.02) (Fig. 4c).
Comparative Ct values ± SD method
Similar to geNorm and BestKeeper, this method also showed miR-16 to be the least stable candidate gene with the highest SD of Ct values of 3.15 respectively. U6 was found to be the most stable gene with the average SD values of 1.06 respectively (Fig. 4d).
On the basis of the calculated geometric mean of ranks, U6 and U48 was found to be at the top rank and thus used for the further quantitative real time PCR experiments (Fig. 4e). miR-16 and miR-425 were found to be on the lowest rank as per combined analysis.
Out of a large number of differentially expressed miRNAs, a total of 12 mature miRNAs were selected (criteria mentioned above) for further validation by using q RT-PCR SYBR green method in both complicated and uncomplicated group of patients. The results of candidate miRNAs were normalized to the geometric mean of two selected housekeeping genes U6 and U48. Out of 12 miRNAs, the expression of 5 miRNAs hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378a-5p, hsa-miR-194-5p, hsa-miR-3667-5p were found to be statistically significantly (p < 0.05) up-regulated in the complicated P. vivax group as compared to the uncomplicated P. vivax group of patients (Fig. 5). No significant difference in the expression was obtained for the rest of the 7 miRNAs.
Evaluation of the diagnostic potential of 5 up-regulated miRNA for complicated Plasmodium vivax infections
To investigate the potential of five significantly up-regulated miRNAs (hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378-5p, hsa-miR-194-5p, hsa-miR-3667-5p) to be used as putative biomarker for P. vivax complicated infections the Receiver operating characteristic (ROC) curve analysis was performed for all five miRNAs individually. The hsa-miR-7977 represented the area under curve (AUC) of 0.7344; 95% CI 0.5434 to 0.9253; p = 0.02378, hsa-miR-28-3p with AUC = 0.7366; 95% CI 0.5579 to 0.9153; p = 0.02762, hsa-miR-378-5p with AUC = 0.7347; 95% CI 0.5606 to 0.9144; p = 0.02429)], hsa-miR-194-5p with AUC = 0.7277; 95% CI 0.5473 to 0.9081; p = 0.03405 and hsa-miR-3667-5p with AUC = 0.7347; 95% CI 0.5491 to 0.9203; p = 0.03460 respectively (Fig. 6).
Gene Ontology analysis of predicted targets of five up-regulated miRNAs
In order to study the role of these five up-regulated miRNAs in complicated P. vivax malaria the putative gene targets of these miRNAs were assessed using TargetScan with a cut off of 0.1 of total context score. A total of 3877, 107, 2307, 396 and 1960 target genes were predicted for hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378-5p, hsa-miR-194-5p and hsa-miR-3667-5p respectively. Then the first 3000 targets in case of hsa-miR-7977 and all gene targets of rest of the miRNAs were then subjected to gene ontology (GO) analysis in DAVID v6.7 tool. At cutoff standard of p < 0.01, a total of 30, 7, 22 and 23 GO terms of biological process, molecular function and cellular component for hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378-5p and hsa-miR-3667-5p was found for the complicated Plasmodium vivax group. No GO term was found for hsa-miR-194-5p using DAVID. A total of 30 GO terms were obtained for combined gene targets, which were common in any of the two miRNAs (Fig. 7a, c, e, g, i).
KEGG and reactome pathway analysis of the predicted targets of the up-regulated microRNAs
To gain an overall knowledge of these miRNA in complicated P. vivax infections, KEGG and Reactome pathway analysis was performed in DAVID v 6.7. The pathways were ranked significant as per their p value (p < 0.01). A total of 18, 2, 3, 4 and 7 pathways were significantly enriched by targeted genes of hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378-5p and hsa-miR-3667-5p and for combined target genes of all five miRNAs. The top pathways in the combined analysis shown here are to be the ones involved in generic transcription pathway, TGF-beta signalling pathway and mRNA-splicing pathways Fig. 7b, d, f, h, j.
Protein-protein interaction (PPI) analysis
The summarized layout for the detailed analysis of PPI networks are shown in Additional file 2: Figure S1. The common gene targets between any of the two miRNAs were selected for the network construction using STRING version: 10.5 (Additional file 3: Figure S2). A total of 1624 nodes with 6445 edges were found in the network with a p value of 0.000398. The data obtained from STRING was uploaded into the cytoscape (Fig. 8). The nodes with combined score of 0.9 were then used to construct the network. The highest degree of 71 was found for UBA52. To construct mRNA–miRNA network the nodes with degree of more than 21 were selected. Out of 344 nodes, 13 nodes (UBA52, POLR2D, CSNK1D, TP53, CDK1, GTF2F1, PAFAH1B1, RBM8A, RPS24, RPS14, HIST1H2BD, MAPK1, UBE2V2) were found to have degree range from 21 to 71. The mRNA–miRNA network was constructed using cytoscape and hsa-miR-7977 was found to have the highest degree of 11 as compared to other miRNAs (Fig. 9). Hence, miR-7977 might be a major regulator of the complications and could be used as a potential biomarker for the prognosis of complications in P. vivax malaria.
In the recent years, the importance of severe P. vivax is being outweighed by P. falciparum which is known to be associated with fatal infections . A number of studies from different geographical regions viz. Thailand, Brazil, Indonesia, Papua New Guinea and India have reported an increase in the number of severe malaria caused by P. vivax malaria [4, 37]. In P. vivax the mechanism initiating the transformation of uncomplicated malaria to severe P. vivax malaria still remains poorly understood. Increased deformability of red blood cells and almost absence of parasite sequestration in contrast to P. falciparum is capable of causing fatal manifestations . An early detection with initiation of correct treatment and close monitoring of the patient would be lifesaving. In the previous studies, the association of various plasma detectable markers (PCT, CRP, ICAM, ANG2, TNFα and IP-10) with mortality in patients with severe or cerebral malaria have been reported . Recently the research has focused mainly on the potential role of miRNA in the underlying normal and pathophysiological conditions of various diseases .
Till date only a handful of available studies have highlighted the probable link of miRNA in the pathogenesis and as a diagnostic marker of P. falciparum and P. vivax malaria. In a case study of acute P. vivax infection, Baro et al.  have reported a distinct series of differentially expressed miRNA associated with erythropoiesis. Another study has reported the levels of miR-451 and miR-16 to be significantly down regulated in the P. vivax infections as compared to normal subjects suggesting to be the probable biomarker of P. vivax malaria infections . LaMonte et al.  have reported the impaired growth of parasites due to mRNA translation inhibition, which was observed due to increased accumulation of miR-451 and miR-223 in both the HbAS and HbSS erythrocytes. Study on P. chabaudi infected C57BL/6 mice has suggested the association of reprogramming of 19 distinct miRNA in female mouse liver with the development of protective immunity against blood stages of P. chabaudi infection . Assaad et al.  have showed a significant change in the expression of let-7i, miR-150, miR-27a in the brain tissue of PbA infected WT mice, suggestive of their critical involvement in the triggering of the neurological syndrome. A recent study for the first time depicted the innate ability of normal RBCs to prevent infection by P. falciparum parasite by the transfer of the hAgo2-miRNA complexes via microparticles (MPs) to infected RBCs and thus down regulating of the expression of PfEMP1 . Wang et al. have also described the potential role of miR-155 in host response to malaria via BBB dysfunction through microvascular leak and regulation of endothelial activation in cerebral malaria (CM) . Another study conducted on controlled human blood-stage infection (CHMI) has shown the dichotomous miRNA expression distinguishing high-miR responders from low-miR responders with significantly enhanced antimalarial antibody responses in the high-miR responders group . At present there is no approved diagnostic modality available to detect the severity of the P. vivax malaria . In the present study we identified the potential whole blood miRNA marker for the probable diagnosis of complicated P. vivax malaria using miRNA microarray.
The present study is first of its type where miRNA microarray was carried on whole blood samples of P. vivax complicated and uncomplicated malaria patients in order to identify putative biomarker of complicated P. vivax malaria. A total of five miRNAs (hsa-miR-7977, hsa-miR-28-3p, hsa-miR-378-5p, hsa-miR-194-5p, hsa-miR-3667-5p) were found to be statistically significantly (p < 0.03) up-regulated in the complicated P. vivax group of patients as compared to the uncomplicated group of patients using qRT-PCR. The diagnostic potential for all these miRNAs were also assessed and was found to be significant with good sensitivity and specificity. Among these five miRNAs, hsa-miR-7977 had the highest degree of interaction in the mRNA–miRNA networks and thus might have a role to play in the pathogenesis of complicated P. vivax malaria. Previous studies have reported the validity of the reference genes being affected by the type of sample and experimental procedures used . Therefore, in this study we have confirmed our selection of reference genes (U6 and U48) using four different softwares and thus used them for further q RT-PCR.
The information on the pathogenesis of severe anaemia in P. vivax remains majorly obscure. Molecular mechanisms underlying the inhibition of the erythroid cell proliferation and differentiation by P. vivax infection is complex and remains unclear . Studies have provided the evidence of the dyserythropoiesis and ineffective erythropoiesis in P. vivax malaria patient’s bone marrow . Ineffective erythropoiesis due to the P. vivax inhibited erythroid development has been reported to play a potential role in severe anaemia caused due to P. vivax . Recently the role of miR-7977 in the hematopoietic dysfunction of mesenchymal stromal cells by poly (rC) binding protein 1supression has been explained by Horiguchi et al. . The significant upregulation of miR-7977 in the complicated group of patients obtained in this study may suggest their putative role in the regulation of erythropoiesis. However, further studies will be required to establish their probable function in complicated P. vivax malaria.
The TGF-β signalling pathway was found to be the most significantly regulated pathway in the present study. Depending on the environment and concentration, the TGF-β is capable of exhibiting both pro- and anti-inflammatory properties . TGF-β is known to supress production of TNF-α and nitric oxide from macrophages and of IFN-γ and TNF-α from NK cells . The simultaneous increase of TNF and IFN-γ in case of P. vivax malaria infections has been linked with severe disease progression . Study by Omer et al. have reported the low circulating levels of TGF-β, and bioactive TGF-β produced by splenocytes in deadly infections with Plasmodium berghei and significant production of TGF- β in P. chabaudi or P. yoelii resolving infections . Study by Perkins et al.  have also demonstrated the association of low levels of TGF- β1 and IL-12 with the malaria disease severity. In our study the TGF- β signalling pathways being the most effected pathway might have some role to play with the severe P. vivax infections.
In-silico analysis was performed for the five up-regulated miRNAs to identify possible role of these miRNAs in complicated P. vivax malaria. The targeted genes of up-regulated miRNA were found to be mainly enriched in the protein binding, nucleic acid binding, proteasome-mediated ubiquitin-dependent protein catabolic process, Poly (A) RNA binding, generic transcription pathway, TGF-beta signalling pathway, regulation of pyruvate dehydrogenase (PDH) and mRNA-splicing pathways. Generic transcription and m-RNA splicing pathways are one of the majorly regulated pathways by almost all the miRNAs in the present study. Our results are in accordance with the previous belief of the possibility, that Plasmodium might utilize the host miRNAs to regulate their gene expression [42, 59, 60]. In the present study, UBA52 was found to have the highest degree of interaction in the PPI network. Among the nuclear and cytoplasmic proteins, UBA52 has a role in the targeted cellular protein degradation, regulation of gene expression and also involved in the stress response. UBA52 is an important ubiquitin supplier of the ubiquitin pool and also is an important regulator of ribosomal-protein complex  clearly depicting the role in ubiquitin-dependent translation mechanism. The higher degree of interaction of this protein in complicated P. vivax group suggests its potential role in the gene regulation mechanism in the severe disease. However, future studies are required to evaluate the functions of the highly connected nodes in the pathogenesis of complicated Plasmodium vivax malaria.
To summarize the present study, we found five miRNA which are up-regulated in the P. vivax complicated group of patients as compared to the uncomplicated group of patients. These miRNAs have a good sensitivity and specificity to be used as a predictor of severity in P. vivax malaria. hsa-miR-7977 was found to be highly up-regulated in complicated P. vivax and may have a role in exacerbating the disease pathology through UBA52 or TGF-beta signalling pathway. Hence, hsa-miR-7977 should be explored as a potential biomarker for differentiating complicated vs uncomplicated P. vivax. However, the findings of the present study need to be further explored in future studies to understand the disease pathogenesis.
area under curve
acute respiratory distress syndrome
acute kidney injury
intercellular adhesion molecule
tumor necrosis factor
receiver operating characteristic
Database for Annotation, Visualization and Integrated Discovery
Kyoto Encyclopedia of Genes and Genomes
transforming growth factor
search tool for the retrieval of interacting genes/proteins
Ubiquitin A-52 Residue Ribosomal Protein Fusion Product 1
RNA polymerase II subunit D
casein kinase 1 Delta
tumor protein P53
cyclin-dependent kinase 1
general transcription factor IIF subunit 1
platelet activating factor acetylhydrolase 1b regulatory subunit 1
RNA binding motif protein 8A
ribosomal protein S24
ribosomal protein S14
histone cluster 1 H2B family member D
mitogen-activated protein kinase 1
ubiquitin conjugating enzyme E2 V2
Wassmer SC, Taylor TE, Rathod PK, Mishra SK, Mohanty S, Arevalo-Herrera M, et al. Investigating the pathogenesis of severe malaria: a multidisciplinary and cross-geographical approach. Am J Trop Med Hyg. 2015;93(3 Suppl):42–56. https://doi.org/10.4269/ajtmh.14-0841.
World Health Organization. World Malaria Report 2017. Geneva: WHO; 2017. http://www.who.int/malaria/publications/world-malaria-report-2017/en/. Accessed 2 Jul 2018.
Im JH, Kwon HY, Baek J, Park SW, Durey A, Lee KH, et al. Severe Plasmodium vivax infection in Korea. Malaria J. 2017;16(1):51. https://doi.org/10.1186/s12936-017-1684-4.
Rahimi BA, Thakkinstian A, White NJ, Sirivichayakul C, Dondorp AM, Chokejindachai W. Severe vivax malaria: a systematic review and meta-analysis of clinical studies since 1900. Malaria J. 2014;13(1):481. https://doi.org/10.1186/1475-2875-13-481.
Kochar A, Kalra P, Sb V, Ukirade V, Chahar A, Kochar DK, et al. Retinopathy of vivax malaria in adults and its relation with severity parameters. Pathog Glob Health. 2016;110(4–5):185–93. https://doi.org/10.1080/20477724.2016.1213948.
White NJ, Pukrittayakamee S, Hien TT, Faiz MA, Mokuolu OA, Dondorp AM. Malaria. Lancet (London, England). 2014;383(9918):723–35. http://linkinghub.elsevier.com/retrieve/pii/S0140673613600240. Accessed 9 Jul 2018.
Kaur H, Sehgal R, Bansal D, Sultan AA, Bhalla A, Singhi SC. Development of visually improved loop mediated isothermal amplification for the diagnosis of Plasmodium vivax malaria in a tertiary hospital in Chandigarh, North India. Am J Trop Med Hyg. 2018;98(5):1374–81. https://doi.org/10.4269/ajtmh.17-0857.
Prakash J, Singh AK, Kumar NS, Saxena RK. Acute renal failure in Plasmodium vivax malaria. J Assoc Physicians India. 2003;51:265–7.
Beg MA, Khan R, Baig SM, Gulzar Z, Hussain R, Smego RA. Cerebral involvement in benign tertian malaria. Am J Trop Med Hyg. 2002;67(3):230–2.
Lomar AV, Vidal JE, Lomar FP, Barbas CV, de Matos GJ, Boulos M. Acute respiratory distress syndrome due to vivax malaria: case report and literature review. Braz J Infect Dis. 2005;9(5):425–30.
Kochar DK, Das A, Kochar SK, Saxena V, Sirohi P, Garg S, et al. Severe Plasmodium vivax malaria: a report on serial cases from Bikaner in northwestern India. Am J Trop Med Hyg. 2009;80(2):194–8.
Kumari M, Ghildiyal R. Clinical profile of Plasmodium vivax malaria in children and study of severity parameters in relation to mortality: a tertiary care centre perspective in Mumbai, India. Malaria Res Treat. 2014;2014:765657. http://www.hindawi.com/journals/mrt/2014/765657/. Accessed 9 Jul 2018.
Chamnanchanunt S, Fucharoen S, Umemura T. Circulating microRNAs in malaria infection: bench to bedside. Malaria J. 2017;16(1):334. https://doi.org/10.1186/s12936-017-1990-x.
Miotto P, Mwangoka G, Valente IC, Norbis L, Sotgiu G, Bosu R, et al. miRNA signatures in sera of patients with active pulmonary tuberculosis. PloS ONE. 2013;8(11):e80149. https://doi.org/10.1371/journal.pone.0080149.
Verma P, Pandey RK, Prajapati P, Prajapati VK. Circulating microRNAs: potential and emerging biomarkers for diagnosis of human infectious diseases. Front Microbiol. 2016;7:1274. https://doi.org/10.3389/fmicb.2016.01274/abstract.
Serafin A, Foco L, Blankenburg H, Picard A, Zanigni S, Zanon A, et al. Identification of a set of endogenous reference genes for miRNA expression studies in Parkinson’s disease blood samples. BMC Res Notes. 2014;7(1):715. https://doi.org/10.1186/1756-0500-7-715.
Pordzik J, Pisarz K, De Rosa S, Jones AD, Eyileten C, Indolfi C, et al. The potential role of platelet-related microRNAs in the development of cardiovascular events in high-risk populations, including diabetic patients: a review. Front Endocrinol. 2018;9:74. https://doi.org/10.3389/fendo.2018.00074/full.
Qi Y, Cui L, Ge Y, Shi Z, Zhao K, Guo X, et al. Altered serum microRNAs as biomarkers for the early diagnosis of pulmonary tuberculosis infection. BMC Infect Dis. 2012;12(1):384. https://doi.org/10.1186/1471-2334-12-384.
Zhao Q, Sun X, Liu C, Li T, Cui J, Qin C. Expression of the microRNA-143/145 cluster is decreased in hepatitis B virus-associated hepatocellular carcinoma and may serve as a biomarker for tumorigenesis in patients with chronic hepatitis B. Oncol Lett. 2018;15(5):6115–22. https://doi.org/10.3892/ol.2018.8117.
Cohen A, Combes V, Grau GE. MicroRNAs and malaria—a dynamic interaction still incompletely understood. J Neuroinfect Dis. 2015;6(1):165.
He X, Sai X, Chen C, Zhang Y, Xu X, Zhang D, et al. Host serum miR-223 is a potential new biomarker for Schistosoma japonicum infection and the response to chemotherapy. Parasit Vectors. 2013;6(1):272. https://doi.org/10.1186/1756-3305-6-272.
Tiwari N, Kumar V, Gedda MR, Singh AK, Singh VK, Gannavaram S, et al. Identification and characterization of miRNAs in response to Leishmania donovani infection: delineation of their roles in macrophage dysfunction. Front Microbiol. 2017;8:314. https://doi.org/10.3389/fmicb.2017.00314/full.
Jia B, Chang Z, Wei X, Lu H, Yin J, Jiang N, et al. Plasma microRNAs are promising novel biomarkers for the early detection of Toxoplasma gondii infection. Parasit Vectors. 2014;7(1):433. https://doi.org/10.1186/1756-3305-7-433.
Judice CC, Bourgard C, Kayano ACAV, Albrecht L, Costa FTM. MicroRNAs in the host-apicomplexan parasites interactions: a review of immunopathological aspects. Front Cell Infect Microbiol. 2016;6:5. https://doi.org/10.3389/fcimb.2016.00005/abstract.
Thirugnanam S, Rout N, Gnanasekar M. Possible role of Toxoplasma gondii in brain cancer through modulation of host microRNAs. Infect Agents Cancer. 2013;8(1):8. https://doi.org/10.1186/1750-9378-8-8.
Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah MS, Rahman HA. A genus- and species-specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. Am J Trop Med Hyg. 1999;60(4):687–92.
Han ETET, Watanabe R, Sattabongkot J, Khuntirat B, Sirichaisinthop J, Iriko H, et al. Detection of four Plasmodium species by genus- and species-specific loop-mediated isothermal amplification for clinical diagnosis. J Clin Microbiol. 2007;45(8):2521–8.
World Health Organization (WHO). Management of severe malaria—a practical handbook, 3th edition. WHO; 2012. http://www.who.int/malaria/publications/atoz/9789241548526/en/. Accessed 9 Jul 2018.
Parsi S, Soltani BM, Hosseini E, Tousi SE, Mowla SJ. Experimental verification of a predicted intronic microRNA in human NGFR gene with a potential pro-apoptotic function. PloS ONE. 2012;7(4):e35561. https://doi.org/10.1371/journal.pone.0035561.
McDermott AM, Kerin MJ, Miller N. Identification and validation of miRNAs as endogenous controls for RQ-PCR in blood specimens for breast cancer studies. PloS ONE. 2013;8(12):e83718. https://doi.org/10.1371/journal.pone.0083718.
Zhao H, Ma T-F, Lin J, Liu L-L, Sun W-J, Guo L-X, et al. Identification of valid reference genes for mRNA and microRNA normalisation in prostate cancer cell lines. Scientific Rep. 2018;8(1):1949. http://www.nature.com/articles/s41598-018-19458-z. Accessed 3 Jul 2018.
Agarwal V, Bell GW, Nam J-W, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. eLife. 2015;4:e05005. https://elifesciences.org/articles/05005. Accessed 3 Jul 2018.
Jiao X, Sherman BT, Huang DW, Stephens R, Baseler MW, Lane HC, et al. DAVID-WS: a stateful web service to facilitate gene/protein list analysis. Bioinformatics. 2012;28(13):1805–6.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8. https://doi.org/10.1093/nar/gkw937.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Price RN, Tjitra E, Guerra CA, Yeung S, White NJ, Anstey NM. Vivax malaria: neglected and not benign. Am J Trop Med Hyg. 2007;77(6 Suppl):79–87.
Naing C, Whittaker MA, Nyunt Wai V, Mak JW. Is Plasmodium vivax malaria a severe malaria? A systematic review and meta-analysis. PLoS Negl Trop Dis. 2014;8(8):e3071. https://doi.org/10.1371/journal.pntd.0003071.
Ray S, Patel SK, Venkatesh A, Bhave A, Kumar V, Singh V, et al. Clinicopathological analysis and multipronged quantitative proteomics reveal oxidative stress and cytoskeletal proteins as possible markers for severe vivax malaria. Scientific Rep. 2016;6(1):24557. http://www.nature.com/articles/srep24557. Accessed 4 Jul 2018.
Rubio M, Bassat Q, Estivill X, Mayor A. Tying malaria and microRNAs: from the biology to future diagnostic perspectives. Malar J. 2016;15(1):167.
Baro B, Deroost K, Raiol T, Brito M, Almeida ACG, de Menezes-Neto A, et al. Plasmodium vivax gametocytes in the bone marrow of an acute malaria patient and changes in the erythroid miRNA profile. PLoS Negl Trop Dis. 2017;11(4):e0005365. https://doi.org/10.1371/journal.pntd.0005365.
Chamnanchanunt S, Kuroki C, Desakorn V, Enomoto M, Thanachartwet V, Sahassananda D, et al. Downregulation of plasma miR-451 and miR-16 in Plasmodium vivax infection. Exp Parasitol. 2015;155:19–25. http://linkinghub.elsevier.com/retrieve/pii/S0014489415001022. Accessed 3 Jul 2018.
LaMonte G, Philip N, Reardon J, Lacsina JR, Majoros W, Chapman L, et al. Translocation of sickle cell erythrocyte microRNAs into Plasmodium falciparum inhibits parasite translation and contributes to malaria resistance. Cell Host Microbe. 2012;12(2):187–99.
Delić D, Dkhil M, Al-Quraishy S, Wunderlich F. Hepatic miRNA expression reprogrammed by Plasmodium chabaudi malaria. Parasitol Res. 2011;108(5):1111–21. https://doi.org/10.1007/s00436-010-2152-z.
El-Assaad F, Hempel C, Combes V, Mitchell AJ, Ball HJ, Kurtzhals JAL, et al. Differential microRNA expression in experimental cerebral and noncerebral malaria. Infect Immun. 2011;79(6):2379–84. https://doi.org/10.1128/iai.01136-10.
Wang Z, Xi J, Hao X, Deng W, Liu J, Wei C, et al. Red blood cells release microparticles containing human argonaute 2 and miRNAs to target genes of Plasmodium falciparum. Emerg Microbes Infect. 2017;6(8):e75.
Barker KR, Lu Z, Kim H, Zheng Y, Chen J, Conroy AL, et al. miR-155 modifies inflammation, endothelial activation and blood-brain barrier dysfunction in cerebral malaria. Mol Med. 2017;23(1):1. http://www.molmed.org/content/pdfstore/16_139_Barker.pdf. Accessed 4 Jul 2018.
Burel JG, Apte SH, Groves PL, Boyle MJ, Langer C, Beeson JG, et al. Dichotomous miR expression and immune responses following primary blood-stage malaria. JCI Insight. 2017;2(15). https://insight.jci.org/articles/view/93434. Accessed 4 Jul 2018.
Chamnanchanunt S, Fucharoen S, Umemura T. Circulating microRNAs in malaria infection: bench to bedside. Malaria J. 2017;16(1):334. https://doi.org/10.1186/s12936-017-1990-x.
Oturai DB, Søndergaard HB, Börnsen L, Sellebjerg F, Christensen JR. Identification of suitable reference genes for peripheral blood mononuclear cell subset studies in multiple sclerosis. Scand J Immunol. 2016;83(1):72–80. https://doi.org/10.1111/sji.12391.
Panichakul T, Ponnikorn S, Roytrakul S, Paemanee A, Kittisenachai S, Hongeng S, et al. Plasmodium vivax inhibits erythroid cell growth through altered phosphorylation of the cytoskeletal protein ezrin. Malaria J. 2015;14(1):138. http://www.malariajournal.com/content/14/1/138. Accessed 4 Jul 2018.
Wickramasinghe SN, Looareesuwan S, Nagachinta B, White NJ. Dyserythropoiesis and ineffective erythropoiesis in Plasmodium vivax malaria. Br J Haematol. 1989;72(1):91–9.
Panichakul T, Payuhakrit W, Panburana P, Wongborisuth C, Hongeng S, Udomsangpetch R. Suppression of erythroid development in vitro by Plasmodium vivax. Malaria J. 2012;11(1):173. https://doi.org/10.1186/1475-2875-11-173.
Horiguchi H, Kobune M, Kikuchi S, Yoshida M, Murata M, Murase K, et al. Extracellular vesicle miR-7977 is involved in hematopoietic dysfunction of mesenchymal stromal cells via poly(rC) binding protein 1 reduction in myeloid neoplasms. Haematologica. 2016;101(4):437–47. https://doi.org/10.3324/haematol.2015.134932.
Wahl SM, McCartney-Francis N, Mergenhagen SE. Inflammatory and immunomodulatory roles of TGF-beta. Immunol Today. 1989;10(8):258–61.
Omer FM, Riley EM. Transforming growth factor beta production is inversely correlated with severity of murine malaria infection. J Exp Med. 1998;188(1):39–48.
Chaves YO, da Costa AG, Pereira MLM, de Lacerda MVG, Coelho-Dos-Reis JG, Martins-Filho OA, et al. Immune response pattern in recurrent Plasmodium vivax malaria. Malaria J. 2016;15(1):445. https://doi.org/10.1186/s12936-016-1501-5.
Wenisch C, Parschalk B, Burgmann H, Looareesuwan S, Graninger W. Decreased serum levels of TGF-beta in patients with acute Plasmodium falciparum malaria. J Clin Immunol. 1995;15(2):69–73.
Perkins DJ, Weinberg JB, Kremsner PG. Reduced interleukin-12 and transforming growth factor-beta1 in severe childhood malaria: relationship of cytokine balance with disease severity. J Infect Dis. 2000;182(3):988–92. https://doi.org/10.1086/315762.
Coulson RMR, Hall N, Ouzounis CA. Comparative genomics of transcriptional control in the human malaria parasite Plasmodium falciparum. Genome Res. 2004;14(8):1548–54.
Hall N, Carlton J. Comparative genomics of malaria parasites. Curr Opin Genet Dev. 2005;15(6):609–13.
Kobayashi M, Oshima S, Maeyashiki C, Nibe Y, Otsubo K, Matsuzawa Y, et al. The ubiquitin hybrid gene UBA52 regulates ubiquitination of ribosome and sustains embryonic development. Scientific Rep. 2016;6(1):36780. http://www.nature.com/articles/srep36780. Accessed 4 Jul 2018.
Venny. Venny 2.1.0. http://bioinfogp.cnb.csic.es/tools/venny/. Accessed 3 Jul 2018.
HK and AK carried out the experiment, analysis and drafted the manuscript; other authors contributed to the critical revision of the manuscript. All authors read and approved the final manuscript.
We thank ICMR for providing Senior Research Fellowship to Hargobinder Kaur (No 80/883/2014-ECD-1). We thank Qatar National Research Fund (a member of Qatar Foundation) for providing partial financial assistance.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the present study results and additional files.
Consent for publication
Ethics approval and consent to participate
The present study was carried out in a tertiary care hospital PGIMER, North India. The ethical approval vide number Histo/14/2926 was obtained from the PGIMER’s ethical committee. The study protocol and research proposal and have been thoroughly reviewed by the committee members. Prior to sample collection the written informed consent from all the subjects was obtained.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differentially expressed miRNAs in 1. Complicated vs uncomplicated P. vivax 2. Complicated P. vivax vs Healthy 3. Complicated P. vivax vs P. falciparum with p-value < 0.05.
Summarized layout for the detailed analysis of PPI networks.
About this article
Cite this article
Kaur, H., Sehgal, R., Kumar, A. et al. Screening and identification of potential novel biomarker for diagnosis of complicated Plasmodium vivax malaria. J Transl Med 16, 272 (2018). https://doi.org/10.1186/s12967-018-1646-9
- Plasmodium vivax
- Complicated P. vivax malaria