Genomic insights into the expansion of carbapenem-resistant Klebsiella pneumoniae within Portuguese hospitals.

Carbapenem-resistant Klebsiella pneumoniae (CR-KP) are a public health concern, causing infections with a high mortality rate, limited therapeutic options and challenging infection control strategies. In Portugal, the CR-KP rate has increased sharply, but the factors associated with this increase are poorly explored. In order to address this question, phylogenetic and resistome analysis were used to compare the draft genomes of 200 CR-KP isolates collected in 2017-2019 from five hospitals in the Lisbon region, Portugal. Most CR-KP belonged to sequence type (ST) 13 (29%), ST17 (15%), ST348 (13%), ST231 (12%) and ST147 (7%). Carbapenem resistance was conferred mostly by the presence of KPC-3 (74%) or OXA-181 (18%), which were associated with IncF/IncN and IncX plasmids, respectively. Almost all isolates were multi-drug resistant, harbouring resistance determinants to aminoglycosides, beta-lactams, trimethoprim, fosfomycin, quinolones and sulphonamides. In addition, 11% of isolates were resistant to colistin. Colonizing and infecting isolates were highly related, and most colonized patients (89%) reported a previous hospitalization. Moreover, among the 171 events of cross-dissemination identified by core genome multi-locus sequence typing data analysis (fewer than five allelic differences), 41 occurred between different hospitals and 130 occurred within the same hospital. The results suggest that CR-KP dissemination in the Lisbon region results from acquisition of carbapenemases in mobile genetic elements, influx of CR-KP into the hospitals by colonized ambulatory patients, and transmission of CR-KP within and between hospitals. Prudent use of carbapenems, patient screening at hospital entry, and improvement of infection control are needed to decrease the burden of CR-KP infection in Portugal.

Insights into the population structure and pan-genome of Haemophilus influenzae.

The human-restricted bacterium Haemophilus influenzae is responsible for respiratory infections in both children and adults. While colonization begins in the upper airways, it can spread throughout the respiratory tract potentially leading to invasive infections. Although the spread of H. influenzae serotype b (Hib) has been prevented by vaccination, the emergence of infections by other serotypes as well as by non-typeable isolates (NTHi) have been observed, prompting the need for novel prevention strategies. Here, we aimed to study the population structure of H. influenzae and to get some insights into its pan-genome. We studied 305H. influenzae strains, enrolling 217 publicly available genomes, as well as 88 newly sequenced H. influenzae invasive strains isolated in Portugal, spanning a 24-year period. NTHi isolates presented a core-SNP-based genetic diversity about 10-fold higher than the one observed for Hib. The analysis of key factors involved in pathogenesis, such as lipooligosaccharides, hemagglutinating pili and High Molecular Weight-adhesins, suggests that NTHi shape its virulence repertoire, either by acquisition and loss of genes or by SNP-based diversification, likely towards host immune evasion and persistence. Discreet NTHi subpopulations structures are proposed based on core-genome supported with 17 candidate genetic markers identified in the accessory genome. Additionally, this study provides two bioinformatics tools for in silico rapid identification of H. influenzae serotypes and NTHi clades previously proposed, obviating laboratory-based demanding procedures. The present study constitutes an important genomic framework that could lay way for future studies on the genetic determinants underlying invasiveness and disease and population structure of H. influenzae.

Cost-effectiveness of an Adjuvanted Recombinant Zoster Vaccine in older adults in the United States.

