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

Mycolyltransferase is important for biofilm formation and pathogenesis of Tsukamurella keratitis

Tsukamurella, a group of multi-drug resistant, Gram-positive, aerobic, and partially acid-fast bacteria, are emerging causes of bacterial conjunctivitis and keratitis. However, the pathogenesis of Tsukamurella keratitis is largely unknown. To address this, we used New Zealand White rabbits to develop the first eye infection model and conducted in vitro tests to study the pathogenesis mechanisms of Tsukamurella. There is increasing evidence that biofilms play a significant role in ocular infections, leading us to hypothesize that biofilm formation is crucial for effective Tsukamurella infection. In order to look for potential candidate genes which are important in biofilm formation and Tsukamurella keratitis. We performed genome sequencing of two ocular isolates, T. pulmonis-PW1004 and T. tyrosinosolvens-PW899, to identify potential virulence factors. Through in vitro and in vivo studies, we characterized their biological roles in mediating Tsukamurella keratitis. Our findings confirmed that Tsukamurella is an ocular pathogen by fulfilling Koch's postulates, and using genome sequence data, we identified tmytC, encoding a mycolyltransferase, as a crucial gene in biofilm formation and causing Tsukamurella keratitis in the rabbit model. This is the first report demonstrating the novel role of mycolyltransferase in causing ocular infections. Overall, our findings contribute to a better understanding of Tsukamurella pathogenesis and provide a potential target for treatment. Specific inhibitors targeting TmytC could serve as an effective treatment option for Tsukamurella infections.

Genomic analysis of Rhodopirellula sp. P2 reveals its role in fucoidan degradation

Fucoidan, the main polysaccharide in various species of brown seaweed, has a high annual production. It is an important source of marine organic carbon and exhibits diverse biological activities and significant application potential. Rhodopirellula sp. P2, a novel marine bacterium of the phylum Planctomycetota, was isolated from intertidal algae samples collected from the Weihai coast, the Yellow Sea, China. The strain P2 is a Gram-negative, aerobic, and pear-shaped bacterium. Here, we report the complete genome sequence of Rhodopirellula sp. P2. The genome of strain P2 consists of a single circular chromosome with 7,291,416 bp and a GC content of 57.38 %, including 5462 protein-coding genes, 2 rRNA genes, and 48 tRNA genes. Genomic analysis revealed that strain P2 possessed 173 CAZymes and 106 sulfatases, indicating that strain P2 has the potential ability to utilize multiple polysaccharides, especially hydrolyze fucoidan to fucose. The genome of strain P2 also encodes a gene cluster related to bacterial microcompartment, suggesting the ability of strain P2 to metabolize fucose. These results enhance the understanding of the diversity and ecological functions of Planctomycetota, and also facilitate the exploitation of Planctomycetota and enzyme resources to utilize fucoidan. This study provides genetic insights into fucoidan catabolism by Planctomycetota, expanding our understanding of fucoidan-degrading microbial groups.

Expression and characterization of monofunctional alcohol dehydrogenase enzymes in Clostridium thermocellum

Clostridium thermocellum is a thermophilic anaerobic bacterium that could be used for cellulosic biofuel production due to its strong native ability to consume cellulose, however its ethanol production ability needs to be improved to enable commercial application. In our previous strain engineering work, we observed a spontaneous mutation in the native adhE gene that reduced ethanol production. Here we attempted to complement this mutation by heterologous expression of 18 different alcohol dehydrogenase (adh) genes. We were able to express all of them successfully in C. thermocellum. Surprisingly, however, none of them increased ethanol production, and several actually decreased it. Our findings contribute to understanding the correlation between C. thermocellum ethanol production and Adh enzyme cofactor preferences. The identification of a set of adh genes that can be successfully expressed in this organism provides a foundation for future investigations into how the properties of Adh enzymes affect ethanol production.

The complete genome sequence of proteases-producing Shewanella sp. H8 isolated from Antarctica

Bacteria of the genus Shewanella in the class Gammaproteobacteria are widely distributed in marine environments. Shewanella sp. H8, was isolated from a red algae sample collected from Nelson Island, Antarctica. Here, we present the complete genome sequence of strain H8, which consists of a single circular chromosome comprising 4,490,743 nucleotides with 40.59 % G + C content without any plasmid. In total, 3983 protein coding genes, 95 tRNA genes, and 25 rRNA genes were obtained. Genomic analysis of strain H8 showed that it contains four cold shock proteins and three fatty acid desaturases and possesses the potential to synthesize hglE-KS, arylpolyene, betalactone and RiPP-like compounds. Through genomic annotation, 91 protease-encoding genes were identified within the genome of strain H8. These proteases are classified into six categories based on their catalytic types. Among these proteases, metalloproteinases and serine proteases are dominant. These proteases may provide carbon and nitrogen sources to H8 by degrading proteins in the environment. This study will provide potential genetic information for the future research and development of cold-adapted proteases.

Comparative genome analysis of endophytic Bacillus amyloliquefaciens MR4: a potential biocontrol agent isolated from wild medicinal plant root tissue

In this investigation, 396 endophytic bacterial strains from six indigenous medicinal plant species within the Xinjiang Tumor Peak National Nature Reserve were subjected to screening. The strain MR4 emerged as a noteworthy contender, demonstrating pronounced biocontrol capabilities coupled with exceptional cold tolerance. Through morphological scrutiny and comprehensive genomic sequencing, MR4 was identified as Bacillus amyloliquefaciens. Antagonistic assays revealed MR4's efficacy in suppressing the causative agents of cotton wilt and verticillium wilt, achieving inhibition rates surpassing 50%. Analyses, underpinned by PCR methodologies, indicated MR4's capacity to biosynthesize a minimum of eight distinct antimicrobial agents. The whole-genome sequencing data indicated that B. amyloliquefaciens MR4 had the genome size and GC content of 4,017,872 bp and 47.14%, respectively, and 4191 coding genes were identified. The genome consists of a single chromosome and one plasmid. Moreover, it was augmented by annotations from various databases, including GO, KEGG, and COG. The pathogenicity of MR4 undergoes evaluation, while predictions concerning MR4's secondary metabolites have disclosed gene clusters for 13 varieties of these compounds, with particular emphasis on surfactins and fengycin. Comparative analyses with four paradigmatic strains shed light on MR4's genomic composition and its phylogenetic lineage within the Bacillus genus. The genomic data pertaining to MR4 have been duly submitted to the NCBI GenBank, bearing the accession numbers CP146236 (Chr1) and CP146237 (plas1). This study endeavors to furnish potent microbial resources for the biocontrol and enhancement of plant growth, thereby providing a theoretical groundwork for MR4's agronomic utilization.

Metagenomic analysis unveiled the response of microbial community and antimicrobial resistome in natural water body to duck farm sewage

Sewages from duck farms are often recognized as a major source of antimicrobial resistance and pathogenic bacteria discharged to natural water bodies, but few studies depicted the dynamic changes in resistome and microbial communities in the rivers under immense exposure of sewage discharge. In this study, we investigated the ecological and environmental risks of duck sewages to the rivers that geographically near to the duck farms with short-distance (<1 km) using 16S rRNA amplicon and metagenomic sequencing. The results showed that a total of 20 ARG types were identified with abundances ranged from 0.61 to 1.33 cpc. Of note, the genes modulate resistances against aminoglycoside, bacitracin and beta-lactam were the most abundant ARGs. Limnohabitans, Fluviibacter and Cyanobium were the top 3 predominant genera in the microbial community. The alpha diversity of overall microbial community decrease while the abundance of pathogen increase during the input of sewage within 200 m. Sul1 and bacA were the dominant ARGs brought from duck farm sewage. The community variations of ARGs and microbiome were primarily driven by pH and temperature. Total phosphorus was significantly correlated to alpha diversity and top 30 ARGs subtype. Stochastic processes was the dominated microbial assembly pattern and did not be altered by sewage. We also highlighted the ecological risk caused by blaGES which possibly could be mitigated by Cyanobacteria, and the natural water body can purify partial ARGs as well as microbiome from duck farms sewage. These findings expanded our knowledge regarding the ecological risks by wastes from the livestock farm, and underscoring the necessity to monitor ARGs in farm-surrounding water bodies.

An Optimized Method for Reconstruction of Transcriptional Regulatory Networks in Bacteria Using ChIP-exo and RNA-seq Datasets

