Literature DB >> 32944567

Detection and quantification of Campylobacter in foods: New analytic approaches to detect and quantify Campylobacter spp. in food samples.

Maria Francesca Peruzy1, Yolande Thérèse Rose Proroga2, Federico Capuano2, Federica Corrado2, Serena Santonicola3, Dario De Medici4, Elisabetta Delibato4, Nicoletta Murru1.   

Abstract

The aim of the present study was to develop rapid qualitative and quantitative methods based on the use of Real-Time PCR and Droplet Digital PCR (ddPCR), in order to have reliable techniques to detect and quantify Campylobacter spp. in food samples. The gene 16S-rRNA was used as specific target for Campylobacter spp. Real- Time PCR evaluation assay and a not competitive internal control was ushered in it. To investigate the selectivity of the method, 26 Campylobacter strains and 40 non-Campylobacter strains were tested and in order to verify the application of Real- Time PCR method, 5 pork meat samples were experimentally inoculated with a Campylobacter jejuni strain. Subsequently, dilutions with a bacterial load of Campylobacter jejuni within 10-106 CFU/mL were chosen for the optimization of the ddPCR assay. Lastly, a total of 54 naturally contaminated foods samples were analyzed through molecular (Real-Time PCR and ddPCR) and traditional methods. The Real-Time PCR protocol demonstrated to amplify only the Campylobacter spp. strains and when Campylobacter jejuni was experimentally inoculated in meat samples the pathogen was always detected. The ddPCRs assay allowed to quantify a level of contamination of 10 CFU/mL, but it was unable to quantify levels of 105 - 106 CFU/mL. Lastly, Campylobacter spp. was never detected in the 54 samples tested. In conclusion, the novel analytic approach proposed, based on an initial screening of the samples with Real-Time PCR and then on quantification of Campylobacter spp. with a ddPCR on those positive, represents a quick monitoring tool and, if used correctly, it would allow the implementation of food safety. ©Copyright: the Author(s).

Entities:  

Keywords:  Campylobacter spp; Droplet digital PCR; Real-Time PCR

Year:  2020        PMID: 32944567      PMCID: PMC7477723          DOI: 10.4081/ijfs.2020.8591

Source DB:  PubMed          Journal:  Ital J Food Saf        ISSN: 2239-7132


Introduction

Since 2005, Campylobacter (C) is the most commonly reported gastrointestinal bacterial pathogen in humans in the European Union (EU) (EFSA, 2018). In particular, according to the European Food Safety Authority (EFSA), in 2017 there were 246,158 confirmed cases of human campylobacteriosis in the EU (EFSA, 2018). Globally, the incidence varies among countries, but the current incidence of campylobacteriosis is probably underestimated because of difficulties with diagnosis, differences in reporting systems and in surveillance (Hansson et al., 2018). Among Campylobacter species, C. jejuni and C. coli have the highest rate of foodborne-related clinical campylobacteriosis (Ricke et al., 2019). Humans can acquire the infection through direct contact with infected animal or through the consumption of contaminated foods. The gastrointestinal illness usually is self-limiting, but in a small number of cases complications such as reactive arthritis and neurological disorders may occur (Fabiani et al., 2019). The most commonly identified sources of human infection are represented by broiler meat and milk, but also dairy product and red meat have been responsible for outbreaks in 2017 (EFSA, 2018). Conventional methods for the qualitative and quantitative detection of Campylobacter involve the use of selective media and elevated incubation temperature (42°C) in microaerophilic atmosphere (ISO 10272-1 and ISO 10272-2) but although inexpensive and simple, these methods are long and timeconsuming, taking up to 6 days (Erdosi et al., 2018; Fabiani et al., 2019). The development of rapid methods is important to conduct a proper control of the contaminated products, to identify quickly outbreak sources and to prevent the spread of illness (Wang and Salazar, 2017.). Over the past few decades, to reduce the detection time from days to hours, molecular methods, such as polymerase chain reaction (PCR), have been introduced in daily laboratory routines to detect pathogenic bacteria in different foods (Peruzy et al., 2019; Ricke et al., 2019). According to the EC regulation 1495/2017 a quantitative determination of Campylobacter is required, thus rapid, reliable and reproducible methods for the quantification of this pathogen in foods are needed (Yu et al., 2019). Recently, Droplet digital PCR (ddPCR) emerged as a useful tool for precise quantification of nucleic acids in different samples and it allows the detection of target DNA in complex matrices in presence of abundant sequences. In this technique, target DNA molecules are distributed across thousands of small droplets. Each droplet constitutes a separate reaction compartment and contains 0 or 1 (or more) template copies (Hindson et al., 2011). After amplification through a conventional PCR, a Poisson distribution is used to extrapolate the fraction of positive end-point reactions. It has been firstly proposed for clinical and environmental studies but more recently it has been introduced also in food microbiology analysis (Gobert et al., 2018). In fact, this technology has been recently used to quantify Bacillus cereus and Salmonella typhimurium in milk (Porcellato et al., 2016; Wang et al., 2018) and noroviruses in oysters (Persson et al., 2018). Moreover, compared with Real-Time PCR, ddPCR is more rapid since does not require standard curve for the analysis and it is less sensitive to inhibitions (Gobert et al., 2018; Wang et al., 2018). The aim of the present study was the development of rapid qualitative and quantitative methods based on the use of Real-Time PCR and Droplet Digital PCR (ddPCR), respectively, in order to have sensitive and reliable techniques to detect and quantify Campylobacter spp. in food samples. Moreover, the occurrence of Campylobacter spp. in different pork meat samples and dairy products was evaluated.

