Expanded microbial genome coverage and improved protein family annotation in the COG database

One high quality genome and two transcriptome datasets for new species of Mantamonas, a deep-branching eukaryote clade

Mantamonads were long considered to represent an "orphan" lineage in the tree of eukaryotes, likely branching near the most frequently assumed position for the root of eukaryotes. Recent phylogenomic analyses have placed them as part of the "CRuMs" supergroup, along with collodictyonids and rigifilids. This supergroup appears to branch at the base of Amorphea, making it of special importance for understanding the deep evolutionary history of eukaryotes. However, the lack of representative species and complete genomic data associated with them has hampered the investigation of their biology and evolution. Here, we isolated and described two new species of mantamonads, Mantamonas vickermani sp. nov. and Mantamonas sphyraenae sp. nov., for each of which we generated transcriptomic sequence data, as well as a high-quality genome for the latter. The estimated size of the M. sphyraenae genome is 25 Mb; our de novo assembly appears to be highly contiguous and complete with 9,416 predicted protein-coding genes. This near-chromosome-scale genome assembly is the first described for the CRuMs supergroup.

Pangenome analysis of Corynebacterium striatum: insights into a neglected multidrug-resistant pathogen

BACKGROUND

Over the past two decades, Corynebacterium striatum has been increasingly isolated from clinical cultures with most isolates showing increased antimicrobial resistance (AMR) to last resort agents. Advances in the field of pan genomics would facilitate the understanding of the clinical significance of such bacterial species previously thought to be among commensals paving the way for identifying new drug targets and control strategies.

METHODS

We constructed a pan-genome using 310 genome sequences of C. striatum. Pan-genome analysis was performed using three tools including Roary, PIRATE, and PEPPAN. AMR genes and virulence factors have been studied in relation to core genome phylogeny. Genomic Islands (GIs), Integrons, and Prophage regions have been explored in detail.

RESULTS

The pan-genome ranges between a total of 5253-5857 genes with 2070 - 1899 core gene clusters. Some antimicrobial resistance genes have been identified in the core genome portion, but most of them were located in the dispensable genome. In addition, some well-known virulence factors described in pathogenic Corynebacterium species were located in the dispensable genome. A total of 115 phage species have been identified with only 44 intact prophage regions.

CONCLUSION

This study presents a detailed comparative pangenome report of C. striatum. The species show a very slowly growing pangenome with relatively high number of genes in the core genome contributing to lower genomic variation. Prophage elements carrying AMR and virulence elements appear to be infrequent in the species. GIs appear to offer a prominent role in mobilizing antibiotic resistance genes in the species and integrons occur at a frequency of 50% in the species. Control strategies should be directed against virulence and resistance determinants carried on the core genome and those frequently occurring in the accessory genome.

Genomic Sequence Resource of Talaromyces albobiverticillius, the Causative Pathogen of Pomegranate Pulp Rot Disease

Talaromyces albobiverticillius, a prominent pathogen responsible for pomegranate pulp rot disease, inflicts significant damage on Punica granatum L. Besides its pathogenicity, this fungus possesses the potential to produce substantial amounts of red pigments, making it promising for industrial applications. This study presents the genome annotation of T. albobiverticillius field strain Tp-2, isolated from pomegranates. The genome assembly, generated through a combination of Oxford Nanopore and Illumina sequencing reads, yielded a high-quality assembly with 14 contigs, featuring an N50 length of 4,594,200 bp. The complete genome of strain Tp-2 spans 38,354,882 bp, with a GC content of 45.78%. Importantly, the assembly exhibits remarkable integrity, with 98.3% of complete Benchmarking Universal Single-Copy Orthologs validating genome completeness. Genome prediction analysis reveals the presence of 10,380 protein-coding genes. To our knowledge, this study is the first report on the genome sequence of T. albobiverticillius, offering valuable insights into its genetic variation and molecular mechanisms of pigment production.

Protaetiibacter mangrovi sp. nov., isolated from mangrove soil

A novel Gram-stain-positive, aerobic, non-flagellated and rod-shaped actinobacterium, designated 10F1B-8-1T, was isolated from mangrove soil sampled at Futian Mangrove Nature Reserve, China. The isolate was able to grow at 10-40 °C (optimum 30-32 °C), at pH 6-8 (optimum 7) and in the presence of 0-6% (w/v) NaCl (optimum 0%). Strain 10F1B-8-1T shared the highest 16S rRNA gene sequence similarity to Protaetiibacter larvae NBRC 113051T (98.3%), followed by Protaetiibacter intestinalis NBRC 113050T (98.2%). Phylogenetic trees based on 16S rRNA gene sequences and the core proteomes exhibited that strain 10F1B-8-1T formed a new phyletic line in the clade of genus Protaetiibacter, indicating that this strain belonged to the genus Protaetiibacter. Strain 10F1B-8-1T showed low average nucleotide identity (<84%) and digital DNA-DNA hybridization values (<27%) with closely related taxa, suggesting that strain 10F1B-8-1T was a hitherto undescribed species of the genus Protaetiibacter. Strain 10F1B-8-1T contained D-2,4-diaminobutyric acid as the diagnostic diamino acid, and the peptidoglycan type was characterized as type B2β. The major fatty acids were iso-C16:0, anteiso-C15:0 and anteiso-C17:0. The major menaquinones were MK-13 and MK-14. The polar lipid profile included diphosphatidylglycerol, phosphatidylglycerol, an unidentified glycolipid and five unidentified lipids. Notably, the ethyl acetate extracts of strain 10F1B-8-1T showed effective antibacterial activity against Bacillus subtilis CPCC 100029 and Escherichia coli △tolC. According to the polyphasic data, strain 10F1B-8-1T should be classified as a novel species of the genus Protaetiibacter, for which the name Protaetiibacter mangrovi sp. nov. is proposed, with the type strain 10F1B-8-1T (=JCM 33142T = CPCC 205428T).

Validation of 16S rRNA gene sequencing and metagenomics for evaluating microbial immigration in a methanogenic bioreactor

To quantitatively evaluate the impact of microbial immigration from an upstream community on the microbial assembly of a downstream community, an ecological genomics (ecogenomics)-based mass balance (EGMB) model coupled with 16S rRNA gene sequencing was previously developed. In this study, a mock community was used to further validate the EGMB models and demonstrate the feasibility of using metagenome-based EGMB model to reveal both microbial activity and function. The mock community consisting of Aeromonas, Escherichia, and Pseudomonas was fed into a lab-scale methanogenic bioreactor together with dissolved organic substrate. Using qPCR, 16S rRNA gene, 16S rRNA gene copy number normalization (GCN), and metagenome, results showed highly comparable community profiles in the feed. In the bioreactor, Aeromonas and Pseudomonas exhibited negative growth rates throughout the experiment by all approaches. Escherichia's growth rate was negative by most biomarkers but was slightly positive by 16S rRNA gene. Still, all approaches showed a decreasing trend toward negative in the growth rate of Escherichia as reactor operation time increased. Uncultivated populations of phyla Desulfobacterota, Chloroflexi, Actinobacteriota, and Spirochaetota were observed to increase in abundance, suggesting their contribution in degrading the feed biomass. Based on metabolic reconstruction of metagenomes, these populations possessed functions of hydrolysis, fermentation, fatty acid degradation, or acetate oxidation. Overall results supported the application of both 16S rRNA gene- and metagenome-based EGMB models to measure the growth rate of microbes in the bioreactor, and the latter had advantage in providing insights into the microbial functions of uncultivated populations.

Kaistella polysaccharea sp. nov., isolated from Antarctic intertidal sediment produces a novel extracellular polymeric substance

An exopolysaccharide-producing bacterial strain GW4-15T, belonging to the genus Kaistella, was isolated from intertidal sediment from King George Island, Antarctic. The strain was Gram-stain-negative, aerobic, rod-shaped, non-motile and yellow-pigmented. The strain was able to grow in the presence of 0-2 % (w/v) NaCl (optimum, 0 %), at 4-30 °C (optimum, 20-28 °C) and at pH 5.0-10.0 (optimum, pH 8.0). A phylogenetic tree based on 16S rRNA gene sequences showed that strain GW4-15T formed a lineage within the genus Kaistella with the closest phylogenetic neighbours Kaistella carnis NCTC 13525T (98.3 %), Kaistella gelatinilytica G5-32T (97.7 %), Kaistella antarctica LMG 24720T (97.4 %) and Kaistella yonginensis HMD1043T (96.9 %). Digital DNA-DNA hybridization values of strain GW4-15T with K. carnis NCTC 13525T, K. antarctica LMG 24720T, K. gelatinilytica G5-32T and K. yonginensis HMD1043T were 22.8, 22.0, 21.7 and 21.6 %, respectively. The average nucleotide identity values between strain GW4-15T and K. carnis NCTC 13525T , K. antarctica LMG 24720T, K. gelatinilytica G5-32T and K. yonginensis HMD1043T were 79.3, 78.6, 77.5 and 77.2 %, respectively. The G+C content of the genome was 36.2 mol%. The major phospholipids were phosphatidylethanolamine and aminophospholipid. The predominant menaquinone was MK-6. The major fatty acids were anteiso-C15 : 0 (28.7 %), iso-C16 : 0 3-OH (15.7 %), iso-C16 : 0 H (10.0 %), iso-C16 : 0 (5.4 %), summed feature 9 (comprising iso-C17 : 1  ω9c and/or 10-methyl C16 : 0; 5.2 %) and iso-C15 : 0 (5.1 %). The monosaccharide composition of the new type of extracellular polymeric of GW4-15T was Glc, GalN, GlcN, Rha, Man and Gal with a molar ratio of 3.14 : 3.83 : 8.38 : 5.16 : 1 : 2.82. Based on phenotypic, phylogenetic and genotypic data, a novel species, Kaistella polysaccharea sp. nov., is proposed with the type strain GW4-15T (=CGMCC 1.19368T=KCTC 92753T).

