Metabolomic changes in cats with renal disease and calcium oxalate uroliths

Introduction There is a significant incidence of cats with renal disease (RD) and calcium oxalate (CaOx) kidney uroliths in domesticated cats. Foods which aid in the management of these diseases may be enhanced through understanding the underlying metabolomic changes. Objective Assess the metabolomic profile with a view to identifying metabolomic targets which could aid in the management of renal disease and CaOx uroliths. Method This is a retrospective investigation of 42 cats: 19 healthy kidney controls, 11 with RD, and 12 that formed CaOx nephroliths. Cats were evaluated as adults (2 through 7 years) and at the end of life for plasma metabolomics, body composition, and markers of renal dysfunction. Kidney sections were assessed by Pizzolato stain at the end of life for detection of CaOx crystals. CaOx stone presence was also assessed by analysis of stones removed from the kidney at the end of life. Results There were 791 metabolites identified with 91 having significant (p < 0.05, q < 0.1) changes between groups. Many changes in metabolite concentrations could be explained by the loss of renal function being most acute in the cats with RD while the cats with CaOx stones were intermediate between control and RD (e.g., urea, creatinine, pseudouridine, dimethylarginines). However, the concentrations of some metabolites differentiated RD from CaOx stone forming cats. These were either increased in the RD cats (e.g., cystathionine, dodecanedioate, 3-(3-amino-3-carboxypropyl) uridine, 5-methyl-2′-deoxycytidine) or comparatively increased in the CaOx stone forming cats (phenylpyruvate, 4-hydroxyphenylpyruvate, alpha-ketobutyrate, retinal). Conclusions The metabolomic changes show specific metabolites which respond generally to both renal diseases while the metabolomic profile still differentiates cats with RD and cats with CaOx uroliths.


Introduction
During the mid-1980's there was a reported increase in CaOx stone formation of domestic cats (Osborne et al., 2009). In humans, CaOx stones are also the most prevalent stone type (Lieske, et al., 2014). The occurrence of CaOx uroliths in cats has compartively declined so that in 2003 the frequency of CaOx uroliths (44.3%) was similar to struvite (44.9%) (Osborne et al., 2009). Calcium oxalate uroliths are unlike struvite uroliths in that although there are a number of possible management options (DeFarges et al., 2013) they are not readily dissolvable (Lulich et al, 2013). There is a known connection between CaOx uroliths and renal disease (RD) (O'Kell et al, 2017). It is likely that the increased loss of renal function seen in some cats with CaOx stones is through damage from urine blockage. This may be because of an increased incidence of infection or other associated damage. Cats with CaOx uroliths were reported to have increased symmetric dimethylarginine (SDMA), creatinine, and reduced urine specific gravity (uSG) as well as having a reduced lifespan as compared to all other cats (Hall et al., 2017b).
The relationship of metabolomic changes in RD has been recently reviewed (Wang et al., 2019). In this review it was reported that RD resulted in plasma concentration changes resulting from changes in the metabolism of: energy, nucleic acids, amino acids, lipids, and carbohydrates. These data were mostly from changes in the plasma of humans with RD.
Of specific interest to this study were the reported increase in arginine and dimethyl arginines, citrulline, allantoin, pseudouridine, and kynurenic acid. A previous study reported the differences in metabolites of healthy cats and cats with RD (Hall et al, 2020). The metabolites were either increased (urea, citrulline, dimethyl arginines, N-acetylarginine, glycine, S-adenosylhomocysteine) or were decreased (malate, γ-tocopherol/βtocopherol) in RD cats as compared to healthy controls. The metabolomic changes associated with CaOx stones are less reported. However, a study reporting CaOx stones in rats (Gao et al., 2016) reported 13 metabolite changes in the serum. Of specific interest are the increased concentrations of: 4-Hydroxy-l-proline, creatine, valerylglycine, and cholic acid. The metabolites that declined with stone formation were: succinic acid, l-carnitine, propionylcarnitine, 2-Methylbutyroylcarnitine, and taurocholic acid.
Some of the likely causes of RD include infectious disease (Pedersen, 2009), and renal lymphosarcoma (DiBartola et al., 1987), as well as nephritis and other inflammatory diseases (Weiss et al, 1996). Early detection of acute renal damage can facilitate intervention and reduce disease severity (Grauer, 2005). Regarding CaOx stone formation supersaturation of calcium oxalate is required, however, healthy normal cats may have urine which is supersaturated for CaOx (Hezel et al, 2007). A review of CaOx stone formation (Bartges, 2016) states that known risk factors include: increased urinary excretion of calcium or oxalate, aciduria, a reduction in urinary inhibitors of stone formation, and reduced urine volume.
Cats with RD benefit from foods with controlled phosphorus and protein (Elliott et al, 2000;Ross et al, 2006). There is also a benefit for cats with early stage renal decline (SDMA > 14 μg/dl) which was not severe enough to be RD to consume foods enhanced with dietary n-3 fatty acids and antioxidants (Hall et al., 2016). An increased water intake and foods with moderate dietary magnesium, phosphorus, or calcium concentrations have been shown to aid in the management of CaOx urolith formation (Lekcharoensuk et al., 2001). Increased dietary fatty acids which are an aid in the management of RD have also been shown to reduce the risk of feline CaOx (Hall et al, 2017a).
The hypothesis evaluated in this study is that each specific disease state will result in metabolomic changes in the plasma associated with that disease. Plasma was used to reduce the metabolite contamination of hemolysis. This metabolomic milieu as well as the markers of renal function, body composition, and circulating concentrations of calcium and phosphorus were evaluated to not only differentiate cats with these diseases but suggest targets for solutions that aid in their management.

