eggNOG 6.0: enabling comparative genomics across 12 535 organisms

Chromosome-scale Elaeis guineensis and E. oleifera assemblies: comparative genomics of oil palm and other Arecaceae

Elaeis guineensis and E. oleifera are the two species of oil palm. E. guineensis is the most widely cultivated commercial species, and introgression of desirable traits from E. oleifera is ongoing. We report an improved E. guineensis genome assembly with substantially increased continuity and completeness, as well as the first chromosome-scale E. oleifera genome assembly. Each assembly was obtained by integration of long-read sequencing, proximity ligation sequencing, optical mapping, and genetic mapping. High interspecific genome conservation is observed between the two species. The study provides the most extensive gene annotation to date, including 46,697 E. guineensis and 38,658 E. oleifera gene predictions. Analyses of repetitive element families further resolve the DNA repeat architecture of both genomes. Comparative genomic analyses identified experimentally validated small structural variants between the oil palm species and resolved the mechanism of chromosomal fusions responsible for the evolutionary descending dysploidy from 18 to 16 chromosomes.

The Spartina alterniflora genome sequence provides insights into the salt-tolerance mechanisms of exo-recretohalophytes

Spartina alterniflora is an exo-recretohalophyte Poaceae species that is able to grow well in seashore, but the genomic basis underlying its adaptation to salt tolerance remains unknown. Here, we report a high-quality, chromosome-level genome assembly of S. alterniflora constructed through PacBio HiFi sequencing, combined with high-throughput chromosome conformation capture (Hi-C) technology and Illumina-based transcriptomic analyses. The final 1.58 Gb genome assembly has a contig N50 size of 46.74 Mb. Phylogenetic analysis suggests that S. alterniflora diverged from Zoysia japonica approximately 21.72 million years ago (MYA). Moreover, whole-genome duplication (WGD) events in S. alterniflora appear to have expanded gene families and transcription factors relevant to salt tolerance and adaptation to saline environments. Comparative genomics analyses identified numerous species-specific genes, significantly expanded genes and positively selected genes that are enriched for 'ion transport' and 'response to salt stress'. RNA-seq analysis identified several ion transporter genes including the high-affinity K+ transporters (HKTs), SaHKT1;2, SaHKT1;3 and SaHKT1;8, and high copy number of Salt Overly Sensitive (SOS) up-regulated under high salt conditions, and the overexpression of SaHKT2;4 in Arabidopsis thaliana conferred salt tolerance to the plant, suggesting specialized roles for S. alterniflora to adapt to saline environments. Integrated metabolomics and transcriptomics analyses revealed that salt stress activate glutathione metabolism, with differential expressions of several genes such as γ-ECS, GSH-S, GPX, GST and PCS in the glutathione metabolism. This study suggests several adaptive mechanisms that could contribute our understanding of evolutional basis of the halophyte.

Comprehensive double-mutant analysis of the Bacillus subtilis envelope using double-CRISPRi

Understanding bacterial gene function remains a major biological challenge. Double-mutant genetic interaction (GI) analysis addresses this challenge by uncovering the functional partners of targeted genes, allowing us to associate genes of unknown function with novel pathways and unravel connections between well-studied pathways, but is difficult to implement at the genome-scale. Here, we develop and use double-CRISPRi to systematically quantify genetic interactions at scale in the Bacillus subtilis envelope, including essential genes. We discover > 1000 known and novel genetic interactions. Our analysis pipeline and experimental follow-ups reveal the distinct roles of paralogous genes such as the mreB and mbl actin homologs, and identify new genes involved in the well-studied process of cell division. Overall, our study provides valuable insights into gene function and demonstrates the utility of double-CRISPRi for high-throughput dissection of bacterial gene networks, providing a blueprint for future studies in diverse bacterial species.

Radiation-Tolerant Fibrivirga spp. from Rhizosphere Soil: Genome Insights and Potential in Agriculture

The rhizosphere of plants contains a wide range of microorganisms that can be cultivated and used for the benefit of agricultural practices. From garden soil near the rhizosphere region, Strain ES10-3-2-2 was isolated, and the cells were Gram-negative, aerobic, non-spore-forming rods that were 0.3-0.8 µm in diameter and 1.5-2.5 µm in length. The neighbor-joining method on 16S rDNA similarity revealed that the strain exhibited the highest sequence similarities with "Fibrivirga algicola JA-25" (99.2%) and Fibrella forsythia HMF5405T (97.3%). To further explore its biotechnological potentialities, we sequenced the complete genome of this strain employing the PacBio RSII sequencing platform. The genome of Strain ES10-3-2-2 comprises a 6,408,035 bp circular chromosome with a 52.8% GC content, including 5038 protein-coding genes and 52 RNA genes. The sequencing also identified three plasmids measuring 212,574 bp, 175,683 bp, and 81,564 bp. Intriguingly, annotations derived from the NCBI-PGAP, eggnog, and KEGG databases indicated the presence of genes affiliated with radiation-resistance pathway genes and plant-growth promotor key/biofertilization-related genes regarding Fe acquisition, K and P assimilation, CO2 fixation, and Fe solubilization, with essential roles in agroecosystems, as well as genes related to siderophore regulation. Additionally, T1SS, T6SS, and T9SS secretion systems are present in this species, like plant-associated bacteria. The inoculation of Strain ES10-3-2-2 to Arabidopsis significantly increases the fresh shoot and root biomass, thereby maintaining the plant quality compared to uninoculated controls. This work represents a link between radiation tolerance and the plant-growth mechanism of Strain ES10-3-2-2 based on in vitro experiments and bioinformatic approaches. Overall, the radiation-tolerant bacteria might enable the development of microbiological preparations that are extremely effective at increasing plant biomass and soil fertility, both of which are crucial for sustainable agriculture.

Exploring a unique class of flavoenzymes: Identification and biochemical characterization of ribosomal RNA dihydrouridine synthase

Dihydrouridine (D), a prevalent and evolutionarily conserved base in the transcriptome, primarily resides in tRNAs and, to a lesser extent, in mRNAs. Notably, this modification is found at position 2449 in the Escherichia coli 23S rRNA, strategically positioned near the ribosome's peptidyl transferase site. Despite the prior identification, in E. coli genome, of three dihydrouridine synthases (DUS), a set of NADPH and FMN-dependent enzymes known for introducing D in tRNAs and mRNAs, characterization of the enzyme responsible for D2449 deposition has remained elusive. This study introduces a rapid method for detecting D in rRNA, involving reverse transcriptase-blockage at the rhodamine-labeled D2449 site, followed by PCR amplification (RhoRT-PCR). Through analysis of rRNA from diverse E. coli strains, harboring chromosomal or single-gene deletions, we pinpoint the yhiN gene as the ribosomal dihydrouridine synthase, now designated as RdsA. Biochemical characterizations uncovered RdsA as a unique class of flavoenzymes, dependent on FAD and NADH, with a complex structural topology. In vitro assays demonstrated that RdsA dihydrouridylates a short rRNA transcript mimicking the local structure of the peptidyl transferase site. This suggests an early introduction of this modification before ribosome assembly. Phylogenetic studies unveiled the widespread distribution of the yhiN gene in the bacterial kingdom, emphasizing the conservation of rRNA dihydrouridylation. In a broader context, these findings underscore nature's preference for utilizing reduced flavin in the reduction of uridines and their derivatives.

Nutrient availability contributes to structural and functional diversity of microbiome in Xinjiang oilfield

Indigenous microbial enhanced oil recovery (IMEOR) is a promising alternative way to promote oil recovery. It activates oil recovery microorganisms in the reservoir by adding nutrients to the injected water, utilizing microbial growth and metabolism to enhance recovery. However, few studies have focused on the impact of injected nutrients on reservoir microbial community composition and potential functions. This limits the further strategic development of IMEOR. In this study, we investigated the effects of nutrition on the composition of the reservoir bacterial community and functions in the Qizhong block of Xinjiang Oilfield, China, by constructing a long core microbial flooding simulation device. The results showed that the microbial community structure of the reservoir changed from aerobic state to anaerobic state after nutrient injection. Reducing the nutrient concentration increased the diversity and network stability of the reservoir bacterial community. At the same time, the nitrogen metabolism function also showed the same change response. Overall, these results indicated that nutrition significantly affected the community structure and function of reservoir microorganisms. Injecting low concentrations of nutrients may be more beneficial to improve oil recovery. This study is of great significance for guiding IMEOR technology and saving costs at the field site.

