Literature DB >> 31886353

Data on the differentiation among Leishmania (Viannia) spp., Leishmania (Leishmania) infantum and Leishmania (Leishmania) amazonensis in Brazilian clinical samples using real-time PCR.

Aurora Diotallevi1, Gloria Buffi1, Marcello Ceccarelli1, Herintha Coeto Neitzke-Abreu2, Laisa Vieira Gnutzmann2, Manoel Sebastião da Costa Lima Junior3, Alice Di Domenico1, Mauro De Santi1, Mauro Magnani1, Luca Galluzzi1.   

Abstract

This article contains the data regarding Leishmania species identification in human and canine clinical samples from a Brazilian region endemic for Leishmania (Viannia) spp., Leishmania (Leishmania) infantum and Leishmania (Leishmania) amazonensis, using a previously developed approach involving two qPCR assays (qPCR-ML and qPCR-ama). The data are related to the article "Real-time PCR to differentiate among Leishmania (Viannia) subgenus, Leishmania (Leishmania) infantum and Leishmania (Leishmania) amazonensis: application on Brazilian clinical samples" [1], and include also details of clinical evaluation/diagnosis of human patients and primer sequences used in the qPCR assays. The Leishmania species has been determined in 27 canine samples and 11 human samples, exploiting HRM analysis of qPCR-ML and Cq values of qPCR-ML and qPCR-ama, as reported previously [2]. The qPCR data were in agreement with the species characterization obtained with other methods such as conventional species-specific PCR, ITS1 PCR-RFLP or DNA sequencing. Despite the limited number of clinical samples, these data are encouraging for a potential application in regions where L. (Viannia) spp., L. (L.) infantum and L. (L.) amazonensis are co-endemic.
© 2019 The Author(s).

Entities:  

Keywords:  HRM; Leishmania amazonensis; Leishmania infantum; Real-time PCR; Viannia; qPCR

Year:  2019        PMID: 31886353      PMCID: PMC6920495          DOI: 10.1016/j.dib.2019.104914

Source DB:  PubMed          Journal:  Data Brief        ISSN: 2352-3409


Specifications Table These data are useful to point out the potential field application of our SYBR Green-based qPCR assays to distinguish among Leishmania (Viannia) subgenus, L. (L.) infantum and L. (L.) amazonensis, exploiting HRM and Cq values. The rapid Leishmania species identification can be particularly useful for diagnosis in regions where Leishmania (Viannia) subgenus, L. (L.) infantum and L. (L.) amazonensis are co-endemic. These data can be considered to further develop qPCR-based assays for other species differentiation in the Viannia subgenus (e.g. L. (V.) braziliensis and L. (V.) panamensis) or in the Leishmania subgenus (e.g. L. (L.) amazonensis and L. (L.) mexicana).

Data

The presented data first describe the Clinical evaluation/diagnosis of human patients (Table 1) and the primers used in the qPCR assays targeting Leishmania minicircle kDNA (qPCR-ML and qPCR-ama) and host genes (Table 2). Then, the Cq and HRM values obtained from qPCR-ML and qPCR-ama used for Leishmania species identification in canine (n = 36) and human (n = 11) clinical samples spotted on filter paper are presented (Table 3 and Table 4, respectively). Single replicate data for canine and human samples are represented in supplementary file 1 and 2, respectively. As a further confirmation of Leishmania infantum species identification, ITS1 amplicons from two human clinical samples were directly sequenced (supplementary file 3 and 4) and a phylogenetic analysis was performed with ITS1 sequences from L. (V.) braziliensis (n = 21) L. (L.) infantum (n = 7) and L. (L.) amazonensis (n = 4) available in genbank database (Fig. 1). Moreover, the data also show the specificity of qPCR-ML and qPCR-ama assays in the conditions used to amplify DNA samples from filter paper (Fig. 2).
Table 1

Diagnosis and clinical evaluation of human patients.

Sample IDDiagnosisPatient clinical evaluation
E1VLAnemia, splenomegaly
E2, E3, E4VLNo clinical information
E5VLHIV patient, fever, pancytopenia, hyporexia
E6, E7VLPancytopenia, hepatosplenomegaly
E8VL5 years old; hepatosplenomegaly, anemia, submandibular ganglia, fever, pancytopenia, abdominal distension
E9n.a.No clinical information
F2CLNasal lesion
F3VLAnemia, splenomegaly
Table 2

