Search

Soybean is a crucial crop worldwide, used as a source of food, feed, and industrial products due to its high protein and oil content. Previously, the rapid accumulation of soybean RNA-seq data in public databases and the computational challenges of processing raw RNA-seq data motivated us to develop the Soybean Expression Atlas, a gene expression database of over a thousand RNA-seq samples. Over the past few years, our database has allowed researchers to explore the expression profiles of important gene families, discover genes associated with agronomic traits, and understand the transcriptional dynamic of cellular processes. Here, we present the Soybean Expression Atlas v2, an updated version of our database with a 4-fold increase in the number of samples, featuring transcript- and gene-level transcript abundance matrices for 5481 publicly available RNA-seq samples. New features in our database include the availability of transcript-level abundance estimates and equivalence classes to explore differential transcript usage, abundance estimates in bias-corrected counts to increase the accuracy of differential gene expression analyses, a new web interface with improved data visualization and user experience, and a reproducible and scalable pipeline available as an R package. The Soybean Expression Atlas v2 is available athttps://soyatlas.venanciogroup.uenf.br/, and it will accelerate soybean research, empowering researchers with high-quality and easily accessible gene expression data.

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.

Enterobacter bugandensis is a recently described species that has been largely associated with nosocomial infections. Here, we report the genome of a non-clinical E. bugandensis strain. We used this and other several publicly available E. bugandensis genomes to obtain the species pangenome, investigate the conservation of important genes, and elucidate general population structure features of the species. Core- and whole-genome multilocus sequence typing (cgMLST and wgMLST, respectively) allowed the detection of five E. bugandensis phylogroups (PG-A to E). We found important antimicrobial resistance and virulence determinants associated with specific PGs, notably PG-A and PG-E. IncFII was the most prevalent plasmid replicon type in this species. We uncovered several extended-spectrum β-lactamases, including blaCTX-M-55 and blaNDM-5, present in an IncX replicon type plasmid, described here for the first time in E. bugandensis. Genetic context analysis of blaNDM-5 revealed the resemblance of this plasmid with other IncX plasmids isolated from other bacteria from the same country. Further, three distinctive siderophore producing operons were found in the E. bugandensis pangenome: enterobactin (ent), aerobactin (iuc/iut), and salmochelin (iro). The latter operon is conserved in all PG-E isolates. Collectively, our findings provide novel insights on the lifestyle, physiology, antimicrobial, and virulence profiles of E. bugandensis.

TheStenotrophomonas maltophiliacomplex (Smc) is a cosmopolitan bacterial group that has been proposed an emergent multidrug-resistant pathogen. Taxonomic studies support the genomic heterogeneity of Smc, which comprises genogroups exhibiting a range of phenotypically distinct strains from different sources. Here, we report the genome sequencing and in-depth analysis ofS. maltophiliaUENF-4GII, isolated from vermicompost. This genome harbors a unique region encoding a penicillin-binding protein (pbpX) that was carried by a transposon, as well as horizontally-transferred genomic islands involved in anti-phage defense via DNA modification, and pili glycosylation. We also analyzed all available Smc genomes to investigate genes associated with resistance and virulence, niche occupation, and population structure.S. maltophiliaUENF-4GII belongs to genogroup 3 (Sm3), which comprises three phylogenetic clusters (PC). Pan-GWAS analysis uncovered 471 environment-associated and 791 PC-associated genes, including antimicrobial resistance (e.g.blaL1 andblaR1) and virulence determinants (e.g.treS andkatG) that provide insights on the resistance and virulence potential of Sm3 strains. Together, the results presented here provide the grounds for more detailed clinical and ecological investigations ofS. maltophilia.

