Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont

Benchmarking of five NGS mapping tools for the reference alignment of bacterial outer membrane vesicles-associated small RNAs

Advances in small RNAs (sRNAs)-related studies have posed a challenge for NGS-related bioinformatics, especially regarding the correct mapping of sRNAs. Depending on the algorithms and scoring matrices on which they are based, aligners are influenced by the characteristics of the dataset and the reference genome. These influences have been studied mainly in eukaryotes and to some extent in prokaryotes. However, in bacteria, the selection of aligners depending on sRNA-seq data associated with outer membrane vesicles (OMVs) and the features of the corresponding bacterial reference genome has not yet been investigated. We selected five aligners: BBmap, Bowtie2, BWA, Minimap2 and Segemehl, known for their generally good performance, to test them in mapping OMV-associated sRNAs from Aliivibrio fischeri to the bacterial reference genome. Significant differences in the performance of the five aligners were observed, resulting in differential recognition of OMV-associated sRNA biotypes in A. fischeri. Our results suggest that aligner(s) should not be arbitrarily selected for this task, which is often done, as this can be detrimental to the biological interpretation of NGS analysis results. Since each aligner has specific advantages and disadvantages, these need to be considered depending on the characteristics of the input OMV sRNAs dataset and the corresponding bacterial reference genome to improve the detection of existing, biologically important OMV sRNAs. Until we learn more about these dependencies, we recommend using at least two, preferably three, aligners that have good metrics for the given dataset/bacterial reference genome. The overlapping results should be considered trustworthy, yet their differences should not be dismissed lightly, but treated carefully in order not to overlook any biologically important OMV sRNA. This can be achieved by applying the intersect-then-combine approach. For the mapping of OMV-associated sRNAs of A. fischeri to the reference genome organized into two circular chromosomes and one circular plasmid, containing copies of sequences with rRNA- and tRNA-related features and no copies of sequences with protein-encoding features, if the aligners are used with their default parameters, we advise avoiding Segemehl, and recommend using the intersect-then-combine approach with BBmap, BWA and Minimap2 to improve the potential for discovery of biologically important OMV-associated sRNAs.

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

Many purported pseudogenes in bacterial genomes are bona fide genes

BACKGROUND

Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly.

RESULTS

Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors.

CONCLUSIONS

Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

A large-scale phylogeny-guided analysis of pseudogenes in Pseudomonas aeruginosa bacterium

Pseudogenes, once considered "junk DNA" based on the incorrect assumption that the absence of full coding potential means a complete lack of functionality, have recently become a subject of significant interest in the scientific community. Concurrently, it is widely assumed that bacterial genomes are compact and have a high density of coding genes with little room for non-coding genes, including pseudogenes. A key aspect of genome annotation is the correct identification of genes and the distinction between coding genes and pseudogenes, as it directly impacts functional and comparative genomics studies. In this study, we analyzed the genomic data of 4,699 strains of the bacterium Pseudomonas aeruginosa (P. aeruginosa) as they exhibit high variability in the number of annotated pseudogenes. In particular, we looked for correlations between the number of pseudogenes and other genomic and meta-features of the strains. We identified clusters of orthologous genes and pseudogenes and compared cluster size distributions and length homogeneity within clusters. We then mapped and examined orthology relationships between genes and pseudogenes. Additionally, we generated a phylogenetic tree of the strains and found that phylogenetically related strains are more homogeneous in the number of pseudogenes and share a significant amount of pseudogenes. Finally, we delved into clusters of orthologous genes and pseudogenes and quantified their phylogenetic neighborhood, classifying pseudogenes into evolutionary preserved pseudogenes, mis-annotated pseudogenes, or pseudogenes formed by failed horizontal transfer events. This in-depth study provides important insights that can be incorporated into pseudogene annotation pipelines in the future. IMPORTANCE Accurate annotation of genes and pseudogenes is vital for comparative genomics analysis. Recent studies have shown that bacterial pseudogenes have an important role in regulatory processes and can provide insight into the evolutionary history of homologous genes or the genome as a whole. Due to pseudogenes' nature as non-functional genes, there is no commonly accepted definition of a pseudogene, which poses difficulties in verifying the annotation through experimental methods and resolving discrepancies among different annotation techniques. Our study introduces an in-depth analysis of annotated genes and pseudogenes and insights that can be incorporated into improved pseudogene annotation pipelines in the future.

Bee foraging preferences, microbiota and pathogens revealed by direct shotgun metagenomics of honey

