Quantifying functional redundancy in polysaccharide-degrading prokaryotic communities.

BACKGROUND: Functional redundancy (FR) is widely present, but there is no consensus on its formation process and influencing factors. Taxonomically distinct microorganisms possessing genes for the same function in a community lead to within-community FR, and distinct assemblies of microorganisms in different communities playing the same functional roles are termed between-community FR. We proposed two formulas to respectively quantify the degree of functional redundancy within and between communities and analyzed the FR degrees of carbohydrate degradation functions in global environment samples using the genetic information of glycoside hydrolases (GHs) encoded by prokaryotes. RESULTS: Our results revealed that GHs are each encoded by multiple taxonomically distinct prokaryotes within a community, and the enzyme-encoding prokaryotes are further distinct between almost any community pairs. The within- and between-FR degrees are primarily affected by the alpha and beta community diversities, respectively, and are also affected by environmental factors (e.g., pH, temperature, and salinity). The FR degree of the prokaryotic community is determined by deterministic factors. CONCLUSIONS: We conclude that the functional redundancy of GHs is a stabilized community characteristic. This study helps to determine the FR formation process and influencing factors and provides new insights into the relationships between prokaryotic community biodiversity and ecosystem functions. Video Abstract.

Nitrate reduction to ammonium: a phylogenetic, physiological, and genetic aspects in Prokaryotes and eukaryotes.

The microbe-mediated conversion of nitrate (NO3-) to ammonium (NH4+) in the nitrogen cycle has strong implications for soil health and crop productivity. The role of prokaryotes, eukaryotes and their phylogeny, physiology, and genetic regulations are essential for understanding the ecological significance of this empirical process. Several prokaryotes (bacteria and archaea), and a few eukaryotes (fungi and algae) are reported as NO3- reducers under certain conditions. This process involves enzymatic reactions which has been catalysed by nitrate reductases, nitrite reductases, and NH4+-assimilating enzymes. Earlier reports emphasised that single-cell prokaryotic or eukaryotic organisms are responsible for this process, which portrayed a prominent gap. Therefore, this study revisits the similarities and uniqueness of mechanism behind NO3- -reduction to NH4+ in both prokaryotes and eukaryotes. Moreover, phylogenetic, physiological, and genetic regulation also shed light on the evolutionary connections between two systems which could help us to better explain the NO3--reduction mechanisms over time. Reports also revealed that certain transcription factors like NtrC/NtrB and Nit2 have shown a major role in coordinating the expression of NO3- assimilation genes in response to NO3- availability. Overall, this review provides a comprehensive information about the complex fermentative and respiratory dissimilatory nitrate reduction to ammonium (DNRA) processes. Uncovering the complexity of this process across various organisms may further give insight into sustainable nitrogen management practices and might contribute to addressing global environmental challenges.

Internal carbon recycling by heterotrophic prokaryotes compensates for mismatches between phytoplankton production and heterotrophic consumption.

Molecular observational tools are useful for characterizing the composition and genetic endowment of microbial communities but cannot measure fluxes, which are critical for the understanding of ecosystems. To overcome these limitations, we used a mechanistic inference approach to estimate dissolved organic carbon (DOC) production and consumption by phytoplankton operational taxonomic units and heterotrophic prokaryotic amplicon sequence variants and inferred carbon fluxes between members of this microbial community from Western English Channel time-series data. Our analyses focused on phytoplankton spring and summer blooms, as well as bacteria summer blooms. In spring blooms, phytoplankton DOC production exceeds heterotrophic prokaryotic consumption, but in bacterial summer blooms heterotrophic prokaryotes consume three times more DOC than produced by the phytoplankton. This mismatch is compensated by heterotrophic prokaryotic DOC release by death, presumably from viral lysis. In both types of summer blooms, large amounts of the DOC liberated by heterotrophic prokaryotes are reused through internal recycling, with fluxes between different heterotrophic prokaryotes being at the same level as those between phytoplankton and heterotrophic prokaryotes. In context, internal recycling accounts for approximately 75% and 30% of the estimated net primary production (0.16 vs 0.22 and 0.08 vs 0.29 µmol l-1 d-1) in bacteria and phytoplankton summer blooms, respectively, and thus represents a major component of the Western English Channel carbon cycle. We have concluded that internal recycling compensates for mismatches between phytoplankton DOC production and heterotrophic prokaryotic consumption, and we encourage future analyses on aquatic carbon cycles to investigate fluxes between heterotrophic prokaryotes, specifically internal recycling.

