key: cord-0004309-51tkew79
authors: Gizachew, Mucheye; Tiruneh, Moges; Moges, Feleke; Adefris, Mulat; Tigabu, Zemene; Tessema, Belay
title: Molecular characterization of Streptococcus agalactiae isolated from pregnant women and newborns at the University of Gondar Comprehensive Specialized Hospital, Northwest Ethiopia
date: 2020-01-13
journal: BMC Infect Dis
DOI: 10.1186/s12879-020-4776-7
sha: 24b1cd9fce054cf3c1aaa7351d7981d9c0cb83d5
doc_id: 4309
cord_uid: 51tkew79

BACKGROUND: Streptococcus agalctiae (Group B Streptococcus, GBS) is a perinatal pathogen and a leading cause of neonatal infections worldwide. Serotype, sequence type, clonality, antibiotic resistance genes and surface protein profiles of GBS are scarce in Ethiopia, a reason that this study was planned to investigate. . METHODS: Sixteen colonizing GBS isolates obtained from recto-vaginal swabs of pregnant women and body surfaces of newborns were further analyzed. Minimum inhibitory concentration (MIC) test, and whole genome sequence (WGS) methods were done for antibiotic susceptibility test, and molecular characterization of the isolates. RESULTS: All the GBS isolates analyzed were belonged to four capsular serotypes: II, 11/16(68.8%), V, 3/16(18.8%), Ia and VI each with 1/16(6.3%) and five sequence type (ST-2, ST-10, ST-14, ST-569 and ST-933). Sequence type-10 was the most predominant ST followed by ST-569. The five STs were grouped into the four clonal complexes (CC - 1, CC-10, CC-19, and CC-23). Different surface proteins and pili families such as ALP1, ALPHA, ALP23, PI-1 / PI-2A1, PI-1 / PI-2B, and Srr1 were detected from WGS data. All isolates were found to be susceptible to the tested antibiotics except for tetracycline in MIC and WGS test methods used. Tetracycline resistant determinant genes such as TETM and TETL / TETM combination were identified. CONCLUSION: Further studies on serotype and molecular epidemiology will provide a comprehensive data of the GBS capsular serotype and clones available in Ethiopia.

Group B Streptococcus (GBS) or Streptococcus agalactiae remains a leading cause of neonatal sepsis, pneumonia and meningitis, often associated with high morbidity and mortality in Europe, USA, and Australia even though there was a wide use of intrapartum antibiotic prophylaxis (IAP) since the mid 1990s. However, molecular profiles of GBS in many countries outside of these regions was not well documented [1, 2] . GBS causes early-onset neonatal disease (EOD) which is characterized by onset of clinical symptoms during the first week of life (0-6 days) following birth; it usually results in vertical transmission from colonized pregnant women to newborns during or just before delivery [3] . The other form of GBS disease appears in infants is the late-onset disease (LOD) -characterized by the onset of symptoms within 7 days to 89 days of age. It can be acquired from the mother or environmental sources [4, 5] . It is also an important cause of preterm delivery, stillbirth, and puerperal sepsis [6, 7] .

To colonize and cause diseases, GBS use a collection of virulence factors; one of the most important being a capsular polysaccharide (CPS) [8] . Human GBS has currently 10 serotypes, Ia, Ib, II to IX, based on serotype specific antigens [9] and the current candidate capsular polysaccharide conjugate vaccines target only a subset of these [10] . A study described that the five most common serotypes, Ia, Ib, II, III, and V accounted for more than 85% of serotypes in global regions that have reported serotype data, including the Americas (96%), Europe (93%), and the Western Pacific (89%) [11] , and Africa (91.8%) [12] . A review demonstrated that among the ten serotypes, Ia, Ib, II, III, and V accounts for 98% of the identified colonizing isolates worldwide. Serotype III, associated with invasive disease, accounts for 25%, but is less frequent in some South American and Asian countries, but serotypes VI-IX are more common in Asia [13] . Another report from Iran showed that serotypes V (19.6%), II (12.5%) and IV (12.5%) were the most frequent followed by serotypes III (10.7%) and VI (10.7%), Ib (8.9%), Ia (7/ 1%), VII (5/3%) and VIII (5/3%); and 7.1% of strains were non-type-able [14] . A study in Ghana showed as serotypes VII (38.5 to 42.9%) and IX(26.9 to 32.1%) were the most common serotypes [15] . Serotyping of GBS is useful to understand the local epidemiology, for monitoring of serotype replacement or capsular switching, and for contributing existing serotype profiles in the area for rational, effective and broad serotype coverage GBS vaccine development [16, 17] . But, there is scarcity of serotype profile data of GBS colonizing pregnant women and newborns in Northwest Ethiopia for the last three decades.