Transcriptional regulatory networks (TRNs) in bacteria are crucial for elucidating the mechanisms that regulate gene expression and cellular responses to environmental stimuli. These networks delineate the interactions between transcription factors (TFs) and their target genes, thereby uncovering the regulatory processes that modulate gene expression under varying environmental conditions. Analyzing TRNs offers valuable insights into bacterial adaptation, stress responses, and metabolic optimization from an evolutionary standpoint. Additionally, understanding TRNs can drive the development of novel antimicrobial therapies and the engineering of microbial strains for biofuel and bioproduct production. This protocol integrates advanced data analysis pipelines, including ChEAP, DEOCSU, and DESeq2, to analyze omics datasets that encompass genome-wide TF binding sites and transcriptome profiles derived from ChIP-exo and RNA-seq experiments. This approach minimizes both the time required and the risk of bias, making it accessible to non-expert users. Key steps in the protocol include preprocessing and peak calling from ChIP-exo data, differential expression analysis of RNA-seq data, and motif and regulon analysis. This method offers a comprehensive and efficient framework for TRN reconstruction across various bacterial strains, enhancing both the accuracy and reliability of the analysis while providing valuable insights for basic and applied research.

Genetic potential for aerobic respiration and denitrification in globally distributed respiratory endosymbionts

The endosymbiont Candidatus Azoamicus ciliaticola was proposed to generate ATP for its eukaryotic host, an anaerobic ciliate of the Plagiopylea class, fulfilling a function analogous to mitochondria in other eukaryotic cells. The discovery of this respiratory endosymbiosis has major implications for both evolutionary history and ecology of microbial eukaryotes. However, with only a single species described, knowledge of its environmental distribution and diversity is limited. Here we report four complete, circular metagenome assembled genomes (cMAGs) representing respiratory endosymbionts inhabiting groundwater in California, Ohio, and Germany. These cMAGs form two lineages comprising a monophyletic clade within the uncharacterized gammaproteobacterial order UBA6186, enabling evolutionary analysis of their key protein complexes. Strikingly, all four cMAGs encode a cytochrome cbb3 oxidase, which indicates that these endosymbionts have the capacity for aerobic respiration. Accordingly, we detect these respiratory endosymbionts in diverse habitats worldwide, thus further expanding the ecological scope of this respiratory symbiosis.

Screening, identification, engineering, and characterization of Bacillus-derived α-amylase for effective tobacco starch degradation

In this study, two high-performing α-amylase-producing strains, CK3-5 and A8-1 were successfully isolated and characterized, which were taxonomically confirmed as Bacillus velezensis through whole-genome sequencing and bioinformatics. Bioinformatic sequence analysis and molecular docking revealed the catalytic triad (Asp173-Glu208-Asp274) essential for α-amylase function. Through metabolic engineering, the recombinant strain BAX-5/PT17amy(A8-1)SP002 was developed, which exhibited the highest α-amylase activity of 1440 U/mL upon fermentation optimization, marking a 9.2-fold enhancement over the wild-type strain A8-1, and it successfully degraded 6 % of the starch in the tobacco leaves within 48 h, while the content of 13 harmful substances, including acetamide, pyridine, and acetonitrile, was reduced by 8.6 % to 25.2 %. This study reveals a novel α-amylase gene from B. velezensis and establishes an efficient expression system in B. amyloliquefaciens, offering valuable insights for industrial α-amylase production.

High lead-tolerant mutant Bacillus tropicus AT31-1 from rhizosphere soil of Pu-erh and its remediation mechanism

In this study, we successfully generated the mutant strain Bacillus tropicus AT31-1 from AT31 through atmospheric room-temperature plasma mutagenesis. This mutant strain AT31-1 demonstrated an impressive 48.6 % removal efficiency in 400 mg/L lead medium. Comparative genomic analysis showed that the mutant strain AT31-1 had three mutation sites, which affect the efflux RND transporter permease subunit, the response regulator transcription factor, and a gene with unknown function. The transcriptional analysis showed a notable upregulation in the expression of 283 genes in AT31-1 as lead concentrations increased from 0 to 200 mg/L and then to 400 mg/L, which include zinc-transporting ATPase, ferrous iron transport protein B, NADH dehydrogenase, and others. The Gene ontology function of the peptide metabolic process, along with the KEGG pathway of carbon metabolism were identified as closely linked to the extreme lead tolerance of AT31-1. This study presents novel insights into the lead tolerance mechanisms of bacteria.

Exploring protein natural diversity in environmental microbiomes with DeepMetagenome

Protein natural diversity offers a vast sequence space for protein engineering, and deep learning enables its detection from metagenomes/proteomes without prior assumptions. DeepMetagenome, a Python-based method, explores protein diversity through modules for training and analyzing sequence datasets. The deep learning model includes Embedding, Conv1D, LSTM, and Dense layers, with sequence feature analysis for data cleaning. Applied to metallothioneins from a database of over 146 million coding features, DeepMetagenome identified over 500 high-confidence metallothionein sequences, outperforming DIAMOND and CNN-based models. It showed stable performance compared to a Transformer-based model over 25 epochs. Among 23 synthesized sequences, 20 exhibited metal resistance. The tool also successfully explored the diversity of three additional protein families and is freely available on GitHub with detailed instructions.

COG database update 2024

The Clusters of Orthologous Genes (COG) database, originally created in 1997, has been updated to reflect the constantly growing collection of completely sequenced prokaryotic genomes. This update increased the genome coverage from 1309 to 2296 species, including 2103 bacteria and 193 archaea, in most cases, with a single representative genome per genus. This set covers all genera of bacteria and archaea that included organisms with 'complete genomes' as per NCBI databases in November 2023. The number of COGs has been expanded from 4877 to 4981, primarily by including protein families involved in bacterial protein secretion. Accordingly, COG pathways and functional groups now include secretion systems of types II through X, as well as Flp/Tad and type IV pili. These groupings allow straightforward identification and examination of the prokaryotic lineages that encompass-or lack-a particular secretion system. Other developments include improved annotations for the rRNA and tRNA modification proteins, multi-domain signal transduction proteins, and some previously uncharacterized protein families. The new version of COGs is available at https://www.ncbi.nlm.nih.gov/research/COG, as well as on the NCBI FTP site https://ftp.ncbi.nlm.nih.gov/pub/COG/, which also provides archived data from previous COG releases.

Genome sequence of Ehrlichia muris from Ixodes ricinus collected in Italy on a migratory bird provides epidemiological and evolutionary insights

Ticks are prominent vectors of several zoonotic diseases. Tick-borne pathogens include the members of the genus Ehrlichia, which are obligate intracellular bacteria infecting immune and hematopoietic cells. Ehrlichia muris predominantly affects rodents, but was also reported to be a human pathogen. The known geographical distribution of this bacterium ranges from Asia, to the USA and eastern Europe. In the present work, we report the finding of E. muris in an Ixodes ricinus tick collected from a migratory bird (Turdus iliacus) in Italy, southern Europe. We sequenced the total DNA from this tick sample, and, thanks to a dedicated bioinformatic pipeline, selectively assembled the genome of the bacterium, which represents the first one for E. muris from Europe. Phylogenetic and comparative genomic analyses were then performed. Accounting for tick species distribution, bird migratory routes, and molecular phylogeny of the bacterium, it is likely that this bird transported the tick to Italy from an endemic area of E. muris, such as eastern Europe. In addition, comparative genomic analyses highlighted that E. muris and other Ehrlichia spp. display copy number variations in two families of membrane proteins, likely due to recent gene duplication, deletion and recombination events. These differences are probably a source of variability for surface antigens to evade host immunity, with a potential role in host adaptation and specificity. The present results underline the impact of migratory birds on the spread of tick-borne pathogens towards non-endemic areas, highlighting the need for further epidemiological surveillance at bird ringing stations in Italy, and advocating further investigations on possible local transmission of E. muris in competent mammalian hosts.

Transcriptomic and proteomic changes associated with cobalamin-dependent propionate production by the rumen bacterium Xylanibacter ruminicola

