Investigations into genome diversity of Haemophilus influenzae using whole genome sequencing of clinical isolates and laboratory transformants

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Multiple interspecies recombination events documented by whole-genome sequencing in multidrug-resistant Haemophilus influenzae clinical isolates

Introduction

Haemophilus influenzae (Hi) was long known as an easy-to-treat bacterium, but increasing resistance against beta-lactams and other critically important antibiotics is now a growing concern. We describe here the whole-genome sequencing (WGS) analysis of three non-typeable Hi isolates received in 2018-2019 by the Belgian National Reference Centre (NRC) for Haemophilus influenzae, as they presented an unusual multi-resistant profile.

Methods

All three isolates were sequenced by WGS and mapped to the reference isolate Hi Rd KW20. Shorten uptake signal sequences (USSs) known to be associated with homologous recombination were sought in ftsI, murE and murF genes, and inner partial sequences were compared against the blast nucleotide database to look for similarity with other Haemophilus species. Their antimicrobial resistance (AMR) genotype was studied. Core-genome multilocus sequence typing (MLST) was performed on the NTHi database pubMLST to place our isolates in the actual worldwide epidemiology.

Results

The isolates also harboured interspecies recombination patterns in the murF-murE-ftsI region involved in cell wall synthesis. The three isolates were multidrug resistant and two of them were also resistant to amoxicillin-clavulanic acid and showed a reduced susceptibility to meropenem. All three isolates belonged to the MLST clonal complex (CC) 422, and WGS revealed that the three were very similar. They harboured mobile genetic elements (carrying blaTEM-1B, mefA and msrD genes associated with resistance), mutations in gyrA and parC linked to fluoroquinolone resistance as well as remodelling events in ompP2 that might be related to lower carbapenem susceptibility.

Conclusion

The Hi evolution towards antimicrobial multiresistance (AMR) is a complex and poorly understood phenomenon, although probably linked to a large degree to the presence of USSs and exchange within the family Pasteurellaceae. To better understand the respective roles of clonal expansion, horizontal gene transfers, spontaneous mutations and interspecies genetic rearrangements in shaping Hi AMR, both analysis of Hi communities over time within individuals and worldwide monitoring of non-typeable Hi causing infections should be conducted.

Identifying Bacterial Airways Infection in Stable Severe Asthma Using Oxford Nanopore Sequencing Technologies

Previous metagenomic studies in asthma have been limited by inadequate sequencing depth for species-level bacterial identification and by heterogeneity in clinical phenotyping. We hypothesize that chronic bacterial airways infection is a key "treatable trait" whose prevalence, clinical phenotype and reliable biomarkers need definition. In this study, we have applied a method for Oxford Nanopore sequencing for the unbiased metagenomic characterization of severe asthma. We optimized methods to compare performance of Illumina MiSeq, Nanopore sequencing, and RT-qPCR on total sputum DNA extracts against culture/MALDI-TOF for analysis of induced sputum samples from highly phenotyped severe asthma during clinical stability. In participants with severe asthma (n = 23) H. influenzae was commonly cultured (n = 8) and identified as the dominant bacterial species by metagenomic sequencing using an optimized method for Illumina MiSeq and Oxford Nanopore. Alongside superior operational characteristics, Oxford Nanopore achieved near complete genome coverage of H. influenzae and demonstrated a high level of agreement with Illumina MiSeq data. Clinically significant infection was confirmed with validated H. influenzae plasmid-based quantitative PCR assay. H. influenzae positive patients were found to have sputum neutrophilia and lower FeNO. In conclusion, using an optimized method of direct sequencing of induced sputum samples, H. influenzae was identified as a clinically relevant pathogen in severe asthma and was identified reliably using metagenomic sequencing. Application of these protocols in ongoing analysis of large patient cohorts will allow full characterization of this clinical phenotype. IMPORTANCE The human airways were once thought sterile in health. Now metagenomic techniques suggest bacteria may be present, but their role in asthma is not understood. Traditional culture lacks sensitivity and current sequencing techniques are limited by operational problems and limited ability to identify pathogens at species level. We optimized a new sequencing technique-Oxford Nanopore technologies (ONT)-for use on human sputum samples and compared it with existing methods. We found ONT was effective for rapidly analyzing samples and could identify bacteria at the species level. We used this to show Haemophilus influenzae was a dominant bacterium in the airways in people with severe asthma. The presence of Haemophilus was associated with a "neutrophilic" form of asthma - a subgroup for which we currently lack specific treatments. Therefore, this technique could be used to target chronic antibiotic therapy and in research to characterize the full breadth of bacteria in the airways.