Metabolomic rearrangement controls the intrinsic microbial response to temperature changes

Temperature is one of the key determinants of microbial behavior and survival, whose impact is typically studied under heat- or cold-shock conditions that elicit specific regulation to combat lethal stress. At intermediate temperatures, cellular growth rate varies according to the Arrhenius law of thermodynamics without stress responses, a behavior whose origins have not yet been elucidated. Using single-cell microscopy during temperature perturbations, we show that bacteria exhibit a highly conserved, gradual response to temperature upshifts with a time scale of ~1.5 doublings at the higher temperature, regardless of initial/final temperature or nutrient source. We find that this behavior is coupled to a temperature memory, which we rule out as being neither transcriptional, translational, nor membrane dependent. Instead, we demonstrate that an autocatalytic enzyme network incorporating temperature-sensitive Michaelis-Menten kinetics recapitulates all temperature-shift dynamics through metabolome rearrangement, which encodes a temperature memory and successfully predicts alterations in the upshift response observed under simple-sugar, low-nutrient conditions, and in fungi. This model also provides a mechanistic framework for both Arrhenius-dependent growth and the classical Monod Equation through temperature-dependent metabolite flux.

Navigating the Landscape: A Comprehensive Review of Current Virus Databases

Viruses are abundant and diverse entities that have important roles in public health, ecology, and agriculture. The identification and surveillance of viruses rely on an understanding of their genome organization, sequences, and replication strategy. Despite technological advancements in sequencing methods, our current understanding of virus diversity remains incomplete, highlighting the need to explore undiscovered viruses. Virus databases play a crucial role in providing access to sequences, annotations and other metadata, and analysis tools for studying viruses. However, there has not been a comprehensive review of virus databases in the last five years. This study aimed to fill this gap by identifying 24 active virus databases and included an extensive evaluation of their content, functionality and compliance with the FAIR principles. In this study, we thoroughly assessed the search capabilities of five database catalogs, which serve as comprehensive repositories housing a diverse array of databases and offering essential metadata. Moreover, we conducted a comprehensive review of different types of errors, encompassing taxonomy, names, missing information, sequences, sequence orientation, and chimeric sequences, with the intention of empowering users to effectively tackle these challenges. We expect this review to aid users in selecting suitable virus databases and other resources, and to help databases in error management and improve their adherence to the FAIR principles. The databases listed here represent the current knowledge of viruses and will help aid users find databases of interest based on content, functionality, and scope. The use of virus databases is integral to gaining new insights into the biology, evolution, and transmission of viruses, and developing new strategies to manage virus outbreaks and preserve global health.

Terabase-Scale Coassembly of a Tropical Soil Microbiome

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.

Screening and genome-wide analysis of lignocellulose-degrading bacteria from humic soil

Crop straw contains huge amounts of exploitable energy, and efficient biomass degradation measures have attracted worldwide attention. Mining strains with high yields of cellulose-degrading enzymes is of great significance for developing clean energy and industrial production of related enzymes. In this study, we reported a high-quality genome sequence of Bacillus velezensis SSF6 strain using high-throughput sequencing technology (Illumina PE150 and PacBio) and assessed its lignocellulose degradation potential. The results demonstrated that the genome of B. velezensis SSF6 was 3.89 Mb and contained 4,015 genes, of which 2,972, 3,831 and 158 genes were annotated in the COGs (Clusters of Orthologous Groups), KEGG (Kyoto Encyclopedia of Genes and Genomes) and CAZyme (Carbohydrate-Active enZymes) databases, respectively, and contained a large number of genes related to carbohydrate metabolism. Furthermore, B. velezensis SSF6 has a high cellulose degradation capacity, with a filter paper assay (FPA) and an exoglucanase activity of 64.48 ± 0.28 and 78.59 ± 0.42 U/mL, respectively. Comparative genomic analysis depicted that B. velezensis SSF6 was richer in carbohydrate hydrolase gene. In conclusion, the cellulose-degrading ability of B. velezensis SSF6 was revealed by genome sequencing and the determination of cellulase activity, which laid a foundation for further cellulose degradation and bioconversion.

Impact of the fungal pathogen Fusarium oxysporum on the taxonomic and functional diversity of the common bean root microbiome

BACKGROUND

Plants rely on their root microbiome as the first line of defense against soil-borne fungal pathogens. The abundance and activities of beneficial root microbial taxa at the time prior to and during fungal infection are key to their protective success. If and how invading fungal root pathogens can disrupt microbiome assembly and gene expression is still largely unknown. Here, we investigated the impact of the fungal pathogen Fusarium oxysporum (fox) on the assembly of rhizosphere and endosphere microbiomes of a fox-susceptible and fox-resistant common bean cultivar.

RESULTS

Integration of 16S-amplicon, shotgun metagenome as well as metatranscriptome sequencing with community ecology analysis showed that fox infections significantly changed the composition and gene expression of the root microbiome in a cultivar-dependent manner. More specifically, fox infection led to increased microbial diversity, network complexity, and a higher proportion of the genera Flavobacterium, Bacillus, and Dyadobacter in the rhizosphere of the fox-resistant cultivar compared to the fox-susceptible cultivar. In the endosphere, root infection also led to changes in community assembly, with a higher abundance of the genera Sinorhizobium and Ensifer in the fox-resistant cultivar. Metagenome and metatranscriptome analyses further revealed the enrichment of terpene biosynthesis genes with a potential role in pathogen suppression in the fox-resistant cultivar upon fungal pathogen invasion.

CONCLUSION

Collectively, these results revealed a cultivar-dependent enrichment of specific bacterial genera and the activation of putative disease-suppressive functions in the rhizosphere and endosphere microbiome of common bean under siege.

Insights into the Ecological Diversification of the Hymenochaetales based on Comparative Genomics and Phylogenomics With an Emphasis on Coltricia

To elucidate the genomic traits of ecological diversification in the Hymenochaetales, we sequenced 15 new genomes, with attention to ectomycorrhizal (EcM) Coltricia species. Together with published data, 32 genomes, including 31 Hymenochaetales and one outgroup, were comparatively analyzed in total. Compared with those of parasitic and saprophytic members, EcM species have significantly reduced number of plant cell wall degrading enzyme genes, and expanded transposable elements, genome sizes, small secreted proteins, and secreted proteases. EcM species still retain some of secreted carbohydrate-active enzymes (CAZymes) and have lost the key secreted CAZymes to degrade lignin and cellulose, while possess a strong capacity to degrade a microbial cell wall containing chitin and peptidoglycan. There were no significant differences in secreted CAZymes between fungi growing on gymnosperms and angiosperms, suggesting that the secreted CAZymes in the Hymenochaetales evolved before differentiation of host trees into gymnosperms and angiosperms. Nevertheless, parasitic and saprophytic species of the Hymenochaetales are very similar in many genome features, which reflect their close phylogenetic relationships both being white rot fungi. Phylogenomic and molecular clock analyses showed that the EcM genus Coltricia formed a clade located at the base of the Hymenochaetaceae and divergence time later than saprophytic species. And Coltricia remains one to two genes of AA2 family. These indicate that the ancestors of Coltricia appear to have originated from saprophytic ancestor with the ability to cause a white rot. This study provides new genomic data for EcM species and insights into the ecological diversification within the Hymenochaetales based on comparative genomics and phylogenomics analyses.

Complete Genome Sequence of Pantoea ananatis Strain LCFJ-001 Isolated from Bacterial Wilt Mulberry