Xylanibacter ruminicola is an abundant rumen bacterium that produces propionate in a cobalamin (vitamin B12)-dependent manner via the succinate pathway. However, the extent to which this occurs across ruminal Xylanibacter and closely related bacteria, and the effect of cobalamin supplementation on the expression of propionate pathway genes and enzymes has yet to be investigated. To assess this, we screened 14 strains and found that almost all strains produced propionate when supplemented with cobalamin. X. ruminicola KHP1 was selected for further study, including complete genome sequencing, and comparative transcriptomics and proteomics of KHP1 cultures grown with and without supplemented cobalamin. The complete KHP1 genome was searched for cobalamin-binding riboswitches and four were predicted, though none were closely located to any of the succinate pathway genes, which were dispersed at numerous genomic loci. Cobalamin supplementation led to the differential expression of 17.5% of genes, including genes encoding the cobalamin-dependent methylmalonyl-CoA mutase and some methylmalonyl-CoA decarboxylase subunits, but most propionate biosynthesis pathway genes were not differentially expressed. The effect of cobalamin supplementation on the KHP1 proteome was much less pronounced, with the only differentially abundant propionate pathway enzyme being methylmalonyl-CoA mutase, which had greater abundance when supplemented with cobalamin. Our results demonstrate that cobalamin supplementation does not result in induction of the entire propionate biosynthesis pathway, but consistently increased expression of methylmalonyl-CoA mutase at transcriptome and proteome levels. The magnitude of the differential expression of propionate pathway genes observed was minor compared to that of genes proximate to predicted cobalamin riboswitches.

IMPORTANCE

In ruminants, the rumen microbial community plays a critical role in nutrition through the fermentation of feed to provide vital energy substrates for the host animal. Propionate is a major rumen fermentation end-product and increasing its production is desirable given its importance in host glucose production and impact on greenhouse gas production. Vitamin B12 (cobalamin) can induce propionate production in the prominent rumen bacterium Xylanibacter ruminicola, but it is not fully understood how cobalamin regulates propionate pathway activity. Contrary to expectation, we found that cobalamin supplementation had little effect on propionate pathway expression at transcriptome and proteome levels, with minor upregulation of genes encoding the cobalamin-dependent enzyme of the pathway. These findings provide new insights into factors that regulate propionate production and suggest that cobalamin-dependent propionate production by X. ruminicola is controlled post-translationally.

Pseudogymnoascus destructans invasion stage impacts the skin microbial functions of highly vulnerable Myotis lucifugus

The role of the skin microbiome in resistance and susceptibility of wildlife to fungal pathogens has been examined from a taxonomic perspective but skin microbial function, in the context of fungal infection, has yet to be studied. Our objective was to understand effects of a bat fungal pathogen site infection status and course of invasion on skin microbial function. We sampled seven hibernating colonies of Myotis lucifugus covering three-time points over the course of Pseudogymnoascus destructans (Pd) invasion and white nose syndrome (pre-invasion, epidemic, and established). Our results support three new hypotheses about Pd and skin functional microbiome: (1) there is an important effect of Pd invasion stage, especially at the epidemic stage; (2) disruption by the fungus at the epidemic stage could decrease anti-fungal functions with potential negative effects on the microbiome and bat health; (3) the collection site might have a larger influence on microbiomes at the pre-invasion stage rather than at epidemic and established stages. Future studies with larger sample sizes and using meta-omics approaches will help confirm these hypotheses, and determine the influence of the microbiome on wildlife survival to fungal disease.

Comprehensive screening of industrially relevant components at genome scale using a high-quality gene overexpression collection of Corynebacterium glutamicum

Development of efficient microbial strains for biomanufacturing requires deep understanding of the biology and functional components responsible for the synthesis, transport, and tolerance of the target compounds. A high-quality controllable gene overexpression strain collection was constructed for the industrial workhorse Corynebacterium glutamicum covering 99.7% of its genes. The collection was then used for comprehensive screening of components relevant to biomanufacturing features. In total, 15 components endowing cells with improved hyperosmotic tolerance and l-lysine productivity were identified, including novel transcriptional factors and DNA repair proteins. Systematic interrogation of a subset of the collection revealed efficient and specific exporters functioning in both C. glutamicum and Escherichia coli. Application of the new exporters was showcased to construct a strain with the highest l-threonine production level reported for C. glutamicum (75.1 g/l and 1.5 g/l·h) thus far. The genome-scale gene overexpression collection will serve as a valuable resource for fundamental biological studies and for developing industrial microorganisms for producing amino acids and other biochemicals.

MVP: a modular viromics pipeline to identify, filter, cluster, annotate, and bin viruses from metagenomes

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.IMPORTANCEThe significance of our work lies in the development of Modular Viromics Pipeline (MVP), an integrated and user-friendly pipeline tailored exclusively for viromics analyses. MVP stands out due to its modular design, which ensures easy installation, execution, and integration of new tools and databases. By combining state-of-the-art tools such as geNomad and CheckV, MVP provides high-quality viral genome recovery and taxonomy and host assignment, and functional annotation, addressing the limitations of existing pipelines. MVP's ability to handle diverse sample types, including environmental, human microbiome, and plant-associated samples, makes it a versatile tool for the broader microbiome research community. By standardizing the analysis process and providing easily interpretable results, MVP enables researchers to perform comprehensive studies of viral communities, significantly advancing our understanding of viral ecology and its impact on various ecosystems.

Total biosynthesis of the medicinal triterpenoid saponin astragalosides

Astragalus membranaceus has been used in traditional Chinese medicine for over 2,000 years. Its major active triterpenoid saponins, astragalosides, have attracted great attention due to their multiple health benefits and applications in medicine. Despite this, the biosynthetic machinery for astragalosides remains enigmatic. Here a chromosome-level genome assembly of A. membranaceus was generated. The identification of two tailoring enzymes required for astragaloside biosynthesis enabled the discovery of a triterpenoid biosynthetic gene cluster, leading to elucidation of the complete astragaloside biosynthetic pathway. This pathway is characterized by a sequence of selective hydroxylation, epoxidation and glycosylation reactions, which are mediated by three cytochrome P450s, one 2-oxoglutarate-dependent dioxygenase and two glycosyltransferases. Reconstitution of this biosynthetic machinery in Nicotiana benthamiana allowed for heterologous production of astragaloside IV. These findings build a solid foundation for addressing the sourcing issues associated with astragalosides and broaden our understanding of the diversity of terpene biosynthetic gene clusters.

Metagenomics unravel distinct taxonomic and functional diversities between terrestrial and aquatic biomes

Microbes in terrestrial and aquatic ecosystems play crucial roles in driving ecosystem functions, but currently, there is a lack of comparison regarding their taxonomic and functional diversities. Here, we conducted a global analysis to investigate the disparities in microbial taxonomy and microbial-mediated biogeochemical cycles between terrestrial and aquatic ecosystems. Results showed a higher relative abundance of bacteria, especially Actinobacteria and Acidobacteria, in soil than water metagenomes, leading to a greater proportion of genes related to membrane transport, regulatory, and cellular signaling. Moreover, there was a higher abundance of genes associated with carbohydrate, sulfur, and potassium metabolisms in the soil, while those involved in nitrogen and iron metabolisms were more prevalent in the water. Thus, both soil and water microbiomes exhibited unique taxonomic and functional properties associated with biogeochemical processes, providing valuable insights into predicting and understanding the adaptation of microbes in different ecosystems in the face of climate change.

Succeed to culture a novel lineage symbiotic bacterium of Mollicutes which widely found in arthropods intestine uncovers the potential double-edged sword ecological function

Symbiotic gut bacteria play crucial role in host health. Symbionts are widely distributed in arthropod intestines, but their ecological functions are poorly understood due to the inability to cultivate them. Members of Candidatus Bacilliplasma (CB) are widely distributed in crustacean intestine and maybe commensals with hosts, but the paucity of pure cultures has limited further insights into their physiologies and functions. Here, four strains of representative CB bacteria in shrimp intestine were successfully isolated and identified as members of a novel Order in the Phylum Mycoplasmatota. Through genome assembly, the circular genome maps of the four strains were obtained, and the number of coding genes ranged from 1,886 to 1,980. Genomic analysis suggested that the bacteria were missing genes for many critical pathways including the TCA cycle and biosynthesis pathways for amino acids and coenzyme factors. The analysis of 16S amplification data showed that Shewanella, Pseudomonas and CB were the dominant at the genera level in the intestine of Penaeus vannamei. Ecological functional experiments revealed that the strains were symbionts and colonized shrimp intestines. Our valued findings can greatly enhance our understanding and provides new insights into the potentially significant role of uncultured symbiotic bacteria in modulating host health.

Dairy farm waste: A potential reservoir of diverse antibiotic resistance and virulence genes in aminoglycoside- and beta-lactam-resistant Escherichia coli in Gansu Province, China

