Next Article in Journal
Differential Transcriptome Profiling Unveils Novel Deregulated Gene Signatures Involved in Pathogenesis of Alzheimer’s Disease
Previous Article in Journal
Checkpoint Inhibitors and Induction of Celiac Disease-like Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neural Network-Based Calculator for Rat Glomerular Filtration Rate

by
Óscar J. Pellicer-Valero
1,†,
Giampiero A. Massaro
2,3,4,5,6,†,
Alfredo G. Casanova
2,3,4,5,6,
María Paniagua-Sancho
2,3,4,5,6,
Isabel Fuentes-Calvo
2,3,5,6,
Mykola Harvat
1,
José D. Martín-Guerrero
1,7,‡,
Carlos Martínez-Salgado
2,3,5,6,7,‡ and
Francisco J. López-Hernández
2,3,4,5,6,7,8,*,‡
1
Intelligent Data Analysis Laboratory (IDAL), Department Electronic Engineering, School of Engineering (ETSE-UV), Universitat de València, 46100 Valencia, Spain
2
Institute of Biomedical Research of Salamanca, 37007 Salamanca, Spain
3
Departmento de Fisiología y Farmacología, Universidad de Salamanca, 37007 Salamanca, Spain
4
Fundación Instituto de Estudios de Ciencias de la Salud de Castilla y León, 42002 Soria, Spain
5
Group of Translational Research on Renal and Cardiovascular Diseases (TRECARD), 37007 Salamanca, Spain
6
National Network for Kidney Research REDINREN, RD016/0009/0025, Instituto de Salud Carlos III, 28029 Madrid, Spain
7
Disease and Theranostic Modelling Consortium (DisMOD), 37007 Salamanca, Spain
8
Group of Biomedical Research on Critical Care (BioCritic), 47003 Valladolid, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
Biomedicines 2022, 10(3), 610; https://doi.org/10.3390/biomedicines10030610
Submission received: 6 February 2022 / Revised: 2 March 2022 / Accepted: 3 March 2022 / Published: 5 March 2022
(This article belongs to the Section Molecular and Translational Medicine)

Abstract

:
Glomerular filtration is a pivotal process of renal physiology, and its alterations are a central pathological event in acute kidney injury and chronic kidney disease. Creatinine clearance (ClCr), a standard method for glomerular filtration rate (GFR) measurement, requires a long and tedious procedure of timed (usually 24 h) urine collection. We have developed a neural network (NN)-based calculator of rat ClCr from plasma creatinine (pCr) and body weight. For this purpose, matched pCr, weight, and ClCr trios from our historical records on male Wistar rats were used. When evaluated on the training (1165 trios), validation (389), and test sets (660), the model committed an average prediction error of 0.196, 0.178, and 0.203 mL/min and had a correlation coefficient of 0.863, 0.902, and 0.856, respectively. More importantly, for all datasets, the NN seemed especially effective at comparing ClCr among groups within individual experiments, providing results that were often more congruent than those measured experimentally. ACLARA, a friendly interface for this calculator, has been made publicly available to ease and expedite experimental procedures and to enhance animal welfare in alignment with the 3Rs principles by avoiding unnecessary stressing metabolic caging for individual urine collection.

Graphical Abstract

1. Introduction