Honeybees (Apis mellifera) continue to succumb to human and environmental pressures despite their crucial role in providing essential ecosystem services. Owing to their foraging and honey production activities, honeybees form complex relationships with species across all domains, such as plants, viruses, bacteria and other hive pests, making honey a valuable biomonitoring tool for assessing their ecological niche. Thus, the application of honey shotgun metagenomics (SM) has paved the way for a detailed description of the species honeybees interact with. Nevertheless, SM bioinformatics tools and DNA extraction methods rely on resources not necessarily optimized for honey. In this study, we compared five widely used taxonomic classifiers using simulated species communities commonly found in honey. We found that Kraken 2 with a threshold of 0.5 performs best in assessing species distribution. We also optimized a simple NaOH-based honey DNA extraction methodology (Direct-SM), which profiled species seasonal variability similarly to an established column-based DNA extraction approach (SM). Both approaches produce results consistent with melissopalinology analysis describing the botanical landscape surrounding the apiary. Interestingly, we detected a strong stability of the bacteria constituting the core and noncore gut microbiome across seasons, pointing to the potential utility of honey for noninvasive assessment of bee microbiota. Finally, the Direct-SM approach to detect Varroa correlates well with the biomonitoring of mite infestation observed in hives. These observations suggest that Direct-SM methodology has the potential to comprehensively describe honeybee ecological niches and can be tested as a building block for large-scale studies to assess bee health in the field.

Phenogenomic Characterization of a Newly Domesticated and Novel Species from the Genus Verrucosispora

The concept of bacterial dark matter stems from our inability to culture most microbes and represents a fundamental gap in our knowledge of microbial diversity. Here, we present the domestication of such an organism: a previously uncultured, novel species from the rare Actinomycetes genus Verrucosispora. Although initial recovery took >4 months, isolation of phenotypically distinct, domesticated generations occurred within weeks. Two isolates were subjected to phenogenomic analyses, revealing domestication correlated with enhanced growth rates in nutrient-rich media but diminished capacity to metabolize diverse amino acids. This is seemingly mediated by genomic atrophy through a mixed approach of pseudogenization and reversion of pseudogenization of amino acid metabolism genes. Conversely, later generational strains had enhanced spore germination rates, potentially through the reversion of a sporulation-associated kinase from pseudogene to true gene status. We observed that our most wild-type isolate had the greatest potential for antibacterial activity, which correlated with extensive mutational attrition of biosynthetic gene clusters in domesticated strains. Comparative analyses revealed wholesale genomic reordering in strains, with widespread single nucleotide polymorphism, indel, and pseudogene-impactful mutations observed. We hypothesize that domestication of this previously unculturable organism resulted from the shedding of genomic flexibility required for life in a dynamic marine environment, parsing out genetic redundancy to allow for a newfound cultivable amenability. IMPORTANCE The majority of environmental bacteria cannot be cultured within the laboratory. Understanding why only certain environmental isolates can be recovered is key to unlocking the abundant microbial dark matter that is widespread on our planet. In this study, we present not only the culturing but domestication of just such an organism. Although initial recovery took >4 months, we were able to isolate distinct, subpassaged offspring from the originating colony within mere weeks. A phenotypic and genotypic analysis of our generational strains revealed that adaptation to life in the lab occurred as a result of wholesale mutational changes. These permitted an enhanced ability for growth in nutrient rich media but came at the expense of reduced genomic flexibility. We suggest that without dynamic natural environmental stressors our domesticated strains effectively underwent genomic atrophy as they adapted to static conditions experienced in the laboratory.

Ecological Divergence Within the Enterobacterial Genus Sodalis: From Insect Symbionts to Inhabitants of Decomposing Deadwood

The bacterial genus Sodalis is represented by insect endosymbionts as well as free-living species. While the former have been studied frequently, the distribution of the latter is not yet clear. Here, we present a description of a free-living strain, Sodalis ligni sp. nov., originating from decomposing deadwood. The favored occurrence of S. ligni in deadwood is confirmed by both 16S rRNA gene distribution and metagenome data. Pangenome analysis of available Sodalis genomes shows at least three groups within the Sodalis genus: deadwood-associated strains, tsetse fly endosymbionts and endosymbionts of other insects. This differentiation is consistent in terms of the gene frequency level, genome similarity and carbohydrate-active enzyme composition of the genomes. Deadwood-associated strains contain genes for active decomposition of biopolymers of plant and fungal origin and can utilize more diverse carbon sources than their symbiotic relatives. Deadwood-associated strains, but not other Sodalis strains, have the genetic potential to fix N2, and the corresponding genes are expressed in deadwood. Nitrogenase genes are located within the genomes of Sodalis, including S. ligni, at multiple loci represented by more gene variants. We show decomposing wood to be a previously undescribed habitat of the genus Sodalis that appears to show striking ecological divergence.

Incipient genome erosion and metabolic streamlining for antibiotic production in a defensive symbiont

