Herpesviridae and Atypical Bacteria Co-Detections in Lower Respiratory Tract Samples of SARS-CoV-2-Positive Patients Admitted to an Intensive Care Unit.

Shortly after the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), cases of viral, bacterial, and fungal coinfections in hospitalized patients became evident. This retrospective study investigates the prevalence of multiple pathogen co-detections in 1472 lower respiratory tract (LRT) samples from 229 SARS-CoV-2-positive patients treated in the largest intensive care unit (ICU) in Slovenia. In addition to SARS-CoV-2, (rt)RT-PCR tests were used to detect cytomegalovirus (CMV), Epstein-Barr virus (EBV), herpes simplex virus 1 (HSV-1), herpes simplex virus 2 (HSV-2), varicella zoster virus (VZV), and atypical bacteria: Chlamydia pneumoniae, Mycoplasma pneumoniae and Legionella pneumophila/spp. At least one co-detection was observed in 89.1% of patients. EBV, HSV-1, and CMV were the most common, with 74.7%, 58.1%, and 38.0% of positive patients, respectively. The median detection time of EBV, HSV-1, and CMV after initial SARS-CoV-2 confirmation was 11 to 20 days. Bronchoalveolar lavage (BAL) and tracheal aspirate (TA) samples showed equivalent performance for the detection of EBV, CMV, and HSV-1 in patients with both available samples. Our results indicate that SARS-CoV-2 infection could be a risk factor for latent herpesvirus reactivation, especially HSV-1, EBV, and CMV. However, additional studies are needed to elucidate the clinical importance of these findings.

Seasonality and Genotype Diversity of Human Rhinoviruses during an Eight-Year Period in Slovenia.

Due to the high socioeconomic burden of rhinoviruses, the development of prevention and treatment strategies is of high importance. Understanding the epidemiological and clinical features of rhinoviruses is essential in order to address these issues. Our study aimed to define the seasonality and molecular epidemiology of rhinoviruses in Slovenia. Over a period of eight years, a total of 20,425 patients from sentinel primary healthcare settings and sentinel hospitals were examined for a panel of respiratory viruses in the national programme for the surveillance of influenza-like illnesses and acute respiratory infections. The patients were from all age groups and had respiratory infections of various severity. Infection with a rhinovirus was confirmed using an RT-rPCR in 1834 patients, and 1480 rhinoviruses were genotyped. The molecular analysis was linked to demographical and meteorological data. We confirmed the year-round circulation of rhinoviruses with clear seasonal cycles, resulting in two seasonal waves with peaks in spring and autumn. High levels of genotype variability and co-circulation were confirmed between and within seasons and were analysed in terms of patient age, the patient source reflecting disease severity, and meteorological factors. Our study provides missing scientific information on the genotype diversity of rhinoviruses in Slovenia. As most previous investigations focused on exclusive segments of the population, such as children or hospitalised patients, and for shorter study periods, our study, with its design, size and length, contributes complementary aspects and new evidence-based knowledge to the regional and global understanding of rhinovirus seasonality and molecular epidemiology.

Enrichment of SARS-CoV-2 sequence from nasopharyngeal swabs whilst identifying the nasal microbiome.

Simultaneously characterising the genomic information of coronaviruses and the underlying nasal microbiome from a single clinical sample would help characterise infection and disease. Metatranscriptomic approaches can be used to sequence SARS-CoV-2 (and other coronaviruses) and identify mRNAs associated with active transcription in the nasal microbiome. However, given the large sequence background, unenriched metatranscriptomic approaches often do not sequence SARS-CoV-2 to sufficient read and coverage depth to obtain a consensus genome, especially with moderate and low viral loads from clinical samples. In this study, various enrichment methods were assessed to detect SARS-CoV-2, identify lineages and define the nasal microbiome. The methods were underpinned by Oxford Nanopore long-read sequencing and variations of sequence independent single primer amplification (SISPA). The utility of the method(s) was also validated on samples from patients infected seasonal coronaviruses. The feasibility of profiling the nasal microbiome using these enrichment methods was explored. The findings shed light on the performance of different enrichment strategies and their applicability in characterising the composition of the nasal microbiome.

Molecular epidemiology of the 2022 monkeypox virus outbreak in Slovenia.