Soybean is one of the major sources of edible protein and oil. Oil content is a quantitative trait that is significantly determined by genetic and environmental factors. Over the past 30 years, a large volume of soybean genetic, genomic, and transcriptomic data have been accumulated. Nevertheless, integrative analyses of such data remain scarce, in spite of their importance for crop improvement. We hypothesized that the co-occurrence of genomic regions for oil-related traits in different studies may reveal more stable regions encompassing important genetic determinants of oil content and quality in soybean. We integrated publicly available data, obtained with distinct techniques, to discover and prioritize candidate genes involved in oil biosynthesis and regulation in soybean. We detected key fatty acid biosynthesis genes (e.g., BCCP and ACCase, FADs, KAS family proteins) and several transcripton factors, which are likely regulators of oil biosynthesis. In addition, we identified new candidates for seed oil accumulation and quality, such as Glyma.03G213300 and Glyma.19G160700, which encode a translocator protein and a histone acetyltransferase, respectively. Further, oil and protein genomic hotspots are strongly associated with breeding and not with domestication, suggesting that soybean domestication prioritized other traits. The genes identified here are promising targets for breeding programs and for the development of soybean lines with increased oil content and quality.

Serratia marcescens is a global spread nosocomial pathogen. This rod-shaped bacterium displays a broad host range and worldwide geographical distribution. Here we analyze an international collection of this multidrug-resistant, opportunistic pathogen from 35 countries to infer its population structure. We show that S. marcescens comprises 12 lineages; Sm1, Sm4, and Sm10 harbor 78.3% of the known environmental strains. Sm5, Sm6, and Sm7 comprise only human-associated strains which harbor smallest pangenomes, genomic fluidity and lowest levels of core recombination, indicating niche specialization. Sm7 and Sm9 lineages exhibit the most concerning resistome; blaKPC-2 plasmid is widespread in Sm7, whereas Sm9, also an anthropogenic-exclusive lineage, presents highest plasmid/lineage size ratio and plasmid-diversity encoding metallo-beta-lactamases comprising blaNDM-1. The heterogeneity of resistance patterns of S. marcescens lineages elucidated herein highlights the relevance of surveillance programs, using whole-genome sequencing, to provide insights into the molecular epidemiology of carbapenemase producing strains of this species.

The Stenotrophomonas maltophilia complex (Smc) is a cosmopolitan bacterial group that has been proposed an emergent multidrug-resistant pathogen. Taxonomic studies support the genomic heterogeneity of Smc, which comprises genogroups exhibiting a range of phenotypically distinct strains from different sources. Here, we report the genome sequencing and in-depth analysis of S. maltophilia UENF-4GII, isolated from vermicompost. This genome harbors a unique region encoding a penicillin-binding protein (pbpX) that was carried by a transposon, as well as horizontally-transferred genomic islands involved in anti-phage defense via DNA modification, and pili glycosylation. We also analyzed all available Smc genomes to investigate genes associated with resistance and virulence, niche occupation, and population structure. S. maltophilia UENF-4GII belongs to genogroup 3 (Sm3), which comprises three phylogenetic clusters (PC). Pan-GWAS analysis uncovered 471 environment-associated and 791 PC-associated genes, including antimicrobial resistance (e.g. blaL1 and blaR1) and virulence determinants (e.g. treS and katG) that provide insights on the resistance and virulence potential of Sm3 strains. Together, the results presented here provide the grounds for more detailed clinical and ecological investigations of S. maltophilia.

