Genomic analyses of ciprofloxacin-resistant Neisseria gonorrhoeae isolates recovered from the largest South American metropolitan area.

We sequenced 13 Neisseria gonorrhoeae isolates exhibiting distinct susceptibility profiles and which were recovered over 12 years in the metropolitan region of São Paulo, Brazil. Whole Genome Sequencing (WGS) was performed on an Illumina MiSeq™ 2 × 300 bp paired-end reads. Bioinformatics analyses were carried out using CGE, PATRIC, and BLAST databases for manual curation of obtained genomes. Multilocus sequence typing (MLST) analysis identified seven STs, namely ST1580, ST1590, ST1901, ST1902, ST8161, ST9363, and ST15640. Moreover, a diversity of mutations was observed in MtrR/G45D-A39T, PIB/G120K-A121S, and PBP1/L421P. Mutations associated with sulfonamides (DHPS/R228S) and rifampicin (RNAP/H552N) were also detected, as well as tetracycline resistance determinants, namely rpsJ/V57M and tet(M). The results presented herein can contribute to the knowledge of N. gonorrhoeae strains circulating in Sao Paulo, Brazil.

Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil.

Almost simultaneously, several studies reported the emergence of novel SARS-CoV-2 lineages characterized by their phylogenetic and genetic distinction (1), (2), (3), (4).….

Characterization of BKC-1 class A carbapenemase from Klebsiella pneumoniae clinical isolates in Brazil.

Three Klebsiella pneumoniae clinical isolates demonstrating carbapenem resistance were recovered from different patients hospitalized at two medical centers in São Paulo, Brazil. Resistance to all ß-lactams, quinolones, and some aminoglycosides was observed for these isolates that were susceptible to polymyxin B. Carbapenem hydrolysis, which was inhibited by clavulanic acid, was observed for all K. pneumoniae isolates that belonged to the same pulsed-field gel electrophoresis (PFGE) type and a novel sequence type (ST), ST1781 (clonal complex 442 [CC442]). A 10-kb nonconjugative incompatibility group Q (IncQ) plasmid, denominated p60136, was transferred to Escherichia coli strain TOP10 cells by electroporation. The full sequencing of p60136 showed that it was composed of a mobilization system, ISKpn23, the phosphotransferase aph3A-VI, and a 941-bp open reading frame (ORF) that codified a 313-amino acid protein. This ORF was named bla BKC-1. Brazilian Klebsiella carbapenemase-1 (BKC-1) showed a pI of 6.0 and possessed the highest identity (63%) with a ß-lactamase of Sinorhizobium meliloti, an environmental bacterium. Hydrolysis studies demonstrated that purified BKC-1 not only hydrolyzed carbapenems but also penicillins, cephalosporins, and monobactams. However, the carbapenems were less efficiently hydrolyzed due to their very low kcat values (0.0016 to 0.031 s(-1)). In fact, oxacillin was the best substrate for BKC-1 (kcat /Km , 53,522.6 mM(-1) s(-1)). Here, we report a new class A carbapenemase, confirming the diversity and rapid evolution of ß-lactamases in K. pneumoniae clinical isolates.

Comparative genomics of the major fungal agents of human and animal Sporotrichosis: Sporothrix schenckii and Sporothrix brasiliensis.