Materials and Methods

Campylobacter strains and other organism

Twenty-six Campylobacter strains, belonging to Campylobacter jejuni, Campylobacter coli, Campylobacter lari, Campylobacter upsaliensis and Campylobacter fetus fetus and 40 non- Campylobacter microorganisms (Table 1) isolated from foods and human cases were included in this study. The strains were supplied from the collection of the Pathogenic Enterobacteria Unit of Istituto Superiore di Sanità, Rome, Italy. Bacterial strains were cultured in Tryptic Soy Agar (TSA; CM0131, OXOID, Basingstoke,UK) for 48-h at 41,5°C in microaerophilic conditions.
Table 1.

Number and species of the Campylobacter strains and other organisms included in the study and Real-Time PCR results.

Bacterial strainsN. of strainsReal-Time PCR
Campylobacter strains
    Campylobacter jejuni10+
    Campylobacter coli10+
    Campylobacter lari2+
    Campylobacter upsaliensis2+
    Campylobacter fetus fetus2+
Other microorganisms
    Aeromonas idrofila1-
    Bacillus subtilis1-
    Citrobacter freundii1-
    Enterobacter aerogene1-
    Enterobacter cloacae1-
    Escherichia coli1-
    Listeria innocua1-
    Listeria monocytogenes2-
    Proteus hauseri1-
    Pseudomonas aerugionosa3-
    Shighella boydii1-
    Shighella flexeneri1-
    Shighella sonnei1-
    Staphyloccocus aureus6-
    Staphyloccocus xylosus1-
    Yersinia enterocolitica5-
    Yersinia intermedia1-
    Yersinia kristensenii1-
    S. Enteritidis2-
    S. Typhimurium2-
    S. Napoli2-
    S. Thompson2-
    S. Veneziana2-

DNA extraction

DNA was extracted using the Chelex- 100-resin method (Bio-Rad, Hercules, CA, USA) whereby one colony of each strain was suspended in 300 μL of 6% Chelex 100 by vortexing, and incubated for 20 min at 56°C and again for 8 min at 100°C. The suspension was immediately chilled on ice for 1 min and centrifuged for 5 min at 10,000 ×g at 4°C (Peruzy et al., 2019).

qPCR