Comparative pangenome analysis of capsulated Haemophilus influenzae serotype f highlights their high genomic stability

Haemophilus influenzae is an opportunistic pathogen adapted to the human respiratory tract. Non-typeable H. influenzae are highly heterogeneous, but few studies have analysed the genomic variability of capsulated strains. This study aims to examine the genetic diversity of 37 serotype f isolates from the Netherlands, Portugal, and Spain, and to compare all capsulated genomes available on public databases. Serotype f isolates belonged to CC124 and shared few single nucleotide polymorphisms (SNPs) (n = 10,999), but a high core genome (> 80%). Three main clades were identified by the presence of 75, 60 and 41 exclusive genes for each clade, respectively. Multi-locus sequence type analysis of all capsulated genomes revealed a reduced number of clonal complexes associated with each serotype. Pangenome analysis showed a large pool of genes (n = 6360), many of which were accessory genome (n = 5323). Phylogenetic analysis revealed that serotypes a, b, and f had greater diversity. The total number of SNPs in serotype f was significantly lower than in serotypes a, b, and e (p < 0.0001), indicating low variability within the serotype f clonal complexes. Capsulated H. influenzae are genetically homogeneous, with few lineages in each serotype. Serotype f has high genetic stability regardless of time and country of isolation.

Learning from -omics strategies applied to uncover Haemophilus influenzae host-pathogen interactions: Current status and perspectives

Haemophilus influenzae has contributed to key bacterial genome sequencing hallmarks, as being not only the first bacterium to be genome-sequenced, but also starring the first genome-wide analysis of chromosomes directly transformed with DNA from a divergent genotype, and pioneering Tn-seq methodologies. Over the years, the phenomenal and constantly evolving development of -omic technologies applied to a whole range of biological questions of clinical relevance in the H. influenzae-host interplay, has greatly moved forward our understanding of this human-adapted pathogen, responsible for multiple acute and chronic infections of the respiratory tract. In this way, essential genes, virulence factors, pathoadaptive traits, and multi-layer gene expression regulatory networks with both genomic and epigenomic complexity levels are being elucidated. Likewise, the unstoppable increasing whole genome sequencing information underpinning H. influenzae great genomic plasticity, mainly when referring to non-capsulated strains, poses major challenges to understand the genomic basis of clinically relevant phenotypes and even more, to clearly highlight potential targets of clinical interest for diagnostic, therapeutic or vaccine development. We review here how genomic, transcriptomic, proteomic and metabolomic-based approaches are great contributors to our current understanding of the interactions between H. influenzae and the human airways, and point possible strategies to maximize their usefulness in the context of biomedical research and clinical needs on this human-adapted bacterial pathogen.

Whole-genome analyses reveal gene content differences between nontypeable Haemophilus influenzae isolates from chronic obstructive pulmonary disease compared to other clinical phenotypes

Nontypeable Haemophilus influenzae (NTHi) colonizes human upper respiratory airways and plays a key role in the course and pathogenesis of acute exacerbations of chronic obstructive pulmonary disease (COPD). Currently, it is not possible to distinguish COPD isolates of NTHi from other clinical isolates of NTHi using conventional genotyping methods. Here, we analysed the core and accessory genome of 568 NTHi isolates, including 40 newly sequenced isolates, to look for genetic distinctions between NTHi isolates from COPD with respect to other illnesses, including otitis media, meningitis and pneumonia. Phylogenies based on polymorphic sites in the core-genome did not show discrimination between NTHi strains collected from different clinical phenotypes. However, pan-genome-wide association studies identified 79 unique NTHi accessory genes that were significantly associated with COPD. Furthermore, many of the COPD-related NTHi genes have known or predicted roles in virulence, transmembrane transport of metal ions and nutrients, cellular respiration and maintenance of redox homeostasis. This indicates that specific genes may be required by NTHi for its survival or virulence in the COPD lung. These results advance our understanding of the pathogenesis of NTHi infection in COPD lungs.

Characterization of biofilm formation and induction of apoptotic DNA fragmentation by nontypeable Haemophilus influenzae on polarized human airway epithelial cells

