Search

The NEREA (Naples Ecological REsearch for Augmented observatories) initiative aims to establish an augmented observatory in the Gulf of Naples (GoN), designed to advance the understanding of marine ecosystems through a holistic approach. Inspired by the Tara Oceans expedition and building on the scientific legacy of the MareChiara Long-Term Ecological Research (LTER-MC) site, NEREA integrates traditional physical, chemical, and biological measurements with state-of-the-art methodologies such as metabarcoding and metagenomics. Here we present the first 10 months of NEREA data, collected from April 2019 to January 2020, encompassing physico-chemical parameters, plankton biodiversity (e.g., microscopy and flow cytometry), prokaryotic and eukaryotic metabarcoding, a prokaryotic gene catalogue, and a collection of 3818 prokaryotic Metagenome-Assembled Genomes (MAGs). NEREA's efforts produce a significant volume of multifaceted data, which enhances our understanding of marine ecosystems and promotes the development of scientific hypotheses and ideas. [ABSTRACT FROM AUTHOR]

Oil reservoirs, being one of the significant subsurface repositories of energy and carbon, host diverse microbial communities affecting energy production and carbon emissions. Viruses play crucial roles in the ecology of microbiomes, however, their distribution and ecological significance in oil reservoirs remain undetermined. Here, we assemble a catalogue encompassing viral and prokaryotic genomes sourced from oil reservoirs. The catalogue comprises 7229 prokaryotic genomes and 3,886 viral Operational Taxonomic Units (vOTUs) from 182 oil reservoir metagenomes. The results show that viruses are widely distributed in oil reservoirs, and 85% vOTUs in oil reservoir are detected in less than 10% of the samples, highlighting the heterogeneous nature of viral communities within oil reservoirs. Through combined microcosm enrichment experiments and bioinformatics analysis, we validate the ecological roles of viruses in regulating the community structure of sulfate reducing microorganisms, primarily through a virulent lifestyle. Taken together, this study uncovers a rich diversity of viruses and their ecological functions within oil reservoirs, offering a comprehensive understanding of the role of viral communities in the biogeochemical cycles of the deep biosphere. Oil reservoirs host diverse microbial communities affecting energy production and carbon emissions. Here, the authors use metagenomics data to construct and characterize a catalogue of viral and prokaryotic genomes from 182 oil reservoirs world-wide, and further validate an ecological role of viruses in regulating the community structure of sulfate reducing microorganisms. [ABSTRACT FROM AUTHOR]

Background: Understanding the interactions and dynamics of microbiotas within biological wastewater treatment systems is essential for ensuring their stability and long-term sustainability. In this study, we developed a systematic framework employing multi-omics and Hi-C sequencing to extensively investigate prokaryotic and phage communities within a hybrid biofilm and activated sludge system. Results: We uncovered distinct distribution patterns, metabolic capabilities, and activities of functional prokaryotes through the analysis of 454 reconstructed prokaryotic genomes. Additionally, we reconstructed a phage catalog comprising 18,645 viral operational taxonomic units (vOTUs) with high length and contiguity using hybrid assembly, and a distinct distribution of phages was depicted between activated sludge (AS) and biofilm. Importantly, 1340 host-phage pairs were established using Hi-C and conventional in silico methods, unveiling the host-determined phage prevalence. The majority of predicted hosts were found to be involved in various crucial metabolic processes, highlighting the potential vital roles of phages in influencing substance metabolism within this system. Moreover, auxiliary metabolic genes (AMGs) related to various categories (e.g., carbohydrate degradation, sulfur metabolism, transporter) were predicted. Subsequent activity analysis emphasized their potential ability to mediate host metabolism during infection. We also profiled the temporal dynamics of phages and their associated hosts using 13-month time-series metagenomic data, further demonstrating their tight interactions. Notably, we observed lineage-specific infection patterns, such as potentially host abundance- or phage/host ratio-driven phage population changes. Conclusions: The insights gained from this research contribute to the growing body of knowledge surrounding interactions and dynamics of host-phage and pave the way for further exploration and potential applications in the field of microbial ecology. 7mTZubWt51_LqpKWPutADa Video Abstract [ABSTRACT FROM AUTHOR]

The knowledge of the different population‐level processes operating within a species, and the genetic variability of the individual prokaryotic genomes, is key to understanding the adaptability of microbial populations. Here, we characterized the flexible genome of ammonia‐oxidizing archaeal (AOA) populations using a metagenomic recruitment approach and long‐read (PacBio HiFi) metagenomic sequencing. In the lower photic zone of the western Mediterranean Sea (75 m deep), the genomes Nitrosopelagicus brevis CN25 and Nitrosopumilus catalinensis SPOT1 had the highest recruitment values among available complete AOA genomes. They were used to analyse the diversity of flexible genes (variable from strain to strain) by examining the long‐reads located within the flexible genomic islands (fGIs) identified by their under‐recruitment. Both AOA genomes had a large fGI involved in the glycosylation of exposed structures, highly variable, and rich in glycosyltransferases. N. brevis had two fGIs related to the transport of phosphorus and ammonium respectively. N. catalinensis had fGIs involved in phosphorus transportation and metal uptake. A fGI5 previously reported as 'unassigned function' in N. brevis could be associated with defense. These findings demonstrate that the microdiversity of marine microbe populations, including AOA, can be effectively characterized using an approach that incorporates third‐generation sequencing metagenomics. [ABSTRACT FROM AUTHOR]

Background: Structural variations play an important role in bacterial genomes. They can mediate genome adaptation quickly in response to the external environment and thus can also play a role in antibiotic resistance. The detection of structural variations in bacteria is challenging, and the recognition of even small rearrangements can be important. Even though most detection tools are aimed at and benchmarked on eukaryotic genomes, they can also be used on prokaryotic genomes. The key features of detection are the ability to detect small rearrangements and support haploid genomes. Because of the limiting performance of a single detection tool, combining the detection abilities of multiple tools can lead to more robust results. There are already available workflows for structural variation detection for long-reads technologies and for the detection of single-nucleotide variation and indels, both aimed at bacteria. Yet we are unaware of structural variations detection workflows for the short-reads sequencing platform. Motivated by this gap we created our workflow. Further, we were interested in increasing the detection performance and providing more robust results. Results: We developed an open-source bioinformatics pipeline, ProcaryaSV, for the detection of structural variations in bacterial isolates from paired-end short sequencing reads. Multiple tools, starting with quality control and trimming of sequencing data, alignment to the reference genome, and multiple structural variation detection tools, are integrated. All the partial results are then processed and merged with an in-house merging algorithm. Compared with a single detection approach, ProcaryaSV has improved detection performance and is a reproducible easy-to-use tool. Conclusions: The ProcaryaSV pipeline provides an integrative approach to structural variation detection from paired-end next-generation sequencing of bacterial samples. It can be easily installed and used on Linux machines. It is publicly available on GitHub at https://github.com/robinjugas/ProcaryaSV. [ABSTRACT FROM AUTHOR]

Although the principles of synthetic biology were initially established in model bacteria, microbial producers, extremophiles and gut microbes have now emerged as valuable prokaryotic chassis for biological engineering. Extending the host range in which designed circuits can function reliably and predictably presents a major challenge for the concept of synthetic biology to materialize. In this work, we systematically characterized the cross‐species universality of two transcriptional regulatory modules—the T7 RNA polymerase activator module and the repressors module—in three non‐model microbes. We found striking linear relationships in circuit activities among different organisms for both modules. Parametrized model fitting revealed host non‐specific parameters defining the universality of both modules. Lastly, a genetic NOT gate and a band‐pass filter circuit were constructed from these modules and tested in non‐model organisms. Combined models employing host non‐specific parameters were successful in quantitatively predicting circuit behaviors, underscoring the potential of universal biological parts and predictive modeling in synthetic bioengineering. [ABSTRACT FROM AUTHOR]

