Mutation Analysis in Regulator DNA-Binding Regions for Antimicrobial Efflux Pumps in 17,000 Pseudomonas aeruginosa Genomes

Mutations leading to upregulation of efflux pumps can produce multiple drug resistance in the pathogen Pseudomonas aeruginosa. Changes in their DNA binding regions, i.e., palindromic operators, can compromise pump depression and subsequently enhance resistance against several antibacterials and biocides. Here, we have identified (pseudo)palindromic repeats close to promoters of genes encoding 13 core drug-efflux pumps of P. aeruginosa. This framework was applied to detect mutations in these repeats in 17,292 genomes. Eighty-nine percent of isolates carried at least one mutation. Eight binary genetic properties potentially related to expression were calculated for mutations. These included palindromicity reduction, mutation type, positioning within the repeat and DNA-bending shift. High-risk ST298, ST308 and ST357 clones commonly carried four conserved mutations while ST175 and the cystic fibrosis-linked ST649 clones showed none. Remarkably, a T-to-C transition in the fourth position of the upstream repeat for mexEF-oprN was nearly exclusive of the high-risk ST111 clone. Other mutations were associated with high-risk sublineages using sample geotemporal metadata. Moreover, 1.5% of isolates carried five or more mutations suggesting they undergo an alternative program for regulation of their effluxome. Overall, P. aeruginosa shows a wide range of operator mutations with a potential effect on efflux pump expression and antibiotic resistance.

StaR Is a Positive Regulator of Topoisomerase I Activity Involved in Supercoiling Maintenance in Streptococcus pneumoniae

The DNA topoisomerases gyrase and topoisomerase I as well as the nucleoid-associated protein HU maintain supercoiling levels in Streptococcus pneumoniae, a main human pathogen. Here, we characterized, for the first time, a topoisomerase I regulator protein (StaR). In the presence of sub-inhibitory novobiocin concentrations, which inhibit gyrase activity, higher doubling times were observed in a strain lacking staR, and in two strains in which StaR was over-expressed either under the control of the ZnSO₄-inducible P_Zn promoter (strain ΔstaRP_ZnstaR) or of the maltose-inducible P_Mal promoter (strain ΔstaRpLS1ROMstaR). These results suggest that StaR has a direct role in novobiocin susceptibility and that the StaR level needs to be maintained within a narrow range. Treatment of ΔstaRP_ZnstaR with inhibitory novobiocin concentrations resulted in a change of the negative DNA supercoiling density (σ) in vivo, which was higher in the absence of StaR (σ = -0.049) than when StaR was overproduced (σ = -0.045). We have located this protein in the nucleoid by using super-resolution confocal microscopy. Through in vitro activity assays, we demonstrated that StaR stimulates TopoI relaxation activity, while it has no effect on gyrase activity. Interaction between TopoI and StaR was detected both in vitro and in vivo by co-immunoprecipitation. No alteration of the transcriptome was associated with StaR amount variation. The results suggest that StaR is a new streptococcal nucleoid-associated protein that activates topoisomerase I activity by direct protein-protein interaction.

The challenges of the genome-based identification of antifungal resistance in the clinical routine

The increasing number of chronic and life-threatening infections caused by antimicrobial resistant fungal isolates is of critical concern. Low DNA sequencing cost may facilitate the identification of the genomic profile leading to resistance, the resistome, to rationally optimize the design of antifungal therapies. However, compared to bacteria, initiatives for resistome detection in eukaryotic pathogens are underdeveloped. Firstly, reported mutations in antifungal targets leading to reduced susceptibility must be extensively collected from the literature to generate comprehensive databases. This information should be complemented with specific laboratory screenings to detect the highest number possible of relevant genetic changes in primary targets and associations between resistance and other genomic markers. Strikingly, some drug resistant strains experience high-level genetic changes such as ploidy variation as much as duplications and reorganizations of specific chromosomes. Such variations involve allelic dominance, gene dosage increments and target expression regime effects that should be explicitly parameterized in antifungal resistome prediction algorithms. Clinical data indicate that predictors need to consider the precise pathogen species and drug levels of detail, instead of just genus and drug class. The concomitant needs for mutation accuracy and assembly quality assurance suggest hybrid sequencing approaches involving third-generation methods will be utilized. Moreover, fatal fast infections, like fungemia and meningitis, will further require both sequencing and analysis facilities are available in-house. Altogether, the complex nature of antifungal resistance demands extensive sequencing, data acquisition and processing, bioinformatic analysis pipelines, and standard protocols to be accomplished prior to genome-based protocols are applied in the clinical setting.

Unsupervised Machine Learning Organization of the Functional Dark Proteome of Gram-Negative "Superbugs": Six Protein Clusters Amenable for Distinct Scientific Applications

Uncharacterized proteins have been underutilized as targets for the development of novel therapeutics for difficult-to-treat bacterial infections. To facilitate the exploration of these proteins, 2819 predicted, uncharacterized proteins (19.1% of the total) from reference strains of multidrug Acinetobacter baumannii, Klebsiella pneumoniae, and Pseudomonas aeruginosa species were organized using an unsupervised k-means machine learning algorithm. Classification using normalized values for protein length, pI, hydrophobicity, degree of conservation, structural disorder, and %AT of the coding gene rendered six natural clusters. Cluster proteins showed different trends regarding operon membership, expression, presence of unknown function domains, and interactomic relevance. Clusters 2, 4, and 5 were enriched with highly disordered proteins, nonworkable membrane proteins, and likely spurious proteins, respectively. Clusters 1, 3, and 6 showed closer distances to known antigens, antibiotic targets, and virulence factors. Up to 21.8% of proteins in these clusters were structurally covered by modeling, which allowed assessment of druggability and discontinuous B-cell epitopes. Five proteins (4 in Cluster 1) were potential druggable targets for antibiotherapy. Eighteen proteins (11 in Cluster 6) were strong B-cell and T-cell immunogen candidates for vaccine development. Conclusively, we provide a feature-based schema to fractionate the functional dark proteome of critical pathogens for fundamental and biomedical purposes.

Association between HLA-C alleles and COVID-19 severity in a pilot study with a Spanish Mediterranean Caucasian cohort

The clinical presentations of COVID-19 may range from an asymptomatic or mild infection to a critical or fatal disease. Several host factors such as elderly age, male gender, and previous comorbidities seem to be involved in the most severe outcomes, but also an impaired immune response that causes a hyperinflammatory state but is unable to clear the infection. In order to get further understanding about this impaired immune response, we aimed to determine the association of specific HLA alleles with different clinical presentations of COVID-19. Therefore, we analyzed HLA Class I and II, as well as KIR gene sequences, in 72 individuals with Spanish Mediterranean Caucasian ethnicity who presented mild, severe, or critical COVID-19, according to their clinical characteristics and management. This cohort was recruited in Madrid (Spain) during the first and second pandemic waves between April and October 2020. There were no significant differences in HLA-A or HLA-B alleles among groups. However, despite the small sample size, we found that HLA-C alleles from group C1 HLA-C*08:02, -C*12:03, or -C*16:01 were more frequently associated in individuals with mild COVID-19 (43.8%) than in individuals with severe (8.3%; p = 0.0030; pc = 0.033) and critical (16.1%; p = 0.0014; pc = 0.0154) disease. C1 alleles are supposed to be highly efficient to present peptides to T cells, and HLA-C*12:03 may present a high number of verified epitopes from abundant SARS-CoV-2 proteins M, N, and S, thereby being allegedly able to trigger an efficient antiviral response. On the contrary, C2 alleles are usually poorly expressed on the cell surface due to low association with β2-microglobulin (β2M) and peptides, which may impede the adequate formation of stable HLA-C/β2M/peptide heterotrimers. Consequently, this pilot study described significant differences in the presence of specific HLA-C1 alleles in individuals with different clinical presentations of COVID-19, thereby suggesting that HLA haplotyping could be valuable to get further understanding in the underlying mechanisms of the impaired immune response during critical COVID-19.