Primer sequences for qPCR. Primer names are in bold.

TargetForward primer (5′-3′)Reverse primer (5′-3′)
B2M canineGTCCCACAGATCCCCCAAAGCTGGTGGATGGAACCCTGAC
GAPDH humanCCATGTTCGTCATGGGTGTGGGTGCTAAGCAGTTGGTGGTG
kDNA (qPCR-ML)MLF - CGTTCTGCGAAAACCGAAAMLR - CGGCCCTATTTTACACCAACC
kDNA (qPCR-ama)LMi-amaF - AAAATGAGTGCAGAAACCCMLR - CGGCCCTATTTTACACCAACC

B2M canine, Canis familiaris Beta-2-Microglobulin.

GAPDH human, Homo sapiens glyceraldehyde-3-phosphate dehydrogenase.

Table 3

qPCR-ML and qPCR-ama data in canine clinical samples.

Sourceclinical sample (ID)qPCR-ML (Cq ± SD)qPCR-ama (Cq ± SD)qPCR-ML (HRM Tm °C)Species identification
Peripheral bloodA1an.d.n.d.n.d.
A2an.d.n.d.n.d.
A3an.d.n.d.n.d.
A4a20.58 ± 0.6231.83 ± 0.8984.08 ± 0.11L. (L.) infantum
A5a26.51 ± 0.40n.d.84.12 ± 0.02L. (L.) infantum
A6an.d.n.d.n.d.
A7an.d.n.d.n.d.
A8an.d.n.d.n.d.
A9a27.60 ± 1.3939.23 ± 0.8184.03 ± 0.04L. (L.) infantum
B1an.d.n.d.n.d.
B2an.d.n.d.n.d.
B3a23.69 ± 0.1432.28 ± 0.6683.64 ± 0.01L. (L.) infantum
B4a18.16 ± 0.3524.77 ± 0.1683.95 ± 0.32L. (L.) infantum
Bone marrowB5a17.96 ± 0.1726.09 ± 0.0983.94 ± 0.08L. (L.) infantum
B6a22.90 ± 0.3030.34 ± 0.0884.08 ± 0.11L. (L.) infantum
B7a17.39 ± 4.5325.93 ± 3.4783.79 ± 0.02L. (L.) infantum
B8a18.18 ± 3.6225.12 ± 2.6484.05 ± 0.14L. (L.) infantum
B9a25.58 ± 1.2033.57 ± 0.8684.90 ± 0.04L. (L.) infantum
C1a17.41 ± 3.3926.27 ± 2.2583.85 ± 0.02L. (L.) infantum
C2a14.80 ± 0.2322.62 ± 0.0183.96 ± 0.08L. (L.) infantum
C3a16.02 ± 0.1323.73 ± 0.0784.12 ± 0.02L. (L.) infantum
C4a17.53 ± 0.0625.20 ± 0.0683.97 ± 0.07L. (L.) infantum
C5a14.57 ± 0.2222.43 ± 0.1383.90 ± 0.07L. (L.) infantum
C6a25.95 ± 10.55n.d.84.86 ± 0.05L. (L.) infantum
C7a20.38 ± 0.5026.40 ± 3.5784.10 ± 0.05L. (L.) infantum
C8a22.52 ± 0.50n.d.84.05 ± 0.05L. (L.) infantum
C9a30.81 ± 0.35n.d.84.10 ± 0.07L. (L.) infantum
D1an.d.n.d.n.d.
D2a28.51 ± 3.35n.d.84.03 ± 0.05L. (L.) infantum
D3b20.04 ± 0.1232.36 ± 0.3580.44 ± 0.04L. (L.) infantum
D4b17.65 ± 0.1828.62 ± 0.9380.08 ± 0.09L. (L.) infantum
D5b12.14 ± 0.0123.68 ± 0.3279.95 ± 0.03L. (L.) infantum
D6b12.17 ± 0.0922.63 ± 0.1780.13 ± 0.11L. (L.) infantum
D7b8.72 ± 0.0118.85 ± 0.1380.06 ± 0.04L. (L.) infantum
D8b8.38 ± 0.1818.50 ± 0.0580.07 ± 0.08L. (L.) infantum
D9a30.75 ± 0.04n.d.84.10 ± 0.08L. (L.) infantum

