Literature DB >> 28852516

Compatibility of preparatory procedures for the analysis of cortisol concentrations and stable isotope (δ13C, δ15N) ratios: a test on brown bear hair.

Agnieszka Sergiel1, Keith A Hobson2,3, David M Janz4, Marc Cattet5,6, Nuria Selva1, Luciene Kapronczai7, Chantel Gryba2, Andreas Zedrosser8,9.   

Abstract

The measurement of naturally occurring glucocorticoids and stable isotopes of several elements has gained importance in wildlife studies in recent decades and opened a myriad of ecological applications. Cortisol and stable isotopes equilibrate in animal tissues over periods of integration related to the growth rate of the tissue, providing information reflecting systemic cortisol secretion and dietary intake. Sample preparation shares the common step of first cleaning the sample of external contamination. However, it is not well understood how different solvents used in sample preparation affect isotopic and cortisol values, and whether it is safe to follow the same procedures for both measures to optimize analyses of the same sample. We conducted an experiment to compare different preparation protocols for the analysis of cortisol concentrations and stable carbon (δ13C) and nitrogen (δ15N) isotope ratios in hair. Hair samples from 12 brown bears (Ursus arctos) were each divided into five aliquots; two aliquots were rinsed with a 2:1 chloroform:methanol (v/v) mixture with one aliquot ground prior to cortisol analysis and the other left intact for stable isotope analyses; two aliquots were washed with methanol with one aliquot ground prior to cortisol analysis and the other left intact for stable isotope analyses; and one aliquot washed with methanol and ground prior to stable isotope analyses. The cortisol, δ13C and δ15N values remained consistent following all treatments. Our results indicate that hair samples rinsed with a 2:1 chloroform:methanol mixture or washed with methanol can be used for both types of analyses. Further, hair that has been ground in a standard hair cortisol procedure can also be used for stable isotope analysis. This information is important for improving laboratory efficiency and compatibility of procedures used for wildlife physiological ecology studies where concurrent measurements of cortisol and stable isotopes in hair are required.

Entities:  

Keywords:  Dietary intake; Ursus arctos; laboratory procedures; stable isotopes; stress; wildlife physiological ecology

Year:  2017        PMID: 28852516      PMCID: PMC5569930          DOI: 10.1093/conphys/cox021

Source DB:  PubMed          Journal:  Conserv Physiol        ISSN: 2051-1434            Impact factor:   3.079


Introduction