Materials and methods
This study was approved by the Institutional Animal Care and Use Committee, Hill's Pet Nutrition, Inc., Topeka, KS. There were 42 cats used in this study (Table 1). All cats in the Hill's Pet Nutrition, Inc. colony at the end of their life were available for inclusion in this study. Cats were assigned to the CaOx stone formers if during their life they had confirmed CaOx stones. Cats were assigned to the RD group if they were diagnosed with renal disease. Cats were chosen for the healthy kidney controls by attempting to match age at death of the two disease groups and having died with no perceived stones or loss of renal function.
During their annual physical examination, plasma was analyzed for the reported analytes except metabolomics and SDMA. Metabolomics and SDMA concentrations were evaluated from banked samples (stored at − 80 °C). Body composition was evaluated at the end of life and during adult life by dual energy X-ray methods previously described (Toll et al., 1994). For the variables where adult cat response was calculated, the mean value of each cat averaged from analysis after 1 year of age and before 8 years of age. All cats were allowed exercise in indoor runs and had access to toys. Cats were housed with access to natural light and enclosed porches with lighting that varied with seasonal changes. In accordance to accepted veterinary practice all pets were cared for and supported during their life to optimize quality of life, even in the presence of CaOx stones or renal disease. Samples were obtained in a postmortem examination to discover the cause of death or the extent of disease. Euthanasia was performed in order to minimize the pets' pain and suffering when it was concluded that the pets' quality of life was degraded to a point where it was no longer humane to do otherwise.

Metabolomics
Samples archived for metabolomic analysis were collected during annual exams or at the end of life. The Healthy Kidney group (n = 19) had 100 samples analyzed, the Renal Disease group (n = 11) had 52 samples analyzed and the CaOx group (n = 12) had 60 samples anlayzed. A commercial laboratory (Metabolon, Morrisville, NC, USA) performed the analysis as described previously (Hall et al., 2018). In brief, samples were extracted with methanol under vigorous shaking for 2 min to precipitate protein and dissociate small molecules bound to protein or trapped in the precipitated protein matrix, followed by centrifugation to recover chemically diverse metabolites. The resulting extract was divided into five fractions: two for analysis by two separate reverse phase methods using positive ion mode electrospray ionization (ESI), one for analysis by RP/UPLC-MS/MS using negative ion mode ESI, one for analysis by HILIC/UPLC-MS/MS using negative ion mode ESI, and one reserved for backup. Proprietary software was used to match ions to an in-house library of standards for metabolite identification and for metabolite quantitation by peak area integration.

