The Soybean Expression Atlas v2: A comprehensive database of over 5000 RNA-seq samples.

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 dynamics of cellular processes. Here, we present the Soybean Expression Atlas v2, an updated version of our database with a fourfold 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 at https://soyatlas.venanciogroup.uenf.br/, and it will accelerate soybean research, empowering researchers with high-quality and easily accessible gene expression data.

Comparative genomics and phylogenomics of Campylobacter unveil potential novel species and provide insights into niche segregation.

Campylobacter is a bacterial genus associated with community outbreaks and gastrointestinal symptoms. Studies on Campylobacter generally focus on specific pathogenic species such as C. coli and C. jejuni. Currently, there are thousands of publicly available Campylobacter genomes, allowing a more complete assessment of the genus diversity. In this work, we report a network-based analysis of all available Campylobacter genomes to explore the genus structure and diversity, revealing potentially new species and elucidating genus features. We also hypothesize that the previously established Clade III of C. coli is in fact a novel species (referred here as Campylobacter spp12). Finally, we found a negative correlation between pangenome fluidity and saturation coefficient, with potential implications to the lifestyles of distinct Campylobacter species. Since pangenome analysis depends on the number of available genomes, this correlation could help estimate pangenome metrics of Campylobacter species with less sequenced genomes, helping understand their lifestyle and niche adaptation. Together, our results indicate that the Campylobacter genus should be re-evaluated, with particular attention to the interplay between genome structure and niche segregation.

BioNERO: an all-in-one R/Bioconductor package for comprehensive and easy biological network reconstruction.

Currently, standard network analysis workflows rely on many different packages, often requiring users to have a solid statistics and programming background. Here, we present BioNERO, an R package that aims to integrate all aspects of network analysis workflows, including expression data preprocessing, gene coexpression and regulatory network inference, functional analyses, and intraspecies and interspecies network comparisons. The state-of-the-art methods implemented in BioNERO ensure that users can perform all analyses with a single package in a simple pipeline, without needing to learn a myriad of package-specific syntaxes. BioNERO offers a user-friendly framework that can be easily incorporated in systems biology pipelines.

Phylogenetic analysis and population structure of Pseudomonas alloputida.

The Pseudomonas putida group comprises strains with biotechnological and clinical relevance. P. alloputida was proposed as a new species and highlighted the misclassification of P. putida. Nevertheless, the population structure of P. alloputida remained unexplored. We retrieved 11,025 Pseudomonas genomes and used P. alloputida Kh7T to delineate the species. The P. alloputida population structure comprises at least 7 clonal complexes (CCs). Clinical isolates are mainly found in CC4 and acquired resistance genes are present at low frequency in plasmids. Virulence profiles support the potential of CC7 members to outcompete other plant or human pathogens through a type VI secretion system. Finally, we found that horizontal gene transfer had an important role in shaping the ability of P. alloputida to bioremediate aromatic compounds such as toluene. Our results provide the grounds to understand P. alloputida genetic diversity and its potential for biotechnological applications.

Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae).

Transcription factors (TFs) are essential for plant growth and development. Several legumes (e.g. soybean) are rich sources of protein and oil and have great economic importance. Here we report a phylogenomic analysis of TF families in legumes and their potential association with important traits (e.g. nitrogen fixation). We used TF DNA-binding domains to systematically screen the genomes of 15 leguminous and five non-leguminous species. Transcription factor orthologous groups (OGs) were used to estimate OG sizes in ancestral nodes using a gene birth-death model, which allowed the identification of lineage-specific expansions. The OG analysis and rate of synonymous substitutions show that major TF expansions are strongly associated with whole-genome duplication (WGD) events in the legume (approximately 58 million years ago) and Glycine (approximately 13 million years ago) lineages, which account for a large fraction of the Phaseolus vulgaris and Glycine max TF repertoires. Of the 3407 G. max TFs, 1808 and 676 have homeologs within single syntenic regions in Phaseolus vulgaris and Vitis vinifera, respectively. We found a trend for TFs expanded in legumes to be preferentially transcribed in roots and nodules, supporting their recruitment early in the evolution of nodulation in the legume clade. Some families also showed count differences between G. max and the wild soybean Glycine soja, including genes located within important quantitative trait loci. Our findings strongly support the roles of two WGDs in shaping the TF repertoires in the legume and Glycine lineages, and these are probably related to important aspects of legume and soybean biology.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas.

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.

