Fisheries and Aquatic Sciences
BioMed Central
Research Article

Establishment of an analytical method for butaphosphan (BTP), a stress-attenuating agent, and its application in the preliminary pharmacokinetic evaluation of residues in olive flounder Paralichthys olivaceus

Ji-Hoon Lee1, Jun Sung Bae2, Chae Won Lee2, Chan Yeong Yang2, Ji-Sung Choi2, Sang-Hoon Choi2, Yue-Jai Kang3, Kwan Ha Park2,*
1Aquatic Disease Control DivisionNational Institute of Fisheries Science216 Gijanghaean-Ro, Gijang-Gun46083GyeongnamKorea
2Department of Aquatic Life Medicine, College of Ocean Science & TechnologyKunsan National University558 Daehak-ro54150Gunsan CityJeonbukKorea
3Department of Aquatic Life and Medical SciencesSun Moon University221 Sunmoon-ro31460Asan CityChungnamKorea

© The Author(s) 2020. Open AccessThis 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

Received: Sep 5, 2019; Accepted: Mar 2, 2020

Published Online: Apr 7, 2020



Butaphosphan (BTP) has recently been introduced into the Korean aquaculture sector as a stress-attenuating agent. In this study, a sensitive chemical analytical method was established for the detection of BTP in the olive flounder (Paralichthys olivaceus) tissues.


Utilizing a method employing liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS), detection sensitivity, specificity, and precision were satisfactorily established. Temporal changes in the BTP plasma and muscle concentrations were assessed after a single intramuscular injection of BTP (50 and 150 mg/kg) to the olive flounder maintained at 13 °C or 22 °C.


High BTP plasma levels were achieved immediately after the injection, and the drug was rapidly eliminated. Additionally, plasma BTP levels were markedly dependent on the elimination rate, which, in turn, seemed dependent on the water temperature, with the drug elimination half-life and mean residence time significantly shorter at 22 °C than 13 °C. Overall, muscle BTP levels were markedly lower than the plasma levels. Notably, muscle levels were not influenced by water temperatures. Muscle BTP concentrations were used to estimate the necessary withdrawal period for drugs used in food fish, with BTP levels maintained far below the possible hazardous limit.


In conclusion, the established LC-MS/MS method can be used for BTP residue detection with high sensitivity and reproducibility.

Keywords: Butaphosphan (BTP); Olive flounder P. olivaceus; Plasma and muscle concentrations; Residue pharmacokinetics


Butaphosphan (1-butylamino-1-methylethyl-phosphonic acid, BTP, Fig. 1) has been used in several domestic animals, usually in combination with vitamin B12, to attenuate the symptoms stress occurring during various situations. After calving in dairy cows, BTP is believed to improve metabolic ketosis and milk production during the postpartum period (Pereira et al. 2013a; Rolling et al. 2010; Tabeleão et al. 2017). Additionally, the drug is helpful in cows recovering from surgical procedures (Fürll et al. 2006; Fürll et al. 2010). In sheep stressed due to drug treatments, BTP improved appetite and metabolic adjustments (Pereira et al. 2013b). In piglets, an assessment of cortisol levels and behavioral changes has indicated reduced social stress induced by co-housing with unfamiliar individuals (Van der Staay et al. 2007). In dogs, BTP prevents abnormalities in hepatic function and digestive enzyme activity caused by an anti-inflammatory steroid (Deniz et al. 2009). Furthermore, BTP pretreatment reportedly improved endurance in rodents (Hasi et al. 2004, 2005a, 2005b).

Fig. 1. Chemical structure of 1-(butylamino)-1-methylethyl]-phosphonic acid in BTP (active component)
Download Original Figure