A Pantoea ananatis strain, named LCFJ-001 (GDMCC: 1.6101), was isolated for the first time from bacterial wilt-diseased roots of mulberry (Morus atropurpurea) in the western part of the Guangxi Zhuang Autonomous Region, China. Moreover, through Koch's postulates, it was proven that LCFJ-001 can cause mulberry wilt, which is one of the pathogens of mulberry bacterial wilt. Here, we report a complete, annotated genome sequence of P. ananatis LCFJ-001. The entire genome sequence of P. ananatis strain LCFJ-001 was a 4,499,350 bp circular chromosome with 53.50% GC content. In total, 3,521 genes were annotated, of which 3,418 were assigned protein-coding genes. In addition, 22 ribosomal RNAs and 81 transfer RNAs were identified. The presented resource will help explore the pathogenetic mechanisms of mulberry wilt disease caused by the genus Pantoea.

A chromosome-level reference genome assembly of the Reeve's moray eel (Gymnothorax reevesii)

Due to potentially hostile behaviors and elusive habitats, moray eels (Muraenidae) as one group of apex predators in coral reefs all across the globe have not been well investigated. Here, we constructed a chromosome-level genome assembly for the representative Reeve's moray eel (Gymnothorax reevesii). This haplotype genome assembly is 2.17 Gb in length, and 97.87% of the sequences are anchored into 21 chromosomes. It contains 56.34% repetitive sequences and 23,812 protein-coding genes, of which 96.77% are functionally annotated. This sequenced marine species in Anguilliformes makes a good complement to the genetic resource of eel genomes. It not only provides a genetic resource for in-depth studies of the Reeve's moray eel, but also enables deep-going genomic comparisons among various eels.

MGnify Genomes: A Resource for Biome-specific Microbial Genome Catalogues

An increasingly common output arising from the analysis of shotgun metagenomic datasets is the generation of metagenome-assembled genomes (MAGs), with tens of thousands of MAGs now described in the literature. However, the discovery and comparison of these MAG collections is hampered by the lack of uniformity in their generation, annotation and storage. To address this, we have developed MGnify Genomes, a growing collection of biome-specific non-redundant microbial genome catalogues generated using MAGs and publicly available isolate genomes. Genomes within a biome-specific catalogue are organised into species clusters. For species that contain multiple conspecific genomes, the highest quality genome is selected as the representative, always prioritising an isolate genome over a MAG. The species representative sequences and annotations can be visualised on the MGnify website and the full catalogue and associated analysis outputs can be downloaded from MGnify servers. A suite of online search tools is provided allowing users to compare their own sequences, ranging from a gene to sets of genomes, against the catalogues. Seven biomes are available currently, comprising over 300,000 genomes that represent 11,048 non-redundant species, and include 36 taxonomic classes not currently represented by cultured genomes. MGnify Genomes is available at https://www.ebi.ac.uk/metagenomics/browse/genomes/.

Genomic comparison between two Inonotus hispidus strains isolated from growing in different tree species

Inonotus hispidus mainly growing in broad-leaved trees, including Morus alba, Fraxinus mandshurica, and Ulmus macrocarpa etc. The fruiting body of I. hispidus growing in M. alba (hereafter as MA) is used as a traditional Chinese medicine "Sanghuang". However, differences between the genetic material basis of I. hispidus growing in other tree species have not been reported. Therefore, in this paper, the genomic comparison between MA and I. hispidus growing in F. mandshurica (hereafter as FM) were studied. The whole genome of MA monokaryon was sequenced by Illumina combined with Pac Bio platform. Next, genome assembly, genome component prediction and genome functional annotation were performed. Comparative genomics analysis was performed between FM monokaryon and MA monokaryon, using MA as the reference. The results showed that, MA had 24 contigs with a N50 length of 2.6 Mb. Specifically, 5,342, 6,564, 1,595, 383 and 123 genes were annotated from GO, KEGG, KOG, CAZymes and CYP450, respectively. Moreover, comparative genomics showed that, the coding genes and total number of genes annotated in different databases of FM were higher than that of MA. This study provides a foundation for the medicinal application of FM as MA from the perspective of genetic composition.

Four Novel Caudoviricetes Bacteriophages Isolated from Baltic Sea Water Infect Colonizers of Aurelia aurita

The moon jellyfish Aurelia aurita is associated with a highly diverse microbiota changing with provenance, tissue, and life stage. While the crucial relevance of bacteria to host fitness is well known, bacteriophages have often been neglected. Here, we aimed to isolate virulent phages targeting bacteria that are part of the A. aurita-associated microbiota. Four phages (Pseudomonas phage BSwM KMM1, Citrobacter phages BSwM KMM2-BSwM KMM4) were isolated from the Baltic Sea water column and characterized. Phages KMM2/3/4 infected representatives of Citrobacter, Shigella, and Escherichia (Enterobacteriaceae), whereas KMM1 showed a remarkably broad host range, infecting Gram-negative Pseudomonas as well as Gram-positive Staphylococcus. All phages showed an up to 99% adsorption to host cells within 5 min, short latent periods (around 30 min), large burst sizes (mean of 128 pfu/cell), and high efficiency of plating (EOP > 0.5), demonstrating decent virulence, efficiency, and infectivity. Transmission electron microscopy and viral genome analysis revealed that all phages are novel species and belong to the class of Caudoviricetes harboring a tail and linear double-stranded DNA (formerly known as Siphovirus-like (KMM3) and Myovirus-like (KMM1/2/4) bacteriophages) with genome sizes between 50 and 138 kbp. In the future, these isolates will allow manipulation of the A. aurita-associated microbiota and provide new insights into phage impact on the multicellular host.

Mining of key genes for cold adaptation from Pseudomonas fragi D12 and analysis of its cold-adaptation mechanism

The psychrotroph Pseudomonas fragi D12, which grew strongly under low temperatures, was screened from tundra soil collected from the permanent alpine zone on Changbai Mountain. To mine the genes critical for cold tolerance and to investigate the cold-adaptation mechanism, whole-genome sequencing, comparative genomic analysis, and transcriptome analysis were performed with P. fragi. A total of 124 potential cold adaptation genes were identified, including nineteen unique cold-adaptive genes were detected in the genome of P. fragi D12. Three unique genes associated with pili protein were significantly upregulated at different degrees of low temperature, which may be the key to the strong low-temperature adaptability of P. fragi D12. Meanwhile, we were pleasantly surprised to find that Pseudomonas fragi D12 exhibited different cold-adaptation mechanisms under different temperature changes. When the temperature declined from 30°C to 15°C, the response included maintenance of the fluidity of cell membranes, increased production of extracellular polymers, elevation in the content of compatibility solutes, and reduction in the content of reactive oxygen species, thereby providing a stable metabolic environment. When the temperature decreased from 15°C to 4°C, the response mainly included increases in the expression of molecular chaperones and transcription factors, enabling the bacteria to restore normal transcription and translation. The response mechanism of P. fragi D12 to low-temperature exposure is discussed. The results provide new ideas for the cold-adaptation mechanism of cold-tolerant microorganisms.

Mariluticola halotolerans gen. nov., sp. nov., a novel member of the family Devosiaceae isolated from South China Sea sediment

A novel member of class Alphaproteobacteria was isolated from marine sediment of the South China Sea. Cells of strain LMO-2^T were Gram-stain negative, greyish in colour, motile, with a single lateral flagellum and short rod in shape with a slight curve. Strain LMO-2^T was positive for oxidase and negative for catalase. The bacterium grew aerobically at 10-40 °C (optimum, 25-30 °C), pH 5.5-10.0 (optimum, pH 7.0) and 0-9 % NaCl (w/v; optimum, 2-3 %). Phylogenetic analysis of the 16S rRNA gene sequence and phylogenomic analysis of the whole genome sequence indicated that strain LMO-2^T represents a new genus and a new species within the family Devosiaceae, class Alphaproteobacteria, phylum Pseudomonadota. Comparisons of the 16S rRNA gene sequences of strain LMO-2^T showed 94.8 % similarity to its closest relative. The genome size is ~3.45 Mbp with a DNA G+C content of 58.17 mol%. The strain possesses potential capability for the degradation of complex organic matter, i.e. fatty acid and benzoate. The predominant cellular fatty acids (>10 %) were C_16 : 0 and C_18 : 1 ω7c 11-methyl. The sole respiratory quinone was ubiquinone-10. The major identified polar lipids were diphosphatidylglycerol, phosphatidylglycerol and phospholipid. Based on the polyphasic taxonomic data, strain LMO-2^T represents a novel genus and a novel species for which the name Mariluticola halotolerans gen. nov., sp. nov., was proposed in the family Devosiaceae. The type strain is LMO-2^T (=CGMCC 1.19273^T=JCM 34934^T).

Description and genomic characterization of Gallaecimonas kandeliae sp. nov., isolated from the sediments of mangrove plant Kandelia obovate