Genome reduction is commonly observed in bacteria of several phyla engaging in obligate nutritional symbioses with insects. In Actinobacteria, however, little is known about the process of genome evolution, despite their importance as prolific producers of antibiotics and their increasingly recognized role as defensive partners of insects and other organisms. Here, we show that “Streptomyces philanthi,” a defensive symbiont of digger wasps, has a G+C-enriched genome in the early stages of erosion, with inactivating mutations in a large proportion of genes, causing dependency on its hosts for certain nutrients, which was validated in axenic symbiont cultures. Additionally, overexpressed catabolic and biosynthetic pathways of the bacteria inside the host indicate host–symbiont metabolic integration for streamlining and control of antibiotic production.

Genome erosion is a frequently observed result of relaxed selection in insect nutritional symbionts, but it has rarely been studied in defensive mutualisms. Solitary beewolf wasps harbor an actinobacterial symbiont of the genus Streptomyces that provides protection to the developing offspring against pathogenic microorganisms. Here, we characterized the genomic architecture and functional gene content of this culturable symbiont using genomics, transcriptomics, and proteomics in combination with in vitro assays. Despite retaining a large linear chromosome (7.3 Mb), the wasp symbiont accumulated frameshift mutations in more than a third of its protein-coding genes, indicative of incipient genome erosion. Although many of the frameshifted genes were still expressed, the encoded proteins were not detected, indicating post-transcriptional regulation. Most pseudogenization events affected accessory genes, regulators, and transporters, but “Streptomyces philanthi” also experienced mutations in central metabolic pathways, resulting in auxotrophies for biotin, proline, and arginine that were confirmed experimentally in axenic culture. In contrast to the strong A+T bias in the genomes of most obligate symbionts, we observed a significant G+C enrichment in regions likely experiencing reduced selection. Differential expression analyses revealed that—compared to in vitro symbiont cultures—“S. philanthi” in beewolf antennae showed overexpression of genes for antibiotic biosynthesis, the uptake of host-provided nutrients and the metabolism of building blocks required for antibiotic production. Our results show unusual traits in the early stage of genome erosion in a defensive symbiont and suggest tight integration of host–symbiont metabolic pathways that effectively grants the host control over the antimicrobial activity of its bacterial partner.

Simulating the evolutionary trajectories of metabolic pathways for insect symbionts in the genus Sodalis

Insect-bacterial symbioses are ubiquitous, but there is still much to uncover about how these relationships establish, persist and evolve. The tsetse endosymbiont Sodalis glossinidius displays intriguing metabolic adaptations to its microenvironment, but the process by which this relationship evolved remains to be elucidated. The recent chance discovery of the free-living species of the genus Sodalis, Sodalis praecaptivus, provides a serendipitous starting point from which to investigate the evolution of this symbiosis. Here, we present a flux balance model for S. praecaptivus and empirically verify its predictions. Metabolic modelling is used in combination with a multi-objective evolutionary algorithm to explore the trajectories that S. glossinidius may have undertaken from this starting point after becoming internalized. The order in which key genes are lost is shown to influence the evolved populations, providing possible targets for future in vitro genetic manipulation. This method provides a detailed perspective on possible evolutionary trajectories for S. glossinidius in this fundamental process of evolutionary and ecological change.

Dynamics of Insect-Microbiome Interaction Influence Host and Microbial Symbiont

Insects share an intimate relationship with their gut microflora and this symbiotic association has developed into an essential evolutionary outcome intended for their survival through extreme environmental conditions. While it has been clearly established that insects, with very few exceptions, associate with several microbes during their life cycle, information regarding several aspects of these associations is yet to be fully unraveled. Acquisition of bacteria by insects marks the onset of microbial symbiosis, which is followed by the adaptation of these bacterial species to the gut environment for prolonged sustenance and successful transmission across generations. Although several insect-microbiome associations have been reported and each with their distinctive features, diversifications and specializations, it is still unclear as to what led to these diversifications. Recent studies have indicated the involvement of various evolutionary processes operating within an insect body that govern the transition of a free-living microbe to an obligate or facultative symbiont and eventually leading to the establishment and diversification of these symbiotic relationships. Data from various studies, summarized in this review, indicate that the symbiotic partners, i.e., the bacteria and the insect undergo several genetic, biochemical and physiological changes that have profound influence on their life cycle and biology. An interesting outcome of the insect-microbe interaction is the compliance of the microbial partner to its eventual genome reduction. Endosymbionts possess a smaller genome as compared to their free-living forms, and thus raising the question what is leading to reductive evolution in the microbial partner. This review attempts to highlight the fate of microbes within an insect body and its implications for both the bacteria and its insect host. While discussion on each specific association would be too voluminous and outside the scope of this review, we present an overview of some recent studies that contribute to a better understanding of the evolutionary trajectory and dynamics of the insect-microbe association and speculate that, in the future, a better understanding of the nature of this interaction could pave the path to a sustainable and environmentally safe way for controlling economically important pests of crop plants.