VOGDB-Database of Virus Orthologous Groups

Computational models of homologous protein groups are essential in sequence bioinformatics. Due to the diversity and rapid evolution of viruses, the grouping of protein sequences from virus genomes is particularly challenging. The low sequence similarities of homologous genes in viruses require specific approaches for sequence- and structure-based clustering. Furthermore, the annotation of virus genomes in public databases is not as consistent and up to date as for many cellular genomes. To tackle these problems, we have developed VOGDB, which is a database of virus orthologous groups. VOGDB is a multi-layer database that progressively groups viral genes into groups connected by increasingly remote similarity. The first layer is based on pair-wise sequence similarities, the second layer is based on the sequence profile alignments, and the third layer uses predicted protein structures to find the most remote similarity. VOGDB groups allow for more sensitive homology searches of novel genes and increase the chance of predicting annotations or inferring phylogeny. VOGD B uses all virus genomes from RefSeq and partially reannotates them. VOGDB is updated with every RefSeq release. The unique feature of VOGDB is the inclusion of both prokaryotic and eukaryotic viruses in the same clustering process, which makes it possible to explore old evolutionary relationships of the two groups. VOGDB is freely available at vogdb.org under the CC BY 4.0 license.

Large-scale investigation of species-specific orphan genes in the human gut microbiome elucidates their evolutionary origins

Species-specific genes, also known as orphans, are ubiquitous across life's domains. In prokaryotes, species-specific orphan genes (SSOGs) are mostly thought to originate in external elements such as viruses followed by horizontal gene transfer, whereas the scenario of native origination, through rapid divergence or de novo, is mostly dismissed. However, quantitative evidence supporting either scenario is lacking. Here, we systematically analyzed genomes from 4644 human gut microbiome species and identified more than 600,000 unique SSOGs, representing an average of 2.6% of a given species' pangenome. These sequences are mostly rare within each species yet show signs of purifying selection. Overall, SSOGs use optimal codons less frequently, and their proteins are more disordered than those of conserved genes (i.e., non-SSOGs). Importantly, across species, the GC content of SSOGs closely matches that of conserved ones. In contrast, the ∼5% of SSOGs that share similarity to known viral sequences have distinct characteristics, including lower GC content. Thus, SSOGs with similarity to viruses differ from the remaining SSOGs, contrasting an external origination scenario for most of them. By examining the orthologous genomic region in closely related species, we show that a small subset of SSOGs likely evolved natively de novo and find that these genes also differ in their properties from the remaining SSOGs. Our results challenge the notion that external elements are the dominant source of prokaryotic genetic novelty and will enable future studies into the biological role and relevance of species-specific genes in the human gut.

The TSC22D, WNK, and NRBP gene families exhibit functional buffering and evolved with Metazoa for cell volume regulation

The ability to sense and respond to osmotic fluctuations is critical for the maintenance of cellular integrity. We used gene co-essentiality analysis to identify an unappreciated relationship between TSC22D2, WNK1, and NRBP1 in regulating cell volume homeostasis. All of these genes have paralogs and are functionally buffered for osmo-sensing and cell volume control. Within seconds of hyperosmotic stress, TSC22D, WNK, and NRBP family members physically associate into biomolecular condensates, a process that is dependent on intrinsically disordered regions (IDRs). A close examination of these protein families across metazoans revealed that TSC22D genes evolved alongside a domain in NRBPs that specifically binds to TSC22D proteins, which we have termed NbrT (NRBP binding region with TSC22D), and this co-evolution is accompanied by rapid IDR length expansion in WNK-family kinases. Our study reveals that TSC22D, WNK, and NRBP genes evolved in metazoans to co-regulate rapid cell volume changes in response to osmolarity.

Effects of vegetation cover and aquaculture pollution on viral assemblages in mangroves sediments

Mangrove forests, a critical coastal ecosystem, face numerous anthropogenic threats, particularly from aquaculture activities. Despite the acknowledged significance of viruses in local and global biogeochemical cycles, there is limited knowledge regarding the community structure, genomic diversity, and ecological roles of viruses in mangrove forests ecosystems, especially regarding their responses to aquaculture. In this study, we identified 17,755 viral operational taxonomic units (vOTUs) from nine sediments viromes across three distinct ecological regions of the mangrove forests ecosystem: mangrove, bare flat, and aquaculture regions. Viral assemblages varied among three regions, and the pathogenic viruses associated with marine animals, such as the white spot syndrome virus (WSSV) from Nimaviridae, were identified in this study. The relative abundance of Nimaviridae in the bare flat region was higher than in other regions. Furthermore, viruses in distinct mangrove forests sediments regions have adapted to their environments by adopting distinct survival strategies and encoding various auxiliary metabolic genes involved in carbon metabolism and antibiotic resistance. These adaptations may have profound impacts on biogeochemical cycles. This study provides the first insights into the effects of vegetation cover and aquaculture on the community structure and ecological roles of viruses in mangrove forests sediments. These findings are crucial for understanding the risks posed by anthropogenic threats to mangrove forests ecosystems and informing effective management strategies.

Lysine 2-hydroxyisobutyrylation proteomics analyses reveal the regulatory mechanism of CaMYB61-CaAFR1 module in regulating stem development in Capsicum annuum L

Plant stems constitute the most abundant renewable resource on earth. The function of lysine (K)-2-hydroxyisobutyrylation (Khib), a novel post-translational modification (PTM), has not yet been elucidated in plant stem development. Here, by assessing typical pepper genotypes with straight stem (SS) and prostrate stem (PS), we report the first large-scale proteomics analysis for protein Khib to date. Khib-modifications influenced central metabolic processes involved in stem development, such as glycolysis/gluconeogenesis and protein translation. The high Khib level regulated gene expression and protein accumulation associated with cell wall formation in the pepper stem. Specially, we found that CaMYB61 knockdown lines that exhibited prostrate stem phenotypes had high Khib levels. Most histone deacetylases (HDACs, e.g., switch-independent 3 associated polypeptide function related 1, AFR1) potentially function as the "erasing enzymes" involved in reversing Khib level. CaMYB61 positively regulated CaAFR1 expression to erase Khib and promote cellulose and hemicellulose accumulation in the stem. Therefore, we propose a bidirectional regulation hypothesis of "Khib modifications" and "Khib erasing" in stem development, and reveal a novel epigenetic regulatory network in which the CaMYB61-CaAFR1 molecular module participating in the regulation of Khib levels and biosynthesis of cellulose and hemicellulose for the first time.

Characterization and Genome Analysis of the Delftia lacustris Strain LzhVag01 Isolated from Vaginal Discharge

Delftia has been separated from freshwater, sludge, and soil and has emerged as a novel opportunistic pathogen in the female vagina. However, the genomic characteristics, pathogenicity, and biotechnological properties still need to be comprehensively investigated. In this study, a Delftia strain was isolated from the vaginal discharge of a 43-year-old female with histologically confirmed cervical intraepithelial neoplasm (CIN III), followed by whole-genome sequencing. Phylogenetic analysis and average nucleotide identity (ANI) analysis demonstrated that it belongs to Delftia lacustris, named D. lacustris strain LzhVag01. LzhVag01 was sensitive to β-lactams, macrolides, and tetracyclines but exhibited resistance to lincoamines, nitroimidazoles, aminoglycosides, and fluoroquinolones. Its genome is a single, circular chromosome of 6,740,460 bp with an average GC content of 66.59%. Whole-genome analysis identified 16 antibiotic resistance-related genes, which match the antimicrobial susceptibility profile of this strain, and 11 potential virulence genes. These pathogenic factors may contribute to its colonization in the vaginal environment and its adaptation and accelerate the progression of cervical cancer. This study sequenced and characterized the whole-genome of Delftia lacustris isolated from vaginal discharge, which provides investigators and clinicians with valuable insights into this uncommon species.

Aureococcus anophagefferens strain CCMP1851: draft genome of a second Kratosvirus quantuckense-susceptible host strain for an emerging host-giant virus model system

Here, we report the draft genome of Aureococcus anophagefferens strain CCMP1851, which is susceptible to the virus Kratosvirus quantuckense. CCMP1851 complements an available genome for a virus-resistant strain (CCMP1850) isolated from the same bloom. Future studies can now use this genome to examine genetic hints of virus resistance and susceptibility.