Aminoglycosides (AGs) and beta-lactams are the most commonly used antimicrobials in animal settings, particularly on dairy farms. Dairy farm waste is an important reservoir of antibiotic resistance genes (ARGs) and virulence genes (VGs) in environmental Escherichia coli, which is an important indicator of environmental contamination and foodborne pathogen that potentially threaten human and animal health. In the present study, we aimed to characterize the ARGs and VGs in AG- and beta-lactam-resistant E. coli from dairy farm waste in Gansu Province, China. The dairy farm waste consisted of fecal (n = 265) and sewage (n = 54) samples processed using standard microbiological techniques and the Clinical & Laboratory Standards Institute guidelines. The total DNA of AG- and beta-lactam-resistant E. coli was extracted, and whole-genome sequencing (WGS) was performed using the Illumina NovaSeq platform and analyzed using various bioinformatics tools. In this study, among 84.3% (269/319) of the E. coli strains, 23.8% (64/269) were identified as AG- and beta-lactam-resistant E. coli. WGS analysis revealed a large pool of ARGs belonging to multiple classes such as AGs, beta-lactams, aminocoumarins, fluoroquinolones, macrolides, phenicol, tetracyclines, phosphonic acid, disinfecting and antiseptic agents, elfamycin, rifamycin, and multidrug resistance genes. Furthermore, virulome analysis of 64 E. coli strains revealed clinically important virulence factors associated with adherence, biofilm, invasion, auto-transportation, siderophores, secretion systems, toxins, anti-phagocytosis, quorum sensing, regulation, metabolism, and motility. We identified dairy farm feces and sewage waste as important reservoirs of antimicrobial resistance and virulence determinants in E. coli in Gansu, China, which can threaten human and animal health through ecological exposure and contamination of food and water. We recommend continuous large-scale surveillance in dairy farm settings to formulate protective guidelines for public health safety.

The marine environmental microbiome mediates physiological outcomes in host nematodes

BACKGROUND

Nematodes are the most abundant metazoans in marine sediments, many of which are bacterivores; however, how habitat bacteria affect physiological outcomes in marine nematodes remains largely unknown.  RESULTS: Here, we used a Litoditis marina inbred line to assess how native bacteria modulate host nematode physiology. We characterized seasonal dynamic bacterial compositions in L. marina habitats and examined the impacts of 448 habitat bacteria isolates on L. marina development, then focused on HQbiome with 73 native bacteria, of which we generated 72 whole genomes sequences. Unexpectedly, we found that the effects of marine native bacteria on the development of L. marina and its terrestrial relative Caenorhabditis elegans were significantly positively correlated. Next, we reconstructed bacterial metabolic networks and identified several bacterial metabolic pathways positively correlated with L. marina development (e.g., ubiquinol and heme b biosynthesis), while pyridoxal 5'-phosphate biosynthesis pathway was negatively associated. Through single metabolite supplementation, we verified CoQ10, heme b, acetyl-CoA, and acetaldehyde promoted L. marina development, while vitamin B6 attenuated growth. Notably, we found that only four development correlated metabolic pathways were shared between L. marina and C. elegans. Furthermore, we identified two bacterial metabolic pathways correlated with L. marina lifespan, while a distinct one in C. elegans. Strikingly, we found that glycerol supplementation significantly extended L. marina but not C. elegans longevity. Moreover, we comparatively demonstrated the distinct gut microbiota characteristics and their effects on L. marina and C. elegans physiology.

CONCLUSIONS

Given that both bacteria and marine nematodes are dominant taxa in sedimentary ecosystems, the resource presented here will provide novel insights to identify mechanisms underpinning how habitat bacteria affect nematode biology in a more natural context. Our integrative approach will provide a microbe-nematodes framework for microbiome mediated effects on host animal fitness.

Hohaiivirga grylli gen. nov., sp. nov., a New Member of the Family Methylobacteriaceae, Isolated from Cricket (Gryllus chinensis)

A Gram-staining negative, non-motile, rod-shaped, oxidase negative and catalase positive strain WL0021T was isolated from cricket (Gryllus chinensis) living in the campus of Hohai University. Strain WL0021T was characterized utilizing a polyphasic taxonomy approach. The major fatty acids (> 5%) for strain WL0021T were C16:0 and summed feature 8, and the major polar lipids were diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, phosphatidylcholine, phospholipid, two aminolipids, and an unidentified polar lipid. Ubiquinone-10 was detected as the predominant respiratory quinone. The results of 16S rRNA gene phylogenetic analyses revealed that strain WL0021T had the highest sequence similarity of 95.3% to Microvirga flavescens c27j1T and strain WL0021T formed a distinct linage within the family Methylobacteriaceae in the phylogenetic trees. Whole genomic DNA G+C content was 48.3%. Combined with the results from this study, strain WL0021T should represent a novel genus in the family Methylobacteriaceae, for which the name Hohaiivirga grylli gen. nov., sp. nov. (type strain WL0021T=GDMCC 1.2420T =JCM 34655T=MCCC 1K05886T) is proposed.

Adlercreutzia wanghongyangiae sp. nov., and Adlercreutzia shanghongiae sp. nov., two new members of the genus Adlercreutzia isolated from plateau pika (Ochotona curzoniae)

Four anaerobic, Gram-stain-positive, non-motile, non-sporulating rod-shaped bacterial strains (R7T, R21, R22 and R25T) were isolated from the intestinal contents of plateau pika (Ochotona curzoniae) collected from the Qinghai-Tibet Plateau, PR China. The four isolates grew at between 25 and 42 °C (optimally at 35-37 °C), and with 0.3-3.3% NaCl (w/v) [optimum, 1.3% (w/v)]. Adding l-arginine to the medium could promote their growth. Strains R7T and R21 were most closely related to Adlercreutzia caecimuris B7T (97.48% 16S rRNA gene sequence similarity). Strains R25T and R22 were most closely related to Adlercreutzia equolifaciens DSM 19450T (98.25% 16S rRNA gene sequence similarity). The genome sequences of R7T and R25T were 2.89 and 2.90 Mb in size with 63.6 and 62.8 mol% DNA G+C contents, respectively. Phylogenetic analysis based on 16S rRNA gene sequences and core genes revealed that R7T and R21 were most closely related to A. caecimuris B7T and Adlercreutzia mucosicola DSM 19490T, whereas R25T and R22 were most closely related to A. equolifaciens DSM 19450T and Adlercreutzia rubneri ResAG-91T. R7T, R25T and the closely related species had average nucleotide identity (ANI) values of 81.9-83.2% as well as digital DNA-DNA hybridisation (dDDH) values between 27.3 and 27.9%, which clearly indicated that they represent two novel species within the genus Adlercreutzia. For R7T and R25T, meso-diaminopimelic acid was the diagnostic diamino acid in the cell-wall peptidoglycan, and the whole cell sugars included galactose, glucose and ribose. On the basis of these results, we propose that strains R7T and R25T represent two novel species of the genus Adlercreutzia, namely Adlercreutzia wanghongyangiae sp. nov. and Adlercreutzia shanghongiae sp. nov., respectively. The type strains are R7T (=GDMCC 1.4459T=KCTC 25860T) and R25T (=GDMCC 1.4458T=KCTC 25861T).

Volatile Organic Compound Metabolism on Early Earth

Biogenic volatile organic compounds (VOCs) constitute a significant portion of gas-phase metabolites in modern ecosystems and have unique roles in moderating atmospheric oxidative capacity, solar radiation balance, and aerosol formation. It has been theorized that VOCs may account for observed geological and evolutionary phenomena during the Archaean, but the direct contribution of biology to early non-methane VOC cycling remains unexplored. Here, we provide an assessment of all potential VOCs metabolized by the last universal common ancestor (LUCA). We identify enzyme functions linked to LUCA orthologous protein groups across eight literature sources and estimate the volatility of all associated substrates to identify ancient volatile metabolites. We hone in on volatile metabolites with confirmed modern emissions that exist in conserved metabolic pathways and produce a curated list of the most likely LUCA VOCs. We introduce volatile organic metabolites associated with early life and discuss their potential influence on early carbon cycling and atmospheric chemistry.

The Car1 Knockout Mice Exhibit Antidepressant-like Behaviors Accompanied with Gut Microbiota Disturbance

