Literature DB >> 30732150

Drug resistance mutations and viral load in human immunodeficiency virus type 2 and dual HIV-1/HIV-2 infected patients in Ghana.

Christopher Z Abana1, Kwamena W C Sagoe2, Evelyn Y Bonney1, Edward K Maina1,3, Ishmael D Aziati1,4, Esinam Agbosu1, Gifty Mawuli1, Linda M Styer5, Koichi Ishikawa6, James A M Brandful1, William K Ampofo1.   

Abstract

Antiretroviral therapy (ART) and drug resistance studies worldwide have focused almost exclusively on human immunodeficiency virus type 1 (HIV-1). As a result, there is limited information on ART and drug resistance in HIV-2 patients. In Ghana, the HIV epidemic is characterized by the domination of HIV-1, with cocirculating HIV-2. We, therefore, sought to determine viral load and drug resistance mutations in HIV-2 patients to inform the clinical management of such individuals in Ghana.We used purposive sampling to collect blood from 16 consented patients, confirmed as HIV-2 or HIV-1/2 dual infections by serology. A 2-step real-time RT-PCR assay was used to determine plasma HIV-2 RNA viral loads. For drug resistance testing, nucleic acids were extracted from plasma and peripheral blood mononuclear cells. The reverse transcriptase and protease genes of HIV-2 were amplified, sequenced and analyzed for drug resistance mutations and HIV-2 group.HIV-2 viral load was detected in 9 of 16 patients. Six of these had quantifiable viral loads (range: 2.62-5.45 log IU/mL) while 3 had viral loads below the limit of quantification. Sequences were generated from 7 out of 16 samples. Five of these were classified as HIV-2 group B and 2 as HIV-2 group A. HIV-2 drug resistance mutations (M184V, K65R, Y115F) were identified in 1 patient.This study is the first to report HIV-2 viral load and drug resistance mutations in HIV-2 strains from Ghana. The results indicate the need for continuous monitoring of drug resistance among HIV-2- infected patients to improve their clinical management.

Entities:  

Mesh:

Substances:

Year:  2019        PMID: 30732150      PMCID: PMC6380870          DOI: 10.1097/MD.0000000000014313

Source DB:  PubMed          Journal:  Medicine (Baltimore)        ISSN: 0025-7974            Impact factor:   1.817


Introduction