Population structure and pangenome analysis of Enterobacter bugandensis uncover the presence of blaCTX-M-55, blaNDM-5 and blaIMI-1, along with sophisticated iron acquisition strategies.

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 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 from other bacteria from the same country. Three distinctive siderophore producing operons were found in E. bugandensis: enterobactin (ent), aerobactin (iuc/iut), and salmochelin (iro). Our findings provide novel insights on the lifestyle, physiology, antimicrobial, and virulence profiles of E. bugandensis.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks.

MAIN CONCLUSION: We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Evolution and multiple roles of the Pancrustacea specific transcription factor zelda in insects.

Gene regulatory networks (GRNs) evolve as a result of the coevolutionary processes acting on transcription factors (TFs) and the cis-regulatory modules they bind. The zinc-finger TF zelda (zld) is essential for the maternal-to-zygotic transition (MZT) in Drosophila melanogaster, where it directly binds over thousand cis-regulatory modules to regulate chromatin accessibility. D. melanogaster displays a long germ type of embryonic development, where all segments are simultaneously generated along the whole egg. However, it remains unclear if zld is also involved in the MZT of short-germ insects (including those from basal lineages) or in other biological processes. Here we show that zld is an innovation of the Pancrustacea lineage, being absent in more distant arthropods (e.g. chelicerates) and other organisms. To better understand zld´s ancestral function, we thoroughly investigated its roles in a short-germ beetle, Tribolium castaneum, using molecular biology and computational approaches. Our results demonstrate roles for zld not only during the MZT, but also in posterior segmentation and patterning of imaginal disc derived structures. Further, we also demonstrate that zld is critical for posterior segmentation in the hemipteran Rhodnius prolixus, indicating this function predates the origin of holometabolous insects and was subsequently lost in long-germ insects. Our results unveil new roles of zld in different biological contexts and suggest that changes in expression of zld (and probably other major TFs) are critical in the evolution of insect GRNs.

Cell wall dynamics and gene expression on soybean embryonic axes during germination.

MAIN CONCLUSION: Identification of the structural changes and cell wall-related genes likely involved in cell wall extension, cellular water balance and cell wall biosynthesis on embryonic axes during germination of soybean seeds. Cell wall is a highly organized and dynamic structure that provides mechanical support for the cell. During seed germination, the cell wall is critical for cell growth and seedling establishment. Although seed germination has been widely studied in several species, key aspects regarding the regulation of cell wall dynamics in germinating embryonic axes remain obscure. Here, we characterize the gene expression patterns of cell wall pathways and investigate their impact on the cell wall dynamics of embryonic axes of germinating soybean seeds. We found 2143 genes involved in cell wall biosynthesis and assembly in the soybean genome. Key cell wall genes were highly expressed at specific germination stages, such as expansins, UDP-Glc epimerases, GT family, cellulose synthases, peroxidases, arabinogalactans, and xyloglucans-related genes. Further, we found that embryonic axes grow through modulation of these specific cell wall genes with no increment in biomass. Cell wall structural analysis revealed a defined pattern of cell expansion and an increase in cellulose content during germination. In addition, we found a clear correlation between these structural changes and expression patterns of cell wall genes during germination. Taken together, our results provide a better understanding of the complex transcriptional regulation of cell wall genes that drive embryonic axes growth and expansion during soybean germination.

Molecular epidemiology of 16S rRNA methyltransferase in Brazil: RmtG in Klebsiella aerogenes ST93 (CC4).

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.

Genome sequencing and assessment of plant growth-promoting properties of a Serratia marcescens strain isolated from vermicompost.

