A strategy for trade monitoring and substitution of the organs of threatened animals

The use of threatened animals as a source of traditional medicines is accelerating the extinction of such species and imposes great challenges to animal conservation. In this study, we propose a feasible strategy for the conservation of threatened medicinal animals that combines trade monitoring and the search for substitutes. First, DNA barcoding provides a powerful technique for monitoring the trade of animal species, which helps in restricting the excessive use and illegal trade of such species. Second, pharmacological tests have been adopted to evaluate the biological equivalence of threatened and domestic animals; based on such testing, potential substitutes are recommended. Based on a review of threatened animal species and their substitutes, we find that the search for substitutes deserves special attention; however, this work is far from complete. These results may be of great value for the conservation of threatened animals and maintaining the heritage of traditional medicine.

D espite legislation from the Convention on International Trade in Endangered Species (CITES) and the establishment of the International Union for Conservation of Nature (IUCN) Red List 1 , threatened animal organs are still widely used in many countries, such as those of East and Southeast Asia 2,3 . The excessive use of animal organs as a source of traditional medicine (TM) is one of the main reasons why a number of species are threatened or endangered 4 . Some TMs are known to include trade-restricted species of animals as ingredients and therefore contravene the CITES legislation. For example, rhinoceros horn powder 5 , saiga antelope horn powder 6 and bear bile crystals 7 are commonly used in traditional Chinese medicines (TCMs) and occupy a special position in traditional Asian medicine; however, the biological origins of species have been decreasing rapidly in recent decades due to their increased use. Several international unions, such as the World Society for the Protection of Animals and the Wildlife Conservation Society, have attached great importance to these problems. Ascertaining the biological origin of these materials, enforcing legislation and prosecuting cases of illegal trade are important ways of conserving wildlife. To our knowledge, ingredients containing animal organs are presented mostly in the form of powders, tablets and capsules, and the biological origins of these materials are difficult to recognize using regular methods, such as chromatographic analysis, after they have been processed 8 . DNA barcoding, with its unique reproducibility, sequence versatility, and comparability among different species, offers a powerful approach to genetically authenticate animal organs [9][10][11] . However, these traditional medicines have been the primary treatment for more than 3,000 years and have been accepted by a large portion of the population; therefore, they are difficult to eliminate 12 . Seeking substitutes that possess equipotent bioactivities from close phylogenetic species has been given much attention worldwide by ecologists, pharmacologists, economists and sociologists 13,14 ; these efforts will greatly alleviate the pressure on threatened species. Liu et al. investigated substitutes by elucidating the active components on specific pharmacological activity, e.g., by purifying novel antioxidant peptides 15 . However, species authentication and comparative pharmacology of endangered animal-derived drugs (ADDs) and their substitutes should be considered to verify the reliability of possible substitutes.
This study aims to propose a strategy for wildlife conservation by combining trade monitoring and the search for substitutes (Fig. 1). Fifteen animal species were investigated, including three CITES-listed species (see the appendices at www.cites.org): saiga antelope (Saiga tatarica, listed in Appendix II), yellow-margined box turtle (Cuora flavomarginata, listed in Appendix II) and Chinese pangolin (Manis pentadactyla, listed in Appendix II); five state-protected species: black rhinoceros (Diceros bicornis, listed as critically endangered in the IUCN Red List), Père David's deer (Elphurus davidianus, listed as extinct in the wild in the IUCN Red List), sika deer (Cervus nippon, listed at the first state protection level in China), red deer (Cervus elaphus, listed at the second state protection level in China), Chinese pond turtle (Chinemys reevesii, listed as endangered in the IUCN Red List); and seven domestic species: sheep (Ovis aries), domestic goat (Capra hircus), Asian water buffalo (Bubalus bubalis), Chinese softshell turtle (Trionyx sinensis), hawksbill (Eretmochelys imbricate), yellow-headed tortoise (Indotestudo elongata), and yellow-bellied slider turtle (Trachemys scripta). The samples genetically audited in this study using DNA barcoding are shown in Fig. 2. DNA barcoding was applied to authenticate the animal species, which involved the DNA sequencing and characterization of mitochondrial DNA (mtDNA) genes. This technique is able to identify threatened species and domestic species with confidence, thereby offering a reliable tool for customs officials. In addition to concerns about trade monitoring, the lack of substitutes has also affected threatened animal species. For this reason, we adopted various pharmacological methods to evaluate the bioactivities of several ADDs using their biological fingerprints and characteristics. We expected to find equivalence regarding the specific functions of ADDs obtained from threatened species to derive substitutes for use in complementary and alternative medicines. The design of the present study is depicted in Fig. 1.