Alterations in GFR are a hallmark of many renal ailments, including acute kidney injury (AKI) [1] and chronic kidney disease (CKD) [2,3], and a gold standard diagnostic parameter [4,5,6]. Methods to measure GFR are based on the clearance of specific probe molecules from the blood. These must be cleared solely by renal excretion, not metabolized, freely filtered at the glomerulus, or secreted or reabsorbed by the tubule, and must not interfere with the GFR [7,8]. Different exogenous molecules with these characteristics have been used [8,9,10], but the clearance of inulin (ClIn), a plant fructan oligosaccharide, is still the gold standard GFR measurement method [7,11,12,13].
These clearance methods are overtly inconvenient for routine clinical practice [14,15] due to the time-consuming and costly procedures, such as the need for exogenous probe administration, timed urine and blood collection and analysis, and the involvement of radioactivity [9]. For these reasons, GFR measurement turned to endogenous probes [16], the most used of which is creatinine (Cr), a metabolite of phosphocreatine and a waste product of muscle activity. Creatinine clearance (ClCr) is not the ideal method for GFR measurement either, because Cr is also secreted in the tubule and, in some circumstances, can also be cleared from the organism with the feces [17]. In fact, 10–40% of the creatinine found in the urine is secreted in normal individuals [18], which leads to an overestimation of the GFR when measured as ClCr, compared to ClIn [19]. Yet endogenous ClCr, introduced in 1937 [20], avoids probe administration and the use of radioactivity because Cr can be measured with a colorimetric reaction known as Jaffe’s method [7,21,22], or with enzymatic assays [23]. Its reduced complexity, cost, and risk rapidly led to the adoption of ClCr [7], which still requires timed urine sampling (usually 24-h collection).
Serum creatinine concentration (sCr) is a gross and inverse reflection of ClCr and, thus, of the GFR because it is produced in the muscle at a more or less constant average daily rate. Variations in sCr are therefore mostly derived from changes in ClCr. Accordingly, with all its limitations [24,25], sCr was progressively adopted as the election proxy to estimate GFR, especially after methods for automatic analysis were developed by the 1960s [21,25]. Still, sCr is dependent on muscle mass, race, and other factors; thus, algorithms have been developed to estimate GFR from sCr and personal and body habitus data [9]. In 1976, the Schwartz formula (adjusted by height) [26] and the Cockroft–Gault equation (adjusted by age and weight) [27] were published. Over 25 formulas have been developed thereafter [28], most notably including the Modification of Diet in Renal Disease equation (MDRD) [29] and the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) group equation [30], both adjusted by age, sex, and ethnicity. Due to its reasonably high accuracy and practicality, estimated GFR has substituted measured GFR (i.e., clearance methods) in most clinical situations [4]. Recently, a novel approach has been proposed to efficiently estimate GFR from sCr through a numerical algorithm based on bivariate numerical modeling [31].
Similarly, in experimental research with animals, ClIn [32,33] is the preferred method for GFR assessment, but endogenous ClCr [34,35] is most widely used due to reduced complexity. Still, individual, laborious 24-hour urine collection is necessary. A transcutaneous exogenous fluorescence-labeled sinistrin clearance (ClSin) method has been developed which avoids individual housing and allows longitudinal studies with minimal invasiveness [36]. However, ClSin is also laborious and technically more demanding than traditional clearance techniques. Due to simplicity, sCr is again the most frequent surrogate for GFR.
In this study we developed and experimentally validated a neural network-based calculator of ClCr for rats from a single analyte (i.e., sCr) and body weight data, useful in both healthy and disease conditions. This calculator may be used in many acute and chronic pathophysiological circumstances to prevent long, tedious, and costly procedures and avoid the stressful conditions for experimental animals associated with ClCr measurement.

2. Materials and Methods

2.1. Data Mining and Database Generation

To train the model, 1554 data trios of measured ClCr (mClCr), pCr, and body weight were collected from the historical records of our laboratory. These datasets came from published and unpublished experiments with AKI and CKD models (and their healthy, normal controls) (Figure 1 and Table 1). A random subset with 75% of the data trios (1165) was used for model training (i.e., the training set), and the remaining 25% (389) was used for validation (i.e., the validation set), which was utilized for making design decisions (e.g., model selection, preprocessing pipeline design, etc.). Once the final model was obtained, it was field-tested against a third dataset with 660 data trios (i.e., the test set) from rats and experiments different from those used for training and validation. In all cases, the data arose from male Wistar rats, with no exclusion criteria. ClCr had been measured as we previously described [37,38]. Briefly, ClCr had been determined with the following formula: ClCr = UF24h × uCr/pCr, where UF24h is the volume of urine collected in 24 h, pCr is the plasma concentration of creatinine, and uCr is the urinary concentration of creatinine. For individual urine collection, rats had been allocated in metabolic cages and allowed free access to water and regular chow. Blood samples had been collected from a small incision in the tail and immediately centrifuged to obtain the plasma, which was frozen at −80 °C until use. Both pCr and uCr had been measured with a commercial colorimetric assay based on Jaffe’s method (Quantichrom Creatinine Assay Kits, BioAssay Systems, Hayward, CA, USA) following the manufacturer’s instructions. A database was generated with rat identifiers, experimental model, ClCr, pCr, and body weight values.

2.2. Algorithm Description

Linear regression (LR), random forest (RF), and feedforward neural network (FFNN) [43,44] models were evaluated comparatively. Among them, LR is the simplest algorithm as it can only learn independent linear relationships between each input feature and the output. Conversely, RF is a non-linear model able to capture complex relationships among input and output features as it uses the compound prediction from many (e.g., 100) decision trees to perform regression. Finally, FFNNs are connectionist models that allow for fitting non-linear relationships by stacking many LR units (i.e., neurons) with non-linear activation functions in between, similar to how simple biological neurons combine to form a complex net: the brain. For the LR and RF models, the Scikit-learn library (version 0.24.1) [45] was employed, while for the FFNN, TensorFlow (version 2.1.0) [46] was used.

2.3. Data Preprocessing

The relationship between pCr and mClCr was linearized by logarithmic transformation (Equation (1) and Supplementary Figure S1) to help the ML algorithms better capture the relationship. The input features (pCr and weight) were then standardized by Equation (2), which centers the data distribution around zero and scales it to unit variance. Standardization is required by LR, and it is helpful for training other ML models too since it equalizes the relative importance of the features.
f ( x ) = log ( 1 x )
x ˜ = x m e a n ( x ) s t d ( x )

2.4. Neural Network Training, Validation, and Testing

