Continued circulation of mpox: an epidemiological and phylogenetic assessment, European Region, 2023 to 2024.

During the summer of 2023, the European Region experienced a limited resurgence of mpox cases following the substantial outbreak in 2022. This increase was characterised by asynchronous and bimodal increases, with countries experiencing peaks at different times. The demographic profile of cases during the resurgence was largely consistent with those reported previously. All available sequences from the European Region belonged to clade IIb. Sustained efforts are crucial to control and eventually eliminate mpox in the European Region.

Low prevalence of influenza viruses and predominance of A(H3N2) virus with respect to SARS-CoV-2 during the 2021-2022 season in Bulgaria.

Social distancing, mask-wearing, and travel restrictions during the COVID-19 pandemic have significantly impacted the spread of influenza viruses. The objectives of this study were to analyze the pattern of influenza virus circulation with respect to that of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Bulgaria during the 2021-2022 season and to perform a phylogenetic/molecular analysis of the hemagglutinin (HA) and neuraminidase (NA) sequences of representative influenza strains. Influenza infection was confirmed using real-time reverse transcription polymerase chain reaction in 93 (4.2%) of the 2193 patients with acute respiratory illness tested wherein all detected viruses were subtyped as A(H3N2). SARS-CoV-2 was identified in 377 (24.3%) of the 1552 patients tested. Significant differences in the incidence of influenza viruses and SARS-CoV-2 were found between individual age groups, outpatients/inpatients, and in the seasonal distribution of cases. Two cases of coinfections were identified. In hospitalized patients, the Ct values of influenza viruses at admission were lower in adults aged ≥65 years (indicating higher viral load) than in children aged 0-14 years (p < 0.05). In SARS-CoV-2-positive inpatients, this association was not statistically significant. HA genes of all A(H3N2) viruses analyzed belonged to subclade 3C.2a1b.2a. The sequenced viruses carried 11 substitutions in HA and 5 in NA, in comparison to the vaccine virus A/Cambodia/e0826360/2020, including several substitutions in the HA antigenic sites B and C. This study revealed extensive changes in the typical epidemiology of influenza infection, including a dramatic reduction in the number of cases, diminished genetic diversity of circulating viruses, changes in age, and seasonal distribution of cases.

Induction of humoral and cellular immune responses to COVID-19 mRNA and vector vaccines: A prospective cohort study in Bulgarian healthcare workers.

Installing efficient protective immunity by anti-SARS-CoV-2 vaccines is the only current means to overcome coronavirus disease 2019 pandemics. The cellular and humoral immune responses induced with an messenger RNA (mRNA) (BNT162b2) or with a vector (ChAdOx1nCoV-19) vaccine among Bulgarian healthcare workers (n = 123, aged 23-71 years) were studied in the course of 16 weeks after priming. Receptor-binding domain (RBD)-blocking Abs and SARS-CoV-2 RBD immunoglobulin A  (IgA) were evaluated in parallel with interferon gamma (IFNγ)-producing virus-specific T cells. Both vaccines induced RBD-blocking Abs in 100% of the participants after complete immunization while the levels of protection after a single dose largely varied (22%-98%). Advanced age had a negative impact on the level and longevity of virus-neutralizing activity induced by one dose mRNA, but not by the vector vaccine. RBD-binding IgA was detected in 100% of tested donors from the mRNA vaccine cohort, and in 67% of tested from the vector vaccine cohort, at least 1 month after completed immunization. One month after completing mRNA immunization, the number of IFNγ-producing T cells correlated significantly with the levels of RBD-specific IgA and virus-neutralizing activity induced after priming. Enumeration of circulating virus-specific IFNγ+ T cells is not recommended for evaluation of protective immunity as their detection may require longer stimulation beyond the firstmonth postimmunization. In conclusion, BNT162B2 and ChAdOx1nCoV-19 induced potent and comparable humoral and cellular anti-SARS-CoV-2 immune responses, peaking between 10 and 30 days after complete immunization. A single dose of any vaccine did not induce adequate protection in a great part of donors, making the shorter interval between mRNA vaccine doses preferable in the settings of increased risk of infection.

Initial introduction and spread of the SARS-CoV-2 AY.4.2.1 Delta variant in Bulgaria, a genomic insight.

