Brevundimonas brasiliensis sp. nov.: a New Multidrug-Resistant Species Isolated from a Patient in Brazil

To increase knowledge on Brevundimonas pathogens, we conducted in-depth genomic and phenotypic characterization of a Brevundimonas strain isolated from the cerebrospinal fluid of a patient admitted in a neonatal intensive care unit. The strain was identified as a member of the genus Brevundimonas based on Vitek 2 system results and 16S rRNA gene sequencing and presented a multidrug resistance profile (MDR). Several molecular and biochemical tests were used to characterize and identify the species for in-depth results. The draft genome assembly of the isolate has a total length of 3,261,074 bp and a G+C of 66.86%, similar to other species of the genus. Multilocus sequence analysis, Type (Strain) Genome Server, digital DNA-DNA hybridization, and average nucleotide identity confirmed that the Brevundimonas sp. studied represents a distinct species, for which we propose the name Brevundimonas brasiliensis sp. nov. In silico analysis detected antimicrobial resistance genes (AMRGs) mediating resistance to β-lactams (penP, blaTEM-16, and blaBKC-1) and aminoglycosides [strA, strB, aac(6')-Ib, and aac(6')-Il]. We also found AMRGs encoding the AcrAB efflux pump that confers resistance to a broad spectrum of antibiotics. Colistin and quinolone resistance can be attributed to mutation in qseC and/or phoP and GyrA/GyrB, respectively. The Brevundimonas brasiliensis sp. nov. genome contained copies of type IV secretion system (T4SS)-type integrative and conjugative elements (ICEs); integrative mobilizable elements (IME); and Tn3-type and IS3, IS6, IS5, and IS1380 families, suggesting an important role in the development and dissemination of antibiotic resistance. The isolate presented a range of virulence-associated genes related to biofilm formation, adhesion, and invasion that can be relevant for its pathogenicity. Our findings provide a wealth of data to hinder the transmission of MDR Brevundimonas and highlight the need for monitoring and identifying new bacterial species in hospital environments. IMPORTANCE Brevundimonas species is considered an opportunistic human pathogen that can cause multiple types of invasive and severe infections in patients with underlying pathologies. Treatment of these pathogens has become a major challenge because many isolates are resistant to most antibiotics used in clinical practice. Furthermore, there are no consistent therapeutic results demonstrating the efficacy of antibacterial agents. Although considered a rare pathogen, recent studies have provided evidence of the emergence of Brevundimonas in clinical settings. Hence, we identified a novel pathogenic bacterium, Brevundimonas brasiliensis sp. nov., that presented a multidrug resistance (MDR) profile and carried diverse genes related to drug resistance, virulence, and mobile genetic elements. Such data can serve as a baseline for understanding the genomic diversity, adaptation, evolution, and pathogenicity of MDR Brevundimonas.

Whole genome sequencing of the multidrug-resistant Chryseobacterium indologenes isolated from a patient in Brazil

Chryseobacterium indologenes is a non-glucose-fermenting Gram-negative bacillus. This emerging multidrug resistant opportunistic nosocomial pathogen can cause severe infections in neonates and immunocompromised patients. This study aimed to present the first detailed draft genome sequence of a multidrug-resistant C. indologenes strain isolated from the cerebrospinal fluid of an infant hospitalized at the Neonatal Intensive Care Unit of Brazilian Tertiary Hospital. We first analyzed the susceptibility of C. indologenes strain to different antibiotics using the VITEK 2 system. The strain demonstrated an outstanding resistance to all the antibiotic classes tested, including β-lactams, aminoglycosides, glycylcycline, and polymyxin. Next, C. indologenes was whole-genome-sequenced, annotated using Prokka and Rapid Annotation using Subsystems Technology (RAST), and screened for orthologous groups (EggNOG), gene ontology (GO), resistance genes, virulence genes, and mobile genetic elements using different software tools. The draft genome contained one circular chromosome of 4,836,765 bp with 37.32% GC content. The genomic features of the chromosome present numerous genes related to cellular processes that are essential to bacteria. The MDR C. indologenes revealed the presence of genes that corresponded to the resistance phenotypes, including genes to β-lactamases (bla_IND–13, bla_CIA–3, bla_TEM–116, bla_OXA–209, bla_VEB–15), quinolone (mcbG), tigecycline (tet(X6)), and genes encoding efflux pumps which confer resistance to aminoglycosides (RanA/RanB), and colistin (HlyD/TolC). Amino acid substitutions related to quinolone resistance were observed in GyrA (S83Y) and GyrB (L425I and K473R). A mutation that may play a role in the development of colistin resistance was detected in lpxA (G68D). Chryseobacterium indologenes isolate harbored 19 virulence factors, most of which were involved in infection pathways. We identified 13 Genomic Islands (GIs) and some elements associated with one integrative and conjugative element (ICEs). Other elements linked to mobile genetic elements (MGEs), such as insertion sequence (ISEIsp1), transposon (Tn5393), and integron (In31), were also present in the C. indologenes genome. Although plasmids were not detected, a ColRNAI replicon type and the most resistance genes detected in singletons were identified in unaligned scaffolds. We provided a wide range of information toward the understanding of the genomic diversity of C. indologenes, which can contribute to controlling the evolution and dissemination of this pathogen in healthcare settings.