BACKGROUND: Plant-bacteria associations have been extensively studied for their potential in increasing crop productivity in a sustainable manner. Serratia marcescens is a species of Enterobacteriaceae found in a wide range of environments, including soil. RESULTS: Here we describe the genome sequencing and assessment of plant growth-promoting abilities of S. marcescens UENF-22GI, a strain isolated from mature cattle manure vermicompost. In vitro, S. marcescens UENF-22GI is able to solubilize P and Zn, to produce indole compounds (likely IAA), to colonize hyphae and counter the growth of two phytopathogenic fungi. Inoculation of maize with this strain remarkably increased seedling growth and biomass under greenhouse conditions. The S. marcescens UENF-22GI genome has 5 Mb, assembled in 17 scaffolds comprising 4662 genes (4528 are protein-coding). No plasmids were identified. S. marcescens UENF-22GI is phylogenetically placed within a clade comprised almost exclusively of non-clinical strains. We identified genes and operons that are likely responsible for the interesting plant-growth promoting features that were experimentally described. The S. marcescens UENF-22GI genome harbors a horizontally-transferred genomic island involved in antibiotic production, antibiotic resistance, and anti-phage defense via a novel ADP-ribosyltransferase-like protein and possible modification of DNA by a deazapurine base, which likely contributes to its competitiveness against other bacteria. CONCLUSIONS: Collectively, our results suggest that S. marcescens UENF-22GI is a strong candidate to be used in the enrichment of substrates for plant growth promotion or as part of bioinoculants for agriculture.

Expansion and diversification of the gibberellin receptor GIBBERELLIN INSENSITIVE DWARF1 (GID1) family in land plants.

KEY MESSAGE: Here we uncover the major evolutionary events shaping the evolution of the GID1 family of gibberellin receptors in land plants at the sequence, structure and gene expression levels. Gibberellic acid (gibberellin, GA) controls key developmental processes in the life cycle of land plants. By interacting with the GIBBERELLIN INSENSITIVE DWARF1 (GID1) receptor, GA regulates the expression of a wide range of genes through different pathways. Here we report the systematic identification and classification of GID1s in 54 plants genomes, encompassing from bryophytes and lycophytes, to several monocots and eudicots. We investigated the evolutionary relationship of GID1s using a comparative genomics framework and found strong support for a previously proposed phylogenetic classification of this family in land plants. We identified lineage-specific expansions of particular subfamilies (i.e. GID1ac and GID1b) in different eudicot lineages (e.g. GID1b in legumes). Further, we found both, shared and divergent structural features between GID1ac and GID1b subgroups in eudicots that provide mechanistic insights on their functions. Gene expression data from several species show that at least one GID1 gene is expressed in every sampled tissue, with a strong bias of GID1b expression towards underground tissues and dry legume seeds (which typically have low GA levels). Taken together, our results indicate that GID1ac retained canonical GA signaling roles, whereas GID1b specialized in conditions of low GA concentrations. We propose that this functional specialization occurred initially at the gene expression level and was later fine-tuned by mutations that conferred greater GA affinity to GID1b, including a Phe residue in the GA-binding pocket. Finally, we discuss the importance of our findings to understand the diversification of GA perception mechanisms in land plants.

Discovering and prioritizing candidate resistance genes against soybean pests by integrating GWAS and gene coexpression networks.