The evolution of the emerging SARS-CoV-2 variants carrying mutations in the spike protein raises concerns about the possibility of accelerated transmission in the ever-evolving COVID-19 pandemic worldwide. AY.4.2, a sublineage of the Delta variant, was considered a variant under investigation (VUI) and also gained the nickname "Delta Plus," due to its extra mutations, Y145H and A222V. In this study, using genomic epidemiology, we provide the first insights into the introduction of AY.4.2 in Bulgaria and the AY.4.2.1 sublineage that found larger dissemination only in Bulgaria and the United Kingdom.

Specialist laboratory networks as preparedness and response tool - the Emerging Viral Diseases-Expert Laboratory Network and the Chikungunya outbreak, Thailand, 2019.

We illustrate the potential for specialist laboratory networks to be used as preparedness and response tool through rapid collection and sharing of data. Here, the Emerging Viral Diseases-Expert Laboratory Network (EVD-LabNet) and a laboratory assessment of chikungunya virus (CHIKV) in returning European travellers related to an ongoing outbreak in Thailand was used for this purpose. EVD-LabNet rapidly collected data on laboratory requests, diagnosed CHIKV imported cases and sequences generated, and shared among its members and with the European Centre for Disease Prevention and Control. Data across the network showed an increase in CHIKV imported cases during 1 October 2018-30 April 2019 vs the same period in 2018 (172 vs 50), particularly an increase in cases known to be related to travel to Thailand (72 vs 1). Moreover, EVD-LabNet showed that strains were imported from Thailand that cluster with strains of the ECSA-IOL E1 A226 variant emerging in Pakistan in 2016 and involved in the 2017 outbreaks in Italy. CHIKV diagnostic requests increased by 23.6% between the two periods. The impact of using EVD-LabNet or similar networks as preparedness and response tool could be improved by standardisation of the collection, quality and mining of data in routine laboratory management systems.

Laboratory management of Crimean-Congo haemorrhagic fever virus infections: perspectives from two European networks.

BackgroundCrimean-Congo haemorrhagic fever virus (CCHFV) is considered an emerging infectious disease threat in the European Union. Since 2000, the incidence and geographic range of confirmed CCHF cases have markedly increased, following changes in the distribution of its main vector, Hyalomma ticks.AimsTo review scientific literature and collect experts' opinion to analyse relevant aspects of the laboratory management of human CCHF cases and any exposed contacts, as well as identify areas for advancement of international collaborative preparedness and laboratory response plans.MethodsWe conducted a literature review on CCHF molecular diagnostics through an online search. Further, we obtained expert opinions on the key laboratory aspects of CCHF diagnosis. Consulted experts were members of two European projects, EMERGE (Efficient response to highly dangerous and emerging pathogens at EU level) and EVD-LabNet (Emerging Viral Diseases-Expert Laboratory Network).ResultsConsensus was reached on relevant and controversial aspects of CCHF disease with implications for laboratory management of human CCHF cases, including biosafety, diagnostic algorithm and advice to improve lab capabilities. Knowledge on the diffusion of CCHF can be obtained by promoting syndromic approach to infectious diseases diagnosis and by including CCHFV infection in the diagnostic algorithm of severe fevers of unknown origin.ConclusionNo effective vaccine and/or therapeutics are available at present so outbreak response relies on rapid identification and appropriate infection control measures. Frontline hospitals and reference laboratories have a crucial role in the response to a CCHF outbreak, which should integrate laboratory, clinical and public health responses.

Improving preparedness to respond to cross-border hepatitis A outbreaks in the European Union/European Economic Area: towards comparable sequencing of hepatitis A virus.

