Conservation of HLA Spike Protein Epitopes Supports T Cell Cross-Protection in SARS-CoV-2 Vaccinated Individuals against the Potentially Zoonotic Coronavirus Khosta-2.

Heterologous vaccines, which induce immunity against several related pathogens, can be a very useful and rapid way to deal with new pandemics. In this study, the potential impact of licensed COVID-19 vaccines on cytotoxic and helper cell immune responses against Khosta-2, a novel sarbecovirus that productively infects human cells, was analyzed for the 567 and 41 most common HLA class I and II alleles, respectively. Computational predictions indicated that most of these 608 alleles, covering more than 90% of the human population, contain sufficient fully conserved T-cell epitopes between the Khosta-2 and SARS-CoV-2 spike-in proteins. Ninety percent of these fully conserved peptides for class I and 93% for class II HLA molecules were verified as epitopes recognized by CD8+ or CD4+ T lymphocytes, respectively. These results show a very high correlation between bioinformatic prediction and experimental assays, which strongly validates this study. This immunoinformatics analysis allowed a broader assessment of the alleles that recognize these peptides, a global approach at the population level that is not possible with experimental assays. In summary, these findings suggest that both cytotoxic and helper cell immune protection elicited by currently licensed COVID-19 vaccines should be effective against Khosta-2 virus infection. Finally, by being rapidly adaptable to future coronavirus pandemics, this study has potential public health 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.

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 ZnSO4-inducible PZn promoter (strain ΔstaRPZnstaR) or of the maltose-inducible PMal 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 ΔstaRPZnstaR 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.

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.

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.

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.

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.

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.

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.

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 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.

Tubulin homolog TubZ in a phage-encoded partition system.

Partition systems are responsible for the process whereby large and essential plasmids are accurately positioned to daughter cells during bacterial division. They are typically made of three components: a centromere-like DNA zone, an adaptor protein, and an assembling protein that is either a Walker-box ATPase (type I) or an actin-like ATPase (type II). A recently described type III segregation system has a tubulin/FtsZ-like protein, called TubZ, for plasmid movement. Here, we present the 2.3 Å structure and dynamic assembly of a TubZ tubulin homolog from a bacteriophage and unravel the Clostridium botulinum phage c-st type III partition system. Using biochemical and biophysical approaches, we prove that a gene upstream from tubZ encodes the partner TubR and localize the centromeric region (tubS), both of which are essential for anchoring phage DNA to the motile TubZ filaments. Finally, we describe a conserved fourth component, TubY, which modulates the TubZ-R-S complex interaction.

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.

Identification of 88 regulatory small RNAs in the TIGR4 strain of the human pathogen Streptococcus pneumoniae.

Streptococcus pneumoniae is the main etiological agent of community-acquired pneumonia and a major cause of mortality and morbidity among children and the elderly. Genome sequencing of several pneumococcal strains revealed valuable information about the potential proteins and genetic diversity of this prevalent human pathogen. However, little is known about its transcriptional regulation and its small regulatory noncoding RNAs. In this study, we performed deep sequencing of the S. pneumoniae TIGR4 strain RNome to identify small regulatory RNA candidates expressed in this pathogen. We discovered 1047 potential small RNAs including intragenic, 5'- and/or 3'-overlapping RNAs and 88 small RNAs encoded in intergenic regions. With this approach, we recovered many of the previously identified intergenic small RNAs and identified 68 novel candidates, most of which are conserved in both sequence and genomic context in other S. pneumoniae strains. We confirmed the independent expression of 17 intergenic small RNAs and predicted putative mRNA targets for six of them using bioinformatics tools. Preliminary results suggest that one of these six is a key player in the regulation of competence development. This study is the biggest catalog of small noncoding RNAs reported to date in S. pneumoniae and provides a highly complete view of the small RNA network in this pathogen.

Evolution of invasive pneumococcal disease by serotype 3 in adults: a Spanish three-decade retrospective study.

Background: Invasive pneumococcal disease due to serotype 3 (S3-IPD) is associated with high mortality rates and long-term adverse effects. The introduction of the 13-valent pneumococcal conjugate vaccine (PCV13) into the Spanish paediatric immunisation programme has not led to a decrease in the adult S3-IPD. We aimed to analyse the incidence, clinical characteristics and genomics of S3-IPD in adults in Spain. Methods: Adult IPD episodes hospitalized in a Southern Barcelona hospital were prospectively collected (1994-2020). For genomic comparison, S3-IPD isolates from six Spanish hospitals (2008-2020) and historical isolates (1989-1993) were analysed by WGS (Illumina and/or MinION). Findings: From 1994 to 2020, 270 S3-IPD episodes were detected. When comparing pre-PCV (1994-2001) and late-PCV13 (2016-2020) periods, only modest changes in S3-IPD were observed (from 1.58 to 1.28 episodes per 100,000 inhabitants year). In this period, the incidence of the two main lineages shifted from 0.38 to 0.67 (CC180-GPSC12) and from 1.18 to 0.55 (CC260-GPSC83). The overall 30-day mortality remained high (24.1%), though a decrease was observed between the pre-PCV (32.4%; 95.0% CI, 22.0-45.0) and the late-PCV13 period (16.7%; 95.0% CI, 7.5-32.0) (p = 0.06). At the same time, comorbidities increased from 77.3% (95.0% CI, 65.0-86.0) to 85.7% (95.0% CI, 71.0-94.0) (p = 0.69). There were no differences in clinical characteristics or 30-day mortality between the two S3 lineages. Although both lineages were genetically homogeneous, the CC180-GPSC12 lineage presented a higher SNP density, a more open pan-genome, and a major presence of prophages and mobile genetic elements carrying resistance genes. Interpretation: Adult S3-IPD remained stable in our area over the study period despite PCV13 introduction in children. However, a clonal shift was observed. The decrease in mortality rates and the increase in comorbidities suggest a change in clinical management and overall population characteristics. The low genetic variability and absence of clinical differences between lineages highlight the role of the S3 capsule in the disease severity. Funding: This study has been funded by Instituto de Salud Carlos III (ISCIII) "PI18/00339", "PI21/01000", "INT22/00096", "FI22/00279", CIBER "CIBERES-CB06/06/0037", "CIBERINFEC-CB21/13/00009" and MSD grant "IISP 60168".