PbGP43 Genotyping Using Paraffin-Embedded Biopsies of Human Paracoccidioidomycosis Reveals a Genetically Distinct Lineage in the Paracoccidioides brasiliensis Complex

Paracoccidioidomycosis (PCM) is a systemic mycosis caused by a group of cryptic species embedded in the Paracoccidioides brasiliensis complex and Paracoccidioides lutzii. Four species were recently inferred to belong to the P. brasiliensis complex, but the high genetic diversity found in both human and environmental samples have suggested that the number of lineages may be higher. This study aimed to assess the 43-kilodalton glycoprotein genotypes (PbGP43) in paraffin-embedded samples from PCM patients to infer the phylogenetic lineages of the P. brasiliensis complex responsible for causing the infection. Formalin-fixed, paraffin-embedded (FFPE) tissue samples from patients with histopathological diagnosis of PCM were analyzed. DNAs were extracted and amplified for a region of the second exon of the PbGP43 gene. Products were sequenced and aligned with other PbGP43 sequences available. A haplotype network and the phylogenetic relationships among sequences were inferred. Amino acid substitutions were investigated regarding the potential to modify physicochemical properties in the proteins. Six phylogenetic lineages were identified as belonging to the P. brasiliensis complex. Two lineages did not group with any of the four recognized species of the complex, and, interestingly, one of them comprised only FFPE samples. A coinfection involving two lineages was found. Five parsimony-informative sites were identified and three of them showed radical non-synonymous substitutions with the potential to promote changes in the protein. This study expands the knowledge regarding the genetic diversity existing in the P. brasiliensis complex and shows the potential of FFPE samples in species identification and in detecting coinfections.

VRILLE shows high divergence among Higher Diptera flies but may retain role as transcriptional repressor of clock

In the circadian system, the clock gene vrille (vri) is an essential component of the second feedback loop, being responsible in Drosophila for the rhythmicity of the Clock (Clk) gene transcription by its repression. Here we studied vri in a fruit fly pest, the Tephritidae Anastrepha fraterculus, aimingtoinvestigate its molecular evolution and expression patterns from whole-head extracts. We used a combination of transcriptomic, genomic and gene walking strategies to sequence and characterize Afravri in male and female head transcriptomes of A. fraterculus and detected two putative isoforms that may correspond to A and D vri isoforms of Drosophila. Both isoforms produced a full-length sequence that translates to 842 amino acids. While the protein sequence showed significant divergence to orthologous sequences from other organisms, the bZIP domain was highly conserved. Molecular evolutionary analyses showed that vri in higher Diptera flies has been evolving under positive selection. A more detailed analysis showed positive selection also in Tephritidae with 29 sites evolving under positive selection in comparison with Drosophilidae. Real time expression analysis in LD and DD conditions showed cyclic expression of Afravri mRNA with oscillation opposite to AfraClk, suggesting that VRI may also behave in Anastrepha as a transcriptional repressor of Clk, providing another indication that higher Diptera might share common interlocked transcript-translation feedback loops (TTFLs) mechanisms that differ from other insects in target genes.

Expression profiles of neotropical termites reveal microbiota-associated, caste-biased genes and biotechnological targets

Termites are well recognized by their complex development trajectories, involving dynamic differentiation process between non-reproductive castes, workers and soldiers. These insects are associated with endosymbiotic microorganisms, which help in lignocellulose digestion and nitrogen metabolism. Aiming to identify genes harbouring biotechnological potential, we analyzed workers and soldiers RNA-Seq data of three neotropical termites: Heterotermes tenuis (Isoptera: Rhinotermitidae), Velocitermes heteropterus (Isoptera: Termitidae) and Cornitermes cumulans (Isoptera: Termitidae). We observed differences in the microbiota associated with each termite family, and found protists' genes in both Termitidae species. We found an opposite pattern of caste-biased gene expression between H. tenuis and the termitids studied. Moreover, the two termitids are considerably different concerning the number of differentially expressed genes (DEGs). Functional annotation indicated considerable differences in caste-biased gene content between V. heteropterus and C. cumulans, even though they share similar diet and biological niche. Among the most DEGs, we highlighted those involved in caste differentiation and cellulose digestion, which are attractive targets for studying more efficient technologies for termite control, biomass digestion and other biotechnological applications.

