Hidden Markov model speed heuristic and iterative HMM search procedure

Structural modeling of ion channels using AlphaFold2, RoseTTAFold2, and ESMFold

Ion channels play key roles in human physiology and are important targets in drug discovery. The atomic-scale structures of ion channels provide invaluable insights into a fundamental understanding of the molecular mechanisms of channel gating and modulation. Recent breakthroughs in deep learning-based computational methods, such as AlphaFold, RoseTTAFold, and ESMFold have transformed research in protein structure prediction and design. We review the application of AlphaFold, RoseTTAFold, and ESMFold to structural modeling of ion channels using representative voltage-gated ion channels, including human voltage-gated sodium (NaV) channel - NaV1.8, human voltage-gated calcium (CaV) channel - CaV1.1, and human voltage-gated potassium (KV) channel - KV1.3. We compared AlphaFold, RoseTTAFold, and ESMFold structural models of NaV1.8, CaV1.1, and KV1.3 with corresponding cryo-EM structures to assess details of their similarities and differences. Our findings shed light on the strengths and limitations of the current state-of-the-art deep learning-based computational methods for modeling ion channel structures, offering valuable insights to guide their future applications for ion channel research.

Mining proteomes for zinc finger persulfidation

Hydrogen sulfide (H2S) is an endogenous gasotransmitter that signals via persulfidation. There is evidence that the cysteine residues of certain zinc finger (ZF) proteins, a common type of cysteine rich protein, are modified to persulfides by H2S. To determine how frequently ZF persulfidation occurs in cells and identify the types of ZFs that are persulfidated, persulfide specific proteomics data were evaluated. 22 datasets from 16 studies were analyzed via a meta-analysis approach. Persulfidated ZFs were identified in a range of eukaryotic species, including Homo sapiens, Mus musculus, Rattus norvegicus, Arabidopsis thaliana, and Emiliania huxley (single-celled phytoplankton). The types of ZFs identified for each species encompassed all three common ZF ligand sets (4-cysteine, 3-cysteine-1-histidine, and 2-cysteine-2-hisitidine), indicating that persulfidation of ZFs is broad. Overlap analysis between different species identified several common ZFs. GO and KEGG analysis identified pathway enrichment for ubiquitin-dependent protein catabolic process and viral carcinogenesis. These collective findings support ZF persulfidation as a wide-ranging PTM that impacts all classes of ZFs.

PII-like signaling proteins: a new paradigm in orchestrating cellular homeostasis

Members of the PII superfamily are versatile, multitasking signaling proteins ubiquitously found in all domains of life. They adeptly monitor and synchronize the cell's carbon, nitrogen, energy, redox, and diurnal states, primarily by binding interdependently to adenyl-nucleotides, including charged nucleotides (ATP, ADP, and AMP) and second messengers such as cyclic adenosine monophosphate (cAMP), cyclic di-adenosine monophosphate (c-di-AMP), and S-adenosylmethionine-AMP (SAM-AMP). These proteins also undergo a variety of posttranslational modifications, such as phosphorylation, adenylation, uridylation, carboxylation, and disulfide bond formation, which further provide cues on the metabolic state of the cell. Serving as precise metabolic sensors, PII superfamily proteins transmit this information to diverse cellular targets, establishing dynamic regulatory assemblies that fine-tune cellular homeostasis. Recently discovered, PII-like proteins are emerging families of signaling proteins that, while related to canonical PII proteins, have evolved to fulfill a diverse range of cellular functions, many of which remain elusive. In this review, we focus on the evolution of PII-like proteins and summarize the molecular mechanisms governing the assembly dynamics of PII complexes, with a special emphasis on the PII-like protein SbtB.

A molecular switch controls assembly of bacterial focal adhesions

Cell motility universally relies on spatial regulation of focal adhesion complexes (FAs) connecting the substrate to cellular motors. In bacterial FAs, the Adventurous gliding motility machinery (Agl-Glt) assembles at the leading cell pole following a Mutual gliding-motility protein (MglA)-guanosine 5'-triphosphate (GTP) gradient along the cell axis. Here, we show that GltJ, a machinery membrane protein, contains cytosolic motifs binding MglA-GTP and AglZ and recruiting the MreB cytoskeleton to initiate movement toward the lagging cell pole. In addition, MglA-GTP binding triggers a conformational shift in an adjacent GltJ zinc-finger domain, facilitating MglB recruitment near the lagging pole. This prompts GTP hydrolysis by MglA, leading to complex disassembly. The GltJ switch thus serves as a sensor for the MglA-GTP gradient, controlling FA activity spatially.

Phylogenomics reveals the evolutionary origins of lichenization in chlorophyte algae

Mutualistic symbioses have contributed to major transitions in the evolution of life. Here, we investigate the evolutionary history and the molecular innovations at the origin of lichens, which are a symbiosis established between fungi and green algae or cyanobacteria. We de novo sequence the genomes or transcriptomes of 12 lichen algal symbiont (LAS) and closely related non-symbiotic algae (NSA) to improve the genomic coverage of Chlorophyte algae. We then perform ancestral state reconstruction and comparative phylogenomics. We identify at least three independent gains of the ability to engage in the lichen symbiosis, one in Trebouxiophyceae and two in Ulvophyceae, confirming the convergent evolution of the lichen symbioses. A carbohydrate-active enzyme from the glycoside hydrolase 8 (GH8) family was identified as a top candidate for the molecular-mechanism underlying lichen symbiosis in Trebouxiophyceae. This GH8 was acquired in lichenizing Trebouxiophyceae by horizontal gene transfer, concomitantly with the ability to associate with lichens fungal symbionts (LFS) and is able to degrade polysaccharides found in the cell wall of LFS. These findings indicate that a combination of gene family expansion and horizontal gene transfer provided the basis for lichenization to evolve in chlorophyte algae.

Novel enzymes for biodegradation of polycyclic aromatic hydrocarbons identified by metagenomics and functional analysis in short-term soil microcosm experiments

Polycyclic aromatic hydrocarbons (PAHs) are highly toxic, carcinogenic substances. On soils contaminated with PAHs, crop cultivation, animal husbandry and even the survival of microflora in the soil are greatly perturbed, depending on the degree of contamination. Most microorganisms cannot tolerate PAH-contaminated soils, however, some microbial strains can adapt to these harsh conditions and survive on contaminated soils. Analysis of the metagenomes of contaminated environmental samples may lead to discovery of PAH-degrading enzymes suitable for green biotechnology methodologies ranging from biocatalysis to pollution control. In the present study, our goal was to apply a metagenomic data search to identify efficient novel enzymes in remediation of PAH-contaminated soils. The metagenomic hits were further analyzed using a set of bioinformatics tools to select protein sequences predicted to encode well-folded soluble enzymes. Three novel enzymes (two dioxygenases and one peroxidase) were cloned and used in soil remediation microcosms experiments. The experimental design of the present study aimed at evaluating the effectiveness of the novel enzymes on short-term PAH degradation in the soil microcosmos model. The novel enzymes were found to be efficient for degradation of naphthalene and phenanthrene. Adding the inorganic oxidant CaO2 further increased the degrading potential of the novel enzymes for anthracene and pyrene. We conclude that metagenome mining paired with bioinformatic predictions, structural modelling and functional assays constitutes a powerful approach towards novel enzymes for soil remediation.