Results
Sequence information, polymerase chain reaction amplification capacity, and sequencing success rate. Polymerase chain reaction (PCR) amplification success and sequencing success are important indices for evaluating DNA barcodes. In the present study, both PCR amplification and sequencing achieved a 100% success rate. The PCR amplification electropherogram was obtained using agarose gel electrophoresis. Using DnaSP version 5.10.01, 72 sequences ( Fig. 3) were analyzed, and 477 sites (i.e., base pairs) were obtained; 463 sites were without gaps or missing data. Of these, 156 sites (32.7%) were polymorphic; the nucleotide diversity (pi) was 10.9%. The average   nucleotide composition of all the sequences was calculated in MEGA 4.0 16 as follows: T 5 31.3%; C 5 24.8%; A 5 28.5%; and G 5 15.4%; the average AT content was 59.8%.
Intra-and interspecific genetic divergence analysis. For successful DNA barcoding, the sequence variation between species needs to be sufficiently high so that they can be distinguished from one another; however, the variation must be sufficiently low within species that a clear threshold between intra-and interspecific genetic variations can be defined. In the present study, the mean intra-and interspecific genetic variations were 0.4% and 15.8%, respectively. All intraspecific variations measured between sequences were less than 1.8%, and all interspecific genetic variations were greater than 2.0%.
Phylogenetic analyses. A cladogram was clustered based on neighborjoining analysis and demonstrated that all ADD haplotypes were strongly grouped with their orthologous mtDNA. We also noted that with eight species, the ADDs were clustered in four subfamilies in the dendrogram (Fig. S2).
To test the accuracy of the barcodes and expand our identification capability, we downloaded all complete mitochondrial cytochrome c oxidase I (COI) genes available in GenBank (November 2011). All 241 sequences (Table S1) were cropped to the 477-base pair barcode region amplified by the primer pair MahF1 and MahR1. A phylogenetic analysis of 313 sequences was conducted, which included the GenBank dataset and our test sample dataset (Fig. 3). The cluster branches indicated that all haplotypes strongly grouped with their orthologous mtDNA and every species clustered in the appropriate branch. Based on this analysis, we were able to deduce that this region of the COI gene is a good candidate for distinguishing between the ADDs and between the taxonomy of related species.
Barcoding gap test and identification capability. Ideal DNA barcodes should contain a conspicuous spacer region that delineates intra-and interspecific variations -the barcoding gap 17,18 . In the present study, a visible barcoding gap was found in only a small portion of the ADDs with overlapping intra-and interspecific variation (Fig. S3); this indicated that all ADDs could be accurately identified.
Comparing species-identification capability relative to other genetic fragments is one of the criteria for evaluating the quality of various barcodes. In this study, the taxonomic affiliation of each sequence (including those of 72 test samples and 241 individuals from GenBank) was determined using BLAST1 against the GenBank database. This analysis indicated that all mitochondrial sequences were successfully identified.
In vivo study of endotoxin-induced rabbit fever. The effects of the ADDs on the rabbit fever models are shown in Fig. S4A and Table S2. As is evident from Fig. S4A, the blank control had a wave range of 0.1uC, which indicates a successful fever model. Compared with negative controls, saiga antelope horn, pangolin scale, goat horn, and aspirin significantly reduced the temperature of fevered rabbits; velvet antler and antler significantly enhanced the temperature of fevered rabbits. Saiga antelope horn exhibited the strongest anti-fever effect, followed by pangolin scale and goat horn.
The endotoxin, PGE1, PGE2, and TNF-a levels in the serum of rabbits were determined to establish the effective mechanisms of the ADDs. Compared with negative controls, saiga antelope horn and goat horn clearly reduced the levels of PGE1, PGE2, and TNF-a, although velvet antler caused significant increases (Figs. S4B-S4D). Blood endotoxin content was higher when using velvet antler compared with the other ADDs (see Table S2). It has been established that endotoxins stimulate monocytes and macrophages to produce proinflammatory cytokines, including TNF-a and IL-1b 19,20 , which subsequently cause fever. We demonstrated that saiga antelope horn  and goat horn inhibited the production of TNF-a and PGE2, thereby relieving fever.
In vitro study of liver microsomes. Thermogenic curves of the growth of liver microsomes recorded at 37uC under the action of ADDs are shown in Fig. S5A. The heat-flow power (HFP) time curve could be delineated using the following equation 21 : P t 5 P 0 exp(kt) or lnP t 5 ln P 0 1 kt, where P 0 and P t represent the HFP at time 0 and t (min), respectively. From this equation, the growth-rate constants (k) of the exponential phase for microsomal metabolism at 37uC in the absence of any substance could be calculated using the data obtained for the highest peak. Then, quantitative thermokinetic parameters, such as the HFP of the highest peak (P), the appearance time of the highest peak (t), and the total heat output (Q t ) were obtained from the HFP time curve of the microsomal growths affected by the various ADDs (Fig. S5B).
The inhibition ratio (I) of ADDs on microsomal growth was calculated from k. Using principal component analysis 22 , P and Q t were calculated as the main components, accounting for 91.89% of the variability. The inhibition ratio was calculated using the following equation: I 5 (P t 2 P 0 )/k 0 , where P 0 and P t represent the HFP in the absence of any substance and in the presence of the ADDs, respectively. By analyzing the main components, we noted two findings. Saiga antelope horn and goat horn inhibited the heat output power and the total heat output, which inhibited microsomal growth. In addition, velvet antler and antler enhanced the heat output power and the total heat output, which indicated increased heat output by the microsomes.
In vitro study of spleen lymphocyte proliferation. The thermogenic curves at 37uC in the presence of ADDs are shown in Fig. S6E. Based on the liver microsomal study, k, P, t, and Q t were obtained; among these, P and Q t were identified as the main components. From an analysis of these main components, we determined that velvet antler and antler increased P and Q t to varying degrees, and saiga antelope horn, pangolin scale, goat horn, and buffalo horn clearly inhibited the heat output of the spleen lymphocytes (Figs. S6F, S6G). These results agree with the theory of the ''nature of heat and cold'' in traditional Chinese medicine 23 . Because lymphocytes are closely related to immunization 24,25 , we believe that ADDs differ in their immunomodulatory effects. To a certain extent, we also consider that ADDs from the horns of threatened animal species can be partially substituted by other substances.
In vivo study of thermotaxis in mice. The nature of heat and cold is an intrinsic property of matter 23 . In this part of the study, animal behavior was investigated to explain the nature of ADDs. Mice were randomly divided into 6 treatments, consisting of 10 replicates. Investigations were carried out under the following temperature gradients: 15-30uC, 15-40uC, 20-30uC, 20-40uC, 25-40uC, and 30-40uC. Under these conditions, the mice exhibited an obvious temperature tropism: they congregated in either the warm zone or the cool zone. At 15-25uC, 15-35uC, 20-35uC, or 25-35uC, the mice showed no obvious temperature tropism: the numbers of animals in the warm and cool zones were similar (Fig. S7). For mice not exposed to ADDs, equal numbers were found in the two zones for the first five hours. However, beyond that time, the numbers of animals in the zones changed significantly, as did the 30-min movement distances and the number of times the animals crossed from one zone to another (Fig. S8). Thus, the results were repeatable within the first five hours.
The results of thermotaxis in mice exposed to ADDs showed that velvet antler and antler reduced the number of animals in the warm zone and increased the 30-min movement distance and the number of times the mice crossed from one zone to another (Figs. S9D -S9F). In contrast, saiga antelope horn, goat horn, and pangolin scale promoted the movement toward heat. There was no significant difference in the effects of saiga antelope horn, goat horn, and pangolin scale on thermotaxis (Figs. S9D -S9F). In addition, velvet antler and antler increased the autonomic action of the mice, whereas saiga antelope horn, goat horn, and pangolin scale decreased this action.
In addition to studying animal behavior, we determined the levels of 5-hydroxytryptamine (5-HT), cyclic adenosine monophosphate (cAMP), and cyclic guanosine monophosphate (cGMP) in the serum and brain tissues of the mice. 5-HT is an inhibitory neurotransmitter in the vertebrate nervous system that plays an important role in controlling mental states, perception, and muscular fatigue. We found that saiga antelope and pangolin scale significantly enhanced the level of 5-HT in both the brain tissues and sera of mice compared with blank controls. However, velvet antler, antler, and goat horn showed no significant effect on 5-HT levels ( Fig. S9G; Fig. S10). cAMP and cGMP are cyclic nucleotides that exert antagonist effects in vivo 26 . The balance between cAMP and cGMP is essential for cellular metabolism and plays a role in many pathological processes 27 . We determined that velvet antler and antler significantly enhanced the level of cAMP in both the brain tissues and sera of mice compared with blank controls (Fig. S9H); velvet antler also clearly enhanced the level of cGMP in both the brain tissues and sera of mice (Fig. S9I). Conversely, saiga antelope horn, goat horn, and pangolin scale decreased cAMP and cGMP levels to various degrees (Figs. S9H and S9I). Moreover, saiga antelope horn significantly increased cAMP/cGMP values, whereas velvet antler decreased these levels (Fig. S9J).
In vitro study of thrombus formation 28 . We evaluated the rabbit plasma recalcification and thrombin times of ADDs; saiga antelope horn, pangolin scale, yak horn, and buffalo horn significantly increased the rabbit plasma recalcification (p , 0.01) and thrombin (p , 0.01) times; in addition, these drugs exhibited good anticoagulant effects. Antler, goat horn, and turtle shell also increased the rabbit plasma recalcification (p , 0.05) and thrombin (p , 0.05) times. Saiga antelope horn and Mongolian gazelle horn had no obvious effect on thrombus formation (p . 0.05).