The gene 16S-rRNA was used as specific target for Campylobacter spp. Real- Time PCR evaluation assay. In particular, 5 μL of DNA extracted were used as template and added to a PCR Mix (Fabiani et al., 2019). The mix contained: 1X Quantitect multiplex PCR NoROX, 440 nM forward primer 5′ CTGCTTAACACAAGTTGAGTAGG 3, 480 nM reverse primer 5′ TTCCTTAGGTACCGTCAGAA 3′, 200 nM 16S-rRNA fluorescent probe 5′ FAMTGTCATCCTCCACGCGGCGTTGCTGCBHQ1 3′and water for molecular biology. Moreover, in order to determine the inhibition in the biological samples and to identify false negative results a not competitive internal control (IAC), primer pUC 18-F (5‘–TGT CGT GCC AGC TGC ATT A-3‘), primer pUC 18-R (5‘–GAG CGA GGA AGC GGA AGA g–3‘) and the probe Tm-pUC18 (5‘–HEX– AAT CGG CCA ACG CGC GG –BHQ1–3) were added. In particular, the optimal IAC concentration was 1000 copies per reaction. The reaction was run online at 95°C for 15 min, followed by 40 cycles at 95°C for 30 s and 58°C for 60 s and an extension phase of 1 cycle at 95°C for 60 s, 60°C for 60 s, and 95°C for 60 s (ramp time, 19.59 min). Amplification reactions were performed with Bio-Rad CFX96 platform, using a 96- well PCR multiplate, from Bio-Rad. To investigate whether the primer set and the amplification protocol amplify DNA from Campylobacter spp. (inclusivity/exclusivity), the DNA extracted with the Chelex-100-resin method of the twenty-six Campylobacter strains and of the 40 non-Campylobacter strains (Table 1) were tested as described above.

Inoculation of the food matrix

In order to verify the application of the above described method, five pork meat samples purchased at different supermarkets, were experimentally inoculated with Campylobacter jejuni strain (ATCC33291). In particular, Campylobacter jejuni was activated twice in Tryptone Soy Broth (TSB; CM0129, OXOID, Basingstoke, UK) and incubated at 41,5 °C for 48-h in a microaerophilic conditions and standardized by count on tryptone soy agar. Each sample was divided into three aliquots, one was used as negative control and the other two were inoculated with a bacterial load within 1-10 and 10-102 CFU/mL, respectively. In particular, 25 grams of each sample were placed in a sterile stomacher bag and homogenized for 3 minutes at 230 rpm using a peristaltic homogenizer (BagMixer®400 P, Interscience, Saint Nom, France) after the addition of 225 mL (1:10, W/W) of sterilized Bolton Selective Enrichment Broth (BSEB), prepared by adding bolton broth selective supplement and laked horse blood in bolton Broth, following technical data sheet. After inoculation, all samples were shaken for 1 min to optimize the distribution of the inoculum and then they were incubated at 37°C for 4-h and 41,5°C for 44-h in microaerophilic conditions. After the incubation, two mL of each incubated homogenate were subjected to the DNA extraction phase using the procedure above described. Subsequently the Real- Time PCR analysis was performed as above described. Moreover, in order to compare the results of the molecular method with the microbiological reference method (EN/ISO 10272-1:2017), after 4 and 44 hours of incubation a sterile loop was used to transfer the homogenate broth on the surface of Modified Charcoal Cefoperazone Desoxycholate Agar (mCCD, CM0739 and SR0155, OXOID, Basingstoke, UK) and Campylobacter agar base (karmali, CM0935, OXOID, Basingstoke, UK) with the Campylobacter selective supplement (karmali supplement, SR0167, OXOID, Basingstoke, UK). Subsequently all the plates were incubated in a microaerophilic conditions at 41,5°C for 44- h. Afterward, suspected colonies were picked, subcultured on Columbia blood agar (CM0331, OXOID, Basingstoke, UK) and incubated in a microaerophilic conditions at 41,5°C for 24/48-h for the phenotypic and biochemical confirmation. Number and species of the Campylobacter strains and other organisms included in the study and Real-Time PCR results.

Digital droplet PCR reaction

Ten-fold serial dilutions of the culture of Campylobacter jejuni were prepared in TSB and the concentration of cells was measured using McFarland turbidity standards. Dilutions with a bacterial load within 10-106 CFU/mL were chosen for the optimization of the Droplet Digital PCR (ddPCR). The same DNA extraction protocol, primer set and probe used for the above described Real- Time PCR assay were applied for the optimization of the ddPCR assay. In particular, ddPCR reaction were performed using 1×ddPCR Supermix for probe (Bio- Rad, Pleasanton, CA) and the final volume of the reaction was 22 μL. Twenty microliter of the ddPCR reaction were used to generate the droplet mix in an 8-well cartridge using the QX100 droplet generator (Bio-Rad). The emulsion (40 μL) was then transferred to a 96-well plate and amplified using the same condition of the previously described Real- Time PCR. The PCR products were denatured at 98°C for 10 min and kept at 4°C until the droplets were read. Ramp rate used in the droplet PCR was 2°C s−1. The 96- well plate was then transferred to the QX100 droplet reader (Bio-Rad) and data acquisition and analysis was performed using QuantaSoft software ver 1.7 (Bio- Rad). The fluorescence amplitude threshold used for the discrimination of the positive and negative droplets in QuantaSoft software was set between 2000 and 2200. The concentration values were calculated by QuantaSoft software (in copies μL−1) and multiplied by 22 (the initial PCR volume) to obtain the absolute number of copies added to the PCR reaction.