In the United States, herpes zoster (HZ) and related complications are estimated to result in approximately $1.3 billion in medical care costs and $1.7 billion in indirect costs annually. In this study, we compared the cost-effectiveness of a new Adjuvanted Recombinant Zoster Vaccine (RZV), containing recombinant varicella-zoster virus glycoprotein E and the AS01B Adjuvant System, versus No Vaccine, as well as versus the live attenuated HZ vaccine (Zoster Vaccine Live (ZVL)) in subjects aged 60+ years of age (YOA) and other age cohorts aged 50+ YOA. A multi-cohort Markov model was developed which follows 1 million individuals over their remaining lifetimes from the year of vaccination with annual cycle lengths. Second dose compliance for RZV was assumed to be 69%. Efficacy and waning parameters were derived from clinical trials for both vaccines. Epidemiological parameters, costs and utility model inputs were derived from US-specific population-based data. Costs and outcomes were discounted at 3% per year. Deterministic and probabilistic sensitivity analysis, along with scenario and threshold analysis were carried out to explore the overall uncertainty in the model. The model estimated that, compared to No Vaccine against HZ, RZV would prevent 103,603 HZ cases, 11,197 postherpetic neuralgia (PHN) cases, and 14,455 other complications, at an incremental cost of $11,863 per quality-adjusted life-year saved from a societal perspective. Compared to ZVL, the model estimated that, RZV would prevent 71,638 additional HZ cases, 6403 PHN cases, and over 10,582 other complications, resulting in net total societal cost savings of over $96 million. The results were robust to a wide range of sensitivity analyses. Vaccination against HZ with RZV is cost-effective compared to No Vaccine and cost-saving compared to ZVL, in the US population aged 60+ YOA. Clinicaltrial.gov. registered#: NA.

A primer on microbial bioinformatics for nonbioinformaticians.

BACKGROUND: Presently, the bottleneck in the deployment of high-throughput sequencing technology is the ability to analyse the increasing amount of data produced in a fit-for-purpose manner. The field of microbial bioinformatics is thriving and quickly adapting to technological changes, which creates difficulties for nonbioinformaticians in following the complexity and increasingly obscure jargon of this field. AIMS: This review is directed towards nonbioinformaticians who wish to gain understanding of the overall microbial bioinformatic processes, from raw data obtained from sequencers to final outputs. SOURCES: The software and analytical strategies reviewed are based on the personal experience of the authors. CONTENT: The bioinformatic processes of transforming raw reads to actionable information in a clinical and epidemiologic context is explained. We review the advantages and limitations of two major strategies currently applied: read mapping, which is the comparison with a predefined reference genome, and de novo assembly, which is the unguided assembly of the raw data. Finally, we discuss the main analytical methodologies and the most frequently used freely available software and its application in the context of bacterial infectious disease management. IMPLICATIONS: High-throughput sequencing technologies are overhauling outbreak investigation and epidemiologic surveillance while creating new challenges due to the amount and complexity of data generated. The continuously evolving field of microbial bioinformatics is required for stakeholders to fully harness the power of these new technologies.

Bioinformatics in bacterial molecular epidemiology and public health: databases, tools and the next-generation sequencing revolution.

Advances in typing methodologies have been the driving force in the field of molecular epidemiology of pathogens. The development of molecular methodologies, and more recently of DNA sequencing methods to complement and improve phenotypic identification methods, was accompanied by the generation of large amounts of data and the need to develop ways of storing and analysing them. Simultaneously, advances in computing allowed the development of specialised algorithms for image analysis, data sharing and integration, and for mining the ever larger amounts of accumulated data. In this review, we will discuss how bioinformatics accompanied the changes in bacterial molecular epidemiology. We will discuss the benefits for public health of specialised online typing databases and algorithms allowing for real-time data analysis and visualisation. The impact of the new and disruptive next-generation sequencing methodologies will be evaluated, and we will look ahead into these novel challenges.

Comparison of molecular typing methods for characterization of Staphylococcus epidermidis: proposal for clone definition.