Major depressive disorder (MDD) is a severe mental disorder with largely unknown mechanisms. Carbonic anhydrases convert CO2 to carbonates and protons, playing roles in various brain functions. Carbonic anhydrase 1 (Car1) is particularly abundant and may be linked to microbiota at interstitial sites. We developed Car1-deficient mice to explore the relationship between depression-like behaviors and gut microbiota. Behavioral tests confirmed depression-like behavior in Car1-/- mice. Fecal samples from Car1-/- and WT mice were collected, and 16S rRNA gene sequencing identified distinct microbiota components between the groups. Car1-/- mice exhibited significantly increased immobility in the tail suspension test (TST) compared to WT mice. The gut microbiota composition differed at the phylum level in p_Bacteroidetes, p_Verrucomicrobia, p_Firmicutes, and p_Tenericutes. At the family level, Car1-/- mice had significantly different abundances in eight microbiota groups compared to WT mice. Car1 deficiency is associated with depressive-like behavior and gut microbiota dysbiosis, potentially linked to depressive-like phenotypes.

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Genome-wide survey of the bipartite structure and pathogenesis-related genes of Neostagonosporella sichuanensis, a causal agent of Fishscale bamboo rhombic-spot disease

Bamboo resources have garnered significant global attention due to their excellent capacity for regeneration and high yield. Rhombic-spot disease, a substantial threat to fishscale bamboo (Phyllostachys heteroclada), is primarily caused by Neostagonosporella sichuanensis. This study first reported the genome assemblies and characteristics of two N. sichuanensis isolates using PacBio and Illumina sequencing platforms. The genomes of N. sichuanensis strain SICAUCC 16-0001 and strain SICAUCC 23-0140, with sizes of 48.0 Mb and 48.4 Mb, respectively, revealed 10,289 and 10,313 protein-coding genes. Additionally, they contained 34.99 and 34.46% repetitive sequences within AT-rich regions, with notable repeat-induced point mutation activity. Comparative genome analysis identified 1,049 contracted and 45 expanded gene families in the genome of N. sichuanensis, including several related to pathogenicity. Several gene families involved in mycotoxin metabolism, secondary metabolism, sterol biosynthesis and transport, and cell wall degradation were contracted. Compared to most analyzed necrotrophic, hemibiotrophic, and phaeosphaeriacous pathogens, the genomes of two N. sichuanensis isolates exhibited fewer secondary metabolite enzymes, carbohydrate-active enzymes, plant cell wall degrading enzymes, secreted proteins, and effectors. Comparative genomics analysis suggested that N. sichuanensis shares more similar characteristics with hemibiotrophic pathogens. Based on single carbon source tests, N. sichuanensis strains demonstrated a higher potential for xylan decomposition than pectin and cellulose. The proportion of cell wall-degrading enzyme effectors occupied a high proportion of the total effectors of the N. sichuanensis genomes. These findings provide valuable insights into uncovering the pathogenesis of N. sichuanensis toward the efficient management of rhombic-spot disease of fishscale bamboo.

Annotation of Functions of Sequences of Concern and Its Relevance to the New Biosecurity Regulatory Framework in the United States

Introduction

Recent regulations from United States Government agencies reshape the screening of synthetic nucleic acids. These take a step away from categorizing hazard on the basis of "bad" taxa and invoke the function of the sequence in pathogenesis or intoxication. Ascertaining functions related to pathogenesis and distinguishing these from other molecular abilities that are unproblematic is not simple. Some have suggested that this information can be readily obtained from existing databases of pathogens.

Objectives

We evaluate how virulence factors are described in current databases of pathogens and their adequacy for biothreat data science. We discuss limitations of how virulence factors have been conceived and propose using the sequence of concern (SoC) term to distinguish sequences with biothreat from those without. We discuss ways in which databases of SoCs might be implemented for research and regulatory purposes. We describe ongoing work improving functional descriptions of SoCs.

Methods

We assess the adequacy of descriptions of virulence factors in pathogen databases following extensive engagement with the literature in microbial pathogenesis.

Results/Conclusions

Descriptions of virulence factors in pathogen databases are inadequate for understanding biothreats. Many are not biothreats and would not be concerning if transferred to another pathogen. New gene ontology terms have been authored, and those specific to pathogenic viral processes are being generalized to make them relevant to other pathogenic taxa. This allows better understanding by humans and better recognition by machines. A database of annotated functions of SoCs could benefit the evolving biosecurity regulatory framework in the United States.

Metagenomic functional profiling: to sketch or not to sketch?

MOTIVATION

Functional profiling of metagenomic samples is essential to decipher the functional capabilities of microbial communities. Traditional and more widely used functional profilers in the context of metagenomics rely on aligning reads against a known reference database. However, aligning sequencing reads against a large and fast-growing database is computationally expensive. In general, k-mer-based sketching techniques have been successfully used in metagenomics to address this bottleneck, notably in taxonomic profiling. In this work, we describe leveraging FracMinHash (implemented in sourmash, a publicly available software), a k-mer-sketching algorithm, to obtain functional profiles of metagenome samples.

RESULTS

We show how pieces of the sourmash software (and the resulting FracMinHash sketches) can be put together in a pipeline to functionally profile a metagenomic sample. We named our pipeline fmh-funprofiler. We report that the functional profiles obtained using this pipeline demonstrate comparable completeness and better purity compared to the profiles obtained using other alignment-based methods when applied to simulated metagenomic data. We also report that fmh-funprofiler is 39-99× faster in wall-clock time, and consumes up to 40-55× less memory. Coupled with the KEGG database, this method not only replicates fundamental biological insights but also highlights novel signals from the Human Microbiome Project datasets.

AVAILABILITY AND IMPLEMENTATION

This fast and lightweight metagenomic functional profiler is freely available and can be accessed here: https://github.com/KoslickiLab/fmh-funprofiler. All scripts of the analyses we present in this manuscript can be found on GitHub.

Integrated analyses of the intestinal microbiome and transcriptome in Ningxiang piglets

Ningxiang (NX) pig has been recognized as one of the most famous Chinese indigenous breeds due to its characteristics in stress resistance. However, intestinal microbial feature and gene profiling in NX piglets have not been studied. Here, we compared the intestinal microbiome and transcriptome between NX and Duroc × Landrace × Large white (DLY) piglets and found the high enrichment of several colonic Bacteroides, Prevotella and Clostridium species in NX piglets. Further functional analyses revealed their predominant function in methane, glycolysis and gluconeogenesis metabolism. Our mRNA-sequencing data unraveled the distinct colonic gene expression between these two breeds. In particular, we showed that the improved intestinal function in NX piglets may be determined by enhanced intestinal barrier gene expression and varied immune gene expression through modulating the composition of the gut microbes. Together, our study revealed the intestinal characteristics of NX piglets, providing their potential application in improving breeding strategies and developing dietary interventions.

Characterization of gut symbionts from wild-caught Drosophila and other Diptera: description of Utexia brackfieldae gen. nov., sp. nov., Orbus sturtevantii sp. nov., Orbus wheelerorum sp. nov, and Orbus mooreae sp. nov

Non-culture based surveys show that the bacterial family Orbaceae is widespread in guts of insects, including wild Drosophila. Relatively few isolates have been described, and none has been described from Drosophila. We present the isolation and characterization of five strains of Orbaceae from wild-caught flies of the genera Drosophila (Diptera: Drosophilidae) and Neogriphoneura (Diptera: Lauxaniidae). Cells are generally rod-shaped, mesophilic, and measure 0.8-2.0 µm long by 0.3-0.5 µm wide. Optimal growth was observed under ambient atmosphere. Reconstruction of phylogenies from the 16S rRNA gene and from single-copy orthologs verify placement of these strains within Orbaceae. Cells exhibited similar fatty acid profiles to those of other Orbaceae. Strain lpD01T shared 74% average nucleotide identity (ANI) with its closest relatives Ca. Schmidhempelia bombi Bimp and Zophobihabitans entericus IPMB12T. Results from multiple genome-wide similarity comparisons indicate lpD01T should be classified as a novel species within a novel genus. The major respiratory quinone for lpD01T is ubiquinone Q-8. lpD02T, lpD03, lpD04T, and BiBT are more closely related to Orbus hercynius CN3T (76, 77, 76, and 77% ANI, respectively) than to other described Orbaceae. Genomic and phylogenetic analyses suggest that lpD03 and lpD04T belong to the same species and that lpD02T, lpD03/lpD04T, and BiBT are each novel species of the genus Orbus. The proposed names of these strains are Utexia brackfieldae gen. nov., sp. nov. (type strain lpD01T =NCIMB 15517T =ATCC TSD-399T), Orbus sturtevantii sp. nov (type strain lpD02T =NCIMB 15518T =ATCC TSD-400T), Orbus wheelerorum sp. nov. (type strain lpD04T =NCIMB 15520T =ATCC TSD-401T), and Orbus mooreae sp. nov (type strain BiBT=NCIMB 15516T =ATCC TSD-402T). The isolation and characterization of these strains expands the repertoire of culturable bacteria naturally associated with insects, including the model organism D. melanogaster.