Campylobacter spp. detection from naturally contaminated samples

In order to evaluate the presence of Campylobacter spp. in naturally contaminated food samples, a total of 28 pork meat samples and 26 dairy products were collected in the Campania region in southern Italy (Table 2) and analyzed through molecular methods (Real-Time PCR and ddPCR) as described above and traditional microbiological methods (EN/ISO 10272-2:2017). In particular, for the traditional microbiological method, from each sample 10 grams were placed in a sterile stomacher bag and homogenized for 3 minutes after the addition of 90 mL (1:10, W/W) of sterilized Buffered Peptone Water (BPW; CM0509, OXOID, Basingstoke, UK). All samples were incubated in a microaerophilic conditions at 41,5°C for 1- h and 48-h. After 1 and after 48-h, 10-fold serial dilutions of the incubated homogenate were prepared in BPW and 0,1 mL of the appropriate 10-fold serial dilutions were spread on the surface of Modified Charcoal mCCD. Subsequently all the plates were incubated in a microaerophilic condition at 41,5°C for 24/48-h. After viable counts, five suspected colonies were picked up from each agar plate, subcultured on Columbia blood agar (CM0331, OXOID, Basingstoke, UK) and incubated in a microaerophilic conditions at 41,5°C for 24/48-h. Afterward, phenotypic and biochemical test were used for the confirmation of the suspected colonies.
Table 2.

Meat samples and dairy products (type and number) collected in Campania region and analyzed through Real-Time PCR, ddPCR and EN/ISO 10272 method for the evaluation of the presence of Campylobacter spp.

MatricesDescription Meat samplesN. of samples
Capocollodry-cured pork neck4
Filone stagionatodry-cured pork meat4
Guanciale stagionato al pepedry-cured pork jowl with pepper4
Pancetta al pepe e finocchiettodry-cured rolled pork belly with pepper and fennel4
Salamedry-cured pork meat4
Salsicciafresh pork meat8
Dairy samples
Cow's milk
    Caciocavallosemi-hard pasta filata cheese6
    Mozzarellasemisoft pasta filata cheese2
    Mozzarella nella mortellasemisoft pasta filata cheese wrapped in a myrtus leaf2
    Ricottaunripened acid-heat coagulated dairy product2
    Scamorzasemi-hard pasta filata cheese6
Goat's milk
    Caciotta frescafresh semi-soft cheese1
    Caciotta stagionatacured semi-soft cheese1
Buffalo's milk
    Mozzarellasemisoft pasta filata cheese2
    Ricottaunripened acid-heat coagulated dairy product2
    Scamorzasemi-hard pasta filata cheese2

Results

Real-time PCR assay and detection of Campylobacter jejuni in experimentally contaminated foods

The optimal amount of internal amplification control to allow target detection, even when present with low copy numbers, was estimated as 1000 copies/reaction yielding a 28-30 Ct value. The Real-Time PCR protocol demonstrated to amplify only the Campylobacter spp. strains (Table 1) and not the DNA of the 40 non-Campylobacter strains tested showing a 100% of inclusivity and exclusivity. Meat samples and dairy products (type and number) collected in Campania region and analyzed through Real-Time PCR, ddPCR and EN/ISO 10272 method for the evaluation of the presence of Campylobacter spp. Moreover, when Campylobacter jejuni was experimentally inoculated in pork meat samples to verify the application of the molecular method, the pathogen was always detected, both when the lowest (1-10 CFU/mL) and the highest (10-102 CFU/mL) concentration were used. Furthermore, the molecular method evaluated showed a concordance rate of 100% with the reference microbiological method (EN/ISO 10272-1:2017). Evaluation of the ddPCR sensitivity, was performed using serial dilution of Campylobacter jejuni DNA. The lowest level of detection was 10 CFU/mL whereas ddPCR was unable to quantify high concentration of the pathogen (105-106 CFU/mL) (Figure 1).
Figure 1.