It is estimated that 36.7 million people worldwide are living with the human immunodeficiency virus (HIV).[ Sub-Saharan Africa is most affected, with 25.8 million people living with HIV in 2015.[ HIV type 1 (HIV-1) and type 2 (HIV-2) infections have been described [ with HIV-1 responsible for over 90% of the global acquired immunodeficiency syndrome (AIDS) pandemic.[ Two million people are known to be infected with HIV-2, most of whom live in West Africa.[ HIV-2 infection has spread to other countries in Europe, Asia, and North America[ with a high prevalence (5.4%) reported in Portugal.[ Coinfections of HIV- 1 and HIV- 2 are common in West Africa representing 0.3% to 1% of HIV-infections.[ HIV-1 and HIV-2 have similar routes of transmission, however, HIV-2 is associated with lower viral loads[ resulting in reduced transmission and slower disease progression.[ The HIV epidemic in Ghana is characterized by cocirculation of both viruses with HIV-1 being predominant.[ In view of this, the national antiretroviral therapy (ART) guidelines for the management of HIV/AIDS patients make provisions specifically for persons with HIV-2 and HIV-1/2 dual infections. For such individuals, the guidelines recommend a protease inhibitor (PI)-based regimen instead of those based on non-nucleoside reverse transcriptase inhibitor (NNRTI).[ In 2017, the estimated adult national HIV prevalence in Ghana was 1.67% with 99% HIV-1, 0.2% HIV-2, and 0.8% HIV-1/2 dual infections.[ The emergence of drug-resistant strains during treatment is a major obstacle to the success of any ART program. Both HIV-1 and HIV-2 are known to develop resistance against antiretroviral drugs.[ The development of drug resistance in HIV-1 infections globally has been examined in several studies[; however, information is scarce on HIV-2 drug resistance.[ Previous studies in Ghana have documented drug resistance mutations in treatment naïve or experienced HIV-1 patients.[ However, there is little drug resistance data on HIV-2 infected patients in Ghana, probably due to the lack of commercially available HIV-2 viral load assays and optimized HIV-2 genotyping protocols. This study therefore optimized available protocols[ and applied it to investigate viral loads and drug resistance mutations in HIV-2 and dual HIV-1/HIV-2 infected patients from Ghana.

Methods

Study design and population

This was a cross-sectional study that recruited patients from ART care centers at the Eastern Regional Hospital, Koforidua and the Korle-Bu Teaching Hospital, Accra from November 2014 to June 2015. The participants were limited to Ghanaian men and women aged 18 years and older. There were 2 cohorts in this study; an ART-naïve group and an ART-experienced group. Patients eligible for the study included patients who were seropositive for HIV-2 or dual HIV-1/2. To be able to identify such patients, we searched for patients who were diagnosed as HIV-2 or dual HIV-1/2 from the hospital folders. Purposive sampling was therefore used to collect blood from 27 consented patients who were seropositive for HIV-2 or dual HIV-1/2 using First Response HIV-1-2-0 card test (Premier Medical Corporation Ltd, India) performed at the hospitals. The status of these patients was confirmed as HIV-2 or dual HIV-1/2 using the innogenetics line immuno assay HIV-I/II score assay (Fujirebio Europe N. V, Belgium). After confirmatory testing, 16 patients; 10 ART-naïve and 6 ART-experienced, who agreed to take part in the study were enrolled after written informed consent was sought. Sample size was determined based on the national prevalence of HIV-2 (0.8%) in Ghana for 2014 using a t value of 1.96 at 95% confidence level and 0.05 confidence interval.[ Clinical histories were retrieved from hospital folders of study patients. The study was conducted in accordance with procedures approved by the Institutional Review Board of the Noguchi Memorial Institute for Medical Research (NMIMR-IRB CPN 063/14-15) and the Ethical and Protocol Review Committee of the University of Ghana Medical School (MS-Et/M.4-P4.3/2014-2015).

Whole blood processing into plasma and peripheral blood mononuclear cells

Blood was processed into plasma and peripheral blood mononuclear cells (PBMC) using a sucrose-gradient based protocol with Histopaque 1077 (Sigma Aldrich Company; Darmstadt, Germany). Plasma was stored at −80oC while PBMCs were stored in freezing medium, consisting of 1% dimethyl sulfoxide (DMSO) (Sigma Aldrich, Germany) in fetal bovine serum (Sigma Aldrich), at −80oC until use.

Nucleic acid extraction and purification

Ribonucleic acid (RNA) and deoxyribonucleic acid (DNA) were extracted from plasma and PBMCs, respectively using the QIAamp viral RNA mini kit and the DNAeasy Blood and Tissue kit (Qiagen, Germany) according to the manufacturer's instructions.

HIV-2 viral load and genotyping

HIV-2 viral load was performed on 0.2 ml of plasma following a published protocol.[ Viral load testing was completed at the Wadsworth Center, New York State Department of Health according to the approved IRB protocol (#17-039). The HIV-2 viral load assay had a lower limit of detection of 32 international units (IU)/ml (1.50 log IU/ml) and a lower limit of quantification (LLOQ) of 225 IU/ml (2.35 log IU/ml). For genotyping, reverse transcriptase (RT) and protease (PR) genes of HIV-2 were amplified separately using specific primers.[ A nested polymerase chain reaction (PCR) for the RT gene from PBMC portions was carried out to amplify a genomic region of 1050 base pairs (bp) encoding the RT gene (Table 1). The reactions were carried out in a total volume of 25 μl containing 5 μl of DNA template, 12.5 μl of Supermix (Life Technologies, Invitrogen; Austin, TX) and 0.5 μl of 20 μM of primers (Table 1). The first round PCR conditions were: 94°C for 2 minutes, 40 cycles of 94°C for 30 seconds, 55°C for 1 minute and 72°C for 1 minute 30 seconds, and then an elongation step at 72°C for 7 minutes. A nested PCR was carried out using 5 μl of the first-round PCR product under the following cycle conditions: 94°C for 2 minutes, 40 cycles of 94°C for 30 seconds, 56°C for 1 minute and 72°C for 1 minute, and an extension at 72°C for 5 minutes. The entire coding region for the PR gene (297 bp) was amplified by nested PCR using primers (Table 1). The cycling conditions for PR gene PCRs were the same as that for RT gene.
Table 1

Details of primers used to amplify HIV-2 sub-genomic fragments previously published.

Details of primers used to amplify HIV-2 sub-genomic fragments previously published. The RT and PR genes in plasma were amplified by nested PCR using QIAGEN OneStep RT-PCR Kit (Qiagen, Germany) for the first round in a total reaction volume of 25 μl containing 5 μl of RNA, 5 μl of 5× buffer, 0.75 μl of 20 μM of primers, 1 μl of enzyme mix and dNTPs followed by nested PCR with Amplitaq Gold master mix kit (Applied Biosystems, USA) in 25 μl reaction volume containing 12.5 μl Amplitaq Gold, 0.5 μl of 20 μM of primers (Table 1). The first round PCR had reverse transcription at 50°C for 30 minutes, and enzyme degradation at 95°C for 15 minutes followed by 40 amplification cycles (94°C for 30 seconds, 55°C for 1 minute, 72°C for 1 minute 30 seconds) and then an extension at 72°C for 7 minutes. The second round PCR conditions were: 94°C for 2 minutes, followed by 40 amplification cycles (94°C for 30 seconds, 56°C for 1 minute and 72°C for 1 minute) and extension at 72°C for 5 minutes. The PCR products for both PBMC and plasma were visualized on 1.5% agarose gel under ultraviolet trans-illumination using the KODAK Gel Logic 100 Imaging system (Eastman Kodak Company; Rochester, New York)

DNA sequencing

Nested PCR products were purified using QIAquick PCR purification system (Qiagen, Hilden, Germany) according to the manufacturer's instructions. Cycle sequencing was performed using the Big Dye terminator v3.1 Cycle Sequencing kit (Applied Biosystems). The DNA fragments were sequenced on both strands with sense primers RT3 and RT5-HIV2 and antisense primers RT4 and RT6-HIV2 for the RT region and primers PR3 and PR2 for the PR region (Table 1). A 10 μl reaction mixture, consisting of 2 μl each of 5X Big Dye Sequencing Buffer, Big Dye Terminator Mix, 2 μM primer, nuclease-free water and purified PCR product, was taken through 94°C for 2 minutes, 25 cycles of 94°C for 30 seconds, 50°C for 15 seconds, and 60°C for 4 minutes. Cycle-sequenced products were purified using AgenCourt CleanSeq kit (Agencourt Bioscience Corporation; Beverly, MA) according to manufacturer's protocol and loaded onto an ABI 3130xl genetic analyzer (Applied Biosystems) to read out sequences.

Sequence data analysis

Sequences were edited using the Sequencher software v4.0 (Gene Codes Corporation; Ann Arbor, MI) and aligned in BioEdit.[ The Fasta-formatted sequences were then analyzed in the HIV basic local alignment search tool programme (http://www.hiv.lanl.gov) for identification of groups. Nucleotide sequences were aligned to the consensus group A HIV-2 ROD and group B HIV-2 EHO. HIV-2 group and drug resistance mutations were obtained by submitting the sequences to an online tool used for HIV-2 drug-resistance interpretation (HIV-2 EUv2).[ Sequences were analyzed for the presence of mutations associated with nucleoside reverse transcriptase inhibitor (NRTI), NNRTI), or PI resistance and were categorized as susceptible, intermediate resistance, or resistant.

Results

Demographic and clinical characteristics

Sixteen HIV-2 and HIV-1/2 dual-infected patients from 2 hospitals in Southern Ghana were enrolled into the study from December 2014 to May 2015. Eleven of these patients were HIV-2 infected and 5 were HIV-1/2 dual infected. Ten were ART-naïve and 6 were ART-experienced. Demographic, clinical and viral load data of study patients are shown in Table 2. Four of the ART-experienced patients were infected with HIV-2 and the remaining 2 were infected with HIV-1/2. Treatment options differed in the 2 ART care centers for both HIV-2 and HIV-1/2 patients. The details of the composition of individual regimen are shown in Figure 1.
Table 2

Demographic, clinical, and molecular characteristics of study participants.

Figure 1

CD4 distribution and treatment regimen of ART-experienced patients. The CD4 at baseline represent the CD4 count (cells/μl) taken at the start of ART. CD4 at enrollment represent the CD4 cell count (cells/μl) taken at the time of enrollment. The duration represents the time interval from the start of ART till enrollment into the study. ART = antiretroviral therapy, NNRTI: EFV = efavirenz, NVP = nevirapine; NRTI: 3TC = lamivudine, AZT = zidovudine, TDF = tenofovir; PI: LPV/r = lopinavir/r.

Demographic, clinical, and molecular characteristics of study participants. CD4 distribution and treatment regimen of ART-experienced patients. The CD4 at baseline represent the CD4 count (cells/μl) taken at the start of ART. CD4 at enrollment represent the CD4 cell count (cells/μl) taken at the time of enrollment. The duration represents the time interval from the start of ART till enrollment into the study. ART = antiretroviral therapy, NNRTI: EFV = efavirenz, NVP = nevirapine; NRTI: 3TC = lamivudine, AZT = zidovudine, TDF = tenofovir; PI: LPV/r = lopinavir/r.

HIV-2 viral load and CD4 cell counts

HIV-2 RNA was detected in the plasma of 9 out of the 16 study patients; 7 were ART-naïve and 2 were on treatment. Three had detectable RNA below the LLOQ and 6 had quantifiable viral loads ranging from 2.62 to 5.45 log IU/ml (Table 2). The highest viral load (5.45 log IU/ml) was recorded in an ART-naive patient; the highest viral load in an ART-experienced patient was 4.62 log IU/ml. Generally, cluster of differentiation 4 (CD4) cell counts ranged from 11 to 916 cells/μl. The least (11 cells/μl) and highest (916 cells/μl) CD4 cell count were recorded in the ART-naïve group (Table 2). Patients CD4 cell count taken at baseline and enrollment are shown in Figure 1.

PCR and sequencing of the RT and PR genes from plasma and PBMC

The RT and PR genes were successfully amplified in 5 and 4 samples, respectively. For the RT gene, 4 of the samples amplified were from PBMC and 1 was from plasma. All 5 of these samples were successfully sequenced. Included in this group were sequences obtained from both plasma and PBMC of 1 HIV-2 ART-naïve patient. For the PR gene, 2 were amplified from PBMC and the other 2 from plasma, all from ART-naïve patients. Sequences were only obtained from the PBMC amplicon (Table 3).
Table 3

HIV-2 groups and drug resistance mutations from the ART-naïve and experienced patients.

HIV-2 groups and drug resistance mutations from the ART-naïve and experienced patients.

HIV-2 grouping and drug resistance mutation analysis

Five sequences from the RT region were identified as HIV-2 group B, including those obtained from the paired plasma and PBMC portions of the blood. Two PR sequences were both HIV-2 group A. One patient (KBAN-14-01) had RT and PR sequences that were classified in different groups; RT as group B and PR as group A (Table 3). There were no drug resistance-conferring mutations in any of the sequences derived from ART-naïve patients. In the ART-experienced patient, 3 major drug resistance mutations (M184V, K65R, and Y115F) were found in both the PBMC and plasma sequences. These mutations were classified as NNRTI resistant (Table 3)

Discussion

Attempts to curb the HIV menace globally are almost exclusively focused on HIV-1. The need for HIV-2 viral load and genotyping protocols to monitor HIV-2 infected patients in regions where HIV-1 and HIV-2 cocirculate is imperative. In this study, we estimated HIV-2 viral loads and detected drug resistance mutations by sequencing the RT and PR genes in HIV-2 strains from HIV-2 and HIV-1/2 dual-infected patients in Ghana. HIV-2 strains in our study were found to belong to groups A and B, with group B predominating thus supporting previous reports on the HIV-2 epidemic.[ In addition, we also identified 1 intergroup (A/B) recombinant HIV-2 strain. This strain was classified as group A in the PR gene and group B in the RT gene. Similar recombinant HIV-2 strains had been documented previously in Cote d’Ivoire and Cameroon.[ Our study identified matching HIV-2 drug resistance mutations in plasma and PBMC from the same patient. A previous study also documented matching drug resistance mutations in paired PBMC and plasma for HIV-1.[ This finding from our study supports the recommendations for the use of PBMC as an alternative to plasma for resistance testing, particularly in patients with low plasma viral loads.[ Major drug resistance mutations M184V, K65R, and Y115F were detected in an ART-experienced patient. These key mutations had been reported in genotyping analysis elsewhere.[ The K65R and Q151M substitutions are known to occur more commonly in HIV-2 than in HIV-1 infections,[ and often appear together with M184V in patients that receive NRTI, as shown in an ART-experienced patient in this study. The presence of K65R is indicative of class-wide nucleoside analogue resistance. The occurrence of Y115F, abacavir-induced mutation, together with M184V, which reduces the susceptibility of lamivudine (3TC), emtricitabine, and abacavir confers possible resistance to the drugs in the NRTIs class. This ART-experienced patient was on a regimen containing 2 NRTIs (lamivudine; tenofovir [TDF]) and 1 NNRTI (nevirapine). Such a patient did probably not derive much therapeutic benefit, because HIV-2 strains are known to be naturally resistant to NNRTIs.[ This patient thus received dual therapy of 3TC and TDF and was thus more prone to accumulate NRTI resistance mutations that will eventually lead to treatment failure. The T cell CD4 count for this patient (76 cells/ul) indicated immunological failure but this information was not available to the clinician at the time of sampling due to nationwide shortage of CD4 reagents; the information only became available because of our study. Ideally, patients with CD4 cell count less than 350 cells/μl should initiate therapy as per the ART guidelines in Ghana but this was observed due to an acute shortage of CD4 reagents. Based on the Ghana National ART guidelines, this HIV-2 patient should not have been on an NNRTI-based regimen. However, the clinician managed the patient with respect to HIV-1 infection; the predominant virus. This patient had the highest HIV-2 viral load (4.62 IU/ml) among ART-experienced patients and thus a high possibility of transmitting drug-resistant HIV-2. Availability of HIV typing or HIV-2 viral load data at the time of enrollment for this patient would have probably informed a more suitable regimen and curtailed the emergence of resistance mutations. A high viral load (5.45 IU/ml) was also detected in an ART-naïve patient indicating potential high risk of HIV-2 transmission. The lack of an HIV-2 viral load assay to monitor the progress of patients on therapy at the time of sampling is a limitation to this study. Generally, the lack of commercially available assay for HIV-2 viral load is a major limitation in monitoring HIV-2 infected patients on therapy. Our study did not find any drug resistance mutations in ART-naïve patients, implying the absence of transmitted drug resistance. This finding is in contrast to observations in previous studies in Burkina Faso, Guinea Bissau, Belgium, and Luxemburg, where transmitted drug-resistant HIV-2 strains were reported in ART-naïve patients.[ The presence of NRTI mutations M184V, K65R, and Y115F found in an ART-experienced person in this study confirms the emergence of ART resistance to the currently available antiretroviral drugs under drug pressure. This is the first report on HIV-2 viral load and drug resistance mutations in patients infected with either HIV-2 or HIV-1/2 in Ghana. The study showed evidence of major mutations in HIV-2 strains of patients receiving HIV-1 targeted ART. This underscores the need for continuous monitoring of HIV-2 resistance mutations to improve clinical management of HIV-2 patients.

Author contributions

Conceptualization: Christopher Z. Abana, Evelyn Y. Bonney, James A. M. Brandful, William K. Ampofo. Funding acquisition: Christopher Z. Abana, William K. Ampofo. Investigation: Christopher Z. Abana. Methodology: Christopher Z. Abana, Evelyn Y. Bonney, Ishmael D. Aziati, Esinam Agbosu, Gifty Mawuli, Linda M. Styer, James A. M. Brandful, William K. Ampofo. Resources: Koichi Ishikawa, James A. M. Brandful, William K. Ampofo. Software: Evelyn Y. Bonney, Ishmael D. Aziati, Esinam Agbosu. Supervision: Kwamena W. C. Sagoe, Evelyn Y. Bonney, Linda M. Styer, James A. M. Brandful, William K. Ampofo. Writing – original draft: Christopher Z. Abana. Writing – review and editing: Kwamena W. C. Sagoe, Evelyn Y. Bonney, Edward K. Maina, Ishmael D. Aziati, Esinam Agbosu, Gifty Mawuli, Linda M. Styer, James A. M. Brandful, William K. Ampofo.
  45 in total

1.  Low level of transmitted HIV Drug resistance at two HIV care centres in Ghana: a threshold survey.

Authors:  E Y Bonney; N A Addo; N A A Ntim; F Addo-Yobo; P Bondzie; K-E Aryee; J Barnor; J Brandful; V Bekoe; S-A Ohene; W Ampofo
Journal:  Ghana Med J       Date:  2013-06

2.  Low plasma human immunodeficiency virus type 2 viral load is independent of proviral load: low virus production in vivo.

Authors:  S J Popper; A D Sarr; A Guèye-Ndiaye; S Mboup; M E Essex; P J Kanki
Journal:  J Virol       Date:  2000-02       Impact factor: 5.103

3.  HIV-1/HIV-2 co-infection among voluntary counselling and testing subjects at a regional hospital in Cameroon.

Authors:  D S Nsagha; A L Njunda; H L F Kamga; J C N Assob; E A Bongkem
Journal:  Afr Health Sci       Date:  2012-09       Impact factor: 0.927

4.  Virological Characterization of Dual HIV-1/HIV-2 Seropositivity and Infections in Southern Ghana.

Authors:  Kwc Sagoe; Jaa Mingle; Rk Affram; S Britton; A Dzokoto; A Sonnerborg
Journal:  Ghana Med J       Date:  2008-03

5.  Clinical, immunological and virological response to different antiretroviral regimens in a cohort of HIV-2-infected patients.

Authors:  Marchina E van der Ende; Jan M Prins; Kees Brinkman; Monique Keuter; Jan Veenstra; Sven A Danner; Hubert G M Niesters; Albert D M E Osterhaus; Martin Schutten
Journal:  AIDS       Date:  2003-07       Impact factor: 4.177

6.  HIV type 2 intergroup recombinant identified in Cameroon.

Authors:  Julie Yamaguchi; Ana Vallari; Nicaise Ndembi; Ruthie Coffey; Charlotte Ngansop; Dora Mbanya; Lazare Kaptué; Lutz G Gürtler; Sushil G Devare; Catherine A Brennan
Journal:  AIDS Res Hum Retroviruses       Date:  2008-01       Impact factor: 2.205

7.  Virological response to highly active antiretroviral therapy in patients infected with human immunodeficiency virus type 2 (HIV-2) and in patients dually infected with HIV-1 and HIV-2 in the Gambia and emergence of drug-resistant variants.

Authors:  Sabelle Jallow; Abraham Alabi; Ramu Sarge-Njie; Kevin Peterson; Hilton Whittle; Tumani Corrah; Assan Jaye; Matthew Cotten; Guido Vanham; Samuel J McConkey; Sarah Rowland-Jones; Wouter Janssens
Journal:  J Clin Microbiol       Date:  2009-05-06       Impact factor: 5.948

8.  Validation and clinical use of a sensitive HIV-2 viral load assay that uses a whole virus internal control.

Authors:  Linda M Styer; Thomas T Miller; Monica M Parker
Journal:  J Clin Virol       Date:  2013-09-13       Impact factor: 3.168

9.  Use of peripheral blood DNA for genotype antiretroviral resistance testing in drug-naive HIV-infected subjects.

Authors:  Ilaria Vicenti; Francesca Razzolini; Francesco Saladini; Laura Romano; Maurizio Zazzi
Journal:  Clin Infect Dis       Date:  2007-05-08       Impact factor: 9.079

10.  Transmitted drug resistance, selection of resistance mutations and moderate antiretroviral efficacy in HIV-2: analysis of the HIV-2 Belgium and Luxembourg database.

Authors:  Jean Ruelle; François Roman; Anne-Thérèse Vandenbroucke; Christine Lambert; Katrien Fransen; Fedoua Echahidi; Denis Piérard; Chris Verhofstede; Kristel Van Laethem; Marie-Luce Delforge; Dolorès Vaira; Jean-Claude Schmit; Patrick Goubau
Journal:  BMC Infect Dis       Date:  2008-02-27       Impact factor: 3.090

View more
  1 in total

1.  HIV-2 Drug Resistance Genotyping from Dried Blood Spots.

Authors:  Dana N Raugi; Robert S Nixon; Sally Leong; Khadim Faye; Jean Phillipe Diatta; Fatima Sall; Robert A Smith; ElHadji Ibrahima Sall; Jean Jacques Malomar; Moussa Seydi; Geoffrey S Gottlieb
Journal:  J Clin Microbiol       Date:  2020-12-17       Impact factor: 5.948

  1 in total

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