An evolutionary optimum amid moderate heritability in prokaryotic cell size.

We investigate the distribution and evolution of prokaryotic cell size based on a compilation of 5,380 species. Size spans four orders of magnitude, from 100 nm (Mycoplasma) to more than 1 cm (Thiomargarita); however, most species congregate heavily around the mean. The distribution approximates but is distinct from log normality. Comparative phylogenetics suggests that size is heritable, yet the phylogenetic signal is moderate, and the degree of heritability is independent of taxonomic scale (i.e., fractal). Evolutionary modeling indicates the presence of an optimal cell size to which most species gravitate. The size is equivalent to a coccus of 0.70 µm in diameter. Analyses of 1,361 species with sequenced genomes show that genomic traits contribute to size evolution moderately and synergistically. Given our results, scaling theory, and empirical evidence, we discuss potential drivers that may expand or shrink cells around the optimum and propose a stability landscape model for prokaryotic cell size.

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.

Annotation and Comparative Genomics of Prokaryotic Transposable Elements.

The data generated in nearly 30 years of bacterial genome sequencing has revealed the abundance of transposable elements (TE) and their importance in genome and transcript remodeling through the mediation of DNA insertions and deletions, structural rearrangements, and regulation of gene expression. Furthermore, what we have learned from studying transposition mechanisms and their regulation in bacterial TE is fundamental to our current understanding of TE in other organisms because much of what has been observed in bacteria is conserved in all domains of life. However, unlike eukaryotic TE, prokaryotic TE sequester and transmit important classes of genes that impact host fitness, such as resistance to antibiotics and heavy metals and virulence factors affecting animals and plants, among other acquired traits. This provides dynamism and plasticity to bacteria, which would otherwise be propagated clonally. The insertion sequences (IS), the simplest form of prokaryotic TE, are autonomous and compact mobile genetic elements. These can be organized into compound transposons, in which two similar IS can flank any DNA segment and render it transposable. Other more complex structures, called unit transposons, can be grouped into four major families (Tn3, Tn7, Tn402, Tn554) with specific genetic characteristics. This chapter will revisit the prominent structural features of these elements, focusing on a genomic annotation framework and comparative analysis. Relevant aspects of TE will also be presented, stressing their key position in genome impact and evolution, especially in the emergence of antimicrobial resistance and other adaptive traits.

Investigating the nature of prokaryotic genomic island locations within a genome.

Horizontal gene transfer (HGT) is a powerful evolutionary force that considerably shapes the structure of prokaryotic genomes and is associated with genomic islands (GIs). A GI is a DNA segment composed of transferred genes that can be found within a prokaryotic genome, obtained through HGT. Much research has focused on detecting GIs in genomes, but here we pursue a new course, which is identifying possible preferred locations of GIs in the prokaryotic genome. Here, we identify the locations of the GIs within prokaryotic genomes to examine patterns in those locations. Prokaryotic GIs were analyzed according to the genome structure that they are located in, whether it be a circular or a linear genome. The analytical investigations employed are: (1) studying the GI locations in relation to the origin of replication (oriC); (2) exploring the distances between GIs; and (3) determining the distribution of GIs across the genomes. For each of the investigations, the analysis was performed on all of the GIs in the data set. Moreover, to void bias caused by the distribution of the genomes represented, the GIs in one genome from each species and the GIs of the most frequent species are also analyzed. Overall, the results showed that there are preferred sites for the GIs in the genome. In the linear genomes, these sites are usually located in the oriC region and terminus region, while in the circular genomes, they are located solely in the terminus region. These results also showed that the distance distribution between the GIs is almost exponential, which proves that GIs have preferred sites within genomes. The oriC and termniuns are preferred sites for the GIs and a possible natural explanation for this could be connected to the content of the oriC region. Moreover, the content of the GIs in terms of its protein families was studied and the results demonstrated that the majority of frequent protein families are close to identical in each section.

Macroalgal virosphere assists with host-microbiome equilibrium regulation and affects prokaryotes in surrounding marine environments.