In the present study we give some direction on the selection of the most appropriate typing method(s) to be used for the characterization of Staphylococcus epidermidis, in view of the most recent findings on the evolution, population structure, and epidemiology of this species. In order to achieve this aim, quantitative assessment of the correlation of the results of three typing methods--pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and staphylococcal chromosomal cassette mec (SCCmec) typing, which target different regions of the chromosome that evolve at different rates--was performed. In order to evaluate the discriminatory ability and the strength and direction of the correlation of the different typing methods, Simpson's index of diversity (SID), the adjusted Rand coefficient (AR), and the Wallace coefficient (W) were calculated. PFGE was the most discriminatory method (SID = 99%), followed by MLST (SID = 90%) and SCCmec typing (SID = 75%). The values of AR and W (0.10 < AR < 0.30; 0.50 < W < 0.75) indicated that the partition of the same isolate collection by PFGE, MLST, and SCCmec typing provided results that had only a poor correlation with each other. However, the information provided by the combination of PFGE and SCCmec enabled the prediction of the results obtained by MLST at the level of the clonal complex with a high degree of precision (W > 0.90). We propose that clones of S. epidermidis be defined by the combination of the PFGE type followed by the SCCmec type, which provides reliable information on the short-term epidemiology and the ability to predict with consistency long-term clonal evolution.

Intracranial vascular anomalies in patients with periorbital lymphatic and lymphaticovenous malformations.

BACKGROUND AND PURPOSE: The purpose of this study was to determine the nature, incidence, and radiologic appearance of intracranial vascular anomalies that occur in association with periorbital lymphatic malformation (LM) and lymphaticovenous malformation (LVM). MATERIALS AND METHODS: We retrospectively reviewed clinical records and imaging studies of 33 patients ranging in age from the neonatal period to 39 years (mean age, 5.1 years; median age, 1.0 year) who were evaluated for orbital LM or LVM at our institution between 1953 and 2002. Imaging studies, including CT, MR imaging, and cerebral angiograms, were evaluated by 2 radiologists to determine morphologic features of orbital LM and to identify associated noncontiguous intracranial vascular and parenchymal anomalies, including arteriovenous malformations (AVM), cerebral cavernous malformations (CCM), developmental venous anomalies (DVA), dural arteriovenous malformations (DAVM), and sinus pericranii (SP). RESULTS: The malformation was left-sided in 70% of patients. Twenty-two patients (70%) had intracranial vascular anomalies: DVA (n = 20; 61%), CCM (n = 2; 6%), DAVM (n = 4; 12%), pial AVM (n = 1; 3%), and SP (n = 1; 3%). Arterial shunts were present in the soft tissues in 2 patients (6%). Three patients had jugular venous anomalies. Three patients (9%) had cerebral hemiatrophy, 2 (6%) had focal cerebral atrophy, and 2 had Chiari I malformation. CONCLUSIONS: Intracranial vascular anomalies, some of which are potentially symptomatic and require treatment, are present in more than two thirds of patients with periorbital LM. Initial imaging of patients with orbital LM should include the brain as well as the orbit.

Illustration of a common framework for relating multiple typing methods by application to macrolide-resistant Streptococcus pyogenes.

The studies that correlate the results obtained by different typing methodologies rely solely on qualitative comparisons of the groups defined by each methodology. We propose a framework of measures for the quantitative assessment of correspondences between different typing methods as a first step to the global mapping of type equivalences. A collection of 325 macrolide-resistant Streptococcus pyogenes isolates associated with pharyngitis cases in Portugal was used to benchmark the proposed measures. All isolates were characterized by macrolide resistance phenotyping, T serotyping, emm sequence typing, and pulsed-field gel electrophoresis (PFGE), using SmaI or Cfr9I and SfiI. A subset of 41 isolates, representing each PFGE cluster, was also characterized by multilocus sequence typing (MLST). The application of Adjusted Rand and Wallace indices allowed the evaluation of the strength and the directionality of the correspondences between the various typing methods and showed that if PFGE or MLST data are available one can confidently predict the emm type (Wallace coefficients of 0.952 for both methods). In contrast, emm typing was a poor predictor of PFGE cluster or MLST sequence type (Wallace coefficients of 0.803 and 0.655, respectively). This was confirmed by the analysis of the larger data set available from http://spyogenes.mlst.net and underscores the necessity of performing PFGE or MLST to unambiguously define clones in S. pyogenes.