Although BTP has long been widely used in various livestock animals against stress-related symptoms, its pharmacological effects in fish species have remained unclear until recently. We recently recognized the usefulness of BTP as a stress-attenuating agent against conditions like hypothermia, hypoxia, or antibiotic drug injections in olive flounder: cortisol levels, reduced appetite, and immune parameters were partially prevented after injection of BTP (Seo et al. 2019). Additionally, BTP attenuates stress-induced immune-suppression and tissue damage (Mechesso et al. 2019) following over-crowding of the olive flounder. In fishes, some earlier preliminary reports have suggested possible beneficial effects, although their work may require further confirmatory studies. For example, it has been reported that maternally injected BTP improved spawning rates and relative fecundity in Nile tilapia (Alkobaby 2008). Van de Vis et al. (2004) examined whether BTP can alleviate hypoxia-stress symptoms in sole fish but failed to definitively conclude its effectiveness. Based on recent reports, BTP is now prescribed in farms for relieving stress symptoms in Korean olive flounder aquaculture (National Institute of Fisheries Sciences 2018).

Despite actual application in veterinarian fields, the mechanism of action of BTP remains unclear. It has been suggested that stress-attenuating effects may be associated with the accelerated production of adenosine-5’-triphosphate (ATP) after BTP administration to mice (Hasi et al. 2004) or fish (Seo et al. 2019).

Information concerning drug’s pharmacokinetic behavior is essential to understand the duration of action and chemical clearance time, particularly when the safety of a drug is unknown. To perform a pharmacokinetic assessment, the establishment of an analytical method is a prerequisite. Analytical procedures for BTP have been developed for piglets (Sun et al. 2016) and horses (Du 2004); however, no attempt has been made to develop a procedure suitable in fish tissues.

This study was performed to delineate the pharmacokinetic profiles of BTP in the olive flounder, a major target for the use of BTP. The analytical method utilized for the olive flounder was established using liquid chromatography coupled with mass spectrometry (LC-MS/MS), which is highly powerful for the analysis of many organic chemical drugs.

Materials and methods

Fish husbandry

Olive flounder P. olivaceus, weighing 105 ± 5 g, were obtained from a culture farm located in Taean County, Chungnam Province, Korea, and maintained in recirculating seawater culture tanks of L3.0 m × W1.0 m × H1.2 m, at a density of about 20 fish/tank. The water level was maintained at 0.8 m, dissolved oxygen 6.0 ppm, and a temperature of 22 ± 3 °C. In some studies, the temperature was set at 13 ± 3 °C to examine the temperature effect. A commercial flounder pellet feed (CJ Feed Inc., Gunsan, Jeonbuk, Korea) was supplied at a rate of approximately 3%, twice daily at 9:00 AM and 6:00 PM. Experimental protocols were approved by the Animal Care and Use Committee of the Kunsan National University, Korea (Approval number, AE 2016-04).

Test drug and analytical reagents

The test drug used was a commercially available combination product of [1-(butylamino)-1-methylethyl]-phosphonous acid 100 mg/ml and cyanocobalamin (vitamin B12) 50 μg/ml as an injectable aqueous solution (TITAN, Korea BNP, Yesan, Chungnam, Korea). In this study, this complex test drug was referred to as BTP (Fig. 1) and was diluted in 0.85% saline before the intramuscular injection at 0.5 ml/kg body weight. Most analytical reagents were obtained from Sigma (St. Louis, MO, USA) unless otherwise specified.

BTP analysis using LC-MS/MS technique

BTP analysis was performed using high-performance liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The analytical system was composed of Agilent-6420 model triple quadrupole detectors with an ionization spray as the ionization source (Agilent Technologies, Santa Clara, CA, USA). LC column used was Xbridge C18 (2.1 × 150 mm, 3.5 μm, Waters), and the mobile phase was composed of two solutions, pure water (A) and 100% methanol (B), both containing 0.1% formic acid. Elution was performed by a linear gradient change at 0.2 mL/min: flow through the LC column was started with 90% solution A for 1 min, 20% A over a 10-min period, and then a brief switch to 90% A at the end of the run to flush for the next 10 min. This mobile phase properly eluted the target molecule BTP under positive ESI conditions. The column temperature was set at 35 °C and 10 μL injection volume was used. Table 1 summarizes the mass detector conditions for all ions involved. In addition, the capillary temperature was set at 450 °C and nitrogen (N2) was used as the collision gas. Calibration curves were obtained with BTP dissolved in 0.1% formic acid at a wide range of concentrations.