GBS remained susceptible to the beta-lactam antibiotics, but, resistance to macrolides, lincosamides, fluoroquinolones and other antibiotics used as alternative therapy has been reported [18] . Macrolide resistance in GBS is represented by two mechanisms: Target Site Modification by Erythromycin ribosomal methylase, mediated by ermB, ermA, ermTR, or ermC genes which confers cross resistance to macrolides, lincosamides, and streptogramin B (MLS B phenotype) [19] . Erm genes encode methylase 23S rRNA, which is responsible for methylation of erythromycin and clindamycin receptor sites in ribosomes [20] . This resistance can be constitutive macrolide lincosamidestreptogramin B (cMLS B ) resistance, as well as inducible macrolide lincosamidestreptogramin B (iMLS B ) resistance [21] . The second resistance mechanism that found in GBS is Macrolide efflux pump, mediated by mef genes, which confers resistance to 14 and 15 member macrolides only. In addition, a novel efflux system distinct from the Mef pump and encoded by mreA (for macrolide resistance efflux) was reported in a unique strain of S. agalactiae COH31 γ/δ [22] . Susceptible GBS isolates also showed to possess the mreA gene, and it might function as a housekeeping gene [23] .

Regarding the study of antimicrobial resistance markers (genes), a study revealed that tetM was the most frequently (84%) identified in all groups. Macrolide resistance genes were found in a small proportion of the isolates (8%) and each of the three relevant genes was represented as 3.3% ermB, 2.5% ermA, and 2.3% mefA/ E. In addition, one isolate contained both ermA and ermB, and one contained both ermA and mefA/E [24] . In another study, macrolide resistance genes were screened and the resistant rate in total isolates was reported as 69.5%. Ribosome methylation genes (ermB, ermA/TR, ermC) were screened in all examined isolates and the most prevalence gene (63.04%) was ermB,and ermA\TR gene accounted 23.9%, but ermC gene was not detected. Findings of the study displayed that the prevalence rate of efflux pumps encoding genes (mefA and mreA) were 8.69 and 69.5% respectively. The study also showed that 14 (30.4%) of the isolates were not recognized as a carrier of any erythromycin resistance genes. One isolate among the screened isolates was harbored four erythromycin resistance genes except erm C gene [25] . Muller et al. reported that resistance to erythromycin and clindamycin was found in 8 and 7%, respectively. Macrolide resistance genes mef(A), erm (TR) and erm(B) were found in one, two, and five isolates, respectively; only 5/8 (62.5%) of the isolates exhibited both genotypic as well as phenotypic resistance. One genotype occurred in 36% of the subset [26] .

A study detailed that GBS has different surface proteins such as alpha, Alp1, Alp2, Alp3, Alp4, and Rib. These proteins are collectively named as alpha-like protein (Alp) family [21, 27] . The biological role of Alp proteins is not well documented though they are targets for protective antibodies demonstrated in the animal model. GBS uses alpha protein to invade the epithelial cells by interacting with host cell glycosaminoglycans [21, 28, 29] . Genetic profiles of GBS strains obtained from around the globe is investigated by different molecular typing methods including multi-locus sequence type (MLST) [21, 27, 30] . However, there is scarcity of data regarding to the antimicrobial resistance gene markers plus sequence types and surface protein profiles of GBS in Northwest Ethiopia. Therefore, this study was aimed to determine the serotypes, antibiotic resistance genes, sequence types and surface protein profiles of GBS recovered from pregnant women and newborns in Northwest Ethiopia.

It was conducted at the University of Gondar Comprehensive Specialized Hospital. Northwest Ethiopia. The hospital serves more than 5, 000,000 people and it has about 450 to 600 pregnant women admission services a month. So far, there is no routine screening of pregnant women and provision of IAP for GBS has been established in the hospital [31] .