INTRODUCTION: Monkeypox virus (MPXV), typically endemic in West and Central Africa, has raised global concern due to the recent outbreak in several non-endemic countries with human-to-human transmission. Here we present a comprehensive analysis of MPXV genomes from Slovenia. METHODS: Two real-time polymerase chain reaction (RT-PCR) assays for Orthopoxvirus (OPV) and MPXV genes were used for laboratory confirmation of mpox. Complete MPXV genomic sequences were obtained using nanopore long reads and Illumina technology. Phylogenetic analyses compared the Slovenian MPXV sequences with the global sequences. RESULTS: A total of 49 laboratory-confirmed mpox cases were diagnosed in Slovenia in 2022, mainly affecting males under 40. In 48 cases, a complete genome sequence was obtained and phylogenetic analysis revealed five distinct lineages (B.1, B.1.14, B.1.2, B.1.3, and A.2.1), with B.1 and B.1.3 dominating, suggesting multiple introductions into Slovenia. Genome analysis revealed significant divergence from the reference MPXV-M5312_HM12_Rivers. CONCLUSIONS: The genetic diversity observed in the Slovenian MPXV sequences sheds light on the complex dynamics of the 2022 mpox outbreak and highlights the need for further research to understand the impact of mutations on MPXV functional characteristics and their role in the evolution and diversification of current lineages.

Changes in HRSV Epidemiology but Not Circulating Variants in Hospitalized Children due to the Emergence of SARS-CoV-2.

This study assesses the circulation of human respiratory syncytial virus (HRSV) genotypes before, during, and toward the end of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in children and determines the influence of the pandemic on HRSV circulation patterns and evolution. Phylogenetic analysis of the hypervariable glycoprotein G gene was performed on 221/261 (84.7%) HRSV-positive samples and shows two separated clusters, one belonging to HRSV-A (129/221) and another to HRSV-B (92/221). All Slovenian HRSV-A strains contained the 72-nucleotide-long duplicated region in the attachment glycoprotein G gene and were classified as lineage GA2.3.5. All Slovenian HRSV-B strains similarly contained a 60-nucleotide-long duplicated region in the attachment glycoprotein G gene and were classified as lineage GB5.0.5a. During the 3-year period (2018-2021) covered by the study, no significant differences were observed within strains detected before the SARS-CoV-2 pandemic, during it, and after the implementation of nonpharmaceutical preventive measures. Slovenian HRSV-A strains seem to be more diverse than HRSV-B strains. Therefore, further whole-genome investigations would be required for better monitoring of the long-term impact of SARS-CoV-2 endemic circulation and the formation of new HRSV lineages and epidemiological patterns.

A simplified nasopharyngeal swab collection procedure for minimizing patient discomfort while retaining sample quality.