Structural and mechanistic insights into Quinolone Synthase to address its functional promiscuity

Quinolone synthase from Aegle marmelos (AmQNS) is a type III polyketide synthase that yields therapeutically effective quinolone and acridone compounds. Addressing the structural and molecular underpinnings of AmQNS and its substrate interaction in terms of its high selectivity and specificity can aid in the development of numerous novel compounds. This paper presents a high-resolution AmQNS crystal structure and explains its mechanistic role in synthetic selectivity. Additionally, we provide a model framework to comprehend structural constraints on ketide insertion and postulate that AmQNS's steric and electrostatic selectivity plays a role in its ability to bind to various core substrates, resulting in its synthetic diversity. AmQNS prefers quinolone synthesis and can accommodate large substrates because of its wide active site entrance. However, our research suggests that acridone is exclusively synthesized in the presence of high malonyl-CoA concentrations. Potential implications of functionally relevant residue mutations were also investigated, which will assist in harnessing the benefits of mutations for targeted polyketide production. The pharmaceutical industry stands to gain from these findings as they expand the pool of potential drug candidates, and these methodologies can also be applied to additional promising enzymes.

Structural Reconstruction of E. coli Ubi Metabolon Using an AlphaFold2-Based Computational Framework

Ubiquinone (UQ) is a redox polyisoprenoid lipid found in the membranes of bacteria and eukaryotes that has important roles, notably one in respiratory metabolism, which sustains cellular bioenergetics. In Escherichia coli, several steps of the UQ biosynthesis take place in the cytosol. To perform these reactions, a supramolecular assembly called Ubi metabolon is involved. This latter is composed of seven proteins (UbiE, UbiG, UbiF, UbiH, UbiI, UbiJ, and UbiK), and its structural organization is unknown as well as its protein stoichiometry. In this study, a computational framework has been designed to predict the structure of this macromolecular assembly. In several successive steps, we explored the possible protein interactions as well as the protein stoichiometry, to finally obtain a structural organization of the complex. The use of AlphaFold2-based methods combined with evolutionary information enabled us to predict several models whose quality and confidence were further analyzed using different metrics and scores. Our work led to the identification of a "core assembly" that will guide functional and structural characterization of the Ubi metabolon.

MFTrans: A multi-feature transformer network for protein secondary structure prediction

In the rapidly evolving field of computational biology, accurate prediction of protein secondary structures is crucial for understanding protein functions, facilitating drug discovery, and advancing disease diagnostics. In this paper, we propose MFTrans, a deep learning-based multi-feature fusion network aimed at enhancing the precision and efficiency of Protein Secondary Structure Prediction (PSSP). This model employs a Multiple Sequence Alignment (MSA) Transformer in combination with a multi-view deep learning architecture to effectively capture both global and local features of protein sequences. MFTrans integrates diverse features generated by protein sequences, including MSA, sequence information, evolutionary information, and hidden state information, using a multi-feature fusion strategy. The MSA Transformer is utilized to interleave row and column attention across the input MSA, while a Transformer encoder and decoder are introduced to enhance the extracted high-level features. A hybrid network architecture, combining a convolutional neural network with a bidirectional Gated Recurrent Unit (BiGRU) network, is used to further extract high-level features after feature fusion. In independent tests, our experimental results show that MFTrans has superior generalization ability, outperforming other state-of-the-art PSSP models by 3 % on average on public benchmarks including CASP12, CASP13, CASP14, TEST2016, TEST2018, and CB513. Case studies further highlight its advanced performance in predicting mutation sites. MFTrans contributes significantly to the protein science field, opening new avenues for drug discovery, disease diagnosis, and protein.

Sphingobacterium oryzagri sp. nov., isolated from rice paddy soil

An aerobic, Gram-stain-negative and short rod-shaped bacterial strain, designated M6-31T, was isolated from rice paddy soil sampled in Miryang, Republic of Korea. Growth was observed at 4-35 °C (optimum, 28 °C), pH 6.0-9.0 (optimum, pH 7.0-8.0) and in the presence of 0-4 % (w/v) NaCl (optimum, 0 % w/v). Phylogenetic analysis based on 16S rRNA gene sequences grouped strain M6-31T with Sphingobacterium bambusae IBFC2009T, Sphingobacterium griseoflavum SCU-B140T and Sphingobacterium solani MLS-26-JM13-11T in the same clade, with the 16S rRNA gene sequence similarities ranging from 95.8 to 96.6 %. A genome-based phylogenetic tree reconstructed by using all publicly available Sphingobacterium genomes placed strain M6-31T with S. bambusae KACC 22910T, 'Sphingobacterium deserti' ACCC 05744T, S. griseoflavum CGMCC 1.12966T and Sphingobacterium paludis CGMCC 1.12801T. Orthologous average nucleotide identity and digital DNA-DNA hybridization values between strain M6-31T and its closely related strains were lower than 74.6 and 22.0 %, respectively. The respiratory quinone was menaquinone-7, and the major polar lipid was phosphatidylethanolamine. The major fatty acids (>10 %) were C15 : 0 iso, C17 : 0 iso 3OH and summed feature 3. The phenotypic, chemotaxonomic and genotypic data obtained in this study showed that strain M6-31T represents a novel species of the genus Sphingobacterium, for which the name Sphingobacterium oryzagri sp. nov. (type strain M6-31T=KACC 22765T=JCM 35893T) is proposed.

Discovering a mitochondrion-localized BAHD acyltransferase involved in calystegine biosynthesis and engineering the production of 3β-tigloyloxytropane

Solanaceous plants produce tropane alkaloids (TAs) via esterification of 3α- and 3β-tropanol. Although littorine synthase is revealed to be responsible for 3α-tropanol esterification that leads to hyoscyamine biosynthesis, the genes associated with 3β-tropanol esterification are unknown. Here, we report that a BAHD acyltransferase from Atropa belladonna, 3β-tigloyloxytropane synthase (TS), catalyzes 3β-tropanol and tigloyl-CoA to form 3β-tigloyloxytropane, the key intermediate in calystegine biosynthesis and a potential drug for treating neurodegenerative disease. Unlike other cytosolic-localized BAHD acyltransferases, TS is localized to mitochondria. The catalytic mechanism of TS is revealed through molecular docking and site-directed mutagenesis. Subsequently, 3β-tigloyloxytropane is synthesized in tobacco. A bacterial CoA ligase (PcICS) is found to synthesize tigloyl-CoA, an acyl donor for 3β-tigloyloxytropane biosynthesis. By expressing TS mutant and PcICS, engineered Escherichia coli synthesizes 3β-tigloyloxytropane from tiglic acid and 3β-tropanol. This study helps to characterize the enzymology and chemodiversity of TAs and provides an approach for producing 3β-tigloyloxytropane.