Droplet distribution obtained from the serial dilution of Campylobacter jejuni genomic DNA.

Detection of Campylobacter in naturally contaminated samples

Campylobacter spp. was never detected in the 54 samples tested, both with molecular (Real-Time PCR and ddPCR) and traditional microbiological methods (EN/ISO 10272-2:2017).

Discussion

Conventional methods based on the culture on microbiological media and identification through the biochemical test are still used in the daily laboratory routine for the isolation and/or enumeration of Campylobacter species in food samples. Although these methods are sensitive and inexpensive, they are long and timeconsuming. Over the past few decades, molecular methods, as the nucleic acidbased technology, have been championed as a promising alternative for the fast and accurate detection and quantification of microorganisms from different foods (Peruzy et al., 2017). In the present work, two molecular platforms (Real-Time PCR and ddPCR), based on nucleic acid amplification, have been developed in order to determine and quantify Campylobacter spp. from food samples. Real-Time PCR, performed on a set of different Campylobacter and non- Campylobacter species, discriminated successfully the target DNA from the other microorganisms (100% inclusivity and exclusivity). Moreover, when C. jejuni was experimentally contaminated in pork samples, this platform resulted to be effective at detecting the pathogen even when present at low concentration (1-10 CFU/mL). Furthermore, since the presence of inhibitor in food matrices may lead to false-negative results, a non-competitive IAC has been away used in the Real-Time PCR assay, in order to have a tool to check on inhibition effect in PCR. Even though the principal purpose of the Real-Time PCR technology is the quantification of target bacteria, its use in qualitative analysis, as performed in the present work, should not be underestimated. In fact, compared with traditional PCR (endpoint PCR), Real-Time PCR is more precise, quicker, as does not require post-PCR processing (es. preparation of gels), and moreover, it is safer, as the ethidium bromide is not needed (Ahmed et al., 2009). Thus, the PCR developed may be a useful tool for initial screening of food samples and only those samples positive for Campylobacter would then be subjected to further analysis. However, according to EU Regulation 2017/1495, a quantitative determination of Campylobacter in broiler carcasses is required. Thus, a rapid molecular method for the precise quantification of this pathogen is also needed. Over the past years, quantitative Real-Time PCR has been successfully used to quantify target DNA in different samples. However, the major limitations of the Real- Time PCR are the request of a standard curve for the calculation of the target concentration and the sensitivity of the amplification to inhibitors (Gobert et al., 2018). Recently, droplet digital PCR has been introduced for the absolute quantification of a DNA target; it does not require a standard curve and, moreover, the PCR reaction partitioning in the ddPCR assay reduces the exposure to PCR inhibitors in the droplet (Porcellato et al., 2016). In the present work, ddPCR allowed the quantification of low numbers of target molecules (10 CFU/mL) but it was unable to quantify high concentrations of Campylobacter jejuni (≥ 105 CFU/ml), probably because the ddPCR droplets became completely saturated (Tang et al., 2016). The two developed molecular platforms (Real-Time PCR and ddPCR) together with the microbiological reference method (EN/ISO 10272:2017) were used in order to evaluate the presence of Campylobacter spp. in different food samples. Although in the present work the number of the samples were limited, our results are in disagreement with the ones of previously published researches where the detection/isolation of Campylobacter in different food samples was higher (El-zamkan et al., 2016;Karki, et al., 2018; Zhao et al., 2001). Droplet distribution obtained from the serial dilution of Campylobacter jejuni genomic DNA.

Conclusions

The novel analytic approach proposed, based on an initial screening of the samples with Real-Time PCR and then on quantification of Campylobacter spp. with a ddPCR on those samples resulted positive, represent a quick monitoring tool and, if used, would allow the implementation of food safety. Moreover, by using these methods could be possible to gain more data for a timely and accurate risk assessment. Furthermore, this proposed strategy, developed in the present study for Campylobacter, could be used, also, for a rapid detection and enumeration of other microorganisms important both in food and in clinical microbiology.
  19 in total

1.  Comparison between digital PCR and real-time PCR in detection of Salmonella typhimurium in milk.