IntroductionSequence-based typing of hepatitis A virus (HAV) is important for outbreak detection, investigation and surveillance. In 2013, sequencing was central to resolving a large European Union (EU)-wide outbreak related to frozen berries. However, as the sequenced HAV genome regions were only partly comparable between countries, results were not always conclusive.AimThe objective was to gather information on HAV surveillance and sequencing in EU/European Economic Area (EEA) countries to find ways to harmonise their procedures, for improvement of cross-border outbreak responses.MethodsIn 2014, the European Centre for Disease Prevention and Control (ECDC) conducted a survey on HAV surveillance practices in EU/EEA countries. The survey enquired whether a referral system for confirming primary diagnostics of hepatitis A existed as well as a central collection/storage of hepatitis A cases' samples for typing. Questions on HAV sequencing procedures were also asked. Based on the results, an expert consultation proposed harmonised procedures for cross-border outbreak response, in particular regarding sequencing. In 2016, a follow-up survey assessed uptake of suggested methods.ResultsOf 31 EU/EEA countries, 23 (2014) and 27 (2016) participated. Numbers of countries with central collection and storage of HAV positive samples and of those performing sequencing increased from 12 to 15 and 12 to 14 respectively in 2016, with all countries typing an overlapping fragment of 218 nt. However, variation existed in the sequenced genomic regions and their lengths.ConclusionsWhile HAV sequences in EU/EEA countries are comparable for surveillance, collaboration in sharing and comparing these can be further strengthened.

WEST NILE VIRUS INFECTION WITH NEUROLOGICAL DISORDERS: A CASE REPORT AND A BRIEF REVIEW OF THE SITUATION IN BULGARIA.

A case of a 66-year-old man with West Nile neuroinvassive disease manifested with fever, weakness, fatigue, consciousness disorders and underlying diabetes mellitus type 2 and cardiovascular diseases is presented. Laboratory data showed elevated erythrocyte sedimentation rate and fibrinogen. Serological tests revealed West Nile virus specific antibodies of class IgM and IgG in serum. West Nile virus RNA was detected in urine sample. Supportive therapy was applied.

High seroprevalence for Crimean-Congo haemorrhagic fever virus in ruminants in the absence of reported human cases in many regions of Bulgaria.

Crimean-Congo haemorrhagic fever (CCHF) is a severe tick-borne zoonotic disease in humans caused by CCHF virus. It has been observed in Bulgaria since 1952 and over the years more than 1600 cases have been reported in the country. Close contact with viraemic livestock was shown as one of the main causes of the infection. Detection of CCHF virus specific antibodies in livestock can be used as an indicator for virus circulation and risk assessment. CCHF seroprevalence was investigated in 843 cattle, 88 goats and 130 sheep, originating from all 28 districts of Bulgaria. CCHF virus-specific IgG antibodies were observed in 165 cattle (19.6, 95% CI 17.0-22.4%), in 20 goats (22.7, 95% CI 15.2-32.5%) and in 10 sheep (7.7, 95% CI 4.2-13.6%). The highest seroprevalence was detected in the district of Kardzhali, South Bulgaria (86.7, 95% CI 73.8-93.7%), a well-known focus of CCHF in humans. The other two districts with human CCHF cases, Blagoevgrad and Burgas, located in southwest and southeast Bulgaria, showed overall seroprevalences in livestock of 41.9% (95% CI 28.4-56.7%) and 31.3% (95% CI 22.2-42.1%), respectively. Seroprevalences in districts with no history of human CCHF cases varied between 55% (95% CI 39.8-69.3%) and 22.5% (95% CI 12.3-37.5%). These results suggest frequent CCHF virus infections even in regions without human CCHF cases and an enhanced risk of infection for humans in close contact with the infected livestock.

Country-wide seroprevalence studies on Crimean-Congo hemorrhagic fever and hantavirus infections in general population of Bulgaria.

Crimean-Congo hemorrhagic fever (CCHF) and hantavirus infections are the two viral hemorrhagic fevers spread in Europe. To test actual circulation of CCHF virus (CCHFV) and hantaviruses in Bulgaria, we conducted country-wide seroepidemiological studies. Serum samples were collected prospectively from 1500 residents of all 28 districts in Bulgaria. CCHFV seroprevalence of 3.7% was revealed. Anamnesis for tick bites, contact with livestock, age over 40 years and residency in Haskovo district were found as risk factors. The highest CCHFV seroprevalence was observed in the known endemic districts in southeastern Bulgaria: Haskovo (28%) and Yambol (12%). Reactive samples were found in residents of 20 of the 28 districts in Bulgaria. In comparison with the previous studies, the data presented indicate that CCHFV increased substantially its circulation in the endemic regions and was introduced in many new areas. Hantavirus seroprevalence was based on results of the immunoblot and estimated as 3.1%. Surprisingly, contrary to all available data, Puumala virus seroprevalence rate was 2.3% versus 0.8% of Dobrava-Belgrade virus. Evidence for hantavirus IgG seropositivity was found in residents of 23 of the 28 districts in the country. The first hantavirus seroprevalence study in Bulgaria showed that Puumala virus is probably more wide-spread in the country than Dobrava-Belgrade virus.