Genomic features of predominant non-PCV13 serotypes responsible for adult invasive pneumococcal disease in Spain

BACKGROUND

Although pneumococcal conjugate vaccines (PCVs) effectively prevent invasive pneumococcal disease (IPD), serotype replacement has occurred.

OBJECTIVES

We studied the pangenome, antibiotic resistance mechanisms and presence of mobile elements in predominant non-PCV13 serotypes causing adult IPD after PCV13 vaccine introduction in Spain.

METHODS

We conducted a multicentre study comparing three periods in six Spanish hospitals and analysed through whole genome sequencing representative strains collected in the pre-PCV13, early-PCV13 and late-PCV13 periods.

RESULTS

Among 2197 cases of adult IPD identified, 110 pneumococci expressing non-PCV13 capsules were sequenced. Seven predominant serotypes accounted for 42.6% of IPD episodes in the late-PCV13 period: serotypes 8 (14.4%), 12F (7.5%), 9N (5.2%), 11A (4.1%), 22F (3.9%), 24F (3.9%) and 16F (3.6%). All predominant non-PCV13 serotypes were highly clonal, comprising one or two clonal complexes (CC). In general, CC538, CC4048, CC3016F, CC43322F and CC669N, related to predominant non-PCV13 serotypes, were antibiotic susceptible. CC15611A was associated with resistance to co-trimoxazole, penicillin and amoxicillin. CC23024F was non-susceptible to penicillin and resistant to erythromycin, clindamycin, and tetracycline. Six composite transposon structures of the Tn5252-family were found in CC23024F, CC98912F and CC3016F carrying different combinations of erm(B), tet(M), and cat. Pangenome analysis revealed differences in accessory genomes among the different CC, with most variety in CC3016F (23.9%) and more conservation in CC15611A (8.5%).

CONCLUSIONS

We identified highly clonal predominant serotypes responsible for IPD in adults. The detection of not only conjugative elements carrying resistance determinants but also clones previously associated with vaccine serotypes (CC15611A and CC23024F) highlights the importance of the accessory genome.

Identification of Promoter Region Markers Associated With Altered Expression of Resistance-Nodulation-Division Antibiotic Efflux Pumps in Acinetobacter baumannii

Genetic alterations leading to the constitutive upregulation of specific efflux pumps contribute to antibacterial resistance in multidrug resistant bacteria. The identification of such resistance markers remains one of the most challenging tasks of genome-level resistance predictors. In this study, 487 non-redundant genetic events were identified in upstream zones of three operons coding for resistance-nodulation-division (RND) efflux pumps of 4,130 Acinetobacter baumannii isolates. These events included insertion sequences, small indels, and single nucleotide polymorphisms. In some cases, alterations explicitly modified the expression motifs described for these operons, such as the promoter boxes, operators, and Shine-Dalgarno sequences. In addition, changes in DNA curvature and mRNA secondary structures, which are structural elements that regulate expression, were also calculated. According to their influence on RND upregulation, the catalog of upstream modifications were associated with “experimentally verified,” “presumed,” and “probably irrelevant” degrees of certainty. For experimental verification, DNA of upstream sequences independently carrying selected markers, three for each RND operon, were fused to a luciferase reporter plasmid system. Five out of the nine selected markers tested showed significant increases in expression with respect to the wild-type sequence control. In particular, a 25-fold expression increase was observed with the ISAba1 insertion sequence upstream the adeABC pump. Next, overexpression of each of the three multi-specific RND pumps was linked to their respective antibacterial substrates by a deep A. baumannii literature screen. Consequently, a data flow framework was then developed to link genomic upregulatory RND determinants to potential antibiotic resistance. Assignment of potential increases in minimal inhibitory concentrations at the “experimentally verified” level was permitted for 42 isolates to 7–8 unrelated antibacterial agents including tigecycline, which is overlooked by conventional resistome predictors. Thus, our protocol may represent a time-saving filter step prior to laborious confirmation experiments for efflux-driven resistance. Altogether, a computational-experimental pipeline containing all components required for identifying the upstream regulatory resistome is proposed. This schema may provide the foundational stone for the elaboration of tools approaching antibiotic efflux that complement routine resistome predictors for preventing antimicrobial therapy failure against difficult-to-threat bacteria.

Cross-Recognition of SARS-CoV-2 B-Cell Epitopes with Other Betacoronavirus Nucleoproteins

The B and T lymphocytes of the adaptive immune system are important for the control of most viral infections, including COVID-19. Identification of epitopes recognized by these cells is fundamental for understanding how the immune system detects and removes pathogens, and for antiviral vaccine design. Intriguingly, several cross-reactive T lymphocyte epitopes from SARS-CoV-2 with other betacoronaviruses responsible for the common cold have been identified. In addition, antibodies that cross-recognize the spike protein, but not the nucleoprotein (N protein), from different betacoronavirus have also been reported. Using a consensus of eight bioinformatic methods for predicting B-cell epitopes and the collection of experimentally detected epitopes for SARS-CoV and SARS-CoV-2, we identified four surface-exposed, conserved, and hypothetical antigenic regions that are exclusive of the N protein. These regions were analyzed using ELISA assays with two cohorts: SARS-CoV-2 infected patients and pre-COVID-19 samples. Here we describe four epitopes from SARS-CoV-2 N protein that are recognized by the humoral response from multiple individuals infected with COVID-19, and are conserved in other human coronaviruses. Three of these linear surface-exposed sequences and their peptide homologs in SARS-CoV-2 and HCoV-OC43 were also recognized by antibodies from pre-COVID-19 serum samples, indicating cross-reactivity of antibodies against coronavirus N proteins. Different conserved human coronaviruses (HCoVs) cross-reactive B epitopes against SARS-CoV-2 N protein are detected in a significant fraction of individuals not exposed to this pandemic virus. These results have potential clinical implications.

Predicted impact of the viral mutational landscape on the cytotoxic response against SARS-CoV-2