Table 1. Multiple reaction monitoring (MRM) mode for detection of butaphosphan (BTP)


Retention time (min)

Molecular weight

Precursorion (m/z)

Production (m/z)

Collision energy (eV)



















Download Excel Table

Validation of LC-MS/MS analytical method

Validation of the BTP analytical technique was performed for muscle tissues as the lowest residue levels have usually been observed in the muscle after most drug administrations. The linearity test was performed with working solutions ranging from 1 to 500 μg/kg. Detection and quantification limits were obtained from signal to noise ratios of 3:1 and 10:1, respectively. The degree of precision was estimated by the intra- and inter-day repeatability with 2.5 and 25 μg/kg BTP. The recovery was estimated by comparing signal magnitudes between standard solutions and spiked samples.

BTP Injection to fish for pharmacokinetic and muscle residue studies

BTP was injected into the dorsal muscle of the olive flounder at 50 and 150 mg/kg doses. Immediately before injection, BTP was diluted to the requisite concentrations in saline. Fish samples were obtained at 0, 0.5, 1, 5, 12, 24, 48, 168 (7 days), and 335 h (14 days) after BTP injection. Blood was withdrawn under MS-222 anesthesia (20 mg/l) using heparinized syringes through the caudal vessels. Next, the fish were sacrificed for muscle dissection. In each group, 10 fish were allocated for each time point. All 10 fish were used for muscle analysis, whereas five random samples were utilized for plasma level analysis. Muscle levels were examined primarily to ascertain food safety.

Sample treatment for BTP extraction

BTP extraction from plasma and muscle tissues was performed as described by Du et al. (2004) with some modifications. Plasma was obtained by blood centrifugation for 30 min at 3000 × g under refrigeration. The plasma (500 μl) was mixed with an equal volume of methanol and centrifuged at 4000 × g to obtain the supernatant. The supernatant was mixed with phosphate-buffered saline (100 mM, pH 3.5), and loaded onto C18 solid-phase extraction cartridges (SepPak®, Waters, Milford, MA, USA), pre-activated with 2 ml methanol and distilled water. The cartridges were cleaned with 2 ml distilled water, and then BTP was eluted with 2 ml methanol. BTP in the methanol eluate was concentrated to dryness under nitrogen at 47 °C. Residues were reconstituted in 1 ml of 0.1% formic acid and filtered using 0.2 μm syringe nylon filters (Sigma-Aldrich). For muscle BTP analysis, tissues (2 g) were homogenized in 5 ml methanol and centrifuged for 10 min at 3000 × g. The supernatant was pooled and the pellet was resuspended in methanol to repeat the centrifugation process. The pooled supernatant was concentrated at 35 °C to dryness. Residues were reconstituted in 1.5 ml phosphate-buffered saline and samples were processed through identical clean-up steps as described for plasma samples.

Data analyses

Data are expressed as mean ± SD when appropriate for multiple determinations. Pharmacokinetic parameters were analyzed using the PKSolver program developed by Zhang et al. (2010). As the plasma BTP concentration did not fit any model, a non-compartmental model was applied to analyze pharmacokinetic parameters (Zhang et al. 2010). Additionally, this model examined whether a withdrawal period for the olive flounder can be determined using an analytical program (Hekman 1996) based on the observed muscle BTP levels. For this estimation, acceptable daily intake data (European Medicines Agency 2013), estimated daily fish product consumption by Koreans (63 g/day, Korea Rural Economic Institute, 2017), and exposure estimation methods (Korea Ministry of Food and Drug Administration 2017) were applied.


Analytical method quality validation

Figure 2 shows multiple residue monitoring (MRM) chromatograms developed for olive flounder tissue analysis and a standard curve of BTP spiked muscles. In the muscles, spiked BTP was clearly discernable from interfering substance signals and a wide range of BTP concentrations maintained linearity in signal magnitudes (r2 > 0.9995).