The genus Gallaecimonas, proposed by Rodríguez-Blanco et al. (Int J Syst Evol Microbiol 60:504-509, 2010), is mainly isolated from marine environments. So far, only three species have been identified and characterized in this genus. In this study, a new Gallaecimonas strain named Q10^T was isolated from the sediments of mangrove plant Kandelia obovate taken from Dapeng district, Shenzhen, China. Strain Q10^T was a Gram-stain-negative, non-motile, strictly aerobic, rod-shaped bacterium, and grew with 0-8.0% (w/v) NaCl, at 10-45 °C and at pH 5.5-8.5. Phylogenetic analysis indicated that strain Q10^T and the three Gallaecimonas species formed a clade in the tree, with 16S rRNA gene sequence similarities ranging from 96.0 to 97.0%. The major respiratory quinone is Q8. The polar lipids comprised aminolipid, aminophospholipid, diphosphatidylglycerol, glycolipid, phosphatidylethanolamine, phosphatidylglycerol, glycophospholipid and phospholipid. The predominant fatty acids are C_16:0, C_17:1ω8c, summed feature 3 (C_16:1ω7c/C_16:1ω6c), and iso-C_16:0. The complete genome of strain Q10^T is 3,836,841 bp with a G+C content of 62.6 mol%. The orthologous proteins analysis revealed 55 unique proteins in strain Q10^T related to important biological processes, especially three frataxins related to iron-sulfur cluster assembly, which may play a pivotal role in environmental adaptability of this species. Based on polyphasic taxonomic data, strain Q10^T is considered to represent a novel species within the genus Gallaecimonas, for which the name Gallaecimonas kandelia sp. nov. is proposed. The type strain is Q10^T (=KCTC 92860^T=MCCC 1K08421^T). These results contribute to a better understanding of general features and taxonomy of the genus Gallaecimonas.

Genome Sequence and Evaluation of Safety and Probiotic Potential of Lactiplantibacillus plantarum LPJZ-658

This study aims to systematically evaluate the safety of a novel L. plantarum LPJZ-658 explored on whole-genome sequence analysis, safety, and probiotic properties assessment. Whole genome sequencing results demonstrated that L. plantarum LPJZ-658 consists of 3.26 Mbp with a GC content of 44.83%. A total of 3254 putative ORFs were identified. Of note, a putative bile saline hydrolase (BSH) (identity 70.4%) was found in its genome. In addition, the secondary metabolites were analyzed, and one secondary metabolite gene cluster was predicted to consist of 51 genes, which verified its safety and probiotic properties at the genome level. Additionally, L. plantarum LPJZ-658 exhibited non-toxic and non-hemolytic activity and was susceptible to various tested antibiotics, indicating that L. plantarum LPJZ-658 was safe for consumption. Moreover, the probiotic properties tests confirm that L. plantarum LPJZ-658 also exhibits tolerance to acid and bile salts, preferably hydrophobicity and auto-aggregation, and excellent antimicrobial activity against both Gram-positive and Gram-negative gastrointestinal pathogens. In conclusion, this study confirmed the safety and probiotic properties of L. plantarum LPJZ-658, suggesting it can be used as a potential probiotic candidate for human and animal applications.

Composition and fate of heat-resistant anaerobic spore-formers in the milk powder production line

Anaerobic spore-forming bacteria are a continuous threat to the dairy industry due to their ability to withstand processing conditions, such as those during heat treatment. These ubiquitous microorganisms have ample opportunity for multiple entry points into the milk chain, creating food quality and safety issues. Certain spore-formers, namely bacilli and clostridia, are more problematic due to their ability to spoil dairy products and pathogenicity. In this study, we investigated how milk treatment and milk powder production influenced the composition and survival of anaerobic spore-formers. Samples were obtained on three different days (replicate blocks) during the production of dairy powders and examined in a culture-dependent manner using the most probable number method coupled with 16S rRNA amplicon sequencing and metagenomic analysis of the enriched samples. Results revealed that the milk separation greatly affected the spore-former presence and composition which were detected along the entire production line from raw material to milk powders. Throughout the various points of the production line, the occurrence of species belonging to the Bacillus cereus sensu lato was higher than that of clostridia. Sequence variants (SVs) belonging to the anaerobic spore-forming genus Clostridium were taxonomically assigned to two SVs groups and were detected in all three replicate blocks. A total of 19 metagenome-assembled genomes were recovered from nine enrichments. Four near-complete and two medium-quality genomes were found in raw milk/milk powder samples and further assigned as Clostridium tyrobutyricum and Clostridium diolis, which may constitute a problem in the finished dairy product. In conclusion, our findings highlight spore-formers' importance on dairy quality and may aid in their intervention and control in the dairy production line.

Whole Genome Sequence Analysis of Cupriavidus necator C39, a Multiple Heavy Metal(loid) and Antibiotic Resistant Bacterium Isolated from a Gold/Copper Mine

Here a multiple heavy metal and antibiotic resistant bacterium Cupriavidus necator C39 (C. necator C39) was isolated from a Gold-Copper mine in Zijin, Fujian, China. C. necator C39 was able to tolerate intermediate concentrations of heavy metal(loid)s in Tris Minimal (TMM) Medium (Cu(II) 2 mM, Zn(II) 2 mM, Ni(II) 0.2 mM, Au(III) 70 μM and As(III) 2.5 mM). In addition, high resistance to multiple antibiotics was experimentally observed. Moreover, strain C39 was able to grow on TMM medium containing aromatic compounds such as benzoate, phenol, indole, p-hydroxybenzoic acid or phloroglucinol anhydrous as the sole carbon sources. The complete genome of this strain revealed 2 circular chromosomes and 1 plasmid, and showed the closest type strain is C. necator N-1^T based on Genome BLAST Distance Phylogeny. The arsenic-resistance (ars) cluster GST-arsR-arsICBR-yciI and a scattered gene encoding the putative arsenite efflux pump ArsB were identified on the genome of strain C39, which thereby may provide the bacterium a robust capability for arsenic resistance. Genes encoding multidrug resistance efflux pump may confer high antibiotic resistance to strain C39. Key genes encoding functions in degradation pathways of benzene compounds, including benzoate, phenol, benzamide, catechol, 3- or 4-fluorobenzoate, 3- or 4-hydroxybenzoate and 3,4-dihydroxybenzoate, indicated its potential for degrading those benzene compounds.

Isoptericola croceus sp. nov., a novel actinobacterium isolated from saline-alkali soil

A novel actinomycete, designated strain q2^T, was isolated from the saline-alkaline soil, collected from Daqing, Heilongjiang province, China. The results of phylogenetic analysis based on the 16S rRNA gene sequences indicated that strain q2^T belongs to the genus Isoptericola, and showed the highest sequence similarity to Isoptericola halotolerans KCTC 19046^T (98.48%) and Isoptericola chiayiensis KCTC 19740^T (98.13%), respectively. The average nucleotide identity values between strain q2^T and other members of the genus Isoptericola were lower than 95% recommended for distinguishing novel prokaryotic species. Cells of strain q2^T were Gram-staining-positive, aerobic, non-motile, rod-shaped and non-spore-forming. Colonies of strain q2^T were golden-yellow pigmented, tidy edged and smooth surfaced. Growth occurred at 15-37 °C (optimum, 29 °C), pH 7.0-10.0 (optimum, pH 8.0). The predominant respiratory quinones were MK-9(H₄) and MK-9(H₂). The main detected polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylinositol, and phosphatidylinositol mannoside. The peptidoglycan compositions were L-alanine, D-aspartic, L-glutamic acid and L-lysine (type A4α). The major cellular fatty acids (> 10%) were anteiso-C_15:0, iso-C_15:0, and anteiso-C_17:0. The G+C content of the genomic DNA was determined to be 69.7%. Based on the phenotypic, physiological, genotypic, and phylogenetic data, strain q2^T represents a novel species of the genus Isoptericola, for which the name Isoptericola croceus sp. nov. is proposed. The type strain is q2^T (= GDMCC 1.2923^T = KCTC 49759^T).

A genome catalog of the early-life human skin microbiome

Metagenome-assembled genomes have greatly expanded the reference genomes for skin microbiome. However, the current reference genomes are largely based on samples from adults in North America and lack representation from infants and individuals from other continents. Here we used ultra-deep shotgun metagenomic sequencing to profile the skin microbiota of 215 infants at age 2-3 months and 12 months who were part of the VITALITY trial in Australia as well as 67 maternally-matched samples. Based on the infant samples, we present the Early-Life Skin Genomes (ELSG) catalog, comprising 9,194 bacterial genomes from 1,029 species, 206 fungal genomes from 13 species, and 39 eukaryotic viral sequences. This genome catalog substantially expands the diversity of species previously known to comprise human skin microbiome and improves the classification rate of sequenced data by 25%. The protein catalog derived from these genomes provides insights into the functional elements such as defense mechanisms that distinguish early-life skin microbiome. We also found evidence for vertical transmission at the microbial community, individual skin bacterial species and strain levels between mothers and infants. Overall, the ELSG catalog uncovers the skin microbiome of a previously underrepresented age group and population and provides a comprehensive view of human skin microbiome diversity, function, and transmission in early life.