A nasopharyngeal swab (NPS) is the most frequently collected sample type when molecular diagnosis of respiratory viruses, including SARS CoV-2, is required. An optimal collection technique would provide sufficient sample quality for the diagnostic process and would minimize the discomfort felt by the patient. This study compares a simplified NPS collection procedure with only one rotation of the swab to a more standard procedure with five rotations. Swabs were collected from 76 healthy volunteers by the same healthcare professional on 2 consecutive days at a similar hour to minimize variability. The number of Ubiquitin C copy number per sample was measured by real-time quantitative PCR and patient discomfort was assessed by questionnaire. No statistically significant difference (p = 0.15) was observed in the Ubiquitin C copy number per sample between a NPS collected with one rotation (5.2 ± 0.6 log UBC number copies/sample) or five rotations (5.3 ± 0.5 log UBC number copies/sample). However, a statistically significant difference was observed in discomfort between these two procedures, the second being much more uncomfortable. Additional analysis of the results showed a weak correlation between discomfort and the number of human cells recovered (Spearman's rho = 0.202) and greater discomfort in younger people. The results of this study show that a NPS collected with one slow rotation has the same quality as a NPS collected with five rotations. However, the collection time is shorter and, most importantly, less unpleasant for patients.

Characterisation of the Antibody Response in Sinopharm (BBIBP-CorV) Recipients and COVID-19 Convalescent Sera from the Republic of Moldova.

The early availability of effective vaccines against SARS-CoV-2, the aetiologic cause of COVID-19, has been at the cornerstone of the global recovery from the pandemic. This study aimed to assess the antispike RBD IgG antibody titres and neutralisation potential of COVID-19 convalescent plasma and the sera of Moldovan adults vaccinated with the Sinopharm BBIBP-CorV vaccine. An IgG ELISA with recombinant SARS-CoV-2 spike RBD and two pseudovirus-based neutralisation assays have been developed to evaluate neutralising antibodies against SARS-CoV-2 in biosafety level 2 containment facilities. A significant moderate correlation was observed between IgG titres and the overall neutralising levels for each neutralisation assay (ρ = 0.64, p < 0.001; ρ = 0.52, p < 0.001). A separate analysis of convalescent and vaccinated individuals showed a higher correlation of neutralising and IgG titres in convalescent individuals (ρ = 0.68, p < 0.001, ρ = 0.45, p < 0.001) compared with vaccinated individuals (ρ = 0.58, p < 0.001; ρ = 0.53, p < 0.001). It can be concluded that individuals who recovered from infection developed higher levels of antispike RBD IgG antibodies. In comparison, the Sinopharm-vaccinated individuals produced higher levels of neutralising antibodies than convalescent plasma.

Molecular Typing of Mastadenoviruses in Simultaneously Collected Nasopharyngeal Swabs and Stool Samples from Children Hospitalized for Acute Bronchiolitis, Acute Gastroenteritis, and Febrile Seizures.

This study determines and compares the frequency of human mastadenovirus (HAdV) presence in children with acute bronchiolitis (AB), acute gastroenteritis (AGE), and febrile seizures (FS), ascertains types of HAdVs associated with each individual syndrome and contrasts the findings with a control group of children. The presence of HAdVs was ascertained in simultaneously collected nasopharyngeal (NP) swabs and stool samples amplifying the hexon gene by RT-PCR; these were sequenced to determine the types of HAdVs. HAdVs were grouped into eight different genotypes. Of these, three (F40, F41, and A31) were found solely in stool samples, whereas the others (B3, C1, C2, C5, and C6) were found in both stool samples and NP swabs. The most common genotypes in NP swabs were C2 (found in children with AGE and FS) and C1 (only in children with FS), whereas in stool samples genotypes F41 (in children with AGE) and C2 (in children with AGE and FS) prevailed, and C2 was simultaneously present in both samples. HAdVs were more often detected in stool samples than in NP swabs in patients (with the highest estimated viral load in stool samples in children with AB and AGE) and healthy controls and were more common in NP swabs in children with AGE than in children with AB. In most patients, the characterized genotypes in NP swabs and stool samples were in concordance.

Postnatal CMV Infection and Antiviral Treatment in Extremely Premature Infants: A 12-Year Retrospective Analysis.

BACKGROUND: The impact and outcomes of postnatal cytomegalovirus (CMV) infection are not entirely clear. We aimed to determine the associations between treatment outcomes of postnatal CMV infection and its antiviral treatment. METHODS: Retrospective study in a tertiary center. Infants of < 29 weeks gestational age who were tested for postnatal CMV infection were included. CMV-infected infants were compared to uninfected infants (control group). CMV-infected infants were either treated with ganciclovir and/or valganciclovir (CMVPT group) or not (CMVPNT group). Demographic, clinical, laboratory, treatment, and outcome data were collected. Primary outcomes were the length of stay, death before discharge and hearing impairment, cognitive and motor development as assessed by the Denver Developmental Screening Test II, and neurologic impairment at the corrected age of 1.5-2 years. RESULTS: We included 103 extremely premature infants. The Median (interquartile range [IQR]) length of stay was 94 (69-112) days in control, 85 (70-102) days in CMVPNT, and 100 (88-137) days in the CMVPT group. Mortality before discharge was 6% in control, 3.8% in CMVPNT, and 3.7% in the CMVPT group. Normal hearing at follow-up was found in 30/37 infants in control (81.1%), 13/13 infants in CMVPNT (100%), and 17/20 infants in the CMVPT group (85%). Denver Developmental Screening Test II results did not differ among the three groups. Neurologic impairment was found in 21/37 infants (56.8%) in control, 9/13 infants in CMVPNT (69.2%), and 14/20 infants in CMVPT group (70%). CONCLUSIONS: The associations between antiviral treatment of postnatal CMV infection and better treatment outcomes were nonsignificant.

Evaluation of Two Broadly Used Commercial Methods for Detection of Respiratory Viruses with a Recently Added New Target for Detection of SARS-CoV-2.

The clinical symptoms caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are nonspecific and can be associated with most other respiratory viruses that cause acute respiratory tract infections (ARI). Because the clinical differentiation of COVID-19 patients from those with other respiratory viruses is difficult, the evaluation of automated methods to detect important respiratory viruses together with SARS-CoV-2 seems necessary. Therefore, this study compares two molecular assays for the detection of respiratory viruses, including SARS-CoV-2: the Respiratory Viruses 16-Well Assay (AusDiagnostics, Pty Ltd., Mascot, Australia) and the Allplex™ RV Essential Assay coupled with the Allplex™-nCoV Assay (Seegene Inc., Seoul, Korea). The two methods (AusDiagnostics and AlplexTM-nCoV Assay SARS-CoV-2) had 98.6% agreement with the reference method, cobas 6800, for the detection of SARS-CoV-2. Agreement between the AusDiagnostics assay and the AlplexTM RV Essential Assay for the detection of seven respiratory viruses was 99%. In our experience, the Respiratory Viruses 16-Well Assay proved to be the most valuable and useful medium-throughput method for simultaneous detection of important respiratory viruses and SARS-CoV-2. The main advantages of the method are high specificity for all targets included and their simultaneous detection and medium throughput with the option of having multiple instruments provide a constant run.

Performance of nasopharyngeal swab and saliva in detecting Delta and Omicron SARS-CoV-2 variants.

A prospective cohort study was conducted during the Delta and Omicron severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) epidemic waves from paired nasopharyngeal swab (NPS or NP swab) and saliva samples taken from 624 participants. The study aimed to assess if any differences among participants from both waves could be observed and if any difference in molecular diagnostic performance could be observed among the two sample types. Samples were transported immediately to the laboratory to ensure the highest possible sample quality without any freezing and thawing steps before processing. Nucleic acids from saliva and NPS were prospectively extracted and SARS-CoV-2 was detected using a real-time reverse-transcription polymerase chain reaction. All observed results were statistically analyzed. Although the results obtained with NP and saliva agreed overall, higher viral loads were observed in NP swabs regardless of the day of specimen collection in both SARS-CoV-2 epidemic waves. No significant difference could be observed between the two epidemic waves characterized by Delta or Omicron SARS-CoV-2. To note, Delta infection resulted in higher viral loads both in NP and saliva and more symptoms, including rhinorrhea, cough, and dyspnea, whereas Omicron wave patients more frequently reported sore throat. An increase in the mean log RNA of SARS-CoV-2 was observed with the number of expressed symptoms in both waves, however, the difference was not significant. Data confirmed that results from saliva were concordant with those from NP swabs, although saliva proved to be a challenging sample with frequent inhibitions that required substantial retesting.