The massive assessment of immune evasion due to viral mutations that increase COVID-19 susceptibility can be computationally facilitated. The adaptive cytotoxic T response is critical during primary infection and the generation of long-term protection. Here, potential HLA class I epitopes in the SARS-CoV-2 proteome were predicted for 2,915 human alleles of 71 families using the netMHCIpan EL algorithm. Allele families showed extreme epitopic differences, underscoring genetic variability of protective capacity between humans. Up to 1,222 epitopes were associated with any of the twelve supertypes, that is, allele clusters covering 90% population. Next, from all mutations identified in ~118,000 viral NCBI isolates, those causing significant epitope score reduction were considered epitope escape mutations. These mutations mainly involved non-conservative substitutions at the second and C-terminal position of the ligand core, or total ligand removal by large recurrent deletions. Escape mutations affected 47% of supertype epitopes, which in 21% of cases concerned isolates from two or more sub-continental areas. Some of these changes were coupled, but never surpassed 15% of evaded epitopes for the same supertype in the same isolate, except for B27. In contrast to most supertypes, eight allele families mostly contained alleles with few SARS-CoV-2 ligands. Isolates harboring cytotoxic escape mutations for these families co-existed geographically within sub-Saharan and Asian populations enriched in these alleles according to the Allele Frequency Net Database. Collectively, our findings indicate that escape mutation events have already occurred for half of HLA class I supertype epitopes. However, it is presently unlikely that, overall, it poses a threat to the global population. In contrast, single and double mutations for susceptible alleles may be associated with viral selective pressure and alarming local outbreaks. The integration of genomic, geographical and immunoinformatic information eases the surveillance of variants potentially affecting the global population, as well as minority subpopulations.

Abundance, Betweenness Centrality, Hydrophobicity, and Isoelectric Points Are Relevant Factors in the Processing of Parental Proteins of the HLA Class II Ligandome

Adaptive cellular and humoral immune responses to infectious agents require previous recognition of pathogenic peptides bound to human leukocyte antigen (HLA) class II molecules exposed on the surface of the professional antigen-presenting cells. Knowledge of how these peptide ligands are generated is essential to understand the basis for CD4⁺ T-cell-mediated immunity and tolerance. In this study, a high-throughput mass spectrometry analysis was used to identify more than 16,000 cell peptides bound to several HLA-DR and -DP class II molecules isolated from large amounts of uninfected and virus-infected human cells (ProteomeXchange accession: PXD028006). The analysis of the 1808 parental proteins containing HLA class II ligands revealed that these cell proteins were more acidic, abundant, and highly connected but less hydrophilic than non-parental proteomes. Therefore, the percentage of acidic residues was increased and hydroxyl and polar residues were decreased in the parental proteins for the HLA class II ligandomes versus the non-parental proteomes. This definition of the properties shared by parental proteins that constitute the source of the HLA class II ligandomes can serve as the basis for the development of bioinformatics tools to predict proteins that are most likely recognized by the immune system through the CD4⁺ helper T lymphocytes in both autoimmunity and infection.

Predicted Epitope Abundance Supports Vaccine-Induced Cytotoxic Protection Against SARS-CoV-2 Variants of Concern

The effect of emerging SARS-CoV-2 variants on vaccine efficacy is of critical importance. In this study, the potential impact of mutations that facilitate escape from the cytotoxic cellular immune response in these new virus variants for the 551 most abundant HLA class I alleles was analyzed. Computational prediction showed that most of these alleles, that cover >90% of the population, contain enough epitopes without escape mutations in the principal SARS-CoV-2 variants. These data suggest that the cytotoxic cellular immune protection elicited by vaccination is not greatly affected by emerging SARS-CoV-2 variants.

Streptococcus pneumoniae: a Plethora of Temperate Bacteriophages With a Role in Host Genome Rearrangement

Bacteriophages (phages) are viruses that infect bacteria. They are the most abundant biological entity on Earth (current estimates suggest there to be perhaps 10³¹ particles) and are found nearly everywhere. Temperate phages can integrate into the chromosome of their host, and prophages have been found in abundance in sequenced bacterial genomes. Prophages may modulate the virulence of their host in different ways, e.g., by the secretion of phage-encoded toxins or by mediating bacterial infectivity. Some 70% of Streptococcus pneumoniae (the pneumococcus)—a frequent cause of otitis media, pneumonia, bacteremia and meningitis—isolates harbor one or more prophages. In the present study, over 4000 S. pneumoniae genomes were examined for the presence of prophages, and nearly 90% were found to contain at least one prophage, either defective (47%) or present in full (43%). More than 7000 complete putative integrases, either of the tyrosine (6243) or serine (957) families, and 1210 full-sized endolysins (among them 1180 enzymes corresponding to 318 amino acid-long N-acetylmuramoyl-L-alanine amidases [LytA_PPH]) were found. Based on their integration site, 26 different pneumococcal prophage groups were documented. Prophages coding for tRNAs, putative virulence factors and different methyltransferases were also detected. The members of one group of diverse prophages (PPH090) were found to integrate into the 3’ end of the host lytA_Spn gene encoding the major S. pneumoniae autolysin without disrupting it. The great similarity of the lytA_Spnand lytA_PPH genes (85–92% identity) allowed them to recombine, via an apparent integrase-independent mechanism, to produce different DNA rearrangements within the pneumococcal chromosome. This study provides a complete dataset that can be used to further analyze pneumococcal prophages, their evolutionary relationships, and their role in the pathogenesis of pneumococcal disease.

Designing Multi-Antigen Vaccines Against Acinetobacter baumannii Using Systemic Approaches

Vaccines and monoclonal antibodies are promising approaches for preventing and treating infections caused by multidrug resistant Acinetobacter baumannii. However, only partial protection has been achieved with many previously tested protein antigens, which suggests that vaccines incorporating multiple antigens may be necessary in order to obtain high levels of protection. Several aspects that use the wealth of omic data available for A. baumannii have not been fully exploited for antigen identification. In this study, the use of fractionated proteomic and computational data from ~4,200 genomes increased the number of proteins potentially accessible to the humoral response to 8,824 non-redundant proteins in the A. baumannii panproteome. Among them, 59% carried predicted B-cell epitopes and T-cell epitopes recognized by two or more alleles of the HLA class II DP supertype. Potential cross-reactivity with human proteins was detected for 8.9% of antigens at the protein level and 2.7% at the B-cell epitope level. Individual antigens were associated with different infection types by genomic, transcriptomic or functional analyses. High intra-clonal genome density permitted the identification of international clone II as a “vaccitype”, in which 20% of identified antigens were specific to this clone. Network-based centrality measurements were used to identify multiple immunologic nodes. Data were formatted, unified and stored in a data warehouse database, which was subsequently used to identify synergistic antigen combinations for different vaccination strategies. This study supports the idea that integration of multi-omic data and fundamental knowledge of the pathobiology of drug-resistant bacteria can facilitate the development of effective multi-antigen vaccines against these challenging infections.

Immunization with SP_1992 (DiiA) Protein of Streptococcus pneumoniae Reduces Nasopharyngeal Colonization and Protects against Invasive Disease in Mice

Knowledge-based vaccinology can reveal uncharacterized antigen candidates for a new generation of protein-based anti-pneumococcal vaccines. DiiA, encoded by the sp_1992 locus, is a surface protein containing either one or two repeats of a 37mer N-terminal motif that exhibits low interstrain variability. DiiA belongs to the core proteome, contains several conserved B-cell epitopes, and is associated with colonization and pathogenesis. Immunization with DiiA protein via the intraperitoneal route induced a strong IgG response, including different IgG subtypes. Vaccination with DiiA increased bacterial clearance and induced protection against sepsis, conferring 70% increased survival at 48 h post-infection when compared to the adjuvant control. The immunogenic response and survival rates in mice immunized with a truncated DiiA version lacking 119 N-terminal residues were remarkably lower, confirming the relevance of the repeat zone in the immunoprotection by DiiA. Intranasal immunization of mice with the entire recombinant protein elicited mucosal IgG and IgA responses that reduced bacterial colonization of the nasopharynx, confirming that this protein might be a vaccine candidate for reducing the carrier rate. DiiA constitutes an example of how functionally unannotated proteins may still represent promising candidates that can be used in prophylactic strategies against the pneumococcal carrier state and invasive disease.