The microbiomes in macroalgal holobionts play vital roles in regulating macroalgal growth and ocean carbon cycling. However, the virospheres in macroalgal holobionts remain largely underexplored, representing a critical knowledge gap. Here we unveil that the holobiont of kelp (Saccharina japonica) harbors highly specific and unique epiphytic/endophytic viral species, with novelty (99.7% unknown) surpassing even extreme marine habitats (e.g. deep-sea and hadal zones), indicating that macroalgal virospheres, despite being closest to us, are among the least understood. These viruses potentially maintain microbiome equilibrium critical for kelp health via lytic-lysogenic infections and the expression of folate biosynthesis genes. In-situ kelp mesocosm cultivation and metagenomic mining revealed that kelp holobiont profoundly reshaped surrounding seawater and sediment virus-prokaryote pairings through changing surrounding environmental conditions and virus-host migrations. Some kelp epiphytic viruses could even infect sediment autochthonous bacteria after deposition. Moreover, the presence of ample viral auxiliary metabolic genes for kelp polysaccharide (e.g. laminarin) degradation underscores the underappreciated viral metabolic influence on macroalgal carbon cycling. This study provides key insights into understanding the previously overlooked ecological significance of viruses within macroalgal holobionts and the macroalgae-prokaryotes-virus tripartite relationship.

Substrate uptake patterns shape niche separation in marine prokaryotic microbiome.

Marine heterotrophic prokaryotes primarily take up ambient substrates using transporters. The patterns of transporters targeting particular substrates shape the ecological role of heterotrophic prokaryotes in marine organic matter cycles. Here, we report a size-fractionated pattern in the expression of prokaryotic transporters throughout the oceanic water column due to taxonomic variations, revealed by a multi-"omics" approach targeting ATP-binding cassette (ABC) transporters and TonB-dependent transporters (TBDTs). Substrate specificity analyses showed that marine SAR11, Rhodobacterales, and Oceanospirillales use ABC transporters to take up organic nitrogenous compounds in the free-living fraction, while Alteromonadales, Bacteroidetes, and Sphingomonadales use TBDTs for carbon-rich organic matter and metal chelates on particles. The expression of transporter proteins also supports distinct lifestyles of deep-sea prokaryotes. Our results suggest that transporter divergency in organic matter assimilation reflects a pronounced niche separation in the prokaryote-mediated organic matter cycles.

Can Sb(III)-oxidizing prokaryote also oxidize As(III) under aerobic and anaerobic conditions, and vice versa?

Sb(III) and As(III) share similar chemical features and coexist in the environment. However, their oxidase enzymes have completely different sequences and structures. This raises an intriguing question: Could Sb(III)-oxidizing prokaryotes (SOPs) also oxidize As(III), and vice versa? Regarding this issue, previous investigations have yielded unclear, incorrect and even conflicting data. This work aims to address this matter. First, we prepared an enriched population of SOPs that comprises 55 different AnoA genes, lacking AioAB and ArxAB genes. We found that these SOPs can oxidize both Sb(III) and As(III) with comparable capabilities. To further confirm this finding, we isolated three cultivable SOP strains that have AnoA gene, but lack AioAB and ArxAB genes. We observed that they also oxidize both Sb(III) and As(III) under both anaerobic and aerobic conditions. Secondly, we obtained an enriched population of As(III)-oxidizing prokaryotes (AOPs) from As-contaminated soils, which comprises 69 different AioA genes, lacking AnoA gene. We observed that the AOP population has significant As(III)-oxidizing activities, but lack detectable Sb(III)-oxidizing activities under both aerobic and anaerobic conditions. Therefore, we convincingly show that SOPs can oxidize As(III), but AOPs cannot oxidize Sb(III). These findings clarify the previous ambiguities, confusion, errors or contradictions regarding how SOPs and AOPs oxidize each other's substrate.

DeepReg: a deep learning hybrid model for predicting transcription factors in eukaryotic and prokaryotic genomes.

Deep learning models (DLMs) have gained importance in predicting, detecting, translating, and classifying a diversity of inputs. In bioinformatics, DLMs have been used to predict protein structures, transcription factor-binding sites, and promoters. In this work, we propose a hybrid model to identify transcription factors (TFs) among prokaryotic and eukaryotic protein sequences, named Deep Regulation (DeepReg) model. Two architectures were used in the DL model: a convolutional neural network (CNN), and a bidirectional long-short-term memory (BiLSTM). DeepReg reached a precision of 0.99, a recall of 0.97, and an F1-score of 0.98. The quality of our predictions, the bias-variance trade-off approach, and the characterization of new TF predictions were evaluated and compared against those produced by DeepTFactor, as well as against experimental data from three model organisms. Predictions based on our DLM tended to exhibit less variance and bias than those from DeepTFactor, thus increasing reliability and decreasing overfitting.

Co-transcriptional gene regulation in eukaryotes and prokaryotes.

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