Ultrastructural insights into cellular organization, energy storage and ribosomal dynamics of an ammonia-oxidizing archaeon from oligotrophic oceans

Introduction

Nitrososphaeria, formerly known as Thaumarchaeota, constitute a diverse and widespread group of ammonia-oxidizing archaea (AOA) inhabiting ubiquitously in marine and terrestrial environments, playing a pivotal role in global nitrogen cycling. Despite their importance in Earth's ecosystems, the cellular organization of AOA remains largely unexplored, leading to a significant unanswered question of how the machinery of these organisms underpins metabolic functions.

Methods

In this study, we combined spherical-chromatic-aberration-corrected cryo-electron tomography (cryo-ET), scanning transmission electron microscopy (STEM), and energy dispersive X-ray spectroscopy (EDS) to unveil the cellular organization and elemental composition of Nitrosopumilus maritimus SCM1, a representative member of marine Nitrososphaeria.

Results and Discussion

Our tomograms show the native ultrastructural morphology of SCM1 and one to several dense storage granules in the cytoplasm. STEM-EDS analysis identifies two types of storage granules: one type is possibly composed of polyphosphate and the other polyhydroxyalkanoate. With precise measurements using cryo-ET, we observed low quantity and density of ribosomes in SCM1 cells, which are in alignment with the documented slow growth of AOA in laboratory cultures. Collectively, these findings provide visual evidence supporting the resilience of AOA in the vast oligotrophic marine environment.

Freshwater genome-reduced bacteria exhibit pervasive episodes of adaptive stasis

The emergence of bacterial species is rooted in their inherent potential for continuous evolution and adaptation to an ever-changing ecological landscape. The adaptive capacity of most species frequently resides within the repertoire of genes encoding the secreted proteome (SP), as it serves as a primary interface used to regulate survival/reproduction strategies. Here, by applying evolutionary genomics approaches to metagenomics data, we show that abundant freshwater bacteria exhibit biphasic adaptation states linked to the eco-evolutionary processes governing their genome sizes. While species with average to large genomes adhere to the dominant paradigm of evolution through niche adaptation by reducing the evolutionary pressure on their SPs (via the augmentation of functionally redundant genes that buffer mutational fitness loss) and increasing the phylogenetic distance of recombination events, most of the genome-reduced species exhibit a nonconforming state. In contrast, their SPs reflect a combination of low functional redundancy and high selection pressure, resulting in significantly higher levels of conservation and invariance. Our findings indicate that although niche adaptation is the principal mechanism driving speciation, freshwater genome-reduced bacteria often experience extended periods of adaptive stasis. Understanding the adaptive state of microbial species will lead to a better comprehension of their spatiotemporal dynamics, biogeography, and resilience to global change.

Can Protein Structure Prediction Methods Capture Alternative Conformations of Membrane Transporters?

Understanding the conformational dynamics of proteins, such as the inward-facing (IF) and outward-facing (OF) transition observed in transporters, is vital for elucidating their functional mechanisms. Despite significant advances in protein structure prediction (PSP) over the past three decades, most efforts have been focused on single-state prediction, leaving multistate or alternative conformation prediction (ACP) relatively unexplored. This discrepancy has led to the development of highly accurate PSP methods such as AlphaFold, yet their capabilities for ACP remain limited. To investigate the performance of current PSP methods in ACP, we curated a data set, named IOMemP, consisting of 32 experimentally determined high-resolution IF and OF structures of 16 membrane proteins with substantial conformational changes. We benchmarked 12 representative PSP methods, along with two recent multistate methods based on AlphaFold, against this data set. Our findings reveal a remarkably consistent preference for specific states across various PSP methods. We elucidated how coevolution information in MSAs influences state preference. Moreover, we showed that AlphaFold, when excluding coevolution information, estimated similar energies between the experimental IF and OF conformations, indicating that the energy model learned by AlphaFold is not biased toward any particular state. Our IOMemP data set and benchmark results are anticipated to advance the development of robust ACP methods.

Whole-genome sequencing of four culturable endophytic bacteria from German hardneck garlic cloves (Allium sativum L.)

We present the whole-genome sequence of four bacterial endophytes associated with German hardneck garlic cloves (Allium sativum L.). Among them, Agrobacterium fabrum and Pantoea agglomerans are associated with plant protection, while Rahnella perminowiae and Stenotrophomonas lactitubi are pathogens. These data will facilitate the identification of genes to improve garlic.

A map of the rubisco biochemical landscape

Rubisco is the primary CO2 fixing enzyme of the biosphere yet has slow kinetics. The roles of evolution and chemical mechanism in constraining the sequence landscape of rubisco remain debated. In order to map sequence to function, we developed a massively parallel assay for rubisco using an engineered E. coli where enzyme function is coupled to growth. By assaying >99% of single amino acid mutants across CO2 concentrations, we inferred enzyme velocity and CO2 affinity for thousands of substitutions. We identified many highly conserved positions that tolerate mutation and rare mutations that improve CO2 affinity. These data suggest that non-trivial kinetic improvements are readily accessible and provide a comprehensive sequence-to-function mapping for enzyme engineering efforts.

Protocol for isolating single species of bacteria with swarming ability from human feces

Understanding the specific movements of bacteria isolated from human feces can serve as a novel diagnostic and therapeutic tool for inflammatory bowel disease. Here, we present a protocol for a microbial swarming assay and to isolate the bacteria responsible for swarming activity. We describe steps for identifying bacteria using MALDI-TOF mass spectrometry and whole-genome sequencing. We then detail procedures for validating findings by observing the same swarming phenotype upon reperforming the swarming assay. For complete details on the use and execution of this protocol, please refer to De et al.1.

Coding genomes with gapped pattern graph convolutional network

MOTIVATION

Genome sequencing technologies reveal a huge amount of genomic sequences. Neural network-based methods can be prime candidates for retrieving insights from these sequences because of their applicability to large and diverse datasets. However, the highly variable lengths of genome sequences severely impair the presentation of sequences as input to the neural network. Genetic variations further complicate tasks that involve sequence comparison or alignment.

RESULTS

Inspired by the theory and applications of "spaced seeds," we propose a graph representation of genome sequences called "gapped pattern graph." These graphs can be transformed through a Graph Convolutional Network to form lower-dimensional embeddings for downstream tasks. On the basis of the gapped pattern graphs, we implemented a neural network model and demonstrated its performance on diverse tasks involving microbe and mammalian genome data. Our method consistently outperformed all the other state-of-the-art methods across various metrics on all tasks, especially for the sequences with limited homology to the training data. In addition, our model was able to identify distinct gapped pattern signatures from the sequences.

AVAILABILITY AND IMPLEMENTATION

The framework is available at https://github.com/deepomicslab/GCNFrame.

Single-sequence protein structure prediction by integrating protein language models