Modulation of Natural HLA-B*27:05 Ligandome by Ankylosing Spondylitis-associated Endoplasmic Reticulum Aminopeptidase 2 (ERAP2)

The HLA-B*27:05 allele and the endoplasmic reticulum-resident aminopeptidases are strongly associated with AS, a chronic inflammatory spondyloarthropathy. This study examined the effect of ERAP2 in the generation of the natural HLA-B*27:05 ligandome in live cells. Complexes of HLA-B*27:05-bound peptide pools were isolated from human ERAP2-edited cell clones, and the peptides were identified using high-throughput mass spectrometry analyses. The relative abundance of a thousand ligands was established by quantitative tandem mass spectrometry and bioinformatics analysis. The residue frequencies at different peptide position, identified in the presence or absence of ERAP2, determined structural features of ligands and their interactions with specific pockets of the antigen-binding site of the HLA-B*27:05 molecule. Sequence alignment of ligands identified with species of bacteria associated with HLA-B*27-dependent reactive arthritis was performed. In the absence of ERAP2, peptides with N-terminal basic residues and minority canonical P2 residues are enriched in the natural ligandome. Further, alterations of residue frequencies and hydrophobicity profile at P3, P7, and PΩ positions were detected. In addition, several ERAP2-dependent cellular peptides were highly similar to protein sequences of arthritogenic bacteria, including one human HLA-B*27:05 ligand fully conserved in a protein from Campylobacter jejuni These findings highlight the pathogenic role of this aminopeptidase in the triggering of AS autoimmune disease.

The Versatility of Opportunistic Infections Caused by Gemella Isolates Is Supported by the Carriage of Virulence Factors From Multiple Origins

The molecular basis of the pathogenesis of the opportunistic invasive infections caused by isolates of the Gemella genus remains largely unknown. Moreover, inconsistencies in the current species assignation were detected after genome-level comparison of 16 public Gemella isolates. A literature search detected that, between the two most pathogenic species, Gemella morbillorum causes about twice the number of cases compared to Gemella haemolysans. These two species shared their mean diseases - sepsis and endocarditis - but differed in causing other syndromes. A number of well-known virulence factors were harbored by all species, such as a manganese transport/adhesin sharing 83% identity from oral endocarditis-causing streptococci. Likewise, all Gemellae carried the genes required for incorporating phosphorylcholine into their cell walls and encoded some choline-binding proteins. In contrast, other proteins were species-specific, which may justify the known epidemiological differences. G. haemolysans, but not G. morbillorum, harbor a gene cluster potentially encoding a polysaccharidic capsule. Species-specific surface determinants also included Rib and MucBP repeats, hemoglobin-binding NEAT domains, peptidases of C5a complement factor and domains that recognize extracellular matrix molecules exposed in damaged heart valves, such as collagen and fibronectin. Surface virulence determinants were associated with several taxonomically dispersed opportunistic genera of the oral microbiota, such as Granulicatella, Parvimonas, and Streptococcus, suggesting the existence of a horizontally transferrable gene reservoir in the oral environment, likely facilitated by close proximity in biofilms and ultimately linked to endocarditis. The identification of the Gemella virulence pool should be implemented in whole genome-based protocols to rationally predict the pathogenic potential in ongoing clinical infections caused by these poorly known bacterial pathogens.

Using Omics Technologies and Systems Biology to Identify Epitope Targets for the Development of Monoclonal Antibodies Against Antibiotic-Resistant Bacteria

Over the past few decades, antimicrobial resistance has emerged as an important threat to public health due to the global dissemination of multidrug-resistant strains from several bacterial species. This worrisome trend, in addition to the paucity of new antibiotics with novel mechanisms of action in the development pipeline, warrants the development of non-antimicrobial approaches to combating infection caused by these isolates. Monoclonal antibodies (mAbs) have emerged as highly effective molecules for the treatment of multiple diseases. However, in spite of the fact that antibodies play an important role in protective immunity against bacteria, only three mAb therapies have been approved for clinical use in the treatment of bacterial infections. In the present review, we briefly outline the therapeutic potential of mAbs in the treatment of bacterial diseases and discuss how their development can be facilitated when assisted by “omics” technologies and interpreted under a systems biology paradigm. Specifically, methods employing large genomic, transcriptomic, structural, and proteomic datasets allow for the rational identification of epitopes. Ideally, these include those that are present in the majority of circulating isolates, highly conserved at the amino acid level, surface-exposed, located on antigens essential for virulence, and expressed during critical stages of infection. Therefore, these knowledge-based approaches can contribute to the identification of high-value epitopes for the development of effective mAbs against challenging bacterial clones.

Evolution of the β-lactam-resistant Streptococcus pneumoniae PMEN3 clone over a 30 year period in Barcelona, Spain

Objectives

To analyse the epidemiology and genetic evolution of PMEN3 (Spain9V-156), a penicillin-non-susceptible clone of Streptococcus pneumoniae, causing invasive pneumococcal disease (IPD) in Barcelona during 1987-2016.

Methods

WGS was performed on 46 representative isolates and the data were used to design additional molecular typing methods including partial MLST, PCR-RFLP and detection of surface-exposed proteins and prophages, to assign the remaining isolates to lineages. The isolates were also subjected to antimicrobial susceptibility testing.

Results

Two hundred and twenty-seven adult cases of IPD caused by PMEN3 were identified. PMEN3 caused mainly pneumonia (84%) and the 30 day mortality rate was 23.1%. Evidence of recombination events was found, mostly in three regions, namely the capsular operon (associated with capsular switching) and adjacent regions containing pbp2x and pbp1a, the murM gene and the pbp2b-ddl region. Some of these genetic changes generated successful new variant serotype lineages, including one of serotype 11A that is not included in the current PCV13 vaccine. Other genetic changes led to increased MICs of β-lactams. Notably, most isolates also harboured prophages coding for PblB-like proteins. Despite these adaptations, the ability of this clone to cause IPD remained unchanged over time, highlighting the importance of its core genetic background.

Conclusions

Our study demonstrated successful adaptation of PMEN3 to persist over time despite the introduction of broader antibiotics and conjugate vaccines. In addition to enhancing understanding of the molecular evolution of PMEN3, these findings highlight the need for the development of non-serotype-based vaccines to fight pneumococcal infection.

An Uncharacterized Member of the Gls24 Protein Superfamily Is a Putative Sensor of Essential Amino Acid Availability in Streptococcus pneumoniae