Characterization of KPC-Producing Serratia marcescens in an Intensive Care Unit of a Brazilian Tertiary Hospital

Serratia marcescens has emerged as an important opportunistic pathogen responsible for nosocomial and severe infections. Here, we determined phenotypic and molecular characteristics of 54 S. marcescens isolates obtained from patient samples from intensive-care-unit (ICU) and neonatal intensive-care-unit (NIUC) of a Brazilian tertiary hospital. All isolates were resistant to beta-lactam group antibiotics, and 92.6% (50/54) were not susceptible to tigecycline. Furthermore, 96.3% showed intrinsic resistance to polymyxin E (colistin), a last-resort antibiotic for the treatment of infections caused by MDR (multidrug-resistant) Gram-negative bacteria. In contrast, high susceptibility to other antibiotics such as fluoroquinolones (81.5%), and to aminoglycosides (as gentamicin 81.5%, and amikacin 85.2%) was found. Of all isolates, 24.1% were classified as MDR. The presence of resistance and virulence genes were examined by PCR and sequencing. All isolates carried KPC-carbapenemase (bla_KPC) and extended spectrum beta-lactamase bla_TEM genes, 14.8% carried bla_OXA–1, and 16.7% carried bla_{CTX–M–1group} genes, suggesting that bacterial resistance to β-lactam antibiotics found may be associated with these genes. The genes SdeB/HasF and SdeY/HasF that are associated with efflux pump mediated drug extrusion to fluoroquinolones and tigecycline, respectively, were found in 88.9%. The aac(6′)-Ib-cr variant gene that can simultaneously induce resistance to aminoglycoside and fluoroquinolone was present in 24.1% of the isolates. Notably, the virulence genes to (i) pore-forming toxin (ShlA); (ii) phospholipase with hemolytic and cytolytic activities (PhlA); (iii) flagellar transcriptional regulator (FlhD); and (iv) positive regulator of prodigiosin and serratamolide production (PigP) were present in 98.2%. The genetic relationship among the isolates determined by ERIC-PCR demonstrated that the vast majority of isolates were grouped in a single cluster with 86.4% genetic similarity. In addition, many isolates showed 100% genetic similarity to each other, suggesting that the S. marcescens that circulate in this ICU are closely related. Our results suggest that the antimicrobial resistance to many drugs currently used to treat ICU and NIUC patients, associated with the high frequency of resistance and virulence genes is a worrisome phenomenon. Our findings emphasize the importance of active surveillance plans for infection control and to prevent dissemination of these strains.

High Prevalence of Multidrug-Resistant Klebsiella pneumoniae Harboring Several Virulence and β-Lactamase Encoding Genes in a Brazilian Intensive Care Unit

Klebsiella pneumoniae is an important opportunistic pathogen that commonly causes nosocomial infections and contributes to substantial morbidity and mortality. We sought to investigate the antibiotic resistance profile, pathogenic potential and the clonal relationships between K. pneumoniae (n = 25) isolated from patients and sources at a tertiary care hospital's intensive care units (ICUs) in the northern region of Brazil. Most of K. pneumoniae isolates (n = 21, 84%) were classified as multidrug resistant (MDR) with high-level resistance to β-lactams, aminoglycosides, quinolones, tigecycline, and colistin. All the 25 isolates presented extended-spectrum beta-lactamase-producing (ESBL), including carbapenemase producers, and carried the bla KPC (100%), bla TEM (100%), bla SHV variants (n = 24, 96%), bla OXA-1 group (n = 21, 84%) and bla CTX-M-1 group (n = 18, 72%) genes. The K2 serotype was found in 4% (n = 1) of the isolates, and the K1 was not detected. The virulence-associated genes found among the 25 isolates were mrkD (n = 24, 96%), fimH-1 (n = 22, 88%), entB (100%), iutA (n = 10, 40%), ybtS (n = 15, 60%). The genes related with efflux pumps and outer membrane porins found were AcrAB (100%), tolC (n = 24, 96%), mdtK (n = 22, 88%), OmpK35 (n = 15, 60%), and OmpK36 (n = 7, 28%). ERIC-PCR was employed to determine the clonal relationship between the different isolated strains. The obtained ERIC-PCR patterns revealed that the similarity between isolates was above 70%. To determine the sequence types (STs) a multilocus sequence typing (MLST) assay was used. The results indicated the presence of high-risk international clones among the isolates. In our study, the wide variety of MDR K. pneumoniae harboring β-lactams and virulence genes strongly suggest a necessity for the implementation of effective strategies to prevent and control the spread of antibiotic resistant infections.

