Impact of UGT1A1 gene polymorphisms on plasma dolutegravir trough concentrations and neuropsychiatric adverse events in Japanese individuals infected with HIV-1
BMC Infectious Diseases volume 17, Article number: 622 (2017)
Dolutegravir (DTG) is metabolized mainly by uridine diphosphate (UDP)-glucuronosyltransferase 1A1 (UGT1A1), and partly by cytochrome P450 3A (CYP3A). Therefore, we focused on UGT1A1 gene polymorphisms (*6 and *28) in Japanese individuals infected with human immunodeficiency virus (HIV)-1 to examine the relationship between their plasma trough concentration of DTG and gene polymorphisms. Recently, neuropsychiatric adverse events (NP-AEs) after the use of DTG have become a concern, so the association between UGT1A1 gene polymorphisms and selected NP-AEs was also investigated.
The study subjects were 107 Japanese patients with HIV-1 infections who were receiving DTG. Five symptoms (dizziness, headache, insomnia, restlessness, and anxiety) were selected as NP-AEs. The subjects were classified by their UGT1A1 gene polymorphisms for the group comparison of DTG trough concentration and the presence or absence of NP-AEs.
The subjects consisted of eight (7%) *6 homozygotes, three (3%) *28 homozygotes, four (4%) for *6/*28 compound heterozygotes, 23 (21%) *6 heterozygotes, 18 (17%) *28 heterozygotes, and 51 (48%) patients carrying the normal allele. The plasma DTG trough concentration of the *6 homozygous patients was significantly higher than that of the patients carrying the normal allele (median, 1.43 and 0.82 μg/mL, respectively, p = 0.0054). The *6 and *28 heterozygous patients also showed significantly higher values than those shown by patients with the normal allele. Multivariate analysis revealed that carrying one or two UGT1A1*6 gene polymorphisms, one UGT1A1*28 polymorphism, and age of < 40 years were independent factors associated with high DTG trough concentrations. The median DTG trough concentration was significantly higher in the patients with NP-AEs (1.31 μg/mL) than in those without NP-AEs (1.01 μg/mL). Consistent with these results, subjects carrying UGT1A1*6, UGT1A1*28, or both alleles showed a higher cumulative incidence of having selected NP-AEs than those carrying the normal alleles (p = 0.0454).
In addition to younger age, carrying UGT1A1*6 and/or UGT1A1*28 was demonstrated to be a factor associated with high DTG trough concentrations. Our results also suggest a relationship between plasma DTG trough concentrations and NP-AEs, and that carrying UGT1A1*6 and/or UGT1A1*28 alleles might be a risk factor for NP-AEs.
Dolutegravir (DTG) is a second-generation human immunodeficiency virus (HIV)-1 integrase inhibitor. In patients without drug-resistant mutations related to integrase inhibitors, once daily administration of DTG shows a good virological effect [1,2,3,4]. Some clinical trials have demonstrated the virological superiority of DTG in comparison with that of other agents [1,2,3,4]. Therefore, major guidelines list DTG as a recommended key drug in the initial therapy of HIV infection [5, 6]. DTG is metabolized mainly by uridine diphosphate (UDP)-glucuronosyltransferase 1A1 (UGT1A1) and partly by cytochrome P450 (CYP) 3A. DTG does not have inhibitory or inducing effects on CYP and UGT, and therefore, has a lower propensity for drug-drug interactions compared to other antiretroviral drugs .
A previous population pharmacokinetic analysis has shown a negative correlation between total bilirubin and DTG’s apparent clearance, which has been attributed to the competition of bilirubin and DTG for UGT1A1. Hence, it is likely that changes in UGT1A1 enzymatic activity and its protein level might affect plasma DTG concentrations .
The UGT1A1 gene has several polymorphisms. The promoter region of UGT1A1 normally has six TA repeats in the TATA box, whereas the number is increased to seven in UGT1A1*28, which lowers the mRNA and protein expression levels of UGT1A1 . In addition, the single nucleotide polymorphism in exon 1 of UGT1A1*6 reduces UGT1A1 activity. The frequency of the UGT1A1 gene polymorphisms has been demonstrated to vary depending on race. UGT1A1*28 is relatively common, whereas UGT1A1*6 is rarely found in Caucasian and African American populations . On the other hand, UGT1A1 *6 frequency is high in Asian populations, including Japanese [11, 12].
In a previous study conducted in Caucasian and African American populations, the area under the DTG plasma concentration-time curve (AUC) and the maximum DTG plasma concentration (Cmax) of patients homozygous for UGT1A1*28 were approximately 1.3–1.4 times higher than those of patients carrying normal alleles . However, the effect of the UGT1A1*6 allele was not fully analyzed because few subjects carrying that allele were found in that study. Therefore, we examined the effects of carrying the UGT1A1*6, *28, or both gene polymorphisms on plasma trough concentrations of DTG in Japanese individuals with HIV-1 infections.
Recently, neuropsychiatric adverse events (NP-AEs) after the use of DTG have become a concern [14,15,16,17,18,19,20]. Although a low discontinuation rate of DTG has been reported in phase III clinical trials [1,2,3,4], NP-AEs are known to be a major cause of DTG discontinuation in real-life setting. Therefore, the association between UGT1A1 gene polymorphisms and selected NP-AEs was also investigated.
This study included Japanese patients infected with HIV-1 who visited the National Hospital Organization Osaka National Hospital, who were 20 years old or older, and who received anti-HIV therapy in combination with DTG. Written consent was acquired prior to sample collection and testing. The subjects were limited to those taking DTG in the morning, because of the time limitations of our institution on the collection and processing of patient samples. We interviewed them about their medication adherence to exclude subjects with poor adherence. Subjects concomitantly using a potent inhibitor of UGT1A1 (atazanavir), potent inducers of UGT1A1, CYP3A, or both (rifampicin, carbamazepine, phenytoin, and phenobarbital), potent inhibitors of CYP3A (itraconazole, cobicistat, and ritonavir), and any other drugs that could affect plasma DTG concentrations (aluminum- or magnesium-containing antacids or both) were excluded.
UGT1A1 Genotyping and plasma trough concentration of DTG
Saliva samples were collected from the patients, absorbed with a filter paper, and subsequently dried. A portion of the filter paper containing the saliva sample was cut and used as a template in the polymerase chain reaction (PCR) used to determine the presence or absence of gene polymorphisms using a previously reported method .
Blood samples were collected after oral administration of DTG for 10 days or longer. Whole blood samples (5 mL) were collected in test tubes containing heparin sodium 22–26 h after an oral administration of DTG. The blood samples were centrifuged at 1200×g for 10 min to isolate 2 mL plasma, which was stored in a freezer at −80°C until assayed. We measured the DTG plasma concentration using mass spectrometry based on our previously developed assay for raltegravir (RAL) , an integrase inhibitor. DTG was detected at a mass-to-charge ratio (m/z) of 420 using the MassLynx analysis software version 4.0 (Waters Corp., Milford, MA).
Examination of the relationship of UGT1A1 genotype and DTG plasma trough concentration to NP-AEs
Five symptoms (dizziness, headache, insomnia, restlessness, and anxiety) were selected as NP-AEs. The information on NP-AEs was obtained upon patient entry into this study, and the subjects were followed-up until March 2017 or the discontinuation of DTG. NP-AEs, which had developed after starting DTG but had ceased by initiation of the study, were included in the analyses. Plasma DTG trough concentrations were compared between the subjects with or without NP-AEs.
For group comparisons, the Wilcoxon rank-sum test was used. If a statistically significant result was obtained for multiple group comparison, Steel’s test was used as the post-hoc test. To search for factors associated with trough concentrations of DTG, multivariate analysis using logistic regression was conducted. The subjects were classified into two groups using the median trough DTG concentration (1.06 μg/mL) as a cut-off value. The variables used were age (over 40 years old), weight (less than 60 kg), CD4 count (less than 500/μL), and the number of *6 and *28 alleles; they were analyzed using the forced input method. The cut-off value for age was set to 40 years old to match that of a previously published study . To search for factors associated with selected NP-AEs, the cumulative incidence of NP-AEs was estimated by the Kaplan-Meier method. The subjects were divided into two groups, stratified by the presence or absence of carrying UGT1A1*6, UGT1A1*28, or both alleles (normal allele group, −/− versus reduced-function allele group, −/*28, −/*6, *28/*28, *6/28, and *6/*6), and these two groups were compared. We adopted a significance level of 5%, and the statistical analysis was conducted using JMP Software version 10.0.0 (SAS Institute Inc., Cary, NC). The p value of χ2 in the analyses of Hardy-Weinberg equilibrium was calculated using R software, version 3.3.3 (R Foundation for Statistical Computing, Vienna, Austria).
Frequency of UGT1A1 gene polymorphisms in Japanese patients with HIV-1 infection
Table 1 shows the general characteristics of the 107 individuals infected with HIV-1 who participated in this study based on the type of UGT1A1 gene polymorphism. The subjects comprised eight (7%) *6 homozygotes, three (3%) *28 homozygotes, four (4%) *6/*28 compound heterozygotes, 23 (21%) *6 heterozygotes, 18 (17%) *28 heterozygotes, and 51 (48%) patients homozygous for the wildtype UGT1A1 genotype. The *6 and *28 allele frequencies were 20% and 13%, respectively. Distribution of the *6 and *28 allele was in the Hardy-Weinberg equilibrium (p = 0.83 and p = 0.92, respectively). Table 1 shows that there were no significant differences among UGT1A1 gene polymorphism frequencies with respect to the listed participant characteristics, except for rilpivirine use that was due to one individual.
Relationship between UGT1A1 gene polymorphisms and plasma DTG concentrations
The median DTG trough concentration of the 107 subjects was 1.06 μg/mL. First, the DTG trough concentrations were classified based on UGT1A1 gene polymorphisms and group comparisons were conducted (Fig. 1a, Wilcoxon rank-sum test, p = 0.0008). The median DTG trough concentration of the *6 homozygous patients (1.43 μg/mL) was significantly higher than that of patients carrying both normal alleles (0.82 μg/mL, p = 0.0054). The median concentrations of the *6 and *28 heterozygous patients (1.29 and 1.2 μg/mL, respectively) were also significantly higher than that of the homozygous patients carrying normal alleles (p = 0.0223 and p = 0.0254, respectively). Moreover, no statistically significant difference was detected in median DTG trough concentrations between patients with both normal alleles and those with both *28 alleles (1.18 μg/mL) or between those with both normal alleles and *6/*28 compound heterozygotes (1.37 μg/mL). However, we could not rule out the possibility that the statistical power was insufficient due to the small numbers of patients.
Next, a multivariate logistic regression analysis was conducted using the median DTG trough concentration, 1.06 μg/mL, as the cut-off value to find factors associated with high DTG trough concentrations. Carrying one or two UGT1A1*6 gene alleles and a UGT1A1*28 allele, as well as an age < 40 years were found to be independent factors associated with high DTG trough concentrations (Table 2).
Examination of the factors associated with NP-AEs
During the follow-up period (255 person-years), selected NP-AEs after the use of DTG were observed in 34 subjects (32%, Table 3). Insomnia was the most frequent (n = 14, 13%) among selected NP-AEs, followed by dizziness (n = 9, 8%), and headache (n = 9, 8%). The severity of all NP-AEs was grade 1 or 2. Six subjects discontinued DTG, and all were because of NP-AEs. Two subjects were lost to follow-up because of personal reasons. As shown in Fig. 1b, plasma DTG trough concentration was significantly higher in the patients with NP-AEs (median, 1.31 μg/mL) than in those without NP-AEs (median, 1.01 μg/mL, Wilcoxon rank-sum test, p = 0.0013). Next, we estimated the factors associated with the cumulative incidence of selected NP-AEs by use of Kaplan-Meier curves (Fig. 2). Plasma DTG trough concentration (≥ 1.06 μg/mL), and carrying UGT1A1*6, UGT1A1*28, or both alleles were associated with higher cumulative incidences of selected NP-AEs (log rank-test, p = 0.0105, and 0.0454, respectively). Age did not show this association (log rank-test, p = 0.2477).
The rate of carriers with the UGT1A1 gene polymorphism has been reported to vary depending on race. The allele frequency of the *28 homozygote has been reported to be approximately 30–40% in Caucasian and African American populations , and 10–13% in Japanese populations [10, 12]. In contrast, the *6 homozygote is rarely detected in Caucasian and African American populations but is specifically found in Asian populations with a frequency of approximately 16% . Similar allele frequencies were found in this study, and, therefore, the subjects are unlikely to deviate from a general Japanese population.
In this study, the median DTG trough concentration of the *6 homozygous patients was approximately 1.7 times higher than that of the subjects carrying both normal alleles. On the other hand, in our previous study on RAL, the median RAL trough concentration of *6 homozygous patients was approximately 10 times higher than that of patients carrying both normal alleles . These observations indicate that the UGT1A1*6 had a smaller effect on the DTG trough concentration than it did on the RAL trough concentration. However, this interpretation is limited by the small number of *6 homozygous patients in both studies. Possible factors for the differential effects include a difference in the mediating metabolic pathway between RAL and DTG. RAL is predominantly metabolized by the UGT1A1 pathway and is not a substrate of CYP enzymes. In contrast, the principal pathway responsible for metabolizing approximately 70% of DTG involves UGT1A1, while CYP3A and other enzymes are involved in the secondary metabolic pathway of the remaining 30% . This secondary metabolic pathway could reduce the effect of UGT1A1*6 on the plasma DTG concentration.
Unexpectedly, age < 40 years was an independent factor associated with high DTG trough concentrations. Drug blood concentrations increase with age. However, a population pharmacokinetic analysis in a previous study also showed the association between older age and decreased DTG plasma trough concentration . Several factors could be considered responsible for this observation, including reduced absorption of DTG in the gastrointestinal tract, increased clearance of DTG due to low levels of plasma albumin resulting in an increased unbound fraction of plasma DTG [8, 24], and altered body fat, which affects the distribution of highly lipid-soluble drugs . However, these reasons seem to be unlikely in this study because these physiological changes are apparent in the elderly [25, 26], and 93% of the subjects in this study were 20–50 years old. Further study is required to clarify the reason for the age effect.
This study demonstrated that plasma DTG trough concentrations were higher in the cases with NP-AEs than in those without NP-AEs. In addition, an association between UGT1A1 gene polymorphism and NP-AEs was suggested. The mechanism of the latter seems to be that reduced-function alleles of UGT1A1 can induce NP-AEs by increasing plasma DTG concentrations. In addition, UGT1A1 gene polymorphism, age, gender, body weight, smoking, and food consumption were reported to affect the pharmacokinetics of DTG . Therefore, it seems reasonable that plasma DTG trough concentration showed a stronger association with NP-AEs than UGT1A1 gene polymorphism. Similar to the effect of DTG in this study, plasma efavirenz (EFV) concentration was reported to be associated with toxicities of the central nervous system . That study suggested that a 1–4 μg/mL range at mid-dosing interval was a suitable target for plasma EFV concentration. However, it was difficult to set a cut-off value for DTG in our study, because we did not include patients who stopped taking DTG within 150 days of initiation. In a retrospective study about cases of DTG discontinuation, the median DTG trough concentration of 12 cases was reported to be 1.72 μg/mL . Considering these observations, more detailed pharmacokinetics about the DTG discontinuation cases may be required to set appropriate cut-off values to predict the risk of NP-AEs after the use of DTG.
This study has some limitations. For example, this was an observational study with a limited number of patients at a single center. Especially, the number of *6 and *28 homozygous and *6/*28 compound heterozygous patients was small. The plasma DTG concentration was evaluated only using trough concentration. Furthermore, our cohort may not be representative of the entire population and the discontinuation rate of DTG in this study may differ from that of a typical setting in Japan. Possible reasons for these differences include that we could not recruit patients who stopped taking DTG before plasma DTG concentrations reached steady state, and that the subjects were selected without use of random sampling techniques. Recently, it was reported that gene polymorphisms of drug transporters influenced plasma RAL peak concentrations . It is likely that gene polymorphisms of drug transporters can affect the pharmacokinetics of DTG as well. Gene polymorphisms CYP3A4*22 and CYP3A5*3 could also alter the pharmacokinetics of DTG, because DTG is partially metabolized by CYP3A and these polymorphisms were reported to change the metabolism of tacrolimus . Further studies should be conducted to examine the pharmacokinetics of DTG and the associated factors in addition to UGT.
We demonstrated that carrying UGT1A1*6, UGT1A1*28, or both alleles, as well as younger age, were factors resulting in high trough concentrations of DTG. In addition, our results suggest that a relationship exists between plasma DTG trough concentrations and the development of selected NP-AEs after the use of DTG. Carrying UGT1A1*6 and/or UGT1A1*28 alleles may also be a risk of developing NP-AEs.
Cytochrome P450 3A
Human immunodeficiency virus
Neuropsychiatric adverse events
Uridine diphosphate (UDP)-glucuronosyltransferase 1A1
Walmsley SL, Antela A, Clumeck N, Duiculescu D, Eberhard A, Gutierrez F, et al. Dolutegravir plus abacavir-lamivudine for the treatment of HIV-1 infection. N Engl J Med. 2013;369:1807–18.
Molina JM, Clotet B, van Lunzen J, Lazzarin A, Cavassini M, Henry K, et al. Once-daily dolutegravir versus darunavir plus ritonavir for treatment-naive adults with HIV-1 infection (FLAMINGO): 96 week results from a randomised, open-label, phase 3b study. Lancet HIV. 2015;2:e127–36.
Raffi F, Jaeger H, Quiros-Roldan E, Albrecht H, Belonosova E, Gatell JM, et al. Once-daily dolutegravir versus twice-daily raltegravir in antiretroviral-naive adults with HIV-1 infection (SPRING-2 study): 96 week results from a randomised, double-blind, non-inferiority trial. Lancet Infect Dis. 2013;13:927–35.
Walmsley S, Baumgarten A, Berenguer J, Felizarta F, Florence E, Khuong-Josses MA, et al. Brief report: Dolutegravir plus Abacavir/Lamivudine for the treatment of HIV-1 infection in antiretroviral therapy-naive patients: week 96 and week 144 results from the SINGLE randomized clinical trial. J Acquir Immune Defic Syndr. 2015;70:515–9.
The Department of Health and Human Services: Guidelines for the use of antiretroviral agents in HIV-1-infected adults and adolescents. 2016 (revised on July 14), http://aidsinfo.nih.gov/.
The European AIDS Clinical Society treatment guidelines version 8.1. 2016 (revised on October), http://www.eacsociety.org/.
Kandel CE, Walmsley SL. Dolutegravir - a review of the pharmacology, efficacy, and safety in the treatment of HIV. Drug Des Devel Ther. 2015;9:3547–55.
Zhang J, Hayes S, Sadler BM, Minto I, Brandt J, Piscitelli S, et al. Population pharmacokinetics of dolutegravir in HIV-infected treatment-naive patients. Br J Clin Pharmacol. 2015;80:502–14.
Calcagno A, Cusato J, D'Avolio A, Bonora S. Genetic polymorphisms affecting the pharmacokinetics of antiretroviral drugs. Clin Pharmacokinet. 2017;56:355–69.
Kaniwa N, Kurose K, Jinno H, Tanaka-Kagawa T, Saito Y, Saeki M, et al. Racial variability in haplotype frequencies of UGT1A1 and glucuronidation activity of a novel single nucleotide polymorphism 686C> T (P229L) found in an African-American. Drug Metab Dispos. 2005;33:458–65.
Ando Y, Fujita K, Sasaki Y, Hasegawa Y. UGT1AI*6 And UGT1A1*27 for individualized irinotecan chemotherapy. Curr Opin Mol Ther. 2007;9:258–62.
Saeki M, Saito Y, Jinno H, Sai K, Ozawa S, Kurose K, et al. Haplotype structures of the UGT1A gene complex in a Japanese population. Pharmacogenomics J. 2006;6:63–75.
Chen S, St Jean P, Borland J, Song I, Yeo AJ, Piscitelli S, et al. Evaluation of the effect of UGT1A1 polymorphisms on dolutegravir pharmacokinetics. Pharmacogenomics. 2014;15:9–16.
de Boer MG, van den Berk GE, van Holten N, Oryszcyn JE, Dorama W, Moha DA, et al. Intolerance of dolutegravir-containing combination antiretroviral therapy regimens in real-life clinical practice. AIDS. 2016;30:2831–4.
Hoffmann C, Welz T, Sabranski M, Kolb M, Wolf E, Stellbrink HJ, et al. Higher rates of neuropsychiatric adverse events leading to dolutegravir discontinuation in women and older patients. HIV Med. 2017;18:56–63.
Bonfanti P, Madeddu G, Gulminetti R, Squillace N, Orofino G, Vitiello P, et al. Discontinuation of treatment and adverse events in an Italian cohort of patients on dolutegravir. AIDS. 2017;31:455–7.
Briand C, Dollfus C, Faye A, Kantor E, Avettand-Fenoel V, Caseris M, et al. Efficacy and tolerance of dolutegravir-based combined ART in perinatally HIV-1-infected adolescents: a French multicentre retrospective study. J Antimicrob Chemother. 2017;72:837–43.
Cattaneo D, Rizzardini G, Gervasoni C. Intolerance of dolutegravir-containing combination antiretroviral therapy: not just a pharmacokinetic drug interaction. AIDS. 2017;31:867–8.
Fettiplace A, Stainsby C, Winston A, Givens N, Puccini S, Vannappagari V, et al. Psychiatric symptoms in patients receiving Dolutegravir. J Acquir Immune Defic Syndr. 2017;74:423–31.
Menard A, Montagnac C, Solas C, Meddeb L, Dhiver C, Tomei C, et al. Neuropsychiatric adverse effects on dolutegravir: an emerging concern in Europe. AIDS. 2017;31:1201–3.
Yagura H, Watanabe D, Ashida M, Kushida H, Hirota K, Ikuma M, et al. Correlation between UGT1A1 polymorphisms and raltegravir plasma trough concentrations in Japanese HIV-1-infected patients. J Infect Chemother. 2015;21:713–7.
Takahashi M, Konishi M, Kudaka Y, Okumura N, Hirano A, Terahata N, et al. A conventional LC-MS method developed for the determination of plasma raltegravir concentrations. Biol Pharm Bull. 2008;31:1601–4.
Castellino S, Moss L, Wagner D, Borland J, Song I, Chen S, et al. Metabolism, excretion, and mass balance of the HIV-1 integrase inhibitor dolutegravir in humans. Antimicrob Agents Chemother. 2013;57:3536–46.
Kobayashi M, Yoshinaga T, Seki T, Wakasa-Morimoto C, Brown KW, Ferris R, et al. In vitro antiretroviral properties of S/GSK1349572, a next-generation HIV integrase inhibitor. Antimicrob Agents Chemother. 2011;55:813–21.
Swift CG, Ewen JM, Clarke P, Stevenson IH. Responsiveness to oral diazepam in the elderly: relationship to total and free plasma concentrations. Br J Clin Pharmacol. 1985;20:111–8.
Greenblatt DJ. Reduced serum albumin concentration in the elderly: a report from the Boston collaborative drug surveillance program. J Am Geriatr Soc. 1979;27:20–2.
Marzolini C, Telenti A, Decosterd LA, Greub G, Biollaz J, Buclin T. Efavirenz plasma levels can predict treatment failure and central nervous system side effects in HIV-1-infected patients. AIDS. 2001;15:71–5.
Tsuchiya K, Hayashida T, Hamada A, Oka S, Gatanaga H. Brief report: high peak level of plasma Raltegravir concentration in patients with ABCB1 and ABCG2 genetic variants. J Acquir Immune Defic Syndr. 2016;72:11–4.
Elens L, Bouamar R, Hesselink DA, Haufroid V, van der Heiden IP, van Gelder T, et al. A new functional CYP3A4 intron 6 polymorphism significantly affects tacrolimus pharmacokinetics in kidney transplant recipients. Clin Chem. 2011;57:1574–83.
We are grateful to all the study participants. We would like to thank Editage (http://www.editage.jp) for English language editing.
This work was supported by the Research Program on HIV/AIDS from the Japan Agency for Medical Research and Development (AMED) to HY and DW, and the Collaborative Research in National Hospital Organization Network to HT.
Availability of data and materials
The data used in this study were obtained under an exclusive data-sharing agreement and are not currently publicly available.
Ethics approval and consent to participate
The study protocol was reviewed and approved by the Institutional Review Board of the Osaka National Hospital (Approval Number: 0838). The procedures were carried out in accordance with relevant guidelines and regulations. All samples were collected after acquiring written consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yagura, H., Watanabe, D., Kushida, H. et al. Impact of UGT1A1 gene polymorphisms on plasma dolutegravir trough concentrations and neuropsychiatric adverse events in Japanese individuals infected with HIV-1. BMC Infect Dis 17, 622 (2017). https://doi.org/10.1186/s12879-017-2717-x