Proteins belonging to the Gls24 superfamily are involved in survival of pathogenic Gram-positive cocci under oligotrophic conditions and other types of stress, by a still unknown molecular mechanism. In Firmicutes, this superfamily includes three different valine-rich orthologal families (Gls24A, B, C) with different potential interactive partners. Whereas the Streptococcus pneumoniae Δgls24A deletion mutant experienced a general long growth delay, the Δgls24B mutant grew as the parental strain in the semisynthetic AGCH medium but failed to grow in the complex Todd-Hewitt medium. Bovine seroalbumin (BSA) was the component responsible for this phenotype. The effect of BSA on growth was concentration-dependent and was maintained when the protein was proteolyzed but not when heat-denatured, suggesting that BSA dependence was related to oligopeptide supplementation. Global transcriptional analyses of the knockout mutant revealed catabolic derepression and induction of chaperone and oligopeptide transport genes. This mutant also showed increased sensibility to cadmium and high temperature. The Δgls24B mutant behaved as a poor colonizer in the nasopharynx of mice and showed 20-fold competence impairment. Experimental data suggest that Gls24B plays a central role as a sensor of amino acid availability and its connection to sugar catabolism. This metabolic rewiring can be compensated in vitro, at the expenses of external oligopeptide supplementation, but reduce important bacteria skills prior to efficiently address systemic virulence traits. This is an example of how metabolic factors conserved in enterococci, streptococci, and staphylococci can be essential for survival in poor oligopeptide environments prior to infection progression.

Proteomics Analysis Reveals That Structural Proteins of the Virion Core and Involved in Gene Expression Are the Main Source for HLA Class II Ligands in Vaccinia Virus-Infected Cells

Protective cellular and humoral immune responses require previous recognition of viral antigenic peptides complexed with human leukocyte antigen (HLA) class II molecules on the surface of the antigen presenting cells. The HLA class II-restricted immune response is important for the control and the clearance of poxvirus infection including vaccinia virus (VACV), the vaccine used in the worldwide eradication of smallpox. In this study, a mass spectrometry analysis was used to identify VACV ligands bound to HLA-DR and -DP class II molecules present on the surface of VACV-infected cells. Twenty-six naturally processed viral ligands among the tens of thousands of cell peptides bound to HLA class II proteins were identified. These viral ligands arose from 19 parental VACV proteins: A4, A5, A18, A35, A38, B5, B13, D1, D5, D7, D12, D13, E3, E8, H5, I2, I3, J2, and K2. The majority of these VACV proteins yielded one HLA ligand and were generated mainly, but not exclusively, by the classical HLA class II antigen processing pathway. Medium-sized and abundant proteins from the virion core and/or involved in the viral gene expression were the major source of VACV ligands bound to HLA-DR and -DP class II molecules. These findings will help to understand the effectiveness of current poxvirus-based vaccines and will be important in the design of new ones.

Bridging Chromosomal Architecture and Pathophysiology of Streptococcus pneumoniae

The chromosome of Streptococcus pneumoniae is organized into topological domains based on its transcriptional response to DNA relaxation: Up-regulated (UP), down-regulated (DOWN), nonregulated (NR), and AT-rich. In the present work, NR genes found to have highly conserved chromosomal locations (17% of the genome) were categorized as members of position-conserved nonregulated (pcNR) domains, while NR genes with a variable position (36% of the genome) were classified as members of position-variable nonregulated (pvNR) domains. On average, pcNR domains showed high transcription rates, optimized codon usage, and were found to contain only a small number of RUP/BOX/SPLICE repeats. They were also poor in exogenous genes but enriched in leading strand genes that code for proteins involved in primary metabolism with central roles within the interactome. In contrast, pvNR genes coding for cell wall proteins, paralogs, virulence factors and immunogenic candidates for protein-based vaccines were found to be overrepresented. DOWN domains were enriched in genes essential for infection. Many UP and DOWN domain genes were seen to be activated during different stages of competence, whereas pcNR genes tended to be repressed until the competence was switched off. Pneumococcal genes appear to be subject to a topology-driven selection pressure that defines the chromosomal location of genes involved in metabolism, virulence and competence. The pcNR domains are interleaved between UP and DOWN domains according to a pattern that suggests the existence of macrodomain entities. The term “topogenomics” is here proposed to describe the study of the topological rules of genomes and their relationship with physiology.

Upregulation of the PatAB Transporter Confers Fluoroquinolone Resistance to Streptococcus pseudopneumoniae

We characterized the mechanism of fluoroquinolone-resistance in two isolates of Streptococcus pseudopneumoniae having fluoroquinolone-efflux as unique mechanism of resistance. Whole genome sequencing and genetic transformation experiments were performed together with phenotypic determinations of the efflux mechanism. The PatAB pump was identified as responsible for efflux of ciprofloxacin (MIC of 4 μg/ml), ethidium bromide (MICs of 8-16 μg/ml) and acriflavine (MICs of 4-8 μg/ml) in both isolates. These MICs were at least 8-fold lower in the presence of the efflux inhibitor reserpine. Complete genome sequencing indicated that the sequence located between the promoter of the patAB operon and the initiation codon of patA, which putatively forms an RNA stem-loop structure, may be responsible for the efflux phenotype. RT-qPCR determinations performed on RNAs of cultures treated or not treated with subinhibitory ciprofloxacin concentrations were performed. While no significant changes were observed in wild-type Streptococcus pneumoniae R6 strain, increases in transcription were detected in the ciprofloxacin-efflux transformants obtained with DNA from efflux-positive isolates, in the ranges of 1.4 to 3.4-fold (patA) and 2.1 to 2.9-fold (patB). Ciprofloxacin-induction was related with a lower predicted free energy for the stem-loop structure in the RNA of S. pseudopneumoniae isolates (-13.81 and -8.58) than for R6 (-15.32 kcal/mol), which may ease transcription. The presence of these regulatory variations in commensal S. pseudopneumoniae isolates, and the possibility of its transfer to Streptococcus pneumoniae by genetic transformation, could increase fluoroquinolone resistance in this important pathogen.

The Transcriptome of Streptococcus pneumoniae Induced by Local and Global Changes in Supercoiling

The bacterial chromosome is compacted in a manner optimal for DNA transactions to occur. The degree of compaction results from the level of DNA-supercoiling and the presence of nucleoid-binding proteins. DNA-supercoiling is homeostatically maintained by the opposing activities of relaxing DNA topoisomerases and negative supercoil-inducing DNA gyrase. DNA-supercoiling acts as a general cis regulator of transcription, which can be superimposed upon other types of more specific trans regulatory mechanism. Transcriptomic studies on the human pathogen Streptococcus pneumoniae, which has a relatively small genome (∼2 Mb) and few nucleoid-binding proteins, have been performed under conditions of local and global changes in supercoiling. The response to local changes induced by fluoroquinolone antibiotics, which target DNA gyrase subunit A and/or topoisomerase IV, involves an increase in oxygen radicals which reduces cell viability, while the induction of global supercoiling changes by novobiocin (a DNA gyrase subunit B inhibitor), or by seconeolitsine (a topoisomerase I inhibitor), has revealed the existence of topological domains that specifically respond to such changes. The control of DNA-supercoiling in S. pneumoniae occurs mainly via the regulation of topoisomerase gene transcription: relaxation triggers the up-regulation of gyrase and the down-regulation of topoisomerases I and IV, while hypernegative supercoiling down-regulates the expression of topoisomerase I. Relaxation affects 13% of the genome, with the majority of the genes affected located in 15 domains. Hypernegative supercoiling affects 10% of the genome, with one quarter of the genes affected located in 12 domains. However, all the above domains overlap, suggesting that the chromosome is organized into topological domains with fixed locations. Based on its response to relaxation, the pneumococcal chromosome can be said to be organized into five types of domain: up-regulated, down-regulated, position-conserved non-regulated, position-variable non-regulated, and AT-rich. The AT content is higher in the up-regulated than in the down-regulated domains. Genes within the different domains share structural and functional characteristics. It would seem that a topology-driven selection pressure has defined the chromosomal location of the metabolism, virulence and competence genes, which suggests the existence of topological rules that aim to improve bacterial fitness.