Authors:  Meng Wang; Junjie Yang; Zhongtao Gai; Shengnan Huo; Jianhua Zhu; Jun Li; Ranran Wang; Sheng Xing; Guosheng Shi; Feng Shi; Lei Zhang
Journal:  Int J Food Microbiol       Date:  2017-12-13       Impact factor: 5.277

2.  Prevalence of Campylobacter spp., Escherichia coli, and Salmonella serovars in retail chicken, turkey, pork, and beef from the Greater Washington, D.C., area.

Authors:  C Zhao; B Ge; J De Villena; R Sudler; E Yeh; S Zhao; D G White; D Wagner; J Meng
Journal:  Appl Environ Microbiol       Date:  2001-12       Impact factor: 4.792

3.  The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2017.

Authors: 
Journal:  EFSA J       Date:  2018-12-12

4.  Droplet digital PCR improves absolute quantification of viable lactic acid bacteria in faecal samples.

Authors:  Guillaume Gobert; Aurélie Cotillard; Candice Fourmestraux; Laurence Pruvost; Jean Miguet; Mickaël Boyer
Journal:  J Microbiol Methods       Date:  2018-03-14       Impact factor: 2.363

5.  Rapid in-house detection method of Campylobacter spp. from food by redox potential monitoring combined with real-time PCR.

Authors:  Orsolya Erdősi; Katalin Szakmár; Zsuzsanna Szili; Géza Szita; Sándor Bernáth; József Sövényi; Péter Laczay
Journal:  Acta Vet Hung       Date:  2018-03       Impact factor: 0.955

6.  Comparison between RT droplet digital PCR and RT real-time PCR for quantification of noroviruses in oysters.

Authors:  Sofia Persson; Ronnie Eriksson; James Lowther; Patrik Ellström; Magnus Simonsson
Journal:  Int J Food Microbiol       Date:  2018-07-02       Impact factor: 5.277

7.  Detection and quantification of Bacillus cereus group in milk by droplet digital PCR.

Authors:  Davide Porcellato; Judith Narvhus; Siv Borghild Skeie
Journal:  J Microbiol Methods       Date:  2016-05-17       Impact factor: 2.363

8.  High-throughput droplet digital PCR system for absolute quantitation of DNA copy number.

Authors:  Benjamin J Hindson; Kevin D Ness; Donald A Masquelier; Phillip Belgrader; Nicholas J Heredia; Anthony J Makarewicz; Isaac J Bright; Michael Y Lucero; Amy L Hiddessen; Tina C Legler; Tyler K Kitano; Michael R Hodel; Jonathan F Petersen; Paul W Wyatt; Erin R Steenblock; Pallavi H Shah; Luc J Bousse; Camille B Troup; Jeffrey C Mellen; Dean K Wittmann; Nicholas G Erndt; Thomas H Cauley; Ryan T Koehler; Austin P So; Simant Dube; Klint A Rose; Luz Montesclaros; Shenglong Wang; David P Stumbo; Shawn P Hodges; Steven Romine; Fred P Milanovich; Helen E White; John F Regan; George A Karlin-Neumann; Christopher M Hindson; Serge Saxonov; Bill W Colston
Journal:  Anal Chem       Date:  2011-10-28       Impact factor: 6.986

9.  Prevalence of Campylobacter jejuni and Campylobacter coli in raw milk and some dairy products.

Authors:  Mona A El-Zamkan; Karima G Abdel Hameed
Journal:  Vet World       Date:  2016-10-26

10.  Campylobacter coli From Retail Liver and Meat Products Is More Aerotolerant Than Campylobacter jejuni.

Authors:  Anand B Karki; Daya Marasini; Clark K Oakey; Kaitlin Mar; Mohamed K Fakhr
Journal:  Front Microbiol       Date:  2018-12-12       Impact factor: 5.640

View more
  1 in total

1.  Rapid Oxford Nanopore Technologies MinION Sequencing Workflow for Campylobacter jejuni Identification in Broilers on Site-A Proof-of-Concept Study.

Authors:  Clara Marin; Francisco Marco-Jiménez; Llucia Martínez-Priego; Griselda De Marco-Romero; Vicente Soriano-Chirona; Laura Lorenzo-Rebenaque; Giuseppe D'Auria
Journal:  Animals (Basel)       Date:  2022-08-13       Impact factor: 3.231

  1 in total

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