Natural history of drug-resistant clones of Streptococcus pneumoniae colonizing healthy children in Portugal.

A total of 3,539 Streptococcus pneumoniae (Pn) were recovered from 4,969 nasopharyngeal samples of children attending 13 day-care centers (DCCs) located in Lisbon, Portugal, during a surveillance study from January, 2001, through March, 2003, integrated in the European intervention project (EURIS, European Resistance Intervention Study). All Pn isolates were tested for anti-biotyping and drug-resistant pneumococci (DRPn) were further tested by serotyping and pulsed-field gel electrophoresis (PFGE). Overall carriage of Pn was very high (71.2%) and 39.9% of the isolates were resistant to antimicrobials (22.5% with decreased susceptibility to penicillin and 17.4% susceptible to penicillin and resistant to other antimicrobials). Serotypes 6B, 14, 23 F, 19F, and 19 A were prevalent among the 1,287 DRPn and 5.8% of the isolates were non-typeable. Eighty PFGE patterns were identified among 1,285 DRPn, and 93.1% of the DRPn belonged to 26 major clonal types that comprised: Pneumococcal Molecular Epidemiology Network (PMEN) clones (76.3%), Portuguese (PT)-DCC clones, previously detected in 1996-1999 (14.3%), and EURIS PT-DCC new clones, identified for the first time in the EURIS study, during 2001-2003 (9.4%). Comparing with previous Portuguese surveillance studies carried out since 1996, we observed that carriage increased from 47% to 71%, but no major changes were detected on the prevalence of pneumococcal serotypes. Moreover, although PMEN clones were predominant in all DCCs, in the present study the majority of them were gradually decreasing in time whereas several PT-DCC and new clones seemed to be increasing.

Assessment of band-based similarity coefficients for automatic type and subtype classification of microbial isolates analyzed by pulsed-field gel electrophoresis.

Pulsed-field gel electrophoresis (PFGE) has been the typing method of choice for strain identification in epidemiological studies of several bacterial species of medical importance. The usual procedure for the comparison of strains and assignment of strain type and subtype relies on visual assessment of band difference number, followed by an incremental assignment to the group hosting the most similar type previously seen. Band-based similarity coefficients, such as the Dice or the Jaccard coefficient, are then used for dendrogram construction, which provides a quantitative assessment of strain similarity. PFGE type assignment is based on the definition of a threshold linkage value, below which strains are assigned to the same group. This is typically performed empirically by inspecting the hierarchical cluster analysis dendrogram containing the strains of interest. This approach has the problem that the threshold value selected is dependent on the linkage method used for dendrogram construction. Furthermore, the use of a linkage method skews the original similarity values between strains. In this paper we assess the goodness of classification of several band-based similarity coefficients by comparing it with the band difference number for PFGE type and subtype classification using receiver operating characteristic curves. The procedure described was applied to a collection of PFGE results for 1,798 isolates of Streptococcus pneumoniae, which documented 96 types and 396 subtypes. The band-based similarity coefficients were found to perform equally well for type classification, but with different proportions of false-positive and false-negative classifications in their minimal false discovery rate when they were used for subtype classification.

Characterization of the genetic lineages responsible for pneumococcal invasive disease in Portugal.