A widespread bacterial protein compartment sequesters and stores elemental sulfur.

Subcellular compartments often serve to store nutrients or sequester labile or toxic compounds. As bacteria mostly do not possess membrane-bound organelles, they often have to rely on protein-based compartments. Encapsulins are one of the most prevalent protein-based compartmentalization strategies found in prokaryotes. Here, we show that desulfurase encapsulins can sequester and store large amounts of crystalline elemental sulfur. We determine the 1.78-angstrom cryo-EM structure of a 24-nanometer desulfurase-loaded encapsulin. Elemental sulfur crystals can be formed inside the encapsulin shell in a desulfurase-dependent manner with l-cysteine as the sulfur donor. Sulfur accumulation can be influenced by the concentration and type of sulfur source in growth medium. The selectively permeable protein shell allows the storage of redox-labile elemental sulfur by excluding cellular reducing agents, while encapsulation substantially improves desulfurase activity and stability. These findings represent an example of a protein compartment able to accumulate and store elemental sulfur.

Intracellular signaling in proto-eukaryotes evolves to alleviate regulatory conflicts of endosymbiosis.

The complex eukaryotic cell resulted from a merger between simpler prokaryotic cells, yet the role of the mitochondrial endosymbiosis with respect to other eukaryotic innovations has remained under dispute. To investigate how the regulatory challenges associated with the endosymbiotic state impacted genome and network evolution during eukaryogenesis, we study a constructive computational model where two simple cells are forced into an obligate endosymbiosis. Across multiple in silico evolutionary replicates, we observe the emergence of different mechanisms for the coordination of host and symbiont cell cycles, stabilizing the endosymbiotic relationship. In most cases, coordination is implicit, without signaling between host and symbiont. Signaling only evolves when there is leakage of regulatory products between host and symbiont. In the fittest evolutionary replicate, the host has taken full control of the symbiont cell cycle through signaling, mimicking the regulatory dominance of the nucleus over the mitochondrion that evolved during eukaryogenesis.

Unraveling the plasticity of translation initiation in prokaryotes: Beyond the invariant Shine-Dalgarno sequence.

Translation initiation in prokaryotes is mainly defined, although not exclusively, by the interaction between the anti-Shine-Dalgarno sequence (antiSD), located at the 3'-terminus of the 16S ribosomal RNA, and a complementary sequence, the ribosome binding site, or Shine-Dalgarno (SD), located upstream of the start codon in prokaryotic mRNAs. The antiSD has a conserved 5'-CCUCC-3' core, but inter-species variations have been found regarding the participation of flanking bases in binding. These variations have been described for certain bacteria and, to a lesser extent, for some archaea. To further analyze these variations, we conducted binding-energy prediction analyses on over 6,400 genomic sequences from both domains. We identified 15 groups of antiSD variants that could be associated with the organisms' phylogenetic origin. Additionally, our findings revealed that certain organisms exhibit variations in the core itself. Importantly, an unaltered core is not necessarily required for the interaction between the 3'-terminus of the rRNA and the region preceding the AUG of the mRNA. In our study, we classified organisms into four distinct categories: i) those possessing a conserved core and demonstrating binding; ii) those with a conserved core but lacking evidence of binding; iii) those exhibiting binding in the absence of a conserved core; and iv) those lacking both a conserved core and evidence of binding. Our results demonstrate the flexibility of organisms in evolving different sequences involved in translation initiation beyond the traditional Shine-Dalgarno sequence. These findings are discussed in terms of the evolution of translation initiation in prokaryotic organisms.

The missing part: the Archaeoglobus fulgidus Argonaute forms a functional heterodimer with an N-L1-L2 domain protein.

Argonaute (Ago) proteins are present in all three domains of life (bacteria, archaea and eukaryotes). They use small (15-30 nucleotides) oligonucleotide guides to bind complementary nucleic acid targets and are responsible for gene expression regulation, mobile genome element silencing, and defence against viruses or plasmids. According to their domain organization, Agos are divided into long and short Agos. Long Agos found in prokaryotes (long-A and long-B pAgos) and eukaryotes (eAgos) comprise four major functional domains (N, PAZ, MID and PIWI) and two structural linker domains L1 and L2. The majority (∼60%) of pAgos are short pAgos, containing only the MID and inactive PIWI domains. Here we focus on the prokaryotic Argonaute AfAgo from Archaeoglobus fulgidus DSM4304. Although phylogenetically classified as a long-B pAgo, AfAgo contains only MID and catalytically inactive PIWI domains, akin to short pAgos. We show that AfAgo forms a heterodimeric complex with a protein encoded upstream in the same operon, which is a structural equivalent of the N-L1-L2 domains of long pAgos. This complex, structurally equivalent to a long PAZ-less pAgo, outperforms standalone AfAgo in guide RNA-mediated target DNA binding. Our findings provide a missing piece to one of the first and the most studied pAgos.