Cytomegalovirus UL44 protein induces a potent T-cell immune response in mice.

Due to the severity of CMV infection in immunocompromised individuals the development of a vaccine has been declared a priority. However, despite the efforts made there is no yet a vaccine available for clinical use. We designed an approach to identify new CMV antigens able to inducing a broad immune response that could be used in future vaccine formulations. We have used serum samples from 28 kidney transplant recipients, with a previously acquired CMV-specific immune response to identify viral proteins that were recognized by the antibodies present in the patient serum samples by Western blot. A band of approximately 45 kDa, identified as UL44, was detected by most serum samples. UL44 immunogenicity was tested in BALB/c mice that received three doses of the UL44-pcDNA DNA vaccine. UL44 elicited both, a strong antibody response and CMV-specific cellular response. Using bioinformatic analysis we demonstrated that UL44 is a highly conserved protein and contains epitopes that are able to activate CD8 lymphocytes of the most common HLA alleles in the world population. We constructed a UL44 ORF deletion mutant virus that produced no viral progeny, suggesting that UL44 is an essential viral protein. In addition, other authors have demonstrated that UL44 is one of the most abundant viral proteins after infection and have suggested an essential role of UL44 in viral replication. Altogether, our data suggests that UL44 is a potent antigen, and favored by its abundance, it may be a good candidate to include in a vaccine formulation.

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.

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.

Bacterial tubulin distinct loop sequences and primitive assembly properties support its origin from a eukaryotic tubulin ancestor.

The structure of the unique bacterial tubulin BtubA/B from Prosthecobacter is very similar to eukaryotic αß-tubulin but, strikingly, BtubA/B fold without eukaryotic chaperones. Our sequence comparisons indicate that BtubA and BtubB do not really correspond to either α- or ß-tubulin but have mosaic sequences with intertwining features from both. Their nucleotide-binding loops are more conserved, and their more divergent sequences correspond to discrete surface zones of tubulin involved in microtubule assembly and binding to eukaryotic cytosolic chaperonin, which is absent from the Prosthecobacter dejongeii draft genome. BtubA/B cooperatively assembles over a wider range of conditions than αß-tubulin, forming pairs of protofilaments that coalesce into bundles instead of microtubules, and it lacks the ability to differentially interact with divalent cations and bind typical tubulin drugs. Assembled BtubA/B contain close to one bound GTP and GDP. Both BtubA and BtubB subunits hydrolyze GTP, leading to disassembly. The mutant BtubA/B-S144G in the tubulin signature motif GGG(T/S)G(S/T)G has strongly inhibited GTPase, but BtubA-T147G/B does not, suggesting that BtubB is a more active GTPase, like ß-tubulin. BtubA/B chimera bearing the ß-tubulin loops M, H1-S2, and S9-S10 in BtubB fold, assemble, and have reduced GTPase activity. However, introduction of the α-tubulin loop S9-S10 with its unique eight-residue insertion impaired folding. From the sequence analyses, its primitive assembly features, and the properties of the chimeras, we propose that BtubA/B were acquired shortly after duplication of a spontaneously folding α- and ß-tubulin ancestor, possibly by horizontal gene transfer from a primitive eukaryotic cell, followed by divergent evolution.

Camps 2.0: exploring the sequence and structure space of prokaryotic, eukaryotic, and viral membrane proteins.

Structural bioinformatics of membrane proteins is still in its infancy, and the picture of their fold space is only beginning to emerge. Because only a handful of three-dimensional structures are available, sequence comparison and structure prediction remain the main tools for investigating sequence-structure relationships in membrane protein families. Here we present a comprehensive analysis of the structural families corresponding to α-helical membrane proteins with at least three transmembrane helices. The new version of our CAMPS database (CAMPS 2.0) covers nearly 1300 eukaryotic, prokaryotic, and viral genomes. Using an advanced classification procedure, which is based on high-order hidden Markov models and considers both sequence similarity as well as the number of transmembrane helices and loop lengths, we identified 1353 structurally homogeneous clusters roughly corresponding to membrane protein folds. Only 53 clusters are associated with experimentally determined three-dimensional structures, and for these clusters CAMPS is in reasonable agreement with structure-based classification approaches such as SCOP and CATH. We therefore estimate that ∼1300 structures would need to be determined to provide a sufficient structural coverage of polytopic membrane proteins. CAMPS 2.0 is available at http://webclu.bio.wzw.tum.de/CAMPS2.0/.