Given their remarkable beneficial effects on plant growth, several Azospirillum isolates currently integrate the formulations of various commercial inoculants. Our research group isolated a new strain, Azospirillum sp. UENF-412522, from passion fruit rhizoplane. This isolate uses carbon sources that are partially distinct from closely-related Azospirillum isolates. Scanning electron microscopy analysis and population counts demonstrate the ability of Azospirillum sp. UENF-412522 to colonize the surface of passion fruit roots. In vitro assays demonstrate the ability of Azospirillum sp. UENF-412522 to fix atmospheric nitrogen, to solubilize phosphate and to produce indole-acetic acid. Passion fruit plantlets inoculated with Azospirillum sp. UENF-41255 showed increased shoot and root fresh matter, as well as root dry matter, further highlighting its biotechnological potential for agriculture. We sequenced the genome of Azospirillum sp. UENF-412522 to investigate the genetic basis of its plant-growth promotion properties. We identified the key nif genes for nitrogen fixation, the complete PQQ operon for phosphate solubilization, the acdS gene that alleviates ethylene effects on plant growth, and the napCAB operon, which produces nitrite under anoxic conditions. We also found several genes conferring resistance to common soil antibiotics, which are critical for Azospirillum sp. UENF-412522 survival in the rhizosphere. Finally, we also assessed the Azospirillum pangenome and highlighted key genes involved in plant growth promotion. A phylogenetic reconstruction of the genus was also conducted. Our results support Azospirillum sp. UENF-412522 as a good candidate for bioinoculant formulations focused on plant growth promotion in sustainable systems.

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past five years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions, and genotypes. In this study, we have collected data from 1,298 publicly available soybean transcriptome samples, processed the raw sequencing reads, and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52,737/56,044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground, and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1,349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70,963/74,490) of the known transcripts overlap with those reported here, whereas 3,256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.

Enterobacter bugandensis is a recently described species that has been largely associated with nosocomial infections. We report the genome of a non-clinical E. bugandensis strain, which was integrated with publicly available genomes to study the pangenome and general population structure of E. bugandensis. Core- and whole-genome multilocus sequence typing allowed the detection of five E. bugandensis phylogroups (PG-A to E), which contain important antimicrobial resistance and virulence determinants. We uncovered several extended-spectrum β-lactamases, including bla

Klebsiella aerogenesis an important pathogen in healthcare-associated infections. Nevertheless, in comparison to other clinically important pathogens,K. aerogenespopulation structure, genetic diversity, and pathogenicity remain poorly understood. Here, we elucidateK. aerogenesclonal complexes (CCs) and genomic features associated with resistance and virulence. We present a detailed description of the population structure ofK. aerogenesbased on 97 publicly available genomes by using both, multilocus sequence typing and single nucleotide polymorphisms extracted from core genome. We also assessed virulence and resistance profiles using VFDB and CARD, respectively. We show thatK. aerogeneshas an open pangenome and a large effective population size, which account for its high genomic diversity and support that negative selection prevents fixation of most deleterious alleles. The population is structured in at least ten CCs, including two novel ones identified here, CC9 and CC10. The repertoires of resistance genes comprise a high number of antibiotic efflux proteins as well as narrow and extended spectrum β-lactamases. Regarding the population structure, we identified two clusters based on virulence profile due to the presence of the toxin-encodingclboperon and the siderophore production genes,irpandybt.Notably, CC3 comprises the majority ofK. aerogenesisolates associated with hospital outbreaks, emphasizing the importance of its constant monitoring. Collectively, our results can be useful in the development of new therapeutic and surveillance strategies worldwide.

Aminoglycosides are a class of antibiotics that play a key role in antimicrobial treatment of Multidrug resistant (MDR) Gram-negative bacilli, typically in combination with β-lactams. Ribosomal 16S RNA modification by methyltransferases (e.g. RmtG) is an aminoglycoside resistance mechanism that, along with the occurrence carbapenem-resistant Enterobacteriaceae (CRE), has become a clinical concern. In Brazil, rmtG genes were initially reported in Klebsiella pneumoniae, and monitoring isolates from other species carrying this gene is critical for epidemiological studies and to prevent dissemination. Here we report the presence of rmtG in Klebisella aerogenes D3 and characterize its genetic context in comparison to isolates from other species. Further, we performed a phylogenetic reconstruction of 900 16S rRNA methyltransferases (16S-RMTases) and methyltransferase-related proteins. We show that, in K. aerogenes D3, rmtG co-occurs with sul2, near a transposon with an IS91-like insertion sequence. Resistome analysis revealed the co-production of RmtG and CTX-M-59. Ongoing surveillance of 16S-RMTases is crucial to delay the dissemination of such multiresistant isolates. Our results also highlight the reduction in treatment options for CRE infections, as well as the need of expanding prevention measures of these pathogens worldwide.