Fig. 2. Representative LC-MS/MS chromatograms and standard curve of BTP. a Blank. b Standard BTP 0.5 pg. c Extracted BTP after spiking to olive flounder muscle, 0.025 pg. d Standard curve of BTP
Download Original Figure

A summary of the analytical precision performed with two different BTP concentrations after spiking muscle tissues is shown in Table 2. The limit of quantification was 2.5 μg/kg, with the recovery rates ranging around 100% in both the intra- and inter-day reproducibility tests.

Table 2. Validation summary for analytical procedure of butaphosphan (BTP) in olive flounder muscle


Quantification limit (mg/kg)

Test concentration (mg/kg)

Intra-day recovery (mean ± SD)

Inter-day recovery (mean ± SD)

Recovery (%)

CV (%)

Recovery (%)

CV (%)




104.2 ± 8.0


105.5 ± 8.5



96.82 ± 3.0


95.3 ± 3.3


CV coefficient of variation, n = 3

Download Excel Table

Residue pharmacokinetic characteristics

Figure 3 shows changes in plasma BTP concentrations following the intramuscular administration of BTP at 50 and 150 mg/kg. In general, the 150 mg/kg BTP dose resulted in higher peak BTP concentrations than the 50 mg/kg dose. Interestingly, however, once peak concentrations were achieved following each dose, lower levels were observed in fish maintained at 22 °C during the elimination phase than those at 13 °C. Another interesting characteristic was the biphasic peak concentration pattern, in which the initial lower level peaks were accompanied by secondary higher peaks. The pharmacokinetic patterns were analyzed to obtain kinetic parameters as shown in Table 3. Peak plasma concentrations (Cmax) were roughly dose-dependent when compared between the two doses. However, the elimination half-lives were shorter at 22 °C than at 13 °C. This faster disappearance at higher temperatures was reflected as shorter elimination half-lives (t1/2E), shorter mean residence time (MRT), and smaller area-under the curve (ACU) of BTP at 22 °C compared to their matching doses at 13 °C.

Fig. 3. Plasma concentration of BTP after intramuscular injection in olive flounder. Mean values were presented without error bars for clarity (n = 5 for each time point). Levels at time 0 (non-detectable levels, n = 5) were omitted to avoid complication, and the first symbols denote levels at 0.5 h after injection
Download Original Figure
Table 3. Pharmacokinetic parameters for plasma BTP level changes after intramuscular injection in olive flounder



13 °C ± 3 °C

22 °C ± 3 °C

50 mg/kg

150 mg/kg

50 mg/kg

150 mg/kg

















































Cmax maximum concentration, Tmax time reaching maximum concentration, t1/2E half-life for elimination phase, Vd volume of distribution, Cl clearance, AUC0-∞ area-under-the-time curve to infinity, AUMC0-∞ area-under-the-moment curve to infinity, MRT0-∞ mean residence time to infinity, n = 5 for each parameter

Download Excel Table
Table 4. Pharmacokinetic parameters for muscle BTP level changes after intramuscular injection in olive flounder



13 ± 3 °C

22 ± 3 °C

50 mg/kg

150 mg/kg

50 mg/kg

150 mg/kg



































































k12 rate constant to central compartment, k21 rate constant to peripheral compartment, Vd volume of distribution, Cmax maximum concentration, Tmax time reaching maximum concentration, t1/2E half-life for elimination, Vd central volume of distribution, Vd2 peripheral volume of Cl, clearance; AUC0-∞ area-under-the-time curve, AUMC area-under-the-moment curve, MRT0-∞ mean residence time from zero to infinity, n = 10 for each parameter

Download Excel Table

The pattern of muscle BTP concentrations is shown in Fig. 4. Muscle concentrations were also dose-dependent at peak concentrations between 50 and 150 mg/kg. However, the peak levels were markedly lower in comparison with the plasma concentration after the equivalent doses and water temperatures. However, patterns between the two temperatures, no differences in the elimination phase were observed. Overall, no marked difference in concentrations was noted between the two temperatures within the same injection dose (Table 4).