n.d. = not detectable.

Samples tested with SYBR green PCR master mix (Diatheva srl).

Samples tested with RT2 SYBR Green ROX FAST Mastermix (Qiagen).

Table 4

qPCR-ML and qPCR-ama data in human clinical samples.

SourceHuman clinical sample (ID)qPCR-ML (Cq ± SD)qPCR-ama (Cq ± SD)qPCR-ML (HRM Tm °C)Species identification
Peripheral bloodE1a28.28 ± 0.7234.08 ± 0.7184.20 ± 0.05L. (L.) infantum
E5b28.81 ± 2.36n.d.83.27 ± 0.05dL. (L.) infantum
83.15 ± 0.07L. (V.) sp
E6a19.12 ± 0.0627.35 ± 0.2384.20 ± 0.01L. (L.) infantum
E7a31.12 ± 0.3134.93 ± 0.7483.94 ± 0.08L. (L.) infantum
E8b23.64 ± 0.46d24.06 ± 0.86d83.77 ± 0.04dL. (L.) infantum
82.75 ± 0.28L. (V.) sp
F2a33.85 ± 4.11n.d.84.18 ± 0.04L. (L.) infantum
F3a29.19 ± 0.18n.d.84.14 ± 0.06L. (L.) infantum
Bone marrowE2c30.32 ± 0.96n.d.80.06 ± 0.08L. (L.) infantum
E3cn.d.30.35 ± 1.55n.d.L. (L.) amazonensis
E4c31.99 ± 1.31n.d.80.65 ± 0.05L. (L.) infantum
E9a27.78 ± 1.34n.d.84.62 ± 0.07L. (L.) infantum

n.d. = not detectable.

Samples tested with SYBR green PCR master mix (Diatheva srl).

Samples tested with SYBR green PCR master mix (Diatheva srl) and TB Green premix ex TaqII Mastermix (Takara).

Samples tested with RT2 SYBR Green ROX FAST Mastermix (Qiagen).

values obtained with TB Green premix ex TaqII Mastermix (Takara).

Fig. 1

Phylogenetic analysis of ITS1 fragments of F2 and F3 samples (red dotted line). The phylogenetic tree was constructed by using the maximum likelihood method and Tamura-Nei model. Numbers below branches represent bootstrap support.

Fig. 2

Specificity of qPCR-ML (A) and qPCR-ama (B). The qPCR amplicons were run on a 2% agarose gel. Both qPCR assays were performed using three different PCR master mix: SYBR green PCR master mix, Diatheva (I); RT2 SYBR Green ROX FAST Mastermix, Qiagen (II); TB Green premix ex TaqII Mastermix, Takara Bio (III). For each master mix, DNA from L. (L.) infantum (0,006 ng/μl), L. (L.) amazonensis (0,15 ng/μl), Trypanosoma cruzi (0,1 ng/μl) and human DNA (30 ng/μl) were tested with the condition described in the manuscript. As negative control, a no template control was used for each qPCR run. M: ΦX174 DNA/BsuRI (HaeIII) Marker, 9 (ThermoFisher Scientific); NTC: no template control.

Diagnosis and clinical evaluation of human patients. Primer sequences for qPCR. Primer names are in bold. B2M canine, Canis familiaris Beta-2-Microglobulin. GAPDH human, Homo sapiens glyceraldehyde-3-phosphate dehydrogenase. qPCR-ML and qPCR-ama data in canine clinical samples. n.d. = not detectable. Samples tested with SYBR green PCR master mix (Diatheva srl). Samples tested with RT2 SYBR Green ROX FAST Mastermix (Qiagen). qPCR-ML and qPCR-ama data in human clinical samples. n.d. = not detectable. Samples tested with SYBR green PCR master mix (Diatheva srl). Samples tested with SYBR green PCR master mix (Diatheva srl) and TB Green premix ex TaqII Mastermix (Takara). Samples tested with RT2 SYBR Green ROX FAST Mastermix (Qiagen). values obtained with TB Green premix ex TaqII Mastermix (Takara). Phylogenetic analysis of ITS1 fragments of F2 and F3 samples (red dotted line). The phylogenetic tree was constructed by using the maximum likelihood method and Tamura-Nei model. Numbers below branches represent bootstrap support. Specificity of qPCR-ML (A) and qPCR-ama (B). The qPCR amplicons were run on a 2% agarose gel. Both qPCR assays were performed using three different PCR master mix: SYBR green PCR master mix, Diatheva (I); RT2 SYBR Green ROX FAST Mastermix, Qiagen (II); TB Green premix ex TaqII Mastermix, Takara Bio (III). For each master mix, DNA from L. (L.) infantum (0,006 ng/μl), L. (L.) amazonensis (0,15 ng/μl), Trypanosoma cruzi (0,1 ng/μl) and human DNA (30 ng/μl) were tested with the condition described in the manuscript. As negative control, a no template control was used for each qPCR run. M: ΦX174 DNA/BsuRI (HaeIII) Marker, 9 (ThermoFisher Scientific); NTC: no template control.