Protein structure prediction has been greatly improved by deep learning in the past few years. However, the most successful methods rely on multiple sequence alignment (MSA) of the sequence homologs of the protein under prediction. In nature, a protein folds in the absence of its sequence homologs and thus, a MSA-free structure prediction method is desired. Here, we develop a single-sequence-based protein structure prediction method RaptorX-Single by integrating several protein language models and a structure generation module and then study its advantage over MSA-based methods. Our experimental results indicate that in addition to running much faster than MSA-based methods such as AlphaFold2, RaptorX-Single outperforms AlphaFold2 and other MSA-free methods in predicting the structure of antibodies (after fine-tuning on antibody data), proteins of very few sequence homologs, and single mutation effects. By comparing different protein language models, our results show that not only the scale but also the training data of protein language models will impact the performance. RaptorX-Single also compares favorably to MSA-based AlphaFold2 when the protein under prediction has a large number of sequence homologs.

Analysis and culturing of the prototypic crAssphage reveals a phage-plasmid lifestyle

The prototypic crAssphage (Carjivirus communis) is one of the most abundant, prevalent, and persistent gut bacteriophages, yet it remains uncultured and its lifestyle uncharacterized. For the last decade, crAssphage has escaped plaque-dependent culturing efforts, leading us to investigate alternative lifestyles that might explain its widespread success. Through genomic analyses and culturing, we find that crAssphage uses a phage-plasmid lifestyle to persist extrachromosomally. Plasmid-related genes are more highly expressed than those implicated in phage maintenance. Leveraging this finding, we use a plaque-free culturing approach to measure crAssphage replication in culture with Phocaeicola vulgatus, Phocaeicola dorei, and Bacteroides stercoris, revealing a broad host range. We demonstrate that crAssphage persists with its hosts in culture without causing major cell lysis events or integrating into host chromosomes. The ability to switch between phage and plasmid lifestyles within a wide range of hosts contributes to the prolific nature of crAssphage in the human gut microbiome.

Sensitive remote homology search by local alignment of small positional embeddings from protein language models

Accurately detecting distant evolutionary relationships between proteins remains an ongoing challenge in bioinformatics. Search methods based on primary sequence struggle to accurately detect homology between sequences with less than 20% amino acid identity. Profile- and structure-based strategies extend sensitive search capabilities into this twilight zone of sequence similarity but require slow pre-processing steps. Recently, whole-protein and positional embeddings from deep neural networks have shown promise for providing sensitive sequence comparison and annotation at long evolutionary distances. Embeddings are generally faster to compute than profiles and predicted structures but still suffer several drawbacks related to the ability of whole-protein embeddings to discriminate domain-level homology, and the database size and search speed of methods using positional embeddings. In this work, we show that low-dimensionality positional embeddings can be used directly in speed-optimized local search algorithms. As a proof of concept, we use the ESM2 3B model to convert primary sequences directly into the 3D interaction (3Di) alphabet or amino acid profiles and use these embeddings as input to the highly optimized Foldseek, HMMER3, and HH-suite search algorithms. Our results suggest that positional embeddings as small as a single byte can provide sufficient information for dramatically improved sensitivity over amino acid sequence searches without sacrificing search speed.

Comparative genomics of N-acetyl-5-methoxytryptamine members in four Prunus species with insights into bud dormancy and abiotic stress responses in Prunus avium

KEY MESSAGE

This study provides novel insights into the evolution, diversification, and functions of melatonin biosynthesis genes in Prunus species, highlighting their potential role in regulating bud dormancy and abiotic stresses. The biosynthesis of melatonin (MEL) in plants is primarily governed by enzymatic reactions involving key enzymes such as serotonin N-acetyltransferase (SNAT), tryptamine 5-hydroxylase (T5H), N-acetylserotonin methyltransferase (ASMT) and tryptophan decarboxylase (TDC). In this study, we analyzed Melatonin genes in four Prunus species such as Prunus avium (Pavi), Prunus pusilliflora (Ppus), Prunus serulata (Pser), and Prunus persica (Pper) based on comparative genomics approach. Among the four Prunus species, a total of 29 TDCs, 998 T5Hs, 16 SNATs, and 115 ASMTs within the genome of four Prunus genomes. A thorough investigation of melatonin-related genes was carried out using systematic biological methods and comparative genomics. Through phylogenetic analysis, orthologous clusters, Go enrichment, syntenic relationship, and gene duplication analysis, we discovered both similarities and variations in Melatonin genes among these Prunus species. Additionally, our study revealed the existence of unique subgroup members in the Melatonin genes of these species, which were distinct from those found in Arabidopsis genes. Furthermore, the transcriptomic expression analysis revealed the potential significance of melatonin genes in bud dormancy regulation and abiotic stresses. Our extensive results offer valuable perspectives on the evolutionary patterns, intricate expansion, and functions of PavMEL genes. Given their promising attributes, PavTDCs, PavT5H, PavNAT, and three PavASMT genes warrant in-depth exploration as prime candidates for manipulating dormancy in sweet cherry. This was done to lay the foundation for future explorations into the structural and functional aspects of these factors in Prunus species. This study offers significant insights into the functions of ASMT, SNAT, T5H, and TDC genes and sheds light on their roles in Prunus avium. Moreover, it established a robust foundation for further exploration functional characterization of melatonin genes in fruit species.

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae

The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

Proteins containing photosynthetic reaction centre domains modulate FtsZ-based archaeal cell division

Cell division in all domains of life requires the orchestration of many proteins, but in Archaea most of the machinery remains poorly characterized. Here we investigate the FtsZ-based cell division mechanism in Haloferax volcanii and find proteins containing photosynthetic reaction centre (PRC) barrel domains that play an essential role in archaeal cell division. We rename these proteins cell division protein B 1 (CdpB1) and CdpB2. Depletions and deletions in their respective genes cause severe cell division defects, generating drastically enlarged cells. Fluorescence microscopy of tagged FtsZ1, FtsZ2 and SepF in CdpB1 and CdpB2 mutant strains revealed an unusually disordered divisome that is not organized into a distinct ring-like structure. Biochemical analysis shows that SepF forms a tripartite complex with CdpB1/2 and crystal structures suggest that these two proteins might form filaments, possibly aligning SepF and the FtsZ2 ring during cell division. Overall our results indicate that PRC-domain proteins play essential roles in FtsZ-based cell division in Archaea.

Natural variation of immune epitopes reveals intrabacterial antagonism

Plants and animals detect biomolecules termed Microbe-Associated Molecular Patterns (MAMPs) and induce immunity. Agricultural production is severely impacted by pathogens which can be controlled by transferring immune receptors. However, most studies use a single MAMP epitope and the impact of diverse multi-copy MAMPs on immune induction is unknown. Here we characterized the epitope landscape from five proteinaceous MAMPs across 4,228 plant-associated bacterial genomes. Despite the diversity sampled, natural variation was constrained and experimentally testable. Immune perception in both Arabidopsis and tomato depended on both epitope sequence and copy number variation. For example, Elongation Factor Tu is predominantly single copy and 92% of its epitopes are immunogenic. Conversely, 99.9% of bacterial genomes contain multiple Cold Shock Proteins and 46% carry a non-immunogenic form. We uncovered a new mechanism for immune evasion, intrabacterial antagonism, where a non-immunogenic Cold Shock Protein blocks perception of immunogenic forms encoded in the same genome. These data will lay the foundation for immune receptor deployment and engineering based on natural variation.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Virus-prokaryote infection pairs associated with prokaryotic production in a freshwater lake