Nontypeable Haemophilus influenzae (NTHi) is a common airway commensal and opportunistic pathogen that persists within biofilm communities in vivo. Biofilm studies so far are mainly based on assays on plastic surfaces. The aim of this work was to investigate the capacity of clinical NTHi strains to form biofilm structures on polarized Calu-3 human airway epithelial cells and primary normal human bronchial epithelial cells and to characterize the biofilm architecture. Formation of adherent NTHi biofilms post colonization of host cells at multiple time-points was evaluated using confocal laser scanning microscopy and electron microscopy. NTHi biofilms were analyzed in terms of biofilm height and presence of extracellular matrix components, and their apoptotic effects on epithelial cells were measured by TUNEL assay. Strain Fi176 was observed to form robust biofilms on airway epithelia over time, while disrupting the integrity of Calu-3 monolayer by 72 h of co-culture. NTHi biofilms were observed to induce apoptotic DNA fragmentation in host cells at 24 h post infection. Biofilm formation on cell monolayers by Fi176ΔpilA strain was markedly reduced compared to WT strain. Biofilm inhibition and disruption assays by crystal violet staining indicated that DNA and proteins are part of NTHi biofilms in vitro. Our findings highlight critical stages of NTHi pathogenesis following host colonization and provide useful biofilm models for future antimicrobial drug discovery investigations.

Synergistic Activity of Mobile Genetic Element Defences in Streptococcus pneumoniae

A diverse set of mobile genetic elements (MGEs) transmit between Streptococcus pneumoniae cells, but many isolates remain uninfected. The best-characterised defences against horizontal transmission of MGEs are restriction-modification systems (RMSs), of which there are two phase-variable examples in S. pneumoniae. Additionally, the transformation machinery has been proposed to limit vertical transmission of chromosomally integrated MGEs. This work describes how these mechanisms can act in concert. Experimental data demonstrate RMS phase variation occurs at a sub-maximal rate. Simulations suggest this may be optimal if MGEs are sometimes vertically inherited, as it reduces the probability that an infected cell will switch between RMS variants while the MGE is invading the population, and thereby undermine the restriction barrier. Such vertically inherited MGEs can be deleted by transformation. The lack of between-strain transformation hotspots at known prophage att sites suggests transformation cannot remove an MGE from a strain in which it is fixed. However, simulations confirmed that transformation was nevertheless effective at preventing the spread of MGEs into a previously uninfected cell population, if a recombination barrier existed between co-colonising strains. Further simulations combining these effects of phase variable RMSs and transformation found they synergistically inhibited MGEs spreading, through limiting both vertical and horizontal transmission.

Molecular Signatures of Non-typeable Haemophilus influenzae Lung Adaptation in Pediatric Chronic Lung Disease

Non-typeable Haemophilus influenzae (NTHi), an opportunistic pathogen of the upper airways of healthy children, can infect the lower airways, driving chronic lung disease. However, the molecular basis underpinning NTHi transition from a commensal to a pathogen is not clearly understood. Here, we performed comparative genomic and transcriptomic analyses of 12 paired, isogenic NTHi strains, isolated from the nasopharynx (NP) and bronchoalveolar lavage (BAL) of 11 children with chronic lung disease, to identify convergent molecular signatures associated with lung adaptation. Comparative genomic analyses of the 12 NP-BAL pairs demonstrated that five were genetically identical, with the remaining seven differing by only 1 to 3 mutations. Within-patient transcriptomic analyses identified between 2 and 58 differentially expressed genes in 8 of the 12 NP-BAL pairs, including pairs with no observable genomic changes. Whilst no convergence was observed at the gene level, functional enrichment analysis revealed significant under-representation of differentially expressed genes belonging to Coenzyme metabolism, Function unknown, Translation, ribosomal structure, and biogenesis Cluster of Orthologous Groups categories. In contrast, Carbohydrate transport and metabolism, Cell motility and secretion, Intracellular trafficking and secretion, and Energy production categories were over-represented. This observed trend amongst genetically unrelated NTHi strains provides evidence of convergent transcriptional adaptation of NTHi to pediatric airways that deserves further exploration. Understanding the pathoadaptative mechanisms that NTHi employs to infect and persist in the lower pediatric airways is essential for devising targeted diagnostics and treatments aimed at minimizing disease severity, and ultimately, preventing NTHi lung infections and subsequent chronic lung disease in children.

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.

Patients with Chronic Obstructive Pulmonary Disease harbour a variation of Haemophilus species