Odorant-binding proteins expression patterns in recently diverged species of Anastrepha fruit flies

We studied two species of closely related South American fruit flies, Anastrepha fraterculus and Anastrepha obliqua which, despite being able to interbreed, still show some ecological and reproductive differences. Because part of these differences, such as host and mate preferences, may be related to olfactory perception, we focused our investigation on the differential expression of Odorant-binding protein (OBP) gene family, which participate in initial steps of the olfactory signal transduction cascade. We investigated patterns of expression of eight OBP genes by qPCR in male and female head tissues of both species. The expression patterns of these OBPs suggest that some OBP genes are more likely involved with the location of food resources, while others seem to be associated with mate and pheromone perception. Furthermore, the expression patterns obtained at different reproductive stages indicate that OBP expression levels changed significantly after mating in males and females of both species. All eight OBP genes analyzed here showed significant levels of differential expression between A. fraterculus and A. obliqua, suggesting that they may hold important roles in their olfactory perception differences, and consequently, may potentially be involved in their differentiation.

Head Transcriptomes of Two Closely Related Species of Fruit Flies of the Anastrepha fraterculus Group Reveals Divergent Genes in Species with Extensive Gene Flow

Several fruit flies species of the Anastrepha fraterculus group are of great economic importance for the damage they cause to a variety of fleshy fruits. Some species in this group have diverged recently, with evidence of introgression, showing similar morphological attributes that render their identification difficult, reinforcing the relevance of identifying new molecular markers that may differentiate species. We investigated genes expressed in head tissues from two closely related species: A. obliqua and A. fraterculus, aiming to identify fixed single nucleotide polymorphisms (SNPs) and highly differentiated transcripts, which, considering that these species still experience some level of gene flow, could indicate potential candidate genes involved in their differentiation process. We generated multiple libraries from head tissues of these two species, at different reproductive stages, for both sexes. Our analyses indicate that the de novo transcriptome assemblies are fairly complete. We also produced a hybrid assembly to map each species' reads, and identified 67,470 SNPs in A. fraterculus, 39,252 in A. obliqua, and 6386 that were common to both species. We identified 164 highly differentiated unigenes that had a mean interspecific index ([Formula: see text]) of at least 0.94. We selected unigenes that had Ka/Ks higher than 0.5, or had at least three or more highly differentiated SNPs as potential candidate genes for species differentiation. Among these candidates, we identified proteases, regulators of redox homeostasis, and an odorant-binding protein (Obp99c), among other genes. The head transcriptomes described here enabled the identification of thousands of genes hitherto unavailable for these species, and generated a set of candidate genes that are potentially important to genetically identify species and understand the speciation process in the presence of gene flow of A. obliqua and A. fraterculus.

Molecular evolution of Odorant-binding proteins gene family in two closely related Anastrepha fruit flies

BACKGROUND

Odorant-binding proteins (OBPs) are of great importance for survival and reproduction since they participate in initial steps of the olfactory signal transduction cascade, solubilizing and transporting chemical signals to the olfactory receptors. A comparative analysis of OBPs between closely related species may help explain how these genes evolve and are maintained under natural selection and how differences in these proteins can affect olfactory responses. We studied OBP genes in the closely related species Anastrepha fraterculus and A. obliqua, which have different host preferences, using data from RNA-seq cDNA libraries of head and reproductive tissues from male and female adults, aiming to understand the speciation process occurred between them.

RESULTS

We identified 23 different OBP sequences from Anastrepha fraterculus and 24 from A. obliqua, which correspond to 20 Drosophila melanogaster OBP genes. Phylogenetic analysis separated Anastrepha OBPs sequences in four branches that represent four subfamilies: classic, minus-C, plus-C and dimer. Both species showed five plus-C members, which is the biggest number found in tephritids until now. We found evidence of positive selection in four genes and at least one duplication event that preceded the speciation of these two species. Inferences on tertiary structures of putative proteins from these genes revealed that at least one positively selected change involves the binding cavity (the odorant binding region) in the plus-C OBP50a.

CONCLUSIONS

A. fraterculus and A. obliqua have a bigger OBP repertoire than the other tephritids studied, though the total number of Anastrepha OBPs may be larger, since we studied only a limited number of tissues. The contrast of these closely related species reveals that there are several amino acid changes between the homologous genes, which might be related to their host preferences. The plus-C OBP that has one amino acid under positive selection located in the binding cavity may be under a selection pressure to recognize and bind a new odorant. The other positively selected sites found may be involved in important structural and functional changes, especially ones in which site-specific changes would radically change amino acid properties.