A nationwide seroprevalence screening for West Nile virus and Tick-borne encephalitis virus in the population of Bulgaria.

To assess local circulation and risk for human infections with West Nile virus (WNV) and Tick-borne encephalitis virus (TBEV) in Bulgaria, a nationwide seroprevalence study was conducted. In total, 1451 residents of all 28 districts in Bulgaria were tested for WNV-specific and TBEV-specific IgG antibodies. The survey found overall seroprevalence of 1.5% and 0.6%, respectively. The highest WNV seroprevalence was found in Sofia Province and districts near the river Danube. TBEV circulation was detected among residents of six districts. The results showed that the two virus infections seem to be more wide-spread in the country as has been described.

Biosafety standards for working with Crimean-Congo hemorrhagic fever virus.

In countries from which Crimean-Congo haemorrhagic fever (CCHF) is absent, the causative virus, CCHF virus (CCHFV), is classified as a hazard group 4 agent and handled in containment level (CL)-4. In contrast, most endemic countries out of necessity have had to perform diagnostic tests under biosafety level (BSL)-2 or -3 conditions. In particular, Turkey and several of the Balkan countries have safely processed more than 100 000 samples over many years in BSL-2 laboratories. It is therefore advocated that biosafety requirements for CCHF diagnostic procedures should be revised, to allow the tests required to be performed under enhanced BSL-2 conditions with appropriate biosafety laboratory equipment and personal protective equipment used according to standardized protocols in the countries affected. Downgrading of CCHFV research work from CL-4, BSL-4 to CL-3, BSL-3 should also be considered.

Molecular epidemiology of Crimean-Congo hemorrhagic fever in Bulgaria--An update.

Crimean-Congo hemorrhagic fever (CCHF) is endemic in Bulgaria. During 2013-2014, 11 confirmed CCHF cases have been reported in the country (seven in 2013 and four in 2014). The present study provides the CCHF molecular epidemiology in Bulgaria based on all currently available S, M, and L RNA segment nucleotide sequences spanning the years 1978-2014. A relatively low genetic difference (0-6%, the maximum seen in the M RNA segment) was seen among the CCHFV sequences suggesting that a slow evolving CCHFV strain belonging to "Europe 1" clade is present in Bulgaria. Although the virus emerged in new foci during the recent years, it is more active in the established endemic foci which seem to offer the most suitable ecosystem and environment. Understanding the CCHF epidemiology and virus evolution is the basis for public health programs and vaccine design.

Clinical Evaluation of Rapid Diagnostic Test Kit Using the Polysaccharide as a Genus-Specific Diagnostic Antigen for Leptospirosis in Korea, Bulgaria, and Argentina.

Leptospirosis, a zoonotic disease that is caused by many serovars which are more than 200 in the world, is an emerging worldwide disease. Accurate and rapid diagnostic tests for leptospirosis are a critical step to diagnose the disease. There are some commercial kits available for diagnosis of leptospirosis, but the obscurity of a species- or genus-specific antigen of pathogenic Leptospira interrogans causes the reduced sensitivity and specificity. In this study, the polysaccharide derived from lipopolysaccharide (LPS) of nonpathogenic Leptospira biflexa serovar patoc was prepared, and the antigenicity was confirmed by immunoblot and enzyme linked immunosorbent assay (ELISA). The performance of the rapid diagnostic test (RDT) kit using the polysaccharide as a diagnostic antigen was evaluated in Korea, Bulgaria and Argentina. The sensitivity was 93.9%, 100%, and 81.0% and the specificity was 97.9%, 100%, and 95.4% in Korea (which is a rare region occurring with 2 serovars mostly), Bulgaria (epidemic region with 3 serovars chiefly) and Argentina (endemic region with 19 serovars mainly) respectively. These results indicate that this RDT is applicable for global diagnosis of leptospirosis. This rapid and effective diagnosis will be helpful for diagnosis and manage of leptospirosis to use and the polysaccharide of Leptospira may be called as genus specific antigen for diagnosis.