Viruses infect and kill prokaryotic populations in a density- or frequency-dependent manner and affect carbon cycling. However, the effects of the stratification transition, including the stratified and de-stratified periods, on the changes in prokaryotic and viral communities and their interactions remain unclear. We conducted a monthly survey of the surface and deep layers of a large and deep freshwater lake (Lake Biwa, Japan) for a year and analyzed the prokaryotic production and prokaryotic and viral community composition. Our analysis revealed that, in the surface layer, 19 prokaryotic species, accounting for approximately 40% of the total prokaryotic abundance, could potentially contribute to the majority of prokaryotic production, which is the highest during the summer and is suppressed by viruses. This suggests that a small fraction of prokaryotes and phages were the key infection pairs during the peak period of prokaryotic activity in the freshwater lake. We also found that approximately 50% of the dominant prokaryotic and viral species in the deep layer were present throughout the study period. This suggests that the "kill the winner" model could explain the viral impact on prokaryotes in the surface layer, but other dynamics may be at play in the deep layer. Furthermore, we found that annual vertical mixing could result in a similar rate of community change between the surface and deep layers. These findings may be valuable in understanding how communities and the interaction among them change when freshwater lake stratification is affected by global warming in the future.IMPORTANCEViral infection associated with prokaryotic production occurs in a density- or frequency-dependent manner and regulates the prokaryotic community. Stratification transition and annual vertical mixing in freshwater lakes are known to affect the prokaryotic community and the interaction between prokaryotes and viruses. By pairing measurements of virome analysis and prokaryotic production of a 1-year survey of the depths of surface and deep layers, we revealed (i) the prokaryotic infection pairs associated with prokaryotic production and (ii) the reset in prokaryotic and viral communities through annual vertical mixing in a freshwater lake. Our results provide a basis for future work into changes in stratification that may impact the biogeochemical cycling in freshwater lakes.

Genome-wide identification, expression and function analysis of the MTP gene family in tulip (Tulipa gesneriana)

Currently, soil heavy metal contamination is a severe issue, particularly with Cd pollution. The metal tolerance protein (MTP) proteins, as plant divalent cation transporters, play a crucial role in the transport and tolerance of heavy metals in plants. This study conducted comprehensive identification and characterization of the MTP gene family in the tulip. A total of 11 TgMTP genes were identified and phylogenetically classified into three subfamilies. Conserved motif and gene structure analyses unveiled commonalities and variations among subfamily members. Expression profiling demonstrated several TgMTPs were markedly upregulated under Cd exposure, including the TgMTP7.1. Heterologous expression in yeast validated that TgMTP7.1 could ameliorate Cd sensitivity and enhance its tolerance. These results provide primary insights into the MTP gene family in tulip. Phylogenetic relationships and functional analyses establish a framework for elucidating the transporters and molecular mechanisms governing Cd accumulation and distribution in tulip. Key TgMTPs identified, exemplified by TgMTP7.1, may illuminate molecular breeding efforts aimed at developing Cd-tolerant cultivars for the remediation of soil Cd contamination.

Reappraisal of the DNA phosphorothioate modification machinery: uncovering neglected functional modalities and identification of new counter-invader defense systems

The DndABCDE systems catalysing the unusual phosphorothioate (PT) DNA backbone modification, and the DndFGH systems, which restrict invasive DNA, have enigmatic and paradoxical features. Using comparative genomics and sequence-structure analyses, we show that the DndABCDE module is commonly functionally decoupled from the DndFGH module. However, the modification gene-neighborhoods encode other nucleases, potentially acting as the actual restriction components or suicide effectors limiting propagation of the selfish elements. The modification module's core consists of a coevolving gene-pair encoding the DNA-scanning apparatus - a DndD/CxC-clade ABC ATPase and DndE with two ribbon-helix-helix (MetJ/Arc) DNA-binding domains. Diversification of DndE's DNA-binding interface suggests a multiplicity of target specificities. Additionally, many systems feature DNA cytosine methylase genes instead of PT modification, indicating the DndDE core can recruit other nucleobase modifications. We show that DndFGH is a distinct counter-invader system with several previously uncharacterized domains, including a nucleotide kinase. These likely trigger its restriction endonuclease domain in response to multiple stimuli, like nucleotides, while blocking protective modifications by invader methylases. Remarkably, different DndH variants contain a HerA/FtsK ATPase domain acquired from multiple sources, including cellular genome-segregation systems and mobile elements. Thus, we uncovered novel HerA/FtsK-dependent defense systems that might intercept invasive DNA during replication, conjugation, or packaging.

Soil Humic Acid Stimulates Potentially Active Dissimilatory Arsenate-Reducing Bacteria in Flooded Paddy Soil as Revealed by Metagenomic Stable Isotope Probing

Dissimilatory arsenate reduction contributes a large proportion of arsenic flux from flooded paddy soil, which is closely linked to soil organic carbon input and efflux. Humic acid (HA) represents a natural ingredient in soil and is shown to enhance microbial arsenate respiration to promote arsenic mobility. However, the community and function profiles of metabolically active arsenate-respiring bacteria and their interactions with HA in paddy soil remain unclear. To probe this linkage, we performed a genome-centric comparison of potentially active arsenate-respiring bacteria in anaerobic microcosms amended with 13C-lactate and HA by combining stable-isotope probing with genome-resolved metagenomics. Indeed, HA greatly accelerated the microbial reduction of arsenate to arsenite. Enrichment of bacteria that harbor arsenate-respiring reductase genes (arrA) in HA-enriched 13C-DNA was confirmed by metagenomic binning, which are affiliated with Firmicutes (mainly Desulfitobacterium, Bacillus, Brevibacillus, and Clostridia) and Acidobacteria. Characterization of reference extracellular electron transfer (EET)-related genes in these arrA-harboring bacteria supports the presence of EET-like genes, with partial electron-transport chain genes identified. This suggests that Gram-positive Firmicutes- and Acidobacteria-related members may harbor unspecified EET-associated genes involved in metal reduction. Our findings highlight the link between soil HA and potentially active arsenate-respiring bacteria, which can be considered when using HA for arsenic removal.

Mining of disease-resistance genes in Crocus sativus based on transcriptome sequencing