Fig. 4. Muscle concentration of BTP after intramuscular injection in olive flounder. Mean values were presented without error bars for clarity (n = 10 for each time point). Levels at time 0 (non-detectable levels, n = 10) were omitted to avoid complication, and the first symbols denote levels at 0.5 h after injection
Download Original Figure

In all sets of experiments, BTP was not detected in plasma or right muscle after test drug injections, excluding the possibility of recent BTP administration.

Estimation of BTP residue standard and withdrawal time

Figure 5 presents the edible tissue (muscle) concentration over time after BTP injection and maximally acceptable tissue concentrations for food uses (European Medicines Agency 2013). As the acceptable daily intake (ADI) of BTP is 0.6 mg/kg, a person can tolerate up to 36 mg/day (60 kg body weight, 60 × 0.6 = 36). In the fish muscle, this can result in as much as 571.4 mg/kg BTP theoretically acceptable for safe consumption (daily consumption of 63 g fishery products, 1000÷63 × 36 = 571.4). Evidently, the muscle levels of BTP were far below the top broken lines (571.4 mg/kg) at any time point, in both 13 °C and 22 °C. This large separation from the limit level and actual muscle concentrations, even at peak levels, rendered withdrawal time setting unnecessary.

Fig. 5. Muscle concentration of BTP after intramuscular injection in olive flounder. Dots represent all muscle residue levels (n = 15 fish for each time point), and the two broken lines indicate upper and lower confidence limits of the highest concentrations possible, respectively. The top horizontal lines represent maximally allowable residue levels for a fish (571.4 mg/kg) calculated from acceptable daily intake information (36 mg/person)
Download Original Figure


This study was designed to develop a sensitive analytical method for the pharmacokinetic analysis of butaphosphan (BTP). Although BTP has widely been prescribed to attenuate stress symptoms in both domestic and aquaculture sectors (Deniz et al. 2009; Fürll et al. 2006; Fürll et al. 2010; Mechesso et al. 2019; Pereira et al. 2013a, 2013b; Rolling et al. 2010; Seo et al. 2019; Tabeleão et al. 2017; Van der Staay et al. 2007), analytical procedures have not well been described. In addition, understanding the pharmacokinetic characteristics of a drug is crucial to estimate the onset, duration, and termination of biological activities (Benez 1984). Furthermore, pharmacokinetic information is essential to assess the consumer safety of a drug if its residues are possibly hazardous to humans.

The analytical method for BTP determination using a singular LC-MS was developed by Du (2004), reporting the quantification limit of 0.01 mg/kg in piglet tissues. More recently, Sun et al. (2016) used an advanced LC-MS/MS method, a tandem MS system, to analyze BTP in piglet plasma samples with a quantifiable level of 0.05 mg/l. In studies by both Du et al. (2004) and Sun et al. (2016), recovery from the above animal tissues was excellent and nearly 100%. In the present study, a level as low as 0.0025 mg/kg was quantifiable, indicating over 4 times the sensitivity than methods reported in previous studies (Du 2004; Sun et al. 2016). In the current procedure and previous studies, recovery rates with the biological samples utilized were equally satisfactory, indicating that BTP recovery from animal tissues was not an obstacle in residue analysis.

We attempted to delineate the pharmacokinetic characteristics of BTP in the olive flounder by employing the analytical procedure established in the present study. One important purpose for characterizing the pharmacokinetic parameters of a drug is its relation to the temporal aspects of pharmacological activity. For example, the duration of activity can be approximately linked to the time a drug stays in the body. In turn, time is dependent on the elimination half-life (t1/2E) from plasma. The plasma t1/2E of BTP was extremely short, indicating that the duration of activity could also be quite brief if the drug works directly to exert its biological effects (Meyer and Fish, 2008). However, in the case of BTP, the duration of action is not likely to be closely linked to t1/2E as it attenuates stress symptoms acting indirectly, i.e., it stimulates adenosine-5’-triphosphate (ATP) production and the elevated ATP amount helps reduce detrimental stress responses (Hasi et al. 2004; Seo et al. 2019). By acting indirectly, the pharmacological actions can outlast tissue levels. In other words, even when BTP had already been eliminated due to incorporation into biomolecules as the phosphorus element (Rolling et al. 2010), biological actions can still occur.