Colonizing GBS isolates were recovered from recto-vaginal area of pregnant women at time of delivery and body surfaces of newborns immediately following delivery and stored in deep freezer (− 80°C). The GBS isolates were randomly picked from a stored collection in a deep freezer. Of these, ten were from recto-vaginal swabs of pregnant women, and six were from the newborns (three nasal, two eye and one umbilical swabs). The isolates were further processed by using streptococci Lancefield grouping GBS latex slide agglutination kits and short-read whole genome sequencing (WGS) as previously described [32] . This study protocol was in accordance with the ethical standards of the responsible regional committee on human experimentation and the Declaration of Helsinki of 1975 (revised in 1983). It was approved by the ethics committees at the University of Gondar. Informed written consent was obtained from the pregnant women and care providers of newborns enrolled in this study.

Sixteen isolates were subjected to broth micro-dilution testing (BDT) for determination of Minimum Inhibitory Concentration (MIC) and compared it with WGS-based characterization of antibiotic susceptibility patterns including resistance genes in the GBS isolate. Antibiotics included six β-lactams (ceftizoxime, cefoxitin, cefotaxime, cefazolin, ampicillin and penicillin), erythromycin, clindamycin, tetracycline, levofloxacin, ciprofloxacin, daptomycin, vancomycin and linezolid) were tested. A well containing both erythromycin and clindamycin detected inducible clindamycin-resistance, as previously described [33] . In addition, these isolates were test for their antimicrobial susceptibility test by using the disc diffusion method as per the CLSI, 2014 guideline [34] .

Whole genome sequencing 'The WGS of GBS was performed at the CDC streptococcus laboratory in the Atlanta, Georgia (USA), and drug resistant genes, serotypes, sequence types, clonal complexity and surface protein profiles were predicted from the WGS data. GBS colonies were cultured on Todd Hewitt Broth with 0.5% yeast and incubated overnight at 37°C in 5%CO 2 . Genomic DNA for WGS was extracted manually using a modified QIAamp DNA mini kit protocol (Qiagen, Inc., Valencia, CA, USA). Nucleic acid concentration was quantified by an Invitrogen Qubit 3.0 Flurometer assay (Thermo Fisher Scientific Inc.,Waltham, MA, USA) and samples were sheared using a CovarisM220 ultrasonicator (Covaris, Inc.,Woburn, MA, USA) programmed to generate 500-bp fragments. Libraries were constructed on theSciCloneG3 (PerkinElmer Inc., Waltham, MA, USA) using a Tru-SeqDNA PCR-Free HT library preparation kit with 96 dual indices (Illumina Inc., San Diego, CA, USA) and quantified by a KAPA qPCR library quantification method (Kapa Biosystems Inc., Wilmington, MA, USA). WGS was generated employing MiSeqinstrument and the MiSeq v2 500 cycle kit (Illumina Inc).WGS-based identification of serotypes employed the query DNA sequences described in the Additional file 1.