The MiiA motif is a common marker present in polytopic surface proteins of oral and urinary tract invasive bacteria

Many surface virulence factors of bacterial pathogens show mosaicism and confounding phylogenetic origin. The Streptococcus gordonii platelet-binding GspB protein, the Streptococcus sanguinis SrpA adhesin and the Streptococcus pneumoniae DiiA protein, share an imperfect 27-residue motif. Given the disparate domain architectures of these proteins and its association to invasive disease, this motif was named MiiA from Multiarchitecture invasion-involved motif A. MiiA is predicted to adopt a beta-sheet folding, probably related to the Ig-like fold, with a symmetrical positioning of two conserved aspartic residues. A specific hidden Markov model profiling MiiA was built, which specifically detected the motif in proteins from 58 species, mainly in cell-wall proteins from Gram-positive bacteria. These proteins contained one to ten MiiA motifs, which were embedded within larger repeat units of 70-82 residues. MiiA motifs combined to other domains and elements such as coiled-coils and low-complexity regions. The species carrying MiiA-proteins included commensals from the urogenital tract and the oral cavity, which can cause opportunistic endocarditis and sepsis. Intra-protein MiiA repeats showed a complex mixture of orthologal, paralogal and inter-species relationships, suggestive of a multistep origin. Presence of these repeats in proteins involved in oligosaccharide recognition and lifestyle of species suggest a putative function for MiiA repeats in sugars binding, probably those present in receptors of epithelial and blood cells. MiiA modules appear to have been transferred horizontally between species co-habiting in the same niche to create their own MiiA-containing determinants. The present work provides a global study and a catalog of potential MiiA virulence factors that should be analyzed experimentally.

An increase in negative supercoiling in bacteria reveals topology-reacting gene clusters and a homeostatic response mediated by the DNA topoisomerase I gene

We studied the transcriptional response to an increase in DNA supercoiling in Streptococcus pneumoniae by using seconeolitsine, a new topoisomerase I inhibitor. A homeostatic response allowing recovery of supercoiling was observed in cells treated with subinhibitory seconeolitsine concentrations. Supercoiling increases of 40.7% (6 μM) and 72.9% (8 μM) were lowered to 8.5% and 44.1%, respectively. Likewise, drug removal facilitated the recovery of cell viability and DNA-supercoiling. Transcription of topoisomerase I depended on the supercoiling level. Also specific binding of topoisomerase I to the gyrase A gene promoter was detected by chromatin-immunoprecipitation. The transcriptomic response to 8 μM seconeolitsine had two stages. An early stage, associated to an increase in supercoiling, affected 10% of the genome. A late stage, manifested by supercoiling recovery, affected 2% of the genome. Nearly 25% of the early responsive genes formed 12 clusters with a coordinated transcription. Clusters were 6.7–31.4 kb in length and included 9–22 responsive genes. These clusters partially overlapped with those observed under DNA relaxation, suggesting that bacteria manage supercoiling stress using pathways with common components. This is the first report of a coordinated global transcriptomic response that is triggered by an increase in DNA supercoiling in bacteria.

Insights into the Evolutionary Relationships of LytA Autolysin and Ply Pneumolysin-Like Genes in Streptococcus pneumoniae and Related Streptococci

Streptococcus pneumoniae (pneumococcus) is a major human pathogen. The main pneumococcal autolysin LytA and the pneumolysin Ply are two of the bacterium's most important virulence factors. The lytA- and ply-related genes are also found in other streptococci of the Mitis group (SMG). The precise characteristics of the lytA-related-but not the ply-related-genes of SMG and their prophages have been previously described. A search of the more than 400 SMG genomic sequences available in public databases (ca. 300 for S. pneumoniae), showed Streptococcus pseudopneumoniae IS7493 to harbor four ply-related genes, two of which (plyA and plyB) have 98% identical nucleotides. The plyA homolog of S. pseudopneumoniae is conserved in all S. pneumoniae strains, and seems to be included in a pathogenicity island together with the lytA gene. However, only nonencapsulated S. pneumoniae strains possess a plyB gene, which is part of an integrative and conjugative element. Notably, the existence of a bacterial lytA-related gene in a genome is linked to the presence of plyA and vice versa. The present analysis also shows there are eight main types of plyA-lytA genomic islands. A possible stepwise scenario for the evolution of the plyA-lytA island in S. pneumoniae is proposed.

Inspecting the potential physiological and biomedical value of 44 conserved uncharacterised proteins of Streptococcus pneumoniae

Background

The major Gram-positive coccoid pathogens cause similar invasive diseases and show high rates of antimicrobial resistance. Uncharacterised proteins shared by these organisms may be involved in virulence or be targets for antimicrobial therapy.

Results

Forty four uncharacterised proteins from Streptococcus pneumoniae with homologues in Enterococcus faecalis and/or Staphylococcus aureus were selected for analysis. These proteins showed differences in terms of sequence conservation and number of interacting partners. Twenty eight of these proteins were monodomain proteins and 16 were modular, involving domain combinations and, in many cases, predicted unstructured regions. The genes coding for four of these 44 proteins were essential. Genomic and structural studies showed one of the four essential genes to code for a promising antibacterial target. The strongest impact of gene removal was on monodomain proteins showing high sequence conservation and/or interactions with many other proteins. Eleven out of 40 knockouts (one for each gene) showed growth delay and 10 knockouts presented a chaining phenotype. Five of these chaining mutants showed a lack of putative DNA-binding proteins. This suggest this phenotype results from a loss of overall transcription regulation. Five knockouts showed defective autolysis in response to penicillin and vancomycin, and attenuated virulence in an animal model of sepsis.

Conclusions

Uncharacterised proteins make up a reservoir of polypeptides of different physiological importance and biomedical potential. A promising antibacterial target was identified. Five of the 44 examined proteins seemed to be virulence factors.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-652) contains supplementary material, which is available to authorized users.

Role of global and local topology in the regulation of gene expression in Streptococcus pneumoniae

The most basic level of transcription regulation in Streptococcus pneumoniae is the organization of its chromosome in topological domains. In response to drugs that caused DNA-relaxation, a global transcriptional response was observed. Several chromosomal domains were identified based on the transcriptional response of their genes: up-regulated (U), down-regulated (D), non-regulated (N), and flanking (F). We show that these distinct domains have different expression and conservation characteristics. Microarray fluorescence units under non-relaxation conditions were used as a measure of gene transcriptional level. Fluorescence units were significantly lower in F genes than in the other domains with a similar AT content. The transcriptional level of the domains categorized them was D>U>F. In addition, a comparison of 12 S. pneumoniae genome sequences showed a conservation of gene composition within U and D domains, and an extensive gene interchange in F domains. We tested the organization of chromosomal domains by measuring the relaxation-mediated transcription of eight insertions of a heterologous P_tccat cassette, two in each type of domain, showing that transcription depended on their chromosomal location. Moreover, transcription from the four promoters directing the five genes involved in supercoiling homeostasis, located either in U (gyrB), D (topA), or N (gyrA and parEC) domains was analyzed both in their chromosomal locations and in a replicating plasmid. Although expression from the chromosomal P_gyrB and P_topA showed the expected domain regulation, their expression was down-regulated in the plasmid, which behaved as a D domain. However, both P_parE and P_gyrA carried their own regulatory signals, their topology-dependent expression being equivalent in the plasmid or in the chromosome. In P_gyrA a DNA bend acted as a DNA supercoiling sensor. These results revealed that DNA topology functions as a general transcriptional regulator, superimposed upon other more specific regulatory mechanisms.