H. haemolyticus is often misidentified as NTHi due to their close phylogenetic relationship. Differentiating between the two is important for correct identification and appropriate treatment of infective organism and to ensure any role of H. haemolyticus in disease is not being overlooked. Speciation however is not completely reliable by culture and PCR methods due to the loss of haemolysis by H. haemolyticus and the heterogeneity of NTHi. Haemophilus isolates from COPD as part of the AERIS study (ClinicalTrials - NCT01360398) were speciated by analysing sequence data for the presence of molecular markers. Further investigation into the genomic relationship was carried out using average nucleotide identity and phylogeny of allelic and genome alignments. Only 6.3% were identified as H. haemolyticus. Multiple in silico methods were able to distinguish H. haemolyticus from NTHi. However, no single gene target was found to be 100% accurate. A group of omp2 negative NTHi were observed to be phylogenetically divergent from H. haemolyticus and remaining NTHi. The presence of an atypical group from a geographically and disease limited set of isolates supports the theory that the heterogeneity of NTHi may provide a genetic continuum between NTHi and H. haemolyticus.

The return of Pfeiffer's bacillus: Rising incidence of ampicillin resistance in Haemophilus influenzae

Haemophilus influenzae, originally named Pfeiffer's bacillus after its discoverer Richard Pfeiffer in 1892, was a major risk for global health at the beginning of the 20th century, causing childhood pneumonia and invasive disease as well as otitis media and other upper respiratory tract infections. The implementation of the Hib vaccine, targeting the major capsule type of H. influenzae, almost eradicated the disease in countries that adapted the vaccination scheme. However, a rising number of infections are caused by non-typeable H. influenzae (NTHi), which has no capsule and against which the vaccine therefore provides no protection, as well as other serotypes equally not recognised by the vaccine. The first line of treatment is ampicillin, but there is a steady rise in ampicillin resistance. This is both through acquired as well as intrinsic mechanisms, and is cause for serious concern and the need for more surveillance. There are also increasing reports of new modifications of the intrinsic ampicillin-resistance mechanism leading to resistance against cephalosporins and carbapenems, the last line of well-tolerated drugs, and ampicillin-resistant H. influenzae was included in the recently released priority list of antibiotic-resistant bacteria by the WHO. This review provides an overview of ampicillin resistance prevalence and mechanisms in the context of our current knowledge about population dynamics of H. influenzae.

Haemophilus influenzae genome evolution during persistence in the human airways in chronic obstructive pulmonary disease

Nontypeable Haemophilus influenzae (NTHi) exclusively colonize and infect humans and are critical to the pathogenesis of chronic obstructive pulmonary disease (COPD). In vitro and animal models do not accurately capture the complex environments encountered by NTHi during human infection. We conducted whole-genome sequencing of 269 longitudinally collected cleared and persistent NTHi from a 15-y prospective study of adults with COPD. Genome sequences were used to elucidate the phylogeny of NTHi isolates, identify genomic changes that occur with persistence in the human airways, and evaluate the effect of selective pressure on 12 candidate vaccine antigens. Strains persisted in individuals with COPD for as long as 1,422 d. Slipped-strand mispairing, mediated by changes in simple sequence repeats in multiple genes during persistence, regulates expression of critical virulence functions, including adherence, nutrient uptake, and modification of surface molecules, and is a major mechanism for survival in the hostile environment of the human airways. A subset of strains underwent a large 400-kb inversion during persistence. NTHi does not undergo significant gene gain or loss during persistence, in contrast to other persistent respiratory tract pathogens. Amino acid sequence changes occurred in 8 of 12 candidate vaccine antigens during persistence, an observation with important implications for vaccine development. These results indicate that NTHi alters its genome during persistence by regulation of critical virulence functions primarily by slipped-strand mispairing, advancing our understanding of how a bacterial pathogen that plays a critical role in COPD adapts to survival in the human respiratory tract.

The adhesins of non-typeable Haemophilus influenzae

INTRODUCTION