For FFNN training, the parameters, i.e., weights, minimizing the model’s prediction error, i.e., the mean squared error (MSE, Equation (3)), between eClCr and mClCr in the training set were found by the stochastic gradient descent (SGD) [47] iterative optimizer.
M S E ( e C l C r ,   m C l C r ) = m e a n ( ( e C l C r m C l C r ) 2 )
Model performance was mainly assessed with the mean absolute error (MAE) and the Pearson product-moment correlation coefficient (Equations (4) and (5), respectively), although p10 and p30 metrics were also included for completeness. The MAE represents the absolute error (in mL/min) that is committed on average in every prediction. The correlation coefficient measures the strength of the relationship between mCl_Cr and eCl_Cr with a value between 0 and 1, where 1 stands for a perfect prediction performance (i.e., eClCr = mClCr for all data trios). Finally, p10 and p30 represent the fraction of predictions (eCl_Cr) that fall within 10% and 30% of the mCl_Cr, respectively.
M A E ( e C l C r ,   m C l C r ) = m e a n ( | e C l C r m C l C r | )
c o r r ( e C l C r ,   m C l C r ) = c o v ( e C l C r ,   m C l C r ) s t d (   m C l C r ) · s t d ( e C l C r )
To avoid potential overfitting, all engineering decisions (FFNN architecture, number of layers, neurons per layer, activation functions, etc.) were made based on the validation MAE and correlation coefficient.
Once the model (NN architecture + weights) was finalized, it was challenged with the test set as a proxy for field performance.

3. Results

3.1. Model Development, Evaluation, and Validation

With the objective of developing a machine learning (ML) model (i.e., a calculator) to estimate rat ClCr (i.e., estimated ClCr; eClCr) from plasma creatinine concentration (pCr) and body weight, three ML algorithms were assessed and compared: LR, RF, and FFNN. Table 2 shows the performance metrics for the three models after training and tuning them. The FFNN outperformed the LR and RF models on the validation set and was therefore selected for further testing. Values of eClCr were compared to the measured (i.e., experimental) values of ClCr (mClCr).
The optimal FFNN architecture was found to have three intermediate (hidden) layers with 30, 20, and 10 neurons respectively, with the hyperbolic tangent as an activation function and using the Adam SGD optimizer with a learning rate of 0.001 and a batch size of 400 with MSE as the loss function. With the aid of the early stopping technique, training was stopped after 200 SGD iterations without improvement of the validation loss. The characteristics of the final FFNN architecture are shown in Supplementary Figure S2 and a decision surface map for the FFNN is shown in Supplementary Figure S3. Agreement between eClCr and mClCr was also studied in the test set by means of a Bland–Altman plot (Supplementary Figure S4).
In addition to analyzing the global performance of the model, the results of eClCr were also compared with their corresponding mClCr values in the context of individual experiments. Figure 2 shows the mClCr and eClCr for four representative experiments (from the training and validation datasets). In general terms, the information provided by measured ClCr (mClCr) and eClCr is almost identical. In particular, in those cases in which mClCr and eClCr are less similar (e.g., groups D2.1 and D2.5 in Figure 2), eClCr is consistently more congruent within its context. Just as an example, the group D2.1 is a control, and groups D2.2 to D2.5 refer to animals treated with cisplatin in several conditions. There is no reason for mClCr to increase in group D2.5 on day 5. On the contrary, following a similar pattern to the other AKI groups, eClCr shows a more logical result for group D2.5 and overall.

3.2. Model Field Testing

As shown in Table 2, the MAE and Pearson product-moment correlation coefficient for the test set were 0.2035 and 0.8564, respectively, which are similar to those obtained in the training and validation sets, hence proving that the model is not overfitted and can generalize to unseen data. Furthermore, a paired t-test was employed to compare the mean of the absolute error committed by the NN and the LR models in the test set; the NN was found to be superior with very high statistical significance (p < 0.001, using SciPy library [48], version 1.5.2). As above, eClCr predictions were compared with their corresponding mClCr values in the context of specific experiments. Figure 3 shows the mClCr and eClCr for four illustrative experiments (from the test data). There is a good match at the whole experiment level between mClCr and eClCr. Again, eClCr is more congruent than mClCr. As exemplifying tokens, the mClCr in group T1.5 by day 2 and group T1.1 by day 8 (a control) shows illogical increments that are normalized or softened by the corresponding eClCr. In experiment T2, the open circles correspond to the control group and the filled circles to a highly nephrotoxic, one-dose treatment. As expected, by day 4, there is a notorious drop in the GFR, which is progressively restored afterwards. The gradual recovery described by eClCr contrasts with the more erratic trajectory depicted by mClCr.
The influence of weight and mClCr on the distribution of the prediction error (eClCr-mClCr) in the whole test set is shown in Figure 4. The variance of the prediction error grows (top panel) with the weight, probably as a result of an increasingly lower amount of data for higher weights. The medians suggest that the model may underestimate the ClCr for high weights. In contrast, the variance of the prediction error seems to be similar for all ClCr values (bottom panel). Yet the medians show that the model tends to slightly overestimate ClCr for small to medium mClCr values (below 0.683 mL/min) and to slightly underestimate ClCr for medium to larger mClCr values (above 1.138 mL/min).