Experimental design, materials, and methods

Canine and human clinical samples

Thirty-six canine clinical samples and 11 human clinical samples were collected in Mato Grosso do Sul (Brazil), an endemic area of leishmaniasis. The canine samples consisted in 13 peripheral blood (A1-A9, B1–B4) and 23 bone marrow (B5–B9, C1–C9, D1-D9), collected by the Control Center of Zoonoses (CCZ) of Campo Grande City. Dogs were diagnosed positive for visceral canine leishmaniasis (VCL) by the immunochromatographic Dual-Path Platform (DPP™, Bio-Manguinhos/Fiocruz, Rio de Janeiro, Brazil) and the ELISA test (EIE™; Bio-Manguinhos/Fiocruz, Rio de Janeiro, Brazil), and direct identification of Leishmania amastigotes from Giemsa-stained smears analyzed by optical microscopy. The human samples consisted in 7 peripheral blood (E1, E5-E8, F2, F3) and 4 bone marrow (E2-E4, E9), collected in Hospital Universitário of Dourados City from patients with a diagnosis of leishmaniasis. After evaluation of clinical signs (Table 1), patients were diagnosed using serology rapid test (rK39) (Kalazar Detect™; InBios, Washington, US) and/or direct identification of Leishmania amastigotes from Giemsa-stained smears analyzed by optical microscopy.

DNA extraction

DNA from canine and human clinical samples, positive for Leishmania infection, was obtained as described in [1]. Briefly, 300 μl volume of sample was added to 500 μl of 20% Sodium Dodecyl Sulfate (Sigma-Aldrich) and homogenized. Then, 400 μl of chloroform and 300 μl of protein precipitation solution (3 M potassium acetate, 11% glacial acetic acid) were added. After centrifugation at 10,000×g for 10 min, the supernatant was transferred to new tubes containing 1 ml of cold absolute ethanol, gently homogenized by inversion and centrifuged for 5 min at 10,000×g. The pellets were washed twice with 1 ml of 70% cold ethanol each time and centrifuged for 2 min at 10,000×g. The supernatant was discarded and the pellet was dried in a dry bath (AccuBlock) at 65 °C for 5 min. The DNA pellet was resuspended in TE buffer (10 mM Tris, 1 mM EDTA, pH 8.0) to a final concentration of 200–500 ng/μl, stored at 4 °C for 24 h, and then frozen at −20 °C. For sample transfer, DNA samples (8 μl) were spotted on filter paper (Macherey-Nagel MN 818), air-dried and stored at room temperature until analysis.

Quantitative PCR (qPCR) assays