Pizzolato staining
Staining for detecting CaOx stones in renal tissues was done using the modified method described (Pizzolato, 1964). In short, both kidneys were removed from the abdominal cavity during the necropsy, excess tissue was removed, and assessed for visual normality. Any visual abnormalities are noted at this time. Formalin-fixed samples (10% formalin) collected from necropsies were transferred to the histopathology lab at Kansas State University for further pathological processing and analysis. The paraffin-embedded blocks (FFPE) are returned to in house for further staining. The FFPE were sectioned at five microns and transferred onto slides for the Pizolatto staining procedure.

Statistical analysis
Statistical analysis was completed using SAS 9.4 or JMP (SAS Institute, Cary, NC, USA). To test the hypothesis that classification (RD, CaOX, or healthy) was a significant component of the observed variation, statistical analysis of the metabolomic data was performed on natural log transformed data and age was used as a covariant. As it is common for metabolomic data to not have normally distributed residuals in the analysis of variance, log transformation was done in order to improve the models by reducing uniquely high or low values to produce normal or near-normal residuals.
Although other transformations can be useful (such as square root or cube root transformations) the log transformation is most common and has been defined as standard practice for metabolomic analysis (Xia et al., 2013). All analytes with p < 0.05 and q < 0.1 were included in Table 2. For the nonmetabolomic data, adult values were calculated by taking the mean of the values before diagnosis from 2 to 7 years of age.
In both metabolomic and non-metabolomic data analysis cat identification was used as a random effect when there was multiple measures for the same response of a given cat. Age was used as a continuous variable in the metabolomic data as it was not binned between adult and end of life analysis. Classification was always used as a fixed effect. The posthoc comparisons for the metabolomic data used Welch's two-sample t-test, while in the non-metabolomic data the PDIFF option was used in PROC Mixed of SAS. Gender distribution was evaluated by the PROC Frequency procedure of SAS. A p-value < 0.05 was used as a significant cutoff for all anlaysis. In order to control the false discovery rate in the metabolomic analysis a q value of < 0.1 was used. This gives an accurate measure of the level of false positives and strikes a balance between accepting a difference that in reality only exists because of random variation and missing a difference that is real but does not rise to the pre-established level of significance. This approach to false discovery rate has been shown to be valuable (Storey, 2002) and recently reviewed in comparison to other methods of controling the false discovery rate (Cheng & Pounds, 2007). A correlation analysis was completed to evaluate the effect of gender and lean body mass on indicators of renal disease (SDMA and creatinine).

Classification of cats
All cats were placed in their treatment classification after the end of life. Cats were assigned to the CaOx stone forming group by the presence of stones in the kidney that were analyzed as CaOx stones. Cats were assigned to the renal group by evaluation of the kidney, and concentrations of SDMA, creatinine, and uSG. A cat was assigned to the renal disease group if SDMA was over 18 µg/dl, creatinine over 1.6 mg/dl, uSG under 1.030 or had structural signs of renal disease. Cats with confirmed CaOx stones were not made available for the RD group. Healthy kidney cats had no signs of stones or RD throughout life or at death. Calcium Oxalate stone formation was also confirmed through histology, as shown in Fig. 1.