Defense-associated sirtuin 2 (DSR2) systems are widely distributed across prokaryotic genomes, providing robust protection against phage infection. DSR2 recognizes phage tail tube proteins and induces abortive infection by depleting intracellular NAD+, a process that is counteracted by another phage-encoded protein, DSR Anti Defense 1 (DSAD1). Here, we present cryo-EM structures of Bacillus subtilis DSR2 in its apo, Tube-bound, and DSAD1-bound states. DSR2 assembles into an elongated tetramer, with four NADase catalytic modules clustered in the center and the regulatory-sensing modules distributed at four distal corners. Interestingly, monomeric Tube protein, rather than its oligomeric states, docks at each corner of the DSR2 tetramer to form a 4:4 DSR2-Tube assembly, which is essential for DSR2 NADase activity. DSAD1 competes with Tube for binding to DSR2 by occupying an overlapping region, thereby inhibiting DSR2 immunity. Thus, our results provide important insights into the assembly, activation and inhibition of the DSR2 anti-phage defense system. The defense-associated sirtuin 2 (DSR2) system protects bacteria from phages by depleting NAD+. Here, authors elucidate the molecular mechanisms underlying DSR2 assembly, activation, and inhibition, providing important insights into bacterial anti-phage defense. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Background: Brucellosis is one of the most common zoonoses caused by Brucella species. Brucella is an intracellular pathogen that causes abortion in domestic animals and brucellosis in humans. The disease has no specific pathognomonic signs; hence a prompt and accurate diagnosis of the disease plays a very critical role in preventing brucella transmission and treating the disease. Methods: Research was carried out at the Department of Veterinary Science of Kazakh National Agrarian Research University and in the Molecular Biology laboratory of the Institute for Plant Biology and Biotechnology under the Committee of Science of the Ministry of Education and Science of the Republic of Kazakhstan. Studies were carried out on a killed culture of the vaccine strain Brucella abortus RB19, genomic DNA of B. abortus RB19, Grapevine virus A, plasmid vectors for DNA cloning, a fluorescent nano marker, and Brucella membrane proteins. Results: The research demonstrated that brucellosis remained one of the most common zoonotic diseases worldwide. This disease not only resulted in significant economic losses in agriculture but also posed a serious threat to public health. Conclusions: The study successfully established prokaryotic expression vectors for efficient protein production, validated through sequencing and Western blotting, and proposed a promising nano-platform for plant-based synthesis of brucellosis proteins, with potential implications for diagnostic kit development and vaccine research in agriculture. [ABSTRACT FROM AUTHOR]

Background: Lactiplantibacillus plantarum is a renowned probiotic known for its ability to enhance health when ingested in sufficient quantities. While L. plantarum typically exhibits unique characteristics, variations in strains and environments can occur. Therefore, it is crucial to conduct genetic identification of the strain before incorporating it into products or administering it for health benefits. Methods: After being isolated from classic artisanal milk cheeses, the bacterial DNA extraction kit from GENEAID, Korea, was used, and strain HA9's genomic DNA was isolated. Then, the highly conserved section of ribosomal RNA was amplified with universal primers. The results were confirmed by a blast analysis of the NCBI database. The genome of HA9 was sequenced using the Illumina and PacBio Sequel II technologies (HN00194138; Psomagen/USA). The Prokaryotic Dynamic Programming Gene-Finding Algorithm software was applied to forecast coding genes. The Pathosystems Resource Integration Center constructed the L. plantarum genome. Result: According to the findings, the genome of the L. plantarum HA9 strain has a length of 3,321,986 bp and a GC percentage of 44.40. This assembled genome contained a total of 1,420 coding proteins and 82 contigs, as discovered by the Comprehensive Genome Analysis. The isolate had been submitted to the NCBI and received the HA9 genome sequence (accession ID: JAUTDL000000000.1). Conclusion Microbes' whole-genome sequencing plays a crucial role in clarifying their probiotic properties and enables accurate identification down to the strain level. The authors argue that further research should be done on probiotics' whole genomes, with a special emphasis on second-generation strains. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod‐infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral‐encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria‐infecting phages resulted to be significantly CpG‐depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG‐rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications. [ABSTRACT FROM AUTHOR]