Nontypeable Haemophilus influenzae (NTHi) is an opportunistic pathogen of the respiratory tract and the greatest contributor to invasive Haemophilus disease. Additionally, in children, NTHi is responsible for the majority of otitis media (OM) which can lead to chronic infection and hearing loss. In adults, NTHi infection in the lungs is responsible for the onset of acute exacerbations in chronic obstructive pulmonary disease (COPD). Unfortunately, there is currently no vaccine available to protect against NTHi infections. Areas covered: NTHi uses an arsenal of adhesins to colonise the respiratory epithelium. The adhesins also have secondary roles that aid in the virulence of NTHi, including mechanisms that avoid immune clearance, adjust pore size to avoid antimicrobial destruction, form micro-colonies and invoke phase variation for protein mediation. Bacterial adhesins can also be ideal antigens for subunit vaccine design due to surface exposure and immunogenic capabilities. Expert commentary: The host-pathogen interactions of the NTHi adhesins are not fully investigated. The relationship between adhesins and the extracellular matrix (ECM) play a part in the success of NTHi colonisation and virulence by immune evasion, migration and biofilm development. Further research into these immunogenic proteins would further our understanding and enable a basis for better combatting NTHi disease.

Simultaneous identification of Haemophilus influenzae and Haemophilus haemolyticus using real-time PCR

AIM

To design a highly specific and sensitive multiplex real-time PCR assay for the differentiation of the pathogen Haemophilus influenzae from its nonpathogenic near-neighbor Haemophilus haemolyticus.

MATERIALS & METHODS

A comparison of 380 Haemophilus spp. genomes was used to identify loci specific for each species. Novel PCR assays targeting H. haemolyticus (hypD) and H. influenzae (siaT) were designed.

RESULTS & DISCUSSION

PCR screening across 143 isolates demonstrated 100% specificity for hypD and siaT. These two assays were multiplexed with the recently described fucP assay for further differentiation among H. influenzae.

CONCLUSION

The triplex assay provides rapid, unambiguous, sensitive and highly specific genotyping results for the simultaneous detection of hypD and siaT, including fucose-positive H. influenzae (fucP), in a single PCR.

Virulence factor activity relationships (VFARs): a bioinformatics perspective

Virulence factor activity relationships (VFARs) - a concept loosely based on quantitative structure-activity relationships (QSARs) for chemicals was proposed as a predictive tool for ranking risks due to microorganisms relevant to water safety. A rapid increase in sequencing capabilities and bioinformatics tools has significantly increased the potential for VFAR-based analyses. This review summarizes more than 20 bioinformatics databases and tools, developed over the last decade, along with their virulence and antimicrobial resistance prediction capabilities. With the number of bacterial whole genome sequences exceeding 241 000 and metagenomic analysis projects exceeding 13 000 and the ability to add additional genome sequences for few hundred dollars, it is evident that further development of VFARs is not limited by the availability of information at least at the genomic level. However, additional information related to co-occurrence, treatment response, modulation of virulence due to environmental and other factors, and economic impact must be gathered and incorporated in a manner that also addresses the associated uncertainties. Of the bioinformatics tools, a majority are either designed exclusively for virulence/resistance determination or equipped with a dedicated module. The remaining have the potential to be employed for evaluating virulence. This review focusing broadly on omics technologies and tools supports the notion that these tools are now sufficiently developed to allow the application of VFAR approaches combined with additional engineering and economic analyses to rank and prioritize organisms important to a given niche. Knowledge gaps do exist but can be filled with focused experimental and theoretical analyses that were unimaginable a decade ago. Further developments should consider the integration of the measurement of activity, risk, and uncertainty to improve the current capabilities.

3D Reconstruction of the Human Airway Mucosa In Vitro as an Experimental Model to Study NTHi Infections

We have established an in vitro 3D system which recapitulates the human tracheo-bronchial mucosa comprehensive of the pseudostratified epithelium and the underlying stromal tissue. In particular, we reported that the mature model, entirely constituted of primary cells of human origin, develops key markers proper of the native tissue such as the mucociliary differentiation of the epithelial sheet and the formation of the basement membrane. The infection of the pseudo-tissue with a strain of NonTypeable Haemophilus influenzae results in bacteria association and crossing of the mucus layer leading to an apparent targeting of the stromal space where they release large amounts of vesicles and form macro-structures. In summary, we propose our in vitro model as a reliable and potentially customizable system to study mid/long term host-pathogen processes.

Targeted single molecule mutation detection with massively parallel sequencing