Introduction: Crocus sativus L. has an important medicinal and economic value in traditional perennial Chinese medicine. However, due to its unique growth characteristics, during cultivation it is highly susceptible to disease. The absence of effective resistance genes restricts us to breed new resistant varieties of C. sativus. Methods: In present study, comprehensive transcriptome sequencing was introduced to explore the disease resistance of the candidate gene in healthy and corm rot-infected C. sativus. Results and discussion: Totally, 43.72 Gb of clean data was obtained from the assembly to generate 65,337 unigenes. By comparing the gene expression levels, 7,575 differentially expressed genes (DEGs) were primarily screened. A majority of the DEGs were completely in charge of defense and metabolism, and 152 of them were annotated as pathogen recognition genes (PRGs) based on the PGRdb dataset. The expression of some transcription factors including NAC, MYB, and WRKY members, changed significantly based on the dataset of transcriptome sequencing. Therefore, this study provides us some valuable information for exploring candidate genes involved in the disease resistance in C. sativus.

Improving deep learning protein monomer and complex structure prediction using DeepMSA2 with huge metagenomics data

Leveraging iterative alignment search through genomic and metagenome sequence databases, we report the DeepMSA2 pipeline for uniform protein single- and multichain multiple-sequence alignment (MSA) construction. Large-scale benchmarks show that DeepMSA2 MSAs can remarkably increase the accuracy of protein tertiary and quaternary structure predictions compared with current state-of-the-art methods. An integrated pipeline with DeepMSA2 participated in the most recent CASP15 experiment and created complex structural models with considerably higher quality than the AlphaFold2-Multimer server (v.2.2.0). Detailed data analyses show that the major advantage of DeepMSA2 lies in its balanced alignment search and effective model selection, and in the power of integrating huge metagenomics databases. These results demonstrate a new avenue to improve deep learning protein structure prediction through advanced MSA construction and provide additional evidence that optimization of input information to deep learning-based structure prediction methods must be considered with as much care as the design of the predictor itself.

Accurate modeling of peptide-MHC structures with AlphaFold

Major histocompatibility complex (MHC) proteins present peptides on the cell surface for T cell surveillance. Reliable in silico prediction of which peptides would be presented and which T cell receptors would recognize them is an important problem in structural immunology. Here, we introduce an AlphaFold-based pipeline for predicting the three-dimensional structures of peptide-MHC complexes for class I and class II MHC molecules. Our method demonstrates high accuracy, outperforming existing tools in class I modeling accuracy and class II peptide register prediction. We validate its performance and utility with new experimental data on a recently described cancer neoantigen/wild-type peptide pair and explore applications toward improving peptide-MHC binding prediction.

Is Novelty Predictable?

Machine learning-based design has gained traction in the sciences, most notably in the design of small molecules, materials, and proteins, with societal applications ranging from drug development and plastic degradation to carbon sequestration. When designing objects to achieve novel property values with machine learning, one faces a fundamental challenge: how to push past the frontier of current knowledge, distilled from the training data into the model, in a manner that rationally controls the risk of failure. If one trusts learned models too much in extrapolation, one is likely to design rubbish. In contrast, if one does not extrapolate, one cannot find novelty. Herein, we ponder how one might strike a useful balance between these two extremes. We focus in particular on designing proteins with novel property values, although much of our discussion is relevant to machine learning-based design more broadly.

Accelerated Molecular Dynamics and AlphaFold Uncover a Missing Conformational State of Transporter Protein OxlT

Transporter proteins change their conformations to carry their substrate across the cell membrane. The conformational dynamics is vital to understanding the transport function. We have studied the oxalate transporter (OxlT), an oxalate:formate antiporter from Oxalobacter formigenes, significant in avoiding kidney stone formation. The atomic structure of OxlT has been recently solved in the outward-open and occluded states. However, the inward-open conformation is still missing, hindering a complete understanding of the transporter. Here, we performed a Gaussian accelerated molecular dynamics simulation to sample the extensive conformational space of OxlT and successfully predicted the inward-open conformation where cytoplasmic substrate formate binding was preferred over oxalate binding. We also identified critical interactions for the inward-open conformation. The results were complemented by an AlphaFold2 structure prediction. Although AlphaFold2 solely predicted OxlT in the outward-open conformation, mutation of the identified critical residues made it partly predict the inward-open conformation, identifying possible state-shifting mutations.

Evolutionary history of the cytochrome P450s from Colletotrichum species and prediction of their putative functional roles during host-pathogen interactions

The genomes of species belonging to the genus Colletotrichum harbor a substantial number of cytochrome P450 monooxygenases (CYPs) encoded by a broad diversity of gene families. However, the biological role of their CYP complement (CYPome) has not been elucidated. Here, we investigated the putative evolutionary scenarios that occurred during the evolution of the CYPome belonging to the Colletotrichum Graminicola species complex (s.c.) and their biological implications. The study revealed that most of the CYPome gene families belonging to the Graminicola s.c. experienced gene contractions. The reductive evolution resulted in species restricted CYPs are predominant in each CYPome of members from the Graminicola s.c., whereas only 18 families are absolutely conserved among these species. However, members of CYP families displayed a notably different phylogenetic relationship at the tertiary structure level, suggesting a putative convergent evolution scenario. Most of the CYP enzymes of the Graminicola s.c. share redundant functions in secondary metabolite biosynthesis and xenobiotic metabolism. Hence, this current work suggests that the presence of a broad CYPome in the genus Colletotrichum plays a critical role in the optimization of the colonization capability and virulence.

Global biogeography and ecological implications of cobamide-producing prokaryotes

Cobamides, a class of essential coenzymes synthesized only by a subset of prokaryotes, are model nutrients in microbial interaction studies and play significant roles in global ecosystems. Yet, their spatial patterns and functional roles remain poorly understood. Herein, we present an in-depth examination of cobamide-producing microorganisms, drawn from a comprehensive analysis of 2862 marine and 2979 soil metagenomic samples. A total of 1934 nonredundant metagenome-assembled genomes (MAGs) potentially capable of producing cobamides de novo were identified. The cobamide-producing MAGs are taxonomically diverse but habitat specific. They constituted only a fraction of all the recovered MAGs, with the majority of MAGs being potential cobamide users. By mapping the distribution of cobamide producers in marine and soil environments, distinct latitudinal gradients were observed: the marine environment showed peak abundance at the equator, whereas soil environments peaked at mid-latitudes. Importantly, significant and positive links between the abundance of cobamide producers and the diversity and functions of microbial communities were observed, as well as their promotional roles in essential biogeochemical cycles. These associations were more pronounced in marine samples than in soil samples, which suggests a heightened propensity for microorganisms to engage in cobamide sharing in fluid environments relative to the more spatially restricted soil environment. These findings shed light on the global patterns and potential ecological roles of cobamide-producing microorganisms in marine and soil ecosystems, enhancing our understanding of large-scale microbial interactions.

Exploring the versatility of sesquiterpene biosynthesis in guava plants: a comparative genome-wide analysis of two cultivars