Soybean is one of the most important legume crops worldwide. Soybean pests have a considerable impact on crop yield. Here, we integrated publicly available genome-wide association studies and transcriptomic data to prioritize candidate resistance genes against the insects Aphis glycines and Spodoptera litura, and the nematode Heterodera glycines. We identified 171, 7, and 228 high-confidence candidate resistance genes against A. glycines, S. litura, and H. glycines, respectively. We found some overlap of candidate genes between insect species, but not between insects and H. glycines. Although 15% of the prioritized candidate genes encode proteins of unknown function, the vast majority of the candidates are related to plant immunity processes, such as transcriptional regulation, signaling, oxidative stress, recognition, and physical defense. Based on the number of resistance alleles, we selected the ten most promising accessions against each pest species in the soybean USDA germplasm. The most resistant accessions do not reach the maximum theoretical resistance potential, indicating that they might be further improved to increase resistance in breeding programs or through genetic engineering. Finally, the coexpression networks we inferred in this work are available in a user-friendly web application (https://soypestgcn.venanciogroup.uenf.br/) and an R/Shiny package (https://github.com/almeidasilvaf/SoyPestGCN) that serve as a public resource to explore soybean-pest interactions at the transcriptional level.

Comparative Genomics Reveals Novel Species and Insights into the Biotechnological Potential, Virulence, and Resistance of Alcaligenes.

Alcaligenes is a cosmopolitan bacterial genus that exhibits diverse properties which are beneficial to plants. However, the genomic versatility of Alcaligenes has also been associated with the ability to cause opportunistic infections in humans, raising concerns about the safety of these microorganisms in biotechnological applications. Here, we report an in-depth comparative analysis of Alcaligenes species using all publicly available genomes to investigate genes associated with species, biotechnological potential, virulence, and resistance to multiple antibiotics. Phylogenomic analysis revealed that Alcaligenes consists of at least seven species, including three novel species. Pan-GWAS analysis uncovered 389 species-associated genes, including cold shock proteins (e.g., cspA) and aquaporins (e.g., aqpZ) found exclusively in the water-isolated species, Alcaligenes aquatilis. Functional annotation of plant-growth-promoting traits revealed enrichment of genes for auxin biosynthesis, siderophores, and organic acids. Genes involved in xenobiotic degradation and toxic metal tolerance were also identified. Virulome and resistome profiles provide insights into selective pressures exerted in clinical settings. Taken together, the results presented here provide the grounds for more detailed clinical and ecological studies of the genus Alcaligenes.

Molecular Characterization of Salmonella Phage Wara Isolated from River Water in Brazil.

Antimicrobial resistance is increasing despite new treatments being employed, so novel strategies are required to ensure that bacterial infections remain treatable. Bacteriophages (phages; bacteria viruses) have the potential to be used as natural antimicrobial methods to control bacterial pathogens such as Salmonella spp. A Salmonella phage, Wara, was isolated from environmental water samples at the Subaé River Basin, Salvador de Bahia, Brazil. The basin has environmental impacts in its main watercourses arising from the dumping of domestic and industrial effluents and agricultural and anthropological activities. The phage genome sequence was determined by Oxford Nanopore Technologies (ONT) MinION and Illumina HiSeq sequencing, and assembly was carried out by Racon (MinION) and Unicycler (Illumina, Illumina + MinION). The genome was annotated and compared to other Salmonella phages using various bioinformatics approaches. MinION DNA sequencing combined with Racon assembly gave the best complete genome sequence. Phylogenetic analysis revealed that Wara is a member of the Tequintavirus genus. A lack of lysogeny genes, antimicrobial resistance, and virulence genes indicated that Wara has therapeutic and biocontrol potential against Salmonella species in healthcare and agriculture.

Pathogenesis-related protein 1 (PR-1) genes in soybean: Genome-wide identification, structural analysis and expression profiling under multiple biotic and abiotic stresses.

Plant pathogenesis-related (PR) proteins are a large group of proteins, classified in 17 families, that are induced by pathological conditions. Here, we characterized the soybean PR-1 (GmPR-1) gene repertoire at the sequence, structural and expression levels. We found 24 GmPR-1 genes, clustered in two phylogenetic groups. GmPR-1 genes are under strong purifying selection, particularly those that emerged by tandem duplications. GmPR-1 promoter regions are abundant in cis-regulatory elements associated with major stress-related transcription factor families, namely WRKY, ERF, HD-Zip, C2H2, NAC, and GATA. We observed that 23 GmPR-1 genes are induced by stress conditions or exclusively expressed upon stress. We explored 1972 transcriptome samples, including 26 stress conditions, revealing that most GmPR-1 genes are differentially expressed in a plethora of biotic and abiotic stresses. Our findings highlight stress-responsive GmPR-1 genes with potential biotechnological applications, such as the development of transgenic lines with increased resistance to biotic and abiotic stresses.

Network analysis of ten thousand genomes shed light on Pseudomonas diversity and classification.

The growth of sequenced bacterial genomes has revolutionized the assessment of microbial diversity. Pseudomonas is a widely diverse genus, containing more than 254 species. Although type strains have been employed to estimate Pseudomonas diversity, they represent a small fraction of the genomic diversity at a genus level. We used 10,035 available Pseudomonas genomes, including 210 type strains, to build a genomic distance network to estimate the number of species through community identification. We identified taxonomic inconsistencies with several type strains and found that 25.65 % of the Pseudomonas genomes deposited on Genbank are misclassified. The phylogenetic tree using single-copy genes from representative genomes in each species cluster in the distance network revealed at least 14 Pseudomonas groups, including the P. alcaligenes group proposed here. We show that Pseudomonas is likely an admixture of different genera and should be further divided. This study provides an overview of Pseudomonas diversity from a network and phylogenomic perspective that may help reduce the propagation of mislabeled Pseudomonas genomes.

Integrating omics approaches to discover and prioritize candidate genes involved in oil biosynthesis in soybean.

Soybean is a major source 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., BCCP2 and ACCase, FADs, KAS family proteins) and several transcription 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 homolog 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.