COG database update: focus on microbial diversity, model organisms, and widespread pathogens

A review of computational tools for generating metagenome-assembled genomes from metagenomic sequencing data

Metagenomic sequencing provides a culture-independent avenue to investigate the complex microbial communities by constructing metagenome-assembled genomes (MAGs). A MAG represents a microbial genome by a group of sequences from genome assembly with similar characteristics. It enables us to identify novel species and understand their potential functions in a dynamic ecosystem. Many computational tools have been developed to construct and annotate MAGs from metagenomic sequencing, however, there is a prominent gap to comprehensively introduce their background and practical performance. In this paper, we have thoroughly investigated the computational tools designed for both upstream and downstream analyses, including metagenome assembly, metagenome binning, gene prediction, functional annotation, taxonomic classification, and profiling. We have categorized the commonly used tools into unique groups based on their functional background and introduced the underlying core algorithms and associated information to demonstrate a comparative outlook. Furthermore, we have emphasized the computational requisition and offered guidance to the users to select the most efficient tools. Finally, we have indicated current limitations, potential solutions, and future perspectives for further improving the tools of MAG construction and annotation. We believe that our work provides a consolidated resource for the current stage of MAG studies and shed light on the future development of more effective MAG analysis tools on metagenomic sequencing.

Discovery of nondiazotrophic Trichodesmium species abundant and widespread in the open ocean

Filamentous and colony-forming cells within the cyanobacterial genus Trichodesmium might account for nearly half of nitrogen fixation in the sunlit ocean, a critical mechanism that sustains plankton's primary productivity. Trichodesmium has long been portrayed as a diazotrophic genus. By means of genome-resolved metagenomics, here we reveal that nondiazotrophic Trichodesmium species not only exist but also are abundant and widespread in the open ocean, benefiting from a previously overlooked functional lifestyle to expand the biogeography of this prominent marine genus. Near-complete environmental genomes for those closely related candidate species reproducibly shared functional features including a lack of genes related to nitrogen fixation, hydrogen recycling, and hopanoid lipid production concomitant with the enrichment of nitrogen assimilation genes. Our results elucidate fieldwork observations of Trichodesmium cells fixing carbon but not nitrogen. The Black Queen hypothesis and burden of low-oxygen concentration requirements provide a rationale to explain gene loss linked to nitrogen fixation among Trichodesmium species. Disconnecting taxonomic signal for this genus from a microbial community's ability to fix nitrogen will help refine our understanding of the marine nitrogen balance. Finally, we are reminded that established links between taxonomic lineages and functional traits do not always hold true.

Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification

Command-line annotation software tools have continuously gained popularity compared to centralized online services due to the worldwide increase of sequenced bacterial genomes. However, results of existing command-line software pipelines heavily depend on taxon-specific databases or sufficiently well annotated reference genomes. Here, we introduce Bakta, a new command-line software tool for the robust, taxon-independent, thorough and, nonetheless, fast annotation of bacterial genomes. Bakta conducts a comprehensive annotation workflow including the detection of small proteins taking into account replicon metadata. The annotation of coding sequences is accelerated via an alignment-free sequence identification approach that in addition facilitates the precise assignment of public database cross-references. Annotation results are exported in GFF3 and International Nucleotide Sequence Database Collaboration (INSDC)-compliant flat files, as well as comprehensive JSON files, facilitating automated downstream analysis. We compared Bakta to other rapid contemporary command-line annotation software tools in both targeted and taxonomically broad benchmarks including isolates and metagenomic-assembled genomes. We demonstrated that Bakta outperforms other tools in terms of functional annotations, the assignment of functional categories and database cross-references, whilst providing comparable wall-clock runtimes. Bakta is implemented in Python 3 and runs on MacOS and Linux systems. It is freely available under a GPLv3 license at https://github.com/oschwengers/bakta. An accompanying web version is available at https://bakta.computational.bio.

Genomic Stability and Genetic Defense Systems in Dolosigranulum pigrum, a Candidate Beneficial Bacterium from the Human Microbiome

Dolosigranulum pigrum is positively associated with indicators of health in multiple epidemiological studies of human nasal microbiota. Knowledge of the basic biology of D. pigrum is a prerequisite for evaluating its potential for future therapeutic use; however, such data are very limited. To gain insight into D. pigrum's chromosomal structure, pangenome, and genomic stability, we compared the genomes of 28 D. pigrum strains that were collected across 20 years. Phylogenomic analysis showed closely related strains circulating over this period and closure of 19 genomes revealed highly conserved chromosomal synteny. Gene clusters involved in the mobilome and in defense against mobile genetic elements (MGEs) were enriched in the accessory genome versus the core genome. A systematic analysis for MGEs identified the first candidate D. pigrum prophage and insertion sequence. A systematic analysis for genetic elements that limit the spread of MGEs, including restriction modification (RM), CRISPR-Cas, and deity-named defense systems, revealed strain-level diversity in host defense systems that localized to specific genomic sites, including one RM system hot spot. Analysis of CRISPR spacers pointed to a wealth of MGEs against which D. pigrum defends itself. These results reveal a role for horizontal gene transfer and mobile genetic elements in strain diversification while highlighting that in D. pigrum this occurs within the context of a highly stable chromosomal organization protected by a variety of defense mechanisms. IMPORTANCE Dolosigranulum pigrum is a candidate beneficial bacterium with potential for future therapeutic use. This is based on its positive associations with characteristics of health in multiple studies of human nasal microbiota across the span of human life. For example, high levels of D. pigrum nasal colonization in adults predicts the absence of Staphylococcus aureus nasal colonization. Also, D. pigrum nasal colonization in young children is associated with healthy control groups in studies of middle ear infections. Our analysis of 28 genomes revealed a remarkable stability of D. pigrum strains colonizing people in the United States across a 20-year span. We subsequently identified factors that can influence this stability, including genomic stability, phage predators, the role of MGEs in strain-level variation, and defenses against MGEs. Finally, these D. pigrum strains also lacked predicted virulence factors. Overall, these findings add additional support to the potential for D. pigrum as a therapeutic bacterium.