Membrane fusion and fission during eukaryogenesis.

All eukaryotes can be traced back to a single shared ancestral lineage that emerged from interactions between different prokaryotic cells. Current models of eukaryogenesis describe various selective forces and evolutionary mechanisms that contributed to the formation of eukaryotic cells. Central to this process were significant changes in cellular structure, resulting in the configuration of a new cell type characterized by internal membrane compartments. Additionally, eukaryogenesis results in a life cycle that relies on cell-cell fusion. We discuss the potential roles of proteins involved in remodeling cellular membranes, highlighting two critical stages in the evolution of eukaryotes: the internalization of symbiotic partners and a scenario wherein the emergence of sexual reproduction is linked to a polyploid ancestor generated by cell-cell fusion.

The nuclease-associated short prokaryotic Argonaute system nonspecifically degrades DNA upon activation by target recognition.

Prokaryotic Argonautes (pAgos) play a vital role in host defense by utilizing short nucleic acid guides to recognize and target complementary nucleic acids. Despite being the majority of pAgos, short pAgos have only recently received attention. Short pAgos are often associated with proteins containing an APAZ domain and a nuclease domain including DUF4365, SMEK, or HNH domain. In contrast to long pAgos that specifically cleave the target DNA, our study demonstrates that the short pAgo from Thermocrispum municipal, along with its associated DUF4365-APAZ protein, forms a heterodimeric complex. Upon RNA-guided target DNA recognition, this complex is activated to nonspecifically cleave DNA. Additionally, we found that the TmuRE-Ago complex shows a preference for 5'-OH guide RNA, specifically requires a uridine nucleotide at the 5' end of the guide RNA, and is sensitive to single-nucleotide mismatches between the guide RNA and target DNA. Based on its catalytic properties, our study has established a novel nucleic acid detection method and demonstrated its feasibility. This study not only expands our understanding of the defense mechanism employed by short pAgo systems but also suggests their potential applications in nucleic acid detection.

Auto-inhibition and activation of a short Argonaute-associated TIR-APAZ defense system.

Short prokaryotic Ago accounts for most prokaryotic Argonaute proteins (pAgos) and is involved in defending bacteria against invading nucleic acids. Short pAgo associated with TIR-APAZ (SPARTA) has been shown to oligomerize and deplete NAD+ upon guide-mediated target DNA recognition. However, the molecular basis of SPARTA inhibition and activation remains unknown. In this study, we determined the cryogenic electron microscopy structures of Crenotalea thermophila SPARTA in its inhibited, transient and activated states. The SPARTA monomer is auto-inhibited by its acidic tail, which occupies the guide-target binding channel. Guide-mediated target binding expels this acidic tail and triggers substantial conformational changes to expose the Ago-Ago dimerization interface. As a result, SPARTA assembles into an active tetramer, where the four TIR domains are rearranged and packed to form NADase active sites. Together with biochemical evidence, our results provide a panoramic vision explaining SPARTA auto-inhibition and activation and expand understanding of pAgo-mediated bacterial defense systems.

Catalytically inactive long prokaryotic Argonaute systems employ distinct effectors to confer immunity via abortive infection.

Argonaute proteins (Agos) bind short nucleic acids as guides and are directed by them to recognize target complementary nucleic acids. Diverse prokaryotic Agos (pAgos) play potential functions in microbial defense. The functions and mechanisms of a group of full-length yet catalytically inactive pAgos, long-B pAgos, remain unclear. Here, we show that most long-B pAgos are functionally connected with distinct associated proteins, including nucleases, Sir2-domain-containing proteins and trans-membrane proteins, respectively. The long-B pAgo-nuclease system (BPAN) is activated by guide RNA-directed target DNA recognition and performs collateral DNA degradation in vitro. In vivo, the system mediates genomic DNA degradation after sensing invading plasmid, which kills the infected cells and results in the depletion of the invader from the cell population. Together, the BPAN system provides immunoprotection via abortive infection. Our data also suggest that the defense strategy is employed by other long-B pAgos equipped with distinct associated proteins.