4. Discussion

Our study reveals that an FFNN fed with pCr and weight values as inputs predicts ClCr (i.e., eClCr) with reasonably high accuracy, regardless of additional variables such as age, comorbidities, or treatment. This calculator may be useful for the estimation of ClCr in those laboratories where metabolic cages are not available, or when avoidance of experimental complexity and resource optimization are sought.
When evaluated with the training, validation, and test datasets, the model showed a relatively low average error. Furthermore, the model proved especially useful for evaluating experiments on the whole. In fact, the behavior and evolution of eClCr for each group was, in general, more homogeneous and congruent than the behavior of mClCr, while the relations among experimental groups resulting from eClCr remained extremely similar to those obtained with mClCr. Our interpretation of this data-equalization effect is that, as opposed to mClCr, eClCr obviates the error introduced by urine collection in metabolic cages, which is unavoidably translated to uCr estimation and, in turn, to mClCr, and thus sometimes yields unexplainable and illogical fluctuations (as argued in the Results section). This advantage adds to the procedural simplicity of eClCr compared to mClCr. Accordingly, this calculator is not aimed at being a method for the precise determination of ClCr but a practical tool providing rapid, easy, and sufficiently accurate results for many purposes of basic research with rats, in the same way that human GFR-estimating formulas are useful in the clinical setting. In our view, this calculator is particularly appropriate for those experiments in which it is not the exact individual value of ClCr that is required but rather its comparison among experimental groups or time points. In addition to facilitating preclinical research, eClCr also provides an alternative to mClCr that improves experimental animal welfare. Individual confinement in metabolic cages [37,38] required for urine collection is known to stress social animals such as rats [49]. Hence, this method is tightly aligned with the 3Rs principles (replacement, reduction, and refinement) concerning animal research [50].
In rats, pCr is mostly constant throughout their lifespan. However, the GFR (and thus ClCr) increases within an initial age range, after which it largely stabilizes [51]. Therefore, age should be taken into account when modeling the ClCr–pCr relationship. Yet because rat age might not be precisely known by the researcher in some cases, we used weight as a surrogate and found that its inclusion significantly improved model performance (even more than age itself). Furthermore, we contend that weight recapitulates other potential sources of disparity related to body habitus. Ethnicity and sex were not included in the model as we did not have sufficient data from other rat strains, nor of female rats. Accordingly, this model is adjusted only for male Wistar rats. Its performance must also be tested, and the model adjusted (if and as necessary), for different rat strains and female sex.
Another important factor potentially uncoupling the relationship between pCr and ClCr is creatinine secretion. Differential affecting of creatinine secretion by different experimental conditions might increase variability and hamper model performance. However, the model was generated on heterogeneous data from etiologically diverse disease models, and there were no deviation tendencies to specific characteristics. This implies that creatinine secretion is not significantly modified by the treatments or conditions included and that the observed data variability mainly reflects experimental variability derived from sample collection and analytic error.

5. Conclusions