To identify the clonality, the sequence types (STs) of GBS isolates were investigated using the seven housekeeping genes: adhp, atr, glck, glna, phes, sdha and tkt by comparing with the standard references available at the MLST 1.8 database (https://cge.cbs.dtu.dk//services/ MLST/). To visualize the possible evolutionary relationships between isolates, STs of the study isolates and the globally reported strains were computed using PHYLO-ViZ software v2.0 based on Global optimal eBURST (goeBURST) algorithm [35] . Table 1 , the genomes size of the GBS ranged from 2,018,343 to 3,891,488 bp. The numbers of contig in the DNA sequences per genome ranged from 29 to 1797.Of the 16 isolates sequenced, four serotypes were identified as, serotype II, 11/16(68.8%), V, 3/ 16(18.8%), Ia and VI each with 1/16(6.3%). Capsular polysaccharide serotypes II and V were accounted for 14/16(87.5%). Capsular serotype II accounted for 5/ 6(83.3%), and 6/10(60.0%) among the neonatal and pregnant women isolates analyzed respectively. Antibiotic susceptibility testing was performed for all 16 GBS isolates by using MIC. Results have shown that, they were sensitive to all antibiotics tested except for tetracycline to which 93.8% of GBS isolates showed resistance Antibiotic susceptibility test results from sequenced data also showed that all isolates were susceptible except for tetracycline to which the 93.8% isolates showed resistance. Tetracycline resistance was predominantly due to tetM, which was detected alone in the 56.3% of the isolates, or in 37.5% of the isolates in association with tetL. However, there are no genes encoding drug resistance detected for other antibiotics tested. In addition, penicillin binding protein (PBP1A and PBP2X) for beta-lactam antibiotic was not detected in our isolates. Of the surface proteins screened from our isolates, the alpha-like protein (ALPH), serine-rich repeat-1 (SRR1),and Pili, were detected. Among the ALPH families, 43.8% of the isolates had ALPHA, followed by ALP1 (25.0%), and all isolates harbored at least one pilus islands (PI), of which, the PI1/PI2A1 combination was the most prevalent subfamily of Pili that was identified among the 14/16 (87.5%) isolates ( Table 1) .

The MLST data reveals that the study isolates showed five sequence types (STs), and four clonal complexes (CCs) which are of different phylogeny and is suggestive of different clones circulating in Northwest Ethiopia. Overall, the most common STs were ST-10 (37.5%) followed by ST-569 (31.3%). Sequence type (ST)-14 was represented by a single GBS isolate. The predominant CCs were CC-10 (68.8%) followed by CC-19 (18.8%). Clonal complex 1 and 23 each was represented by a single isolate. goeBURST reveals about 68 clonal complexes from the available datasets globally. The goeBURST diagram shows that the sequence types STs-10 and STs- 569, a double locus variant (DLV) of ST-10 ( Fig. 1b) , from the study isolates was belonging to clonal complex (CC) 10 with ST-10 as founder ST, whereas, ST-2, a subgroup founder, was belonging to the CC-19, ST-14 was grouped into the CC-1, and ST-933 was categorized into the CC-23 (Table 1 , and Fig. 1a ).

In addition, Fig. 2 showed that ST-569 and ST-10 were in the double locus variant (DLV), where as the rest of the STs were with the single locus variants (SLV), the STs those which were directly connected to the founder and they differ from the founder by only one single locus or allele. The size of the black dots indicate the number of the isolates included in the specified STs in which the larger the size, the higher the number of the isolates included. ST-14, and ST-933 were 1/16(6.25%) each, but STs-10 were 6/16 (37.5%).

As it is shown in Table 2 , the 16 GBS isolates were tested for about 15 different antibiotics by using either disc diffusion, minimum inhibitory concentration (MIC) or short-read whole genome sequencing methods. Certain isolates showed resistance to different antibiotics tested ranging between 1(6.25%) for penicillin to 14(87.5%) for tetracycline in the disc diffusion test method.

In this study, serotype II represented by ST-10 and ST-569, belongs to CC-10 was the most prevalent colonizer of newborns. It is in congruent with that of a study conducted in Taiwan where serotype II was the most colonizer,4/11 (33.3%) [36] . However, it is different from a study conducted before three decades in the same hospital where 60% were type Ib followed by type Ia (16%) [37] ; in Poland, serotype III (35%) was the most prevalent types followed by Ia (20%), V (17%), II (15%), Ib (8%) and IV (5%) [38] . Similarly, a study from Iran showed that serotype III (50%) was the most common serotype, followed by serotype II (25%), Ia (12%), V (11%), and Ib (2%) [39] .

Regarding the MLST analysis, we found that the most predominant ST was ST-10 which is quite different from [42] . The variation observed might be attributed to the geographical variations, the type of study participants involved in which certain studies isolated invasive GBS while in our case is the colonizing GBS isolated from recto-vaginal swabs of pregnant women and body surfaces of newborns. This variation also might be attributed to the limited number of GBS isolates sequenced in our study. The most predominant CCs we found was CC-10 followed by CC-19, which is slightly congruent with the reports from Poland where CC-23 (49.0%), CC-19 (17.0%), and CC-10 and CC-17 each with 10.0% [43] , and in Taiwan, CC- [44] . In our study, ALPH, SRR and Pilli were identified. Of the ALPH families, ALPHA was found to be the predominant surface protein followed by ALP1. All our isolates carried SRR1, a glycoprotein that may play an important role in crossing the blood brain-barrier (BBB) and for the development of streptococcal meningitis [45, 46] . Serine-rich repeat proteins also help to promote initial contact of GBS with host epithelium, cervical and vaginal attachment [47] which result in GBS cervicovaginal colonization and provide a survival benefit in the vaginal mucosa. In addition, three forms of pilus proteins (PI-1 / PI-2A1, PI-1 / PI-2B, and PI-1 / PI-2A3) were identified in our study. A study showed that pilus vaccine containing the PI-1, PI-2a, and PI-2b elicited good opsonophagocytic antibodies that confer protection in mice and this suggested that pilus components of GBS are highly immunogenic [48, 49] . It implies that increasing our knowledge about the profiles of the pilus island (PI) among GBS isolates obtained from different geographic locations is useful for future efforts that aimed at the development of maternal GBS vaccines from bacterial protein components (pilusbased), which were effective in animal models. This variable presence of specific PIs has considerable implications for the development of maternal GBS vaccines which target these bacterial proteins.

The drug susceptibility patterns of the 16 GBS isolates tested by BDI and WGS seem to resemble what is seen in other countries, such as studies conducted in Iran [40] and Poland [43] showed that, all GBS isolates were susceptible to penicillin, vancomycin, and linezolid, Improving the laboratory facilities for routine antibiotic susceptibility tests, mainly in developing countries, should be mandatory.

The prevalent tetracycline resistance was caused by the presence of the TETL.TETM genes. In agreement with our results, a study described that tetracycline resistance was predominantly due to tetM (83.7%) or along with tetL (12.2%) [50] . Another study conducted in Italy also revealed that the tetracycline resistance rate was 62/91 (68.1%), all resistant isolates harbored the tet(M) gene [21] . On the same way, resistant to tetracycline was observed in all 19 (100%) strains and 18/19 (98.0%) was correlated with presence of the tetM gene in a study carried out in Iran [40] . Tetracycline resistance is due to acquisition of tet determinant that encodes for antibiotic efflux or ribosomal protection in gram positive cocci [51, 52] .

The limitation of this study was that the numbers of the GBS isolates sequenced were small.

In this study, the MLST data reveals that the study isolates are of different phylogeny and is suggestive of different clones circulating in Ethiopia. WGS data showed that tetM and tetL / tetM were found as the tetracycline resistant determinant genes in 93.8% of the GBS isolates. All other antibiotics tested for by using MIC and WGS methods were susceptible. We found that penicillin is still the drug of choice for maternal GBS IAP. Vaccine formulation would include serotype II, V, Ia, and VI, besides serotype Ib, and III to give broader coverage of GBS prevention across different geographical locations including the study area particularly and in the globe generally since serotype VI was also often reported from the Asian studies. Further studies on serotype and molecular epidemiology will provide a comprehensive data of the GBS capsular serotypes, sequence types, clones, surface proteins and drug resistance determinant genes available in Ethiopia.

Supplementary information accompanies this paper at https://doi.org/10. 1186/s12879-020-4776-7.

Additional file 1: Whole genome sequencing (WGS) and Minimum inhibitory concentration (MIC) by using broth dilution method test results of colonizing GBS. 

Group B streptococcal disease in infants aged younger than 3 months: systematic review and meta-analysis

Early onset neonatal sepsis: the burden of group B streptococcal and E coli disease continues

An overview of the natural history of early onset group B streptococcal disease in the UK

Late and ultra late onset Streptococcus B meningitis: clinical and bacteriological data over 6 years in France

Evaluation of universal antenatal screening for group B streptococcus

An update on vaccination against group B streptococcus

Group B streptococcal disease: its importance in the developing world and prospect for prevention with vaccines

Prevention of C3 deposition by capsular polysaccharide is a virulence mechanism of type III group B streptococci

An overview of global GBS epidemiology

Real-time PCR-based serotyping of Streptococcus agalactiae

Group B streptococcal epidemiology and vaccine needs in developed countries

Streptococcus agalactiae maternal colonization, antibiotic resistance and serotype proiles in Africa: a meta-analysis

Maternal colonization with group B Streptococcus and serotype distribution worldwide: systematic review and meta-analyses

Capsular serotype and antibiotic resistance of group B streptococci isolated from pregnant women in Ardabil, Iran

Frimodt-Møller N. Carriage and serotype distribution of Streptococcus agalactiae in third trimester pregnancy in southern Ghana

Evidence for rare capsular switching in Streptococcus agalactiae

Capsular switching in group B Streptococcus CC17 hypervirulent clone: a future challenge for polysaccharide vaccine development

Antimicrobial resistance profiles and genetic characterisation of macrolide resistant isolates of Streptococcus agalactiae

Three new macrolide efflux (mef) gene variants in Streptococcus agalactiae

Genetic characterization and diversity of Streptococcus agalactiae isolates with macrolideresistance

Molecular epidemiology and distribution of serotypes, surface proteins, and antibiotic resistance among group B streptococci in Italy

Cloning and characterization of a novel macrolide efflux gene, mreA, from Streptococcus agalactiae

Phenotype and genotype of some clinical group B Streptococcus isolates resistant to erythromycin in Egypt

Distribution of genotypes and antibiotic resistance genes among invasive Streptococcus agalactiae (group B streptococcus) isolates from Australasian patients belonging to different age groups

Prevalence of macrolide resistance genes among group B streptococci in pregnant women

Low rate of carriage of macrolide-resistant group B streptococci in pregnant women in The Netherlands

Surface proteins of Streptococcus agalactiae and related proteins in other bacterial pathogens

Alpha C protein of group B Streptococcus binds host cell surface glycosaminoglycan and enters cells by an actin-dependent mechanism

Identification of a glycosaminoglycan-binding region of the alpha C protein that mediates entry of group B streptococci into host cells

Comparison of a 3-set genotyping system with multilocus sequence typing for Streptococcus agalactiae (group B streptococcus)

Newborn colonization and antibiotic susceptibility patterns of Streptococcus agalactiae at the University of Gondar Referral Hospital

Shortread whole genome sequencing for determination of antimicrobial resistance mechanisms and capsular serotypes of current invasive Streptococcus agalactiae recovered in the USA

Evaluation of disk approximation and single-well broth tests for detection of inducible clindamycin resistance in Streptococcus pneumoniae

Performance Standards for Antimicrobial Susceptibility Testing; Twenty-Fourth Informational Supplement. CLSI document M100-S24 Wayne

PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods

Characterization of group B Streptococcus colonization in full-term and lete-preterm neonates in Taiwan

Colonization of pregnant women and their newborn infants with group B streptococci in the Gondar College of Medical Sciences

Molecular characterization of capsular polysaccharides and surface protein genes in relation to genetic similarity of group B streptococci isolated from polish pregnant women

Molecular characterization of Streptococcus agalactiae isolates from pregnant and non-pregnant women at Yazd University hospital

Characterization of virulence factors, antimicrobial resistance pattern and clonal complexes of group B streptococci isolated from neonates

Associations between capsular serotype, multilocus sequence type, and macrolide resistance in Streptococcus agalactiae isolates from Japanese infants with invasive infections

Multilocus Sequence Typing of Swedish Invasive Group B Streptococcus Isolates Indicates a Neonatally Associated Genetic Lineage and Capsule Switching

Multilocus sequence types of invasive and colonizing neonatal group B streptococci in Poland

Multilocus sequence typing of invasive group B Streptococcus in central area of Taiwan

The alpha C protein mediates internalization of group B Streptococcus within human cervical epithelial cells

Protection against experimental infection with group B streptococcus by immunization with a bivalent protein vaccine

Serine-rich repeat proteins and pili promote Streptococcus agalactiae colonization of the vaginal tract

Preventing bacterial infections with pilus-based vaccines: the group B Streptococcus paradigm

Distribution of pilus islands in Streptococcus agalactiae that cause human infections: insights into evolution and implication for vaccine development

Macrolide-and tetracycline-resistance determinants of colonizing group B streptococcus in women in Egypt

Sixmonth multicenter study on invasive infections due to group B streptococci in Argentina

Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria That Grow Aerobically

Approved Standard, 9th. CLSI document M07-A9

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The authors would like to acknowledge the financial and material support from the College of Medicine and Health Sciences, University of Gondar, Northwest Ethiopia, the German Academic Exchange Service (DAAD), Germany, and the Institute of Virology, Leipzig University, Germany. We extend our appreciation to CDC Streptococcus laboratory team including Lesley McGee for their valuable contribution in sequencing the isolates. The authors are also grateful to the study participants' data collectors.

Authors' contributions MG participated in conception and design of the study, acquisition of data, analysis and Interpretation of findings. BT participated in conception and design of the study, and interpretation of data; he also participated in revision of the manuscript. MT participated in design of the study, acquisition and interpretation data; in revision of the manuscript. FM participated in conception of the study; and revision of the manuscript. MA participated in design of the study, and revision of the manuscript. ZT participated in interpretation of data, drafting and revision of the manuscript. All authors read and approved the final manuscript.

Financial and material support were from the University of Gondar, Ethiopia, the German Academic Exchange Service (DAAD) In-Country scholarship and the Institute of Virology, Leipzig University, Germany for salary, material procurement, data collection, analysis, and interpretation and in writing the manuscript, printing and copying services, and CDC for whole genome sequencing.