Next-generation sequencing (NGS) technologies have transformed genomic research and have the potential to revolutionize clinical medicine. However, the background error rates of sequencing instruments and limitations in targeted read coverage have precluded the detection of rare DNA sequence variants by NGS. Here we describe a method, termed CypherSeq, which combines double-stranded barcoding error correction and rolling circle amplification (RCA)-based target enrichment to vastly improve NGS-based rare variant detection. The CypherSeq methodology involves the ligation of sample DNA into circular vectors, which contain double-stranded barcodes for computational error correction and adapters for library preparation and sequencing. CypherSeq is capable of detecting rare mutations genome-wide as well as those within specific target genes via RCA-based enrichment. We demonstrate that CypherSeq is capable of correcting errors incurred during library preparation and sequencing to reproducibly detect mutations down to a frequency of 2.4 × 10⁻⁷ per base pair, and report the frequency and spectra of spontaneous and ethyl methanesulfonate-induced mutations across the Saccharomycescerevisiae genome.

Expression of IgA Proteases by Haemophilus influenzae in the Respiratory Tract of Adults With Chronic Obstructive Pulmonary Disease

BACKGROUND

Immunoglobulin (Ig)A proteases of Haemophilus influenzae are highly specific endopeptidases that cleave the hinge region of human IgA1 and also mediate invasion and trafficking in human respiratory epithelial cells, facilitating persistence of H. influenzae. Little is known about the expression of IgA proteases in clinical settings of H. influenzae infection.

METHODS

We identified and characterized IgA protease genes in H. influenzae and studied their expression and proteolytic specificity, in vitro and in vivo in 169 independent strains of H. influenzae collected longitudinally over 10 years from adults with chronic obstructive pulmonary disease.

RESULTS

The H. influenzae pangenome has 2 alleles of IgA protease genes; all strains have igaA, and 40% of strains have igaB. Each allele has 2 variants with differing proteolytic specificities for human IgA1. A total of 88% of 169 strains express IgA protease activity. Expression of the 4 forms of IgA protease varies among strains. Based on the presence of IgA1 fragments in sputum samples, each of the different forms of IgA protease is selectively expressed in the human airways during infection.

CONCLUSIONS

Four variants of IgA proteases are variably expressed by H. influenzae during infection of the human airways.

Comparative whole-genome analysis of clinical isolates reveals characteristic architecture of Mycobacterium tuberculosis pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.

Vaccines for Nontypeable Haemophilus influenzae: the Future Is Now

Infections due to nontypeable Haemophilus influenzae result in enormous global morbidity in two clinical settings: otitis media in children and respiratory tract infections in adults with chronic obstructive pulmonary disease (COPD). Recurrent otitis media affects up to 20% of children and results in hearing loss, delays in speech and language development and, in developing countries, chronic suppurative otitis media. Infections in people with COPD result in clinic and emergency room visits, hospital admissions, and respiratory failure. An effective vaccine would prevent morbidity, help control health care costs, and reduce antibiotic use, a major contributor to the global crisis in bacterial antibiotic resistance. The widespread use of the pneumococcal conjugate vaccines is causing a relative increase in H. influenzae otitis media. The partial protection against H. influenzae otitis media induced by the pneumococcal H. influenzae protein D conjugate vaccine represents a proof of principle of the feasibility of a vaccine for nontypeable H. influenzae. An ideal vaccine antigen should be conserved among strains, have abundant epitopes on the bacterial surface, be immunogenic, and induce protective immune responses. Several surface proteins of H. influenzae have been identified as potential vaccine candidates and are in various stages of development. With continued research, progress toward a broadly effective vaccine to prevent infections caused by nontypeable H. influenzae is expected over the next several years.

Haemophilus influenzae: recent advances in the understanding of molecular pathogenesis and polymicrobial infections

PURPOSE OF REVIEW

Non-typeable Haemophilus influenzae (NTHi) is a human-specific mucosal pathogen and one of the most common causes of bacterial infections in children and patients with chronic obstructive pulmonary disease. It is also frequently found in polymicrobial superinfections. Great strides have recently been made in the understanding of the molecular mechanisms underlying NTHi pathogenesis.

RECENT FINDINGS

By using new methodology, such as experimental human colonization models and whole-genome approaches, investigators have shed light upon the various strategies of NTHi that are involved in pathogenesis. These include the escape of the mucociliary elevator, evasion of host immunity, survival in environments with scarce nutrients, and finally participation in polymicrobial infections. Lipooligosaccharide branching, proteinous adhesins, metabolic adaption to nutrient availability and many scavenging systems are implicated in these processes. Interestingly, genome-based studies comparing virulent and commensal strains have identified many hypothetical proteins as virulence determinants, suggesting that much regarding the molecular pathogenesis of NTHi remains to be solved.

SUMMARY