Prokaryotic genomes are dynamic tapestries that are strongly influenced by mobile genetic elements (MGEs), including transposons (Tn's), plasmids, and bacteriophages. Of these, miniature inverted-repeat transposable elements (MITEs) are undoubtedly the least studied MGEs in bacteria and archaea. This review explores the diversity and distribution of MITEs in prokaryotes and describes what is known about their functional roles in the host and involvement in genomic plasticity and evolution. [ABSTRACT FROM AUTHOR]

Background: One of the stranger phenomena that can occur during gene translation is where, as a ribosome reads along the mRNA, various cellular and molecular properties contribute to stalling the ribosome on a slippery sequence and shifting the ribosome into one of the other two alternate reading frames. The alternate frame has different codons, so different amino acids are added to the peptide chain. More importantly, the original stop codon is no longer in-frame, so the ribosome can bypass the stop codon and continue to translate the codons past it. This produces a longer version of the protein, a fusion of the original in-frame amino acids, followed by all the alternate frame amino acids. There is currently no automated software to predict the occurrence of these programmed ribosomal frameshifts (PRF), and they are currently only identified by manual curation. Results: Here we present PRFect, an innovative machine-learning method for the detection and prediction of PRFs in coding genes of various types. PRFect combines advanced machine learning techniques with the integration of multiple complex cellular properties, such as secondary structure, codon usage, ribosomal binding site interference, direction, and slippery site motif. Calculating and incorporating these diverse properties posed significant challenges, but through extensive research and development, we have achieved a user-friendly approach. The code for PRFect is freely available, open-source, and can be easily installed via a single command in the terminal. Our comprehensive evaluations on diverse organisms, including bacteria, archaea, and phages, demonstrate PRFect's strong performance, achieving high sensitivity, specificity, and an accuracy exceeding 90%. The code for PRFect is freely available and installs with a single terminal command. Conclusion: PRFect represents a significant advancement in the field of PRF detection and prediction, offering a powerful tool for researchers and scientists to unravel the intricacies of programmed ribosomal frameshifting in coding genes. [ABSTRACT FROM AUTHOR]

The South China Sea (SCS) is a marginal sea characterized by strong land-sea biogeochemical interactions. SCS has a distinctive landscape with a multitude of seamounts in its basin. Seamounts create "seamount effects" that influence the diversity and distribution of planktonic microorganisms in the surrounding oligotrophic waters. Although the vertical distribution and community structure of marine microorganisms have been explored in certain regions of the global ocean, there is a lack of comprehensive microbial genomic surveys for uncultured microorganisms in SCS, particularly in the seamount regions. Here, we employed a metagenomic approach to study the uncultured microbial communities sampled from the Xianbei seamount region to the North Coast waters of SCS. A total of 1887 non-redundant prokaryotic metagenome-assembled genomes (MAGs) were reconstructed, of which, 153 MAGs were classified as high-quality MAGs based on the MIMAG standards. The community structure and genomic information provided by this dataset could be used to analyze microbial distribution and metabolism in the SCS. [ABSTRACT FROM AUTHOR]

Toxin-antitoxin systems (TAs) are generally two-component genetic modules present in almost every prokaryotic genome. The production of the free and active toxin is able to disrupt key cellular processes leading to the growth inhibition or death of its host organism in absence of its cognate antitoxin. The functions attributed to TAs rely on this lethal phenotype ranging from mobile genetic elements stabilization to phage defense. Their abundance in prokaryotic genomes as well as their lethal potential make them attractive targets for new antibacterial strategies. The hijacking of TAs requires a deep understanding of their regulation to be able to design such approach. In this review, we summarize the accumulated knowledge on how bacteria cope with these toxic genes in their genome. The characterized TAs can be grouped based on the way they prevent toxicity. Some systems rely on a tight control of the expression to prevent the production of the toxin while others control the activity of the toxin at the post-translational level. [ABSTRACT FROM AUTHOR]

Background: The anaerobic digestion process degrades organic matter into simpler compounds and occurs in strictly anaerobic and microaerophilic environments. The process is carried out by a diverse community of microorganisms where each species has a unique role and it has relevant biotechnological applications since it is used for biogas production. Some aspects of the microbiome, including its interaction with phages, remains still unclear: a better comprehension of the community composition and role of each species is crucial for a cured understanding of the carbon cycle in anaerobic systems and improving biogas production. Results: The primary objective of this study was to expand our understanding on the anaerobic digestion microbiome by jointly analyzing its prokaryotic and viral components. By integrating 192 additional datasets into a previous metagenomic database, the binning process generated 11,831 metagenome-assembled genomes from 314 metagenome samples published between 2014 and 2022, belonging to 4,568 non-redundant species based on ANI calculation and quality verification. CRISPR analysis on these genomes identified 76 archaeal genomes with active phage interactions. Moreover, single-nucleotide variants further pointed to archaea as the most critical members of the community. Among the MAGs, two methanogenic archaea, Methanothrix sp. 43zhSC_152 and Methanoculleus sp. 52maCN_3230, had the highest number of SNVs, with the latter having almost double the density of most other MAGs. Conclusions: This study offers a more comprehensive understanding of microbial community structures that thrive at different temperatures. The findings revealed that the fraction of archaeal species characterized at the genome level and reported in public databases is higher than that of bacteria, although still quite limited. The identification of shared spacers between phages and microbes implies a history of phage-bacterial interactions, and specifically lysogenic infections. A significant number of SNVs were identified, primarily comprising synonymous and nonsynonymous variants. Together, the findings indicate that methanogenic archaea are subject to intense selective pressure and suggest that genomic variants play a critical role in the anaerobic digestion process. Overall, this study provides a more balanced and diverse representation of the anaerobic digestion microbiota in terms of geographic location, temperature range and feedstock utilization. [ABSTRACT FROM AUTHOR]

Simple Summary: Transposable elements (TEs) are mobile genetic elements that play vital role in defining and contributing to the size, shape and structure of both prokaryotic and eukaryotic genomes in nature. The piggyBac (PB), a superfamily of DNA transposons, has been isolated over the years from various organisms including insects, fungi and plants. These piggyBac transposon systems have high efficiency with a wide usage in the study of gene therapies, mutagenesis and transgenesis. Currently, there is limited information available on DNA transposons in small (compact) genomes of animals. Therefore, this study aims to annotate the PB transposons in small genomes of animals, revealing their evolution profiles in both vertebrate and invertebrate genomes. TEs, including DNA transposons, are major contributors of genome expansions, and have played a very significant role in shaping the evolution of animal genomes, due to their capacity to jump from one genomic position to the other. In this study, we investigated the evolution landscapes of PB transposons, including their distribution, diversity, activity and structure organization in 79 species of small (compact) genomes of animals comprising both vertebrate and invertebrates. Overall, 212 PB transposon types were detected from almost half (37) of the total number of the small genome species (79) investigated. The detected PB transposon types, which were unevenly distributed in various genera and phyla, have been classified into seven distinct clades or families with good bootstrap support (>80%). The PB transposon types that were identified have a length ranging from 1.23 kb to 9.51 kb. They encode transposases of approximately ≥500 amino acids in length, and possess terminal inverted repeats (TIRs) ranging from 4 bp to 24 bp. Though some of the transposon types have long TIRs (528 bp), they still maintain the consistent and reliable 4 bp target site duplication (TSD) of TTAA. However, PiggyBac-2_Cvir transposon originating from the Crassostrea virginica species exhibits a unique TSD of TATG. The TIRs of the transposons in all the seven families display high divergence, with a highly conserved 5′ end motif. The core transposase domains (DDD) were better conserved among the seven different families compared to the other protein domains, which were less prevalent in the vertebrate genome. The divergent evolution dynamics analysis also indicated that the majority of the PB transposon types identified in this study are either relatively young or old, with some being active. Additionally, numerous invasions of PB transposons were found in the genomes of both vertebrate and invertebrate animals. The data reveals that the PB superfamily is widely distributed in these species. PB transposons exhibit high diversity and activity in the small genomes of animals, and might play a crucial role in shaping the evolution of these small genomes of animals. [ABSTRACT FROM AUTHOR]

Aim: Glycoside hydrolases (GHs) encoded by prokaryotes play a crucial role in the degradation of carbohydrates, which is a key step in the carbon cycle. However, their global distribution remains poorly understood. Here, we conducted a comprehensive investigation of the global biogeographic patterns of prokaryotic‐encoded GHs. Location: Global. Time period: Present. Major taxa studied: Prokaryotes. Methods: We conducted a global survey of GH distribution based on 10,000 samples from the Earth Microbiome Project (EMP) and 17,923 completely sequenced prokaryotic genomes. The samples were retrieved from seven continents, covering all 17 EMP Ontology environments. We evaluated the assembly process and influencing factors of GH functional genes within prokaryotic communities, and analysed the rank abundance distributions (RADs), occupancy frequency distributions (OFDs), and taxonomic sources of GH families across various environments. Results: The median number of GH genes per thousand genes encoded by global prokaryotic communities was 9.11 (7.16–11.58), with a higher GH abundance in host‐associated communities than in free‐living communities. The median number of GH families in prokaryotic communities was 105 (86–121), with fewer families but tighter functional connections between families in host‐associated communities compared to free‐living communities. In any environment, the RAD of GH families followed a lognormal model, and their OFD was bimodal, with a strong positive correlation between the two. In addition, most GH families had complex and diverse taxonomic sources across different environments. Main conclusions: This study provides a panoramic view of the functional potentials for carbohydrate degradation across prokaryotic communities under diverse environmental conditions. The composition of GHs encoded by prokaryotic communities exhibits clear habitat specificity, which is driven by a combination of biotic and abiotic factors and generated through highly deterministic assembly processes. [ABSTRACT FROM AUTHOR]

The guanine/cytosine (GC) content of prokaryotic genomes is species‐specific, taking values from 16% to 77%. This diversity of selection for GC content remains contentious. We analyse the correlations between GC content and a range of phenotypic and genotypic data in thousands of prokaryotes. GC content integrates well with these traits into r/K selection theory when phenotypic plasticity is considered. High GC‐content prokaryotes are r‐strategists with cheaper descendants thanks to a lower average amino acid metabolic cost, colonize unstable environments thanks to flagella and a bacillus form and are generalists in terms of resource opportunism and their defence mechanisms. Low GC content prokaryotes are K‐strategists specialized for stable environments that maintain homeostasis via a high‐cost outer cell membrane and endospore formation as a response to nutrient deprivation, and attain a higher nutrient‐to‐biomass yield. The lower proteome cost of high GC content prokaryotes is driven by the association between GC‐rich codons and cheaper amino acids in the genetic code, while the correlation between GC content and genome size may be partly due to functional diversity driven by r/K selection. In all, molecular diversity in the GC content of prokaryotes may be a consequence of ecological r/K selection. [ABSTRACT FROM AUTHOR]

We define two new computational problems in the domain of perfect genome rearrangements, and propose three algorithms to solve them. The rearrangement scenarios modeled by the problems consider Reversal and Block Interchange operations, and a PQ-tree is utilized to guide the allowed operations and to compute their weights. In the first problem, Constrained TreeToString Divergence ( CTTSD ), we define the basic structure-informed rearrangement measure. Here, we assume that the gene order members of the gene cluster from which the PQ-tree is constructed are permutations. The PQ-tree representing the gene cluster is ordered such that the series of gene IDs spelled by its leaves is equivalent to that of the reference gene order. Then, a structure-informed genome rearrangement distance is computed between the ordered PQ-tree and the target gene order. The second problem, TreeToString Divergence ( TTSD ), generalizes CTTSD , where the gene order members are not necessarily permutations and the structure informed rearrangement measure is extended to also consider up to d S and d T gene insertion and deletion operations, respectively, when modelling the PQ-tree informed divergence process from the reference gene order to the target gene order. The first algorithm solves CTTSD in O (n γ 2 · (m p · 1. 381 γ + m q)) time and O (n 2) space, where γ is the maximum number of children of a node, n is the length of the string and the number of leaves in the tree, and m p and m q are the number of P-nodes and Q-nodes in the tree, respectively. If one of the penalties of CTTSD is 0, then the algorithm runs in O (n m γ 2) time and O (n 2) space. The second algorithm solves TTSD in O (n 2 γ 2 d T 2 d S 2 m 2 (m p · 5 γ γ + m q)) time and O (d T d S m (m n + 5 γ)) space, where γ is the maximum number of children of a node, n is the length of the string, m is the number of leaves in the tree, m p and m q are the number of P-nodes and Q-nodes in the tree, respectively, and allowing up to d T deletions from the tree and up to d S deletions from the string. The third algorithm is intended to reduce the space complexity of the second algorithm. It solves a variant of the problem (where one of the penalties of TTSD is 0) in O (n γ 2 d T 2 d S 2 m 2 (m p · 4 γ γ 2 n (d T + d S + m + n) + m q)) time and O (γ 2 n m 2 d T d S (d T + d S + m + n)) space. The algorithm is implemented as a software tool, denoted MEM-Rearrange, and applied to the comparative and evolutionary analysis of 59 chromosomal gene clusters extracted from a dataset of 1487 prokaryotic genomes. [ABSTRACT FROM AUTHOR]

Background: The RCA (Roseobacter clade affiliated) cluster belongs to the family Roseobacteracea and represents a major Roseobacter lineage in temperate to polar oceans. Despite its prevalence and abundance, only a few genomes and one described species, Planktomarina temperata, exist. To gain more insights into our limited understanding of this cluster and its taxonomic and functional diversity and biogeography, we screened metagenomic datasets from the global oceans and reconstructed metagenome-assembled genomes (MAG) affiliated to this cluster. Results: The total of 82 MAGs, plus five genomes of isolates, reveal an unexpected diversity and novel insights into the genomic features, the functional diversity, and greatly refined biogeographic patterns of the RCA cluster. This cluster is subdivided into three genera: Planktomarina, Pseudoplanktomarina, and the most deeply branching Candidatus Paraplanktomarina. Six of the eight Planktomarina species have larger genome sizes (2.44–3.12 Mbp) and higher G + C contents (46.36–53.70%) than the four Pseudoplanktomarina species (2.26–2.72 Mbp, 42.22–43.72 G + C%). Cand. Paraplanktomarina is represented only by one species with a genome size of 2.40 Mbp and a G + C content of 45.85%. Three novel species of the genera Planktomarina and Pseudoplanktomarina are validly described according to the SeqCode nomenclature for prokaryotic genomes. Aerobic anoxygenic photosynthesis (AAP) is encoded in three Planktomarina species. Unexpectedly, proteorhodopsin (PR) is encoded in the other Planktomarina and all Pseudoplanktomarina species, suggesting that this light-driven proton pump is the most important mode of acquiring complementary energy of the RCA cluster. The Pseudoplanktomarina species exhibit differences in functional traits compared to Planktomarina species and adaptations to more resource-limited conditions. An assessment of the global biogeography of the different species greatly expands the range of occurrence and shows that the different species exhibit distinct biogeographic patterns. They partially reflect the genomic features of the species. Conclusions: Our detailed MAG-based analyses shed new light on the diversification, environmental adaptation, and global biogeography of a major lineage of pelagic bacteria. The taxonomic delineation and validation by the SeqCode nomenclature of prominent genera and species of the RCA cluster may be a promising way for a refined taxonomic identification of major prokaryotic lineages and sublineages in marine and other prokaryotic communities assessed by metagenomics approaches. 4nWA41D1uUohyPVze8XWvD Video Abstract [ABSTRACT FROM AUTHOR]

Prokaryotic genomes are generally organized in haploid. In synthetic biological research, efficient chassis cells must be constructed to produce bio‐based products. Here, the essential division of the ftsZ gene to create functional polyploid E. coli is regulated. The artificial polyploid E. coli containing 2–4 chromosomes is confirmed through PCR amplification, terminator localization, and flow cytometry. The polyploid E. coli exhibits a larger cell size, and its low pH tolerance and acetate resistance are stronger than those of haploid E. coli. Transcriptome analysis shows that the genes of the cell's main functional pathways are significantly upregulated in the polyploid E. coli. These advantages of the polyploid E. coli results in the highest reported L‐threonine yield (160.3 g L−1) in fed‐batch fermentation to date. In summary, an easy and convenient method for constructing polyploid E. coli and demonstrated its application in L‐threonine production is developed. This work provides a new approach for creating an excellent host strain for biochemical production and studying the evolution of prokaryotes and their chromosome functions. [ABSTRACT FROM AUTHOR]

Microbial genome recovery from metagenomes can further explain microbial ecosystem structures, functions and dynamics. Thus, this study developed the Additional Clustering Refiner (ACR) to enhance high-purity prokaryotic and eukaryotic metagenome-assembled genome (MAGs) recovery. ACR refines low-quality MAGs by subjecting them to iterative k -means clustering predicated on contig abundance and increasing bin purity through validated universal marker genes. Synthetic and real-world metagenomic datasets, including short- and long-read sequences, evaluated ACR's effectiveness. The results demonstrated improved MAG purity and a significant increase in high- and medium-quality MAG recovery rates. In addition, ACR seamlessly integrates with various binning algorithms, augmenting their strengths without modifying core features. Furthermore, its multiple sequencing technology compatibilities expand its applicability. By efficiently recovering high-quality prokaryotic and eukaryotic genomes, ACR is a promising tool for deepening our understanding of microbial communities through genome-centric metagenomics. [ABSTRACT FROM AUTHOR]

The human gastrointestinal (GI) tract harbors eukaryotic and prokaryotic viruses and their genomes, metabolites, and proteins, collectively known as the "gut virome". This complex community of viruses colonizing the enteric mucosa is pivotal in regulating host immunity. The mechanisms involved in cross communication between mucosal immunity and the gut virome, as well as their relationship in health and disease, remain largely unknown. Herein, we review the literature on the human gut virome's composition and evolution and the interplay between the gut virome and enteric mucosal immunity and their molecular mechanisms. Our review suggests that future research efforts should focus on unraveling the mechanisms of gut viruses in human homeostasis and pathophysiology and on developing virus-prompted precision therapies. [ABSTRACT FROM AUTHOR]

The ever-growing repertoire of genomic techniques continues to expand our understanding of the true diversity and richness of prokaryotic genomes. Riboproteogenomics laid the foundation for dynamic studies of previously overlooked genomic elements. Most strikingly, bacterial genomes were revealed to harbor robust repertoires of small open reading frames (sORFs) encoding a diverse and broadly expressed range of small proteins, or sORF-encoded polypeptides (SEPs). In recent years, continuous efforts led to great improvements in the annotation and characterization of such proteins, yet many challenges remain to fully comprehend the pervasive nature of small proteins and their impact on bacterial biology. In this work, we review the recent developments in the dynamic field of bacterial genome reannotation, catalog the important biological roles carried out by small proteins and identify challenges obstructing the way to full understanding of these elusive proteins. [ABSTRACT FROM AUTHOR]

The global sulphur cycle has implications for human health, climate change, biogeochemistry and bioremediation. The organosulphur compounds that participate in this cycle not only represent a vast reservoir of sulphur but are also used by prokaryotes as sources of energy and/or carbon. Closely linked to the inorganic sulphur cycle, it involves the interaction of prokaryotes, eukaryotes and chemical processes. However, ecological and evolutionary studies of the conversion of organic sulphur compounds are hampered by the poor conservation of the relevant pathways and their variation even within strains of the same species. In addition, several proteins involved in the conversion of sulphonated compounds are related to proteins involved in sulphur dissimilation or turnover of other compounds. Therefore, the enzymes involved in the metabolism of organic sulphur compounds are usually not correctly annotated in public databases. To address this challenge, we have developed HMSS2, a profiled Hidden Markov Model‐based tool for rapid annotation and synteny analysis of organic and inorganic sulphur cycle proteins in prokaryotic genomes. Compared to its previous version (HMS‐S‐S), HMSS2 includes several new features. HMM‐based annotation is now supported by nonhomology criteria and covers the metabolic pathways of important organosulphur compounds, including dimethylsulphoniopropionate, taurine, isethionate, and sulphoquinovose. In addition, the calculation speed has been increased by a factor of four and the available output formats have been extended to include iTol compatible data sets, and customized sequence FASTA files. [ABSTRACT FROM AUTHOR]

Many gut microorganisms critical to human health rely on nutrients produced by each other for survival; however, these cross-feeding interactions are still challenging to quantify and remain poorly characterized. Here, we introduce a Metabolite Exchange Score (MES) to quantify those interactions. Using metabolic models of prokaryotic metagenome-assembled genomes from over 1600 individuals, MES allows us to identify and rank metabolic interactions that are significantly affected by a loss of cross-feeding partners in 10 out of 11 diseases. When applied to a Crohn's disease case-control study, our approach identifies a lack of species with the ability to consume hydrogen sulfide as the main distinguishing microbiome feature of disease. We propose that our conceptual framework will help prioritize in-depth analyses, experiments and clinical targets, and that targeting the restoration of microbial cross-feeding interactions is a promising mechanism-informed strategy to reconstruct a healthy gut ecosystem. Gut microbes rely on nutrient exchange for survival, but these cross-feeding interactions remain poorly characterized. Here, Marcelino et al. present a metabolite-exchange scoring system derived from metagenome-scale metabolic models, designed to identify the potential microbial cross-feeding interactions most affected in human diseases. [ABSTRACT FROM AUTHOR]

Mobile genetic elements (MGEs) are relevant agents in bacterial adaptation and evolutionary diversification. Stable appropriation of these DNA elements depends on host factors, among which are the nucleoid-associated proteins (NAPs). NAPs are highly abundant proteins that bind and bend DNA, altering its topology and folding, thus affecting all known cellular DNA processes from replication to expression. Even though NAP coding genes are found in most prokaryotic genomes, their functions in host chromosome biology and xenogeneic silencing are only known for a few NAP families. Less is known about the occurrence, abundance, and roles of MGE-encoded NAPs in foreign elements establishment and mobility. In this study, we used a combination of comparative genomics and phylogenetic strategies to gain insights into the diversity, distribution, and functional roles of NAPs within the class Acidithiobacillia with a special focus on their role in MGE biology. Acidithiobacillia class members are aerobic, chemolithoautotrophic, acidophilic sulfur-oxidizers, encompassing substantial genotypic diversity attributable to MGEs. Our search for NAP protein families (PFs) in more than 90 genomes of the different species that conform the class, revealed the presence of 1,197 proteins pertaining to 12 different NAP families, with differential occurrence and conservation across species. Pangenome-level analysis revealed 6 core NAP PFs that were highly conserved across the class, some of which also existed as variant forms of scattered occurrence, in addition to NAPs of taxa-restricted distribution. Core NAPs identified are reckoned as essential based on the conservation of genomic context and phylogenetic signals. In turn, various highly diversified NAPs pertaining to the flexible gene complement of the class, were found to be encoded in known plasmids or, larger integrated MGEs or, present in genomic loci associated with MGE-hallmark genes, pointing to their role in the stabilization/maintenance of these elements in strains and species with larger genomes. Both core and flexible NAPs identified proved valuable as markers, the former accurately recapitulating the phylogeny of the class, and the later, as seed in the bioinformatic identification of novel episomal and integrated mobile elements. [ABSTRACT FROM AUTHOR]

Type-V-F Cas12f proteins, also known as Cas14, have drawn significant interest within the diverse CRISPR-Cas nucleases due to their compact size. This study involves analyzing and comparing Cas14-homology proteins in prokaryotic genomes through mining, sequence comparisons, a phylogenetic analysis, and an array/repeat analysis. In our analysis, we identified and mined a total of 93 Cas14-homology proteins that ranged in size from 344 aa to 843 aa. The majority of the Cas14-homology proteins discovered in this analysis were found within the Firmicutes group, which contained 37 species, representing 42% of all the Cas14-homology proteins identified. In archaea, the DPANN group had the highest number of species containing Cas14-homology proteins, a total of three species. The phylogenetic analysis results demonstrate the division of Cas14-homology proteins into three clades: Cas14-A, Cas14-B, and Cas14-U. Extensive similarity was observed at the C-terminal end (CTD) through a domain comparison of the three clades, suggesting a potentially shared mechanism of action due to the presence of cutting domains in that region. Additionally, a sequence similarity analysis of all the identified Cas14 sequences indicated a low level of similarity (18%) between the protein variants. The analysis of repeats/arrays in the extended nucleotide sequences of the identified Cas14-homology proteins highlighted that 44 out of the total mined proteins possessed CRISPR-associated repeats, with 20 of them being specific to Cas14. Our study contributes to the increased understanding of Cas14 proteins across prokaryotic genomes. These homologous proteins have the potential for future applications in the mining and engineering of Cas14 proteins. [ABSTRACT FROM AUTHOR]

Prokaryotes, the earliest forms of life on Earth, play crucial roles in global biogeochemical processes in virtually all ecosystems. The ever-increasing amount of prokaryotic genome sequencing data provides a wealth of information to examine fundamental and applied questions through systematic genome comparison. Genomic features, such as genome size and GC content, and taxonomy-centric genomic features of complete prokaryotic genomes (CPGs) are crucial for various fields of microbial research and education, yet they are often overlooked. Additionally, creating systematically curated datasets that align with research concerns is an essential yet challenging task for wet-lab researchers. In this study, we introduce CPGminer, a user-friendly tool that allows researchers to quickly and easily examine the genomic features and taxonomy of CPGs and curate genome datasets. We also provide several examples to demonstrate its practical utility in addressing descriptive questions. [ABSTRACT FROM AUTHOR]

The availability of an ever-increasing diversity of prokaryotic genomes and metagenomes represents a major opportunity to understand and decipher the mechanisms behind the functional diversification of microbial biosynthetic pathways. However, it remains unclear to what extent a pathway producing a specific molecule from a specific precursor can diversify. In this study, we focus on the biosynthesis of ubiquinone (UQ), a crucial coenzyme that is central to the bioenergetics and to the functioning of a wide variety of enzymes in Eukarya and Pseudomonadota (a subgroup of the formerly named Proteobacteria). UQ biosynthesis involves three hydroxylation reactions on contiguous carbon atoms. We and others have previously shown that these reactions are catalyzed by different sets of UQ-hydroxylases that belong either to the iron-dependent Coq7 family or to the more widespread flavin monooxygenase (FMO) family. Here, we combine an experimental approach with comparative genomics and phylogenetics to reveal how UQ-hydroxylases evolved different selectivities within the constrained framework of the UQ pathway. It is shown that the UQ-FMOs diversified via at least three duplication events associated with two cases of neofunctionalization and one case of subfunctionalization, leading to six subfamilies with distinct hydroxylation selectivity. We also demonstrate multiple transfers of the UbiM enzyme and the convergent evolution of UQ-FMOs toward the same function, which resulted in two independent losses of the Coq7 ancestral enzyme. Diversification of this crucial biosynthetic pathway has therefore occurred via a combination of parallel evolution, gene duplications, transfers, and losses. [ABSTRACT FROM AUTHOR]

Tuberculosis, caused by Mycobacterium tuberculosis, is a lethal infectious disease of significant public health concern. The rise of multidrug-resistant and drug-tolerant strains has necessitated novel approaches to combat the disease. Toxin–antitoxin (TA) systems, key players in bacterial adaptive responses, are prevalent in prokaryotic genomes and have been linked to tuberculosis. The genome of M. tuberculosis strains harbors an unusually high number of TA systems, prompting questions about their biological roles. The VapBC family, a representative type II TA system, is characterized by the VapC toxin, featuring a PilT N-terminal domain with nuclease activity. Its counterpart, VapB, functions as an antitoxin, inhibiting VapC's activity. Additionally, we explore peptide mimics designed to replicate protein helical structures in this review. Investigating these synthetic peptides offers fresh insights into molecular interactions, potentially leading to therapeutic applications. These synthetic peptides show promise as versatile tools for modulating cellular processes and protein–protein interactions. We examine the rational design strategies employed to mimic helical motifs, their biophysical properties, and potential applications in drug development and bioengineering. This review aims to provide an in-depth understanding of TA systems by introducing known complex structures, with a focus on both structural aspects and functional and molecular details associated with each system. [ABSTRACT FROM AUTHOR]

Toxin–antitoxin (TA) systems are entities found in the prokaryotic genomes, with eight reported types. Type II, the best characterized, is comprised of two genes organized as an operon. Whereas toxins impair growth, the cognate antitoxin neutralizes its activity. TAs appeared to be involved in plasmid maintenance, persistence, virulence, and defence against bacteriophages. Most Type II toxins target the bacterial translational machinery. They seem to be antecessors of Higher Eukaryotes and Prokaryotes Nucleotide-binding (HEPN) RNases, minimal nucleotidyltransferase domains, or CRISPR–Cas systems. A total of four TAs encoded by Streptococcus pneumoniae , RelBE, YefMYoeB, Phd-Doc, and HicAB, belong to HEPN-RNases. The fifth is represented by PezAT/Epsilon–Zeta. PezT/Zeta toxins phosphorylate the peptidoglycan precursors, thereby blocking cell wall synthesis. We explore the body of knowledge (facts) and hypotheses procured for Type II TAs and analyse the data accumulated on the PezAT family. Bioinformatics analyses showed that homologues of PezT/Zeta toxin are abundantly distributed among 14 bacterial phyla mostly in Proteobacteria (48%), Firmicutes (27%), and Actinobacteria (18%), showing the widespread distribution of this TA. The pezAT locus was found to be mainly chromosomally encoded whereas its homologue, the tripartite omega–epsilon–zeta locus, was found mostly on plasmids. We found several orphan pezT / zeta toxins, unaccompanied by a cognate antitoxin. [ABSTRACT FROM AUTHOR]

CRISPR-Cas immunity systems safeguard prokaryotic genomes by inhibiting the invasion of mobile genetic elements. Here, we screened prokaryotic genomic sequences and identified multiple natural transpositions of insertion sequences (ISs) into cas genes, thus inactivating CRISPR-Cas defenses. We then generated an IS-trapping system, using Escherichia coli strains with various ISs and an inducible cas nuclease, to monitor IS insertions into cas genes following the induction of double-strand DNA breakage as a physiological host stress. We identified multiple events mediated by different ISs, especially IS1 and IS10, displaying substantial relaxed target specificity. IS transposition into cas was maintained in the presence of DNA repair machinery, and transposition into other host defense systems was also detected. Our findings highlight the potential of ISs to counter CRISPR activity, thus increasing bacterial susceptibility to foreign DNA invasion. CRISPR-Cas immunity systems safeguard prokaryotic genomes by inhibiting the invasion of mobile genetic elements. Here, the authors show that insertion sequences can efficiently insert into cas genes, thus inactivating CRISPR defenses and increasing bacterial susceptibility to foreign DNA invasion. [ABSTRACT FROM AUTHOR]

As most eukaryotic genomes are yet to be sequenced, the mechanisms underlying their contribution to different ecosystem processes remain untapped. Although approaches to recovering Prokaryotic genomes have become common in genome biology, few studies have tackled the recovery of eukaryotic genomes from metagenomes. This study assessed the reconstruction of microbial eukaryotic genomes using 6000 metagenomes from terrestrial and some transition environments using the EukRep pipeline. Only 215 metagenomic libraries yielded eukaryotic bins. From a total of 447 eukaryotic bins recovered 197 were classified at the phylum level. Streptophytes and fungi were the most represented clades with 83 and 73 bins, respectively. More than 78% of the obtained eukaryotic bins were recovered from samples whose biomes were classified as host‐associated, aquatic, and anthropogenic terrestrial. However, only 93 bins were taxonomically assigned at the genus level and 17 bins at the species level. Completeness and contamination estimates were obtained for a total of 193 bins and consisted of 44.64% (σ = 27.41%) and 3.97% (σ = 6.53%), respectively. Micromonas commoda was the most frequent taxon found while Saccharomyces cerevisiae presented the highest completeness, probably because more reference genomes are available. Current measures of completeness are based on the presence of single‐copy genes. However, mapping of the contigs from the recovered eukaryotic bins to the chromosomes of the reference genomes showed many gaps, suggesting that completeness measures should also include chromosome coverage. Recovering eukaryotic genomes will benefit significantly from long‐read sequencing, development of tools for dealing with repeat‐rich genomes, and improved reference genomes databases. [ABSTRACT FROM AUTHOR]

Background: Marine prokaryotes are a rich source of novel bioactive secondary metabolites for drug discovery. Recent genome mining studies have revealed their great potential to bio-synthesize novel secondary metabolites. However, the exact biosynthetic chemical space encoded by the marine prokaryotes has yet to be systematically evaluated. Results: We first investigated the secondary metabolic potential of marine prokaryotes by analyzing the diversity and novelty of the biosynthetic gene clusters (BGCs) in 7541 prokaryotic genomes from cultivated and single cells, along with 26,363 newly assembled medium-to-high-quality genomes from marine environmental samples. To quantitatively evaluate the unexplored biosynthetic chemical space of marine prokaryotes, the clustering thresholds for constructing the biosynthetic gene cluster and molecular networks were optimized to reach a similar level of the chemical similarity between the gene cluster family (GCF)-encoded metabolites and molecular family (MF) scaffolds using the MIBiG database. The global genome mining analysis demonstrated that the predicted 70,011 BGCs were organized into 24,536 mostly new (99.5%) GCFs, while the reported marine prokaryotic natural products were only classified into 778 MFs at the optimized clustering thresholds. The number of MF scaffolds is only 3.2% of the number of GCF-encoded scaffolds, suggesting that at least 96.8% of the secondary metabolic potential in marine prokaryotes is untapped. The unexplored biosynthetic chemical space of marine prokaryotes was illustrated by the 88 potential novel antimicrobial peptides encoded by ribosomally synthesized and post-translationally modified peptide BGCs. Furthermore, a sea-water-derived Aquimarina strain was selected to illustrate the diverse biosynthetic chemical space through untargeted metabolomics and genomics approaches, which identified the potential biosynthetic pathways of a group of novel polyketides and two known compounds (didemnilactone B and macrolactin A 15-ketone). Conclusions: The present bioinformatics and cheminformatics analyses highlight the promising potential to explore the biosynthetic chemical diversity of marine prokaryotes and provide valuable knowledge for the targeted discovery and biosynthesis of novel marine prokaryotic natural products. 8qkWfGqcrPbn_o7-XLBWPj Video Abstract [ABSTRACT FROM AUTHOR]

Simple Summary: Transposons are mobile genetic elements that can mobilize to other locations in the host genomes. Generally, the mapping of natural transposons is a meticulous endeavor, since the repetitive nature of these elements impedes on the accurate localization of the insertions. In this study, we outlined the landscape of the DNA transposons of a natural population of Drosophila melanogaster from Romania. Harnessing a set of bioinformatics tools designed for transposon mapping, we identified multiple insertions affecting genes with a potential adaptive role. One focus of our work was a detailed analysis of a recently invading transposon, known as the P-element. Another target was to map the DNA natural transposons located in various heterochromatic regions, known as preferential targets for insertions. Our research contributes to a better understanding of the dynamics of DNA transposons and of their roles in the structure and evolution of the D. melanogaster genome. Natural transposons (NTs) represent mobile DNA sequences found in both prokaryotic and eukaryotic genomes. Drosophila melanogaster (the fruit fly) is a eukaryotic model organism with NTs standing for about 20% of its genome and has contributed significantly to the understanding of various aspects of transposon biology. Our study describes an accurate approach designed to map class II transposons (DNA transposons) in the genome of the Horezu_LaPeri fruit fly strain, consecutive to Oxford Nanopore Technology sequencing. A whole genome bioinformatics analysis was conducted using Genome ARTIST_v2, LoRTE and RepeatMasker tools to identify DNA transposons insertions. Then, a gene ontology enrichment analysis was performed in order to evaluate the potential adaptive role of some DNA transposons insertions. Herein, we describe DNA transposon insertions specific for the Horezu_LaPeri genome and a predictive functional analysis of some insertional alleles. The PCR validation of P-element insertions specific for this fruit fly strain, along with a putative consensus sequence for the KP element, is also reported. Overall, the genome of the Horezu_LaPeri strain contains several insertions of DNA transposons associated with genes known to be involved in adaptive processes. For some of these genes, insertional alleles obtained via mobilization of the artificial transposons were previously reported. This is a very alluring aspect, as it suggests that insertional mutagenesis experiments conducting adaptive predictions for laboratory strains may be confirmed by mirroring insertions which are expected to be found at least in some natural fruit fly strains. [ABSTRACT FROM AUTHOR]

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the “dark matter” of the eukaryotic genome and discuss some of the outstanding questions in the field. [ABSTRACT FROM AUTHOR]

One of the main characteristics of prokaryotic genomes is the ratio in which guanine-cytosine bases are used in their DNA sequences. This is known as the genomic GC content and varies widely, from values below 20% to values greater than 74%. It has been demonstrated that the genomic GC content varies in accordance with the phylogenetic distribution of organisms and influences the amino acid composition of their corresponding proteomes. This bias is particularly important for amino acids that are coded by GC content-rich codons such as alanine, glycine, and proline, as well as amino acids that are coded by AT-rich codons, such as lysine, asparagine, and isoleucine. In our study, we extend these results by considering the effect of the genomic GC content on the secondary structure of proteins. On a set of 192 representative prokaryotic genomes and proteome sequences, we identified through a bioinformatic study that the composition of the secondary structures of the proteomes varies in relation to the genomic GC content; random coils increase as the genomic GC content increases, while alpha-helices and beta-sheets present an inverse relationship. In addition, we found that the tendency of an amino acid to form part of a secondary structure of proteins is not ubiquitous, as previously expected, but varies according to the genomic GC content. Finally, we discovered that for some specific groups of orthologous proteins, the GC content of genes biases the composition of secondary structures of the proteins for which they code. [ABSTRACT FROM AUTHOR]

Small open reading frames (smORFs) encoding proteins less than 100 amino acids (aa) are known to be important regulators of key cellular processes. However, their computational identification remains a challenge. Based on a comprehensive analysis of known prokaryotic small ORFs, we have developed the ProsmORF-pred resource which uses a machine learning (ML)-based method for prediction of smORFs in the prokaryotic genome sequences. ProsmORF-pred consists of two ML models, one for initiation site recognition in nucleic acid sequences upstream of putative start codons and the other uses translated amino acid sequences to decipher functional protein like sequences. The nucleotide sequence-based initiation site recognition model has been trained using longer ORFs (>100 aa) in the same genome while the ML model for identification of protein like sequences has been trained using annotated smORFs from Escherichia coli. Comprehensive benchmarking of ProsmORF-pred reveals that its performance is comparable to other state-of-the-art approaches on the annotated smORF set derived from 32 prokaryotic genomes. Its performance is distinctly superior to other tools like PRODIGAL and RANSEPS for prediction of newly identified smORFs which have a length range of 10–30 aa, where prediction of smORFs has been a major challenge. Apart from identification of smORFs in genomic sequences, ProsmORF-pred can also aid in functional annotation of the predicted smORFs based on sequence similarity and genomic neighbourhood similarity searches in ProsmORFDB, a well-curated database of known smORFs. ProsmORF-pred along with its backend database ProsmORFDB is available as a user-friendly web server (http://www.nii.ac.in/prosmorfpred.html). [ABSTRACT FROM AUTHOR]

The discovery of cytidine deaminases as part of toxin–antitoxin systems in bacteria suggests that self-DNA targeting cytidine deaminases of vertebrate immune systems have prokaryotic origins. Comparative analyses of RAG-like transposases of invertebrates and vertebrate RAG recombinase suggest a molecular basis for the loss of integrase activity during the domestication process. The lymphocyte lineage is uniquely dependent on DNA maintenance methylase, which represents an orphan descendant of prokaryotic restriction/methylation systems. Conflicts between genetic parasites and host influence the antigen-recognition properties of adaptive immune systems in vertebrates. Higher order conflicts at the cellular level arise from unscheduled and erroneous activities of immune-related host genome editors. Domestication of ancient toxin/antitoxin systems may have supported the developmental stability of the lymphoid lineage and the emergence of somatic diversification mechanisms of antigen receptor genes. Genetic conflicts shape the genomes of prokaryotic and eukaryotic organisms. Here, we argue that some of the key evolutionary novelties of adaptive immune systems of vertebrates are descendants of prokaryotic toxin–antitoxin (TA) systems. Cytidine deaminases and RAG recombinase have evolved from genotoxic enzymes to programmable editors of host genomes, supporting the astounding discriminatory capability of variable lymphocyte receptors of jawless vertebrates, as well as immunoglobulins and T cell receptors of jawed vertebrates. The evolutionarily recent lymphoid lineage is uniquely sensitive to mutations of the DNA maintenance methylase, which is an orphaned distant relative of prokaryotic restriction-modification systems. We discuss how the emergence of adaptive immunity gave rise to higher order genetic conflicts between genetic parasites and their vertebrate host. [ABSTRACT FROM AUTHOR]

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e. a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.] [ABSTRACT FROM AUTHOR]

Background: CRISPR-Cas-Docker is a web server for in silico docking experiments with CRISPR RNAs (crRNAs) and Cas proteins. This web server aims at providing experimentalists with the optimal crRNA-Cas pair predicted computationally when prokaryotic genomes have multiple CRISPR arrays and Cas systems, as frequently observed in metagenomic data. Results: CRISPR-Cas-Docker provides two methods to predict the optimal Cas protein given a particular crRNA sequence: a structure-based method (in silico docking) and a sequence-based method (machine learning classification). For the structure-based method, users can either provide experimentally determined 3D structures of these macromolecules or use an integrated pipeline to generate 3D-predicted structures for in silico docking experiments. Conclusion: CRISPR-Cas-Docker addresses the need of the CRISPR-Cas community to predict RNA–protein interactions in silico by optimizing multiple stages of computation and evaluation, specifically for CRISPR-Cas systems. CRISPR-Cas-Docker is available at www.crisprcasdocker.org as a web server, and at https://github.com/hshimlab/CRISPR-Cas-Docker as an open-source tool. [ABSTRACT FROM AUTHOR]

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life. [ABSTRACT FROM AUTHOR]

Simple Summary: Bees are one of the mega-diverse class of insects in Hymenoptera. As predators and the main crops pollinator, bees play an essential role in improving crop yield and the quality whiles providing substantial economic importance. Transposable elements (TEs) exist widely in nature and they make up a significant portion of prokaryotic and eukaryotic genomes. piggyBac (PB) transposons, belonging to DNA transposons, have been identified in plants, fungi, and animals. The advantages of the modified PB transposon system or high transposable efficiency and broad host ranges have led to their widespread application in transgenesis, mutagenesis, and gene therapy. At present, there is no systematic report on PB transposon in bees. Therefore, this study focused on annotating the transposons of the PB superfamily in bees to reveal the distribution, diversity, structural characteristics, evolutionary pattern and activity of PB transposons in the available genomes of bees. In this study, we investigated the presence of piggyBac (PB) transposons in 44 bee genomes from the Apoidea order, which is a superfamily within the Hymenoptera, which includes a large number of bee species crucial for pollination. We annotated the PB transposons in these 44 bee genomes and examined their evolution profiles, including structural characteristics, distribution, diversity, activity, and abundance. The mined PB transposons were divided into three clades, with uneven distribution in each genus of PB transposons in Apoidea. The complete PB transposons we discovered are around 2.23–3.52 kb in length and encode transposases of approximately 580 aa, with terminal inverted repeats (TIRs) of about 14 bp and 4 bp (TTAA) target-site duplications. Long TIRs (200 bp, 201 bp, and 493 bp) were also detected in some species of bees. The DDD domains of the three transposon types were more conserved, while the other protein domains were less conserved. Generally, most PB transposons showed low abundance in the genomes of Apoidea. Divergent evolution dynamics of PB were observed in the genomes of Apoidea. PB transposons in some identified species were relatively young, whiles others were older and with some either active or inactive. In addition, multiple invasions of PB were also detected in some genomes of Apoidea. Our findings highlight the contribution of PB transposons to genomic variation in these species and suggest their potential as candidates for future gene transfer tools. [ABSTRACT FROM AUTHOR]

Motivation Several genome annotation tools standardize annotation outputs for comparability. During standardization, these tools do not allow user-friendly customization of annotation databases; limiting their flexibility and applicability in downstream analysis. Results StandEnA is a user-friendly command-line tool for Linux that facilitates the generation of custom databases by retrieving protein sequences from multiple databases. Directed by a user-defined list of standard names, StandEnA retrieves synonyms to search for corresponding sequences in a set of public databases. Custom databases are used in prokaryotic genome annotation to generate standardized presence–absence matrices and reference files containing standard database identifiers. To showcase StandEnA, we applied it to six metagenome-assembled genomes to analyze three different pathways. Availability and implementation StandEnA is an open-source software available at https://github.com/mdsufz/StandEnA. Supplementary information Supplementary data are available at Bioinformatics Advances online. [ABSTRACT FROM AUTHOR]

Background: Wastewater treatment plants (WWTPs) are one of the largest biotechnology applications in the world and are of critical importance to modern urban societies. An accurate evaluation of the microbial dark matter (MDM, microorganisms whose genomes remain uncharacterized) proportions in WWTPs is of great value, while there is no such research yet. This study conducted a global meta-analysis of MDM in WWTPs with 317,542 prokaryotic genomes from the Genome Taxonomy Database and proposed a "wanted list" for priority targets in further investigations of activated sludge. Results: Compared with the Earth Microbiome Project data, WWTPs had relatively lower genome-sequenced proportions of prokaryotes than other ecosystems, such as the animal related environments. Analysis showed that the median proportions of the genome-sequenced cells and taxa (100% identity and 100% coverage in 16S rRNA gene region) in WWTPs reached 56.3% and 34.5% for activated sludge, 48.6% and 28.5% for aerobic biofilm, and 48.3% and 28.5% for anaerobic digestion sludge, respectively. This result meant MDM had high proportions in WWTPs. Besides, all of the samples were occupied by a few predominant taxa, and the majority of the sequenced genomes were from pure cultures. The global-scale "wanted list" for activated sludge contained four phyla that have few representatives and 71 operational taxonomic units with the majority of them having no genome or isolate yet. Finally, several genome mining methods were verified to successfully recover genomes from activated sludge such as hybrid assembly of the second- and third-generation sequencing. Conclusions: This work elucidated the proportion of MDM in WWTPs, defined the "wanted list" of activated sludge for future investigations, and certified potential genome recovery methods. The proposed methodology of this study can be applied to other ecosystems and improve understanding of ecosystem structure across diverse habitats. Ed3xU4yrM6Da7r9rFgJoVi Video Abstract [ABSTRACT FROM AUTHOR]

Motivation The rank distance model represents genome rearrangements in multi-chromosomal genomes as matrix operations, which allows the reconstruction of parsimonious histories of evolution by rearrangements. We seek to generalize this model by allowing for genomes with different gene content, to accommodate a broader range of biological contexts. We approach this generalization by using a matrix representation of genomes. This leads to simple distance formulas and sorting algorithms for genomes with different gene contents, but without duplications. Results We generalize the rank distance to genomes with different gene content in two different ways. The first approach adds insertions, deletions and the substitution of a single extremity to the basic operations. We show how to efficiently compute this distance. To avoid genomes with incomplete markers, our alternative distance, the rank-indel distance, only uses insertions and deletions of entire chromosomes. We construct phylogenetic trees with our distances and the DCJ-Indel distance for simulated data and real prokaryotic genomes, and compare them against reference trees. For simulated data, our distances outperform the DCJ-Indel distance using the Quartet metric as baseline. This suggests that rank distances are more robust for comparing distantly related species. For real prokaryotic genomes, all rearrangement-based distances yield phylogenetic trees that are topologically distant from the reference (65% similarity with Quartet metric), but are able to cluster related species within their respective clades and distinguish the Shigella strains as the farthest relative of the Escherichia coli strains, a feature not seen in the reference tree. Availability and implementation Code and instructions are available at https://github.com/meidanis-lab/rank-indel. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Background: Microsatellites are a ubiquitous occurrence in prokaryotic and eukaryotic genomes. Microsatellites have become one of the most popular classes of genetic markers due to their high reproducibility, multi-allelic nature, co-dominant mode of inheritance, abundance and wide genome coverage. We characterised microsatellites in the genomes and genes of two bat species, Pteropus vampyrus and Miniopterus natalensis. This characterisation was used for gene ontology analysis and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment of coding sequences (CDS). Results: Compared to M. natalensis, the genome size of P. vampyrus is larger and contains more microsatellites, but the total diversity of both species is similar. Mononucleotide and dinucleotide repeats were the most diverse in the genome of the two species. In each bat species, the microsatellite bias was obvious. The microsatellites with the largest number of repeat motifs in P. vampyrus from mononucleotide to hexanucleotide were (A)n, (AC)n, (CAA)n, (AAAC)n, (AACAA)n and (AAACAA)n, with frequencies of 97.94%, 58.75%, 30.53%, 22.82%, 54.68% and 22.87%, respectively, while in M. natalensis were (A)n, (AC)n, (TAT)n, (TTTA)n, (AACAA)n and (GAGAGG)n, with of 92.00%, 34.08%, 40.36%, 21.83%, 25.42% and 12.79%, respectively. In both species, the diversity of microsatellites was highest in intergenic regions, followed by intronic, untranslated and exonic regions and lowest in coding regions. Location analysis indicated that microsatellites were mainly concentrated at both ends of the genes. Microsatellites in the CDS are thus subject to higher selective pressure. In the GO analysis, two unique GO terms were found only in P. vampyrus and M. natalensis, respectively. In KEGG enriched pathway, the biosynthesis of other secondary metabolites and metabolism of other amino acids in metabolism pathways were present only in M. natalensis. The combined biological process, cellular components and molecular function ontology are reflected in the GO analysis and six functional enrichments in KEGG annotation, suggesting advantageous mutations during species evolution. Conclusions: Our study gives a comparative characterization of the genomes of microsatellites composition in the two bat species. And also allow further study on the effect of microsatellites on gene function as well as provide an insight into the molecular basis for species adaptation to new and changing environments. [ABSTRACT FROM AUTHOR]