Crimean-Congo hemorrhagic fever: CXCL10 correlates with the viral load.

Crimean-Congo hemorrhagic fever (CCHF) is a human disease with high fatality rate. Although its pathogenesis is not elucidated yet, it is considered that cytokines play a significant role in the progression and outcome of the disease. Serum CXCL10 levels were estimated in 35 patients with acute CCHF and were correlated with the viral load, and various demographic and clinical parameters. The mean CXCL10 concentration in the patients' group was higher compared to the respective value in the control group (4421.74 pg/ml vs. 28.47 pg/ml, P < 0.05). A strong positive correlation between CXCL10 and viral load was seen (rs = 0.57, P < 0.001), while the outcome of the disease was related with the viral load (rs = 0.47, P = 0.004) and the presence of hemorrhagic manifestations (P < 0.001). The study provides an insight into the strong correlation between CXCL10 and viral load in acute CCHF cases suggesting that it plays an important role in CCHF pathogenesis.

Detection of Dobrava hantavirus RNA in Apodemus mice in Bulgaria.

Several Hantaviruses cause hemorrhagic fever with renal syndrome (HFRS) in Europe: Dobrava-Belgrade virus (DOBV), Puumala, Saaremaa, Sochi, and Seoul virus. Although HFRS is endemic in Bulgaria, genome sequences of hantaviruses have never been detected in wild rodents. To identify rodent reservoirs, a total of 691 rodents from three endemic regions were trapped in 2011-2012 and screened by TaqMan RT-PCR for detection of hantaviral genomic RNA. Partial small (S) and/or large (L)-segment sequences were recovered from six Apodemus mice: five of the species A. flavicollis and one A. agrarius. Phylogenetic analysis revealed that all recovered sequences belonged to DOBV. On the phylogenetic trees, the novel Bulgarian hantavirus sequences clustered together with sequences of established previously DOBV variants recovered from Bulgarian HFRS patients and also with variants found in wild rodents trapped in Slovenia, Greece, and Slovakia. One of the novel Bulgarian DOBV S-sequences from A. agrarius was related closely to DOBV sequences recovered from A. flavicollis, suggesting a spillover of DOBV from its natural host to A. agrarius mice. The results of this study confirmed the circulation of DOBV in wild rodents in Bulgaria. The complexity of the epidemiological situation in the Balkans requires further studies of hantaviruses in rodent hosts and human HFRS cases.

Seroprevalance of Crimean-Congo haemorrhagic fever in Bulgarian livestock.

Crimean-Congo haemorrhagic fever (CCHF) is a tick-borne zoonotic disease. Over the past decade, CCHF cases in humans have emerged in Turkey and reemerged in the Balkan countries, Ukraine and Tajikistan. Occupational contact with infected livestock has been recognized as a common cause of the disease. A cross-sectional seroprevalence study in livestock was conducted in farming communities of an endemic area in Bulgaria, southeastern Europe. Overall, 72% of the tested animals were positive for IgG antibodies to CCHF virus. By the time the animals were one-year old almost 50% had serologic evidence of CCHF infection, and by two years already 80% of them had been infected. The data obtained in this study reflect current situation of CCHF virus infection among livestock in Bulgaria. The results showed active CCHF virus circulation that poses risk for humans to be infected during contacts with animals and requires public health awareness.

Resurgence of influenza with increased genetic diversity of circulating viruses during the 2022-2023 season.