An interesting phenomenon observed in the pharmacokinetic study was that the plasma BTP initially reached a high level after injection and then declined transiently. Subsequently, the level increased to reach higher peaks than the first one. The rationale underlying the occurrence of this peculiar pattern was unclear. One possible hypothesis explaining this biphasic pattern is the possibility of two discrete mobilization routes for BTP from the muscle to plasma via routes similar to multiple absorption sites model (Gabrielsson and Weiner 2000). Based on this model, some BTP molecules will immediately enter capillary vessels on arrival, and the remaining slower molecules will diffuse out to a wide tissue range, before entering into the vascular compartments. Muscle BTP concentrations did not differ between the two bath temperatures, implying that distribution from the site of injection was not markedly influenced by body temperatures.

The patterns of overall BTP plasma concentrations demonstrated dual peak characteristics as briefly described above, with such a case interpreted as a drug with either multiple absorption sites (Gabrielsson and Eeiner 2000) or a drug that exhibits enterohepatic recirculation (Funaki 1999). In the case of such drugs, standard compartmental models, either single- or multiple-compartment, are not applicable. However, in the muscle residue levels, the overall pattern exhibited profiles similar to a two-compartment model. Hence, the non-compartmental analysis was applied to analyze both plasma and muscle levels. Non-compartmental analysis employed for the present BTP data was developed for easier application when there is uncertainty in true underlying compartments while compromising some accuracy of estimation (Zhang et al. 2010). Thus, although it is possible to figure out general tissue patterns of BTP residues, the present results can be considered as a further validation of developed analytical methods rather than definitive estimation of pharmacokinetic parameters.

Plasma pharmacokinetic parameters indicated that BTP quickly reached peak levels and disappeared. The rapid disappearance of BTP was characterized by short MRT (< 11 days) since MRT is an estimate of how long a drug molecule stays in the body before complete elimination (Durisova 2014). Although only limited information is currently available, the pharmacokinetic behavior of BTP in the olive flounder appeared qualitatively similar to that in other animal species. In piglets, for example, BTP achieves peak levels as quickly as 18.6 min, with a reasonably efficacious absorption rate of 74.7% (Sun et al. 2016), and rapid plasma. Furthermore, rapid clearance from the body was observed in horses, demonstrating a decline in plasma BTP from 87.8 μg/ml peaks to undetectable levels in 24 h (Du et al. 2004).

Furthermore, experiments were performed at two different water temperatures: 22 °C as the optimal temperature, and 13 °C representing the low temperature at which the fish eats feed. Notably, higher water temperature conferred lower BTP peak concentrations and lower bioavailability compared with the lower temperature. This temperature dependence could be linked to a more efficient elimination when the fish is maintained at a higher temperature. In fishes, similar temperature dependency has been reported with the antimicrobial agent, oxolinic acid (Kleinow et al. 1994).

Another important objective of pharmacokinetic studies is consumer safety, particularly when the drug is administered to a food animal. Undoubtedly, olive founder is one of the most important aquaculture species in Korea, and BTP is currently employed with this fish (National Institute of Fisheries Sciences 2018). According to the related regulations (Korea Ministry of Food and Drug Administration 2017), the maximally allowable drug residue levels need to be identified in edible tissues, the muscle in the case of fish. At harvest for human consumption, fish should not contain drugs exceeding the allowable concentrations (maximum residue limit, MRL). Withdrawal time, the time after which harvesting is possible after the final drug use, will then depends on residue levels. Due to these food safety reasons, determination of both MRL and withdrawal time are crucial for all animal drugs used in aquaculture. Residue levels of an aquaculture drug are monitored for the purpose of establishing consumer safety. According to the analysis by the European Agency for the Evaluation of Medicinal Products (European Medicines Agency 2013), daily intake of BTP 600 μg/kg is safe for lifetime use in humans. The estimation of a withdrawal period was not possible as peak tissue levels were too low compared with the MRL.