Levels of stress hormones and stable isotope ratios in animal tissues are increasingly recognized as important tools in studies on the physiology and ecology of wildlife (e.g. Kempster ; Barger and Kitaysky, 2011; Deschner ; Bryan ). Keratinous tissues, such as hair and feathers, offer the longest record of an animal's circulating glucocorticoid (GC) concentrations (Bortolotti ; Macbeth ; Jenni-Eiermann ). Steroids are incorporated into the hair shaft and feather during their growth, and therefore, GC levels in those matrices are thought to reflect average systemic levels over the respective growth phase by integrating baseline levels and elevated adrenocortical secretion (Cattet ; Jenni-Eiermann ). Hair or feather GC concentrations are increasingly applied to evaluate chronic exposure to various stressors or potentially stressful conditions (Koren ; Davenport ; Accorsi ; Macbeth ; Comin ; Fairhurst ). Additionally, GCs in those sample matrices are stable over time and resistant to environmental exposure (Bortolotti ; Macbeth ). Therefore, measuring GC concentration provides a useful tool in understanding the role of stress in the life history, conservation physiology, health, and ecology of wildlife species (e.g. Wikelski and Cooke, 2006; Macbeth ; Cattet ). Stable isotopes are also incorporated over time into hair and feathers (Wassenaar, 2008; Bontempo ). Stable nitrogen (δ15N) and carbon (δ13C) isotope ratios are widely employed to reconstruct and assess temporal and spatial variation in diet using extant or archived materials (Newsome ; Crawford ). These isotopes are also used to characterize trophic niche and community structure (Boecklen ), elucidate animal migration and movements (Hobson and Wassenaar, 2008), and determine individual specialization and habitat selection (Newsome ). Stable isotope analyses are also employed in wildlife forensics and ecotoxicological studies (Bowen ). Combining concurrent measures of long-term stress and stable isotopes within the same study may open up a wide selection of testable ecological theories, including relationships between diet and stress (e.g. Bryan ; Fairhurst ; Lafferty ). Compared to blood, saliva, urine or feces, hair and feathers are media that can be transported and stored at room temperature. The use of hair and feathers for quantifying stress offers other desirable features, such as non-invasive collection and retrospective analyses of stress using archived samples (Cattet ). Additionally, some studies have made use of hair collected from museum and archeological specimens to determine cortisol levels (Webb ; Bechshøft ). However, the preparatory methods for the measurement of cortisol concentration differ from the methods used for the determination of stable isotope values, although both sets of measurements may be desired within the same study. For cortisol analysis, hair is commonly washed and extracted using methanol, and ground into a fine powder prior to extraction (Macbeth ). For stable isotope analysis (SIA), hair samples are typically washed in a 2:1 chloroform:methanol (v/v) mixture (e.g. Riofrío-Lazo and Páez-Rosas, 2015), and then cut into small pieces (Newsome ) or powdered (Darimont ; Riofrío-Lazo and Páez-Rosas, 2015). Because cortisol concentrations in blood and sweat are often far greater than those found in hair, failure to adequately wash the external surface of hair shafts contaminated with these substances will falsely elevate measurement of the internal hair shaft (medullar) cortisol concentration (e.g. Stalder and Kirschbaum, 2012). Similarly, the single wash that is often used for stable isotope analyses may not be sufficient to remove external contamination of the hair shaft with blood and/or perspiration (Russel ). On the other hand, if a solvent is too aggressive or if hair is immersed in a solvent for prolonged time, hair shaft cortisol may be lost if the solvent fully penetrates the hair cuticle (Eser ). Stable isotope values may also be affected by sample processing procedures (Font ). Carbon composition of surface oils and waxes can be significantly different from pure tissue composition and thus may influence stable isotope values. Therefore, in cleaning protocols for fixed tissue samples, volatile solvent mixtures are used to remove surface lipids (Wassenaar, 2008). An important question that remains is whether these preparatory procedures used for washing and/or grinding hair can be merged into a single protocol to increase the efficiency of sample preparations, and to reduce labor input and costs for research projects. It is unknown whether the preparation process for δ13C and δ15N analysis affects the reliable measurement of cortisol concentration in hair and vice versa. In this study, we compared commonly used, but different, wash and sample processing procedures to determine their effect on cortisol and stable isotope values measured in hair. We used hair samples collected from brown bears (Ursus arctos) as a model system. Here, we evaluate whether different preparatory procedures yielded consistent results. Specifically, we test the null hypotheses that (1) there are no differences in cortisol concentrations or (2) stable isotope ratios (δ13C and δ15N) between brown bear hair samples prepared using standard washing procedures in cortisol (methanol wash) or stable isotope ratio (2:1 chloroform:methanol rinse) measurement protocols, and that (3) there are no differences in stable isotope ratios between hair samples washed following the procedure in cortisol concentration measurement protocol and left intact or ground into powder before stable carbon and nitrogen ratio measurement.

Material and methods

Sample selection