Comparative genomics reveals ample evidence to Ganoderma sinense cultivars for molecular identification and new FIP exploration

The first dikaryotic genome of Ganoderma cultivar Zizhi S2 (56.76 Mb, 16,681 genes) has been sequenced recently. 98.15% of complete BUSCOs were recovered in this genome assembly and high-confidence annotation rate improved to 91.41%. Collinearity analysis displayed the nuclear genome were 80.2% and 93.84% similar to reference genome of G. sinense at nucleotide and amino acid levels, which presented 8,521 core genes and 880 unique orthologous gene groups. Among that, at least six functional genes (tef1-α, β-tubulin, rpb2, CaM, Mn-SOD and VeA) and a newly discovered fip gene were highly similar 99.27% ∼100% to those in reference genome. And the mt-LSU, mt-SSU and 13 PCGs in their mitogenome were also highly conserved with 99.27%-99.87% and 99.08%-100% identity, respectively. So that, this cultivar Zizhi S2 is confirmed conspecific with Ganoderma sinense (NCBI: txid1077348). The new fip gene (MN635280.1_336bp) existing a novel mutation which can be reflected on the phylogenetic tree and 3-dimensional model topology structure.

Spatially structured exchange of metabolites enhances bacterial survival and resilience in biofilms

Biofilm formation enhances bacterial survival and antibiotic tolerance, but the underlying mechanisms are incompletely understood. Here, we show that biofilm growth is accompanied by a reduction in bacterial energy metabolism and membrane potential, together with metabolic exchanges between the inner and outer regions in biofilms. More specifically, nutrient-starved cells in the interior supply amino acids to cells in the periphery, while peripheral cells experience a decrease in membrane potential and provide fatty acids to interior cells. Fatty acids facilitate the repair of starvation-induced membrane damage in inner cells and enhance their survival in the presence of antibiotics. Thus, metabolic exchanges between inner and outer cells contribute to survival of the nutrient-starved inner cells and contribute to antibiotic tolerance within the biofilm.

Polysaccharides induce deep-sea Lentisphaerae strains to release chronic bacteriophages

Viruses are ubiquitous in nature and play key roles in various ecosystems. Notably, some viruses (e.g. bacteriophage) exhibit alternative life cycles, such as chronic infections without cell lysis. However, the impact of chronic infections and their interactions with the host organisms remains largely unknown. Here, we found for the first time that polysaccharides induced the production of multiple temperate phages infecting two deep-sea Lentisphaerae strains (WC36 and zth2). Through physiological assays, genomic analysis, and transcriptomics assays, we found these bacteriophages were released via a chronic style without host cell lysis, which might reprogram host polysaccharide metabolism through the potential auxiliary metabolic genes. The findings presented here, together with recent discoveries made on the reprogramming of host energy-generating metabolisms by chronic bacteriophages, shed light on the poorly explored marine virus-host interaction and bring us closer to understanding the potential role of chronic viruses in marine ecosystems.

Staphylococcal exoribonuclease YhaM destabilizes ribosomes by targeting the mRNA of a hibernation factor

The hibernation-promoting factor (Hpf) in Staphylococcus aureus binds to 70S ribosomes and induces the formation of the 100S complex (70S dimer), leading to translational avoidance and occlusion of ribosomes from RNase R-mediated degradation. Here, we show that the 3'-5' exoribonuclease YhaM plays a previously unrecognized role in modulating ribosome stability. Unlike RNase R, which directly degrades the 16S rRNA of ribosomes in S. aureus cells lacking Hpf, YhaM destabilizes ribosomes by indirectly degrading the 3'-hpf mRNA that carries an intrinsic terminator. YhaM adopts an active hexameric assembly and robustly cleaves ssRNA in a manganese-dependent manner. In vivo, YhaM appears to be a low-processive enzyme, trimming the hpf mRNA by only 1 nucleotide. Deletion of yhaM delays cell growth. These findings substantiate the physiological significance of this cryptic enzyme and the protective role of Hpf in ribosome integrity, providing a mechanistic understanding of bacterial ribosome turnover.

CemR atypical response regulator impacts energy conversion in Campylobacteria

Campylobacter jejuni and Arcobacter butzleri are microaerobic food-borne human gastrointestinal pathogens that mainly cause diarrheal disease. These related species of the Campylobacteria class face variable atmospheric environments during infection and transmission, ranging from nearly anaerobic to aerobic conditions. Consequently, their lifestyles require that both pathogens need to adjust their metabolism and respiration to the changing oxygen concentrations of the colonization sites. Our transcriptomic and proteomic studies revealed that C. jejuni and A. butzleri, lacking a Campylobacteria-specific regulatory protein, C. jejuni Cj1608, or a homolog, A. butzleri Abu0127, are unable to reprogram tricarboxylic acid cycle or respiration pathways, respectively, to produce ATP efficiently and, in consequence, adjust growth to changing oxygen supply. We propose that these Campylobacteria energy and metabolism regulators (CemRs) are long-sought transcription factors controlling the metabolic shift related to oxygen availability, essential for these bacteria's survival and adaptation to the niches they inhabit. Besides their significant universal role in Campylobacteria, CemRs, as pleiotropic regulators, control the transcription of many genes, often specific to the species, under microaerophilic conditions and in response to oxidative stress.

IMPORTANCE

C. jejuni and A. butzleri are closely related pathogens that infect the human gastrointestinal tract. In order to infect humans successfully, they need to change their metabolism as nutrient and respiratory conditions change. A regulator called CemR has been identified, which helps them adapt their metabolism to changing conditions, particularly oxygen availability in the gastrointestinal tract so that they can produce enough energy for survival and spread. Without CemR, these bacteria, as well as a related species, Helicobacter pylori, produce less energy, grow more slowly, or, in the case of C. jejuni, do not grow at all. Furthermore, CemR is a global regulator that controls the synthesis of many genes in each species, potentially allowing them to adapt to their ecological niches as well as establish infection. Therefore, the identification of CemR opens new possibilities for studying the pathogenicity of C. jejuni and A. butzleri.

Differential Strategies of Ectomycorrhizal Development between Suillus luteus and Pinus massoniana in Response to Nutrient Changes

Ectomycorrhizal fungi employ different strategies for mycelial growth and host colonization under varying nutrient conditions. However, key genes associated with mycorrhizal interaction should be influenced solely by the inoculation treatment and not by nutrient variations. To utilize subtle nutrient differences and rapidly screen for key genes related to the interaction between Suillus luteus and Pinus massoniana, we performed an inoculation experiment using culture bottles containing high- and low-nutrient media. Interestingly, S. luteus LS88 promoted the growth of P. massoniana seedlings without mature ectomycorrhiza, and the impact of LS88 inoculation on P. massoniana roots was greater than that of nutrient changes. In this study, the resequenced genome of the LS88 strain was utilized for transcriptome analysis of the strain. The analysis indicated that a unique gene encoding glutathione S-transferase (GST) in LS88 is likely involved in colonizing P. massoniana roots. In this study, the GST gene expression was independent of nutrient levels. It was probably induced by P. massoniana and could be used as a marker for S. luteus colonization degree.

Transcriptome responses of Lactococcus paracarnosus to different gas compositions and co-culture with Brochothrix thermosphacta