Genetic Characterisation and Comparison of Three Human Coronaviruses (HKU1, OC43, 229E) from Patients and Bovine Coronavirus (BCoV) from Cattle with Respiratory Disease in Slovenia.

Coronaviruses (CoV) are widely distributed pathogens of human and animals and can cause mild or severe respiratory and gastrointestinal disease. Antigenic and genetic similarity of some CoVs within the Betacoronavirus genus is evident. Therefore, for the first time in Slovenia, we investigated the genetic diversity of partial 390-nucleotides of RNA-dependent-RNA polymerase gene (RdRp) for 66 human (HCoV) and 24 bovine CoV (BCoV) positive samples, collected between 2010 and 2016 from human patients and cattle with respiratory disease. The characterized CoV strains belong to four different clusters, in three separate human clusters HCoV-HKU1 (n = 34), HCoV-OC43 (n = 31) and HCoV 229E (n = 1) and bovine grouping only as BCoVs (n = 24). BCoVs from cattle and HCoV-OC43 were genetically the most closely related and share 96.4-97.1% nucleotide and 96.9-98.5% amino acid identity.

Seroprevalence of severe acute respiratory syndrome coronavirus 2 in Slovenia: results of two rounds of a nationwide population study on a probability-based sample, challenges and lessons learned.