BGCFlow: systematic pangenome workflow for the analysis of biosynthetic gene clusters across large genomic datasets

Genome mining is revolutionizing natural products discovery efforts. The rapid increase in available genomes demands comprehensive computational platforms to effectively extract biosynthetic knowledge encoded across bacterial pangenomes. Here, we present BGCFlow, a novel systematic workflow integrating analytics for large-scale genome mining of bacterial pangenomes. BGCFlow incorporates several genome analytics and mining tools grouped into five common stages of analysis such as: (i) data selection, (ii) functional annotation, (iii) phylogenetic analysis, (iv) genome mining, and (v) comparative analysis. Furthermore, BGCFlow provides easy configuration of different projects, parallel distribution, scheduled job monitoring, an interactive database to visualize tables, exploratory Jupyter Notebooks, and customized reports. Here, we demonstrate the application of BGCFlow by investigating the phylogenetic distribution of various biosynthetic gene clusters detected across 42 genomes of the Saccharopolyspora genus, known to produce industrially important secondary/specialized metabolites. The BGCFlow-guided analysis predicted more accurate dereplication of BGCs and guided the targeted comparative analysis of selected RiPPs. The scalable, interoperable, adaptable, re-entrant, and reproducible nature of the BGCFlow will provide an effective novel way to extract the biosynthetic knowledge from the ever-growing genomic datasets of biotechnologically relevant bacterial species.

A long-read draft assembly of the Chinese mantis (Mantodea: Mantidae: Tenodera sinensis) genome reveals patterns of ion channel gain and loss across Arthropoda

Praying mantids (Mantodea: Mantidae) are iconic insects that have captivated biologists for decades, especially the species with cannibalistic copulatory behavior. This behavior has been cited as evidence that insects lack nociceptive capacities and cannot feel pain; however, this behaviorally driven hypothesis has never been rigorously tested at the genetic or functional level. To enable future studies of nociceptive capabilities in mantids, we sequenced and assembled a draft genome of the Chinese praying mantis (Tenodera sinensis) and identified multiple classes of nociceptive ion channels by comparison to orthologous gene families in Arthropoda. Our assembly-produced using PacBio HiFi reads-is fragmented (total size = 3.03 Gb; N50 = 1.8 Mb; 4,966 contigs), but is highly complete with respect to gene content (BUSCO complete = 98.7% [odb10_insecta]). The size of our assembly is substantially larger than that of most other insects, but is consistent with the size of other mantid genomes. We found that most families of nociceptive ion channels are present in the T. sinensis genome; that they are most closely related to those found in the damp-wood termite (Zootermopsis nevadensis); and that some families have expanded in T. sinensis while others have contracted relative to nearby lineages. Our findings suggest that mantids are likely to possess nociceptive capabilities and provide a foundation for future experimentation regarding ion channel functions and their consequences for insect behavior.

Neuronal differentiation requires BRAT1 complex to remove REST from chromatin

Repressor element-1 silencing transcription factor (REST) is required for the formation of mature neurons. REST dysregulation underlies a key mechanism of neurodegeneration associated with neurological disorders. However, the mechanisms leading to alterations of REST-mediated silencing of key neurogenesis genes are not known. Here, we show that BRCA1 Associated ATM Activator 1 (BRAT1), a gene linked to neurodegenerative diseases, is required for the activation of REST-responsive genes during neuronal differentiation. We find that INTS11 and INTS9 subunits of Integrator complex interact with BRAT1 as a distinct trimeric complex to activate critical neuronal genes during differentiation. BRAT1 depletion results in persistence of REST residence on critical neuronal genes disrupting the differentiation of NT2 cells into astrocytes and neuronal cells. We identified BRAT1 and INTS11 co-occupying the promoter region of these genes and pinpoint a role for BRAT1 in recruiting INTS11 to their promoters. Disease-causing mutations in BRAT1 diminish its association with INTS11/INTS9, linking the manifestation of disease phenotypes with a defect in transcriptional activation of key neuronal genes by BRAT1/INTS11/INTS9 complex. Finally, loss of Brat1 in mouse embryonic stem cells leads to a defect in neuronal differentiation assay. Importantly, while reconstitution with wild-type BRAT1 restores neuronal differentiation, the addition of a BRAT1 mutant is unable to associate with INTS11/INTS9 and fails to rescue the neuronal phenotype. Taken together, our study highlights the importance of BRAT1 association with INTS11 and INTS9 in the development of the nervous system.

Utilizing profile hidden Markov model databases for discovering viruses from metagenomic data: a comprehensive review

Profile hidden Markov models (pHMMs) are able to achieve high sensitivity in remote homology search, making them popular choices for detecting novel or highly diverged viruses in metagenomic data. However, many existing pHMM databases have different design focuses, making it difficult for users to decide the proper one to use. In this review, we provide a thorough evaluation and comparison for multiple commonly used profile HMM databases for viral sequence discovery in metagenomic data. We characterized the databases by comparing their sizes, their taxonomic coverage, and the properties of their models using quantitative metrics. Subsequently, we assessed their performance in virus identification across multiple application scenarios, utilizing both simulated and real metagenomic data. We aim to offer researchers a thorough and critical assessment of the strengths and limitations of different databases. Furthermore, based on the experimental results obtained from the simulated and real metagenomic data, we provided practical suggestions for users to optimize their use of pHMM databases, thus enhancing the quality and reliability of their findings in the field of viral metagenomics.

A high-quality chromosome-level genome assembly of Ficus hirta

Ficus species (Moraceae) play pivotal roles in tropical and subtropical ecosystems. Thriving across diverse habitats, from rainforests to deserts, they harbor a multitude of mutualistic and antagonistic interactions with insects, nematodes, and pathogens. Despite their ecological significance, knowledge about the genomic background of Ficus remains limited. In this study, we report a chromosome-level reference genome of F. hirta, with a total size of 297.27 Mb, containing 28,625 protein-coding genes and 44.67% repeat sequences. These findings illuminate the genetic basis of Ficus responses to environmental challenges, offering valuable genomic resources for understanding genome size, adaptive evolution, and co-evolution with natural enemies and mutualists within the genus.

Paracoccus benzoatiresistens sp. nov., a benzoate resistance and selenite reduction bacterium isolated from wetland

A Gram-stain-negative, aerobic, non-motile, catalase- and oxidase-positive, pale orange, rod-shaped strain EF6T, was isolated from a natural wetland reserve in Hebei province, China. The strain grew at 25-37 °C (optimum, 30 °C), pH 5-9 (optimum, pH 7), and in the presence of 1.0-4.0% (w/v) NaCl (optimum, 2%). A phylogenetic analysis based on 16S rRNA gene sequence revealed that strain EF6T belongs to the genus Paracoccus, and the closest members were Paracoccus shandongensis wg2T with 98.1% similarity, Paracoccus fontiphilus MVW-1 T (97.9%), Paracoccus everestensis S8-55 T (97.7%), Paracoccus subflavus GY0581T (97.6%), Paracoccus sediminis CMB17T (97.3%), Paracoccus caeni MJ17T (97.0%), and Paracoccus angustae E6T (97.0%). The genome size of strain EF6T was 4.88 Mb, and the DNA G + C content was 65.3%. The digital DNA-DNA hybridization, average nucleotide identity, and average amino acid identity values between strain EF6T and the reference strains were all below the threshold limit for species delineation (< 32.8%, < 88.0%, and < 86.7%, respectively). The major fatty acids (≥ 5.0%) were summed feature 8 (86.3%, C18:1 ω6c and/or C18:1 ω7c) and C18:1 (5.0%) and the only isoprenoid quinone was Q-10. The polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol, two unidentified glycolipids, five unidentified phospholipids, and an unidentified aminolipid. Strain EF6T displays notable resistance to benzoate and selenite, with higher tolerance levels (25 g/L for benzoate and 150 mM for selenite) compared to the closely related species. Genomic analysis identified six benzoate resistance genes (acdA, pcaF, fadA, pcaC, purB, and catA) and twenty selenite resistance and reduction-related genes (iscR, ssuB, ssuD, selA, selD and so on). Additionally, EF6T possesses unique genes (catA, ssuB, and ssuC) absent in the closely related species for benzoate and selenite resistance. Its robust resistance to benzoate and selenite, coupled with its genomic makeup, make EF6T a promising candidate for the remediation of both organic and inorganic pollutants. It is worth noting that the specific resistance phenotypes described above were not reported in other novel species in Paracoccus. Based on the results of biochemical, physiological, phylogenetic, and chemotaxonomic analyses, combined with comparisons of the 16S rRNA gene sequence and the whole genome sequence, strain EF6T is considered to represent a novel species of the genus Paracoccus within the family Rhodobacteraceae, for which the name Paracoccus benzoatiresistens sp. nov. is proposed. The type strain is EF6T (= GDMCC 1.3400 T = JCM 35642 T = MCCC 1K08702T).