Mapping flexibility and the assembly switch of cell division protein FtsZ by computational and mutational approaches

The molecular switch for nucleotide-regulated assembly and disassembly of the main prokaryotic cell division protein FtsZ is unknown despite the numerous crystal structures that are available. We have characterized the functional motions in FtsZ with a computational consensus of essential dynamics, structural comparisons, sequence conservation, and networks of co-evolving residues. Employing this information, we have constructed 17 mutants, which alter the FtsZ functional cycle at different stages, to modify FtsZ flexibility. The mutant phenotypes ranged from benign to total inactivation and included increased GTPase, reduced assembly, and stabilized assembly. Six mutations clustering at the long cleft between the C-terminal beta-sheet and core helix H7 deviated FtsZ assembly into curved filaments with inhibited GTPase, which still polymerize cooperatively. These mutations may perturb the predicted closure of the C-terminal domain onto H7 required for switching between curved and straight association modes and for GTPase activation. By mapping the FtsZ assembly switch, this work also gives insight into FtsZ druggability because the curved mutations delineate the putative binding site of the promising antibacterial FtsZ inhibitor PC190723.

The antibacterial cell division inhibitor PC190723 is an FtsZ polymer-stabilizing agent that induces filament assembly and condensation

Cell division protein FtsZ can form single-stranded filaments with a cooperative behavior by self-switching assembly. Subsequent condensation and bending of FtsZ filaments are important for the formation and constriction of the cytokinetic ring. PC190723 is an effective bactericidal cell division inhibitor that targets FtsZ in the pathogen Staphylococcus aureus and Bacillus subtilis and does not affect Escherichia coli cells, which apparently binds to a zone equivalent to the binding site of the antitumor drug taxol in tubulin (Haydon, D. J., Stokes, N. R., Ure, R., Galbraith, G., Bennett, J. M., Brown, D. R., Baker, P. J., Barynin, V. V., Rice, D. W., Sedelnikova, S. E., Heal, J. R., Sheridan, J. M., Aiwale, S. T., Chauhan, P. K., Srivastava, A., Taneja, A., Collins, I., Errington, J., and Czaplewski, L. G. (2008) Science 312, 1673-1675). We have found that the benzamide derivative PC190723 is an FtsZ polymer-stabilizing agent. PC190723 induced nucleated assembly of Bs-FtsZ into single-stranded coiled protofilaments and polymorphic condensates, including bundles, coils, and toroids, whose formation could be modulated with different solution conditions. Under conditions for reversible assembly of Bs-FtsZ, PC190723 binding reduced the GTPase activity and induced the formation of straight bundles and ribbons, which was also observed with Sa-FtsZ but not with nonsusceptible Ec-FtsZ. The fragment 2,6-difluoro-3-methoxybenzamide also induced Bs-FtsZ bundling. We propose that polymer stabilization by PC190723 suppresses in vivo FtsZ polymer dynamics and bacterial division. The biochemical action of PC190723 on FtsZ parallels that of the microtubule-stabilizing agent taxol on the eukaryotic structural homologue tubulin. Both taxol and PC190723 stabilize polymers against disassembly by preferential binding to each assembled protein. It is yet to be investigated whether both ligands target structurally related assembly switches.

The genome of Streptococcus pneumoniae is organized in topology-reacting gene clusters

The transcriptional response of Streptococcus pneumoniae was examined after exposure to the GyrB-inhibitor novobiocin. Topoisomer distributions of an internal plasmid confirmed DNA relaxation and recovery of the native level of supercoiling at low novobiocin concentrations. This was due to the up-regulation of DNA gyrase and the down-regulation of topoisomerases I and IV. In addition, >13% of the genome exhibited relaxation-dependent transcription. The majority of the responsive genes (>68%) fell into 15 physical clusters (14.6–85.6 kb) that underwent coordinated regulation, independently of operon organization. These genomic clusters correlated with AT content and codon composition, showing the chromosome to be organized into topology-reacting gene clusters that respond to DNA supercoiling. In particular, down-regulated clusters were flanked by 11–40 kb AT-rich zones that might have a putative structural function. This is the first case where genes responding to changes in the level of supercoiling in a coordinated manner were found organized as functional clusters. Such an organization revealed DNA supercoiling as a general feature that controls gene expression superimposed on other kinds of more specific regulatory mechanisms.

Molecular characterization of disease-associated streptococci of the mitis group that are optochin susceptible

Eight optochin-susceptible (Opt(s)) alpha-hemolytic (viridans) streptococcus isolates were characterized at the molecular level. These isolates showed phenotypic characteristics typical of both viridans streptococci and Streptococcus pneumoniae. Comparison of the sequence of housekeeping genes from these isolates with those of S. pneumoniae, Streptococcus mitis, Streptococcus oralis, and Streptococcus pseudopneumoniae suggested that the Opt(s) isolates corresponded to streptococci of the mitis group. Besides, the Opt(s) streptococci were negative by a Gen-Probe AccuProbe pneumococcus test and hybridized with specific pneumococcal probes (lytA and ply) but also with ant, a gene not present in most S. pneumoniae strains. Moreover, the isolates were insoluble in 1% sodium deoxycholate but completely dissolved in 0.1% deoxycholate. Sequence analysis of the lytA gene revealed that the Opt(s) streptococci carried lytA alleles characteristic of those present in nonpneumococcal streptococci of the mitis group. The determination of the partial nucleotide sequence embracing the atp operon encoding the F(o)F(1) H(+)-ATPase indicated that the optochin susceptibility of the isolates was due to the acquisition of atpC, atpA, and part of atpB from S. pneumoniae by horizontal gene transfer.

Transcriptional analysis of the acid tolerance response in Streptococcus pneumoniae

Streptococcus pneumoniae, one of the major causes of morbidity and mortality in humans, faces a range of potentially acidic conditions in the middle and late stages of growthin vitro, in diverse human fluids during the infection process, and in biofilms present in the nasopharynx of carriers.S. pneumoniaewas shown to develop a weak acid tolerance response (ATR), where cells previously exposed to sublethal pHs (5·8–6·6) showed an increased survival rate of up to one order of magnitude after challenge at the lethal pH (4·4, survival rate of 10−4). Moreover, the survival after challenge of stationary phase cells at pH 4·4 was three orders of magnitude higher than that of cells taken from the exponential phase, due to the production of lactic acid during growth and increasing acidification of the growth medium until stationary phase. Global expression analysis after short-term (5, 15 and 30 min, the adaptation phase) and long-term (the maintenance phase) acidic shock (pH 6·0) was performed by microarray experiments, and the results were validated by real-time RT-PCR. Out of a total of 126 genes responding to acidification, 59 and 37 were specific to the adaptation phase and maintenance phase, respectively, and 30 were common to both periods. In the adaptation phase, both up- and down-regulation of gene transcripts was observed (38 and 21 genes, respectively), whereas in the maintenance phase most of the affected genes were down-regulated (34 out of 37). Genes involved in protein fate (including those involved in the protection of the protein native structure) and transport (including transporters of manganese and iron) were overrepresented among the genes affected by acidification, 8·7 and 24·6 % of the acid-responsive genes compared to 2·8 % and 9·6 % of the genome complement, respectively. Cross-regulation with the response to oxidative and osmotic stress was observed. Potential regulatory motifs involved in the ATR were identified in the promoter regions of some of the regulated genes.