We selected guard hair samples from 12 brown bear individuals from a larger pool of archived samples available, based on the criteria that the samples (1) were of sufficient quantity for further subdivision, (2) originated from the same body region, and (3) were free of visible external debris (e.g. blood, dirt). Samples were collected by plucking from the top of the shoulders (the hump) of bears captured in Sweden (1993–1996), or by cutting close to skin at the top of the neck for bears captured in Poland (2014–2015), and were stored dry in paper envelopes in the dark at room temperature. All samples were collected for ongoing long-term monitoring projects. Bear capture, immobilization and handling were carried out as described in Arnemo , and approved by the appropriate authorities (Sweden: Swedish Environmental Protection Agency, Stockholm, #NV-0758-14, and the Swedish Board of Agriculture, #31-11 102/12; Poland: General Directorate for Environmental Protection, Warsaw, #DOP-OZ.6401.08.2.2013.ls.1) and Ethical Committees (Sweden: Swedish Ethical Committee on Animal Research, Uppsala, #C18/15; Poland: I Local Ethical Committee in Krakow, #21/2013).

Sample preparation and procedures

Preparation of hair samples

We selected two subsamples (≥50 mg per subsample) that appeared similar with respect to length, color, and amount of external contamination from each of 12 samples of guard hair. The hair follicles were excised from all plucked hair samples to ensure that cortisol and stable isotope measurements were based on hair shafts only. Each subsample was further subdivided into two aliquots; one to be used for cortisol analysis and the other for SIA. Aliquots were washed with methanol (CORT procedure) or rinsed with 2:1 chloroform:methanol (SIA procedure) and processed accordingly for cortisol and stable isotope analyses as shown in Fig. 1. Hair aliquots for cortisol measurement were ground prior to analysis using a Retsch MM 301 Mixer Mill (Retsch Inc., Newtown, Pennsylvania, USA), whereas hair for stable isotope analyses were left as intact hair shafts, with the exception of one 6-mg aliquot that was ground and set aside. This aliquot was used to test the effect of the full CORT procedure (both methanol wash and grinding) on stable isotope ratios.
Figure 1:

The breakdown of a single hair sample (≥100 mg) collected from an individual bear into 50 mg subsamples, followed by smaller aliquots, and then processed accordingly following preparation procedures for cortisol (CORT) and stable isotope analyses (SIA; as indicated by arrows).

The breakdown of a single hair sample (≥100 mg) collected from an individual bear into 50 mg subsamples, followed by smaller aliquots, and then processed accordingly following preparation procedures for cortisol (CORT) and stable isotope analyses (SIA; as indicated by arrows).

CORT preparatory procedure

The CORT preparatory procedure is a standard protocol used for hair cortisol analysis that involves washing hair samples with HPLC-grade methanol in three cycles of 3 min each (Macbeth ). The washed hair is then air-dried thoroughly for 24 h at room temperature. Three aliquots per bear were prepared following this protocol (Fig. 1). A 10-mg aliquot was left intact and used for SIA, and the rest was ground into powder, according to the protocol for cortisol analysis. From the powdered hair, a 25-mg aliquot was used for cortisol analysis, and a 6-mg aliquot for SIA (Fig. 1).

SIA preparatory procedure

The SIA preparatory procedure is a standard protocol for SIA that involves rinsing hair samples with a 2:1 (v/v) mixture of chloroform:methanol (Hobson ). The hair is then removed from the solvent and air-dried thoroughly for 24 h at room temperature. Two 25 mg aliquots per bear were prepared following this protocol with one aliquot left intact and used for SIA and the other was ground and used for cortisol analysis (Fig. 1).

Hair cortisol analysis

For cortisol analysis, we used 24 aliquots of hair (2 × 25 mg aliquots per bear) that were washed using either the SIA (chloroform:methanol rinse; one per bear) or CORT (methanol wash; one per bear) procedure. Aliquots of 25 mg were combined with 0.5 ml of HPLC grade methanol for 24 h on a slow rotator to extract cortisol. Subsequent procedures and cortisol analysis were conducted in accordance with the protocol described by Macbeth . The extract was centrifuged for 15 min at 2150g, the methanol was removed, evaporated until dryness under nitrogen gas (38°C), and reconstituted in phosphate buffer (0.2 ml). Cortisol was quantified using a commercially available enzyme immunoassay (EIA) kit previously validated for use with brown bear hair (Oxford Biomedical, Lansing, Michigan, USA; Macbeth ). Extracts were run in duplicates with intra-assay coefficient of variability (CV) of 4.6%, calculated from the individual CVs for all of the duplicates.

Stable isotope analysis

For SIA, we used 24 aliquots of intact hair that were washed using either the SIA (10 mg per bear) or CORT (25 mg per bear) procedures, and additionally 12 aliquots (6 mg per bear) prepared according to the CORT procedure and ground into powder. Subsequent δ13C and δ15N analyses were conducted in accordance with the protocol described in Hobson . Portions (1 mg) of each aliquot were weighed into tin capsules and combusted in a Carlo Erba elemental analyzer (Milan, Italy) interfaced with a Europa 20:20 continuous flow mass spectrometer. Results are reported in delta (δ) notation as parts per thousand (‰) deviation from international standards, Vienna Pee dee Belemnite (VPDB) for δ13C and atmospheric (AIR) for δ15N. Based on within-run (n = 12) replicate measurements of an egg albumen standard, we estimated measurement error to be ±0.1‰ and ±0.2‰ for δ13C and δ15N, respectively.

Statistical analysis

Data distribution was checked using the Shapiro–Wilk test and the assumption of normality was not fulfilled. Therefore, a non-parametric pairwise test (Wilcoxon signed-rank test) was used to compare the cortisol concentration in hair from the same individual (n = 12; Supplementary table 1) prepared by either CORT or SIA washing protocols (Fig. 1). We used the same non-parametric test to compare δ13C and δ15N values obtained by preparing hair samples of the same individual (n = 12; Supplementary table 1) by the SIA washing and the CORT washing procedure, as well as the CORT washing procedure followed by grinding of the hair sample (referred to as CORT/grind in the results). We did not apply sequential Bonferroni corrections in the pairwise tests, because this method is unnecessarily conservative (e.g. Narum, 2006). In addition, all of our test results had P-values far above 0.05 (see Results), which we considered as statistically significant. We present results for all analyses as median (range) as well as mean ± SD. All tests were carried out in the statistical software R 3.3.2 (R Core Team, 2016). A comparison of δ13C, δ15N (‰) and cortisol concentrations (pg/mg) values measured in hair samples of all individuals, descriptive statistics, concordance correlation coefficients, as well as further analyses comparing the effect of cleaning procedures on carbon and nitrogen fraction (%) and carbon:nitrogen ratio are included as Supplementary material.

Results

We found no significant difference (V = 36, P = 0.884) between the pairwise cortisol concentrations of hair washed using the SIA procedure (1.19 (0.66, 6.69); 1.62 ± 1.63 pg/mg) compared to the CORT procedure (1.29 (0.82, 4.60); 1.50 ± 1.03 pg/mg; Fig. 2).
Figure 2:

Box plot of median cortisol level (in picograms per milligram) in hair samples of 12 brown bear individuals across different washing protocols used in preparatory stage. The cortisol concentrations in aliquots rinsed with 2:1 chloroform:methanol mixture, as typical in SIA, were not significantly different from aliquots washed in methanol, as typically used in cortisol concentration analysis (CORT), by the Wilcoxon test (V = 36, P = 0.884).

Box plot of median cortisol level (in picograms per milligram) in hair samples of 12 brown bear individuals across different washing protocols used in preparatory stage. The cortisol concentrations in aliquots rinsed with 2:1 chloroform:methanol mixture, as typical in SIA, were not significantly different from aliquots washed in methanol, as typically used in cortisol concentration analysis (CORT), by the Wilcoxon test (V = 36, P = 0.884). We found no significant differences in δ15N values between CORT and SIA procedures (V = 43, P = 0.791; CORT: 4.8‰ (1.9‰, 6.6‰), 4.9 ± 1.3‰; SIA: 4.7‰ (2.9‰, 6.4‰), 4.9 ± 1.0‰), CORT and CORT/grind procedures (V = 30, P = 0.505; CORT/grind: 4.8‰ (2.6‰, 6.4‰), 5.0 ± 1.2‰), or CORT/grind and SIA procedures (V = 42, P = 0.850; Fig. 3). Similarly, we found also no significant differences in δ13C values between CORT and SIA procedures (V = 45.5, P = 0.638; CORT: −22.3‰ (−23.9‰, −14.7‰), −21.4 ± 2.4‰; SIA: −22.3‰ (−23.8‰, −14.8‰), −21.5 ± 2.5‰), CORT and CORT/grind procedures (V = 29, P = 0.756; CORT/grind: −22.3‰ (−23.9‰, −15.8‰), −21.3 ± 2.2‰) and CORT/grind and SIA procedures (V = 39, P = 0.999; Fig. 3).
Figure 3:

Box plot of median stable nitrogen δ15N and carbon δ13C values in hair samples of 12 brown bear individuals across different washing protocols used in preparatory stage. The δ15N values in aliquots rinsed with 2:1 chloroform:methanol mixture (SIA), were not significantly different from aliquots washed in methanol (CORT), and from the aliquots washed with methanol and ground into powder (CORT/grind), by the Wilcoxon test (V = 43, P = 0.791 and V = 42, P = 0.850, respectively). No significant difference was found in δ15N values between CORT and CORT/grind procedures (V = 30, P = 0.505). The δ13C values were not significantly different between CORT and SIA procedures (V = 45.5, P = 0.638), CORT and CORT/grind procedures (V = 29, P = 0.756) and CORT/grind and SIA procedures (V = 39, P = 0.999).

Box plot of median stable nitrogen δ15N and carbon δ13C values in hair samples of 12 brown bear individuals across different washing protocols used in preparatory stage. The δ15N values in aliquots rinsed with 2:1 chloroform:methanol mixture (SIA), were not significantly different from aliquots washed in methanol (CORT), and from the aliquots washed with methanol and ground into powder (CORT/grind), by the Wilcoxon test (V = 43, P = 0.791 and V = 42, P = 0.850, respectively). No significant difference was found in δ15N values between CORT and CORT/grind procedures (V = 30, P = 0.505). The δ13C values were not significantly different between CORT and SIA procedures (V = 45.5, P = 0.638), CORT and CORT/grind procedures (V = 29, P = 0.756) and CORT/grind and SIA procedures (V = 39, P = 0.999).

Discussion and conclusions

We failed to reject our null hypotheses and found no differences in cortisol concentrations or in stable isotope values for brown bear hair washed with methanol, used in sample preparation for cortisol measurement, compared with the 2:1 chloroform:methanol rinse used for stable isotope preparations. This suggests that the chloroform:methanol rinse was adequate at cleaning the external surface of hair shafts without extracting endogenous cortisol. An important caution, however, is that hair samples used in this study were visibly clean, i.e. no obvious blood or dirt was observed on the hair shafts. It is possible that a single rinse with chloroform:methanol may not have been adequate to remove external contamination had we used hair that was significantly contaminated with blood, sebum, or dirt. There were also no differences found in stable isotope values for ground hair vs. intact hair, indicating that the standard full protocol for hair cortisol sample preparation was sufficient for obtaining reliable stable isotope ratio data. As shown, initial preparatory procedures used for either hair cortisol or stable isotope determinations provide reliable and reproducible measurements. Furthermore, while the methods share a common basis (washing and/or grinding), the overall time in the lab would be reduced with samples prepared and processed in one laboratory and aliquots shared for the final stage of the other analysis. Shorter overall processing and thus laboratory time decrease the labor input, which also lowers personnel costs. With respect to analytical costs, we found a wide range of prices in North America and Europe (price offers available online as for mid-August 2016, converted to US$). The cost of SIA per unprepared sample ranged from $32.50 to $78 (excluding shipping), and decreased to $8.50–20 (excluding shipping) if the sample is sent fully prepared (washed, dried, weighed, cut or ground or homogenized, and packed in tin capsules). For hair cortisol analysis, the costs ranged from $35 to $55 per sample. The analytical costs can therefore decrease by about 50% when samples are sent fully prepared (washed and powdered) for final analytical procedures. Though inter-laboratory calibrations are routine, the exercise to unify preparatory methods presented in this study is worth considering, as it could also lower the ecological footprint of the projects. The additional value comes in using archived or stored materials prepared by one protocol for the other, and pursuing hypotheses related to both assays within same study. Combining cortisol and stable isotope measurements can advance questions relevant to conservation physiology (e.g. relations of stress and nutritional status; Kempster ; Deschner ), and identify mechanisms underlying the patterns and trends observed in wildlife populations (e.g. Albano, 2012). In conclusion, the findings of this study suggest that the common washing procedures and preprocessing steps used for the analyses of cortisol and stable δ13C and δ15N in animal hair do not lead to significant discrepancies in the final results. In other words, one procedure, using either a 2:1 chloroform:methanol mixture or methanol alone, can be used to wash hair in advance of determining cortisol, δ13C, and δ15N values. Further, stable isotope values can be measured using either intact hair shafts or ground hair, the latter being required for cortisol analysis. This greatly reduces the costs and labor of sample preparation, and expands the set of samples available for forensic, physiological and ecological studies.

Author contributions

Conceived and designed the experiment: A.S., K.A.H., D.M.J., N.S., A.Z. Performed the experiments: A.S., L.K., K.A.H., C.G. Analyzed the data: A.S., M.C., A.Z. Contributed reagents/materials/analysis tools: D.M.J., K.A.H., M.C. Wrote the paper: A.S., M.C., D.M.J., A.Z., K.A.H., N.S., L.K., C.G. Click here for additional data file.
  17 in total

Review 1.  Conservation physiology.

Authors:  Martin Wikelski; Steven J Cooke
Journal:  Trends Ecol Evol       Date:  2005-11-11       Impact factor: 17.712

2.  Analysis of endogenous cortisol concentrations in the hair of rhesus macaques.

Authors:  Matthew D Davenport; Stefan Tiefenbacher; Corrine K Lutz; Melinda A Novak; Jerrold S Meyer
Journal:  Gen Comp Endocrinol       Date:  2006-02-17       Impact factor: 2.822

3.  Tracking stress: localisation, deposition and stability of corticosterone in feathers.

Authors:  Gary R Bortolotti; Tracy Marchant; Julio Blas; Sonia Cabezas
Journal:  J Exp Biol       Date:  2009-05       Impact factor: 3.312

4.  Influence of sample preparation on analytical results: drug analysis [GC/MS] on hair snippets versus hair powder using various extraction methods.

Authors:  H P Eser; L Pötsch; G Skopp; M R Moeller
Journal:  Forensic Sci Int       Date:  1997-01-17       Impact factor: 2.395

Review 5.  Analysis of cortisol in hair--state of the art and future directions.

Authors:  Tobias Stalder; Clemens Kirschbaum
Journal:  Brain Behav Immun       Date:  2012-02-15       Impact factor: 7.217

6.  Identification of energy consumption and nutritional stress by isotopic and elemental analysis of urine in bonobos (Pan paniscus).

Authors:  Tobias Deschner; Benjamin T Fuller; Vicky M Oelze; Christophe Boesch; Jean-Jacques Hublin; Roger Mundry; Michael P Richards; Sylvia Ortmann; Gottfried Hohmann
Journal:  Rapid Commun Mass Spectrom       Date:  2012-01-15       Impact factor: 2.419

7.  Isotopic segregation between sympatric seabird species increases with nutritional stress.

Authors:  Christopher P Barger; Alexander S Kitaysky
Journal:  Biol Lett       Date:  2011-12-14       Impact factor: 3.703

8.  Measuring environmental stress in East Greenland polar bears, 1892-1927 and 1988-2009: what does hair cortisol tell us?

Authors:  T Ø Bechshøft; F F Rigét; C Sonne; R J Letcher; D C G Muir; M A Novak; E Henchey; J S Meyer; I Eulaers; V L B Jaspers; M Eens; A Covaci; R Dietz
Journal:  Environ Int       Date:  2012-05-07       Impact factor: 9.621

9.  Do stable isotopes reflect nutritional stress? Results from a laboratory experiment on song sparrows.

Authors:  Bethany Kempster; Liana Zanette; Fred J Longstaffe; Scott A MacDougall-Shackleton; John C Wingfield; Michael Clinchy
Journal:  Oecologia       Date:  2006-11-11       Impact factor: 3.298

10.  Stress and reproductive hormones in grizzly bears reflect nutritional benefits and social consequences of a salmon foraging niche.

Authors:  Heather M Bryan; Chris T Darimont; Paul C Paquet; Katherine E Wynne-Edwards; Judit E G Smits
Journal:  PLoS One       Date:  2013-11-27       Impact factor: 3.240

View more
  3 in total

1.  Whiskers as a novel tissue for tracking reproductive and stress-related hormones in North Pacific otariid pinnipeds.

Authors:  Mandy J Keogh; Patrick Charapata; Brian S Fadely; Tonya Zeppelin; Lorrie Rea; Jason N Waite; Vladimir Burkanov; Chris Marshall; Aubree Jones; Caitlin Sprowls; Matthew J Wooller
Journal:  Conserv Physiol       Date:  2021-01-12       Impact factor: 3.079

2.  Do follicles matter? Testing the effect of follicles on hair cortisol levels.

Authors:  Agnieszka Sergiel; Marc Cattet; Luciene Kapronczai; David M Janz; Nuria Selva; Kamil A Bartoń; Jon E Swenson; Andreas Zedrosser
Journal:  Conserv Physiol       Date:  2020-02-01       Impact factor: 3.079

3.  Can concentrations of steroid hormones in brown bear hair reveal age class?

Authors:  Marc Cattet; Gordon B Stenhouse; John Boulanger; David M Janz; Luciene Kapronczai; Jon E Swenson; Andreas Zedrosser
Journal:  Conserv Physiol       Date:  2018-01-29       Impact factor: 3.079

  3 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.