NTHi is an opportunistic pathogen and highly specialized colonizer of the human respiratory tract that has developed intricate mechanisms to establish growth and survival in the human host. Continued research is needed to further elucidate NTHi host-pathogen and pathogen-pathogen interactions.

Multilocus sequence typing and ftsI sequencing: a powerful tool for surveillance of penicillin-binding protein 3-mediated beta-lactam resistance in nontypeable Haemophilus influenzae

Background

Beta-lactam resistance in Haemophilus influenzae due to ftsI mutations causing altered penicillin-binding protein 3 (PBP3) is increasing worldwide. Low-level resistant isolates with the N526K substitution (group II low-rPBP3) predominate in most geographical regions, while high-level resistant isolates with the additional S385T substitution (group III high-rPBP3) are common in Japan and South Korea.

Knowledge about the molecular epidemiology of rPBP3 strains is limited. We combined multilocus sequence typing (MLST) and ftsI/PBP3 typing to study the emergence and spread of rPBP3 in nontypeable H. influenzae (NTHi) in Norway.

Results

The prevalence of rPBP3 in a population of 795 eye, ear and respiratory isolates (99% NTHi) from 2007 was 15%. The prevalence of clinical PBP3-mediated resistance to ampicillin was 9%, compared to 2.5% three years earlier. Group II low-rPBP3 predominated (96%), with significant proportions of isolates non-susceptible to cefotaxime (6%) and meropenem (20%). Group III high-rPBP3 was identified for the first time in Northern Europe.

Four MLST sequence types (ST) with characteristic, highly diverging ftsI alleles accounted for 61% of the rPBP3 isolates. The most prevalent substitution pattern (PBP3 type A) was present in 41% of rPBP3 isolates, mainly carried by ST367 and ST14. Several unrelated STs possessed identical copies of the ftsI allele encoding PBP3 type A.

Infection sites, age groups, hospitalization rates and rPBP3 frequencies differed between STs and phylogenetic groups.

Conclusions

This study is the first to link ftsI alleles to STs in H. influenzae. The results indicate that horizontal gene transfer contributes to the emergence of rPBP3 by phylogeny restricted transformation.

Clonally related virulent rPBP3 strains are widely disseminated and high-level resistant isolates emerge in new geographical regions, threatening current empiric antibiotic treatment. The need of continuous monitoring of beta-lactam susceptibility and a global system for molecular surveillance of rPBP3 strains is underlined. Combining MLST and ftsI/PBP3 typing is a powerful tool for this purpose.

Extensive cotransformation of natural variation into chromosomes of naturally competent Haemophilus influenzae

Naturally competent bacterial species actively take up environmental DNA and can incorporate it into their chromosomes by homologous recombination. This can bring genetic variation from environmental DNA to recipient chromosomes, often in multiple long “donor” segments. Here, we report the results of genome sequencing 96 colonies of a laboratory Haemophilus influenzae strain, which had been experimentally transformed by DNA from a diverged clinical isolate. Donor segments averaged 6.9 kb (spanning several genes) and were clustered into recombination tracts of ~19.5 kb. Individual colonies had replaced from 0.1 to 3.2% of their chromosomes, and ~1/3 of all donor-specific single-nucleotide variants were present in at least one recombinant. We found that nucleotide divergence did not obviously limit the locations of recombination tracts, although there were small but significant reductions in divergence at recombination breakpoints. Although indels occasionally transformed as parts of longer recombination tracts, they were common at breakpoints, suggesting that indels typically block progression of strand exchange. Some colonies had recombination tracts in which variant positions contained mixtures of both donor and recipient alleles. These tracts were clustered around the origin of replication and were interpreted as the result of heteroduplex segregation in the original transformed cell. Finally, a pilot experiment demonstrated the utility of natural transformation for genetically dissecting natural phenotypic variation. We discuss our results in the context of the potential to merge experimental and population genetic approaches, giving a more holistic understanding of bacterial gene transfer.

Genome sequencing of disease and carriage isolates of nontypeable Haemophilus influenzae identifies discrete population structure