Lactococcus (Lc.) paracarnosus and the phylogenetically closely related Lc. carnosus species are common members of the microbiota in meat stored under modified atmosphere and at low temperature. The effect of these strains on meat spoilage is controversially discussed. While some strains are known to cause spoilage, others are being studied for their potential to suppress the growth of spoilage and pathogenic bacteria. In this study, Lc. paracarnosus DSM 111017T was selected based on a previous study for its ability to suppress the growth of meat spoilers, including Brochothrix thermosphacta. The mechanism by which this bioprotective strain inhibits competing bacteria and how it contributes to spoilage are not yet known. To answer these two questions, we investigated the effect of four different headspace gas mixtures (simulated air (21 % O2/79 % N2); HiOx-MAP (70 % O2/30 % CO2); nonOx-MAP (70 % N2/ 30 % CO2); simulated vacuum (100 % N2) and the presence of Brochothrix (B.) thermosphacta TMW 2.2101 on the growth and transcriptional response of Lc. paracarnosus DSM 111017T when cultured on a meat simulation agar surface at 4 °C. Analysis of genes specifically upregulated by the gas mixtures used revealed metabolic pathways that may lead to different levels of spoilage metabolites production. We propose that under elevated oxygen levels, Lc. paracarnosus preferentially converts pyruvate from glucose and glycerol to uncharged acetoin/diacetyl instead of lactate to counteract acid stress. Due to the potential production of a buttery off-flavour, the strain may not be suitable as a protective culture in meat packaged under high‑oxygen conditions. 70 % N2/ 30 % CO2, simulated vacuum- and the presence of Lc. paracarnosus inhibited the growth of B. thermosphacta TMW 2.2101. However, B. thermosphacta did not affect gene regulation of metabolic pathways in Lc. paracarnosus, and genes previously predicted to be involved in B. thermosphacta growth suppression were not regulated at the transcriptional level. In conclusion, the study indicates that the gas mixture used in packaging significantly affects the metabolism and spoilage potential of Lc. paracarnosus and its ability to inhibit B. thermosphacta growth.

RNA-Seq reveals that Pseudomonas aeruginosa mounts growth medium-dependent competitive responses when sensing diffusible cues from Burkholderia cenocepacia

Most habitats host diverse bacterial communities, offering opportunities for inter-species interactions. While competition might often dominate such interactions, little is known about whether bacteria can sense competitors and mount adequate responses. The competition sensing hypothesis proposes that bacteria can use cues such as nutrient stress and cell damage to prepare for battle. Here, we tested this hypothesis by measuring transcriptome changes in Pseudomonas aeruginosa exposed to the supernatant of its competitor Burkholderia cenocepacia. We found that P. aeruginosa exhibited significant growth-medium-dependent transcriptome changes in response to competition. In an iron-rich medium, P. aeruginosa upregulated genes encoding the type-VI secretion system and the siderophore pyoverdine, whereas genes encoding phenazine toxins and hydrogen cyanide were upregulated under iron-limited conditions. Moreover, general stress response and quorum sensing regulators were upregulated upon supernatant exposure. Altogether, our results reveal nuanced competitive responses of P. aeruginosa when confronted with B. cenocepacia supernatant, integrating both environmental and social cues.

Persistent glucose consumption under antibiotic treatment protects bacterial community

Antibiotics typically induce major physiological changes in bacteria. However, their effect on nutrient consumption remains unclear. Here we found that Escherichia coli communities can sustain normal levels of glucose consumption under a broad range of antibiotics. The community-living resulted in a low membrane potential in the bacteria, allowing slow antibiotic accumulation on treatment and better adaptation. Through multi-omics analysis, we identified a prevalent adaptive response characterized by the upregulation of lipid synthesis, which substantially contributes to sustained glucose consumption. The consumption was maintained by the periphery region of the community, thereby restricting glucose penetration into the community interior. The resulting spatial heterogeneity in glucose availability protected the interior from antibiotic accumulation in a membrane potential-dependent manner, ensuring rapid recovery of the community postantibiotic treatment. Our findings unveiled a community-level antibiotic response through spatial regulation of metabolism and suggested new strategies for antibiotic therapies.

Characterization, whole-genome sequence analysis, and protease production of a new thermophilic Bacillus licheniformis strain isolated from Debagh hot spring, Algeria

A new thermophilic strain, designated as Bacillus sp. LMB3902, was isolated from Hammam Debagh, the hottest spring in Algeria (up to 98 °C). This isolate showed high protease production in skim milk media at 55 °C and exhibited significant specific protease activity by using azocasein as a substrate (157.50 U/mg). Through conventional methods, chemotaxonomic characteristics, 16S rRNA gene sequencing, and comparative genomic analysis with the closely related strain Bacillus licheniformis DSM 13 (ATCC 14580 T), the isolate Bacillus sp. LMB3902 was identified as a potentially new strain of Bacillus licheniformis. In addition, the gene functions of Bacillus sp. LMB3902 strain were predicted using the Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Clusters of Orthologous Groups, Non-Redundant Protein Sequence Database, Swiss-Prot, and Pfam databases. The results showed that the genome size of Bacillus sp. LMB3902 was 4.279.557 bp, with an average GC content of 46%. The genome contained 4.760 predicted genes, including 8 rRNAs, 78 tRNAs, and 24 sRNAs. A total of 235 protease genes were annotated including 50 proteases with transmembrane helix structures and eight secreted proteases with signal peptides. Additionally, the majority of secondary metabolites found by antiSMASH platform showed low similarity to identified natural products, such as fengicin (53%), lichenysin (57%), and surfactin (34%), suggesting that this strain may encode for novel uncharacterized natural products which can be useful for biotechnological applications. This study is the first report that describes the complete genome sequence, taxono-genomics, and gene annotation as well as protease production of the Bacillus genus in this hydrothermal vent.

Evolution and related pathogenic genes of Pseudodiploöspora longispora on Morchella based on genomic characterization and comparative genomic analysis

True morels (Morchella) are globally renowned medicinal and edible mushrooms. White mold disease caused by fungi is the main disease of Morchella, which has the characteristics of wide incidence and strong destructiveness. The disparities observed in the isolation rates of different pathogens indicate their varying degrees of host adaptability and competitive survival abilities. In order to elucidate its potential mechanism, this study, the pathogen of white mold disease from Dafang county, Guizhou Province was isolated and purified, identified as Pseudodiploöspora longispora by morphological, molecular biological and pathogenicity tests. Furthermore, high-quality genome of P. longisporus (40.846 Mb) was assembled N50 of 3.09 Mb, predicts 7381 protein-coding genes. Phylogenetic analysis of single-copy homologous genes showed that P. longispora and Zelopaecilomyces penicillatus have the closest evolutionary relationship, diverging into two branches approximately 50 (44.3-61.4) MYA. Additionally, compared with the other two pathogens causing Morchella disease, Z. penicillatus and Cladobotryum protrusum, it was found that they had similar proportions of carbohydrate enzyme types and encoded abundant cell wall degrading enzymes, such as chitinase and glucanase, indicating their important role in disease development. Moreover, the secondary metabolite gene clusters of P. longispora and Z. penicillatus show a high degree of similarity to leucinostatin A and leucinostatin B (peptaibols). Furthermore, a gene cluster with synthetic toxic substance Ochratoxin A was also identified in P. longispora and C. protrusum, indicating that they may pose a potential threat to food safety. This study provides valuable insights into the genome of P. longispora, contributing to pathogenicity research.

Chromosome-level genome assembly and annotation of the Rhabdophis nuchalis (Hubei keelback)

Rhabdophis nuchalis, a snake widely distributed in China, possesses a unique trait: glands beneath the skin on its neck and back, known as nucho-dorsal glands. These features make it a valuable subject for studying genetic diversity and the evolution of complex traits. In this study, we obtained a high-quality chromosome-level reference genome of R. nuchalis using MGI short-read sequencing, PacBio Revio long-read sequencing, and Hi-C sequencing techniques. The final assembly comprised 1.92 Gb of the R. nuchalis genome, anchored to 20 chromosomes (including 9 macrochromosomes and 11 microchromosomes), with a contig N50 of 104.79 Mb, a scaffold N50 of 204.96 Mb, and a BUSCO completeness of 97.50%. Additionally, we annotated a total of 1.09 Gb of repetitive sequences (which constitute 56.51% of the entire genome) and identified 22,057 protein-coding genes. This high-quality reference genome of R. nuchalis furnishes essential genomic data for comprehending the genetic diversity and evolutionary history of the species, as well as for facilitating species conservation efforts and comparative genomics studies.

Assessing marine ecosystem health using multi-omic analysis of blue mussel liquid biopsies: A case study within a national marine park

Marine ecosystems are under escalating threats from myriad environmental stressors, necessitating a deeper understanding of their impact on biodiversity and the health of sentinel organisms. In this study, we carried out a spatiotemporal multi-omic analysis of liquid biopsies collected from mussels (Mytilus spp.) in marine ecosystems of a national park. We delved into the epigenomic, transcriptomic, glycomic, proteomic, and microbiomic profiles to unravel the intricate interplay between ecosystem biodiversity and mussels' biological response to their environments. Our analysis revealed temporal fluctuations in the alpha diversity of the circulating microbiome associated with human activities. Analysis of the hemolymphatic circulating cell-free DNA (ccfDNA) provided information on the biodiversity and the presence of potential pathogens. Epigenomic analysis revealed widespread hypomethylation sites within the mitochondrial (mtDNA). Comparative transcriptomic and glycomic analyses highlighted differences in metabolic pathways and genes associated with immune and wound healing functions. This study demonstrates the potential of multi-omic analysis of liquid biopsy in sentinel to provide a holistic view of human activities' environmental impacts on marine coastal ecosystems. Overall, this approach has the potential to enhance the effectiveness and efficiency of various conservation efforts, leading to more informed decision-making and better outcomes for biodiversity and ecosystem conservation.

First regional reference database of northern Adriatic diatom transcriptomes

Marine microbial communities form the basis for the functioning of marine ecosystems and the conservation of biodiversity. With the application of metagenomics and metatranscriptomics in marine environmental studies, significant progress has been made in analysing the functioning of microbial communities as a whole. These molecular techniques are highly dependent on reliable, well-characterised, comprehensive and taxonomically diverse sequenced reference transcriptomes of microbial organisms. Here we present a set of 12 individual transcriptome assemblies derived from 6 representative diatom species from the northern Adriatic Sea grown under 2 environmentally relevant growth conditions (phosphate replete vs. phosphate deprived). After filtering the reads and assembly, an average number of 64,932 transcripts per assembly was obtained, of which an average of 8856 were assigned to functionally known proteins. Of all assigned transcripts, an average of 6483 proteins were taxonomically assigned to diatoms (Bacillariophyta). On average, a higher number of assigned proteins was detected in the transcriptome assemblies of diatoms grown under replete media condition. On average, 50% of the mapped proteins were shared between the two growth conditions. All recorded proteins in the dataset were classified into 24 COG categories, with approximately 25% belonging to the unknown function and the remaining 75% belonging to all other categories. The resulting diatom reference database for the northern Adriatic, focussing on the response to nutrient limitation as characteristic for the region and predicted for the future world oceans, provides a valuable resource for analysing environmental metatranscriptome and metagenome data. Each northern Adriatic transcriptome can also be used by itself as a reference database for the (meta)transcriptomes and gene expression studies of the associated species that will be generated in the future.

Bioinformatic analysis of type III CRISPR systems reveals key properties and new effector families

Recognition of RNA from invading mobile genetic elements (MGE) prompts type III CRISPR systems to activate an HD nuclease domain and/or a nucleotide cyclase domain in the Cas10 subunit, eliciting an immune response. The cyclase domain can generate a range of nucleotide second messengers, which in turn activate a diverse family of ancillary effector proteins. These provide immunity by non-specific degradation of host and MGE nucleic acids or proteins, perturbation of membrane potentials, transcriptional responses, or the arrest of translation. The wide range of nucleotide activators and downstream effectors generates a complex picture that is gradually being resolved. Here, we carry out a global bioinformatic analysis of type III CRISPR loci in prokaryotic genomes, defining the relationships of Cas10 proteins and their ancillary effectors. Our study reveals that cyclic tetra-adenylate is by far the most common signalling molecule used and that many loci have multiple effectors. These typically share the same activator and may work synergistically to combat MGE. We propose four new candidate effector protein families and confirm experimentally that the Csm6-2 protein, a highly diverged, fused Csm6 effector, is a ribonuclease activated by cyclic hexa-adenylate.

The Development of a Fluorescent Microsatellite Marker Assay for the Pitaya Canker Pathogen (Neoscytalidium dimidiatum)

Pitaya canker, caused by Neoscytalidium dimidiatum, is a destructive disease that significantly threatens the safety of the pitaya industry. The authors of previous studies have mainly focused on its biological characteristics and chemical control. However, there are no molecular markers available thus far that can be used for the population genetics study of this pathogen. In the present study, a draft genome of N. dimidiatum with a total length of 41.46 MB was assembled in which 9863 coding genes were predicted and annotated. In particular, the microsatellite sequences in the draft genome were investigated. To improve the successful screening rate of potentially polymorphic microsatellite makers, another five N. dimidiatum isolates were resequenced and assembled. A total of eight pairs of polymorphic microsatellite primers were screened out based on the polymorphic microsatellite loci after investigating the sequencing and resequencing assemblies of the six isolates. A total of thirteen representative isolates sampled from different pitaya plantations were genotyped in order to validate the polymorphism of the resulting eight markers. The results indicated that these markers were able to distinguish the isolates well. Lastly, a neighbor-joining tree of 35 isolates, sampled from different pitaya plantations located in different regions, was constructed according to the genotypes of the eight molecular markers. The developed tree indicated that these molecular markers had sufficient genotyping capabilities for our test panel of isolates. In summary, we developed a set of polymorphic microsatellite markers in the following study that can effectively genotype and distinguish N. dimidiatum isolates and be utilized in the population genetics study of N. dimidiatum.

Unraveling the genomic secrets of Tritonibacter mobilis AK171: a plant growth-promoting bacterium isolated from Avicennia marina

The scarcity of freshwater resources resulting in a significant yield loss presents a pressing challenge in agriculture. To address this issue, utilizing abundantly available saline water could offer a smart solution. In this study, we demonstrate that the genome sequence rhizosphere bacterium Tritonibacter mobilis AK171, a halophilic marine bacterium recognized for its ability to thrive in saline and waterlogged environments, isolated from mangroves, has the remarkable ability to enable plant growth using saline irrigation. AK171 is characterized as rod-shaped cells, displays agile movement in free-living conditions, and adopts a rosette arrangement in static media. Moreover, The qualitative evaluation of PGP traits showed that AK171 could produce siderophores and IAA but could not solubilize phosphate nor produce hydrolytic enzymes it exhibits a remarkable tolerance to high temperatures and salinity. In this study, we conducted a comprehensive genome sequence analysis of T. mobilis AK171 to unravel the genetic mechanisms underlying its plant growth-promoting abilities in such challenging conditions. Our analysis revealed diverse genes and pathways involved in the bacterium's adaptation to salinity and waterlogging stress. Notably, T. mobilis AK171 exhibited a high level of tolerance to salinity and waterlogging through the activation of stress-responsive genes and the production of specific enzymes and metabolites. Additionally, we identified genes associated with biofilm formation, indicating its potential role in establishing symbiotic relationships with host plants. Furthermore, our analysis unveiled the presence of genes responsible for synthesizing antimicrobial compounds, including tropodithietic acid (TDA), which can effectively control phytopathogens. This genomic insight into T. mobilis AK171 provides valuable information for understanding the molecular basis of plant-microbial interactions in saline and waterlogged environments. It offers potential applications for sustainable agriculture in challenging conditions.

Correlation analysis of whole genome sequencing of a pathogenic Escherichia coli strain of Inner Mongolian origin

Anal swabs of 1-month-old Holstein calves with diarrhea were collected from an intensive cattle farm, and a highly pathogenic Escherichia coli strain was obtained by isolation and purification. To study the virulence and resistance genes of pathogenic E. coli that cause diarrhea in calves, a strain of E. coli E12 isolated from calf diarrhea samples was used as experimental material in this experiment, and the virulence of the E12 strain were identified by the mouse infection test, and the whole genome map of the E12 strain were obtained by whole-genome sequencing and analyzed for genome characterization. The results showed that the lethality of strain E12 was 100%, the total length of E12-encoded genes was 4,294,530 bp, Cluster of Orthologous Groups of proteins (COG) annotated to 4,194 functional genes, and the virulence genes of sequenced strain E12 were compared with the virulence genes of sequenced strain E12 from the Virulence Factors of Pathogenic Bacteria (VFDB), which contained a total of 366 virulence genes in sequenced strain E12. The analysis of virulence genes of E12 revealed a total of 52 virulence genes in the iron transferrin system, 56 virulence genes in the secretory system, 41 virulence genes in bacterial toxins, and a total of 217 virulence genes in the Adhesin and Invasins group. The antibiotic resistance genes of sequenced strain E12 were identified through the Antibiotic Resistance Genes Database (ARDB) and Comprehensive Antibiotic Research Database, and it was found that its chromosome and plasmid included a total of 127 antibiotic resistance genes in four classes, and that E12 carried 71 genes related to the antibiotic efflux pumps, 36 genes related to antibiotic inactivation, and 14 antibiotic target alteration and reduced penetration into antibiotics, and 6 antibiotic resistance genes, and the resistance phenotypes were consistent with the genotypes. The pathogenic E. coli that causes diarrhea in calves on this ranch contains a large number of virulence and resistance genes. The results provide a theoretical basis for the prevention and treatment of diarrhea and other diseases caused by E. coli disease.