Following the use of BTP in farms, the residue levels in olive flounder muscles may concern consumers regarding its safety in humans. The safety estimation based on muscle concentrations and the theoretical tolerable level indicates that BTP is extremely safe and that such concern is not legitimate.


In conclusion, a sensitive quantification technique for BTP was established using LC-MS/MS analytical procedures. Based on this method, established particularly for olive flounder plasma and muscles, it was possible to delineate the pharmacokinetic behavior of BTP in the fish.


The work was supported by a Grant from the National Institute of Fisheries Science, Republic of Korea (Grant no. R2019058).

Authors’ contributions

JHL, JSB, CWL, CYY, and JSC performed animal experiments and chemical analysis. SHC and YJK designed the study and guided chemical analysis procedures. KHP formulated the whole study. The author(s) read and approved the final manuscript.


The study was supported by the National Institute of Fisheries Science, Republic of Korea (Grant no. R2019058).

Availability of data and materials

Not applicable.

Ethics approval and consent to participate

Experimental protocols were approved by the Animal Care and Use Committee of the Kunsan National University, Korea (Approval number, AE 2016-04).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.



Alkobaby AI. Effects of maternal injection with organic phosphorus and vitamin B12 on reproductive performance and newly hatched offspring of Nile tilapia (Oreochromis niloticus). Proceedings, The 8th International Symposium on Tilapia in Aquaculture, October 12-14, 2008, Cairo, Egypt; Accessed 20 Jan 2010.


Benez LZ. Pharmacokinetic parameters: which are necessary to define a drug substance?. Eur J Respir Dis Supp. 1984; 134:45-61.


Deniz A, Spiecker-Hauser U, Rehagen M. Efficacy of a butafosfan and vitamin B12 combination (Catosal®) on biochemical and hematological blood parameters in dogs treated with dexamethasone. Intern J Appl Res Vet Med. 2009;7:116–29.


Du XY. Studies on pharmacy and pharmacokinetics of compound butafosfan. Ph.D Thesis, China Agricultural University, 2004. (2004). Accessed 20 Jan 2020.


Durisova M. A physiological review and structure of mean residence times. Gen Physiol Biophys. 2014; 33:75-80


European Medicines Agency. Opinion of the committee for medicinal products for veterinary use on the establishment of maximum residue limits: name of the substance-butafosfan. EMEA/V/MRL/003610/EXTN/0003 (2013). Accessed 20 Jan 2020.


Funaki T. Enterohepatic circulation model for population pharmacokinetic analysis. J Pharm Pharmacol. 1999; 51:1143-1148


Fürll M, Deniz A, Westphal B, Illing C, Constable PD. Effect of multiple intravenous injections of butaphosphan and cyanocobalamin on the metabolism of periparturient dairy cows. J Dairy Sci. 2010; 93:4155-4164


Fürll M, Wittek T, Gegenbach S, Schmidt B. Effects of preoperative application of butaphosphan and cyanocobalamin on reconvalescence, clinic-chemical parameters, and antioxidative metabolism and postoperative abomasal emptying in cows with abnormal dislocation. Tieraertl Prax. 2006; 34:351-356.


Gabrielsson J, Weiner D. Pharmacokinetic and pharmacodynamic data analysis: concepts and applications, 3rd eds. 2000; Sweden: Swedish Pharamceutical Press, Stockholm. p. 645-659.


Hasi S, Du X, Zhu B, Jiang J. Studies on effects of compound butaphosfan solution on endurance capability and energy metabolism in mice. Chin J Anim Vet Sci. 2004; 35:290-294.


Hasi S, Jiang J, Du X, Zhu B. Anti-cold stress effects and mechanisms of compound butaphosphan solution. Prog Vet Med. 2005; 26:59-62.


Hasi S, Jiang J, Zhu B, Du X. Studies on anti-heat stress effects and mechanisms of compound butaphosfan solution. Acta Vet Zootech Sin. 2005; 36:1334-1338.


Hekman P. Withdrawal-time calculation program WT1.4. 1996 London, UK: European agency for the evaluation of medicinal products (EMEA). .


