1
|
Pocar P, Scarpa P, Berrini A, Cagnardi P, Rizzi R, Borromeo V. Diagnostic potential of simplified methods for measuring glomerular filtration rate to detect chronic kidney disease in dogs. J Vet Intern Med 2019; 33:2105-2116. [PMID: 31444986 PMCID: PMC6766551 DOI: 10.1111/jvim.15573] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Accepted: 07/11/2019] [Indexed: 12/04/2022] Open
Abstract
Background Glomerular filtration rate (GFR) is the most sensitive indicator of initial renal function decline during chronic kidney disease (CKD), but conventional protocols for measuring GFR are labor‐intensive and stressful for the dog. Objectives To assess the diagnostic potential for detecting CKD with simplified GFR protocols based on iohexol plasma clearance. Animals Seventeen CKD‐positive and 23 CKD‐negative dogs of different breeds and sex. Methods Prospective nonrandomized study. Plasma iohexol was measured 5, 15, 60, 90, and 180 minutes after injection. Glomerular filtration rate was calculated using 5 samples (GFR5) or simplified protocols based on 1, 2, or 3 samples. The GFR5 and simplified GFR were compared by Bland‐Altmann and concordance correlation coefficient (CCC) analysis, and diagnostic accuracy for CKD by receiver operating characteristic curves. A gray zone for each protocol was bounded by the fourth quartile of the CKD‐positive population (lower cutoff) and the first quartile of the CKD‐negative population (upper cutoff). Results All simplified protocols gave reliable GFR measurements, comparable to reference GFR5 (CCC >0.92). Simplified protocols which included the 180‐minutes sampling granted the best GFR measure (CCC: 0.98), with strong diagnostic potential for CKD (area under the receiver operating characteristic curve ± SE: 0.98 ± 0.01). A double cutoff including a zone of CKD uncertainty guaranteed reliable diagnosis outside the gray area and identified borderline dogs inside it. Conclusions The simplified GFR protocols offer an accurate, hands‐on tool for CKD diagnosis in dogs. The gray zone might help decision‐making in the management of early kidney dysfunction.
Collapse
Affiliation(s)
- Paola Pocar
- Department of Veterinary Medicine, Università degli Studi di Milano, Milan, Italy
| | - Paola Scarpa
- Department of Veterinary Medicine, Università degli Studi di Milano, Milan, Italy
| | - Anna Berrini
- Department of Veterinary Medicine, Università degli Studi di Milano, Milan, Italy
| | - Petra Cagnardi
- Department of Health, Animal Science and Food Safety, Università degli Studi di Milano, Milan, Italy
| | - Rita Rizzi
- Department of Veterinary Medicine, Università degli Studi di Milano, Milan, Italy
| | - Vitaliano Borromeo
- Department of Veterinary Medicine, Università degli Studi di Milano, Milan, Italy.,Department of Health, Animal Science and Food Safety, Università degli Studi di Milano, Milan, Italy
| |
Collapse
|
2
|
Paepe D, Lefebvre HP, Concordet D, van Hoek I, Croubels S, Daminet S. Simplified methods for estimating glomerular filtration rate in cats and for detection of cats with low or borderline glomerular filtration rate. J Feline Med Surg 2014; 17:889-900. [PMID: 25518848 DOI: 10.1177/1098612x14561106] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
OBJECTIVES Diagnosis of early feline chronic kidney disease (CKD) is challenging. Glomerular filtration rate (GFR) is the best overall indicator of kidney function, but multisample plasma clearance methods to determine GFR are labour intensive, time consuming and stressful for feline patients. This study aimed to develop simplified methods to detect decreased GFR in cats. METHODS Data from a nine-sample combined plasma exogenous creatinine-iohexol clearance test of 73 cats were used. Limited sampling strategies were developed by comparing all sampling time combinations with the complete nine sampling times set and selecting the best sampling time combinations based on maximum relative error. By regression analysis, the ability of routine blood (serum creatinine, serum urea) and urine (urine specific gravity, urinary protein:creatinine ratio) variables to predict GFR or identify cats with low or borderline GFR was examined. Cut-off clearance marker concentrations to predict low or borderline GFR was determined at three time points after marker injection. All procedures were analysed for three clearance markers (exo-iohexol, creatinine, endo-iohexol). RESULTS For reliable estimation of GFR, at least three blood samples for clinical purposes and five blood samples for research purposes are required. Regression formulae based on routine variables did not reliably predict GFR, but accurately identified cats with low (sensitivity 96.5-98.2%; specificity 60-91.3%) or borderline (sensitivity 91.1-96%; specificity 76.5-81.8%) GFR. Clearance marker concentrations exceeding given marker cut-off concentrations also identified cats with low or borderline GFR with high sensitivities and specificities. CONCLUSIONS AND RELEVANCE These simplified methods will facilitate the detection of early kidney dysfunction in cats. Early diagnosis allows timely therapeutic intervention, and future studies must reveal whether this improves the long-term outcome of cats with CKD.
Collapse
Affiliation(s)
- Dominique Paepe
- Department of Medicine and Clinical Biology of Small Animals, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
| | - Hervé P Lefebvre
- University of Toulouse, INP, National Veterinary School of Toulouse, Clinical Research Unit, Toulouse, France
| | - Didier Concordet
- INRA, UMR 1331, Toxalim, National Veterinary School of Toulouse, 23 Chemin des Capelles, Toulouse, France
| | - Ingrid van Hoek
- Department of Medicine and Clinical Biology of Small Animals, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium Current address: Royal Canin SAS, BP 4 - 650, Avenue de la Petite Camargue, 30470, Aimargues, France
| | - Siska Croubels
- Department of Pharmacology, Toxicology and Biochemistry, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
| | - Sylvie Daminet
- Department of Medicine and Clinical Biology of Small Animals, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
| |
Collapse
|
3
|
Reynolds BS, Massal MR, Nguyen P, Grégoire LL, Périgaud AE, Concordet D, Biourge V, Lefebvre HP. Plasma exogenous creatinine clearance in clinically healthy cats: comparison with urinary exogenous creatinine clearance, tentative reference intervals and indexation to bodyweight. Vet J 2014; 202:157-65. [PMID: 25193408 DOI: 10.1016/j.tvjl.2014.07.022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2013] [Revised: 04/30/2014] [Accepted: 07/26/2014] [Indexed: 10/25/2022]
Abstract
Glomerular filtration rate (GFR) is considered to be the best indicator of overall kidney function. The major objectives of this study were to compare plasma exogenous creatinine clearance (PECC) with a reference method, to establish reference intervals (RIs) for PECC and to assess the effects of indexation of GFR to bodyweight (BW) in cats. PECC was compared with urinary clearance of exogenous creatinine (UECC) in six clinically healthy domestic shorthair cats (experiment 1). Tentative RIs were determined according to current guidelines and the effects of indexation to BW and of covariables on GFR were assessed in 43 clinically healthy cats of various breeds (experiment 2). PECC was 15% higher than UECC (P <0.01), but the two estimates were strongly correlated (r(2)=0.97, P = 0.001). RIs for PECC were 6.4-21.3 mL/min or 1.2-4.9 mL/min/kg. The absolute (i.e. non-indexed) GFR value was not dependent on BW. Thus, indexation of GFR to BW in cats would not standardize the GFR value, but could introduce bias in clinical interpretation. Significant effects of breed, plasma protein concentration and plasma albumin concentration on GFR were demonstrated. Plasma concentrations of urea and creatinine, when assessed separately, were also weakly correlated with GFR in healthy cats. These combined findings contribute to a better understanding of renal function assessment in cats.
Collapse
Affiliation(s)
- B S Reynolds
- Unité de Recherche Clinique, Université de Toulouse, INP, Ecole Nationale Vétérinaire de Toulouse, F-31076 cedex 03 Toulouse, France.
| | - M R Massal
- Unité de Recherche Clinique, Université de Toulouse, INP, Ecole Nationale Vétérinaire de Toulouse, F-31076 cedex 03 Toulouse, France
| | - P Nguyen
- Unité de Nutrition et d'Endocrinologie, Oniris, F-44307, Nantes Cedex, France
| | - L L Grégoire
- Unité de Recherche Clinique, Université de Toulouse, INP, Ecole Nationale Vétérinaire de Toulouse, F-31076 cedex 03 Toulouse, France
| | - A E Périgaud
- Unité de Recherche Clinique, Université de Toulouse, INP, Ecole Nationale Vétérinaire de Toulouse, F-31076 cedex 03 Toulouse, France
| | - D Concordet
- UMR 1331 Toxalim, INRA, Université de Toulouse, INP, Ecole Nationale Vétérinaire de Toulouse, F-31076 cedex 03 Toulouse, France
| | - V Biourge
- Royal Canin SAS, Centre de Recherches, F-30470, Aimargues, France
| | - H P Lefebvre
- Unité de Recherche Clinique, Université de Toulouse, INP, Ecole Nationale Vétérinaire de Toulouse, F-31076 cedex 03 Toulouse, France
| |
Collapse
|
4
|
Abstract
Practical relevance: Routinely used markers of renal function in clinical practice include urea and creatinine. However, these are insensitive markers, particularly in the early stages of kidney disease. Measurement of glomerular filtration rate (GFR) is regarded as the most sensitive index of functioning renal mass. It may be useful for feline patients in varying clinical scenarios; for example, where a more accurate measurement of renal function may aid diagnosis, to enable response to therapeutic interventions to be more closely monitored, or to evaluate renal function prior to the use of nephrotoxic or renally cleared drugs. Clinical challenges: Traditional methods of measuring GFR, such as renal clearance or multisample plasma clearance techniques, are generally impractical for clinical use. Limited sampling and single sample plasma clearance methods using the filtration marker iohexol have been validated in cats. These have the advantages of reduced stress to cats associated with repeated sampling and reduced costs of analysis, and therefore offer greater clinical utility. Attempts to develop an estimated GFR (eGFR) formula similar to that used in human patients have been made in cats, although currently an accurate and reliable formula is not available. Audience: This review presents the basis for the theoretical understanding and practical measurement of GFR for any veterinary practitioner wishing to obtain a more accurate and sensitive assessment of renal function than routinely used markers provide. Evidence base: The review draws evidence from peer-reviewed publications, the author’s PhD thesis and also clinical experience.
Collapse
Affiliation(s)
- Natalie Finch
- School of Veterinary Sciences, University of Bristol, Langford, Bristol BS40 7UD, UK
| |
Collapse
|
5
|
Finch NC, Heiene R, Bird NJ, Elliott J, Peters AM. Evaluation of non-polynomial equations for one-compartment correction of slope-intercept GFR: theoretical prediction and experimental measurement. Scandinavian Journal of Clinical and Laboratory Investigation 2014; 74:611-9. [PMID: 25017421 DOI: 10.3109/00365513.2014.928941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
BACKGROUND Polynomial equations for one-compartment correction of slope-intercept glomerular filtration rate (GFR) will underestimate values at high clearance rates. Non-polynomial correction equations that are independent of patient size and renal function would be advantageous and may have cross-species use. MATERIALS AND METHODS The study explored the theoretical basis of firstly the Jodal and Brochner-Mortensen one-compartment correction equation, replacing plasma volume with extracellular fluid volume, and secondly an equation described by Peters. One-compartment correction factors (a which is related to plasma volume and v which is related to extracellular fluid volume) which avoided the need for scaling to body size were developed. Both factors were determined from the biexponential clearance curve of the markers iohexol and (51)Cr-EDTA in humans and iohexol in cats and dogs. Relationships between a and v and filtration function and body size were then determined using data from humans, cats and dogs to assess their validity and compare this with theoretical predictions. RESULTS In all species, v was higher than a, as theoretically predicted. Both were significantly higher in humans than cats and dogs, ruling out cross-species use. Significant relationships were present between v and measures of filtration function in humans, but were weak with respect to a. Neither a nor v showed significant relationships with filtration function in animals or with body size in any species. CONCLUSIONS a and v (which are factors independent of body size) can be used interchangeably for correcting slope-intercept clearance. However values of both for humans are higher compared to cats and dogs. Therefore a single cross-species factor cannot be used.
Collapse
|
6
|
Paepe D, Daminet S. Feline CKD: Diagnosis, staging and screening - what is recommended? J Feline Med Surg 2014; 15 Suppl 1:15-27. [PMID: 23999183 DOI: 10.1177/1098612x13495235] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
PRACTICAL RELEVANCE Feline chronic kidney disease (CKD) is frequently encountered by veterinarians. Timely diagnosis and staging may facilitate the initiation of adequate therapy and improve the prognosis for patients. CLINICAL CHALLENGES Feline CKD is diagnosed based on the presence of compatible clinical signs and renal azotaemia, which implies that urinalysis (particularly urine specific gravity) is mandatory to confirm the diagnosis. Although the diagnosis of advanced feline CKD and associated complications is usually straightforward, based on complete blood and urine examination, all routine blood and urine tests have their limitations in detecting early CKD. Therefore, diagnosing early or non-azotaemic CKD is much more challenging. Although determination of glomerular filtration rate (GFR) would be ideal to identify early kidney dysfunction, practical limitations hamper its routine use in clinical practice. PATIENT GROUP CKD is typically a disease of aged cats, but may affect cats of all ages. Conclusive breed and sex predispositions for feline CKD are not reported. AUDIENCE This review is directed at practising veterinarians and provides an overview of the required diagnostic tests, the classification system established by the International Renal Interest Society, and the importance of and possible techniques for early detection of CKD. EVIDENCE BASE Staging of cats with CKD is essential as it directs management and provides a prognostic guide. Given that diagnosis at early disease stages is associated with more prolonged survival times, simple, inexpensive and accurate methods for early CKD diagnosis are needed. Techniques currently under investigation include limited sampling strategies to estimate GFR, clearance marker cut-off concentrations to identify cats with low GFR, new indirect GFR markers and urinary biomarkers.
Collapse
Affiliation(s)
- Dominique Paepe
- Department of Small Animal Medicine and Clinical Biology, Faculty of Veterinary Medicine, Ghent University, Salisburylaan 133, B-9820 Merelbeke, Belgium.
| | | |
Collapse
|
7
|
Finch N, Heiene R, Elliott J, Syme H, Peters A. A Single Sample Method for Estimating Glomerular Filtration Rate in Cats. J Vet Intern Med 2013; 27:782-90. [DOI: 10.1111/jvim.12089] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2011] [Revised: 02/14/2013] [Accepted: 03/11/2013] [Indexed: 11/29/2022] Open
Affiliation(s)
| | - R. Heiene
- Norwegian School of Veterinary Science; Oslo Norway
| | | | | | | |
Collapse
|
8
|
Debruyn K, Vandermeulen E, Saunders JH, Dobbeleir AA, Ham HR, Peremans K. Effect of background region of interest and time-interval selection on glomerular filtration ratio estimation by percentage dose uptake of (99m)Tc-DTPA in comparison with (51)Cr-EDTA clearance in healthy cats. J Feline Med Surg 2013; 15:698-705. [PMID: 23349527 DOI: 10.1177/1098612x13475618] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Evaluation of glomerular function is a useful part of the diagnostic approach in animals suspected of having renal disease. Time-interval and background region of interest (bg ROI) selection are determining factors when calculating the glomerular filtration ratio (GFR) based on percentage uptake of (99m)technetium-labelled diethylene triamine penta-acetic acid ((99m)Tc-DTPA). Therefore, three different time intervals (60-120 s, 120-180 s, 60-180 s) and three different bg ROIs (C-shape, caudolateral, cranial + caudal) were investigated. In addition, global GFRs based on percentage dose uptake of (99m)Tc-DTPA for the different time-intervals and bg ROIs were compared with the global GFR based on (51)chromium-ethylene diaminic tetra-acetic acid ((51)Cr-EDTA) plasma clearance in nine healthy European domestic shorthair cats. Paired Student's t-tests and linear regression analysis were used to analyse the data. Different time intervals seemed to cause significant variation (P <0.01) in absolute GFR values, regardless of the choice of bg ROI. Significant differences (P <0.01) between bg ROIs were only observed in the 120-180s time interval between the C-shape and cranial + caudal bg ROI, and between the caudolateral and cranial + caudal bg ROI. The caudolateral bg ROI in the 60-180 s time interval showed the highest correlation coefficient (r = 0.882) between (99m)Tc-DTPA and (51)Cr-EDTA, although a significant difference (P <0.05) was present between both techniques.
Collapse
Affiliation(s)
- Katrien Debruyn
- Department of Veterinary Medical Imaging and Small Animal Orthopaedics, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium.
| | | | | | | | | | | |
Collapse
|
9
|
Von Hendy-Willson VE, Pressler BM. An overview of glomerular filtration rate testing in dogs and cats. Vet J 2011; 188:156-65. [PMID: 20541957 PMCID: PMC2980797 DOI: 10.1016/j.tvjl.2010.05.006] [Citation(s) in RCA: 79] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2009] [Revised: 04/21/2010] [Accepted: 05/07/2010] [Indexed: 12/25/2022]
Abstract
Determination of glomerular filtration rate (GFR) is a valuable, yet underused, diagnostic tool for evaluating renal function in dogs and cats. This article first reviews the hormonal and hemodynamic factors which contribute to GFR, followed by a description of considerations when selecting a pharmacokinetic model and methods of animal-to-animal standardization. The best-characterized existing GFR markers, including creatinine, radiolabeled markers, and iohexol, are reviewed in depth, as well as alternative but lesser used techniques. A weighted means analysis of reported GFR measurements in healthy dogs and cats and a review of selected studies that have examined GFR alterations in animals with naturally occurring and experimental diseases provide the reader with preliminary guidelines on expected GFR results in these species and disease conditions.
Collapse
Affiliation(s)
| | - Barrak M. Pressler
- Department of Veterinary Clinical Sciences, Purdue University, West Lafayette, IN 47907 USA
| |
Collapse
|
10
|
Vandermeulen E, Ham HR, Dobbeleir AA, De Sadeleer C, Piepsz A, Waelbers T, Vermeire ST, Slegers G, Peremans KY. (99m)Tc-DMSA absolute and relative renal uptake in cats: procedure and normal values. J Feline Med Surg 2011; 13:387-92. [PMID: 21440474 DOI: 10.1016/j.jfms.2010.12.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2010] [Revised: 12/02/2010] [Accepted: 12/06/2010] [Indexed: 10/18/2022]
Abstract
In this study we investigated the influence of technical factors (positioning, background (BG) correction and attenuation correction) on qualitative and quantitative (absolute (AU) and relative (RU) uptake) assessment of feline kidneys with (99m)technetium labelled dimercaptosuccinic acid ((99m)Tc-DMSA). Eleven healthy adult cats were included. Influence of BG and depth correction on quantitative assessment was evaluated. Depth correction was based on the geometric mean method (using dorsal and ventral images) and the use of two standards placed over each individual kidney. Visual evaluation showed superiority of dorsal and ventral over lateral positioning due to increased separation of the kidneys permitting region of interest (ROI) placement without overlap. No apparent influence of BG correction was found for RU. However, AU was systematically overestimated without BG correction. Depth correction did not seem to affect RU in most cases, however, in some cats the differences were not negligible. The values for AU without depth correction were lower compared to depth corrected values.
Collapse
Affiliation(s)
- Eva Vandermeulen
- Department of Veterinary Medical Imaging and Small Animal Orthopaedics, Faculty of Veterinary Medicine, Ghent University, Salisburylaan 133, B-9820 Merelbeke, Belgium.
| | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Determination of optimal sampling times for a two blood sample clearance method using 51Cr-EDTA in cats. J Feline Med Surg 2010; 12:577-83. [DOI: 10.1016/j.jfms.2010.01.011] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/21/2010] [Indexed: 11/22/2022]
Abstract
Estimation of the glomerular filtration rate (GFR) is a useful tool in the evaluation of kidney function in feline medicine. GFR can be determined by measuring the rate of tracer disappearance from the blood, and although these measurements are generally performed by multi-sampling techniques, simplified methods are more convenient in clinical practice. The optimal times for a simplified sampling strategy with two blood samples (2BS) for GFR measurement in cats using plasma 51chromium ethylene diamine tetra-acetic acid (51Cr-EDTA) clearance were investigated. After intravenous administration of 51Cr-EDTA, seven blood samples were obtained in 46 cats (19 euthyroid and 27 hyperthyroid cats, none with previously diagnosed chronic kidney disease (CKD)). The plasma clearance was then calculated from the seven point blood kinetics (7BS) and used for comparison to define the optimal sampling strategy by correlating different pairs of time points to the reference method. Mean GFR estimation for the reference method was 3.7±2.5 ml/min/kg (mean±standard deviation (SD)). Several pairs of sampling times were highly correlated with this reference method ( r2≥0.980), with the best results when the first sample was taken 30 min after tracer injection and the second sample between 198 and 222 min after injection; or with the first sample at 36 min and the second at 234 or 240 min ( r2 for both combinations=0.984). Because of the similarity of GFR values obtained with the 2BS method in comparison to the values obtained with the 7BS reference method, the simplified method may offer an alternative for GFR estimation. Although a wide range of GFR values was found in the included group of cats, the applicability should be confirmed in cats suspected of renal disease and with confirmed CKD. Furthermore, although no indications of age-related effect were found in this study, a possible influence of age should be included in future studies.
Collapse
|
12
|
MIYAGAWA Y, TAKEMURA N, HIROSE H. Evaluation of a Single Sampling Method for Estimation of Plasma Iohexol Clearance in Dogs and Cats with Various Kidney Functions. J Vet Med Sci 2010; 72:271-8. [DOI: 10.1292/jvms.09-0294] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Affiliation(s)
- Yuichi MIYAGAWA
- Laboratory of Veterinary Integrative Medicine, School of Veterinary Medicine, Faculty of Veterinary Science, Nippon Veterinary and Life Science University
| | - Naoyuki TAKEMURA
- Laboratory of Veterinary Integrative Medicine, School of Veterinary Medicine, Faculty of Veterinary Science, Nippon Veterinary and Life Science University
| | - Hisashi HIROSE
- Laboratory of Veterinary Internal Medicine, School of Veterinary Medicine, Faculty of Veterinary Science, Nippon Veterinary and Life Science University
| |
Collapse
|