One of the main hurdles for the development of an effective and broadly protective vaccine against nonencapsulated isolates of Haemophilus influenzae (NTHi) lies in the genetic diversity of the species, which renders extremely difficult the identification of cross-protective candidate antigens. To assess whether a population structure of NTHi could be defined, we performed genome sequencing of a collection of diverse clinical isolates representative of both carriage and disease and of the diversity of the natural population. Analysis of the distribution of polymorphic sites in the core genome and of the composition of the accessory genome defined distinct evolutionary clades and supported a predominantly clonal evolution of NTHi, with the majority of genetic information transmitted vertically within lineages. A correlation between the population structure and the presence of selected surface-associated proteins and lipooligosaccharide structure, known to contribute to virulence, was found. This high-resolution, genome-based population structure of NTHi provides the foundation to obtain a better understanding, of NTHi adaptation to the host as well as its commensal and virulence behavior, that could facilitate intervention strategies against disease caused by this important human pathogen.

Comparative genomic analysis reveals distinct genotypic features of the emerging pathogen Haemophilus influenzae type f

BACKGROUND

The incidence of invasive disease caused by encapsulated Haemophilus influenzae type f (Hif) has increased in the post-H. influenzae type b (Hib) vaccine era. We previously annotated the first complete Hif genome from a clinical isolate (KR494) that caused septic shock and necrotizing myositis. Here, the full genome of Hif KR494 was compared to sequenced reference strains Hib 10810, capsule type d (Hid) Rd Kw20, and finally nontypeable H. influenzae 3655. The goal was to identify possible genomic characteristics that may shed light upon the pathogenesis of Hif.

RESULTS

The Hif KR494 genome exhibited large regions of synteny with other H. influenzae, but also distinct genome rearrangements. A predicted Hif core genome of 1390 genes was shared with the reference strains, and 6 unique genomic regions comprising half of the 191 unique coding sequences were revealed. The majority of these regions were inserted genetic fragments, most likely derived from the closely-related Haemophilus spp. including H. aegyptius, H. haemolyticus and H. parainfluenzae. Importantly, the KR494 genome possessed several putative virulence genes that were distinct from non-type f strains. These included the sap2 operon, aef3 fimbriae, and genes for kanamycin nucleotidyltranserase, iron-utilization proteins, and putative YadA-like trimeric autotransporters that may increase the bacterial virulence. Furthermore, Hif KR494 lacked a hisABCDEFGH operon for de novo histidine biosynthesis, hmg locus for lipooligosaccharide biosynthesis and biofilm formation, the Haemophilus antibiotic resistance island and a Haemophilus secondary molybdate transport system. We confirmed the histidine auxotrophy and kanamycin resistance in Hif by functional experiments. Moreover, the pattern of unique or missing genes of Hif KR494 was similar in 20 Hif clinical isolates obtained from different years and geographical areas. A cross-species comparison revealed that the Hif genome shared more characteristics with H. aegyptius than Hid and NTHi.

CONCLUSIONS

The genomic comparative analyses facilitated identification of genotypic characteristics that may be related to the specific virulence of Hif. In relation to non-type f H. influenzae strains, the Hif genome contains differences in components involved in metabolism and survival that may contribute to its invasiveness.

Haemophilus parainfluenzae expresses diverse lipopolysaccharide O-antigens using ABC transporter and Wzy polymerase-dependent mechanisms

Lipopolysaccharide O-antigens are the basis of serotyping schemes for Gram negative bacteria and help to determine the nature of host–bacterial interactions. Haemophilus parainfluenzae is a normal commensal of humans but is also an occasional pathogen. The prevalence, diversity and biosynthesis of O-antigens were investigated in this species for the first time. 18/18 commensal H. parainfluenzae isolates contain a O-antigen biosynthesis gene cluster flanked by glnA and pepB, the same position as the hmg locus for tetrasaccharide biosynthesis in Haemophilus influenzae. The O-antigen loci show diverse restriction digest patterns but fall into two main groups: (1) those encoding enzymes for the synthesis and transfer of FucNAc4N in addition to the Wzy-dependent mechanism of O-antigen synthesis and transport and (2) those encoding galactofuranose synthesis/transfer enzymes and an ABC transporter. The other glycosyltransferase genes differ between isolates. Three H. parainfluenzae isolates fell outside these groups and are predicted to synthesise O-antigens containing ribitol phosphate or deoxytalose. Isolates using the ABC transporter system encode a putative O-antigen ligase, required for the synthesis of O-antigen-containing LPS glycoforms, at a separate genomic location. The presence of an O-antigen contributes significantly to H. parainfluenzae resistance to the killing effect of human serum in vitro. The discovery of O-antigens in H. parainfluenzae is striking, as its close relative H. influenzae lacks this cell surface component.