BACKGROUND: The fungal genus Sporothrix includes at least four human pathogenic species. One of these species, S. brasiliensis, is the causal agent of a major ongoing zoonotic outbreak of sporotrichosis in Brazil. Elsewhere, sapronoses are caused by S. schenckii and S. globosa. The major aims on this comparative genomic study are: 1) to explore the presence of virulence factors in S. schenckii and S. brasiliensis; 2) to compare S. brasiliensis, which is cat-transmitted and infects both humans and cats with S. schenckii, mainly a human pathogen; 3) to compare these two species to other human pathogens (Onygenales) with similar thermo-dimorphic behavior and to other plant-associated Sordariomycetes. RESULTS: The genomes of S. schenckii and S. brasiliensis were pyrosequenced to 17x and 20x coverage comprising a total of 32.3 Mb and 33.2 Mb, respectively. Pair-wise genome alignments revealed that the two species are highly syntenic showing 97.5% average sequence identity. Phylogenomic analysis reveals that both species diverged about 3.8-4.9 MYA suggesting a recent event of speciation. Transposable elements comprise respectively 0.34% and 0.62% of the S. schenckii and S. brasiliensis genomes and expansions of Gypsy-like elements was observed reflecting the accumulation of repetitive elements in the S. brasiliensis genome. Mitochondrial genomic comparisons showed the presence of group-I intron encoding homing endonucleases (HE's) exclusively in S. brasiliensis. Analysis of protein family expansions and contractions in the Sporothrix lineage revealed expansion of LysM domain-containing proteins, small GTPases, PKS type1 and leucin-rich proteins. In contrast, a lack of polysaccharide lyase genes that are associated with decay of plants was observed when compared to other Sordariomycetes and dimorphic fungal pathogens, suggesting evolutionary adaptations from a plant pathogenic or saprobic to an animal pathogenic life style. CONCLUSIONS: Comparative genomic data suggest a unique ecological shift in the Sporothrix lineage from plant-association to mammalian parasitism, which contributes to the understanding of how environmental interactions may shape fungal virulence. . Moreover, the striking differences found in comparison with other dimorphic fungi revealed that dimorphism in these close relatives of plant-associated Sordariomycetes is a case of convergent evolution, stressing the importance of this morphogenetic change in fungal pathogenesis.

The Omicron Lineages BA.1 and BA.2 (Betacoronavirus SARS-CoV-2) Have Repeatedly Entered Brazil through a Single Dispersal Hub.

Brazil currently ranks second in absolute deaths by COVID-19, even though most of its population has completed the vaccination protocol. With the introduction of Omicron in late 2021, the number of COVID-19 cases soared once again in the country. We investigated in this work how lineages BA.1 and BA.2 entered and spread in the country by sequencing 2173 new SARS-CoV-2 genomes collected between October 2021 and April 2022 and analyzing them in addition to more than 18,000 publicly available sequences with phylodynamic methods. We registered that Omicron was present in Brazil as early as 16 November 2021 and by January 2022 was already more than 99% of samples. More importantly, we detected that Omicron has been mostly imported through the state of São Paulo, which in turn dispersed the lineages to other states and regions of Brazil. This knowledge can be used to implement more efficient non-pharmaceutical interventions against the introduction of new SARS-CoV variants focused on surveillance of airports and ground transportation.

Kinetics Analysis of ß-Lactams Hydrolysis by OXA-50 Variants of Pseudomonas aeruginosa.

Pseudomonas aeruginosa is an opportunist pathogen usually associated with life threatening infections and exhibits a set of intrinsic and acquired antimicrobial mechanisms. Although resistance to penicillins-like compounds is commonly associated with the chromosomal Pseudomonas-derived cephalosporinases ß-lactamase, the real contribution of OXA-50, a second chromosomally encoded ß-lactamase, remains unclear. In this study, we characterized the biochemical properties of OXA-50, OXA-488, and OXA-494. Both oxacilinases differ from OXA-50 in two amino acids each. The blaOXA-50, blaOXA-488, and blaOXA-494 were cloned into pET26b+ that was transformed into Escherichia coli DH5α strain, expressed in E. coli BL21 strain, and then purified for obtaining the hydrolytic parameters. Benzylpenicillin was the preferential substrate instead of oxacillin. Besides, OXA-488 showed a threefold increase in catalytic efficiency for benzylpenicillin, and it was twofold more efficient in hydrolyzing imipenem, compared with OXA-50, although such carbapenemase activity was considered weak. In addition, OXA-488 and OXA-494 showed an increased affinity for penicillins, which contributed to the increased catalytic efficiency against ampicillin, especially OXA-488. Chromosomally encoded resistance mechanisms are usually overshadowed by acquired mechanisms. However, understanding their real contribution is essential to comprehend the versatile profiles verified in P. aeruginosa isolates. Such information can help to choose the best therapy in a scenario of limited options.

The Role of Lebanon in the COVID-19 Butterfly Effect: The B.1.398 Example.

In the present study, we provide a retrospective genomic surveillance of the SARS-CoV-2 pandemic in Lebanon; we newly sequence the viral genomes of 200 nasopharyngeal samples collected between July 2020 and February 2021 from patients in different regions of Lebanon and from travelers crossing the Lebanese-Syrian border, and we also analyze the Lebanese genomic dataset available at GISAID. Our results show that SARS-CoV-2 infections in Lebanon during this period were shaped by the turnovers of four dominant SARS-CoV-2 lineages, with B.1.398 being the first to thoroughly dominate. Lebanon acted as a dispersal center of B.1.398 to other countries, with intercontinental transmissions being more common than within-continent. Within the country, the district of Tripoli, which was the source of 43% of the total B.1.398 sequences in our study, was identified as being an important source of dispersal in the country. In conclusion, our findings exemplify the butterfly effect, by which a lineage that emerges in a small area can be spread around the world, and highlight the potential role of developing countries in the emergence of new variants.

Emergence of Within-Host SARS-CoV-2 Recombinant Genome After Coinfection by Gamma and Delta Variants: A Case Report.

In this study, we report the first case of intra-host SARS-CoV-2 recombination during a coinfection by the variants of concern (VOC) AY.33 (Delta) and P.1 (Gamma) supported by sequencing reads harboring a mosaic of lineage-defining mutations. By using next-generation sequencing reads intersecting regions that simultaneously overlap lineage-defining mutations from Gamma and Delta, we were able to identify a total of six recombinant regions across the SARS-CoV-2 genome within a sample. Four of them mapped in the spike gene and two in the nucleocapsid gene. We detected mosaic reads harboring a combination of lineage-defining mutations from each VOC. To our knowledge, this is the first report of intra-host RNA-RNA recombination between two lineages of SARS-CoV-2, which can represent a threat to public health management during the COVID-19 pandemic due to the possibility of the emergence of viruses with recombinant phenotypes.

Phylodynamic analysis of SARS-CoV-2 spread in Rio de Janeiro, Brazil, highlights how metropolitan areas act as dispersal hubs for new variants.

During the first semester of 2021, all of Brazil has suffered an intense wave of COVID-19 associated with the Gamma variant. In July, the first cases of Delta variant were detected in the state of Rio de Janeiro. In this work, we have employed phylodynamic methods to analyse more than 1 600 genomic sequences of Delta variant collected until September in Rio de Janeiro to reconstruct how this variant has surpassed Gamma and dispersed throughout the state. After the introduction of Delta, it has initially spread mostly in the homonymous city of Rio de Janeiro, the most populous of the state. In a second stage, dispersal occurred to mid- and long-range cities, which acted as new close-range hubs for spread. We observed that the substitution of Gamma by Delta was possibly caused by its higher viral load, a proxy for transmissibility. This variant turnover prompted a new surge in cases, but with lower lethality than was observed during the peak caused by Gamma. We reason that high vaccination rates in the state of Rio de Janeiro were possibly what prevented a higher number of deaths.

Genomic insights into the versatility of the plant growth-promoting bacterium Azospirillum amazonense.

BACKGROUND: The species Azospirillum amazonense belongs to a well-known genus of plant growth-promoting bacteria. This bacterium is found in association with several crops of economic importance; however, there is a lack of information on its physiology. In this work, we present a comprehensive analysis of the genomic features of this species. RESULTS: Genes of A. amazonense related to nitrogen/carbon metabolism, energy production, phytohormone production, transport, quorum sensing, antibiotic resistance, chemotaxis/motility and bacteriophytochrome biosynthesis were identified. Noteworthy genes were the nitrogen fixation genes and the nitrilase gene, which could be directly implicated in plant growth promotion, and the carbon fixation genes, which had previously been poorly investigated in this genus. One important finding was that some A. amazonense genes, like the nitrogenase genes and RubisCO genes, were closer phylogenetically to Rhizobiales members than to species of its own order. CONCLUSION: The species A. amazonense presents a versatile repertoire of genes crucial for its plant-associated lifestyle.

Genomic analysis of carbapenem-resistant Pseudomonas aeruginosa ST143 clone showing susceptibility to broad-spectrum cephalosporins.

OBJECTIVES: Using whole-genome sequencing (WGS), we aimed to characterise a Pseudomonas aeruginosa ST143 clinical strain (Pb9) that presented resistance to meropenem and imipenem and susceptibility to piperacillin/tazobactam and broad-spectrum cephalosporins. METHODS: The antimicrobial susceptibility profile was confirmed by broth microdilution. WGS was performed using an Illumina MiSeq platform to identify possible genetic determinants of ß-lactam resistance. Transcription levels of chromosomally encoded efflux systems and oprD were evaluated by RT-qPCR. RESULTS: WGS analysis showed that no acquired carbapenemase-encoding gene was found in isolate Pb9, although mutations in the chromosomally encoded ß-lactamase genes blaOXA-488, blaPIB-1 and blaPDC-5 were observed. In addition, we detected a premature stop codon in the major porin-encoding gene oprD coupled with hyperexpression of MexAB-OprM and MexEF-OprN. CONCLUSION: Our results suggest that the ß-lactam resistance phenotype presented by strain Pb9 might be related to an association of OprD loss with hyperexpression of the efflux pump systems MexAB-OprM and MexEF-OprN. However, the contribution of OXA-488, PDC-5 and PIB-1 to this phenotype remains unclear and warrants further investigation.

Dynamic of High-Risk Acinetobacter baumannii Major Clones in a Brazilian Tertiary Hospital During a Short Time Period.

We characterized by whole-genome sequencing (WGS) six carbapenem-resistant Acinetobacter baumannii strains isolated from a Brazilian tertiary hospital during a 14-day period. The ISAba1-blaOXA-23 structure was found in the chromosome of five isolates, whereas blaOXA-72 was inserted in a 16.6-kb plasmid in two isolates. The presence of ISAba1-blaADC-like justified the high broad-spectrum cephalosporins minimal inhibitory concentrations (MICs) (MIC50, > 512 mg/L) verified in all isolates. Only minocycline (MIC50, ≤ 0.5 µg/mL), polymyxin B (MIC50, 0.5 µg/mL), and tigecycline (MIC50, 0.5 µg/mL) were in vitro active against such isolates. A diversity of other antimicrobial resistance determinants (aph(3')-VIa, aadA1, aac(3')-IIa, strA, strB, sul2, drfA1, mph(E), msr(E), tetB, and floR) was also observed, which may confer resistance to at last six distinct antimicrobial classes. Four distinct pulsed-field gel electrophoresis (PFGE) profiles were observed during the study period, which belonged to ST79/ST258 (n = 2; IC5), ST25/ST229 (n = 2; IC7), ST1 (n = 1; IC1), and ST162/ST235 (n = 1; IC4). Although the ST1 isolate that carried blaOXA-23 and blaOXA-72 was introduced in this hospital setting by a transferred patient, two clonally related ST79/ST258 isolates carrying either one of these carbapenemase encoding genes were recovered from two patients who were hospitalized within the same period of time in the same hospital unit. Finally, a good correlation between PFGE/MLST, blaOXA-51 variant, and single nucleotide polymorphisms was also observed. Here we demonstrated that distinct extensively drug-resistant A. baumannii clones can circulate in the same hospital setting during a short time period, illustrating a very complex epidemiological scenario for this priority pathogen.

Reductive evolution of virulence repertoire to drive the divergence between community- and hospital-associated methicillin-resistant Staphylococcus aureus of the ST1 lineage.

Methicillin-resistant Staphylococcus aureus (MRSA) of the ST1-SCCmecIV lineage has been associated with community-acquired (CA) infections in North America and Australia. In Brazil, multi-drug resistant ST1-SCCmecIV MRSA has emerged in hospital-associated (HA) diseases in Rio de Janeiro. To understand these epidemiological differences, genomic and phylogenetic analyses were performed. In addition, virulence assays were done for representative CA - and HA-MRSA strains. Despite the conservation of the virulence repertoire, some genes were missing in Brazilian ST1-SCCmecIV including lukSF-PV, fnbB, and several superantigen-encoded genes. Additionally, CA-MRSA lost the splDE while HA-MRSA strains conserved the complete operon. Most of these variable genes were located in mobile genetic elements (MGE). However, conservation and maintenance of MGEs were often observed despite the absence of their associated virulence markers. A Bayesian phylogenetic tree revealed the occurrence of more than one entrance of ST1 strains in Rio de Janeiro. The tree shape and chronology allowed us to infer that the hospital-associated ST1-SCCmecIV from Brazil and the community-acquired USA400 from North America are not closely related and that they might have originated from different MSSA strains that independently acquired SCCmecIV cassettes. As expected, representatives of ST1 strains from Brazil showed lower cytotoxicity and a greater ability to survive inside human host cells. We suggest that Brazilian ST1-SCCmecIV strains have adapted to the hospital setting by reducing virulence and gaining the ability to persist and survive inside host cells. Possibly, these evolutionary strategies may balance the biologic cost of retaining multiple antibiotic resistance genes.

Genomic surveillance of SARS-CoV-2 tracks early interstate transmission of P.1 lineage and diversification within P.2 clade in Brazil.

The sharp increase of COVID-19 cases in late 2020 has made Brazil the new epicenter of the ongoing SARS-CoV-2 pandemic. The novel viral lineages P.1 (Variant of Concern Gamma) and P.2, respectively identified in the Brazilian states of Amazonas and Rio de Janeiro, have been associated with potentially higher transmission rates and antibody neutralization escape. In this study, we performed the whole-genome sequencing of 185 samples isolated from three out of the five Brazilian regions, including Amazonas (North region), Rio Grande do Norte, Paraíba and Bahia (Northeast region), and Rio de Janeiro (Southeast region) in order to monitor the spread of SARS-CoV-2 lineages in Brazil in the first months of 2021. Here, we showed a widespread dispersal of P.1 and P.2 across Brazilian regions and, except for Amazonas, P.2 was the predominant lineage identified in the sampled states. We estimated the origin of P.2 lineage to have happened in February, 2020 and identified that it has differentiated into new clades. Interstate transmission of P.2 was detected since March, but reached its peak in December, 2020 and January, 2021. Transmission of P.1 was also high in December and its origin was inferred to have happened in August 2020. We also confirmed the presence of lineage P.7, recently described in the southernmost region of Brazil, to have spread across the Northeastern states. P.1, P.2 and P.7 are descended from the ancient B.1.1.28 strain, which co-dominated the first phase of the pandemic in Brazil with the B.1.1.33 strain. We also identified the occurrence of a new lineage descending from B.1.1.33 that convergently carries the E484K mutation, N.9. Indeed, the recurrent report of many novel SARS-CoV-2 genetic variants in Brazil could be due to the absence of effective control measures resulting in high SARS-CoV2 transmission rates. Altogether, our findings provided a landscape of the critical state of SARS-CoV-2 across Brazil and confirm the need to sustain continuous sequencing of the SARS-CoV-2 isolates worldwide in order to identify novel variants of interest and monitor for vaccine effectiveness.

The Emergence of the New P.4 Lineage of SARS-CoV-2 With Spike L452R Mutation in Brazil.

The emergence of several SARS-CoV-2 lineages presenting adaptive mutations is a matter of concern worldwide due to their potential ability to increase transmission and/or evade the immune response. While performing epidemiological and genomic surveillance of SARS-CoV-2 in samples from Porto Ferreira-São Paulo-Brazil, we identified sequences classified by pangolin as B.1.1.28 harboring Spike L452R mutation, in the RBD region. Phylogenetic analysis revealed that these sequences grouped into a monophyletic branch, with others from Brazil, mainly from the state of São Paulo. The sequences had a set of 15 clade defining amino acid mutations, of which six were in the Spike protein. A new lineage was proposed to Pango and it was accepted and designated P.4. In samples from the city of Porto Ferreira, P.4 lineage has been increasing in frequency since it was first detected in March 2021, corresponding to 34.7% of the samples sequenced in June, the second in prevalence after P.1. Also, it is circulating in 30 cities from the state of São Paulo, and it was also detected in one sample from the state of Sergipe and two from the state of Rio de Janeiro. Further studies are needed to understand whether P.4 should be considered a new threat.

Intra-host evolution during SARS-CoV-2 prolonged infection.

Long-term infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represents a challenge to virus dispersion and the control of coronavirus disease 2019 (COVID-19) pandemic. The reason why some people have prolonged infection and how the virus persists for so long are still not fully understood. Recent studies suggested that the accumulation of intra-host single nucleotide variants (iSNVs) over the course of the infection might play an important role in persistence as well as emergence of mutations of concern. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during prolonged infection. Thirty-three patients who remained reverse transcription polymerase chain reaction (RT-PCR) positive in the nasopharynx for on average 18 days from the symptoms onset were included in this study. Whole-genome sequences were obtained for each patient at two different time points. Phylogenetic, populational, and computational analyses of viral sequences were consistent with prolonged infection without evidence of coinfection in our cohort. We observed an elevated within-host genomic diversity at the second time point samples positively correlated with cycle threshold (Ct) values (lower viral load). Direct transmission was also confirmed in a small cluster of healthcare professionals that shared the same workplace by the presence of common iSNVs. A differential accumulation of missense variants between the time points was detected targeting crucial structural and non-structural proteins such as Spike and helicase. Interestingly, longitudinal acquisition of iSNVs in Spike protein coincided in many cases with SARS-CoV-2 reactive and predicted T cell epitopes. We observed a distinguishing pattern of mutations over the course of the infection mainly driven by increasing A→U and decreasing G→A signatures. G→A mutations may be associated with RNA-editing enzyme activities; therefore, the mutational profiles observed in our analysis were suggestive of innate immune mechanisms of the host cell defense. Therefore, we unveiled a dynamic and complex landscape of host and pathogen interaction during prolonged infection of SARS-CoV-2, suggesting that the host's innate immunity shapes the increase of intra-host diversity. Our findings may also shed light on possible mechanisms underlying the emergence and spread of new variants resistant to the host immune response as recently observed in COVID-19 pandemic.

Pervasive transmission of E484K and emergence of VUI-NP13L with evidence of SARS-CoV-2 co-infection events by two different lineages in Rio Grande do Sul, Brazil.

Emergence of novel SARS-CoV-2 lineages are under the spotlight of the media, scientific community and governments. Recent reports of novel variants in the United Kingdom, South Africa and Brazil (B.1.1.28-E484K) have raised intense interest because of a possible higher transmission rate or resistance to the novel vaccines. Nevertheless, the spread of B.1.1.28 (E484K) and other variants in Brazil is still unknown. In this work, we investigated the population structure and genomic complexity of SARS-CoV-2 in Rio Grande do Sul, the southernmost state in Brazil. Most samples sequenced belonged to the B.1.1.28 (E484K) lineage, demonstrating its widespread dispersion. We were the first to identify two independent events of co-infection caused by the occurrence of B.1.1.28 (E484K) with either B.1.1.248 or B.1.91 lineages. Also, clustering analysis revealed the occurrence of a novel cluster of samples circulating in the state (named VUI-NP13L) characterized by 12 lineage-defining mutations. In light of the evidence for E484K dispersion, co-infection and emergence of VUI-NP13 L in Rio Grande do Sul, we reaffirm the importance of establishing strict and effective social distancing measures to counter the spread of potentially more hazardous SARS-CoV-2 strains.

Maternal SARS-CoV-2 Infection Associated to Systemic Inflammatory Response and Pericardial Effusion in the Newborn: A Case Report.

Vertical transmission of SARS-CoV-2 has already been described, while clinical consequences to the fetus are still under investigation. This article reports a case of systemic fetal inflammatory response and pericardial effusion. As far as is known, this is the first case of fetal/neonatal cardiac complications related to SARS-CoV-2 infection.