Discussion
Threatened animal conservation is a special field within applied ecology. Currently, the objectives of wildlife management fall into the following categories: protecting the natural habitats of animals through controlled exploitation; maintaining protected areas containing rare species; establishing biosphere reserves for endangered species; protecting wildlife through legislation; and imposing specific restrictions on the export of endangered animals or their products. However, it is undeniable that many drugs derived from threatened animals have remarkable curative effects in traditional medicine. We consider that the most effective way of protecting medicinal wild animals is to combine trade monitoring and seeking substitutes.
Regarding trade monitoring, the techniques applied include macroscopic and microscopic morphology, histology and molecular biology; however, these methods fail to identify ADDs with confidence. In recent decades, molecular techniques have promoted the identification of ADDs; these techniques include protein electrophoresis, restriction fragment length polymorphism analysis, and using specific primers such as those for 12S rDNA, 16S rDNA, valine tRNA and Cyt B to amplify specific regions in the mitochondrial genome. However, these markers have no universal primers for broad animal species. Despite the great potential of genetics to assist in the identification of ADDs and the phylogeny of some animal clades, such as the Bovidae, Cervidae, and Moschidae 29-31 , little consensus exists regarding which gene region is the most suitable for studying a wide range of animal species. Recently, the mitochondrial cytochrome c oxidase I (COI) gene has attracted global attention as a DNA barcode www.nature.com/scientificreports SCIENTIFIC REPORTS | 3 : 3108 | DOI: 10.1038/srep03108 for animals because the consortium for the barcode of life (CBOL) uses universal primers to amplify a region of the COI gene approximately 650 bp in length to identify broad animal taxonomy. Miller 32 considered that the DNA barcode technique might usher taxonomic research into a new era. The COI gene is sequenced to yield a DNA barcode for the examined specimen for use in trade monitoring; this barcode is compared to barcodes from reference specimens to identify the species. Our results indicate that the COI gene is a useful barcode for the trade monitoring of threatened animal species, in agreement with a previous report 33 demonstrating that even short COI fragments can efficiently identify antelopes, buffalo and domestic Bovidae; thus, this technique exhibits great potential for use in wildlife conservation. DNA barcoding can serve as a powerful tool in wildlife forensics 34,35 and might prove to be a vital aid in conserving organisms that are threatened by the illegal wildlife trade, such as turtles 36 . The COI gene has also enhanced our understanding of mammalian diversity in Southeast Asia and, in this way, aided conservation planning 37 . The results obtained here confirm that the COI gene exhibits good universality and a powerful ability to identify animal species; therefore, this gene should play an important role in monitoring the trade in endangered animals. However, the heavy reliance on one gene for identification has its shortcomings for wildlife forensic enquiries, so that we suggest that if COI gene could not be obtained, other protein coding gene may be a good backup.
Animal products play an important role in traditional medicine within China 38 , Japan 39 , Brazil 40 and many other countries. Although some animal species are becoming threatened or endangered due to their excessive use as complementary medicine, few measures have been proposed for the sustainable development of traditional medicine. Emphasizing conservation biology while disregarding the traditional use of wildlife in medicine will therefore encounter great resistance. One feasible solution is to substitute the organs of closely related domestic animals by examining their biological properties and ranking the most relevant indices regarding their medicinal values. In the present study, several pharmacological tests were performed on ADDs obtained from both threatened and domestic animals, and the resulting bioactivities (both in vitro and in vivo) were compared in several ways. Fig. 4 presents several regularities. 1) Velvet antler and antler show opposite tendency compared with other ADDs, e.g., saiga antelope, domestic goat, Chinese pangolin and Chinese softshell turtle reduce the temperature of fevered rabbits while velvet antler and antler enhance that. This regularity is validated by the tests of liver microsome, spleen lymphocytes and thermotaxis. 2) The levels of the bioactivities among ADDs vary in different tests. For example, saiga antelope and domestic goat share the same level of activity on temperature and thermotaxis tests; while the activity of saiga antelope is better than that of domestic goat on liver microsome (p , 0.05) and spleen lymphocytes tests (p , 0.05). Another example is that velvet antler and antler share the same level of activity on spleen lymphocytes test while different regarding other tests. 3) Two new bioactivities are found including the antipyretic activity of Chinese pangolin and the anticoagulant activity of velvet antler. 4) ADDs derived from threatened species that may be replaced by ADDs obtained from domestic animals (designated by different colors in the figure) would otherwise be difficult to substitute. We conclude that domestic goat may be the potential substitute of saiga antelope regarding its immunological and antipyretic activities, and Chinese softshell turtle may be the potential substitute of Chinese pangolin regarding its immunological activity. These potential substitutes could alleviate the burden on threatened species and enrich the traditional medicines. 5) Although the bioactivities of antler are not so good as velvet antler, the increase dosage of the antler may achieve the same potency. This substitutive use can reduce harm to deer since antler www.nature.com/scientificreports SCIENTIFIC REPORTS | 3 : 3108 | DOI: 10.1038/srep03108 is ossific angle that will shed naturally once a year while velvet antler is obtained by sawing.
Pharmacological studies indicate that each ADD usually has several therapeutic benefits and that only certain actions of threatened-ADDs are interchangeable. Based on the results of this study, we are able to establish a criterion for the use of ADD substitutes: if a substitute shows the same level of bioactivity as that of a threatened ADD for specific functions, it may be adopted as a candidate for potential clinical use. We believe that further study in this area along the suggested lines would contribute to the use of substitutes. However, it should be noted that the replacement of one particular function does not equate to complete replacement. Thus, when examining possible candidates for substitution, all of the main functions or bioactivities of the ADDs need to be examined and clinical trials should be conducted where appropriate 41 . Additionally, the mechanisms of the ADD bioactivities have to be clarified to aid in identifying possible substitutes 42,43 : ADDs derived from domestic species that share the same mechanisms as those obtained from threatened species should be given the highest priority.
But et al. pointed out an important aspect for the remedies of endangered animal species by using medicinal plants, which are easier to obtain 5 . If so, it should be good news for conservation of wild animals. However, to our knowledge, medicinal plants were only adopted by few cases to replace animal-based remedies. Feng et al. reported that plant substitutes were only suggested in literature without solid evidence to support, while some animal substitutes were superficially studied 44 . Even if efforts on finding the alternative domestic animals for the endangered animals are also not many, we added Table S4 to summarize current investigation. In addition, when it is not possible to replace the medicinal value of ADDs from threatened species with those from domestic species, efforts should be made to examine the relevant active components in the ADDs and investigate any interactions among the various constituents; this would be useful in the attempt to synthesize alternatives. Otherwise, it would be necessary to direct funding toward those species that are at the greatest risk of extinction. However, this strategy is not the most efficient way to promote the recovery of such species or reduce global extinction rates; some of the species at the greatest risk of extinction would require a massive recovery effort, most likely with a small chance of success 45 .
We have summarized the trade monitoring and substitution of threatened animal species (Table S4), which has attracted the strongest concern worldwide; species at risk include cetaceans, carnivora, proboscidea, pholidota, perissodactyla, artiodactyla, the reptiles testudines and the fish pleurotremata. The threatened species are listed in order according to their position on the IUCN Red List of threatened species. Among the listed species, 42 are included in the CITES appendices, which afford different levels of protection from successive use 46 . Table S4 demonstrates that studies of the substitution of threatened animal species have just begun, and many critically endangered species, such as Père David's deer, saiga antelope and hawksbill turtle, cannot be substituted and have no alternative uses. Many articles have been written about the use of threatened animal species in traditional medicines [47][48][49][50][51] , and the illegal trade in these species for medicinal use will possibly destroy biodiversity [52][53][54] . Therefore, the development of substitutes for drugs derived from threatened species, critically endangered species in particular, is of great significance for wildlife conservation.
The use of substitutes not only relates to the bioactivities of these drugs, it also involves traditional and cultural concerns. The functions of organs threatened animals are often exaggerated based on value orientations, religious beliefs, superstitions and totem culture, and the blind faith of some high-end consumers, with complete disregard for the exploitation of wild animal resources.