Psidium guajava L., a fruit crop belonging to the Myrtaceae family, is highly valued for its nutritional and medicinal properties. The family exhibits a diverse chemical profile of essential oils and serves as a valuable resource due to its ecological interactions, adaptability, and dispersal capacity. The Myrtaceae family has been extensively studied for its terpenoids. Genetic studies have focused on foliar terpene yield in species from the Eucalypteae and Melaleucaceae tribes. To understand the evolutionary trends in guava breeding, this study predicted terpene synthase genes (TPS) from different cultivars. Through this analysis, 43 full-length TPS genes were identified, and approximately 77% of them exhibited relative expression in at least one of the five investigated plant tissues (root, leaf, bud, flower, and fruit) of two guava cultivars. We identified intra-species variation in the terpene profile and single nucleotide polymorphisms (SNPs) in twelve TPS genes, resulting in the clustering of 62 genotypes according to their essential oil chemotypes. The high concentration of sesquiterpenes is supported by the higher number of TPS-a genes and their expression. The expansion for TPS sub-families in P. guajava occurred after the expansion of other rosids species. Providing insight into the origin of structural diversification and expansion in each clade of the TPS gene family within Myrtaceae. This study can provide insights into the diversity of genes for specialized metabolites such as terpenes, and their regulation, which can lead to a diverse chemotype of essential oil in different tissues and genotypes. This suggests a mode of enzymatic evolution that could lead to high sesquiterpene production, act as a chemical defense and contribute to the adaptive capacity of this species to different habitats.

Boron homeostasis affects Longan yield: a study of NIP and BOR boron transporter of two cultivars

BACKGROUND

Essential micronutrient Boron (B) plays crucial roles in plant survival and reproduction but becomes toxic in higher quantities. Although plant cells have different B transport systems, B homeostasis is mainly maintained by two transporter protein families: B exporters (BOR) and nodulin-26-like intrinsic proteins (NIP). Their diversity and differential expression are responsible for varied B tolerance among plant varieties and species. Longan is a highly admired subtropical fruit with a rising market in China and beyond. In the present study, we cultured Shixia (SX) and Yiduo (YD), two differently characterized Longan cultivars, with foliar B spray. We analyzed their leaf physiology, fruit setting, B content, and boron transporter gene expression of various tissue samples. We also traced some of these genes' subcellular localization and overexpression effects.

RESULTS

YD and SX foliage share similar microstructures, except the mesophyll cell wall thickness is double in YD. The B spray differently influenced their cellular constituents and growth regulators. Gene expression analysis showed reduced BOR genes expression and NIP genes differential spatiotemporal expression. Using green fluorescent protein, two high-expressing NIPs, NIP1 and NIP19, were found to translocate in the transformed tobacco leaves' cell membrane. NIPs transformation of SX pollen was confirmed using magnetic beads and quantified using a fluorescence microscope and polymerase chain reaction. An increased seed-setting rate was observed when YD was pollinated using these pollens. Between the DlNIP1 and DlNIP19 transformed SX pollen, the former germinated better with increasing B concentrations and, compared to naturally pollinated plants, had a better seed-setting rate in YD♀ × SX♂.

CONCLUSION

SX and YD Longan have different cell wall structures and react differently to foliar B spray, indicating distinct B tolerance and management. Two B transporter NIP genes were traced to localize in the plasma membrane. However, under high B concentrations, their differential expression resulted in differences in Jasmonic acid content, leading to differences in germination rate. Pollination of YD using these NIPs transformed SX pollen also showed NIP1 overexpression might overcome the unilateral cross incompatibility between YD♀ × SX♂ and can be used to increase Longan production.

Quantifying microbial guilds

The ecological role of microorganisms is of utmost importance due to their multiple interactions with the environment. However, assessing the contribution of individual taxonomic groups has proven difficult despite the availability of high throughput data, hindering our understanding of such complex systems. Here, we propose a quantitative definition of guild that is readily applicable to metagenomic data. Our framework focuses on the functional character of protein sequences, as well as their diversifying nature. First, we discriminate functional sequences from the whole sequence space corresponding to a gene annotation to then quantify their contribution to the guild composition across environments. In addition, we identify and distinguish functional implementations, which are sequence spaces that have different ways of carrying out the function. In contrast, we found that orthology delineation did not consistently align with ecologically (or functionally) distinct implementations of the function. We demonstrate the value of our approach with two case studies: the ammonia oxidation and polyamine uptake guilds from the Malaspina circumnavigation cruise, revealing novel ecological dynamics of the latter in marine ecosystems. Thus, the quantification of guilds helps us to assess the functional role of different taxonomic groups with profound implications on the study of microbial communities.

Structure and dispersion of the conjugative mobilome in surface ocean bacterioplankton

Mobile genetic elements (MGEs), collectively referred to as the "mobilome", can have a significant impact on the fitness of microbial communities and therefore on ecological processes. Marine MGEs have mainly been associated with wide geographical and phylogenetic dispersal of adaptative traits. However, whether the structure of this mobilome exhibits deterministic patterns in the natural community is still an open question. The aim of this study was to characterize the structure of the conjugative mobilome in the ocean surface bacterioplankton by searching the publicly available marine metagenomes from the TARA Oceans survey, together with molecular markers, such as relaxases and type IV coupling proteins of the type IV secretion system (T4SS). The T4SS machinery was retrieved in more abundance than relaxases in the surface marine bacterioplankton. Moreover, among the identified MGEs, mobilizable elements were the most abundant, outnumbering self-conjugative sequences. Detection of a high number of incomplete T4SSs provides insight into possible strategies related to trans-acting activity between MGEs, and accessory functions of the T4SS (e.g. protein secretion), allowing the host to maintain a lower metabolic burden in the highly dynamic marine system. Additionally, the results demonstrate a wide geographical dispersion of MGEs throughout oceanic regions, while the Southern Ocean appears segregated from other regions. The marine mobilome also showed a high similarity of functions present in known plasmid databases. Moreover, cargo genes were mostly related to DNA processing, but scarcely associated with antibiotic resistance. Finally, within the MGEs, integrative and conjugative elements showed wider marine geographic dispersion than plasmids.

Limitation of amino acid availability by bacterial populations during enhanced colitis in IBD mouse model

IMPORTANCE

Inflammatory bowel disease is associated with an increase in Enterobacteriaceae and Enterococcus species; however, the specific mechanisms are unclear. Previous research has reported the associations between microbiota and inflammation, here we investigate potential pathways that specific bacteria populations use to drive gut inflammation. Richie et al. show that these bacterial populations utilize an alternate sulfur metabolism and are tolerant of host-derived immune-response products. These metabolic pathways drive host gut inflammation and fuel over colonization of these pathobionts in the dysbiotic colon. Cultured isolates from dysbiotic mice indicated faster growth supplemented with L-cysteine, showing these microbes can utilize essential host nutrients.

With a pinch of salt: metagenomic insights into Namib Desert salt pan microbial mats and halites reveal functionally adapted and competitive communities

IMPORTANCE

The hyperarid Namib Desert is one of the oldest deserts on Earth. It contains multiple clusters of playas which are saline-rich springs surrounded by halite evaporites. Playas are of great ecological importance, and their indigenous (poly)extremophilic microorganisms are potentially involved in the precipitation of minerals such as carbonates and sulfates and have been of great biotechnological importance. While there has been a considerable amount of microbial ecology research performed on various Namib Desert edaphic microbiomes, little is known about the microbial communities inhabiting its multiple playas. In this work, we provide a comprehensive taxonomic and functional potential characterization of the microbial, including viral, communities of sediment mats and halites from two distant salt pans of the Namib Desert, contributing toward a better understanding of the ecology of this biome.