The availability of a conjugate vaccine has the potential to reduce the disease burden of pneumococci and to alter the serotype frequency in the disease-causing population through immunoselection. These changes will probably be reflected in the distributions of individual genetic lineages within the population. We present a characterization of a collection of recent (1999 to 2002) invasive isolates from Portugal (n = 465) by macrorestriction profiling with pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing. During this time, serotypes 14, 1, 3, 4, 8, 9V, 23F, 7F, 19A, and 12B were the 10 most prevalent overall by decreasing rank order. By combining the PFGE data with the sequence types (STs) of 104 isolates, we were able to identify the genetic lineages of the majority of the isolates. We found 66 STs, including 20 novel STs, corresponding to 47 different lineages by e-BURST analysis. We found in our collection a number of previously identified internationally disseminated lineages, especially among macrolide-resistant and penicillin-resistant isolates, and these accounted for most of the isolates. Most of the major lineages (17 of 25) were identified in all years of the study, suggesting that the pneumococcal population associated with invasive disease was stable. This study provides a characterization of the pneumococcal population associated with invasive disease that will be useful for detecting potential selective effects of the novel conjugate vaccine.

Trends in drug resistance, serotypes, and molecular types of Streptococcus pneumoniae colonizing preschool-age children attending day care centers in Lisbon, Portugal: a summary of 4 years of annual surveillance.

Of the nasopharyngeal cultures recovered from 942 day care center (DCC) attendees in Lisbon, Portugal, 591 (62%) yielded Streptococcus pneumoniae during a surveillance performed in February and March of 1999. Forty percent of the isolates were resistant to one or more antimicrobial agents. In particular, 2% were penicillin resistant and 20% had intermediate penicillin resistance. Multidrug resistance to macrolides, lincosamides, and tetracycline was the most frequent antibiotype (17% of all isolates). Serotyping and molecular typing by pulsed-field gel electrophoresis were performed for 202 out of 237 drug-resistant pneumococci (DRPn). The most frequent serotypes were 6B (26%), 14 (22%), 19F (16%), 23F (10%), and nontypeable (12%). The majority (67%) of the DRPn strains were representatives of nine international clones included in the Pneumococcal Molecular Epidemiology Network; eight of them had been detected in previous studies. Fourteen novel clones were identified, corresponding to 26% of the DRPn strains. The remaining 7% of the strains were local clones detected in our previous studies. Comparison with studies conducted since 1996 in Portuguese DCCs identified several trends: (i) the rate of DRPn frequency has fluctuated between 40 and 50%; (ii) the serotypes most frequently recovered have remained the same; (iii) nontypeable strains appear to be increasing in frequency; and (iv) a clone of serotype 33F emerged in 1999. Together, our observations highlight that the nasopharynxes of children in DCCs are a melting pot of successful DRPn clones that are important to study and monitor if we aim to gain a better understanding on the epidemiology of this pathogen.

Analysis of genomic sequences by Chaos Game Representation.

MOTIVATION: Chaos Game Representation (CGR) is an iterative mapping technique that processes sequences of units, such as nucleotides in a DNA sequence or amino acids in a protein, in order to find the coordinates for their position in a continuous space. This distribution of positions has two properties: it is unique, and the source sequence can be recovered from the coordinates such that distance between positions measures similarity between the corresponding sequences. The possibility of using the latter property to identify succession schemes have been entirely overlooked in previous studies which raises the possibility that CGR may be upgraded from a mere representation technique to a sequence modeling tool. RESULTS: The distribution of positions in the CGR plane were shown to be a generalization of Markov chain probability tables that accommodates non-integer orders. Therefore, Markov models are particular cases of CGR models rather than the reverse, as currently accepted. In addition, the CGR generalization has both practical (computational efficiency) and fundamental (scale independence) advantages. These results are illustrated by using Escherichia coli K-12 as a test data-set, in particular, the genes thrA, thrB and thrC of the threonine operon.

Simultaneous transurethral prostatectomy and herniorrhaphy. A study of 63 cases.

Transurethral prostatectomy and herniorrhaphy were performed as a combined procedure for 63 patients during the period 1974 to 1978. The intraoperative and postoperative courses are compared with those of a control group of 82 patients who underwent herniorrhapies as an isolated procedure during the same period. The simultaneous performance of transurethral prostatectomy and herniorrhaphy in this group of patients was safe as well as cost-effective.