OBJECTIVES: Seroprevalence surveys provide crucial information on cumulative severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exposure. This Slovenian nationwide population study is the first longitudinal 6-month serosurvey using probability-based samples across all age categories. METHODS: Each participant supplied two blood samples: 1316 samples in April 2020 (first round) and 1211 in October/November 2020 (second round). The first-round sera were tested using Euroimmun Anti-SARS-CoV-2 ELISA IgG (ELISA) and, because of uncertain estimates, were retested using Elecsys Anti-SARS-CoV-2 (Elecsys-N) and Elecsys Anti-SARS-CoV-2 S (Elecsys-S). The second-round sera were concomitantly tested using Elecsys-N/Elecsys-S. RESULTS: The populations of both rounds matched the overall population (n = 3000), with minor settlement type and age differences. The first-round seroprevalence corrected for the ELISA manufacturer's specificity was 2.78% (95% highest density interval [HDI] 1.81%-3.80%), corrected using pooled ELISA specificity calculated from published data 0.93% (95% CI 0.00%-2.65%), and based on Elecsys-N/Elecsys-S results 0.87% (95% HDI 0.40%-1.38%). The second-round unadjusted lower limit of seroprevalence on 11 November 2020 was 4.06% (95% HDI 2.97%-5.16%) and on 3 October 2020, unadjusted upper limit was 4.29% (95% HDI 3.18%-5.47%). CONCLUSIONS: SARS-CoV-2 seroprevalence in Slovenia increased four-fold from late April to October/November 2020, mainly due to a devastating second wave. Significant logistic/methodological challenges accompanied both rounds. The main lessons learned were a need for caution when relying on manufacturer-generated assay evaluation data, the importance of multiple manufacturer-independent assay performance assessments, the need for concomitant use of highly-specific serological assays targeting different SARS-CoV-2 proteins in serosurveys conducted in low-prevalence settings or during epidemic exponential growth and the usefulness of a Bayesian approach for overcoming complex methodological challenges.

Characterization of Tularemia Cases in Slovenia with Multiple-Locus Variable-Number Tandem Repeat Analysis.

Francisella tularensis is the etiologic agent of tularemia, a bacterial zoonotic disease. The genome of F. tularensis shows a recent evolutionary change, especially in reservoirs. Variable number of tandem repeats (VNTR) is described as a high-speed molecular clock and can thus be used as a high-resolution typing system. The main objective of our study was to investigate the molecular diversity of F. tularensis strains and reveal possible sources of infection. Using real-time PCR targeting the ISFtu2 region, we successfully amplified targeted DNA in 13/31 Slovenian patients with a clinical diagnosis of tularemia, and with PCR targeting the fopA gene, we obtained 11/13 PCR products. Sequencing revealed that all samples were identified as F. tularensis subsp. holarctica. We successfully obtained one F. tularensis isolate from a lymph node aspirate by culture on chocolate agar. Our isolate was clustered into major clade B12 (subclade B43). We optimized VNTR typing to be used directly on clinical samples. Multiple-locus VNTR analysis (MLVA) revealed five unique MLVA types; 45.5% samples had the same MLVA type, another 27.3% shared a different MLVA type, and each of the remaining had a unique MLVA type. Most samples differed at only two VNTR markers (Ft-M03 and Ft-M06). Additionally, we investigated samples from small mammals (n = 532) and Ixodes ricinus ticks (n = 232) captured in the same geographical area in which patients with tularemia were found. No F. tularensis DNA was detected in samples of small mammals or I. ricinus ticks. The diversity of MLVA types in Slovenia was high, despite the small region, but most of the samples from the same region shared the same MLVA type. Our results suggest that MLVA is a useful tool for quick molecular characterization of F. tularensis directly from patient samples, especially when investigating geographically localized outbreaks.

Comparison of Lymphocyte Populations in Patients With Dobrava or Puumala orthohantavirus Infection.

Hemorrhagic fever with renal syndrome (HFRS), caused by Dobrava (DOBV) and Puumala (PUUV) orthohantaviruses, is an endemic disease in Slovenia. DOBV is mainly responsible for a more severe disease, whereas PUUV usually causes a milder form. Therefore, the aim of our study was to determine whether any differences in lymphocyte population in patients infected with these two viruses exist. Mononuclear cells from peripheral blood (PBMCs) were isolated from DOBV or PUUV infected patients and different lymphocyte subpopulations were analyzed with flow cytometry. Decreased concentrations of lymphocyte subpopulation were observed in DOBV and in PUUV infected patients compared with a healthy control, which was especially evident in DOBV infected patients. The lower values of T cells are likely due to the extravasation of the activated cells from the circulation to the infected tissue. Higher percentage of NK cells were detected in DOBV infected patients in comparison to PUUV infected patients, which could be associated with a more severe HFRS caused by DOBV. PUUV infected patients had a significantly higher concentration of activated T cell subsets, expressing markers CD25, CD69, and HLA-DR in comparison to DOBV infected patients. Higher activation of T cell subsets in PUUV infected patients could be a contributor to a milder HFRS. Further studies are necessary to elucidate the relation between the protective and the harmful role of activated lymphocytes subsets in HFRS pathogenesis.