Methods
Taxon sampling. A total of 72 animals from 15 species that are used as sources of ADDs were collected in six provinces and regions of China (Beijing, Jilin, Sichuan, Qinghai, Xizang, and Hubei). The sample information is depicted in Fig. 2. Among the samples, those derived from wild animals, including saiga antelope horn, Cervus nippon velvet antler and Cervus nippon antler, were provided by the Animal Research Institute of the Chinese Academy of Sciences and China's Forestry Sciences ( Table 1). The distribution information was derived from investments made in TM in China (Fig. S1). We also analyzed 241 animal sequences within four related families from GenBank. We received a license from the State Forestry Administration of China to study wild animals, and the license number is Lam Woo Approval [2005] 627; this license is listed on the official website of the State Forestry Administration of China (http://www.forestry.gov.cn). Throughout this study, none of the authors handled any animals. Although we investigated animal products, we had no direct contact with the animals themselves. The samples were collected by administrators in the National Nature Reserve and at aquaculture bases.
DNA extraction. ADDs was frozen at 24uC and subsequently ground in a DNA extraction beveller (Retsch MM400, Germany) for 1 min at 1,800 rpm. The total DNA was then extracted using DNAout kits. Specifically, 200 mg of each ADD sample was dispersed in 1 ml DNAout solution, and 10 ml RNase A and a little silicon were subsequently added. After incubating for 5 minutes at 65uC, the solution was centrifuged at 15000 3 g for 5 minutes, and the supernatant was extracted twice with chloroform, each time 0.2 ml. Then, an equal volume of isopropyl alcohol was added into the supernatant, shaken up, and centrifuged at 15000 3 g for 5 minutes. The precipitation was washed twice with ethanol, each time 1 ml, and then dissolved in buffer solution for use. The PCR primers were designed from a ClustalW alignment of COI sequences available from GenBank for 54 closely related species within the Bovidae, Cervidae, Geoemydidae and Trionychoidae families. The forward primer (MahF1 5'-GCAGGAACAGGCTGAACCGT-3') was anchored at site 364 (the numbering used is relative to the complete CO1 sequence from the above four families); the reverse primer (MahR1 5'-AATATGTGGTGGGCTCATAC-3') was anchored at site 863. The performance of MahF1 and MahR1 in PCR amplification was tested using genomic DNA isolated from 10 Bovidae, Cervidae, Geoemydidae, and Trionychoidae species. Preliminary PCR amplifications using MahF1 and MahR1 achieved a 100% success rate. The animals were raised separately by gender and had unlimited access to food and water in an environmentally controlled breeding room (temperature 22 6 2uC, humidity 60 -80%). The breeding room was illuminated by an artificial light cycle with 12 hours of light and 12 hours of darkness every day and was disinfected regularly.
Sample preparation for the pharmacological study. Because of the limitation of the amount of some ADDs samples, only six ADDs derived from a total of eight species were tested for pharmacological studies. Each ADD was weighed (approximately 300 g) and immersed in cold water for 1 hour. Then, 1,500 ml of water was added; this was decocted three times, once an hour. Next, the decoction solutions were combined and concentrated to 1,000 ml to achieve a final solution equivalent to 0.3 g ingredient ml 21 . The concentrate was then sealed and sterilized using high-pressure steam at 121uC for 30 min and stored at 4uC. In the pharmacological study, all ADDs underwent the same preparation process, and various concentrations were obtained by diluting the mother liquid with distilled water.
In vivo study of endotoxin-induced rabbit fever. Healthy New Zealand pure big-ear white rabbits were subjected to frequent temperature measurements. Animals were selected whose temperatures were 38.6-39.5uC and did not vary more than 0.3uC. The selected rabbits were randomly allocated into nine groups with six rabbits in each group (including six experimental groups, a positive control and blank control). Escherichia coli endotoxin (O111B4) was slowly injected into the ear vein of rabbits with the concentration of 12 EU/kg to build fevered rabbits models. After one hour, ADDs (5.0 ml mother liquid kg 21 ) were administered to rabbits and the change in the rabbits' body temperature was measured to monitor the effects of ADDs. Positive controls were administered 200 mg kg 21 aspirin, and negative controls were administered the same volume of water and endotoxin in the ADD groups; blank controls were given only the same volume of water in the ADD groups. ADDs were administered to rabbits 1 hour before fever was induced, and the initial temperatures were immediately measured. The body temperature of the rabbits was measured at 0.5, 1.0, 1.5, 2.0, 3.0, 4.0, 5.0, and 6.0 h after fever was induced. Ear-vein blood samples were taken from the rabbits, and PGE1, PGE2, and TNF-a levels were measured by radioimmunology. The endotoxin content was determined using a dynamic nephelometric method. In pre-interference experiments, the same amount of endotoxin was added, and the recovery rates were calculated (Table S2).
In vitro study on spleen lymphocyte proliferation. Ten healthy Kunming mice were weighed and immediately sacrificed under anesthesia. The spleens were obtained, and the residual blood and polyp tissues attached to the spleen were washed away using phosphate-buffered saline; dispersed cells were prepared by passing the tissue through a cell sieve (mesh size, 75 mm). The cell suspension was then centrifuged for 5 min at 1,000 3 g, and the pellet was suspended in Dulbecco's modified Eagle's medium culture to a density of 5 3 10 6 cell ml 21 . Then, 100 ml of the cell suspension was added to each well of 96-well culture plates, and 100 ml aliquots of Dulbecco's modified Eagle's medium containing various doses of the ADDs were added to each well. Finally, the 96-well culture plates were incubated for 72 h at 37uC in 5% CO 2 , and the cytopathogenic effect was microscopically observed. Preliminary study showed that ADDs did not cause damage to spleen lymphocytes at concentrations lower than 40 mg ml 21 (Table S3), and 10 mg ml 21 was finally selected for the present experiment. Microcalorimetry 58,59 was used to determine the HFP curves for spleen lymphocyte proliferation. To achieve a cell density of 5 3 10 4 cell ml 21 , 5 ml of RPMI-1640 culture incubated in mouse spleen cells was added quantitatively to a 20-ml ampoule under sterile conditions. Then, 500 ml of RPMI-1640 culture containing one type of ADD (10 mg ml 21 ) was added to the ampoule. Ampoules containing various ADD and spleen cell suspensions were sealed and lowered into the measuring position of the calorimeter block. After approximately 30 min (by which time the temperature of the ampoules had reached 37uC), HFP time curves were recorded until the readings had returned to the baseline level. All data were continuously recorded using a dedicated software package (PicoLog TC-80, TA Corporation, USA).
In vitro study using liver microsomes. The liver tissues of the ten killed mice in above test were taken and placed into iced water. The fat tissue were quickly removed and discarded; a sufficient amount of the tissue was then homogenized with approximately 10 strokes using a glass Dounce homogenizer. A volume 10 times that of the mitochondrial extract (including 0.4 M mannitol, 50 mM tris-hydroxymethyl aminomethane-2 hydrochloric acid, 1 mM disodium ethylenediamine tetraacetate, 5 mM potassium chloride, 0.1% bovine serum albumin and 2 mM mercaptoethanol at pH 7.5) was required to ensure extraction efficiency. The homogenate was transferred to aseptic centrifuge tubes and centrifuged at 1,000 3 g at 4uC for 5 min. The supernatants were transferred to new centrifuge tubes and centrifuged at 12,000 3 g at 4uC for 10 minutes; the pellets were then collected and dissolved with phosphate-buffered saline for use as the mitochondrial suspension. The subsequent tests were determined using microcalorimetry, and the experimental protocol, groups and dosage were the same as the test of spleen lymphocyte.
In vivo study on thermotaxis in mice. Kunming mice were fed and monitored for 3 days before the experiment, and those with no obvious temperature change were used for further training so that they could memorize the various experimental zones. Afterwards, the mice were randomly allocated into seven groups, with ten mice in each group. The experimental groups were administered one type of ADD for 7 days at a dose of 5.0 ml of mother liquid kg 21 day 21 before recording behavioral observations; the negative control was given the same volume of water. After the traces had been recorded, the mice were sacrificed, and blood samples were obtained after enucleating the eyeballs. Immediately thereafter, brain tissue was extracted.
The samples of blood were centrifuged for 5 min at 1,000 3 g, and the resulting serum was stored at -20uC. The level of 5-HT in the blood was directly determined using a 5-HT ELISA Kit (Cusabio Biotech Co., USA). A mixture of 0.1 ml of serum with 2.0 ml of ethanol was evenly shaken and incubated at room temperature for 5 min before being centrifuged for 10 min at 1,500 3 g. The supernatant was transferred to a new centrifuge tube, and the pellet was added to 1 ml 75% ethanol, evenly shaken, and centrifuged for 10 min at 1,500 3 g; this second supernatant was collected and mixed with the first supernatant. The supernatant was dried at 60uC and dissolved with 1 ml of acetate buffer, which was then analyzed by radioimmunoassay.
Two milliliters of ethanol was added to 2 ml of brain homogenate (containing 50 mg of brain tissue), mixed well, incubated for 5 min, and centrifuged for 15 min at 1,500 3 g. The supernatant was collected, dried at 60uC, and stored at 4uC. Half of the residue was dissolved in acetate buffer, and the cAMP and cGMP levels were determined by radioimmunoassay. The remaining half was dissolved in normal saline, evenly shaken, centrifuged for 10 min at 1,000 3 g, and the 5-HT level was determined.