Genome assembly of Melilotus officinalis provides a new reference genome for functional genomics

BACKGROUND

Sweet yellow clover (Melilotus officinalis) is a diploid plant (2n = 16) that is native to Europe. It is an excellent legume forage. It can both fix nitrogen and serve as a medicine. A genome assembly of Melilotus officinalis that was collected from Best corporation in Beijing is available based on Nanopore sequencing. The genome of Melilotus officinalis was sequenced, assembled, and annotated.

RESULTS

The latest PacBio third generation HiFi assembly and sequencing strategies were used to produce a Melilotus officinalis genome assembly size of 1,066 Mbp, contig N50 = 5 Mbp, scaffold N50 = 130 Mbp, and complete benchmarking universal single-copy orthologs (BUSCOs) = 96.4%. This annotation produced 47,873 high-confidence gene models, which will substantially aid in our research on molecular breeding. A collinear analysis showed that Melilotus officinalis and Medicago truncatula shared conserved synteny. The expansion and contraction of gene families showed that Melilotus officinalis expanded by 565 gene families and shrank by 56 gene families. The contacted gene families were associated with response to stimulus, nucleotide binding, and small molecule binding. Thus, it is related to a family of genes associated with peptidase activity, which could lead to better stress tolerance in plants.

CONCLUSIONS

In this study, the latest PacBio technology was used to assemble and sequence the genome of the Melilotus officinalis and annotate its protein-coding genes. These results will expand the genomic resources available for Melilotus officinalis and should assist in subsequent research on sweet yellow clover plants.

A fast comparative genome browser for diverse bacteria and archaea

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

Microbial Diversity and Open Questions about the Deep Tree of Life

In this perspective, we explore the transformative impact and inherent limitations of metagenomics and single-cell genomics on our understanding of microbial diversity and their integration into the Tree of Life. We delve into the key challenges associated with incorporating new microbial lineages into the Tree of Life through advanced phylogenomic approaches. Additionally, we shed light on enduring debates surrounding various aspects of the microbial Tree of Life, focusing on recent advances in some of its deepest nodes, such as the roots of bacteria, archaea, and eukaryotes. We also bring forth current limitations in genome recovery and phylogenomic methodology, as well as new avenues of research to uncover additional key microbial lineages and resolve the shape of the Tree of Life.

Regional centromere configuration in the fungal pathogens of the Pneumocystis genus