Correction to: The absence of the drhm gene is not a marker for human-pathogenicity in European Anaplasma phagocytophilum strains.

An amendment to this paper has been published and can be accessed via the original article.

Five-Day vs 10-Day Postexposure Chemoprophylaxis With Oseltamivir to Prevent Hospital Transmission of Influenza: A Noninferiority Randomized Open-Label Study.

BACKGROUND: The efficacy and optimal duration of postexposure influenza prophylaxis with oseltamivir are undetermined in hospital settings, where immediate separation from index cases is not feasible. METHODS: In an open-label noninferiority randomized clinical trial in a single-center university hospital, the efficacy of 5-day vs 10-day postexposure prophylaxis with oseltamivir was compared in adult patients exposed to influenza who could not be immediately separated from index influenza cases. Influenza incidence was assessed for 10 days after discontinuing prophylaxis. RESULTS: Among 222 exposed patients (median age, 75 years; male 119; median Charlson Comorbidity Index, 5), 110 patients were assigned to 5 days of postexposure prophylaxis with oseltamivir, and 112 patients were assigned to the 10-day group. The median duration of exposure to influenza (interquartile range) was 2 (1-3) days. In the intention-to-treat analysis, the incidence of influenza was 2/110 (1.8%) in the 5-day group and 0/112 (0%) in the 10-day group (difference, 1.8 percentage points; 1-sided 95% CI, -1 to 4.9 percentage points; P = .77). CONCLUSIONS: For patients exposed to influenza in a hospital setting and who were not immediately separated from index cases, postexposure prophylaxis with oseltamivir resulted in low incidence of nosocomial influenza transmission. Five-day postexposure prophylaxis was noninferior to 10-day regimen. CLINICALTRIALSGOV REGISTRATION: NCT03899571.

Outbreak of Human Metapneumovirus Infection in Zoo, Slovenia.

We report a case of human metapneumovirus infection that spread from humans to chimpanzees and back to humans. Bronchopneumonia developed in 4 of 6 members of a chimpanzee family, and 2 subsequently died. The chimpanzees' keeper also became ill. Sequencing showed 100% identity between virus sequences from chimpanzees and the keeper.

The absence of the drhm gene is not a marker for human-pathogenicity in European Anaplasma phagocytophilum strains.

BACKGROUND: Anaplasma phagocytophilum is a Gram-negative obligate intracellular bacterium that replicates in neutrophil granulocytes. It is transmitted by ticks of the Ixodes ricinus complex and causes febrile illness in humans and animals. The geographical distribution of A. phagocytophilum spans the Americas, Europe, Africa and Asia. However, human disease predominantly occurs in North America but is infrequently reported from Europe and Asia. In North American strains, the absence of the drhm gene has been proposed as marker for pathogenicity in humans whereas no information on the presence or absence of the drhm gene was available for A. phagocytophilum strains circulating in Europe. Therefore, we tested 511 European and 21 North American strains for the presence of drhm and compared the results to two other typing methods: multilocus sequence typing (MLST) and ankA-based typing. RESULTS: Altogether, 99% (478/484) of the analyzable European and 19% (4/21) of the North American samples from different hosts were drhm-positive. Regarding the strains from human granulocytic anaplasmosis cases, 100% (35/35) of European origin were drhm-positive and 100% (14/14) of North American origin were drhm-negative. Human strains from North America and Europe were both part of MLST cluster 1. North American strains from humans belonged to ankA gene clusters 11 and 12 whereas European strains from humans were found in ankA gene cluster 1. However, the North American ankA gene clusters 11 and 12 were highly identical at the nucleotide level to the European cluster 1 with 97.4% and 95.2% of identity, respectively. CONCLUSIONS: The absence of the drhm gene in A. phagocytophilum does not seem to be associated with pathogenicity for humans per se, because all 35 European strains of human origin were drhm-positive. The epidemiological differences between North America and Europe concerning the incidence of human A. phagocytophilum infection are not explained by strain divergence based on MLST and ankA gene-based typing.