Plasma metabolomics
There were a significant number of metabolites that increased concentrations in RD and CaOx cats when compared to cats with healthy kidneys. Most of these responded similarly (with an elevated circulating concentrations in the RD and CaOx cats) with the RD cats usually having a more pronounced response. The ten analytes with decreased concentrations in the RD cats (as compared to the healthy cats) were never changed in a significant way in the CaOx cats. This then resulted in eight of the analytes that differed between the CaOx and RD cats being comparatively elevated in the RD cats and seven being comparatively elevated in the CaOxo cats.
a Values are ratios of means between groups (i.e., RD/Healthy is concentration in cats with renal disease/concentration in cats with healthy kidneys; CaOx/Healthy is concentration in cats with CaOx stones/concentration in cats with healthy kidneys; RD/CaOx is concentration in cats with renal disease/concentration in cats with CaOx stones. Values in bold are above 1 (p < 0.05). Values in italics are lower than 1 (p < 0.05)

Body composition
At the end of life, cats with CaOx stones had less lean than healthy controls, while cats with RD were not different from either group. During adult life, cats with CaOx stones had a lower fat perentage and mass than cats with RD, healthy controls were not different from either group (Table 3).

Specific blood and urine response variables
At the end of life, the cats with RD had a lower uSG, as well as higher circulating concentrations of SDMA, creatinine, and urea (Fig. 2). Urine specific gravity was reduced and circulating phophosphorus concentrations were increased when cats with RD were compared to the healthy controls (Table 3). Cats with CaOx stones were not different from the other groups in circulating concentrations, but were different from both groups in uSG (Table 4). There was an effect of gender on creatinine at the end of life (with males having a higher value than females). However, this is likely due to the positive relationship between lean body mass and circulating creatinine as there was no effect of gender when lean was used as a covariant. Gender or lean body mass had no effect on SDMA concentrations. Neutrophils, lymphocytes, and the neutrophil to lymphocyte ratio (NLR) were not different in the cat's adult years (NLR is an important measure of systemic inflammation). At the end of life, cats with RD had a greater percentage of neutrophils than cats that developed CaOx stones, with healthy controls being not

Discussion
There was a general increase in concentrations of those metabolites dependent on renal clearance in both the RD and CaOx cats (e.g., urea, creatinine, dimethylarginines, pseudouridine). This was shown in that there was an elevation in both cats forming CaOx and RD cats as compared to those with healthy kidneys while there was not an observed difference in creatinine, methylated arginines, pseudouridine or urea between RD and CaOx forming cats. These changes associated with renal function were expected while the other metabolite changes allow for a deeper understanding of these diseases and possible changes in the nutrition to aid in their management. These data show the expected increase in creatinine as a result of the renal function decline. In a recent review (Kovarikova, 2018), the relationship between creatinine and renal function was elucidated. Increasing lean body mass, and therefore increased creatine, was associated with increased circulating creatinine concentrations (Hall et al., 2014). This explains the gender effect as males, which had higher lean body mass, would be expected to have higher creatinine concentrations. In this study creatinine was highest in the cats with RD, intermediate in the cats with CaOx stones and lowest in cats with healthy kidneys. These changes were not the result of changing lean body mass. These data support the conclusion that some obstruction resulted in renal function decline in the cats with CaOx stones. However, this loss of renal function as a group did not equal the function loss of the cats with RD.
The cats with RD also had increased circulating SDMA concentrations as compared to healthy controls, while cats with CaOx stones were intermediate between these groups. The circulating concentrations of SDMA was elevated in cats with CaOx stones as compared to healthy controls (Hall et al., 2017b), but it to our knowledge has not been shown in comparison to cats with RD. Because changes in lean body mass did not account for changes in circulating concentrations of SDMA, there is an improved understanding of renal function through the use of circulating SDMA concentration as compared to circulating creatinine. The lower threshold in renal decline associated with increasing circulating SDMA may contribute to the elevation of circulating SDMA being observed before the elevation of creatinine in cats (Hall et al, 2014).
Urinary specific gravity was different between all three groups at the time of death. This again supports the conclusion that although there was renal decline associated with CaOx stones, although it was not as severe in the RD group. Increasing solute load, specifically urinary calcium, is associated with an increased risk of calcium oxalate stones (Bartges, 2016). However, the pre-disease adult concentrations of solutes, as shown by uSG was not different between groups suggesting that this was not the stimulus for stone formation in these cats.
Phosphorus concentrations were elevated with RD as compared to healthy controls. Elevated circulating phosphorus is a significant risk factor for continued renal decline and death and a direct result of renal decline (Laflamme et al, 2020). Increased circulating calcium concentrations have also been associated with renal decline (Tang et al, 2021) in phosphate restricted foods. The lack of an increase in circulating calcium concentrations (there was a numeric decline) in this study in the cats with RD may be because the cats when diagnosed were transitioned to a different food designed to aid in the management of renal disease than that fed to the cats where circulating calcium concentration increased. There were metabolomic changes resulting in decreased concentrations in the RD cats while remaining unchanged in the CaOx stone forming cats. The metabolite N-acetylasparagine was reduced in RD cats while being unchanged in the healthy controls or cats with CaOx stones. This elevation has previously been reported in humans (Boccard, et al., 2021) and may be the result of a change in the toxification method associated with RD (Sekula, et al., 2016). The relative reduction in phenylpyruvate in RD cats is similar to what has been reported in humans (Kopple, 2007) and may be the result of impaired phenylalanine hydroxylase activity resulting in a reduced phenylalanine to tyrosine conversion (Tourian et al., 1969). There was a reduction in the circulating concentrations of alpha-ketobutyrate with increased concentrations of cystathionine in the RD cats and not in the CaOx cats. This is likely the direct result of homocysteine being irreversibly degraded through the transsulfuration pathway (condensing homocysteine and serine to form cystathionine which is split into alpha-ketobutyrate and cysteine). This could have a significant effect on health as alpha ketobutyrate has been reported to have protective effects in vitro in neuronal cells when stressed with hydrogen peroxide (Jacewicz et al., 2014). These changes support the conclusion that homocysteine degradation is disrupted in cats with RD. Because high circulating concentrations of homocysteine are known to be risk factors in RD in humans (Ostrakhovitch & Tabibzadeh, 2015) and that betaine has been shown to reduce homocysteine and increase methionine in cats (Jewell et al., 2021;Verbrugghe & Bakovic, 2013) these results suggest that dietary betaine could be a benefit to aid in the management of RD in cats. These data then provide understanding for the positive effect of a food containing increased dietary betaine on lean body mass in cats with RD (Ephraim-Gebreselassie et al., 2018) as well as the reduction in oxidized glutathione when dietary betaine was combined with short chain fructooligosaccharide (Hall et al., 2020).
The metabolite 1-palmitoyl-GPI, was elevated in the cats with CaOx stones both as compared to cats with RD and cats with healthy kidneys. The elevation of 1-palmitoyl-GPI has been associated with an increased inflammatory condition in humans where it was correlated to the pro-inflammatory cytokine IL-8 (Gao et al., 2021) and was increased in the inflammation associated with periodontitis in mice (Ilievski et al, 2016). Increased concentrations of inflammatory glycosylphosphatidylinositols were observed in the cats eating the more inflammatory foods as compared to initial values or those cats eating a less inflammatory (high fish oil) food . It is an intriguing basis for further research to better understand if dietary changes which change the microbiota in cats with CaOx stones would aid in the management of this disease. The increased inflammatory markers in cats with CaOx stones may explain the benefit of a food with a reduced inflammation fatty acid profile in cats with CaOx stones (Hall et al, 2017a(Hall et al, , 2017b. The higher NLR in the cats with RD is likely a function of increased inflammation at end of life. This ratio has been associated with increased inflammation in both people (McMillan, 2009) and cats (Neumann, 2021). This change at the end of life shows the inflammatory whole-body response happening in RD. Although CaOx stones are known to have an inflammatory component (O'kell et al, 2017), this was not reflected in a significant elevation of the NLR. This could result from a more localized inflammation in the cats with CaOx stones compared to those with RD and account for the changing inflammatory metabolomics. This lymphopenia observed in the RD cats was also observed by Kralova- Kovarikova et al. (2016) who found that the reduction was negatively correlated to circulating concentrations of creatinine and urea and resulted in immunosuppression. The reduced immune response was shown by a reduced proliferation in unstimulated cells or cells stimulated with phytohaemaglutinin, concanavalin A, or pokeweed mitogen. The reduced immune response associated with the reduction in NLR may be a significant target for enhancing the length and quality of life in cats with RD.
At the end of life there was a reduction in body lean in the cats with CaOx stones which may have resulted from reduced food intake during the last phase of life. The reduced body fat (and fat percentage) during adult life with body lean being unchanged (similar to healthy controls) suggests a partitioning of available energy away from body fat. It has been reported that lean is preferentially maintained over body fat in cases of reduced food intake (Floerchinger et al, 2015). Although that research was done with food designed to aid in the management of weight loss, it may be that cats with limited calories during life preferentially support body lean, and this is what is observed in the CaOx forming cats body composition.

Conclusion
Cats with both RD and CaOx stones have metabolomic signs of renal disease and specific metabolomic changes which both give insight into successful management and targets for influence to further aid in the management of these diseases. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.