Heterotrophic bacterial diazotrophs are more abundant than their cyanobacterial counterparts in metagenomes covering most of the sunlit ocean

Biological nitrogen fixation contributes significantly to marine primary productivity. The current view depicts few cyanobacterial diazotrophs as the main marine nitrogen fixers. Here, we used 891 Tara Oceans metagenomes derived from surface waters of five oceans and two seas to generate a manually curated genomic database corresponding to free-living, filamentous, colony-forming, particle-attached, and symbiotic bacterial and archaeal populations. The database provides the genomic content of eight cyanobacterial diazotrophs including a newly discovered population related to known heterocystous symbionts of diatoms, as well as 40 heterotrophic bacterial diazotrophs that considerably expand the known diversity of abundant marine nitrogen fixers. These 48 populations encapsulate 92% of metagenomic signal for known nifH genes in the sunlit ocean, suggesting that the genomic characterization of the most abundant marine diazotrophs may be nearing completion. Newly identified heterotrophic bacterial diazotrophs are widespread, express their nifH genes in situ, and also occur in large planktonic size fractions where they might form aggregates that provide the low-oxygen microenvironments required for nitrogen fixation. Critically, we found heterotrophic bacterial diazotrophs to be more abundant than cyanobacterial diazotrophs in most metagenomes from the open oceans and seas, emphasizing the importance of a wide range of heterotrophic populations in the marine nitrogen balance.

The current state of SubtiWiki, the database for the model organism Bacillus subtilis