Kleinow KM, Jarboe HH, Shomaker HE. Comparative pharmacokinetics and bioavailability of oxolinic acid in channel catfish (Ictalurus punctatus) and rainbow trout (Oncorrhynchus mykiss). Can J Fish Aquat Sci. 1994; 51:1205-1211


Korea Ministry of Food and Drug Administration (2017) Principles of Standard Establishment for Foods. p. 155. Accessed 20 Jan 2020.


Korea Rural Economic Institute. 2017-Food Balance Sheet, 2017.


Mechesso AF, Kim Y, Park S. Effects of butaphosphan and cyanocovalamin combination on plasma immune and biochemical parameters of olive flounder (Paralichthys olivaceus) subjected to crowding stress. Aquacult Res. 2019; 50:1611-1617


Meyer RE, Fish RE. Pharmacology of injectable anesthetics, sedatives, and tranquilizers. In: Anesthesia and analgesia in laboratory animals, 2nd eds. San Diego: Academic Press; 2008. 27–82.


National Institute of Fisheries Sciences. Aquatic medicine catalogue. (2018) Accessed 20 Jan 2010.


Pereira RA, Fensterseifer S, Barcelos VB, Martins CF, Schneider A, Schmitt E, Pfeifer LFM, Del Pino FAB, Corrêa MN. (2013b) Metabolic parameters and dry matter intake of ewes treated with butaphosphan and cyanocobalamin in the early postpartum period. Small Rum Res. 2013; 114:140-145


Pereira RA, Montagner P, Schneider A, Schmitt E, Rabassa VR, Pfeifer LFM, Del Pino FAB, Pulga ME, Corrêa MN. Effect of butaphosphan and cyanocobalamin on postpartum metabolism and milk production in dairy cows. Animal. 2013; 7:1143-1147


Rolling E, Berghaus RD, Rapnicki P, Godden SM, Overton MW. The effect of injectable butaphosphan and cyanocobalamin on postpartum serum beta-hydroxybutyrate, calcium, and phosphorus concentrations in dairy cattle. J Dairy Sci. 2010; 93:978-987


Seo JS, Lee JH, Park JJ, Choi JS, Bae JS, Lee CW, Yang CY, Kang YJ, Choi SH, Park KH. Biochemical and stress-attenuating effects of butaphosphan-cycnocobalamin combination drug in olive flounder Paralichthys olivaceus. Fish Sci. 2019 (published online December 16).


Sun F, Wang J, Yang S, Zhang S, Shen J, Xingyuan C. Pharmacokinetics of butafosfan after intravenous and intramuscular administration in piglets. J Vet Pharmacol Ther. 2016; 40:203-205


Tabeleão VC, Schwegler E, Pereira RA, Krause ART, Montagner P, Feijó JO, Schneider A, Schmitt E, Brauner CC, Rabassa VR, Del Pino FAB, Corrêa MN. Combined butaphosphan and cyanocobalamin on the glucose metabolism of dairy cows after calving. Arq Bras Med Vet Zootec. 2017; 69:317-324


Van de Vis H, Kruijt-Kloosterboer K, Veldman M. Test for potential protective effects of Catosal® on mortality of fish due to temporary hypoxia. Internal Report No. C061/04, RIVO-Netherlands Institute for Fisheries Research, Ymuiden, Netherlands. (2004). Accessed 20 Jan 2020.


Van der Staay DE, Groot J, Van Reenen CC, Hoving-Bolink AH, Schuurman T, Schmidt BH. Effects of butaphosphan on salivary cortisol and behavioral response to social stress in piglets. J Vet Pharmacol Ther. 2007; 30:410-416


Zhang Y, Hou M, Zhou J, Xie S. PKSolver: an add-in program for pharmacokinetic and pharmacodynamics data analysis in Microsoft Excel. Computer Methods Programs Biomed. 2010; 99:306-314





Dissolve oxygen


Liq2uid chromatography coupled with tandem mass spectrometry


Multiple residue monitoring


Maximum residue limit