Efficient SARS-CoV-2 variant detection and monitoring with Spike Screen next-generation sequencing.

The emergence and rapid spread of SARS-CoV-2 prompted the global community to identify innovative approaches to diagnose infection and sequence the viral genome because at several points in the pandemic positive case numbers exceeded the laboratory capacity to characterize sufficient samples to adequately respond to the spread of emerging variants. From week 10, 2020, to week 13, 2023, Slovenian routine complete genome sequencing (CGS) surveillance network yielded 41 537 complete genomes and revealed a typical molecular epidemiology with early lineages gradually being replaced by Alpha, Delta, and finally Omicron. We developed a targeted next-generation sequencing based variant surveillance strategy dubbed Spike Screen through sample pooling and selective SARS-CoV-2 spike gene amplification in conjunction with CGS of individual cases to increase throughput and cost-effectiveness. Spike Screen identifies variant of concern (VOC) and variant of interest (VOI) signature mutations, analyses their frequencies in sample pools, and calculates the number of VOCs/VOIs at the population level. The strategy was successfully applied for detection of specific VOC/VOI mutations prior to their confirmation by CGS. Spike Screen complemented CGS efforts with an additional 22 897 samples sequenced in two time periods: between week 42, 2020, and week 24, 2021, and between week 37, 2021, and week 2, 2022. The results showed that Spike Screen can be applied to monitor VOC/VOI mutations among large volumes of samples in settings with limited sequencing capacity through reliable and rapid detection of novel variants at the population level and can serve as a basis for public health policy planning.

Influenza A, Influenza B, human respiratory syncytial virus and SARSCoV-2 molecular diagnostics and epidemiology in the post COVID-19 era.

BACKGROUND: The concurrent circulation of SARS-CoV-2 with other respiratory viruses is unstoppable and represents a new diagnostic reality for clinicians and clinical microbiology laboratories. Multiplexed molecular testing on automated platforms that focus on the simultaneous detection of multiple respiratory viruses in a single tube is a useful approach for current and future diagnosis of respiratory infections in the clinical setting. METHODS: Two time periods were included in the study: from February to April 2022, an early 2022 period, during the gradual lifting of COVID-19 prevention measures in the country, and from October 2022 to April 2023, the 2022/23 respiratory infections season. We analysed a total of 1,918 samples in the first period and 18,131 respiratory samples in the second period using a multiplex molecular assay for the simultaneous detection of Influenza A (Flu-A), Influenza B (Flu-B), Human Respiratory Syncytial Virus (HRSV) and SARS-CoV-2. RESULTS: The results from early 2022 showed a strong dominance of SARS-CoV-2 infections with 1,267/1,918 (66.1%) cases. Flu-A was detected in 30/1,918 (1.6%) samples, HRSV in 14/1,918 (0.7%) samples, and Flu-B in 2/1,918 (0.1%) samples. Flu-A/SARS-CoV-2 co-detections were observed in 11/1,267 (0.9%) samples, and HRSV/SARS-CoV-2 co-detection in 5/1,267 (0.4%) samples. During the 2022/23 winter respiratory season, SARS-CoV-2 was detected in 1,738/18,131 (9.6%), Flu-A in 628/18,131 (3.5%), Flu-B in 106/18,131 (0.6%), and HRSV in 505/18,131 (2.8%) samples. Interestingly, co-detections were present to a similar extent as in early 2022. CONCLUSION: The results show that the multiplex molecular approach is a valuable tool for the simultaneous laboratory diagnosis of SARS-CoV-2, Flu-A/B, and HRSV in hospitalized and outpatients. Infections with Flu-A/B, and HRSV occurred shortly after the COVID-19 control measures were lifted, so a strong reoccurrence of various respiratory infections and co-detections in the post COVID-19 period was to be expected.

Sequencing analysis of SARS-CoV-2 cases in Slovenian long-term care facilities to support outbreak control.

Introduction: Residents of long-term care facilities (LTCFs) are at high risk of morbidity and mortality due to COVID-19, especially when new variants of concern (VOC) emerge. To provide intradisciplinary data in order to tailor public health interventions during future epidemics, available epidemiologic and genomic data from Slovenian LTCFs during the initial phases of the COVID-19 pandemic was analyzed. Methods: The first part of the study included SARS-CoV-2 reverse-transcription Real-Time PCR (rtRT-PCR) positive LTCF residents, from 21 facilities with COVID-19 outbreaks occurring in October 2020. The second part of the study included SARS-CoV-2 rtRT-PCR positive LTCF residents and staff between January and April 2021, when VOC Alpha emerged in Slovenia. Next-generation sequencing (NGS) was used to acquire SARS-CoV-2 genomes, and lineage determination. In-depth phylogenetic and mutational profile analysis were performed and coupled with available field epidemiological data to assess the dynamics of SARS-CoV-2 introduction and transmission. Results: 370/498 SARS-CoV-2 positive residents as well as 558/699 SARS-CoV-2 positive residents and 301/358 staff were successfully sequenced in the first and second part of the study, respectively. In October 2020, COVID-19 outbreaks in the 21 LTCFs were caused by intra-facility transmission as well as multiple independent SARS-CoV-2 introductions. The Alpha variant was confirmed in the first LTCF resident approximately 1.5 months after the first Alpha case was identified in Slovenia. The data also showed a slower replacement of existing variants by Alpha in residents compared to staff and the general population. Discussion: Multiple SARS CoV-2 introductions as well as intra-facility spreading impacted disease transmission in Slovenian LTCFs. Timely implementation of control measures aimed at limiting new introductions while controlling in-facility transmission are of paramount importance, especially as new VOCs emerge. Sequencing, in conjunction with epidemiological data, can facilitate the determination of the need for future improvements in control measures to protect LTCF residents from COVID-19 or other respiratory infections.

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.

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.

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.

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.

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.

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.

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.