The approach described in Ceccarelli et al. [2], consisting in running two qPCR reactions in parallel (qPCR-ML and qPCR-ama) to amplify different classes of minicircles, has been applied on sample DNA spotted on filter paper. The workflow of this approach involves first a discrimination between subgenera Viannia and Leishmania based on qPCR-ML HRM analysis, performed with primers MLF and MLR (Table 2); if this analysis indicates subgenus Leishmania, the discrimination between L. (L.) infantum and L. (L.) amazonensis is performed through comparison of qPCR-ML and qPCR-ama Cq values. The qPCR-ama was performed using the forward primer LMi-amaF and the reverse primer MLR (Table 2). All samples were tested blindly. To evaluate the DNA integrity and amplifiability in canine and human samples, canine beta-2-microglobulin (B2M) and human glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were amplified under the same conditions as the qPCR-ML and qPCR-ama (described below). To ensure applicability with samples on filter paper, a pre-amplification step was introduced as follows. A punch of filter paper (2 mm in diameter) was placed in 40 μl SYBR green reaction mixture containing 200 nM of each primer. Three different PCR master mix were tested: SYBR green PCR master mix (Diatheva srl, Fano, Italy), RT2 SYBR Green ROX FAST Mastermix (Qiagen, Hilden, Germany), TB Green premix ex TaqII Mastermix (Takara Bio Europe, France). Tubes were placed in a thermal cycler (GeneAmp PCR System 2700), and pre-amplified under the following conditions: 94 °C for 5 min, 10 cycles at 94 °C for 30 s, 60 °C for 20 s and 72 °C for 20 s. At the end of this pre-amplification step, the tubes were centrifuged for few seconds and placed in ice; the filter paper was removed and the reaction was split into two PCR tubes (20 μl each tube). Then, the tubes were placed in the Rotor-Gene 6000 instrument and amplified as follows: 45 cycles at 94 °C for 30 s, 60 °C for 20 s and 72 °C for 20 s. As reference, PCR mixtures containing template DNA isolated from L. (L.) infantum MHOM/FR/78/LEM75, L. (L.) amazonensis MHOM/BR/00/LTB0016 and L. (V.) braziliensis MHOM/BR/75/M2904 were included in each run. A no template control (NTC) was included for each primer pair reaction. To monitor non-specific products or primer dimers, a melting analysis was performed from 79 to 95 °C at the end of each run, with a slope of 1 °C/s, and 5 s at each temperature. The Cq values were evaluated using the quantification analysis of the RotorGene 6000 software. Trypanosoma cruzi and human DNA were amplified using the conditions described above to confirm qPCR_ML and qPCR-ama specificity with the three master mix used and including the pre-amplification step. The absence of specific amplicons was confirmed by visualization on 2% agarose gel electrophoresis.

High-resolution melt (HRM) analysis

The high-resolution melt (HRM) analysis was performed immediately after the amplification reactions in the Rotor-Gene 6000 instrument. HRM analysis was conducted over the range from 79 °C to 90 °C, rising at 0.1 °C/s and waiting for 2 s at each temperature. Each sample was run in duplicate, and the gain was optimized before melting on all tubes. HRM curve analysis was performed with the derivative of the intensity of fluorescence at different temperatures (dF/dT), after smoothing, with the Rotor-Gene 6000 software. Template DNA isolated from L. (L.) infantum MHOM/FR/78/LEM75 and L. (L.) amazonensis MHOM/BR/00/LTB0016 were used as reference for Leishmania subgenus, while template DNA from L. (V.) braziliensis MHOM/BR/75/M2904 was used as reference for Viannia subgenus. Bins were set to define Tm of amplicons for each species. Automated classification of genotypes (i.e. subgenus Leishmania or Viannia) of unknown samples was performed by the Rotor-Gene software according to the presence of a derivative peak located within a defined temperature bin.

ITS1 DNA sequencing and phylogenetic analysis

To confirm Leishmania species identification, ITS1 fragments, obtained by ITS1-PCR RFLP [3], were sequenced in two human samples (F2 and F3) as follows. The ITS1 amplicons were directly digested with 10 U HaeIII (Thermo Fisher Scientifc) at 37 °C for 3 h. The restriction fragments were visualized on a 3.5% high-resolution MetaPhor (Cambrex) agarose gel stained with GelRed (Biotium, Hayward, CA). The larger fragment was excised from gel, purified using MinElute Gel Extraction kit (Qiagen) and directly sequenced. DNA sequencing was performed using the BigDye Terminator v. 1.1 Cycle Sequencing Kit on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems). The F2 and F3 sequences were aligned with ITS1 sequences from L. (L.) infantum (n = 7), L. (V.) braziliensis (n = 21) and L. (L.) amazonensis (n = 4) available in GenBank database, using BioEdit Sequence Alignment Editor using default options. The phylogenetic analysis of aligned ITS1 fragments was conducted with MEGA 6 software. Phylogenetic relationships were inferred by using the maximum likelihood method and Tamura-Nei model. Bootstrap values were calculated from 100 replications.

Specifications Table