Relationship between codon biased genes, microarray expression values and physiological characteristics of Streptococcus pneumoniae

A codon-profile strategy was used to predict gene expression levels in Streptococcus pneumoniae. Predicted highly expressed (PHE) genes included those encoding glycolytic and fermentative enzymes, sugar-conversion systems and carbohydrate-transporters. Additionally, some genes required for infection that are involved in oxidative metabolism and hydrogen peroxide production were PHE. Low expression values were predicted for genes encoding specific regulatory proteins like two-component systems and competence genes. Correspondence analysis localized 484 ORFs which shared a distinctive codon profile in the right horn. These genes had a mean G+C content (33·4 %) that was lower than the bulk of the genome coding sequences (39·7 %), suggesting that many of them were acquired by horizontal transfer. Half of these genes (242) were pseudogenes, ORFs shorter than 80 codons or without assigned function. The remaining genes included several virulence factors, such as capsular genes, iga, lytB, nanB, pspA, choline-binding proteins, and functions related to DNA acquisition, such as restriction-modification systems and comDE. In order to compare predicted translation rate with the relative amounts of mRNA for each gene, the codon adaptation index (CAI) values were compared with microarray fluorescence intensity values following hybridization of labelled RNA from laboratory-grown cultures. High mRNA amounts were observed in 32·5 % of PHE genes and in 64 % of the 25 genes with the highest CAI values. However, high relative amounts of RNA were also detected in 10·4 % of non-PHE genes, such as those encoding fatty acid metabolism enzymes and proteases, suggesting that their expression might also be regulated at the level of transcription or mRNA stability under the conditions tested. The effects of codon bias and mRNA amount on different gene groups in S. pneumoniae are discussed.

Genetic characterization of optochin-susceptible viridans group streptococci

Two clinical isolates of viridans group streptococci (VS) with different degrees of susceptibility to optochin (OPT), i.e., fully OPT-susceptible (Opt(s)) VS strain 1162/99 (for which the MIC was equal to that for Streptococcus pneumoniae, 0.75 micro g/ml) and intermediate Opt(s) VS strain 1174/97 (MIC, 6 micro g/ml) were studied. Besides being OPT susceptible, they showed characteristics typical of VS, such as bile insolubility; lack of reaction with pneumococcal capsular antibodies; and lack of hybridization with rRNA (AccuProbe)-, lytA-, and pnl-specific pneumococcal probes. However, these VS Opt(s) strains and VS type strains hybridized with ant, a gene not present in S. pneumoniae. A detailed characterization of the genes encoding the 16S rRNA and SodA classified isolates 1162/99 and 1174/97 as Streptococcus mitis. Analysis of the atpCAB region, which encodes the c, a, and b subunits of the F(0)F(1) H(+)-ATPase, the target of optochin, revealed high degrees of similarity between S. mitis 1162/99 and S. pneumoniae in atpC, atpA, and the N terminus of atpB. Moreover, amino acid identity between S. mitis 1174/97 and S. pneumoniae was found in alpha helix 5 of the a subunit. The organization of the chromosomal region containing the atp operon of the two Opt(s) VS and VS type strains was spr1284-atpC, with spr1284 being located 296 to 556 bp from atpC, whereas in S. pneumoniae this distance was longer than 68 kb. In addition, the gene order in S. pneumoniae was IS1239-74 bp-atpC. The results suggest that the full OPT susceptibility of S. mitis 1162/99 is due to the acquisition of atpC, atpA, and part of atpB from S. pneumoniae and that the intermediate OPT susceptibility of S. mitis 1174/97 correlates with the amino acid composition of its a subunit.

High-efficiency generation of antibiotic-resistant strains of Streptococcus pneumoniae by PCR and transformation

We designed a method by which to generate antibiotic-resistant strains of Streptococcus pneumoniae at frequencies 4 orders of magnitude greater than the spontaneous mutation rate. The method is based on the natural ability of this organism to be genetically transformed with PCR products carrying sequences homologous to its chromosome. The genes encoding the targets of ciprofloxacin (parC, encoding the ParC subunit of DNA topoisomerase IV), rifampin (rpoB, encoding the beta subunit of RNA polymerase), and streptomycin (rpsL, encoding the S12 ribosomal protein) from susceptible laboratory strain R6 were amplified by PCR and used to transform the same strain. Resistant mutants were obtained with a frequency of 10(-4) to 10(-5), depending on the fidelity of the DNA polymerase used for PCR amplifications. Ciprofloxacin-resistant mutants, for which the MICs were four-to eightfold higher than that for R6, carried a single mutation of a residue in the quinolone resistance-determining region: S79 (change to A, F, or Y) or D83 (change to N or V). Rifampin-resistant strains, for which the MICs were at least 133-fold higher than that for R6, contained a single mutation within cluster I of rpoB: S482 (change to P), Q486 (change to L), D489 (change to V), or H499 (change to L or Y). Streptomycin-resistant mutants, for which the MICs were at least 64-fold higher than that for R6, carried a mutation at either K56 (change to I, R, or T) or K101 (change to E). PCR products obtained from the mutants were able to transform R6 to resistance with high efficiency (>10(4)). This method could be used to efficiently obtain resistant mutants for any drug whose target is known.

The promoter of the operon encoding the F0F1 ATPase of Streptococcus pneumoniae is inducible by pH

The genes encoding the subunits of the F0F1 membrane ATPase of Streptococcus pneumoniae were cloned and sequenced. The eight genes, transcribed to one mRNA, are organized in an operon encoding the c, a, b, delta, alpha, gamma, beta and epsilon subunits of 66, 238, 165, 178, 501, 292, 471 and 139 amino acid residues, respectively, that were expressed in an Escherichia coli system. To investigate the role of the ATPase in the regulation of the intracellular pH, the expression of the operon between pH 5.7 and 7.5 was studied. An increase in both the ATPase activity and the amount of the alpha and beta F1 subunits as shown by Western blot analysis was observed as the pH decreased. These increases were accompanied by an increase in the atp-specific mRNA, as shown by Northern blot and slot-blot analysis. Primer extension experiments and transcriptional fusions between the atp promoter and the reporter cat gene demonstrated that this pH-dependent increase in the mRNA was regulated at the level of initiation of transcription. Transcription of the operon occurs from a promoter with a consensus -35 box (TTGACA) and a -10 box (TACACT) that differs from the consensus (TATAAT). A point mutation at the -10 box of the promoter (change to TGCACT) avoided this increase, suggesting a role for this sequence in the pH-inducible regulation.