Endophytic bacteria associated with wild-type banana seed (Musa balbisiana): whole genome sequencing

We present the whole-genome sequences of five endophytic bacteria isolated from Musa balbisiana seeds. These strains represent five different genera: Bacillus, Brachybacterium, Enterobacter, Enterococcus, and Pantoea. Among these, three genera (Bacillus, Pantoea, and Enterobacter) were previously recognized for their antagonistic effects against Fusarium wilt, a highly destructive disease that affects banana plants.

Whole genome structural predictions reveal hidden diversity in putative oxidative enzymes of the lignocellulose-degrading ascomycete Parascedosporium putredinis NO1

IMPORTANCE

An annotated reference genome has revealed P. putredinis NO1 as a useful resource for the identification of new lignocellulose-degrading enzymes for biorefining of woody plant biomass. Utilizing a "structure-omics"-based searching strategy, we identified new potentially lignocellulose-active sequences that would have been missed by traditional sequence searching methods. These new identifications, alongside the discovery of novel enzymatic functions from this underexplored lineage with the recent discovery of a new phenol oxidase that cleaves the main structural β-O-4 linkage in lignin from P. putredinis NO1, highlight the underexplored and poorly represented family Microascaceae as a particularly interesting candidate worthy of further exploration toward the valorization of high value biorenewable products.

Broad-range capsule-dependent lytic Sugarlandvirus against Klebsiella sp

IMPORTANCE

The emergence of multi-drug resistant bacteria is a global health problem. Among them, Klebsiella pneumoniae is considered a high-priority pathogen, making it necessary to develop new therapeutic tools to reduce the bacterial burden in an effective and sustainable manner. Phages, bacterial viruses, are very promising tools. However, phages are highy specific, rendering large-scale therapeutics costly to implement. This is especially certain in Klebsiella, a capsular bacterium in which phages have been shown to be capsular type dependent, infecting one or a few capsular types through specific enzymes called depolymerases. In this study, we have isolated and characterized novel phages with lytic ability against bacteria from a wide variety of capsular types, representing the Klebsiella phages with the widest range of infection described. Remarkably, these broad-range phages showed capsule dependency, despite the absence of depolymerases in their genomes, implying that infectivity could be governed by alternative mechanisms yet to be uncovered.

Gene Family Expansion during the Adaptation of Colletotrichum gloeosporioides to Woody Plants

Gene gains/losses during evolution are critical for the adaptation of organisms to new environments or hosts. However, it remains unknown whether gene family expansions facilitated the adaptation of phytopathogenic fungi to woody plants. In this study, we compared the newly sequenced genome of the Colletotrichum gloeosporioides strain CFCC80308 with the genomes of two other C. gloeosporioides strains, Cg-14 and Lc-1, isolated from Persea americana and Liriodendron leaves, respectively. The genes in the expanded families, which were associated with plant surface signal recognition, encoded various proteins, including glycosyde hydrolases (GHs) and cytochrome P450. Interestingly, there was a substantial increase in the number of GH family genes in CFCC80308. Specifically, there were 368 enriched genes in the GH families (e.g., GH1, GH3, GH10, GH12, GH15, GH16, GH17, GH18, GH25, GH32, GH53, GH61, GH76, and GH81); the expression levels of these genes were highly up-regulated during the infection of poplar trees. Additionally, the GH17 family was larger in CFCC80308 than in C. gloeosporioides strains Cg-14 and Lc-1. Furthermore, the expansion of the MP65-encoding gene family during the adaptation of Colletotrichum species to woody plants was consistent with the importance of gene gains/losses for the adaptation of organisms to their environments. This study has clarified how C. gloeosporioides adapted to woody plants during evolution.

SPLASH: A statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), which directly analyzes raw sequencing data, using a statistical test to detect a signature of regulation: sample-specific sequence variation. SPLASH detects many types of variation and can be efficiently run at scale. We show that SPLASH identifies complex mutation patterns in SARS-CoV-2, discovers regulated RNA isoforms at the single-cell level, detects the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a unifying approach to genomic analysis that enables expansive discovery without metadata or references.

Diversification of division mechanisms in endospore-forming bacteria revealed by analyses of peptidoglycan synthesis in Clostridioides difficile

The bacterial enzymes FtsW and FtsI, encoded in the highly conserved dcw gene cluster, are considered to be universally essential for the synthesis of septal peptidoglycan (PG) during cell division. Here, we show that the pathogen Clostridioides difficile lacks a canonical FtsW/FtsI pair, and its dcw-encoded PG synthases have undergone a specialization to fulfill sporulation-specific roles, including synthesizing septal PG during the sporulation-specific mode of cell division. Although these enzymes are directly regulated by canonical divisome components during this process, dcw-encoded PG synthases and their divisome regulators are dispensable for cell division during normal growth. Instead, C. difficile uses a bifunctional class A penicillin-binding protein as the core divisome PG synthase, revealing a previously unreported role for this class of enzymes. Our findings support that the emergence of endosporulation in the Firmicutes phylum facilitated the functional repurposing of cell division factors. Moreover, they indicate that C. difficile, and likely other clostridia, assemble a distinct divisome that therefore may represent a unique target for therapeutic interventions.

Brevibacillus ruminantium sp. nov., isolated from cow faeces

An aerobic, Gram-stain-positive, rod-shaped, endospore-forming bacterial strain, designated BB3-R1T, was isolated from cow faeces sampled in Daejeon, Republic of Korea. Growth was observed at 25-45 °C (optimum, 35-40 °C) and pH 7.0-9.0 (optimum, pH 8.0), with up to 3 % (w/v) NaCl (optimum, 0 % NaCl). blast analysis of 16S rRNA gene sequences revealed the highest sequence similarity of strain BB3-R1T to Brevibacillus borstelensis NRRL NRS-818T (98.8 %) followed by Brevibacillus panacihumi JCM 15085T (97.5 %). According to 16S rRNA gene and whole-genome based phylogenetic trees, strain BB3-R1T clustered with Brevibacillus composti FJAT-54423T and B. borstelensis NRRL NRS-818T. OrthoANI and dDDH values of strain BB3-R1T with the closely related strains were lower than 77.5 and 26.8 %, respectively. The major menaquinones and polar lipids of the strain were MK-7 and phosphatidylmonomethylethanolamine, diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine, respectively. The major fatty acids (>10 %) were C14 : 0 iso, C15 : 0 iso, C15 : 0 anteiso and C16 : 1  ω7c alcohol. The cell-wall peptidoglycan contained cross-linked meso-diaminopimelic acid (type A1 gamma). The phenotypic, chemotaxonomic and genotypic data obtained in this study showed that the strain represents a novel species of the genus Brevibacillus, for which the name Brevibacillus ruminantium sp. nov. (type strain BB3-R1T=KACC 22663T=NBRC 115962T) is proposed.