SubjectParasitology
Specific subject areaMolecular diagnostics
Type of dataTableFigure
How data were acquiredThe qPCR assays were performed using Rotor-Gene 6000 instrument (Corbett life science). The qPCR runs were analyzed with Rotor-Gene software version 1.7 to obtain Cq values and High resolution melting temperatures.
Data formatRawAnalyzed
Parameters for data collectionDNA extracted from clinical samples and strains was spotted on filter paper for storage and laboratory transfer. The qPCR assays (qPCR-ML and qPCR-ama), as well as ITS1 PCR, were performed directly from a punch of filter paper, with a pre-amplification step.
Description of data collectionHRM data obtained with qPCR-ML allowed to discriminate between subgenera Viannia and Leishmania. In case of subgenus Leishmania, the discrimination between L. (L.) infantum and L. (L.) amazonensis is performed through comparison of qPCR-ML and qPCR-ama Cq values. Both HRM and Cq data were collected by amplification of DNA extracted from clinical samples and spotted on filter paper.
Data source locationCanine samplesInstitution: Control Center of Zoonoses (CCZ)City/Town/Region: Campo Grande City/Mato Grosso do SulCountry: BrazilHuman samplesInstitution: Hospital UniversitárioCity/Town/Region: Dourados City/Mato Grosso do SulCountry: Brazil
Data accessibilityWith the article
Related research articleAurora Diotallevi, Gloria Buffi, Marcello Ceccarelli, Herintha Coeto Neitzke-Abreu, Laisa Vieira Gnutzmann, Manoel Sebastião da Costa Lima Junior, Alice Di Domenico, Mauro De Santi, Mauro Magnani, Luca Galluzzi. Real-time PCR to differentiate among Leishmania (Viannia) subgenus, Leishmania (Leishmania) infantum and Leishmania (Leishmania) amazonensis: application on Brazilian clinical samples.Acta Tropica, volume 201, January 2020, 105178
Value of the Data

These data are useful to point out the potential field application of our SYBR Green-based qPCR assays to distinguish among Leishmania (Viannia) subgenus, L. (L.) infantum and L. (L.) amazonensis, exploiting HRM and Cq values.

The rapid Leishmania species identification can be particularly useful for diagnosis in regions where Leishmania (Viannia) subgenus, L. (L.) infantum and L. (L.) amazonensis are co-endemic.

These data can be considered to further develop qPCR-based assays for other species differentiation in the Viannia subgenus (e.g. L. (V.) braziliensis and L. (V.) panamensis) or in the Leishmania subgenus (e.g. L. (L.) amazonensis and L. (L.) mexicana).

  3 in total

1.  Real-time PCR to differentiate among Leishmania (Viannia) subgenus, Leishmania (Leishmania) infantum and Leishmania (Leishmania) amazonensis: Application on Brazilian clinical samples.

Authors:  Aurora Diotallevi; Gloria Buffi; Marcello Ceccarelli; Herintha Coeto Neitzke-Abreu; Laisa Vieira Gnutzmann; Manoel Sebastião da Costa Lima; Alice Di Domenico; Mauro De Santi; Mauro Magnani; Luca Galluzzi
Journal:  Acta Trop       Date:  2019-10-10       Impact factor: 3.112

2.  PCR diagnosis and characterization of Leishmania in local and imported clinical samples.

Authors:  Gabriele Schönian; Abedelmajeed Nasereddin; Nicole Dinse; Carola Schweynoch; Henk D F H Schallig; Wolfgang Presber; Charles L Jaffe
Journal:  Diagn Microbiol Infect Dis       Date:  2003-09       Impact factor: 2.803

3.  The use of kDNA minicircle subclass relative abundance to differentiate between Leishmania (L.) infantum and Leishmania (L.) amazonensis.

Authors:  Marcello Ceccarelli; Luca Galluzzi; Aurora Diotallevi; Francesca Andreoni; Hailie Fowler; Christine Petersen; Fabrizio Vitale; Mauro Magnani
Journal:  Parasit Vectors       Date:  2017-05-16       Impact factor: 3.876

  3 in total
  1 in total

1.  Establish an allele-specific real-time PCR for Leishmania species identification.

Authors:  Yun Wu; Mengyuan Jiang; Shaogang Li; Nicholas R Waterfield; Guowei Yang
Journal:  Infect Dis Poverty       Date:  2022-06-02       Impact factor: 10.485

  1 in total

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