Polyphase taxonomy and genome analysis reveal the adaptability of Luteolibacter rhizosphaerae sp. nov. to the rhizosphere soil of Ulmus pumila L

A Gram-stain-negative, rod-shaped, non-flagellated, pale-yellow bacterium, designated GHJ8^T, was isolated from the rhizosphere soil of Ulmus pumila L., Shanxi Province, China. Growth occurred at 20-37 °C (optimum, 28 °C), pH 6.0-11.0 (optimum, pH 8.0), and 0-1% NaCl (optimum, 0%). Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain GHJ8^T was related to members of the genus Luteolibacter, and close to Luteolibacter flavescens GKX^T (98.5%), Luteolibacter luteus G-1-1-1^T (97.3%), Luteolibacter arcticus MC 3726^T (97.2%), and Luteolibacter marinus NBU1238^T (96.0%). The genome size of strain GHJ8^T was 6.2 Mbp, with a G + C content of 62.5%. Genomic mining revealed that the strain contained antibiotic resistance genes and secondary metabolic gene clusters, indicating that it had adaptation mechanisms to environmental stress. Comparative genomic analyses clearly separated strain GHJ8^T from the recognized species of the genus Luteolibacter based on average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values below the thresholds for species delineation. The major cellular fatty acids were iso-C_14:0 (30.8%), C_16:1 ω9c (23.0%), C_16:0 (17.3%), and C_14:0 (13.4%). The quinone system was composed of the major menaquinones MK-8, MK-9, and MK-10, and the principal polar lipids were diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, an unidentified aminophospholipid, an unidentified glycolipid, two unidentified phospholipids, and three unidentified lipids. Based on its phenotypic and genotypic properties and phylogenetic inference, strain GHJ8^T is a novel species of the genus Luteolibacter, for which the name Luteolibacter rhizosphaerae sp. nov. is proposed. The type strain is GHJ8^T (= GDMCC 1.2160^T = KCTC 82452^T = JCM 34400^T).

Genomic and in-vitro characteristics of a novel strain Lacticaseibacillus chiayiensis AACE3 isolated from fermented blueberry

Numerous different species of LAB are used in different fields due to their unique characteristics. However, Lacticaseibacillus chiayiensis, a newly established species in 2018, has limited microorganism resources, and lacks comprehensive evaluations of its properties. In this study, L. chiayiensis AACE3, isolated from fermented blueberry, was evaluated by genomic analysis and in vitro assays of the properties. The genome identified genes associated with biofilm formation (luxS, ccpA, brpA), resistance to oxidative stress (tpx, trxA, trxB, hslO), tolerance to acidic conditions (dltA, dltC), resistance to unfavorable osmotic pressure (opuBB, gbuA, gbuB, gbuC), and adhesion (luxS, dltA, dltC). The AACE3 showed 112 unique genes, relative to the other three L. chiayiensis strains. Among them, the presence of genes such as clpP, pepO, and feoA suggests a possible advantage of AACE3 over other L. chiayiensis in terms of environmental adaptation. In vitro evaluation of the properties revealed that AACE3 had robust antibacterial activity against eight common pathogens: Streptococcus agalactiae, Staphylococcus aureus, Escherichia coli, Salmonella enteritidis, Salmonella choleraesuis, Shigella flexneri, Pseudomonas aeruginosa, and Klebsiella pneumoniae. In addition, AACE3 showed more than 80% survival rate in all tests simulating gastrointestinal fluid, and it exhibited high antioxidant capacity. Interestingly, the cell culture supernatant was superior to intact organisms and ultrasonically crushed bacterial extracts in all tests of antioxidant capacity. These results suggested that the antioxidant capacity may originate from certain metabolites and extracellular enzymes produced by AACE3. Moreover, AACE3 was a moderate biofilm producer due to the self-agglomeration effect. Taken together, L. chiayiensis AACE3 appears to be a candidate strain for combating the growing incidence of pathogen infections and antioxidant production.

Sulfur metabolism in subtropical marine mangrove sediments fundamentally differs from other habitats as revealed by SMDB