The proposed model is a reasonably accurate surrogate for mClCr, solves some of its faults, and even provides more congruent results. eGFR in humans is known to be less accurate in specific populations and circumstances, such as individuals with abnormal muscle mass or body surface area; during changes in metabolism and pregnancy; and when GFR is rapidly changing, such as during growth, AKI, or high-protein ingestion [16,52]. In contrast, laboratory rats present a lower phenotypic heterogeneity as compared to humans, so a higher accuracy is therefore expected. Lastly, we have released ACLARA, a web-based, openly available eClCr calculator (http://idal.uv.es/aclara, accessed on 5 February 2022), to facilitate its use and encourage further research on the topic.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10030610/s1, Figure S1: Data transformation. Figure S2: Graphical description of the FFNN architecture. Figure S3: FFNN decision surface map. Figure S4: Bland–Altman plot of agreement between eClCr and mClCr in the test set.

Author Contributions

G.A.M., A.G.C., M.P.-S. and I.F.-C. carried out data collection and databases. Ó.J.P.-V., M.H. and J.D.M.-G. developed the mathematical model. A.G.C., Ó.J.P.-V., G.A.M. and M.P.-S. made the figures and tables. C.M.-S., J.D.M.-G. and F.J.L.-H. designed the study. Ó.J.P.-V., C.M.-S., J.D.M.-G. and F.J.L.-H. drafted the manuscript. All authors critically revised and approved the final text. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the Government of Spain (Instituto de Salud Carlos III) (PI18/00996, PI21/01226, and Retic RD021/0005/0004, REDINREN), co-funded by FEDER funds, and from the Ministry of Education, Government of Castile and León (IES160P20), co-funded by FEDER funds.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Basile, D.P.; Anderson, M.D.; Sutton, T.A. Pathophysiology of acute kidney injury. Compr. Physiol. 2012, 2, 1303–1353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. López-Novoa, J.M.; Martínez-Salgado, C.; Rodríguez-Peña, A.B.; Hernández, F.J.L. Common pathophysiological mechanisms of chronic kidney disease: Therapeutic perspectives. Pharmacol. Ther. 2010, 128, 61–81. [Google Scholar] [CrossRef] [PubMed]
  3. López-Novoa, J.M.; Rodríguez-Peña, A.B.; Ortiz, A.; Martínez-Salgado, C.; López Hernández, F.J. Etiopathology of chronic tubular, glomerular and renovascular nephropathies: Clinical implications. J. Transl. Med. 2011, 9, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Murphy, D.P.; Hsu, C.Y. Estimating glomerular filtration rate: Is it good enough? and is it time to move on? Curr. Opin. Nephrol. Hypertens. 2013, 22, 310–315. [Google Scholar] [CrossRef] [PubMed]
  5. Levey, A.S.; Becker, C.; Inker, L.A. Glomerular filtration rate and albuminuria for detection and staging of acute and chronic kidney disease in adults: A systematic review. JAMA J. Am. Med. Assoc. 2015, 313, 837–846. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. KDIGO. Clinical Practice Guideline for Acute Kidney Injury. Kidney Int. Suppl. 2012. Available online: https://kdigo.org/guidelines/ckd-evaluation-and-management/ (accessed on 5 February 2022).
  7. Diskin, C.J. Creatinine and glomerular filtration rate: Evolution of an accommodation. Ann. Clin. Biochem. 2007, 44, 16–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Aurell, M. Accurate and feasible measurements of GFR-is the iohexol clearance the answer? Nephrol. Dial. Transplant. 1994, 9, 1222–1224. [Google Scholar] [CrossRef] [PubMed]
  9. Schwartz, G.J.; Furth, S.L. Glomerular filtration rate measurement and estimation in chronic kidney disease. Pediatr. Nephrol. 2007, 22, 1839–1848. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Delanaye, P.; Ebert, N.; Melsom, T.; Gaspari, F.; Mariat, C.; Cavalier, E.; Björk, J.; Christensson, A.; Nyman, U.; Porrini, E.; et al. Iohexol plasma clearance for measuring glomerular filtration rate in clinical practice and research: A review. Part 1: How to measure glomerular filtration rate with iohexol? Clin. Kidney J. 2016, 9, 682–699. [Google Scholar] [CrossRef] [PubMed]
  11. Smith, H. The inulin clearance. In The Kidney; Smith, H., Ed.; Oxford University Press: New York, NY, USA, 1951; pp. 47–55. [Google Scholar]
  12. Filler, G.; Yasin, A.; Medeiros, M. Methods of assessing renal function. Pediatr. Nephrol. 2014, 29, 183–192. [Google Scholar] [CrossRef] [PubMed]
  13. Soveri, I.; Berg, U.B.; Björk, J.; Elinder, C.G.; Grubb, A.; Mejare, I.; Sterner, G.; Bäck, S.E. Measuring GFR: A systematic review. Am. J. Kidney Dis. 2014, 64, 411–424. [Google Scholar] [CrossRef] [PubMed]
  14. Inker, L.A.; Levey, A.S.; Coresh, J. Estimated Glomerular Filtration Rate from a Panel of Filtration Markers—Hope for Increased Accuracy Beyond Measured Glomerular Filtration Rate? Adv. Chronic Kidney Dis. 2018, 25, 67–75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Prasad, N.; Barai, S.; Gambhir, S.; Parasar, D.; Ora, M.; Gupta, A.; Sharma, R. Comparison of glomerular filtration rate estimated by plasma clearance method with modification of diet in renal disease prediction equation and Gates method. Indian J. Nephrol. 2012, 22, 103–107. [Google Scholar] [CrossRef]
  16. Rizk, D.V.; Meier, D.; Sandoval, R.M.; Chacana, T.; Reilly, E.S.; Seegmiller, J.C.; DeNoia, E.; Strickland, J.S.; Muldoon, J.; Molitoris, B.A. A novel method for rapid bedside measurement of GFR. J. Am. Soc. Nephrol. 2018, 29, 1609–1613. [Google Scholar] [CrossRef] [PubMed]
  17. Wixom, R.L.; Davis, G.E.; Flynn, M.A.; Tsutakawa, R.T.; Hentges, D.J. Excretion of Creatine and Creatinine in Feces of Man. Proc. Soc. Exp. Biol. Med. 1979, 161, 452–457. [Google Scholar] [CrossRef] [PubMed]
  18. Levey, A.S.; Perrone, R.D.; Madias, N.E. Serum creatinine and renal function. Annu. Rev. Med. 1988, 39, 465–490. [Google Scholar] [CrossRef]
  19. Zhang, X.; Rule, A.D.; McCulloch, C.E.; Lieske, J.C.; Ku, E.; Hsu, C.Y. Tubular secretion of creatinine and kidney function: An observational study. BMC Nephrol. 2020, 21, 108. [Google Scholar] [CrossRef]
  20. Popper, H.; Mandel, E. Filtrations and reabsorptions leitung in dernierenpathologie. Erg. Inn. Med. Kinder 1937, 53, 685–695. [Google Scholar]
  21. Chasson, A.; Grady, H.; Stanley, M. Determination of creatinine by means of automatic chemical analysis. Tech. Bull. Regist Med. Technol. 1960, 30, 207–212. [Google Scholar] [CrossRef]
  22. Husdan, H.; Rapoport, A. Estimation of creatinine by the Jaffe reaction. A comparison of three methods. Clin. Chem. 1968, 14, 222–238. [Google Scholar] [CrossRef]
  23. Brück, K.; Jager, K.J.; Dounousi, E.; Kainz, A.; Nitsch, D.; Ärnlöv, J.; Rothenbacher, D.; Browne, G.; Capuano, V.; Ferraro, P.M.; et al. Methodology used in studies reporting chronic kidney disease prevalence: A systematic literature review. Nephrol. Dial. Transplant. 2015, 30, iv6–iv16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Delanaye, P.; Cavalier, E.; Pottel, H. Serum Creatinine: Not so Simple! Nephron 2017, 136, 302–308. [Google Scholar] [CrossRef] [PubMed]
  25. Dennen, P.; Parikh, C.R. Biomarkers of acute kidney injury: Can we replace serum creatinine? Clin. Nephrol. 2007, 68, 269–278. [Google Scholar] [CrossRef] [PubMed]
  26. Schwartz, G.; Haycock, G.; Edelmann, C., Jr.; Spitzer, A. A simple estimate of glomerular filtration rate in children derived from body length and plasma creatinine. Pediatrics 1976, 58, 259–263. [Google Scholar] [CrossRef]
  27. Cockcroft, D.W.; Gault, M.H. Prediction of creatinine clearance from serum creatinine. Nephron 1976, 16, 31–41. [Google Scholar] [CrossRef]
  28. Lamb, E.J.; Tomson, C.R.V.; Roderick, P.J. Estimating kidney function in adults using formulae. Ann. Clin. Biochem. 2005, 42, 321–345. [Google Scholar] [CrossRef]
  29. Levey, A.S.; Bosch, J.P.; Lewis, J.B.; Greene, T.; Rogers, N.; Roth, D. A more accurate method to estimate glomerular filtration rate from serum creatinine: A new prediction equation. Ann. Intern. Med. 1999, 130, 461–470. [Google Scholar] [CrossRef]
  30. Levey, A.S.; Stevens, L.A.; Schmid, C.H.; Zhang, Y.; Castro, A.F.; Feldman, H.I.; Kusek, J.W.; Eggers, P.; Van Lente, F.; Greene, T.; et al. A new equation to estimate glomerular filtration rate. Ann. Intern. Med. 2009, 150, 604–612. [Google Scholar] [CrossRef]
  31. Giles, S.N.; Fiori, S. Glomerular Filtration Rate Estimation by a Novel Numerical Binning-Less Isotonic Statistical Bivariate Numerical Modeling Method. Information 2019, 10, 100. [Google Scholar] [CrossRef] [Green Version]
  32. Passos, M.T.; Nishida, S.K.; Câmara, N.O.S.; Shimizu, M.H.; Mastroianni-Kirsztajn, G. Iohexol clearance for determination of glomerular filtration rate in rats induced to acute renal failure. PLoS ONE 2015, 10, e0123753. [Google Scholar] [CrossRef] [Green Version]
  33. Zhu, X.R.; Ge, N.; Wang, Y.; Zhai, J.L.; Liu, C. Corticosteroids significantly increase cystatin C levels in the plasma by promoting cystatin C production in rats. Ren. Fail. 2019, 41, 698–703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zager, R.A. Exogenous Creatinine Clearance Accurately Assesses Filtration Failure in Rat Experimental Nephropathies. Am. J. Kidney Dis. 1987, 10, 427–430. [Google Scholar] [CrossRef]
  35. Van Liew, J.B.; Zamlauski-Tucker, M.J.; Feld, L.G. Endogenous creatinine clearance in the rat: Strain variation. Life Sci. 1993, 53, 1015–1021. [Google Scholar] [CrossRef]
  36. Ellery, S.J.; Cai, X.; Walker, D.D.; Dickinson, H.; Kett, M.M. Transcutaneous measurement of glomerular filtration rate in small rodents: Through the skin for the win? Nephrology 2015, 20, 117–123. [Google Scholar] [CrossRef]
  37. Sancho-Martínez, S.M.; Blanco-Gozalo, V.; Quiros, Y.; Prieto-García, L.; Montero-Gómez, M.J.; Docherty, N.G.; Martínez-Salgado, C.; Morales, A.I.; López-Novoa, J.M.; López-Hernández, F.J. Impaired Tubular Reabsorption Is the Main Mechanism Explaining Increases in Urinary NGAL Excretion Following Acute Kidney Injury in Rats. Toxicol. Sci. 2020, 175, 75–86. [Google Scholar] [CrossRef] [PubMed]
  38. Sancho-Martínez, S.M.; Sánchez-Juanes, F.; Blanco-Gozalo, V.; Fontecha-Barriuso, M.; Prieto-García, L.; Fuentes-Calvo, I.; González-Buitrago, J.M.; Morales, A.I.; Martínez-Salgado, C.; Ramos-Barron, M.A.; et al. Urinary TCP1-eta: A Cortical Damage Marker for the Pathophysiological Diagnosis and Prognosis of Acute Kidney Injury. Toxicol. Sci. 2020, 174, 3–15. [Google Scholar] [CrossRef]
  39. Quiros, Y.; Ferreira, L.; Sancho-Martínez, S.M.; González-Buitrago, J.M.; López-Novoa, J.M.; López-Hernández, F.J. Sub-nephrotoxic doses of gentamicin predispose animals to developing acute kidney injury and to excrete ganglioside M2 activator protein. Kidney Int. 2010, 78, 1006–1015. [Google Scholar] [CrossRef] [Green Version]
  40. Casanova, A.G.; Vicente-Vicente, L.; Hernández-Sánchez, M.T.; Prieto, M.; Rihuete, M.I.; Ramis, L.M.; del Barco, E.; Cruz, J.J.; Ortiz, A.; Cruz-González, I.; et al. Urinary transferrin pre-emptively identifies the risk of renal damage posed by subclinical tubular alterations. Biomed. Pharmacother. 2020, 121, 109684. [Google Scholar] [CrossRef]
  41. Vicente-Vicente, L.; Sánchez-Juanes, F.; García-Sánchez, O.; Blanco-Gozalo, V.; Pescador, M.; Sevilla, M.A.; González-Buitrago, J.M.; López-Hernández, F.J.; López-Novoa, J.M.; Morales, A.I. Sub-nephrotoxic cisplatin sensitizes rats to acute renal failure and increases urinary excretion of fumarylacetoacetase. Toxicol. Lett. 2015, 234, 99–109. [Google Scholar] [CrossRef]
  42. Casanova, A.G.; Fuentes-Calvo, I.; Hernández-Sánchez, M.T.; Quintero, M.; Toral, P.; Caballero, M.T.; Martínez-Salgado, C.; Morales, A.I.; Layton, A.T.; Eleno, N.; et al. The furosemide stress test and computational modeling identify renal damage sites associated with predisposition to acute kidney injury in rats. Transl. Res. 2020, 231, 76–91. [Google Scholar] [CrossRef]
  43. Alpaydin, E. Introduction to Machine Learning, 4th ed.; MT Press: Cambridge, MA, USA, 2020. [Google Scholar]
  44. Breinman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  45. Pedregosa, F.; Michel, V.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Vanderplas, J.; Cournapeau, D.; Pedregosa, F.; Varoquaux, G.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  46. Abadi, M.; Agarwal, A.; Barham, P.; Brevdo, E.; Chen, Z.; Citro, C.; Corrado, G.S.; Davis, A.; Dean, J.; Devin, M.; et al. TensorFlow: Large-Scale Machine Learning on Heterogeneous Distributed Systems. arXiv 2016, arXiv:1603.04467. [Google Scholar]
  47. Kiefer, J.; Wolfowitz, J. Stochastic Estimation of the Maximum of a Regression Function. Ann. Math. Stat. 1952, 23, 462–466. [Google Scholar] [CrossRef]
  48. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
  49. Whittaker, A.L.; Lymn, K.A.; Howarth, G.S. Effects of Metabolic Cage Housing on Rat Behavior and Performance in the Social Interaction Test. J. Appl. Anim. Welf. Sci. 2016, 19, 363–374. [Google Scholar] [CrossRef]
  50. Sneddon, L.U.; Halsey, L.G.; Bury, N.R. Considering aspects of the 3Rs principles within experimental animal biology. J. Exp. Biol. 2017, 220, 3007–3016. [Google Scholar] [CrossRef] [Green Version]
  51. Potter, D.; Jarrah, A.; Sakai, T.; Harrah, J.; Holliday, M.A. Character of function and size in kidney during normal growth of rats. Pediatr. Res. 1969, 3, 51–59. [Google Scholar] [CrossRef] [Green Version]
  52. Filler, G.; Lee, M. Educational review: Measurement of GFR in special populations. Pediatr. Nephrol. 2018, 33, 2037–2046. [Google Scholar] [CrossRef]
Figure 1. Summary of the datasets used for model development and field testing. Data trios were composed of matched pCr, ClCr, and body weight data.
Figure 1. Summary of the datasets used for model development and field testing. Data trios were composed of matched pCr, ClCr, and body weight data.
Biomedicines 10 00610 g001
Figure 2. Comparison of the measured ClCr and estimated ClCr, corresponding to four exemplifying actual experiments using data from the training and validation sets. ClCr, creatinine clearance. CORR, Pearson product-moment correlation coefficient. MAE, mean average error.
Figure 2. Comparison of the measured ClCr and estimated ClCr, corresponding to four exemplifying actual experiments using data from the training and validation sets. ClCr, creatinine clearance. CORR, Pearson product-moment correlation coefficient. MAE, mean average error.
Biomedicines 10 00610 g002
Figure 3. Comparison of the measured ClCr and estimated ClCr, corresponding to four exemplifying actual experiments, using data from the test set. ClCr, creatinine clearance. CORR, Pearson product-moment correlation coefficient. MAE, mean average error.
Figure 3. Comparison of the measured ClCr and estimated ClCr, corresponding to four exemplifying actual experiments, using data from the test set. ClCr, creatinine clearance. CORR, Pearson product-moment correlation coefficient. MAE, mean average error.
Biomedicines 10 00610 g003
Figure 4. Prediction error analysis. Prediction error box plots according to the rat’s body weight range (top) and the mClCr range (bottom). Error bars represent the 2.5th and 97.5th percentiles. A scatter plot of prediction errors is also shown superimposed.
Figure 4. Prediction error analysis. Prediction error box plots according to the rat’s body weight range (top) and the mClCr range (bottom). Error bars represent the 2.5th and 97.5th percentiles. A scatter plot of prediction errors is also shown superimposed.
Biomedicines 10 00610 g004
Table 1. Characteristics of the datasets used for model development and field testing. CTRL, control. T-AKI, toxic acute kidney injury. I/R, ischemia/reperfusion. M, male. RMR, 5/6 renal mass reduction. PRD, predisposed-to-AKI (as in [39,40,41,42]).
Table 1. Characteristics of the datasets used for model development and field testing. CTRL, control. T-AKI, toxic acute kidney injury. I/R, ischemia/reperfusion. M, male. RMR, 5/6 renal mass reduction. PRD, predisposed-to-AKI (as in [39,40,41,42]).
Modeln
(Rats)
Age
at Start
(Weeks)
Duration
(Weeks)
n
(Data Trios)
Model Development (Training + Validation Datasets)
CTRL98151
3482–28194
T-AKI10981583
1982–898
6482–28381
I/R1582–863
RMR3582–28184
Test Dataset
CTRL138159
28510
38721
PRD6181279
T-AKI98538
78749
5081204
Table 2. Performance metrics for the three considered models in the training and validation sets. LR, linear regression. RF, random forest. FFNN, feed forward neural network. MAE, mean average error. Correlation, Pearson product-moment correlation coefficient. P10/P30, fraction of predictions with an error within the 10%/30% threshold. Best results for each set are highlighted in bold.
Table 2. Performance metrics for the three considered models in the training and validation sets. LR, linear regression. RF, random forest. FFNN, feed forward neural network. MAE, mean average error. Correlation, Pearson product-moment correlation coefficient. P10/P30, fraction of predictions with an error within the 10%/30% threshold. Best results for each set are highlighted in bold.
Training SetValidation SetTest Set
LRRFFFNNLRRFFFNNLRRFFFNN
MAE0.22100.18010.19560.19750.18940.17800.25150.22570.2035
Correlation0.84020.88690.86320.87320.88710.90280.79220.74500.8563
P100.26350.33220.30990.26990.28280.28790.25300.28030.3136
P300.69010.77510.73910.70440.71890.75320.63940.67270.6955
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pellicer-Valero, Ó.J.; Massaro, G.A.; Casanova, A.G.; Paniagua-Sancho, M.; Fuentes-Calvo, I.; Harvat, M.; Martín-Guerrero, J.D.; Martínez-Salgado, C.; López-Hernández, F.J. Neural Network-Based Calculator for Rat Glomerular Filtration Rate. Biomedicines 2022, 10, 610. https://doi.org/10.3390/biomedicines10030610

AMA Style

Pellicer-Valero ÓJ, Massaro GA, Casanova AG, Paniagua-Sancho M, Fuentes-Calvo I, Harvat M, Martín-Guerrero JD, Martínez-Salgado C, López-Hernández FJ. Neural Network-Based Calculator for Rat Glomerular Filtration Rate. Biomedicines. 2022; 10(3):610. https://doi.org/10.3390/biomedicines10030610

Chicago/Turabian Style

Pellicer-Valero, Óscar J., Giampiero A. Massaro, Alfredo G. Casanova, María Paniagua-Sancho, Isabel Fuentes-Calvo, Mykola Harvat, José D. Martín-Guerrero, Carlos Martínez-Salgado, and Francisco J. López-Hernández. 2022. "Neural Network-Based Calculator for Rat Glomerular Filtration Rate" Biomedicines 10, no. 3: 610. https://doi.org/10.3390/biomedicines10030610

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop