- 1Department of Veterinary Medical Sciences, Alma Mater Studiorum-University of Bologna, Ozzano Emilia, Italy
- 2Department of Veterinary Sciences, University of Parma, Parma, Italy
- 3Department of Statistical Sciences, Alma Mater Studiorum-University of Bologna, Bologna, Italy
Introduction: Hypothyroidism can be easily misdiagnosed in dogs, and prediction models can support clinical decision-making, avoiding unnecessary testing and treatment. The aim of this study is to develop and internally validate diagnostic prediction models for hypothyroidism in dogs by applying machine-learning algorithms.
Methods: A single-institutional cross-sectional study was designed searching the electronic database of a Veterinary Teaching Hospital for dogs tested for hypothyroidism. Hypothyroidism was diagnosed based on suggestive clinical signs and thyroid function tests. Dogs were excluded if medical records were incomplete or a definitive diagnosis was lacking. Predictors identified after data processing were dermatological signs, alopecia, lethargy, hematocrit, serum concentrations of cholesterol, creatinine, total thyroxine (tT4), and thyrotropin (cTSH). Four models were created by combining clinical signs and clinicopathological variables expressed as quantitative (models 1 and 2) and qualitative variables (models 3 and 4). Models 2 and 4 included tT4 and cTSH, models 1 and 3 did not. Six different algorithms were applied to each model. Internal validation was performed using a 10-fold cross-validation. Apparent performance was evaluated by calculating the area under the receiver operating characteristic curve (AUROC).
Results: Eighty-two hypothyroid and 233 euthyroid client-owned dogs were included. The best performing algorithms were naive Bayes in model 1 (AUROC = 0.85; 95% confidence interval [CI] = 0.83–0.86) and in model 2 (AUROC = 0.98; 95% CI = 0.97–0.99), logistic regression in model 3 (AUROC = 0.88; 95% CI = 0.86–0.89), and random forest in model 4 (AUROC = 0.99; 95% CI = 0.98–0.99). Positive predictive value was 0.76, 0.84, 0.93, and 0.97 in model 1, 2, 3, and 4, respectively. Negative predictive value was 0.89, 0.89, 0.99, and 0.99 in model 1, 2, 3, and 4, respectively.
Discussion: Machine learning-based prediction models were accurate in predicting and quantifying the likelihood of hypothyroidism in dogs based on internal validation performed in a single-institution, but external validation is required to support the clinical applicability of these models.
Introduction
Hypothyroidism is considered a common endocrine disease in dogs, although its actual prevalence remains unknown. The most typical clinical signs include dermatological abnormalities (e.g., alopecia, poor-quality coat, pyoderma, seborrhea) and non-specific metabolic signs (e.g., decreased appetite, weight gain, lethargy, asthenia, and cold intolerance) (1). The most common clinicopathological abnormalities are hypercholesterolemia and mild to moderate normocytic normochromic non-regenerative anemia. Hypothyroid dogs less commonly show azotemia, increased liver enzymes and creatine kinase activity, and increased serum fructosamine concentrations (1–4). The diagnosis is easily confirmed in dogs with a clinical suspicion of hypothyroidism and concurrent low serum total thyroxine (tT4) concentration and high serum thyroid-stimulating hormone (TSH) concentrations. However, serum tT4 can be decreased in euthyroid dogs due to non-thyroidal illness syndrome (NTIS) or previous/ongoing drug treatment (e.g., glucocorticoids, phenobarbital, sulphonamides, tyrosine kinase inhibitors) (5–10). Serum TSH can be normal in 20 to 40% of hypothyroid dogs (11, 12). Moreover, breed-specific reference intervals have been reported in different breeds (e.g., Greyhounds and other Sighthounds, Basenji, Dogue de Bordeaux) (13–16). For all these reasons, the correct assessment of thyroid function in dogs can often be challenging, possibly leading to overdiagnosis of hypothyroidism. The vast majority of hypothyroid dogs obtain a complete resolution of clinical signs and clinicopathological abnormalities with appropriate levothyroxine supplementation. On the contrary, euthyroid dogs misdiagnosed with hypothyroidism receive inappropriate treatment which fails to improve clinical signs and leads to a possible delay in achieving the correct diagnosis, unjustified costs for the owner, and the risk of iatrogenic hyperthyroidism. When in doubt, additional testing to confirm the diagnosis is a better approach than a therapeutic trial with levothyroxine. Free T4 measured by equilibrium dialysis has a higher specificity than tT4; however, it is laborious and expensive (11). A recombinant human TSH (rhTSH) stimulation test (TSHst) or a thyroid scintigraphy are currently considered to be the gold standards for confirming hypothyroidism; however, both are expensive for the owners and not readily available for primary care practitioners (4, 17–20). Recently, prediction models aimed to assist the diagnostic process have been formulated both in human and veterinary medicine, using regression or machine learning methods (21–23). A prediction model could improve the blood testing predictive values, thus supporting or discouraging the therapeutic choice, and suggesting when it is or is not appropriate to carry out more demanding tests, such as gold standard testing. This approach could improve the overall diagnostic accuracy and decrease the rate of misdiagnosis. To be clinically useful in primary care practice, an easy-to-use tool must be developed from a model.
The aim of this study was to develop an easy-to-use prediction tool to assist in clinical decision-making when evaluating dogs with suspected hypothyroidism. It was hypothesized that a machine learning-based model built on clinical signs and clinicopathological data would help in defining the likelihood of hypothyroidism in dogs in which it was suspected.
Materials and methods
Study design
The digital patient management system (Fenice, Zaksoft Software Technology) of the Veterinary Teaching Hospital of the University of Bologna was searched for dogs in which a joint measurement of basal tT4 and endogenous TSH serum concentrations, or an rhTSH stimulation test were available in the period between January 2006 and September 2020. The medical records of the hospital were searched for information regarding signalment (i.e., sex, age, and breed), body weight, body condition score (BCS), presence/absence of clinical signs (i.e., asthenia, lethargy, polyuria/polydipsia, changes in appetite, obesity, alopecia, dermatological signs, and neurological alterations), complete blood count (i.e., hematocrit value [HCT], hemoglobin concentration [Hgb], red blood cell [RBC] count, mean corpuscular volume [MCV] and mean corpuscular hemoglobin concentration [MCHC]), serum chemistry (i.e., cholesterol, triglycerides, and creatinine concentrations, alanine aminotransferase [ALT], and aspartate aminotransferase [AST] activities), urinalysis (specific gravity and urine protein to creatinine ratio), and thyroid function evaluation (i.e., basal and post-rhTSH stimulation serum tT4 concentrations and TSH concentrations). The category of dermatological signs included all dermatological abnormalities other than alopecia, such as dry/poor quality coat, skin hyperpigmentation, pyoderma, seborrhea and recurrent otitis. When available, a diagnosis and follow-up of the patients were obtained. Dogs were excluded from the study based on the following criteria: (a) complete lack of anamnestic, clinical or clinicopathological data, (b) thyroid function testing carried out for research purposes or to monitor hypothyroid dogs receiving treatment, (c) the presence of congenital or suspect secondary hypothyroidism, and (d) treatment with levothyroxine in the month before the hormonal tests. Congenital hypothyroidism was defined based on age of presentation and typical clinical signs, while secondary hypothyroidism was suspected in dogs with acquired hypothyroidism and concurrent pituitary macrotumor. Dogs with concurrent diseases and dogs that were receiving medications known to affect serum tT4 and TSH concentration were included, providing their diagnosis was based on rhTSH stimulation results. Cases in which a clear confirmation or exclusion of hypothyroidism could not be obtained were also excluded from the study. Dogs were classified as euthyroid based on a serum post-rhTSH stimulation tT4 > 1.7 μg/dL or if they showed (a) a normal tT4 and TSH serum concentration and (b) lack of clinical or clinicopathological abnormalities consistent with hypothyroidism, or a diagnosis other than hypothyroidism was reached, or if clinical signs improved without treatment with levothyroxine. Dogs were classified as hypothyroid based on a serum post-TSH stimulation tT4 < 1.3 μg/dL or if they showed an increased TSH serum concentration associated with a decreased basal tT4 serum concentration, with clinical or clinicopathological signs consistent with hypothyroidism (4). Two authors (AC, FF) reviewed all the case records.
Clinicopathological analysis
All the analysis was performed in the Clinical Pathology Laboratory of our Institution. Fasting blood samples collected for biochemistry analysis and for hormonal assays were processed on a routine basis, according to quality standard procedures. Urine samples were prepared using low-speed centrifugation. The urine supernatant was separated and specific gravity was analyzed using a spectrometer. The serum and urine chemistry analysis were carried out using two different analyzers, which replaced each other during the inclusion period: Olympus AU400 (Two Corporate Center Drive, Melville, New York, United States) from 2006 to 2016, and Beckman Coulter-Olympus AU 480 (Brea, California, United States) from 2016 to 2020; both analyzers used the same methods for serum and urine chemistry analysis. Blood samples for the complete blood count were collected in K3 ethylene diamine tetra-acetic acid (EDTA) tubes and analyzed using two different analyzers, which replaced each other during the inclusion period: Abbott Cell-Dyn 3,500 (Abbott Laboratories, Green Oaks, Illinois, United States) from 2006 to 2010, and ADVIA 2120 (Siemens Healthcare Diagnostics, Tarrytown, New York, United States) from 2010 to 2020. Hormonal analysis on the serum was carried out using two different analyzers, which replaced each other during the inclusion period: Immulite One (Medical Systems SpA, Genova, Italy) from 2006 to 2017, and Immulite 2000 (Siemens Healthineers, Flanders, New Jersey, United States) from 2017 to 2020. The serum tT4 concentration was determined using commercially available chemiluminescent enzyme immunometric assays (Immulite Canine Total T4, Diagnostic Products Corporation; Immulite 2000 Canine Total T4, Diagnostic Products Corporation, Los Angeles, California, United States) validated for use in dogs (reference range 1 to 3.98 μg/dL [12.8 to 51.2 nmol/L]) (24, 25). The serum TSH concentrations were measured using chemiluminescent enzyme immunometric assays (Immulite Canine TSH, Diagnostic Products Corporation; Immulite 2000 Canine TSH, Diagnostic Products Corporation, Los Angeles, California, United States) validated for use in dogs (upper limit of the reference range 0.38 ng/mL) (26). Based on product data sheets, both the method comparison between Immulite 2000 canine total T4 and Immulite One canine total T4, and the method comparison between Immulite 2000 canine TSH and Immulite One canine TSH showed strong correlations (R = 0.991 and R = 0.988, respectively) (PIL2KCT-5 Immulite 2000 Canine Total T4 [June 27, 2005]; PIL2KKT-15 Immulite 2000 Canine TSH [March 6, 2017]). The results of the biochemistry analysis, the hormonal assays, and the complete blood count were compared with reference intervals calculated internally in the Clinical Pathology Laboratory, according to previously published guidelines (27). The TSH stimulation test was carried out using a dose of 75 μg/dog of rhTSH (Thyrogen, Genzyme Corporation, Suffolk, UK) (4). Blood samples were collected immediately before and 6 h after the intravenous administration of rhTSH.
Statistical analysis
Machine learning models were built to describe how the variability of the response variable (i.e., hypothyroidism, no hypothyroidism) was generated by the relationship with the explanatory variables (i.e., clinical and clinicopathological parameters) using an algorithm approach.
The dataset obtained from the medical records was analyzed to identify and exclude variables which had >5% of missing data. Then, feature selection was carried out to remove irrelevant and redundant information, and define a subset of variables which would provide good prediction results. The feature selection was conducted in two phases. In the first phase, the significantly related variables were identified by means of an analysis of mutual dependence and the variables which had a significant relationship with the objective variable were identified using a filter method, such as mutual information (MI); a permutation test was applied to the MI to verify its significance. Based on these results, 4 different models, which included different combinations of variables, were created. Specifically, the models differed based on the inclusion of quantitative or qualitative variables and the inclusion or exclusion of thyroid hormones concentrations. Model 1 and 2 were built using quantitative variables, without and with thyroid hormones concentrations, respectively; Model 3 and 4 were built using qualitative variables, without and with thyroid hormones concentrations, respectively. In the qualitative models, the categories were defined as normal if the value was within the reference interval; increased if the value was above the RI; decreased if the value was below the RI; markedly increased (cholesterol) if the value was 2 times higher than the upper reference limit and markedly decreased (tT4) if the value was below the lower limit of detection of the assay (6 mmol/L). In the second phase, a Wrapper method, in which the influence of the learning algorithm was considered, was applied to these 4 models. Six different learning algorithms, namely Classification Tree (CT), Random Forest (RF), Gradient Boosting Machine (GBM), Support Vector Machine (SVM), Logistic Regression (LR) and naive Bayes (NB), were applied to cover a broad spectrum of existing tools. The evaluation of the predictive performance was performed on the original dataset in a 10-fold cross-validation setup (i.e., the initial sample was randomly divided into 10 sub-samples and, from time to time, each of these 10 sub-samples was used as a validation set as compared to the other 9 sub-samples, used as training sets), repeated 10 times, to prevent the results from being affected by the selection of the training and test samples. For each setup, the set of classification rules which best predicted the data of the training set was created; then, the classification rules were applied to the test set and, finally, within the test set, the predicted values were compared with the real observed values, and the performance of the model was evaluated by calculating the normalized Matthews correlation coefficient (nMCC). The nMCC ranges from 0 to 1, with 0 indicating perfect misclassification, 0.5 indicating random classification, and 1 indicating perfect classification (28). The predictive performance estimated by this process was defined as the apparent performance because it was calculated on the same dataset used to develop the prediction model. For each prediction model, the optimal model was defined as the one with the highest apparent performance, defined as the highest nMCC. The sensitivity, specificity, positive predictive values (PPV), negative predictive values (NPV), accuracy, and area under the receiver operating characteristic curve (AUROC) of the optimal models were reported. The statistical analysis was carried out using Cran-R software, and the statistical significance was set at p < 0.05.
User interface building-up
A graphic user interface was created to facilitate the implementation of the prediction models in clinical practice. In order to make the application more user-friendly, an attempt was made to standardize the interface by identifying what was the common combination of variables of the different optimal models which presented the best predictive capacity, and then verifying whether these modified models were significantly different from the optimal ones. The nMCC of the optimal and the modified models were compared using t-test.
Results
Dataset
The medical records of 619 dogs were recovered from the digital patient management system; of these, 304 (49%) cases presented at least one exclusion criterion. Three hundred and fifteen (51%) cases were ultimately included in the study. Eighty-two (26%) dogs were classified as hypothyroid and 233 (74%) were classified as euthyroid.
In the hypothyroid group, there were 40 (49%) males, of which 6 were neutered, and 42 (51%) females, of which 28 were spayed. The median age was 8.0 years (range, 1.8–15.6), and the median body weight was 30.0 kg (range, 7.3–69.0). The most represented breeds were mixed breeds (n = 26; 32%), Dobermann Pinschers (n = 7; 9%), and English Setters (n = 5; 6%). Serum TSH concentration was within reference interval in 20 (24%) hypothyroid dogs, all of which were diagnosed based on rhTSHst. Overall, in 39 (48%) cases an rhTSHst was carried out to confirm the clinical suspicion of hypothyroidism.
In the euthyroid group, there were 123 (53%) males, 31 of which were neutered, and 110 (47%) females, of which 66 were spayed. The median age was 9.2 years (range, 1.3–17.3), and the median body weight was 25.5 kg (range, 1.4–84.0). Most represented breeds were mixed breeds (n = 59; 25%), Labrador Retrievers (n = 22; 9%), Dobermann Pinschers (n = 18; 8%), and Golden Retrievers (n = 13; 6%). In 67 (29%) cases, a TSHst was required to exclude hypothyroidism due to a clinical (n = 45; 67%), clinicopathological (n = 12; 18%), or clinical and clinicopathological (n = 10; 15%) suspicion of the disease.
For both groups, the results of quantitative and qualitative variables are listed in Tables 1, 2, respectively. Concurrent diseases were overall reported in 198/315 (63%) dogs, 21 (26%) hypothyroid and 177 (76%) euthyroid dogs.
Table 1. Descriptive statistics for quantitative variables in hypothyroid dogs (n = 82) and euthyroid dogs (n = 233).
Table 2. Descriptive statistics for qualitative variables in hypothyroid dogs (n = 82) and euthyroid dogs (n = 233).
Statistical model
Analysis of the missing data is reported in Supplementary Table S1. The result of the analysis of the correlation performed between the continuous variables is represented in Supplementary Table S2. Based on this, the Hgb and RBC counts were removed because of the overlapping with HCT. The results of the feature selection identified 11 variables: breed, HCT, serum concentrations of cholesterol, creatinine, serum tT4 and TSH concentrations, and the presence of alopecia, dermatopathy, lethargy/depression, asthenia, and/or obesity (Table 3). These variables were grouped into 4 different models, as previously described (Table 4). The results of the performance evaluation for the 4 different models are reported in Figure 1. Based on the values of the nMCC, the best-performing algorithms were NB for models 1 and 2, LR for model 3 and RF for model 4 (Table 5). The sensitivity, specificity, PPV, NPV, accuracy, and AUROC of the optimal models were described in Table 6.
Figure 1. Box and whisker plots showing the apparent performances of the 4 models. For each model, the predictive performance of the 6 machine-learning algorithms (i.e., Classification Tree [CT], Random Forest [RF], Gradient Boosting Machine [GBM], Support Vector Machine [SVM], Logistic Regression [LR] and naïve Bayers [NB] were sorted using a normalized Matthews Correlation Coefficient (nMCC). The boxes represent the interquartile range from the 25th to the 75th percentile. The horizontal bar in each box represents the median value. The whiskers represent the interquartile range from the 2.5th to the 97.5th percentile, with the outliers represented as dots. The dotted lines represent the limits of the reference interval.
Table 5. Comparison between the explanatory variables and the predictive performances of the best performing machine learning models (optimal models) and of the machine learning models selected for implementation in the user-friendly prediction tool (modified models).
Table 6. Main indicators of the predictive performance of the best performing machine learning models (optimal models) and machine learning models selected for implementation in the user-friendly prediction tool (modified models).
User interface
The modified models, which included a set of variables common among different models, did not show a significant loss in predictive ability as compared to the optimal models (Table 5). The interface was designed using a sequential approach. First, the veterinarian had to select the preferred model; then, the clinical and clinicopathological variables required could be added to the model using input boxes. At that point, the interface displayed a percent chance that the patient had hypothyroidism (Figures 2, 3).
Figure 2. Graphical user interface for Model 3. The user selects the appropriate category for each parameter from the drop-down menus and clicks the ‘Calculate Probability’ button. The algorithm displays the probability of the dog being hypothyroid as a percentage.
Figure 3. Graphical user interface for model 4. The user enters medical record numbers using the menu and clicks the ‘Calculate Probability’ button. The algorithm displays the probability of the dog being hypothyroid as a percentage.
Discussion
The Authors demonstrated that the prediction models described in the present study, based on clinical and clinicopathological parameters obtained from a single-institution and built using the algorithm modeling approach, have good accuracy in determining whether a dog is hypothyroid or not. A graphic user interface was also created which translated the results of these models into a predicted likelihood of the disease, thus enabling their application in clinical practice, as long as external validation supports the diagnostic accuracy described herein.
Estimating the probability of a disease is inherently a multivariable-based process; every clinician naturally integrates signalment, clinical signs and test results to assess the probability of disease which can very rarely be defined on the basis of a single predictor. A multivariable diagnostic prediction model is a mathematical equation which relates multiple predictors for a particular individual with the probability of the presence (diagnosis) of a particular outcome (29, 30). There are no validated diagnostic prediction models for hypothyroidism in dogs; however, they could prove to be extremely useful, considering that hypothyroidism can be difficult to diagnose especially when concurrent diseases are present, serum TSH concentrations are not increased, and gold-standard testing (i.e., rhTSH stimulation test, thyroid scintigraphy) is not available or possible. Several studies have applied machine learning to create diagnostic prediction models in dogs and cats for different diseases, including kidney diseases in cats, canine leishmaniasis and endocrine diseases, such as hypoadrenocorticism and Cushing’s Syndrome in dogs (21, 22, 31–33). All the models presented in these studies were, for the most part, built on clinical data and the results of easy-to-perform laboratory testing results with good predictive performance. The predictors assessed in this study were identified a priori using the current knowledge of the disease, based on the existing scientific literature. All the predictors were demographic factors, clinical signs or clinicopathological data commonly assessed in dogs with suspected hypothyroidism which were easy and affordable to obtain in a clinical setting. The 9 predictors included in the final models were: alopecia, dermatological signs, lethargy/depression, asthenia, HCT, serum cholesterol concentration, and serum creatinine concentration plus serum tT4 and TSH concentrations. These results were consistent with the existing literature. Alopecia, other dermatological signs and lethargy are clinical signs commonly reported in hypothyroid dogs. None of them is specific; however, taken together, they are fairly indicative since the majority of hypothyroid dogs show at least one of them (2, 4, 34). Hypercholesterolemia and mild to moderate normocytic normochromic non-regenerative anemia are the most common clinicopathological abnormalities in hypothyroid dogs. Hypercholesterolemia, caused by an altered lipid metabolism, is reported in 70–80% of cases, and it has been suggested that the larger the increase in cholesterol the more likely hypothyroidism rather than non-thyroidal illness (2, 4, 34). Anemia is less common but still described in approximately 30–40% of cases, presumably as a consequence of the decreased production of erythropoietin (2, 4, 34). Azotemia results from decreased glomerular filtration rate due to lack of thyroid hormones and has been variably reported in hypothyroid dogs; the majority of studies have reported azotemia in approximately 10–15% of cases; however, it reached up to 33% in a recent study (1, 3, 4).
All four models presented in the present study showed good to excellent apparent performance in predicting the presence of hypothyroidism. The term apparent performance defined the predictive ability of a model quantified on the same data from which the model was built; it could result in an optimistic estimate of performance, due to overfitting. For this reason, an internal validation is necessary to adjust the model for overfitting which was carried out using a 10-fold cross-validation.
Based on the inclusion criteria applied (i.e., measurement of thyroid hormones), the original dataset, included for the most part, but not only, dogs in which hypothyroidism was at least deemed possible by the Authors or by the referring veterinarian. If dogs in which hypothyroidism was not suspected (i.e., sick dogs without clinical signs or clinicopathological abnormalities suggestive of hypothyroidism) were included, it is likely that the predictive performance of the present models would have been different. However, to optimize the predictive performance and its clinical utility, the model must be developed using a sample representative of the target population, namely dogs in which hypothyroidism is suspected; in fact, it is senseless and misleading to apply the prediction models if hypothyroidism is not even suspected at first.
As expected, the performance of the prediction models markedly improved when tT4 and TSH concentrations were included in the model. Even if models 1 and 3 were the least performing since they did not include tT4 and TSH, they were important since they could be applied to dogs when only clinical findings and a minimum clinicopathological database were available, and before thyroid hormone evaluation, thus helping in the assessment of the pre-test probability. Moreover, the models were built including both continuous variables (models 1 and 2) and mixed continuous and categorical variables (models 3 and 4). The continuous variables allowed for a fine-tuned assessment of the variable of interest; however, the reference interval for any variable usually depends on the laboratory which carries out the analysis. Thus, models 1 and 2 which incorporated categorical variables could be better suited for application in clinical practice, considering that the reference intervals for specific variables (e.g., serum cholesterol concentration) could vary between different centers, sometimes greatly, and that the reference intervals reported in the present study are not universally valid.
The prediction models presented in this study were neither meant to be used as substitutes for clinical reasoning nor as a gold-standard for diagnosis, but rather to support clinicians in their decision-making, helping them to define the likelihood of disease in an individual dog. In addition, the present prediction tool could be helpful in everyday practice in many different ways. At the time of the first consultation, the tool could give the clinicians a pre-test probability helping them to decide whether additional testing (e.g., serum tT4 and TSH measurement) was warranted or not. Even more, if thyroid testing is subsequently carried out, the tool would quantify the actual predictive value of the results obtained. A classic example is the approach to tT4 assessment in a dog with low pre-test-probability (e.g., no alopecia, no dermatological signs, no lethargy, no anemia and no hypercholesterolemia); in fact, mildly/moderately decreased tT4 has a low positive predictive value in this setting and, at the least, additional testing is required before starting treatment. On the contrary, if the dog had a high pre-test probability, moderately decreased tT4 could be enough to merit treatment with levothyroxine. Thus, obtaining a quantitative value of the predicted likelihood of hypothyroidism could aid in avoiding a misdiagnosis and help communication with the owners regarding the reasons behind the decision to discard hypothyroidism, start treatment, or suggest gold-standard testing.
Of note, any prediction model, even if proven accurate, should be trusted and routinely applied in clinical practice only if built on solid data and, as much as possible, free from bias. For this reason, some limitations should be assessed. The main limitation of the present study was the number of cases included; in fact, larger sample sizes would lead to the development of more robust prediction models. A general rule of thumb for the sample size required is to include 10 events (i.e., hypothyroid dogs) for each candidate predictor variable considered in the model (35). Based on this, the ideal sample for this study should have included approximately 200 hypothyroid dogs. However, models including less than 10 events per parameter should not be systemically disregarded and could still be clinically useful, especially if based on strong predictors, as is the case of the present models (36). The Authors carried out a case by case review of the medical records, and stringent inclusion criteria were applied to the original dataset. This approach was appropriate in order to minimize the risk of misclassification bias; however, some hypothyroid dogs were likely excluded from the study due to the inability to confirm the diagnosis. Furthermore, specific clinical signs not clearly stated as ‘present’ or ‘not present’ in the medical records were recorded as ‘not present’ in the present dataset, based on the assumption that clinicians usually do not record absent clinical signs. It is possible that some clinical signs remained unnoticed by the owners and were misclassified by the attending clinician; however, this was likely the same in both groups. The use of tT4 and TSH concentrations as part of the inclusion criteria may have overestimated the performance of model 2 and 4, but hormone concentrations were not used as sole criteria for classifying the dogs and, most importantly, the algorithm was designed to specifically consider the degree of serum tT4 decrease and serum TSH elevation. Finally, before applying these results to very different populations, an external validation should be performed. The Authors aim to carry out an external validation of their models, which can be done either by using independent data collected by the same authors but sampled from a later period of time or by using data collected by different investigators in different clinical settings (i.e., primary vs. secondary care) or different centers. This approach strongly supports the clinical applicability of prediction models and allows for updating and adjustments in the case that poor performance is detected (29, 30).
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
Ethical approval was not required for the studies involving animals in accordance with the local legislation and institutional requirements because it is not required by the local animal ethicals for studies using medical records and retrospectively collected data. No procedure was performed on animals for this study. Written informed consent was not obtained from the owners for the participation of their animals in this study because it is not required by the local animal ethicals for studies using medical records and retrospectively collected data. No procedure was performed on animals for this study.
Author contributions
AC: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. FL: Conceptualization, Data curation, Formal analysis, Investigation, Resources, Visualization, Writing – original draft, Writing – review & editing. FA: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. ID: Formal analysis, Methodology, Software, Supervision, Writing – original draft. EF: Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing. FF: Conceptualization, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
Part of this study was presented as an oral abstract at the online ECVIM Congress 2021.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1292988/full#supplementary-material
References
1. Mooney, CT . Canine hypothyroidism: a review of aetiology and diagnosis. N Z Vet J. (2011) 59:105–14. doi: 10.1080/00480169.2011.563729
2. Dixon, RM , Reid, SW , and Mooney, CT . Epidemiological, clinical, haematological and biochemical characteristics of canine hypothyroidism. Vet Rec. (1999) 145:481–7. doi: 10.1136/vr.145.17.481
3. di Paola, A , Carotenuto, G , Dondi, F , Corsini, A , Corradini, S , and Fracassi, F . Symmetric dimethylarginine concentrations in dogs with hypothyroidism before and after treatement with levothyroxine. J Small Anim Pract. (2021) 62:89–96. doi: 10.1111/jsap.13212
4. Corsini, A , Faroni, E , Lunetta, F , and Fracassi, F . Recombinant human thyrotropin stimulation test in 114 dogs with suspected hypothyroidism: a cross-sectional study. J Small Anim Pract. (2021) 62:257–64. doi: 10.1111/jsap.13290
5. Gaskill, CL , Burton, SA , Gelens, HCJ , Ihle, SL , Miller, JB , Shaw, DH, et al. Changes in serum thyroxine and thyroid-stimulating hormone concentrations in epileptic dogs receiving phenobarbital for one year. J Vet Pharmacol Ther. (2000) 23:243–9. doi: 10.1046/j.1365-2885.2000.00278.x
6. Kantrowitz, LB , Peterson, ME , Melián, C , and Nichols, R . Serum total thyroxine, total triiodothyronine, free thyroxine, and thyrotropin concentrations in dogs with nonthyroidal disease. J Am Vet Med Assoc. (2001) 219:765–9. doi: 10.2460/javma.2001.219.765
7. Mooney, CT , Shiel, RE , and Dixon, RM . Thyroid hormone abnormalities and outcome in dogs with non-thyroidal illness. J Small Anim Pract. (2008) 49:11–6. doi: 10.1111/j.1748-5827.2007.00418.x
8. Hume, KR , Rizzo, VL , Cawley, JR , and Balkman, CE . Effects of Toceranib phosphate on the hypothalamic-pituitary-thyroid Axis in tumor-bearing dogs. J Vet Intern Med. (2018) 32:377–83. doi: 10.1111/jvim.14882
9. Harper, A , Blackwood, L , and Mason, S . Investigation of thyroid function in dogs treated with the tyrosine kinase inhibitor toceranib. Vet Comp Oncol. (2020) 18:433–7. doi: 10.1111/vco.12538
10. NISHII, N , OKADA, R , MATSUBA, M , TAKASHIMA, S , KOBATAKE, Y , and KITAGAWA, H . Risk factors for low plasma thyroxine and high plasma thyroid-stimulating hormone concentrations in dogs with non-thyroidal diseases. J Vet Med Sci. (2019) 81:1097–103. doi: 10.1292/jvms.19-0169
11. Dixon, RM , and Mooney, CT . Evaluation of serum free thyroxine and thyrotropin concentrations in the diagnosis of canine hypothyroidism. J Small Anim Pract. (1999) 40:72–8. doi: 10.1111/j.1748-5827.1999.tb03040.x
12. Boretti, FS , and Reusch, CE . Endogenous TSH in the diagnosis of hypothyroidism in dogs. Schweiz Arch Tierheilkd. (2004) 146:183–8. doi: 10.1024/0036-7281.146.4.183
13. Shiel, RE , Brennan, SF , Omodo-Eluk, AJ , and Mooney, CT . Thyroid hormone concentrations in young, healthy, pretraining greyhounds. Vet Rec. (2007) 161:616–9. doi: 10.1136/vr.161.18.616
14. Seavers, A , Snow, D , Mason, K , and Malik, R . Evaluation of the thyroid status of basenji dogs in Australia. Aust Vet J. (2008) 86:429–34. doi: 10.1111/j.1751-0813.2008.00357.x
15. Lavoué, R , Geffré, A , Braun, JP , Peeters, D , and Trumel, C . Breed-specific biochemical reference intervals for the adult Dogue de Bordeaux. Vet Clin Pathol. (2013) 42:346–59. doi: 10.1111/vcp.12067
16. Uhríková, I , Lačňáková, A , Tandlerová, K , Kuchařová, V , Řeháková, K , Jánová, E, et al. Haematological and biochemical variations among eight sighthound breeds. Aust Vet J. (2013) 91:452–9. doi: 10.1111/avj.12117
17. Boretti, FS , Sieber-Ruckstuhl, NS , Favrot, C , Lutz, H , Hofmann-Lehmann, R , and Reusch, CE . Evaluation of recombinant human thyroid-stimulating hormone to test thyroid function in dogs suspected of having hypothyroidism. Am J Vet Res. (2006) 67:2012–6. doi: 10.2460/ajvr.67.12.2012
18. Diaz Espineira, MM , Mol, JA , Peeters, ME , Pollak, YWEA , Iversen, L , van Dijk, JE, et al. Assessment of thyroid function in dogs with low plasma thyroxine concentration. J Vet Intern Med. (2007) 21:25–32. doi: 10.1892/0891-6640(2007)21[25:AOTFID]2.0.CO;2
19. Boretti, FS , Sieber-Ruckstuhl, NS , Wenger-Riggenbach, B , Gerber, B , Lutz, H , Hofmann-Lehmann, R, et al. Comparison of 2 doses of recombinant human thyrotropin for thyroid function testing in healthy and suspected hypothyroid dogs. J Vet Intern Med. (2009) 23:856–61. doi: 10.1111/j.1939-1676.2009.0336.x
20. Shiel, RE , Pinilla, M , McAllister, H , and Mooney, CT . Assessment of the value of quantitative thyroid scintigraphy for determination of thyroid function in dogs. J Small Anim Pract. (2012) 53:278–85. doi: 10.1111/j.1748-5827.2011.01205.x
21. Schofield, I , Brodbelt, DC , Niessen, SJM , Church, DB , Geddes, RF , Kennedy, N, et al. Development and internal validation of a prediction tool to aid the diagnosis of Cushing’s syndrome in dogs attending primary-care practice. J Vet Intern Med. (2020) 34:2306–18. doi: 10.1111/jvim.15851
22. Reagan, KL , Reagan, BA , and Gilor, C . Machine learning algorithm as a diagnostic tool for hypoadrenocorticism in dogs. Domest Anim Endocrinol. (2020) 72:106396. doi: 10.1016/j.domaniend.2019.106396
23. Saberi-Karimian, M , Khorasanchi, Z , Ghazizadeh, H , Tayefi, M , Saffar, S , Ferns, GA, et al. Potential value and impact of data mining and machine learning in clinical diagnostics. Crit Rev Clin Lab Sci. (2021) 58:275–96. doi: 10.1080/10408363.2020.1857681
24. Bruner, JM , Scott-Moncrieff, JCR , and Williams, DA . Effect of time of sample collection on serum thyroid-stimulating hormone concentrations in euthyroid and hypothyroid dogs. J Am Vet Med Assoc. (1998) 212:1572–5.
25. Wolff, EDS , Bilbrough, G , Moore, G , Guptill, L , and Scott-Moncrieff, JC . Comparison of 2 assays for measuring serum total thyroxine concentration in dogs and cats. J Vet Intern Med. (2020) 34:607–15. doi: 10.1111/jvim.15703
26. Kooistra, HS , Diaz-Espineira, M , Mol, JA , van den Brom, WE , and Rijnberk, A . Secretion pattern of thyroid-stimulating hormone in dogs during euthyroidism and hypothyroidism. Domest Anim Endocrinol. (2000) 18:19–29. doi: 10.1016/S0739-7240(99)00060-0
27. Friedrichs, KR , Harr, KE , Freeman, KP , Szladovits, B , Walton, RM , Barnhart, KF, et al. ASVCP reference interval guidelines: determination of de novo reference intervals in veterinary species and other related topics. Vet Clin Pathol. (2012) 41:441–53. doi: 10.1111/vcp.12006
28. Chicco, D , Starovoitov, V , and Jurman, G . The benefits of the Matthews correlation coefficient (MCC) over the diagnostic odds ratio (DOR) in binary classification assessment. IEEE Access. (2021) 9:47112–24. doi: 10.1109/ACCESS.2021.3068614
29. Moons, KGM , Altman, DG , Reitsma, JB , Ioannidis, JPA , Macaskill, P , Steyerberg, EW, et al. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): explanation and elaboration. Ann Intern Med. (2015) 162:W1–W73. doi: 10.7326/M14-0698
30. Collins, GS , Reitsma, JB , Altman, DG , and Moons, KGM . Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ. (2015) 350:g7594. doi: 10.1136/bmj.g7594
31. Bradley, R , Tagkopoulos, I , Kim, M , Kokkinos, Y , Panagiotakos, T , Kennedy, J, et al. Predicting early risk of chronic kidney disease in cats using routine clinical laboratory tests and machine learning. J Vet Intern Med. (2019) 33:2644–56. doi: 10.1111/jvim.15623
32. Renard, J , Faucher, MR , Combes, A , Concordet, D , and Reynolds, BS . Machine-learning algorithm as a prognostic tool in non-obstructive acute-on-chronic kidney disease in the cat. J Feline Med Surg. (2021) 23:1140–8. doi: 10.1177/1098612X211001273
33. Ferreira, TS , Santana, EEC , Jacob Junior, AFL , Silva Junior, PF , Bastos, LS , Silva, ALA, et al. Diagnostic classification of cases of canine Leishmaniasis using machine learning. Sensors. (2022) 22:1–13. doi: 10.3390/s22093128
34. Panciera, DL . Hypothyroidism in dogs: 66 cases (1987-1992). J Am Vet Med Assoc. (1994) 204:761–7.
35. Riley, RD , Ensor, J , Snell, KIE , Harrell, FE Jr, Martin, GP , Reitsma, JB, et al. Calculating the sample size required for developing a clinical prediction model. BMJ. (2020) 368:m441. doi: 10.1136/bmj.m441
36. Vittinghoff, E , and McCulloch, CE . Relaxing the rule of ten events per variable in logistic and cox regression. Am J Epidemiol. (2007) 165:710–8. doi: 10.1093/aje/kwk052
Glossary
Keywords: diagnosis, canine, thyroid, thyroxine, endocrinology, logistic regression, machine learning, artificial intelligence
Citation: Corsini A, Lunetta F, Alboni F, Drudi I, Faroni E and Fracassi F (2023) Development and internal validation of diagnostic prediction models using machine-learning algorithms in dogs with hypothyroidism. Front. Vet. Sci. 10:1292988. doi: 10.3389/fvets.2023.1292988
Edited by:
Carmel T. Mooney, University College Dublin, IrelandReviewed by:
Barbara Contiero, University of Padova, ItalyTommaso Banzato, University of Padua, Italy
Copyright © 2023 Corsini, Lunetta, Alboni, Drudi, Faroni and Fracassi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Andrea Corsini, YW5kcmVhLmNvcnNpbmlAdW5pcHIuaXQ=