Shotgun metagenome sequencing provides the opportunity to recover underexplored rare populations and identify difficult-to-elucidate biochemical pathways. However, information on sulfur genes, including their sequences, is scattered in public databases. Here, we introduce SMDB ( https://smdb.gxu.edu.cn/ )-a manually curated database of sulfur genes based on an in-depth review of the scientific literature and orthology database. The SMDB contained a total of 175 genes and covered 11 sulfur metabolism processes with 395,737 representative sequences affiliated with 110 phyla and 2340 genera of bacteria/archaea. The SMDB was applied to characterize the sulfur cycle from five habitats and compared the microbial diversity of mangrove sediments with that of other habitats. The structure and composition of microorganism communities and sulfur genes were significantly different among the five habitats. Our results show that microorganism alpha diversity in mangrove sediments was significantly higher than in other habitats. Genes involved in dissimilatory sulfate reduction were abundant in subtropical marine mangroves and deep-sea sediments. The neutral community model results showed that microbial dispersal was higher in the marine mangrove ecosystem than in others habitats. The Flavilitoribacter of sulfur-metabolizing microorganism becomes a reliable biomarker in the five habitats. SMDB will assist researchers to analyze genes of sulfur cycle from the metagenomic efficiently.

Comparative Genomic Analysis of Cold-Water Coral-Derived Sulfitobacter faviae: Insights into Their Habitat Adaptation and Metabolism

Sulfitobacter is one of the major sulfite-oxidizing alphaproteobacterial groups and is often associated with marine algae and corals. Their association with the eukaryotic host cell may have important ecological contexts due to their complex lifestyle and metabolism. However, the role of Sulfitobacter in cold-water corals remains largely unexplored. In this study, we explored the metabolism and mobile genetic elements (MGEs) in two closely related Sulfitobacter faviae strains isolated from cold-water black corals at a depth of ~1000 m by comparative genomic analysis. The two strains shared high sequence similarity in chromosomes, including two megaplasmids and two prophages, while both contained several distinct MGEs, including prophages and megaplasmids. Additionally, several toxin-antitoxin systems and other types of antiphage elements were also identified in both strains, potentially helping Sulfitobacter faviae overcome the threat of diverse lytic phages. Furthermore, the two strains shared similar secondary metabolite biosynthetic gene clusters and genes involved in dimethylsulfoniopropionate (DMSP) degradation pathways. Our results provide insight into the adaptive strategy of Sulfitobacter strains to thrive in ecological niches such as cold-water corals at the genomic level.

Functional insight into Cordyceps militaris sugar transporters by structure modeling, network analysis and allosteric regulation

Insights into the structures, functions and dynamics of Cordyceps militaris (C. militaris) sugar transporters are necessary for understanding their versatile metabolic capability for fungal growth. The sequence-function relationship study of 85 C. militaris sugar transporters showed that there is a gap between phylogenetic-based subfamily classification and their functions. Beyond protein sequences, structural modeling and principal component analysis of the structural ensemble revealed the different folds of the Car and Org subfamilies. Performing channel detection and network analysis found that the Alp and Hex subfamilies can be specifically distinguished from others by the betweenness of channel residues. Signature dynamics analysis further suggested that the Hex subfamily demonstrates different dynamics, with high flexibility at the H1 region in TM11. Furthermore, the H1 region as an allosteric site was examined by network parameter calculations that guided allosteric pathways between this region and the channel cavity. Together with gene expression data of C. militaris, e.g., Hex06741 in the Hex subfamily, it was promisingly expressed when sugar utilization was altered. This work demonstrates an in silico framework for investigating C. militaris sugar transporters as an example case study of the allosteric activity of the Hex subfamily and can facilitate sugar transporter engineering design that can further optimize the preferable sugar utilization and fermentation process of C. militaris.

Comparative genomics of Stutzerimonas balearica (Pseudomonas balearica): diversity, habitats, and biodegradation of aromatic compounds

Stutzerimonas balearica (Pseudomonas balearica) has been found principally in oil-polluted environments. The capability of S. balearica to thrive from the degradation of pollutant compounds makes it a species of interest for potential bioremediation applications. However, little has been reported about the diversity of S. balearica. In this study, genome sequences of S. balearica strains from different origins were analyzed, revealing that it is a diverse species with an open pan-genome that will continue revealing new genes and functionalities as the genomes of more strains are sequenced. The nucleotide signatures and intra- and inter-species variation of the 16S rRNA genes of S. balearica were reevaluated. A strategy of screening 16S rRNA gene sequences in public databases enabled the detection of 158 additional strains, of which only 23% were described as S. balearica. The species was detected from a wide range of environments, although mostly from aquatic and polluted environments, predominantly related to petroleum oil. Genomic and phenotypic analyses confirmed that S. balearica possesses varied inherent capabilities for aromatic compounds degradation. This study increases the knowledge of the biology and diversity of S. balearica and will serve as a basis for future work with the species.

Dissimilar gene repertoires of Dickeya solani involved in the colonization of lesions and roots of Solanum tuberosum

Dickeya and Pectobacterium species are necrotrophic pathogens that macerate stems (blackleg disease) and tubers (soft rot disease) of Solanum tuberosum. They proliferate by exploiting plant cell remains. They also colonize roots, even if no symptoms are observed. The genes involved in pre-symptomatic root colonization are poorly understood. Here, transposon-sequencing (Tn-seq) analysis of Dickeya solani living in macerated tissues revealed 126 genes important for competitive colonization of tuber lesions and 207 for stem lesions, including 96 genes common to both conditions. Common genes included acr genes involved in the detoxification of plant defense phytoalexins and kduD, kduI, eda (=kdgA), gudD, garK, garL, and garR genes involved in the assimilation of pectin and galactarate. In root colonization, Tn-seq highlighted 83 genes, all different from those in stem and tuber lesion conditions. They encode the exploitation of organic and mineral nutrients (dpp, ddp, dctA, and pst) including glucuronate (kdgK and yeiQ) and synthesis of metabolites: cellulose (celY and bcs), aryl polyene (ape), and oocydin (ooc). We constructed in-frame deletion mutants of bcsA, ddpA, apeH, and pstA genes. All mutants were virulent in stem infection assays, but they were impaired in the competitive colonization of roots. In addition, the ΔpstA mutant was impaired in its capacity to colonize progeny tubers. Overall, this work distinguished two metabolic networks supporting either an oligotrophic lifestyle on roots or a copiotrophic lifestyle in lesions. This work revealed novel traits and pathways important for understanding how the D. solani pathogen efficiently survives on roots, persists in the environment, and colonizes progeny tubers.

Whole-genome sequencing of a biocontrol Myxococcus xanthus R31 isolate and comparative genomic analysis

Tomato bacterial wilt (TBW) caused by Ralstonia solanacearum is one of the most destructive soil-borne diseases. Myxococcus xanthus R31, isolated from healthy tomato rhizosphere soil using the R. solanacearum baiting method, exhibiting good biocontrol efficacy against TBW. However, the genomic information and evolutionary features of R31 are largely unclear. Here, the high-quality genome assembly of R31 was presented. Using Nanopore sequencing technology, we assembled the 9.25 Mb complete genome of R31 and identified several extracellular enzyme proteins, including carbohydrate-active enzymes (CAZymes) and peptidases. We also performed a comparative genome analysis of R31 and 17 other strains of M. xanthus with genome sequences in the NCBI database to gain insights into myxobacteria predation and genome size expansion. Average nucleotide identity and digital DNA-DNA hybridization calculation and phylogenetic analysis indicated that R31 was closely related to the species M. xanthus. Further comparative genomics analysis suggested that, in addition to characteristics of predatory microorganisms, R31 contains many strain-specific genes, which may provide a genetic basis for its proficient predatory ability. This study provides new insights into R31 and other closely related species and facilitates studies using genetic approaches to further elucidate the predation mechanism of myxobacteria.

Assembly and analytical validation of a metagenomic reference catalog of human gut microbiota based on co-barcoding sequencing

Human gut microbiota is associated with human health and disease, and is known to have the second-largest genome in the human body. The microbiota genome is important for their functions and metabolites; however, accurate genomic access to the microbiota of the human gut is hindered due to the difficulty of cultivating and the shortcomings of sequencing technology. Therefore, we applied the stLFR library construction method to assemble the microbiota genomes and demonstrated that assembly property outperformed standard metagenome sequencing. Using the assembled genomes as references, SNP, INDEL, and HGT gene analyses were performed. The results demonstrated significant differences in the number of SNPs and INDELs among different individuals. The individual displayed a unique species variation spectrum, and the similarity of strains within individuals decreased over time. In addition, the coverage depth analysis of the stLFR method shows that a sequencing depth of 60X is sufficient for SNP calling. HGT analysis revealed that the genes involved in replication, recombination and repair, mobilome prophages, and transposons were the most transferred genes among different bacterial species in individuals. A preliminary framework for human gut microbiome studies was established using the stLFR library construction method.

Curtanaerobium respiraculi gen. nov., sp. nov., a novel anaerobic bacterium isolated from human bronchoalveolar lavage fluid

Two related anaerobic strains, designated as SWB101512^T and SWB19611, were isolated from the bronchoalveolar lavage fluid of two lung cancer patients. Cells were Gram-stain-positive, non-motile and non-spore-forming. Growth could be observed at 26-45 °C (optimum, 37 °C), pH 5.0-8.5 (optimum, pH 7.0) and with 0.5-2.0 % (v/w) NaCl (optimum, 1.0%). The 16S rRNA gene sequences of SWB101512^T and SWB19611 showed the highest similarities to Denitrobacterium detoxificans DSM 21843^T (91.1 and 91.3 %, respectively). The phylogenetic tree based on the 16S rRNA gene sequences and the core genome sequences demonstrated that the two strains clustered together and formed a distinct lineage within the family Eggerthellaceae. The DNA G+C contents of strains SWB101512^T and SWB19611 were 62.0 and 61.9 mol%, respectively. The predominant cellular fatty acids of strains SWB101512^T and SWB19611 were C_16 : 0 DMA (27.8 and 28.8 %, respectively). The respiratory menaquinone in both strains was menaquinone 6 and the polar lipid profile consisted of diphosphatidylglycerol, phosphatidylglycerol, two phospholipids, three glycolipids and three unidentified lipids. Based on evidence from phenotypic, chemotaxonomic and genomic analyses, a new genus and species belonging to the family Eggerthellaceae, named Curtanaerobium respiraculi gen. nov., sp. nov. is proposed. The type strain is SWB101512^T (=GDMCC 1.2991^T=JCM 35330^T).

A phylogenomic and comparative genomic analysis of Commensalibacter, a versatile insect symbiont

BACKGROUND

To understand mechanisms of adaptation and plasticity of pollinators and other insects a better understanding of diversity and function of their key symbionts is required. Commensalibacter is a genus of acetic acid bacterial symbionts in the gut of honey bees and other insect species, yet little information is available on the diversity and function of Commensalibacter bacteria. In the present study, whole-genome sequences of 12 Commensalibacter isolates from bumble bees, butterflies, Asian hornets and rowan berries were determined, and publicly available genome assemblies of 14 Commensalibacter strains were used in a phylogenomic and comparative genomic analysis.

RESULTS

The phylogenomic analysis revealed that the 26 Commensalibacter isolates represented four species, i.e. Commensalibacter intestini and three novel species for which we propose the names Commensalibacter melissae sp. nov., Commensalibacter communis sp. nov. and Commensalibacter papalotli sp. nov. Comparative genomic analysis revealed that the four Commensalibacter species had similar genetic pathways for central metabolism characterized by a complete tricarboxylic acid cycle and pentose phosphate pathway, but their genomes differed in size, G + C content, amino acid metabolism and carbohydrate-utilizing enzymes. The reduced genome size, the large number of species-specific gene clusters, and the small number of gene clusters shared between C. melissae and other Commensalibacter species suggested a unique evolutionary process in C. melissae, the Western honey bee symbiont.

CONCLUSION

The genus Commensalibacter is a widely distributed insect symbiont that consists of multiple species, each contributing in a species specific manner to the physiology of the holobiont host.

Identification, classification, and functional characterization of novel sponge-associated acidimicrobiial species

Sponges are known to harbour an exceptional diversity of uncultured microorganisms, including members of the phylum Actinobacteriota. While members of the actinobacteriotal class Actinomycetia have been studied intensively due to their potential for secondary metabolite production, the sister class of Acidimicrobiia is often more abundant in sponges. However, the taxonomy, functions, and ecological roles of sponge-associated Acidimicrobiia are largely unknown. Here, we reconstructed and characterized 22 metagenome-assembled genomes (MAGs) of Acidimicrobiia from three sponge species. These MAGs represented six novel species, belonging to five genera, four families, and two orders, which are all uncharacterized (except the order Acidimicrobiales) and for which we propose nomenclature. These six uncultured species have either only been found in sponges and/or corals and have varying degrees of specificity to their host species. Functional gene profiling indicated that these six species shared a similar potential to non-symbiotic Acidimicrobiia with respect to amino acid biosynthesis and utilization of sulfur compounds. However, sponge-associated Acidimicrobiia differed from their non-symbiotic counterparts by relying predominantly on organic rather than inorganic sources of energy, and their predicted capacity to synthesise bioactive compounds or their precursors implicated in host defence. Additionally, the species possess the genetic capacity to degrade aromatic compounds that are frequently found in sponges. The novel Acidimicrobiia may also potentially mediate host development by modulating Hedgehog signalling and by the production of serotonin, which can affect host body contractions and digestion. These results highlight unique genomic and metabolic features of six new acidimicrobiial species that potentially support a sponge-associated lifestyle.

Quantitative model for genome-wide cyclic AMP receptor protein binding site identification and characteristic analysis

Cyclic AMP receptor proteins (CRPs) are important transcription regulators in many species. The prediction of CRP-binding sites was mainly based on position-weighted matrixes (PWMs). Traditional prediction methods only considered known binding motifs, and their ability to discover inflexible binding patterns was limited. Thus, a novel CRP-binding site prediction model called CRPBSFinder was developed in this research, which combined the hidden Markov model, knowledge-based PWMs and structure-based binding affinity matrixes. We trained this model using validated CRP-binding data from Escherichia coli and evaluated it with computational and experimental methods. The result shows that the model not only can provide higher prediction performance than a classic method but also quantitatively indicates the binding affinity of transcription factor binding sites by prediction scores. The prediction result included not only the most knowns regulated genes but also 1089 novel CRP-regulated genes. The major regulatory roles of CRPs were divided into four classes: carbohydrate metabolism, organic acid metabolism, nitrogen compound metabolism and cellular transport. Several novel functions were also discovered, including heterocycle metabolic and response to stimulus. Based on the functional similarity of homologous CRPs, we applied the model to 35 other species. The prediction tool and the prediction results are online and are available at: https://awi.cuhk.edu.cn/∼CRPBSFinder.

Annelid adult cell type diversity and their pluripotent cellular origins

Annelids are a broadly distributed, highly diverse, economically and environmentally important group of animals. Most species can regenerate missing body parts, and many are able to reproduce asexually. Therefore, many annelids can generate all adult cell types in adult stages. However, the putative adult stem cell populations involved in these processes, as well as the diversity of adult cell types generated by them, are still unknown. Here, we recover 75,218 single cell transcriptomes of Pristina leidyi, a highly regenerative and asexually-reproducing freshwater annelid. We characterise all major annelid adult cell types, and validate many of our observations by HCR in situ hybridisation. Our results uncover complex patterns of regionally expressed genes in the annelid gut, as well as neuronal, muscle and epidermal specific genes. We also characterise annelid-specific cell types such as the chaetal sacs and globin+ cells, and novel cell types of enigmatic affinity, including a vigilin+ cell type, a lumbrokinase+ cell type, and a diverse set of metabolic cells. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. In these piwi+ cells, we also find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors. Finally, lineage reconstruction analyses reveal the existence of differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids for the first time and serve as a resource for studying annelid cell types and their evolution. On the other hand, our characterisation of a piwi+ cell population with a pluripotent stem cell signature will serve as a platform for the study of annelid stem cells and their role in regeneration.

The nickel-pincer coenzyme of lactate racemase: A case study of uncovering cofactor structure and biosynthesis

Cofactors are essential components of numerous enzymes, therefore their characterization by structural, biophysical, and biochemical approaches is crucial for understanding the resulting catalytic and regulatory mechanisms. In this chapter, we present a case study of a recently discovered cofactor, the nickel-pincer nucleotide (NPN), by demonstrating how we identified and thoroughly characterized this unprecedented nickel-containing coenzyme that is tethered to lactase racemase from Lactiplantibacillus plantarum. In addition, we describe how the NPN cofactor is biosynthesized by a panel of proteins encoded in the lar operon and describe the properties of these novel enzymes. Comprehensive protocols for conducting functional and mechanistic studies of NPN-containing lactate racemase (LarA) and the carboxylase/hydrolase (LarB), sulfur transferase (LarE), and metal insertase (LarC) used for NPN biosynthesis are provided for potential applications towards characterizing enzymes in the same or homologous families.

Pseudomonas aeruginosa responds to altered membrane phospholipid composition by adjusting the production of two-component systems, proteases and iron uptake proteins

Membrane protein and phospholipid (PL) composition changes in response to environmental cues and during infections. To achieve these, bacteria use adaptation mechanisms involving covalent modification and remodelling of the acyl chain length of PLs. However, little is known about bacterial pathways regulated by PLs. Here, we investigated proteomic changes in the biofilm of P. aeruginosa phospholipase mutant (∆plaF) with altered membrane PL composition. The results revealed profound alterations in the abundance of many biofilm-related two-component systems (TCSs), including accumulation of PprAB, a key regulator of the transition to biofilm. Furthermore, a unique phosphorylation pattern of transcriptional regulators, transporters and metabolic enzymes, as well as differential production of several proteases, in ∆plaF, indicate that PlaF-mediated virulence adaptation involves complex transcriptional and posttranscriptional response. Moreover, proteomics and biochemical assays revealed the depletion of pyoverdine-mediated iron uptake pathway proteins in ∆plaF, while proteins from alternative iron-uptake systems were accumulated. These suggest that PlaF may function as a switch between different iron-acquisition pathways. The observation that PL-acyl chain modifying and PL synthesis enzymes were overproduced in ∆plaF reveals the interconnection of degradation, synthesis and modification of PLs for proper membrane homeostasis. Although the precise mechanism by which PlaF simultaneously affects multiple pathways remains to be elucidated, we suggest that alteration of PL composition in ∆plaF plays a role for the global adaptive response in P. aeruginosa mediated by TCSs and proteases. Our study revealed the global regulation of virulence and biofilm by PlaF and suggests that targeting this enzyme may have therapeutic potential.

Genome-resolved metagenomics reveals depth-related patterns of microbial community structure and functions in a highly stratified, AMD overlaying mine tailings

Acid mine drainage (AMD) is a worldwide environmental problem, yet bioremediation is hampered by a limited knowledge of the reductive microbial processes in the AMD ecosystem. Here, we generate extensive metagenome and geochemical datasets to investigate how microbial populations and metabolic capacities driving major element cycles are structured in a highly stratified, AMD overlaying tailings environment. The results demonstrated an explicit depth-dependent differentiation of microbial community composition and function profiles between the surface and deeper tailings layers, paralleling the dramatic shifts in major physical and geochemical properties. Specifically, key genes involved in sulfur and iron oxidation were significantly enriched in the surface tailings, whereas those associated with reductive nitrogen, sulfur, and iron processes were enriched in the deeper layers. Genome-resolved metagenomics retrieved 406 intermediate or high-quality genomes spanning 26 phyla, including major new groups (e.g., Patescibacteria and DPANN). Metabolic models involving nitrogen, sulfur, iron, and carbon cycles were proposed based on the functional potentials of the abundant microbial genomes, emphasizing syntrophy and the importance of lesser-known taxa in the degradation of complex carbon compounds. These results have implications for in situ AMD bioremediation.

Microbial conversion of carbon dioxide and hydrogen into the fine chemicals hydroxyectoine and ectoine

This study explores a novel conversion of CO₂ into the chemicals hydroxyectoine and ectoine, which are compounds with high retail values in the pharmaceutical industry. Firstly, 11 species of microbes able to use CO₂ and H₂ and that have the genes for ectoines synthesis (ectABCD) were identified through literature search and genomic mining. Laboratory tests were then conducted to ascertain the capacity of these microbes to produce ectoines from CO₂. Results showed that the most promising bacteria for CO₂ to ectoines bioconversion areHydrogenovibrio marinus, Rhodococcus opacus, and Hydrogenibacillus schlegelii.Upon salinity and H₂/CO₂/O₂ ratio optimization,H. marinus accumulated 85 mg of ectoine g biomass^-1. Interestingly, R.opacusand H. schlegelii mainly produced hydroxyectoine (53 and 62 mg g biomass^-1), which has a higher commercial value. Overall, these results constitute the first proof of a novel valorization platform of CO₂ and lay the foundation for a new economic niche aimed at CO₂ recircularization into pharmaceuticals.

Ancylobacter mangrovi sp. nov., a novel endophytic bacterium isolated form mangrove plant☆

Two novel strains GSK1Z-4-2^T and MQZ15Z-1 were isolated from branches of mangrove plants collected from Guangxi Zhuang Autonomous Region, China. Both strains were Gram-negative, aerobic, non-flagellated and non-spore-forming bacteria. The comparison of 16S rRNA gene sequences initially indicated that the two strains were assigned to the genus Ancylobacter with sharing the highest similarity to Ancylobacter pratisalsi DSM 102029^T (97.3%). The 16S rRNA gene sequence similarity, average nucleotide identity (ANI) and in silico DNA-DNA hybridization (isDDH) values between strains GSK1Z-4-2^T and MQZ15Z-1 were 99.9%, 97.4% and 77.4%, respectively, which revealed that the two strains belonged to the same species. Phylogenetic analyses based on 16S rRNA gene sequences and the core proteome showed that the two strains formed a well-supported cluster with A. pratisalsi DSM 102029^T. Moreover, the ANI and isDDH values between strain GSK1Z-4-2^T and A. pratisalsi DSM 102029^T were 83.0% and 25.8%, respectively, demonstrating that strain GSK1Z-4-2^T was a previously undescribed species. Meanwhile, strains GSK1Z-4-2^T and MQZ15Z-1 exhibited most of chemotaxonomic and phenotypic features consistent with the description of the genus Ancylobacter. Based on the polyphasic data, strains GSK1Z-4-2^T and MQZ15Z-1 should represent a novel species of the genus Ancylobacter, for which the name Ancylobacter mangrovi sp. nov. is proposed. The type strain is GSK1Z-4-2^T (=MCCC 1K07181^T = JCM 34924^T).

Reduction, evolutionary pattern and positive selection of genes encoding formate dehydrogenase in Wood-Ljungdahl pathway of gastrointestinal acetogens suggests their adaptation to formate-rich habitats

Acetogens are anaerobes using Wood-Ljungdahl pathway (WLP) as the terminal electron acceptor for both assimilation and dissimilation of CO₂ and widely distributed in diverse habitats. However, their habitat adaptation is often unclear. Given that bacterial genome evolution is often the result of environmental selective pressure, hereby we analysed gene copy number, phylogeny and selective pressure of genes involved in WLP within known genomes of 43 species to study the habitat adaption of gastrointestinal acetogens. The gene copy number of formate dehydrogenase (FDH) in gastrointestinal acetogens was much lower than that of non-gastrointestinal acetogens, and in five cases, no FDH genes were found in the genomes of five gastrointestinal acetogens, but that of the other WLP genes showed no difference. The evolutionary pattern of FDH genes was significantly different from that of the other enzymes. Additionally, seven positively selected sites were only identified in the fdhF genes, which means fdhF mutations favoured their adaptation. Collectively, reduction or loss of FDH genes and their evolutionary pattern as well as positive selection in gastrointestinal acetogens indicated their adaptation to formate-rich habitats, implying that FDH genes catalysing CO₂ reduction to formate as the first step of methyl branch of WLP may have evolved independently.

Description and genomic characterization of Nocardioides bruguierae sp. nov., isolated from Bruguiera gymnorhiza

Strains BSK12Z-3^T and BSK12Z-4, two Gram-stain-positive, aerobic, non-spore-forming strains, were isolated from Shankou Mangrove Nature Reserve, Guangxi Zhuang Autonomous Region, China. The diagnostic diamino acid in the cell-wall peptidoglycan of strain BSK12Z-3^T was LL-diaminopimelic acid and MK-8(H₄) was the predominant menaquinone. The polar lipids comprised diphosphatidylglycerol (DPG), phosphatidylglycerol (PG) and phospholipid (PL). The major fatty acids was iso-C_16:0. Phylogenetic analysis based on 16S rRNA gene sequences suggested that the two strains fell within the genus Nocardioides, appearing most closely related to Nocardioides ginkgobilobae KCTC 39594^T (97.5-97.6 % sequence similarity) and Nocardioides marinus DSM 18248^T (97.4-97.6 %). Genome-based phylogenetic analysis confirmed that strains BSK12Z-3^T and BSK12Z-4 formed a distinct phylogenetic cluster within the genus Nocardioides. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values of strains BSK12Z-3^T, BSK12Z-4 with their most related species N. marinus DSM18248^T were within the ranges of 77.2-77.3 % and 21.3-21.4 %, respectively, clearly indicated that strains BSK12Z-3^T, BSK12Z-4 represented novel species. Strains BSK12Z-3^T and BSK12Z-4 exhibited 99.9 % 16S rRNA gene sequence similarity. The ANI and dDDH values between the two strains were 97.8 % and 81.1 %, respectively, suggesting that they belong to the same species. However, DNA fingerprinting discriminated that they were not from one clonal origin. Based on phylogenomic and phylogenetic analyses coupled with phenotypic and chemotaxonomic characterizatons, strains BSK12Z-3^T and BSK12Z-4 could be classified as a novel species of the genus Nocardioides, for which the name Nocardioides bruguierae sp. nov., is proposed. The type strain is BSK12Z-3^T (=CGMCC 4.7709^T = JCM 34554^T).

Functional characterization and taxonomic classification of novel gammaproteobacterial diversity in sponges

Sponges harbour exceptionally diverse microbial communities, whose members are largely uncultured. The class Gammaproteobacteria often dominates the microbial communities of various sponge species, but most of its diversity remains functional and taxonomically uncharacterised. Here we reconstructed and characterised 32 metagenome-assembled genomes (MAGs) derived from three sponge species. These MAGs represent ten novel species and belong to seven orders, of which one is new. We propose nomenclature for all these taxa. These new species comprise sponge-specific bacteria with varying levels of host specificity. Functional gene profiling highlights significant differences in metabolic capabilities across the ten species, though each also often exhibited a large degree of metabolic diversity involving various nitrogen- and sulfur-based compounds. The genomic features of the ten species suggest they have evolved to form symbiotic interaction with their hosts or are well-adapted to survive within the sponge environment. These Gammaproteobacteria are proposed to scavenge substrates from the host environment, including metabolites or cellular components of the sponge. Their diverse metabolic capabilities may allow for efficient cycling of organic matter in the sponge environment, potentially to the benefit of the host and other symbionts.

Characterization and Genome Analysis of Cladobotryum mycophilum, the Causal Agent of Cobweb Disease of Morchella sextelata in China

Cobweb disease is a fungal disease that can cause serious damage to edible mushrooms worldwide. To investigate cobweb disease in Morchella sextelata in Guizhou Province, China, we isolated and purified the pathogen responsible for the disease. Through morphological and molecular identification and pathogenicity testing on infected M. sextelata, we identified Cladobotryum mycophilum as the cause of cobweb disease in this region. This is the first known occurrence of this pathogen causing cobweb disease in M. sextelata anywhere in the world. We then obtained the genome of C. mycophilum BJWN07 using the HiFi sequencing platform, resulting in a high-quality genome assembly with a size of 38.56 Mb, 10 contigs, and a GC content of 47.84%. We annotated 8428 protein-coding genes in the genome, including many secreted proteins, host interaction-related genes, and carbohydrate-active enzymes (CAZymes) related to the pathogenesis of the disease. Our findings shed new light on the pathogenesis of C. mycophilum and provide a theoretical basis for developing potential prevention and control strategies for cobweb disease.

Comparative physiological and transcriptome analysis in cultivated and wild sugarcane species in response to hydrogen peroxide-induced oxidative stress

BACKGROUND

Sugarcane is an important energy crop grown worldwide,supplementing various renewable energy sources. Cultivated and wild sugarcane species respond differently to biotic and abiotic stresses. Generally, wild species are tolerant to various abiotic stresses. In the present study, the physiological and molecular responses of cultivated and wild sugarcane species to oxidative stress at the transcriptional levels were compared. Transcriptional responses were determined using RNAseq. The representative RNA-seq transcript values were validated by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and confirmed through physiological responses.

RESULTS

Oxidative stress causes leaf-rolling and -tip drying in cultivated sugarcane, but the wild species are tolerant. Higher chlorophyll fluorescence was observed in the wild species than that in the cultivated varieties under stress. Wild species can maintain a higher chlorophyll stability index than the cultivated species, which was confirmed by the lower transcripts of the chlorophyllase gene in the wild species than that in the cultivated variety. Transcription factor genes (NAC, MYB, and WRKY) were markedly expressed in response to oxidative stress, revealing their involvement in stress tolerance. The analysis revealed synchronized expression of acetyl-transferase, histone2A, cellulose synthase, and secondary cell wall biosynthetic genes in the wild species. The validation of selected genes and 15 NAC transcription factors using RT-qPCR revealed that their expression profiles were strongly correlated with RNA-seq. To the best of our knowledge, this is the first report on the oxidative stress response in cultivated and wild sugarcane species.

CONCLUSION

Physiological and biochemical changes in response to oxidative stress markedly differ between cultivated and wild sugarcane species. The differentially expressed stress-responsive genes are grouped intothe response to oxidative stress, heme-binding, peroxidase activity, and metal ion binding categories. Chlorophyll maintenance is a stress tolerance response enhanced by the differential regulation of the chlorophyllase gene.There is a considerable difference in the chlorophyll stability index between wild and cultivated varieties. We observed a substantial regulation of secondary wall biosynthesis genes in the wild species compared with that in the cultivated variety, suggesting differences in stress tolerance mechanisms.

Whole-Genome Sequence of Lactiplantibacillus plantarum Mut-3, Isolated from Indonesian Fermented Soybean (Tempeh)

Lactiplantibacillus plantarum Mut-3 was isolated from tempeh. After whole-genome sequencing, analysis of its possibility as a probiotic candidate was performed using subsystem analysis with RAST with the SEED viewer.