Bacillus subtilis is a Gram-positive model bacterium with extensive documented annotation. However, with the rise of high-throughput techniques, the amount of complex data being generated every year has been increasing at a fast pace. Thus, having platforms ready to integrate and give a representation to these data becomes a priority. To address it, SubtiWiki (http://subtiwiki.uni-goettingen.de/) was created in 2008 and has been growing in data and viewership ever since. With millions of requests every year, it is the most visited B. subtilis database, providing scientists all over the world with curated information about its genes and proteins, as well as intricate protein–protein interactions, regulatory elements, expression data and metabolic pathways. However, there is still a large portion of annotation to be unveiled for some biological elements. Thus, to facilitate the development of new hypotheses for research, we have added a Homology section covering potential protein homologs in other organisms. Here, we present the recent developments of SubtiWiki and give a guided tour of our database and the current state of the data for this organism.

Complete Genome Sequence of Weissella confusa LM1 and Comparative Genomic Analysis

The genus Weissella is attracting an increasing amount of attention because of its multiple functions and probiotic potential. In particular, the species Weissella confusa is known to have great potential in industrial applications and exhibits numerous biological functions. However, the knowledge on this bacterium in insects is not investigated. Here, we isolated and identified W. confusa as the dominant lactic acid bacteria in the gut of the migratory locust. We named this strain W. confusa LM1, which is the first genome of an insect-derived W. confusa strain with one complete chromosome and one complete plasmid. Among all W. confusa strains, W. confusa LM1 had the largest genome. Its genome was the closest to that of W. confusa 1001271B_151109_G12, a strain from human feces. Our results provided accurate evolutionary relationships of known Weissella species and W. confusa strains. Based on genomic analysis, the pan-genome of W. confusa is in an open state. Most strains of W. confusa had the unique genes, indicating that these strains can adapt to different ecological niches and organisms. However, the variation of strain-specific genes did represent significant correlations with their hosts and ecological niches. These strains were predicted to have low potential to produce secondary metabolites. Furthermore, no antibiotic resistance genes were identified. At the same time, virulence factors associated with toxin production and secretion system were not found, indicating that W. confusa strains were not sufficient to perform virulence. Our study facilitated the discovery of the functions of W. confusa LM1 in locust biology and their potential application to locust management.

Quantitative elucidation of associations between nucleotide identity and physicochemical properties of amino acids and the functional insight

Studies on codon property would deepen our understanding of the origin of primitive life and enlighten biotechnical application. Here, we proposed a quantitative measurement of codon-amino acid association and found that seven out of 13 physicochemical properties have stronger associations with the nucleotide identity at the second codon position, indicating that protein structure and function may associate more closely with it than the other two sites. When extending the effect of codon-amino acid association to protein level, it was found that the correlation between the second codon position (measured by the relative frequencies of nucleobase T and A at this codon site) and hydrophobicity (by the form of GRAVY value) became stronger with 96% genomes having R > 0.90 and p < 1e-60. Furthermore, we revealed that informational genes encoding proteins have lower GRAVY values than operational proteins (p < 3e-37) in both prokaryotic and eukaryotic genomes. The above results reveal a complete link from codon identity (A₂ versus T₂) to amino acid property (hydrophilic versus hydrophobic) and then to protein functions (informational versus operational). Hence, our work may help to understand how the nucleotide sequence determines protein function.

Isolation and Characterization of Outer Membrane Vesicles of Pectobacterium brasiliense 1692

Pectobacterium brasiliense (Pbr) 1692 is an aggressive phytopathogen affecting a broad host range of crops and ornamental plants, including potatoes. Previous research on animal pathogens, and a few plant pathogens, revealed that Outer Membrane Vesicles (OMVs) are part of Gram-negative bacteria's (GNB) adaptive toolkit. For this reason, OMV production and subsequent release from bacteria is a conserved process. Therefore, we hypothesized that OMVs might transport proteins that play a critical role in causing soft rot disease and in the survival and fitness of Pbr1692. Here, we show that the potato pathogen, Pbr1692, releases OMVs of various morphologies in Luria Bertani media at 31 °C. Scanning Electron Microscopy (SEM) and Transmission Electron Microscopy (TEM) confirmed the production of OMVs by Pbr1692 cells. Transmission Electron Microscopy showed that these exist as chain-, single-, and double-membrane morphologies. Mass spectrometry followed by Gene Ontology, Clusters of Orthologous Groups, Virulence Factor, CAZymes, Antibiotic Resistance Ontology, and Bastion6 T6SE annotations identified 129 OMV-associated proteins with diverse annotated roles, including antibiotic stress response, virulence, and competition. Pbr1692 OMVs contributed to virulence in potato tubers and elicited a hypersensitive response in Nicotiana benthamiana leaves. Furthermore, Pbr1692 OMVs demonstrated antibacterial activity against Dickeya dadantii.

Gulosibacter hominis sp. nov.: a novel human microbiome bacterium that may cause opportunistic infections

We present genomic, phylogenomic, and phenotypic taxonomic data to demonstrate that three human ear isolates represent a novel species within the genus Gulosibacter. These isolates could not be identified reliably using MALDI-TOF mass spectrometry during routine diagnostic work, but partial 16S rRNA gene sequence analysis revealed that they belonged to the genus Gulosibacter. Overall genomic relatedness indices between the draft genome sequences of the three isolates and of the type strains of established Gulosibacter species confirmed that the three isolates represented a single novel Gulosibacter species. A biochemical characterisation yielded differential tests between the novel and established Gulosibacter species, which could also be differentiated using MALDI-TOF mass spectrometry. We propose to formally classify these three isolates into Gulosibacter hominis sp. nov., with 401352-2018 ^T (= LMG 31778 ^T, CCUG 74795 ^T) as the type strain. The whole-genome sequence of strain 401352-2018 ^T has a size of 2,340,181 bp and a G+C content of 62.04 mol%. A Gulosibacter pangenome analysis revealed 467 gene clusters that were exclusively present in G. hominis genomes. While these G. hominis specific gene clusters were enriched in several COG functional categories, this analysis did not reveal functions that suggested a role in the human microbiome, nor did it explain the occurrence of G. hominis in ear infections. The absence of acquired antimicrobial resistance determinants and virulence factors in the G. hominis genomes, and an analysis of publicly available 16S rRNA gene sequences and 16S rRNA amplicon sequencing data sets suggested that G. hominis is a member of the human skin microbiota that may occasionally be involved in opportunistic infections.

Paraburkholderia dioscoreae sp. nov., a novel plant associated growth promotor

A novel bacterium, designated strain Msb3T, was recently isolated from leaves of the yam family plant Dioscorea bulbifera (Dioscoreaceae). Phylogenetic analysis based on the 16S rRNA gene sequence indicated that this strain belonged to the genus Paraburkholderia with Paraburkholderia xenovorans as nearest validly named neighbour taxon (99.3 % sequence similarity towards the P. xenovorans type strain). Earlier genome sequence analysis revealed a genome of 8.35 Mb in size with a G+C content of 62.5 mol%, which was distributed over two chromosomes and three plasmids. Here, we confirm that strain Msb3T represents a novel Paraburkholderia species. In silico DNA-DNA hybridization and average nucleotide identity (OrthoANIu) analyses towards P. xenovorans LB400T yielded 58.4 % dDDH and 94.5 % orthoANIu. Phenotypic and metabolic characterization revealed growth at 15 °C on tryptic soy agar, growth in the presence of 1 % NaCl and the lack of assimilation of phenylacetic acid as distinctive features. Together, these data demonstrate that strain Msb3T represents a novel species of the genus Paraburkholderia, for which we propose the name Paraburkholderia dioscoreae sp. nov. The type strain is Msb3T (=LMG 31881T, DSM 111632T, CECT 30342T).

Forty Years of Oxalobacter formigenes, a Gutsy Oxalate-Degrading Specialist

Oxalobacter formigenes, a unique anaerobic bacterium that relies solely on oxalate for growth, is a key oxalate-degrading bacterium in the mammalian intestinal tract. Degradation of oxalate in the gut by O. formigenes plays a critical role in preventing renal toxicity in animals that feed on oxalate-rich plants. The role of O. formigenes in reducing the risk of calcium oxalate kidney stone disease and oxalate nephropathy in humans is less clear, in part due to difficulties in culturing this organism and the lack of studies which have utilized diets in which the oxalate content is controlled. Herein, we review the literature on the 40th anniversary of the discovery of O. formigenes, with a focus on its biology, its role in gut oxalate metabolism and calcium oxalate kidney stone disease, and potential areas of future research. Results from ongoing clinical trials utilizing O. formigenes in healthy volunteers and in patients with primary hyperoxaluria type 1 (PH1), a rare but severe form of calcium oxalate kidney stone disease, are also discussed. Information has been consolidated on O. formigenes strains and best practices to culture this bacterium, which should serve as a good resource for researchers.

Genome Analysis of Streptomyces nojiriensis JCM 3382 and Distribution of Gene Clusters for Three Antibiotics and an Azasugar across the Genus Streptomyces

Streptomyces spp. have been major contributors of novel natural products that are used in many application areas. We found that the nojirimycin (NJ) producer JCM 3382 has antimicrobial activity against Staphylococcus aureus via cellular degradation. Genome analysis revealed 30 biosynthetic gene clusters, including those responsible for producing antibiotics, including an azasugar NJ. In-depth MS/MS analysis confirmed the production of 1-deoxynojirimycin (DNJ) along with NJ. In addition, the production of tambromycins, setomimycin, and linearmycins was verified by spectroscopic analyses, including LC-MS and NMR. The distribution of the clusters of genes coding for antibiotics in 2061 Streptomyces genomes suggested potential producers of tambromycin, setomimycin, and linearmycin. For a DNJ gene cluster, homologs of gabT1 and gutB1 were commonly found; however, yktC1 was identified in only 112 genomes. The presence of several types of clusters suggests that different strains may produce different types of azasugars. Chemical-profile-inspired comparative genome analysis may facilitate a more accurate assessment of the biosynthetic potential to produce secondary metabolites.

Genomic Analysis and Secondary Metabolites Production of the Endophytic Bacillus velezensis Bvel1: A Biocontrol Agent against Botrytis cinerea Causing Bunch Rot in Post-Harvest Table Grapes

Botrytis bunch rot caused by Botrytis cinerea is one of the most economically significant post-harvest diseases of grapes. In the present study, we showed that the bacterial strain Bvel1 is phylogenetically affiliated to Bacillus velezensis species. The strain Bvel1 and its secreted metabolites exerted an antifungal activity, under in vitro conditions, against B. cinerea. UHPLC-HRMS chemical analysis revealed that iturin A2, surfactin-C13 and -C15, oxydifficidin, bacillibactin, L-dihydroanticapsin, and azelaic acid were among the metabolites secreted by Bvel1. Treatment of wounded grape berries with Bacillus sp. Bvel1 cell culture was effective for controlling grey mold ingress and expansion in vivo. The effectiveness of this biological control agent was a function of the cell culture concentration of the antagonist applied, while preventive treatment proved to be more effective compared to curative. The strain Bvel1 exhibited an adequate colonization efficiency in wounded grapes. The whole-genome phylogeny, combined with ANI and dDDH analyses, provided compelling evidence that the strain Bvel1 should be taxonomically classified as Bacillus velezensis. Genome mining approaches showed that the strain Bvel1 harbors 13 antimicrobial biosynthetic gene clusters, including iturin A, fengycin, surfactin, bacilysin, difficidin, bacillaene, and bacillibactin. The results provide new insights into the understanding of the endophytic Bacillus velezensis Bvel1 biocontrol mechanism against post-harvest fungal pathogens, including bunch rot disease in grape berries.

Production of purple Ma bamboo (Dendrocalamus latiflorus Munro) with enhanced drought and cold stress tolerance by engineering anthocyanin biosynthesis

Overexpression of the leaf color (Lc) gene in Ma bamboo substantially increased the accumulation level of anthocyanin, and improved plant tolerance to cold and drought stresses, probably due to the increased antioxidant capacity. Most bamboos, including Ma bamboo (Dendrocalamus latiflorus Munro), are naturally evergreen and sensitive to cold and drought stresses, while it's nearly impossible to make improvements through conventual breeding due to their long and irregular flowering habit. Moreover, few studies have reported bamboo germplasm innovation through genetic engineering as bamboo genetic transformation remains difficult. In this study, we have upregulated anthocyanin biosynthesis in Ma bamboo, to generate non-green Ma bamboo with increased abiotic stress tolerance. By overexpressing the maize Lc gene, a bHLH transcription activator involved in the anthocyanin biosynthesis in Ma bamboo, we generated purple bamboos with increased anthocyanin levels including cyanidin-3-O-rutinoside, peonidin 3-O-rutinoside, and an unknown cyanidin pentaglycoside derivative. The expression levels of 9 anthocyanin biosynthesis genes were up-regulated. Overexpression of the Lc gene improved the plant tolerance to cold and drought stress, probably due to increased antioxidant capacity. The levels of the cold- and drought-related phytohormone jasmonic acid in the transgenic plants were also enhanced, which may also contribute to the plant stress-tolerant phenotypes. High anthocyanin accumulation level did not affect plant growth. Transcriptomic analysis showed higher expressions of genes involved in the flavonoid pathway in Lc transgenic bamboos compared with those in wild-type ones. The anthocyanin-rich bamboos generated here provide an example of ornamental and multiple agronomic trait improvements by genetic engineering in this important grass species.

Capabilities of bioinformatics tools for optimizing physicochemical features of proteins used in Nano biosensors: A short overview of the tools related to bioinformatics

Protein-protein ligand is one of the most detection methods used in Nano biosensors. Based on the advantage of specific docking between two special 3D structures, they have become a potent candidate in bioanalysis and Nanodiagnostic tools. These tools lease users to do a simple, fast, cost-effective, sensitive, and specific detection of molecular biomarkers in real samples. Recent advantages of using protein-protein ligand Nano-biosensors application is remarkable due to its special docking that refers to each protein unique 3D conformation. However, it challenges different problems such as low rate of docking and hard process for fixation on the basic layer. These challenges make developers to optimize the structure and functions of proteins. The process has different Nano scale calculation that could be done with algorithms and solutions are available as bioinformatics tools. This article aimed to have a short overview of the abilities of bioinformatics tools for modeling and optimization of physiochemical features of proteins in Nano scale.

•

Nano biosensors use different strategies which based on docking between two molecules to detect and identify different proteins.

•

Molecular docking between transducer in Nano biosensors and proteins rely on physicochemical features of transducer, protein and docking strategy.

•

Nano bioinformatics use bioinformatics tools and algorithms as a collective solution for developing functional structure in Nano scale.

•

Nano bioinformatics use different bioinformatics tools to optimize physicochemical features of proteins as a new approach in Nano biosensors and drug discovery.

Genetic Responses of Metabolically Active Limnospira indica Strain PCC 8005 Exposed to γ-Radiation during Its Lifecycle

Two morphotypes of the cyanobacterial Limnospira&nbsp;indica (formerly Arthrospira sp.) strain PCC 8005, denoted as P2 (straight trichomes) and P6 (helical trichomes), were subjected to chronic gamma radiation from spent nuclear fuel (SNF) rods at a dose rate of ca. 80 Gy·h^-1 for one mass doubling period (approximately 3 days) under continuous light with photoautotrophic metabolism fully active. Samples were taken for post-irradiation growth recovery and RNA-Seq transcriptional analysis at time intervals of 15, 40, and 71.5 h corresponding to cumulative doses of ca. 1450, 3200, and 5700 Gy, respectively. Both morphotypes, which were previously reported by us to display different antioxidant capacities and differ at the genomic level in 168 SNPs, 48 indels and 4 large insertions, recovered equally well from 1450 and 3200 Gy. However, while the P2 straight type recovered from 5700 Gy by regaining normal growth within 6 days, the P6 helical type took about 13 days to recover from this dose, indicating differences in their radiation tolerance and response. To investigate these differences, P2 and P6 cells exposed to the intermediate dose of gamma radiation (3200 Gy) were analyzed for differential gene expression by RNA-Seq analysis. Prior to batch normalization, a total of 1553 genes (887 and 666 of P2 and P6, respectively, with 352 genes in common) were selected based on a two-fold change in expression and a false discovery rate FDR smaller or equal to 0.05. About 85% of these 1553 genes encoded products of yet unknown function. Of the 229 remaining genes, 171 had a defined function while 58 genes were transcribed into non-coding RNA including 21 tRNAs (all downregulated). Batch normalization resulted in 660 differentially expressed genes with 98 having a function and 32 encoding RNA. From PCC 8005-P2 and PCC 8005-P6 expression patterns, it emerges that although the cellular routes used by the two substrains to cope with ionizing radiation do overlap to a large extent, both strains displayed a distinct preference of priorities.

Major histocompatibility complex class II genetic diversity and the genetic influence on gut microbiota in Guizhou minipigs

Major histocompatibility complex (MHC) is an important complex that presents antigen to T cells. The second exon of swine MHC (SLA) class II genes has antigen binding sites that bind with extracellular antigen. Populations with high MHC gene diversity result in low gut microbiota, and individuals with MHC gene heterozygote have lower gut microbiota diversity than that of homozygote. The pig is an animal with organs physiologically and anatomically similar to humans than any other mammal, and the pig is also suitably developed as a laboratory animal to establish the animal models of human disease. However, the relationship between SLA genetic diversity and the gut microbes of the pig is ambiguous. We studied the characterization of SLA class II genes and calculated the genetic diversity, and then we selected experimental animal groups divided by different SLA genotypes to investigate the gut microbiota composition by sequencing V3 to V4 hypervariable regions of bacterial 16 s rRNA from fecal samples. Our results showed that Guizhou minipigs had a low SLA genetic diversity, which may be due to the small founder population. The Guizhou minipig population deviated from neutral selection and balancing selection, which shows that Guizhou minipigs experience a strong artificial selection in recent years. We observed that the sex differences influenced gut microbiota much more deeply than that of genetics. Our results also showed that the individual with heterozygote of genes at the SLA class II region had much higher abundant gut microbiota than that of the homozygote.

Functional Redundancy and Specialization of the Conserved Cold Shock Proteins in Bacillus subtilis

Many bacteria encode so-called cold shock proteins. These proteins are characterized by a conserved protein domain. Often, the bacteria have multiple cold shock proteins that are expressed either constitutively or at low temperatures. In the Gram-positive model bacterium Bacillussubtilis, two of three cold shock proteins, CspB and CspD, belong to the most abundant proteins suggesting a very important function. To get insights into the role of these highly abundant proteins, we analyzed the phenotypes of single and double mutants, tested the expression of the csp genes and the impact of CspB and CspD on global gene expression in B. subtilis. We demonstrate that the simultaneous loss of both CspB and CspD results in a severe growth defect, in the loss of genetic competence, and the appearance of suppressor mutations. Overexpression of the third cold shock protein CspC could compensate for the loss of CspB and CspD. The transcriptome analysis revealed that the lack of CspB and CspD affects the expression of about 20% of all genes. In several cases, the lack of the cold shock proteins results in an increased read-through at transcription terminators suggesting that CspB and CspD might be involved in the control of transcription termination.

High-quality chromosome-level genomes of Cucumis metuliferus and Cucumis melo provide insight into Cucumis genome evolution

Cucumis metuliferus (African horned cucumber), a wild relative of Cucumis sativus (cucumber) and Cucumis melo (melon), displays high-level resistance to several important plant pathogens (e.g., root-knot nematodes and several viruses). Here, we report a chromosome-level genome assembly for C. metuliferus, with a 316 Mb genome sequence comprising 29 039 genes. Phylogenetic analysis of related species in family Cucurbitaceae indicated that the divergence time between C. metuliferus and melon was 17.8 million years ago. Comparisons between the C. metuliferus and melon genomes revealed large structural variations (inversions and translocations >1 Mb) in eight chromosomes of these two species. Gene family comparison showed that C. metuliferus has the largest number of resistance-related nucleotide-binding site leucine-rich repeat (NBS-LRR) genes in Cucurbitaceae. The loss of NBS-LRR loci caused by large insertions or deletions (indels) and pseudogenization caused by small indels explained the loss of NBS-LRR genes in Cucurbitaceae. Population structure analysis suggested that C. metuliferus originated in Zimbabwe, then spread to other southern African regions where it likely underwent similar domestic selection as melon. This C. metuliferus reference sequence will accelerate the understanding of the molecular evolution of resistance-related genes and enhance cucurbit crop improvement efforts.

Comparative Genomic Insights Into the Taxonomic Classification, Diversity, and Secondary Metabolic Potentials of Kitasatospora, a Genus Closely Related to Streptomyces

While the genus Streptomyces (family Streptomycetaceae) has been studied as a model for bacterial secondary metabolism and genetics, its close relatives have been less studied. The genus Kitasatospora is the second largest genus in the family Streptomycetaceae. However, its taxonomic position within the family remains under debate and the secondary metabolic potential remains largely unclear. Here, we performed systematic comparative genomic and phylogenomic analyses of Kitasatospora. Firstly, the three genera within the family Streptomycetaceae (Kitasatospora, Streptomyces, and Streptacidiphilus) showed common genomic features, including high G + C contents, high secondary metabolic potentials, and high recombination frequencies. Secondly, phylogenomic and comparative genomic analyses revealed phylogenetic distinctions and genome content differences among these three genera, supporting Kitasatospora as a separate genus within the family. Lastly, the pan-genome analysis revealed extensive genetic diversity within the genus Kitasatospora, while functional annotation and genome content comparison suggested genomic differentiation among lineages. This study provided new insights into genomic characteristics of the genus Kitasatospora, and also uncovered its previously underestimated and complex secondary metabolism.

Non-essential ribosomal proteins in bacteria and archaea identified using COGs

Ribosomal proteins (RPs) are highly conserved across the bacterial and archaeal domains. Although many RPs are essential for survival, genome analysis demonstrates the absence of some RP genes in many bacterial and archaeal genomes. Furthermore, global transposon mutagenesis and/or targeted deletion showed that elimination of some RP genes had only a moderate effect on the bacterial growth rate. Here, we systematically analyze the evolutionary conservation of RPs in prokaryotes by compiling the list of the ribosomal genes that are missing from one or more genomes in the recently updated version of the Clusters of Orthologous Genes (COG) database. Some of these absences occurred because the respective genes carried frameshifts, presumably, resulting from sequencing errors, while others were overlooked and not translated during genome annotation. Apart from these annotation errors, we identified multiple genuine losses of RP genes in a variety of bacteria and archaea. Some of these losses are clade-specific, whereas others occur in symbionts and parasites with dramatically reduced genomes. The lists of computationally and experimentally defined non-essential ribosomal genes show a substantial overlap, revealing a common trend in prokaryote ribosome evolution that could be linked to the architecture and assembly of the ribosomes. Thus, RPs that are located at the surface of the ribosome and/or are incorporated at a late stage of ribosome assembly are more likely to be non-essential and to be lost during microbial evolution, particularly, in the course of genome compaction.IMPORTANCEIn many prokaryote genomes, one or more ribosomal protein (RP) genes are missing. Analysis of 1,309 prokaryote genomes included in the COG database shows that only about half of the RPs are universally conserved in bacteria and archaea. In contrast, up to 16 other RPs are missing in some genomes, primarily, tiny (<1 Mb) genomes of host-associated bacteria and archaea. Ten universal and nine archaea-specific ribosomal proteins show clear patterns of lineage-specific gene loss. Most of the RPs that are frequently lost from bacterial genomes are located on the ribosome periphery and are non-essential in Escherichia coli and Bacillus subtilis These results reveal general trends and common constraints in the architecture and evolution of ribosomes in prokaryotes.

Assessment of phylo-functional coherence along the bacterial phylogeny and taxonomy

In this report we use available curated phylogenies, taxonomy, and genome annotations to assess the phylogenetic and gene content similarity associated with each different taxon and taxonomic rank. Subsequently, we employ the same data to assess the frontiers of functional coherence along the bacterial phylogeny. Our results show that within-group phylogenetic and gene content similarity of taxa in the same rank are not homogenous, and that these values show extensive overlap between ranks. Functional coherence along the 16S rRNA gene-based phylogeny was limited to 44 particular nodes presenting large variations in phylogenetic depth. For instance, the deep subtree affiliated to class Actinobacteria presented functional coherence, while the shallower family Enterobacteriaceae-affiliated subtree did not. On the other hand, functional coherence along the genome-based phylogeny delimited deep subtrees affiliated to phyla Actinobacteriota, Deinococcota, Chloroflexota, Firmicutes, and a subtree containing the rest of the bacterial phyla. The results presented here can be used to guide the exploration of results in many microbial ecology and evolution research scenarios. Moreover, we provide dedicated scripts and files that can be used to continue the exploration of functional coherence along the bacterial phylogeny employing different parameters or input data (https://git.io/Jec5U).

Genomic and Metabolomic Investigation of a Rhizosphere Isolate Streptomyces netropsis WLXQSS-4 Associated with a Traditional Chinese Medicine

Rhizosphere microorganisms play important ecological roles in promoting herb growth and producing abundant secondary metabolites. Studies on the rhizosphere microbes of traditional Chinese medicines (TCMs) are limited, especially on the genomic and metabolic levels. In this study, we reported the isolation and characterization of a Steptomyces netropsis WLXQSS-4 strain from the rhizospheric soil of Clematis manshurica Rupr. Genomic sequencing revealed an impressive total of 40 predicted biosynthetic gene clusters (BGCs), whereas metabolomic profiling revealed 13 secondary metabolites under current laboratory conditions. Particularly, medium screening activated the production of alloaureothin, whereas brominated and chlorinated pimprinine derivatives were identified through precursor-directed feeding. Moreover, antiproliferative activities against Hela and A549 cancer cell lines were observed for five compounds, of which two also elicited potent growth inhibition in Enterococcus faecalis and Staphylococcus aureus, respectively. Our results demonstrated the robust secondary metabolism of S. netropsis WLXQSS-4, which may serve as a biocontrol agent upon further investigation.

An Introduction to Next Generation Sequencing Bioinformatic Analysis in Gut Microbiome Studies

The gut microbiome is a microbial ecosystem which expresses 100 times more genes than the human host and plays an essential role in human health and disease pathogenesis. Since most intestinal microbial species are difficult to culture, next generation sequencing technologies have been widely applied to study the gut microbiome, including 16S rRNA, 18S rRNA, internal transcribed spacer (ITS) sequencing, shotgun metagenomic sequencing, metatranscriptomic sequencing and viromic sequencing. Various software tools were developed to analyze different sequencing data. In this review, we summarize commonly used computational tools for gut microbiome data analysis, which extended our understanding of the gut microbiome in health and diseases.

Density Peak clustering of protein sequences associated to a Pfam clan reveals clear similarities and interesting differences with respect to manual family annotation

BACKGROUND

The identification of protein families is of outstanding practical importance for in silico protein annotation and is at the basis of several bioinformatic resources. Pfam is possibly the most well known protein family database, built in many years of work by domain experts with extensive use of manual curation. This approach is generally very accurate, but it is quite time consuming and it may suffer from a bias generated from the hand-curation itself, which is often guided by the available experimental evidence.

RESULTS

We introduce a procedure that aims to identify automatically putative protein families. The procedure is based on Density Peak Clustering and uses as input only local pairwise alignments between protein sequences. In the experiment we present here, we ran the algorithm on about 4000 full-length proteins with at least one domain classified by Pfam as belonging to the Pseudouridine synthase and Archaeosine transglycosylase (PUA) clan. We obtained 71 automatically-generated sequence clusters with at least 100 members. While our clusters were largely consistent with the Pfam classification, showing good overlap with either single or multi-domain Pfam family architectures, we also observed some inconsistencies. The latter were inspected using structural and sequence based evidence, which suggested that the automatic classification captured evolutionary signals reflecting non-trivial features of protein family architectures. Based on this analysis we identified a putative novel pre-PUA domain as well as alternative boundaries for a few PUA or PUA-associated families. As a first indication that our approach was unlikely to be clan-specific, we performed the same analysis on the P53 clan, obtaining comparable results.

CONCLUSIONS

The clustering procedure described in this work takes advantage of the information contained in a large set of pairwise alignments and successfully identifies a set of putative families and family architectures in an unsupervised manner. Comparison with the Pfam classification highlights significant overlap and points to interesting differences, suggesting that our new algorithm could have potential in applications related to automatic protein classification. Testing this hypothesis, however, will require further experiments on large and diverse sequence datasets.

Structural, Metabolic and Evolutionary Comparison of Bacterial Endospore and Exospore Formation

Sporulation is a specialized developmental program employed by a diverse set of bacteria which culminates in the formation of dormant cells displaying increased resilience to stressors. This represents a major survival strategy for bacteria facing harsh environmental conditions, including nutrient limitation, heat, desiccation, and exposure to antimicrobial compounds. Through dispersal to new environments via biotic or abiotic factors, sporulation provides a means for disseminating genetic material and promotes encounters with preferable environments thus promoting environmental selection. Several types of bacterial sporulation have been characterized, each involving numerous morphological changes regulated and performed by non-homologous pathways. Despite their likely independent evolutionary origins, all known modes of sporulation are typically triggered by limited nutrients and require extensive membrane and peptidoglycan remodeling. While distinct modes of sporulation have been observed in diverse species, two major types are at the forefront of understanding the role of sporulation in human health, and microbial population dynamics and survival. Here, we outline endospore and exospore formation by members of the phyla Firmicutes and Actinobacteria, respectively. Using recent advances in molecular and structural biology, we point to the regulatory, genetic, and morphological differences unique to endo- and exospore formation, discuss shared characteristics that contribute to the enhanced environmental survival of spores and, finally, cover the evolutionary aspects of sporulation that contribute to bacterial species diversification.

No endospore formation confirmed in members of the phylum Proteobacteria

Endospore formation is used by members of the phylum Firmicutes to withstand extreme environmental conditions. Several recent studies have proposed endospore formation in species outside of Firmicutes, particularly in Rhodobacter johrii and Serratia marcescens, members of the phylum Proteobacteria. Here, we aimed to investigate endospore formation in these two species by using advanced imaging and analytical approaches. Examination of the phase-bright structures observed in R. johrii and S. marcescens using cryo-electron tomography failed to identify endospores or stages of endospore formation. We determined that the phase-bright objects in R. johrii cells were triacylglycerol storage granules and those in S. marcescens were aggregates of cellular debris. In addition, R. johrii and S. marcescens containing phase-bright objects do not possess phenotypic and genetic features of endospores, including enhanced resistance to heat, presence of dipicolinic acid, or the presence of many of the genes associated with endospore formation. Our results support the hypothesis that endospore formation is restricted to the phylum Firmicutes.Importance: Bacterial endospore formation is an important process that allows the formation of dormant life forms called spores. As such, organisms able to sporulate can survive harsh environmental conditions for hundreds of years. Here, we follow up on previous claims that two members of Proteobacteria, Serratia marcescens and Rhodobacter johrii, are able to form spores. We conclude that those claims were incorrect and show that the putative spores in R. johrii and S. marcescens are storage granules and cellular debris, respectively. This study concludes that endospore formation is still unique to the phylum Firmicutes.

Post-Transcriptional Modifications of Conserved Nucleotides in the T-Loop of tRNA: A Tale of Functional Convergent Evolution

The high conservation of nucleotides of the T-loop, including their chemical identity, are hallmarks of tRNAs from organisms belonging to the three Domains of Life. These structural characteristics allow the T-loop to adopt a peculiar intraloop conformation able to interact specifically with other conserved residues of the D-loop, which ultimately folds the mature tRNA in a unique functional canonical L-shaped architecture. Paradoxically, despite the high conservation of modified nucleotides in the T-loop, enzymes catalyzing their formation depend mostly on the considered organism, attesting for an independent but convergent evolution of the post-transcriptional modification processes. The driving force behind this is the preservation of a native conformation of the tRNA elbow that underlies the various interactions of tRNA molecules with different cellular components.

The 2021 Nucleic Acids Research database issue and the online molecular biology database collection

The 2021 Nucleic Acids Research database Issue contains 189 papers spanning a wide range of biological fields and investigation. It includes 89 papers reporting on new databases and 90 covering recent changes to resources previously published in the Issue. A further ten are updates on databases most recently published elsewhere. Seven new databases focus on COVID-19 and SARS-CoV-2 and many others offer resources for studying the virus. Major returning nucleic acid databases include NONCODE, Rfam and RNAcentral. Protein family and domain databases include COG, Pfam, SMART and Panther. Protein structures are covered by RCSB PDB and dispersed proteins by PED and MobiDB. In metabolism and signalling, STRING, KEGG and WikiPathways are featured, along with returning KLIFS and new DKK and KinaseMD, all focused on kinases. IMG/M and IMG/VR update in the microbial and viral genome resources section, while human and model organism genomics resources include Flybase, Ensembl and UCSC Genome Browser. Cancer studies are covered by updates from canSAR and PINA, as well as newcomers CNCdatabase and Oncovar for cancer drivers. Plant comparative genomics is catered for by updates from Gramene and GreenPhylDB. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). The NAR online Molecular Biology Database Collection has been substantially updated, revisiting nearly 1000 entries, adding 90 new resources and eliminating 86 obsolete databases, bringing the current total to 1641 databases. It is available at https://www.oxfordjournals.org/nar/database/c/.