Introduction. After two seasons of absence and low circulation, influenza activity increased significantly in the winter of 2022-2023. This study aims to characterize virological and epidemiological aspects of influenza infection in Bulgaria during the 2022-2023 season and perform a phylogenetic/molecular analysis of the hemagglutinin (HA) and neuraminidase (NA) sequences of representative influenza strains.Hypothesis/Gap Statement. Influenza A and B viruses generate new genetic groups/clades each season, replacing previously circulating variants. This results in increased antigenic distances from current vaccine strains. Strengthening existing influenza surveillance is essential to meet the challenges posed by the co-circulation of influenza and SARS-CoV-2.Methodology. We tested 2713 clinical samples from patients with acute respiratory illnesses using a multiplex real-time RT-PCR kit (FluSC2) to detect influenza A/B and Severe acute respiratory syndrome coronavirus-2(SARS-CoV-2) simultaneously. Representative Bulgarian influenza strains were sequenced at the WHO Collaborating Centres in London, UK, and Atlanta, USA.Results. Influenza virus was detected in 694 (25.6 %) patients. Of these, 364 (52.4 %), 213 (30.7 %) and 117 (16.9 %) were positive for influenza A(H1N1)pdm09, A(H3N2) and B/Victoria lineage virus, respectively. HA genes of the 47 influenza A(H1N1)pdm09 viruses fell into clades 5a.2. and 5a.2a.1 within the 6B.5A.1A.5a.2 group. Twenty-seven A(H3N2) viruses belonging to subclades 2b, 2a.1, 2a.1b and 2a.3a.1 within the 3C.2a1b.2a.2 group were analysed. All 23 sequenced B/Victoria lineage viruses were classified into the V1A.3a.2 group. We identified amino acid substitutions in HA and NA compared with the vaccine strains, including several substitutions in the HA antigenic sites.Conclusion. The study's findings showed genetic diversity among the influenza A viruses and, to a lesser extent, among B viruses, circulating in the first season after the lifting of anti-COVID-19 measures.

Monkeypox in Bulgaria: Significance of Various Clinical Samples, Clinical Manifestation, and Molecular Detection.

Background/Objectives: Monkeypox (mpox) is currently the most common orthopoxvirus (OPXV) zoonotic disease, and, since 2022, there has been atypical person-to-person transmission observed in non-endemic countries. The present study aimed to investigate the frequency of monkeypox virus (MPXV) and OPXV DNA detection in recommended and alternative clinical materials taken during the acute and convalescent phases of infection in Bulgarian patients. Methods: The study included laboratory investigation by real time PCR of 181 clinical samples from 42 Bulgarian patients with possible mpox infections. Results: MPXV DNA was detected in 23/181 (12.71%), and OPXV DNA in 20/181 (11.05%) clinical samples. There were six mpox-confirmed patients aged 23 to 44. At the highest frequency, MPXV and OPXV DNA were detected in samples of vesicular contents (6/6) and nasal/oropharyngeal secretions (5/6 and 4/6) during the first three days from the appearance of clinical symptoms. We demonstrated MPXV and OPXV DNA in alternative samples (urine, feces, ejaculate, and saliva), and in follow-up patient samples, taken two weeks after mpox confirmation in the convalescent phase (vesicular contentsand urine). Conclusions: Our findings suggested that MPXV may be detected in a larger set of clinical materials, including alternatives, where the virus can persist for more than two weeks.

Epidemiological and Genetic Characteristics of Respiratory Viral Coinfections with Different Variants of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2).

This study aimed to determine the incidence and etiological, seasonal, and genetic characteristics of respiratory viral coinfections involving severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Between October 2020 and January 2024, nasopharyngeal samples were collected from 2277 SARS-CoV-2-positive patients. Two multiplex approaches were used to detect and sequence SARS-CoV-2, influenza A/B viruses, and other seasonal respiratory viruses: multiplex real-time polymerase chain reaction (PCR) and multiplex next-generation sequencing. Coinfections of SARS-CoV-2 with other respiratory viruses were detected in 164 (7.2%) patients. The most common co-infecting virus was respiratory syncytial virus (RSV) (38 cases, 1.7%), followed by bocavirus (BoV) (1.2%) and rhinovirus (RV) (1.1%). Patients ≤ 16 years of age had the highest rate (15%) of mixed infections. Whole-genome sequencing produced 19 complete genomes of seasonal respiratory viral co-pathogens, which were subjected to phylogenetic and amino acid analyses. The detected influenza viruses were classified into the genetic groups 6B.1A.5a.2a and 6B.1A.5a.2a.1 for A(H1N1)pdm09, 3C.2a1b.2a.2a.1 and 3C.2a.2b for A(H3N2), and V1A.3a.2 for the B/Victoria lineage. The RSV-B sequences belonged to the genetic group GB5.0.5a, with HAdV-C belonging to type 1, BoV to genotype VP1, and PIV3 to lineage 1a(i). Multiple amino acid substitutions were identified, including at the antibody-binding sites. This study provides insights into respiratory viral coinfections involving SARS-CoV-2 and reinforces the importance of genetic characterization of co-pathogens in the development of therapeutic and preventive strategies.