Centromeres are constricted chromosomal regions that are essential for cell division. In eukaryotes, centromeres display a remarkable architectural and genetic diversity. The basis of centromere-accelerated evolution remains elusive. Here, we focused on Pneumocystis species, a group of mammalian-specific fungal pathogens that form a sister taxon with that of the Schizosaccharomyces pombe, an important genetic model for centromere biology research. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of S. pombe. Using organisms from a short-term in vitro culture or infected animal models and chromatin immunoprecipitation (ChIP)-Seq, we identified CENP-A bound regions in two Pneumocystis species that diverged ~35 million years ago. Each species has a unique short regional centromere (<10 kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. These features suggest an epigenetic specification of centromere function. Analysis of centromeric DNA across multiple Pneumocystis species suggests a vertical transmission at least 100 million years ago. The common ancestry of Pneumocystis and S. pombe centromeres is untraceable at the DNA level, but the overall architectural similarity could be the result of functional constraint for successful chromosomal segregation.IMPORTANCEPneumocystis species offer a suitable genetic system to study centromere evolution in pathogens because of their phylogenetic proximity with the non-pathogenic yeast S. pombe, a popular model for cell biology. We used this system to explore how centromeres have evolved after the divergence of the two clades ~ 460 million years ago. To address this question, we established a protocol combining short-term culture and ChIP-Seq to characterize centromeres in multiple Pneumocystis species. We show that Pneumocystis have short epigenetic centromeres that function differently from those in S. pombe.

MetaCerberus: distributed highly parallelized HMM-based processing for robust functional annotation across the tree of life

MOTIVATION

MetaCerberus is a massively parallel, fast, low memory, scalable annotation tool for inference gene function across genomes to metacommunities. MetaCerberus provides an elusive HMM/HMMER-based tool at a rapid scale with low memory. It offers scalable gene elucidation to major public databases, including KEGG (KO), COGs, CAZy, FOAM, and specific databases for viruses, including VOGs and PHROGs, from single genomes to metacommunities.

RESULTS

MetaCerberus is 1.3× as fast on a single node than eggNOG-mapper v2 on 5× less memory using an exclusively HMM/HMMER mode. In a direct comparison, MetaCerberus provides better annotation of viruses, phages, and archaeal viruses than DRAM, Prokka, or InterProScan. MetaCerberus annotates more KOs across domains when compared to DRAM, with a 186× smaller database, and with 63× less memory. MetaCerberus is fully integrated for automatic analysis of statistics and pathways using differential statistic tools (i.e. DESeq2 and edgeR), pathway enrichment (GAGE R), and pathview R. MetaCerberus provides a novel tool for unlocking the biosphere across the tree of life at scale.

AVAILABILITY AND IMPLEMENTATION

MetaCerberus is written in Python and distributed under a BSD-3 license. The source code of MetaCerberus is freely available at https://github.com/raw-lab/metacerberus compatible with Python 3 and works on both Mac OS X and Linux. MetaCerberus can also be easily installed using bioconda: mamba create -n metacerberus -c bioconda -c conda-forge metacerberus.

Unveiling a high-risk epidemic clone (ST 357) of 'Difficult to Treat Extensively Drug-Resistant' (DT-XDR) Pseudomonas aeruginosa from a burn patient in Bangladesh: A resilient beast revealing coexistence of four classes of beta lactamases

OBJECTIVES

Pseudomonas aeruginosa (P. aeruginosa) stands out as a key culprit in the colonization of burn wounds, instigating grave infections of heightened severity. In this study, we have performed comparative whole genome analysis of a difficult to treat extensively drug resistant P. aeruginosa isolated from a burn patient in order to elucidate genomic diversity, molecular patterns, mechanisms and genes responsible for conferring antimicrobial resistance and virulence.

METHOD

P. aeruginosa SHNIBPS206 was isolated from an infected burn wound of a critically injured burn patient. Whole genome sequencing was carried out and annotated with Prokka. Sequence type, serotype, antimicrobial resistance genes and mechanisms, virulence genes, metal resistance genes and CRISPR/Cas systems were investigated. Later, pangenome analysis was carried out to find out genomic diversity.

RESULT

P. aeruginosa SHNIBPS206 (MLST 357, Serotype O11) was resistant to 14 antibiotics including carbapenems and harboured all four classes of beta lactamase producing genes: Class A (blaPME-1, blaVEB-9), Class B (blaNDM-1), Class C (blaPDC-11) and Class D (blaOXA-846). Mutational analysis of Porin D gave valuable insights. Several efflux pump, virulence and metal resistance genes were also detected. Pangenome analysis revealed high genomic diversity among different strains of P. aeruginosa.

CONCLUSION

To our knowledge, this is the first report of an extensively drug resistant ST 357 P. aeruginosa from Bangladesh, which is an epidemic high-risk P. aeruginosa clone. Further research and in-depth comprehensive studies are required to investigate the prevalence of such high-risk clone of P. aeruginosa in Bangladesh.

Beyond blast: enabling microbiologists to better extract literature, taxonomic distributions and gene neighbourhood information for protein families

Capturing the published corpus of information on all members of a given protein family should be an essential step in any study focusing on specific members of that family. Using a previously gathered dataset of more than 280 references mentioning a member of the DUF34 (NIF3/Ngg1-interacting Factor 3) family, we evaluated the efficiency of different databases and search tools, and devised a workflow that experimentalists can use to capture the most information published on members of a protein family in the least amount of time. To complement this workflow, web-based platforms allowing for the exploration of protein family members across sequenced genomes or for the analysis of gene neighbourhood information were reviewed for their versatility and ease of use. Recommendations that can be used for experimentalist users, as well as educators, are provided and integrated within a customized, publicly accessible Wiki.

Unraveling the multiplicity of geranylgeranyl reductases in Archaea: potential roles in saturation of terpenoids

The enzymology of the key steps in the archaeal phospholipid biosynthetic pathway has been elucidated in recent years. In contrast, the complete biosynthetic pathways for proposed membrane regulators consisting of polyterpenes, such as carotenoids, respiratory quinones, and polyprenols remain unknown. Notably, the multiplicity of geranylgeranyl reductases (GGRs) in archaeal genomes has been correlated with the saturation of polyterpenes. Although GGRs, which are responsible for saturation of the isoprene chains of phospholipids, have been identified and studied in detail, there is little information regarding the structure and function of the paralogs. Here, we discuss the diversity of archaeal membrane-associated polyterpenes which is correlated with the genomic loci, structural and sequence-based analyses of GGR paralogs.

Genomic landscape of NDM-1 producing multidrug-resistant Providencia stuartii causing burn wound infections in Bangladesh

The increasing antimicrobial resistance in Providencia stuartii (P. stuartii) worldwide, particularly concerning for immunocompromised and burn patients, has raised concern in Bangladesh, where the significance of this infectious opportunistic pathogen had been previously overlooked, prompting a need for investigation. The two strains of P. stuartii (P. stuartii SHNIBPS63 and P. stuartii SHNIBPS71) isolated from wound swab of two critically injured burn patients were found to be multidrug-resistant and P. stuartii SHNIBPS63 showed resistance to all the 22 antibiotics tested as well as revealed the co-existence of blaVEB-6 (Class A), blaNDM-1 (Class B), blaOXA-10 (Class D) beta lactamase genes. Complete resistance to carbapenems through the production of NDM-1, is indicative of an alarming situation as carbapenems are considered to be the last line antibiotic to combat this pathogen. Both isolates displayed strong biofilm-forming abilities and exhibited resistance to copper, zinc, and iron, in addition to carrying multiple genes associated with metal resistance and the formation of biofilms. The study also encompassed a pangenome analysis utilizing a dataset of eighty-six publicly available P. stuartii genomes (n = 86), revealing evidence of an open or expanding pangenome for P. stuartii. Also, an extensive genome-wide analysis of all the P. stuartii genomes revealed a concerning global prevalence of diverse antimicrobial resistance genes, with a particular alarm raised over the abundance of carbapenem resistance gene blaNDM-1. Additionally, this study highlighted the notable genetic diversity within P. stuartii, significant informations about phylogenomic relationships and ancestry, as well as potential for cross-species transmission, raising important implications for public health and microbial adaptation across different environments.

Proteomic differences between extracellular vesicles and extracellular vesicle-depleted excretory/secretory products of barber's pole worm

BACKGROUND

Components of excretory/secretory products (ESPs) of helminths have been proposed as vaccine targets and shown to play a role in modulating host immune responses for decades. Such research interest is further increased by the discovery of extracellular vesicles (EVs) in the ESPs of parasitic worms. Although efforts have been made to reveal the cargos of EVs, little is known about the proteomic differences between EVs and canonical ESPs released by parasitic worms from animals.

METHODS

The total ESPs of Haemonchus contortus (barber's pole worm) were obtained by short-term in vitro culturing of young adult worms, and small EVs were isolated from ESPs using an ultracentrifugation method. Data-dependent acquisition (DDA) label-free Nano-LC-MS/MS was used to quantify the proteomic difference between small EVs and EV-depleted ESPs of H. contortus. Functional annotation and enrichment of the differential proteins were performed regarding cellular components, molecular functions, pathways, and/or biological processes.

RESULTS

A total of 1697 proteins were identified in small EVs and EV-depleted ESPs of H. contortus adult worms, with 706 unique proteins detected in the former and 597 unique proteins in the latter. It was revealed that proteins in small EVs are dominantly cytoplasmic, whereas proteins in EV-depleted ESPs are mainly extracellular; canonical ESPs such as proteases and small GTPases were abundantly detected in small EVs, and SCP/TAP-, DUF-, and GLOBIN domain-containing proteins were mainly found in EV-depleted ESPs. Compared with well-characterised proteins in small EVs, about 50% of the proteins detected in EV-depleted ESPs were poorly characterised.

CONCLUSIONS

There are remarkable differences between small EVs and EV-depleted ESPs of H. contortus in terms of protein composition. Immune modulatory effects caused by nematode ESPs are possibly contributed mainly by the proteins in small EVs.

Comparative Bioinformatic Analysis of the Proteomes of Rabbit and Human Sex Chromosomes

Studying proteins associated with sex chromosomes can provide insights into sex-specific proteins. Membrane proteins accessible through the cell surface may serve as excellent targets for diagnostic, therapeutic, or even technological purposes, such as sperm sexing technologies. In this context, proteins encoded by sex chromosomes have the potential to become targets for X- or Y-chromosome-bearing spermatozoa. Due to the limited availability of proteomic studies on rabbit spermatozoa and poorly annotated databases for rabbits compared to humans, a bioinformatic analysis of the available rabbit X chromosome proteome (RX), as well as the human X (HX) and Y (HY) chromosomes proteome, was conducted to identify potential targets that could be accessible from the cell surface and predict which of the potential targets identified in humans might also exist in rabbits. We identified 100, 211, and 3 proteins associated with the plasma membrane or cell surface for RX, HX, and HY, respectively, of which 61, 132, and 3 proteins exhibit potential as targets as they were predicted to be accessible from the cell surface. Cross-referencing the potential HX targets with the rabbit proteome revealed an additional 60 proteins with the potential to be RX targets, resulting in a total of 121 potential RX targets. In addition, at least 53 possible common HX and RX targets have been previously identified in human spermatozoa, emphasizing their potential as targets of X-chromosome-bearing spermatozoa. Further proteomic studies on rabbit sperm will be essential to identify and validate the usefulness of these proteins for application in rabbit sperm sorting techniques as targets of X-chromosome-bearing spermatozoa.

IMP: bridging the gap for medicinal plant genomics

Medicinal plants have garnered significant attention in ethnomedicine and traditional medicine due to their potential antitumor, anti-inflammatory and antioxidant properties. Recent advancements in genome sequencing and synthetic biology have revitalized interest in natural products. Despite the availability of sequenced genomes and transcriptomes of these plants, the absence of publicly accessible gene annotations and tabular formatted gene expression data has hindered their effective utilization. To address this pressing issue, we have developed IMP (Integrated Medicinal Plantomics), a freely accessible platform at https://www.bic.ac.cn/IMP. IMP curated a total of 8 565 672 genes for 84 high-quality genome assemblies, and 2156 transcriptome sequencing samples encompassing various organs, tissues, developmental stages and stimulations. With the integrated 10 analysis modules, users could simply examine gene annotations, sequences, functions, distributions and expressions in IMP in a one-stop mode. We firmly believe that IMP will play a vital role in enhancing the understanding of molecular metabolic pathways in medicinal plants or plants with medicinal benefits, thereby driving advancements in synthetic biology, and facilitating the exploration of natural sources for valuable chemical constituents like drug discovery and drug production.

OrthoMaM v12: a database of curated single-copy ortholog alignments and trees to study mammalian evolutionary genomics

To date, the databases built to gather information on gene orthology do not provide end-users with descriptors of the molecular evolution information and phylogenetic pattern of these orthologues. In this context, we developed OrthoMaM, a database of ORTHOlogous MAmmalian Markers describing the evolutionary dynamics of coding sequences in mammalian genomes. OrthoMaM version 12 includes 15,868 alignments of orthologous coding sequences (CDS) from the 190 complete mammalian genomes currently available. All annotations and 1-to-1 orthology assignments are based on NCBI. Orthologous CDS can be mined for potential informative markers at the different taxonomic levels of the mammalian tree. To this end, several evolutionary descriptors of DNA sequences are provided for querying purposes (e.g. base composition and relative substitution rate). The graphical web interface allows the user to easily browse and sort the results of combined queries. The corresponding multiple sequence alignments and ML trees, inferred using state-of-the art approaches, are available for download both at the nucleotide and amino acid levels. OrthoMaM v12 can be used by researchers interested either in reconstructing the phylogenetic relationships of mammalian taxa or in understanding the evolutionary dynamics of coding sequences in their genomes. OrthoMaM is available for browsing, querying and complete or filtered download at https://orthomam.mbb.cnrs.fr/.

GDPF: a data resource for the distribution of prokaryotic protein families across the global biosphere

Microorganisms encode most of the functions of life on Earth. However, conventional research has primarily focused on specific environments such as humans, soil and oceans, leaving the distribution of functional families throughout the global biosphere poorly comprehended. Here, we present the database of the global distribution of prokaryotic protein families (GDPF, http://bioinfo.qd.sdu.edu.cn/GDPF/), a data resource on the distribution of functional families across the global biosphere. GDPF provides global distribution information for 36 334 protein families, 19 734 superfamilies and 12 089 KEGG (Kyoto Encyclopedia of Genes and Genomes) orthologs from multiple source databases, covering typical environments such as soil, oceans, animals, plants and sediments. Users can browse, search and download the distribution data of each entry in 10 000 global microbial communities, as well as conduct comparative analysis of distribution disparities among multiple entries across various environments. The GDPF data resource contributes to uncovering the geographical distribution patterns, key influencing factors and macroecological principles of microbial functions at a global level, thereby promoting research in Earth ecology and human health.

SPIRE: a Searchable, Planetary-scale mIcrobiome REsource

Meta'omic data on microbial diversity and function accrue exponentially in public repositories, but derived information is often siloed according to data type, study or sampled microbial environment. Here we present SPIRE, a Searchable Planetary-scale mIcrobiome REsource that integrates various consistently processed metagenome-derived microbial data modalities across habitats, geography and phylogeny. SPIRE encompasses 99 146 metagenomic samples from 739 studies covering a wide array of microbial environments and augmented with manually-curated contextual data. Across a total metagenomic assembly of 16 Tbp, SPIRE comprises 35 billion predicted protein sequences and 1.16 million newly constructed metagenome-assembled genomes (MAGs) of medium or high quality. Beyond mapping to the high-quality genome reference provided by proGenomes3 (http://progenomes.embl.de), these novel MAGs form 92 134 novel species-level clusters, the majority of which are unclassified at species level using current tools. SPIRE enables taxonomic profiling of these species clusters via an updated, custom mOTUs database (https://motu-tool.org/) and includes several layers of functional annotation, as well as crosslinks to several (micro-)biological databases. The resource is accessible, searchable and browsable via http://spire.embl.de.

Step-by-Step Bacterial Genome Comparison

Thanks to advancements in genome sequencing and bioinformatics, thousands of bacterial genome sequences are available in public databases. This presents an opportunity to study bacterial diversity in unprecedented detail. This chapter describes a complete bioinformatics workflow for comparative genomics of bacterial genomes, including genome annotation, pangenome reconstruction and visualization, phylogenetic analysis, and identification of sequences of interest such as antimicrobial-resistance genes, virulence factors, and phage sequences. The workflow uses state-of-the-art, open-source tools. The workflow is presented by means of a comparative analysis of Salmonella enterica serovar Typhimurium genomes. The workflow is based on Linux commands and scripts, and result visualization relies on the R environment. The chapter provides a step-by-step protocol that researchers with basic expertise in bioinformatics can easily follow to conduct investigations on their own genome datasets.

Protein-Coding Gene Families in Prokaryote Genome Comparisons

The identification of orthologous genes is relevant for comparative genomics, phylogenetic analysis, and functional annotation. There are many computational tools for the prediction of orthologous groups as well as web-based resources that offer orthology datasets for download and online analysis. This chapter presents a simple and practical guide to the process of orthologous group prediction, using a dataset of 10 prokaryotic proteomes as example. The orthology methods covered are OrthoMCL, COGtriangles, OrthoFinder2, and OMA. The authors compare the number of orthologous groups predicted by these various methods, and present a brief workflow for the functional annotation and reconstruction of phylogenies from inferred single-copy orthologous genes. The chapter also demonstrates how to explore two orthology databases: eggNOG6 and OrthoDB.

Molecular docking and dynamics of a dextranase derived from Penicillium cyclopium CICC-4022

This study aimed to investigate the recognition mechanism of dextranase (PC-Edex) produced by Penicillium cyclopium CICC-4022 on dextran. Whole genome information of P. cyclopium CICC-4022 was obtained through genome sequencing technology. The coding information of PC-Edex was determined based on the annotation of the protein-coding genes using protein databases. The three-dimensional structure of PC-Edex was obtained via homology modelling. The active site and binding free energy between PC-Edex and dextran were calculated by molecular docking and molecular dynamics techniques. The results showed that the total sequence length and GC content of P. cyclopium CICC-4022 were 29,710,801 bp and 47.02 %, respectively. The annotation of protein-encoding genes showed that P. cyclopium CICC-4022 is highly active and has many carbohydrate transport and metabolic functions, and most of its proteases are glycolytic anhydrases. Furthermore, the gene encoding PC-Edex was successfully annotated. Molecular dynamics simulations indicated that van der Waals interaction was the main driving force of interaction. Residues Ile114, Asp115, Tyr332, Lys344, and Gln403 significantly promoted the binding between dextran and PC-Edex. In summary, this study explored the active site catalyzed by PC-Edex based on the binding pattern of PC-Edex and dextran. Therefore, this study provides genomic information on dextranase and data supporting the rational modification and enhancement of PC-Edex.

Functional Genome Analyses Reveal the Molecular Basis of Oil Accumulation in Developing Seeds of Castor Beans

Castor (Ricinus communis L.) seeds produce abundant ricinoleic acid during seed maturation, which is important for plant development and human demands. Ricinoleic acid, as a unique hydroxy fatty acid (HFA), possesses a distinct bond structure that could be used as a substitute for fossil fuels. Here, we identified all homologous genes related to glycolysis, hydroxy fatty acid biosynthesis, and triacylglycerol (TAG) accumulation in castor seeds. Furthermore, we investigated their expression patterns globally during five seed development stages. We characterized a total of 66 genes involved in the glycolysis pathway, with the majority exhibiting higher expression levels during the early stage of castor bean seed development. This metabolic process provided abundant acetyl-CoA for fatty acid (FA) biosynthesis. Subsequently, we identified 82 genes involved in the processes of de novo FA biosynthesis and TAG assembly, with the majority exhibiting high expression levels during the middle or late stages. In addition, we examined the expression patterns of the transcription factors involved in carbohydrate and oil metabolism. For instance, RcMYB73 and RcERF72 exhibited high expression levels during the early stage, whereas RcWRI1, RcABI3, and RcbZIP67 showed relatively higher expression levels during the middle and late stages, indicating their crucial roles in seed development and oil accumulation. Our study suggests that the high HFA production in castor seeds is attributed to the interaction of multiple genes from sugar transportation to lipid droplet packaging. Therefore, this research comprehensively characterizes all the genes related to glycolysis, fatty acid biosynthesis, and triacylglycerol (TAG) accumulation in the castor and provides novel insight into exploring the genetic mechanisms underlying seed oil accumulation in the endosperm of castor beans.

Comparative Transcriptome Analysis Reveals the Molecular Mechanisms of Acetic Acid Reduction by Adding NaHSO3 in Actinobacillus succinogenes GXAS137

Acetic acid (AC) is a major by-product from fermentation processes for producing succinic acid (SA) using Actinobacillus succinogenes. Previous experiments have demonstrated that sodium bisulfate (NaHSO3) can significantly decrease AC production by A. succinogenes GXAS137 during SA fermentation. However, the mechanism of AC reduction is poorly understood. In this study, the transcriptional profiles of the strain were compared through Illumina RNA-seq to identify differentially expressed genes (DEGs). A total of 210 DEGs were identified by expression analysis: 83 and 127 genes up-regulated and down-regulated, respectively, in response to NaHSO3 treatment. The functional annotation analysis of DEGs showed that the genes were mainly involved in carbohydrates, inorganic ions, amino acid transport, metabolism, and energy production and conversion. The mechanisms of AC reduction might be related to two aspects: (i) the lipoic acid synthesis pathway (LipA, LipB) was significantly down-regulated, which blocked the pathway catalyzed by pyruvate dehydrogenase complex to synthesize acetyl-coenzyme A (acetyl-CoA) from pyruvate; (ii) the expression level of the gene encoding bifunctional acetaldehyde-alcohol dehydrogenase was significantly up-regulated, and this effect facilitated the synthesis of ethanol from acetyl-CoA. However, the reaction of NaHSO3 with the intermediate metabolite acetaldehyde blocked the production of ethanol and consumed acetyl-CoA, thereby decreasing AC production. Thus, our study provides new insights into the molecular mechanism of AC decreased underlying the treatment of NaHSO3 and will deepen the understanding of the complex regulatory mechanisms of A. succinogenes.

Antimicrobial resistance and biotechnological potential of plastic-associated bacteria isolated from an urban estuary

Plastics have quickly become one of the major pollutants in aquatic environments worldwide and solving the plastic pollution crisis is considered a central goal of modern society. In this study, 10 different plastic samples, including high- and low-density polyethylene and polypropylene, were collected from a deeply polluted urban estuary in Brazil. By employing different isolation and analysis approaches to investigate plastic-associated bacteria, a predominance of potentially pathogenic bacteria such as Acinetobacter, Aeromonas, and Vibrio was observed throughout all plastic samples. Bacteria typically found in the aquatic environment harboured clinically relevant genes encoding resistance to carbapenems (blaKPC ) and colistin (such as mcr-3 and mcr-4), along with genetic determinants associated with potentially active gene mobilization. Whole genome sequencing and annotation of three plastic-associated Vibrio strains further demonstrated the carriage of mobile genetic elements and antimicrobial resistance and virulence genes. On the other hand, bacteria isolated from the same samples were also able to produce esterases, lipases, and bioemulsifiers, thus highlighting that the plastisphere could also be of special interest from a biotechnological perspective.

FAM171B stabilizes vimentin and enhances CCL2-mediated TAM infiltration to promote bladder cancer progression

BACKGROUND

Invasion and metastasis are the main causes of unfavourable prognosis in patients diagnosed with bladder cancer. The efficacy of immunotherapy in bladder cancer remains suboptimal due to the presence of an immunosuppressive microenvironment. The novel protein family with sequence similarity 171B (FAM171B) has been identified, but its precise role and mechanism in bladder cancer remain unclear.

METHODS

In this study, we conducted an analysis to investigate the associations between FAM171B expression and the prognosis and clinicopathological stage of bladder cancer. To this end, we utilized RNA sequencing data from the TCGA and GEO databases, as well as tumor tissue specimens obtained from our clinical centre. RNA sequencing analysis allowed us to examine the biological function of FAM171B at the transcriptional level in bladder cancer cells. Additionally, we used immunoprecipitation and mass spectrometry to identify the protein that interacts with FAM171B in bladder cancer cells. The effects of FAM171B on modulating tumor-associated macrophages (TAMs) and vimentin-mediated tumor progression, as well as the underlying mechanisms, were clarified by phalloidin staining, immunofluorescence staining, ELISA, RNA immunoprecipitation, flow cytometry and a bladder cancer graft model.

RESULTS

FAM171B expression exhibits strong positive correlation with poor survival outcomes and advanced clinicopathological stages in patients with bladder cancer. FAM171B significantly promoted bladder cancer growth and metastasis, accompanied by TAM accumulation in the microenvironment, in vivo and in vitro. Through studies of the molecular mechanism, we found that FAM171B contributes to tumor progression by stabilizing vimentin in the cytoplasm. Additionally, our research revealed that FAM171B enhances the splicing of CCL2 mRNA by interacting with heterogeneous nuclear ribonucleoprotein U (HNRNPU), ultimately leading to increased recruitment and M2 polarization of TAMs.

CONCLUSIONS

In this study, we identified FAM171B as a potent factor that promotes the progression of bladder cancer. These findings establish a solid theoretical foundation for considering FAM171B as a potential diagnostic and therapeutic biomarker for bladder cancer.

A High-Quality Reference Genome Assembly of Prinsepia uniflora (Rosaceae)

This study introduces a meticulously constructed genome assembly at the chromosome level for the Rosaceae family species Prinsepia uniflora, a traditional Chinese medicinal herb. The final assembly encompasses 1272.71 megabases (Mb) distributed across 16 pseudochromosomes, boasting contig and super-scaffold N50 values of 2.77 and 79.32 Mb, respectively. Annotated within this genome is a substantial 875.99 Mb of repetitive sequences, with transposable elements occupying 777.28 Mb, constituting 61.07% of the entire genome. Our predictive efforts identified 49,261 protein-coding genes within the repeat-masked assembly, with 45,256 (91.87%) having functional annotations, 5127 (10.41%) demonstrating tandem duplication, and 2373 (4.82%) classified as transcription factor genes. Additionally, our investigation unveiled 3080 non-coding RNAs spanning 0.51 Mb of the genome sequences. According to our evolutionary study, P. uniflora underwent recent whole-genome duplication following its separation from Prunus salicina. The presented reference-level genome assembly and annotation for P. uniflora will significantly facilitate the in-depth exploration of genomic information pertaining to this species, offering substantial utility in comparative genomics and evolutionary analyses involving Rosaceae species.

oggmap: a Python package to extract gene ages per orthogroup and link them with single-cell RNA data

SUMMARY

For model species, single-cell RNA-based cell atlases are available. A good cell atlas includes all major stages in a species' ontogeny, and soon, they will be standard even for nonmodel species. Here, we propose a Python package called oggmap, which allows for the easy extraction of an orthomap (gene ages per orthogroup) for any given query species from OrthoFinder and other gene family data resources, like homologous groups from eggNOG or PLAZA. oggmap provides extracted gene ages for more than thousand eukaryotic species which can be further used to calculate gene age-weighted expression data from scRNA sequencing objects using the Python Scanpy toolkit. Not limited to one transcriptome evolutionary index, oggmap can visualize the individual gene category (e.g. age class, nucleotide diversity bin) and their corresponding expression profiles to investigate scRNA-based cell type assignments in an evolutionary context.

AVAILABILITY AND IMPLEMENTATION

oggmap source code is available at https://github.com/kullrich/oggmap, documentation is available at https://oggmap.readthedocs.io/en/latest/. oggmap can be installed via PyPi or directly used via a docker container.

Well-hidden methanogenesis in deep, organic-rich sediments of Guaymas Basin

Deep marine sediments (>1mbsf) harbor ~26% of microbial biomass and are the largest reservoir of methane on Earth. Yet, the deep subsurface biosphere and controls on its contribution to methane production remain underexplored. Here, we use a multidisciplinary approach to examine methanogenesis in sediments (down to 295 mbsf) from sites with varying degrees of thermal alteration (none, past, current) at Guaymas Basin (Gulf of California) for the first time. Traditional (13C/12C and D/H) and multiply substituted (13CH3D and 12CH2D2) methane isotope measurements reveal significant proportions of microbial methane at all sites, with the largest signal at the site with past alteration. With depth, relative microbial methane decreases at differing rates between sites. Gibbs energy calculations confirm methanogenesis is exergonic in Guaymas sediments, with methylotrophic pathways consistently yielding more energy than the canonical hydrogenotrophic and acetoclastic pathways. Yet, metagenomic sequencing and cultivation attempts indicate that methanogens are present in low abundance. We find only one methyl-coenzyme M (mcrA) sequence within the entire sequencing dataset. Also, we identify a wide diversity of methyltransferases (mtaB, mttB), but only a few sequences phylogenetically cluster with methylotrophic methanogens. Our results suggest that the microbial methane in the Guaymas subsurface was produced over geologic time by relatively small methanogen populations, which have been variably influenced by thermal sediment alteration. Higher resolution metagenomic sampling may clarify the modern methanogen community. This study highlights the importance of using a multidisciplinary approach to capture microbial influences in dynamic, deep subsurface settings like Guaymas Basin.

Robiginitalea aestuariiviva sp. nov. isolated from sediment of tidal flat located in Zhejiang, PR China

A Gram-stain-negative, strictly aerobic and rod- to coccoid-shaped bacterium, designated as strain M366T, was isolated from coastal sediment of Jiaoshanjiao, Zhejiang Province, PR China (121°54' E 29 °38' N). The draft genome of strain M366T was 3 225 479 bp long (with 55.6 mol% G+C content) and assembled into four contigs. The N50 value was 563 270 bp and the genomic completeness and contamination were estimated to be 99.34 and 0.05 %, respectively. Colonies of strain M366T were yellow-orange, 1 mm in diameter, round, opaque, smooth and convex after incubation on marine agar at 30 °C for 3 days. Cells were catalase-positive but oxidase-negative. Strain M366T was observed to grow at 20-40 °C (optimum, 30 °C), pH 5.5-9.0 (optimum, pH 6.5-7.0) and with 0.5-8.0 % (w/v) NaCl (optimum, 2.5 %). Strain M366T shown highest 16S rRNA gene sequence similarity of 98.1 % to Robiginitalea sediminis O458T, 95.6-95.9 % to other type strains of the genus Robiginitalea and below 93 % to other genera. The average nucleotide identity and digital DNA-DNA hybridization values between strain M366T and its closely related Robiginitalea species were 71.1-75.9 % and 17.5-19.0 %. Menaquinone-6 was the only respiratory quinone. The major fatty acids (>10 %) were iso-C15 : 0, iso-C17 : 0 3-OH and summed feature 1 (iso-C15 : 1 h and/or C13 : 0 3-OH). The main polar lipids included phosphatidylethanolamine, two unidentified phospholipid, two unidentified aminophospholipid, one unidentified glycolipid and five unidentified lipids. According to the above results, Robiginitalea aestuariiviva sp. nov. is proposed and the type strain is M366T (=KCTC 92866T=MCCC 1K04524T=CGMCC 1.61708T).

ORFanID: A web-based search engine for the discovery and identification of orphan and taxonomically restricted genes

With the numerous genomes sequenced today, it has been revealed that a noteworthy percentage of genes in a given taxon of organisms in the phylogenetic tree of life do not have orthologous sequences in other taxa. These sequences are commonly referred to as "orphans" or "ORFans" if found as single occurrences in a single species or as "taxonomically restricted genes" (TRGs) when found at higher taxonomic levels. Quantitative and collective studies of these genes are necessary for understanding their biological origins. However, the current software for identifying orphan genes is limited in its functionality, database search range, and very complex algorithmically. Thus, researchers studying orphan genes must harvest their data from many disparate sources. ORFanID is a graphical web-based search engine that facilitates the efficient identification of both orphan genes and TRGs at all taxonomic levels, from DNA or amino acid sequences in the NCBI database cluster and other large bioinformatics repositories. The software allows users to identify genes that are unique to any taxonomic rank, from species to domain, using NCBI systematic classifiers. It provides control over NCBI database search parameters, and the results are presented in a spreadsheet as well as a graphical display. The tables in the software are sortable, and results can be filtered using the fuzzy search functionality. The visual presentation can be expanded and collapsed by the taxonomic tree to its various branches. Example results from searches on five species and gene expression data from specific orphan genes are provided in the Supplementary Information.

Genomic Characteristics and Comparative Genomics Analysis of the Endophytic Fungus Paraphoma chrysanthemicola DS-84 Isolated from Codonopsis pilosula Root

Paraphoma chrysanthemicola is a newly identified endophytic fungus. The focus of most studies on P. chrysanthemicola has been on its isolation, identification and effects on plants. However, the limited genomic information is a barrier to further research. Therefore, in addition to studying the morphological and physiological characteristics of P. chrysanthemicola, we sequenced its genome and compared it with that of Paraphoma sp. The results showed that sucrose, peptone and calcium phosphate were suitable sources of carbon, nitrogen and phosphorus for this strain. The activities of amylase, cellulase, chitosanase, lipase and alkaline protease were also detected. Sequencing analysis revealed that the genome of P. chrysanthemicola was 44.1 Mb, with a scaffold N50 of 36.1 Mb and 37,077 protein-coding genes. Gene Ontology (GO) annotation showed that mannose-modified glycosylation was predominant in monosaccharide utilisation. The percentage of glycoside hydrolase (GH) modules was the highest in the carbohydrate-active enzymes database (CAZy) analysis. Secondary metabolite-associated gene cluster analysis identified melanin, dimethylcoprogen and phyllostictine A biosynthetic gene clusters (>60% similarity). The results indicated that P. chrysanthemicola had a mannose preference in monosaccharide utilisation and that melanin, dimethylcoprogen and phyllostictine A were important secondary metabolites for P. chrysanthemicola as an endophytic fungus.

Genome Variation Map of Domestic Qinghai-Tibet Plateau Yaks by SLAF-Seq Reveals Genetic Footprint during Artificial Selection

The yak (Bos grunniens) was domesticated in the high-altitude QTP. Research about their genetic diversity and population structure is limited. In this study, we resequenced the genome of 494 domestic yaks using Specific-Locus Amplified Fragment Sequencing (SLAF-seq). The survey was conducted on six populations sampled from isolated locations in China in order to analyze their structure and genetic diversity. These six domestic populations were clearly grouped into two independent clusters, with Jinchuan, Changtai, and Jiulong showing a tight genetic relationship with the wild yak. Nerve development pathways were enriched with GO enrichment analysis of 334 domesticated genes. Major genomic regions associated with the differentiation of domestic yaks were detected. These findings provide preliminary information on the yak genome variability, useful to understand the genomic characteristics of different populations in QTP.

Wide-scale identification of novel/eliminated genes responsible for evolutionary transformations

BACKGROUND

It is generally accepted that most evolutionary transformations at the phenotype level are associated either with rearrangements of genomic regulatory elements, which control the activity of gene networks, or with changes in the amino acid contents of proteins. Recently, evidence has accumulated that significant evolutionary transformations could also be associated with the loss/emergence of whole genes. The targeted identification of such genes is a challenging problem for both bioinformatics and evo-devo research.

RESULTS

To solve this problem we propose the WINEGRET method, named after the first letters of the title. Its main idea is to search for genes that satisfy two requirements: first, the desired genes were lost/emerged at the same evolutionary stage at which the phenotypic trait of interest was lost/emerged, and second, the expression of these genes changes significantly during the development of the trait of interest in the model organism. To verify the first requirement, we do not use existing databases of orthologs, but rely purely on gene homology and local synteny by using some novel quickly computable conditions. Genes satisfying the second requirement are found by deep RNA sequencing. As a proof of principle, we used our method to find genes absent in extant amniotes (reptiles, birds, mammals) but present in anamniotes (fish and amphibians), in which these genes are involved in the regeneration of large body appendages. As a result, 57 genes were identified. For three of them, c-c motif chemokine 4, eotaxin-like, and a previously unknown gene called here sod4, essential roles for tail regeneration were demonstrated. Noteworthy, we established that the latter gene belongs to a novel family of Cu/Zn-superoxide dismutases lost by amniotes, SOD4.

CONCLUSIONS

We present a method for targeted identification of genes whose loss/emergence in evolution could be associated with the loss/emergence of a phenotypic trait of interest. In a proof-of-principle study, we identified genes absent in amniotes that participate in body appendage regeneration in anamniotes. Our method provides a wide range of opportunities for studying the relationship between the loss/emergence of phenotypic traits and the loss/emergence of specific genes in evolution.

BRAT1 associates with INTS11/INTS9 heterodimer to regulate key neurodevelopmental genes

Integrator is a multi-subunits protein complex involved in regulation of gene expression. Several Integrator subunits have been found to be mutated in human neurodevelopmental disorders, suggesting a key role for the complex in the development of nervous system. BRAT1 is similarly linked with neurodegenerative diseases and neurodevelopmental disorders such as rigidity and multifocal-seizure syndrome. Here, we show that INTS11 and INTS9 subunits of Integrator complex interact with BRAT1 and form a trimeric complex in human HEK293T cells as well as in pluripotent human embryonal carcinoma cell line (NT2). We find that BRAT1 depletion disrupts the differentiation of NT2 cells into astrocytes and neural cells. Loss of BRAT1 results in inability to activate many neuronal genes that are targets of REST, a neuronal silencer. We identified BRAT1 and INTS11 co-occupying the promoter region of these genes and pinpoint a role for BRAT1 in recruiting INTS11 to their promoters. Disease-causing mutations in BRAT1 diminish its association with INTS11/INTS9, linking the manifestation of disease phenotypes with a defect in transcriptional activation of key neuronal genes by BRAT1/INTS11/INTS9 complex.

Highlights

Integrator subunits INTS9 and INTS11 tightly interact with BRAT1 Depletion of BRAT1 causes a dramatic delay in human neural differentiation BRAT1 and INTS11 module targets the promoters of neural marker genes and co-regulates their expression. The recruitment of INTS11 to these sites is BRAT1-dependent. Pathogenic E522K mutation in BRAT1 disrupts its interaction with INTS11/INTS9 heterodimer.