Pseudofinder: detection of pseudogenes in prokaryotic genomes

Pseudogenes and host specialization in the emergent bacterial plant pathogen Xylella fastidiosa

Pseudogenes are regarded as "junk" DNA, representing vestigial functions no longer needed for fitness. Accordingly, a higher number of pseudogenes in a bacterial human pathogen was previously hypothesized to be a hallmark of host specialists. In this study, we tested this hypothesis on the emergent bacterial plant pathogen Xylella fastidiosa (Xf) to link pseudogene makeup and host specificity. Xf is an ideal subject for these studies by being a xylem-limited pathogen that underwent extensive genome reduction. Using natural host range data of 151 strains and Pseudofinder analysis on Xf whole genome sequences, we observed that Xf subsp. sandyi had the highest pseudogene content, followed by subsp. morus, while subsp. pauca, fastidiosa, and multiplex had the lowest. The first two subspecies are known to have a limited host range compared to the others, aligning with the hypothesis of a greater number of pseudogenes corresponding to narrower host range. Weed isolates are presumably host specialists because they had the highest pseudogene content. Using a thorough pseudogene map across genomes and empirical pathogenicity data on blueberries, we screened for genes potentially involved in blueberry specialization. Targets were identified by selecting sequences pseudogenized (i) in strains infecting hosts different from blueberry and (ii) only in blueberry asymptomatic strains. Six sequences were identified with a potential role in blueberry infection, including one that was common between the two criteria. Here, we generated hypotheses on host range and specificity of Xf strains that need to be tested experimentally to help understand this devastating plant pathogen.IMPORTANCEXylella fastidiosa is a highly destructive plant pathogen that infects hundreds of landscape and agriculturally important plant species mainly in Europe and the Americas. Nevertheless, the host range of specific genotypes and underlying mechanisms of host specificity remain unclear. These are important aspects to determine the potential risk of infection in specific areas depending on the genetic makeup of the pathogen population and hosts present. This study offers valuable insights into the role of pseudogenization in the genomes of different X. fastidiosa strains, linking it to host specialization. Despite the limited information available for the host range of different strains of this pathogen, this research proposes a relationship between the abundance of pseudogenes and host specificity. These findings are essential for predicting potential host shifts by this pathogen, aiding in the development of strategies to prevent its spread. Additionally, the identification of candidate genes putatively important for symptom development in blueberries offers targets for prevention and control efforts.

A genomic and phenotypic investigation of pigeon-adaptive Salmonella

Salmonella, a significant threat to public safety, inflicts substantial economic losses on the poultry industry. The unique "parental feeding" breeding model of pigeon farms, against the "all-in & all-out" biosecurity strategy, makes them susceptible to Salmonella infections and subsequent outbreaks of pigeon paratyphoid. This study initially studied three pigeon paratyphoid outbreak incidents in Henan, China, in which 53 strains of pigeon-origin Salmonella Typhimurium (STM) were identified. Whole-genome sequencing (WGS) and antimicrobial-resistant profile analysis revealed that the three outbreaks were caused by distinct STM clones (ST128-DT2, ST19-DT99). Global phylogenetic analysis suggested that the United States is a possible origin, indicating a risk of intercontinental transmission via pigeon eggs. Further bacterial virulence and invasion assays, including in vitro and in vivo assays, revealed that pigeon-host-adaptive STM, compared to broad-host-range STM, carried fewer resistance genes, exhibited higher invasion indices and pseudogene levels, displayed a non-rdar (red dry and rough) phenotype, and had strong biofilm formation capability. Additionally, they showed reduced virulence and invasiveness in mice but a pigeon-adaptive feature in cogent models. The collective results support the host adaptation for pigeons among DT2 and DT99 phage-type isolates.

The nitrogen-fixing fern Azolla has a complex microbiome characterized by varying degrees of cophylogenetic signal

PREMISE

Azolla is a genus of floating ferns that has closely evolved with a vertically transmitted obligate cyanobacterium endosymbiont-Anabaena azollae-that fixes nitrogen. There are also other lesser-known Azolla symbionts whose role and mode of transmission are unknown.

METHODS

We sequenced 112 Azolla specimens collected across the state of California and characterized their metagenomes to identify the common bacterial endosymbionts and assess their patterns of interaction.

RESULTS

Four genera were found across all samples, establishing that multiple Azolla endosymbionts were consistently present. We found varying degrees of cophylogenetic signal across these taxa as well as varying degrees of isolation by distance and of pseudogenation, which demonstrates that multiple processes underlie how this endosymbiotic community is constituted. We also characterized the entire Azolla leaf pocket microbiome.

CONCLUSIONS

These results show that the Azolla symbiotic community is complex and features members at potentially different stages of symbiosis evolution, further supporting the utility of the Azolla microcosm as a system for studying the evolution of symbioses.

From eggs to guts: Symbiotic association of Sodalis nezarae sp. nov. with the Southern green shield bug Nezara viridula

Phytophagous insects engage in symbiotic relationships with bacteria that contribute to digestion, nutrient supplementation, and development of the host. The analysis of shield bug microbiomes has been mainly focused on the gut intestinal tract predominantly colonized by Pantoea symbionts and other microbial community members in the gut or other organs have hardly been investigated. In this study, we reveal that the Southern green shield bug Nezara viridula harbours a Sodalis symbiont in several organs, with a notable prevalence in salivary glands, and anterior regions of the midgut. Removing external egg microbiota via sterilization profoundly impacted insect viability but did not disrupt the vertical transmission of Sodalis and Pantoea symbionts. Based on the dominance of Sodalis in testes, we deduce that N. viridula males could be involved in symbiont vertical transmission. Genomic analyses comparing Sodalis species revealed that Sodalis sp. Nvir shares characteristics with both free-living and obligate insect-associated Sodalis spp. Sodalis sp. Nvir also displays genome instability typical of endosymbiont lineages, which suggests ongoing speciation to an obligate endosymbiont. Together, our study reveals that shield bugs harbour unrecognized symbionts that might be paternally transmitted.

mettannotator: a comprehensive and scalable Nextflow annotation pipeline for prokaryotic assemblies

SUMMARY

In recent years, there has been a surge in prokaryotic genome assemblies, coming from both isolated organisms and environmental samples. These assemblies often include novel species that are poorly represented in reference databases creating a need for a tool that can annotate both well-described and novel taxa, and can run at scale. Here, we present mettannotator-a comprehensive, scalable Nextflow pipeline for prokaryotic genome annotation that identifies coding and noncoding regions, predicts protein functions, including antimicrobial resistance, and delineates gene clusters. The pipeline summarizes these results in a GFF (General Feature Format) file that can be easily utilized in downstream analysis or visualized using common genome browsers. Here, we show how it works on 200 genomes from 29 prokaryotic phyla, including isolate genomes and known and novel metagenome-assembled genomes, and present metrics on its performance in comparison to other tools.

AVAILABILITY AND IMPLEMENTATION

The pipeline is written in Nextflow and Python and published under an open source Apache 2.0 licence. Instructions and source code can be accessed at https://github.com/EBI-Metagenomics/mettannotator. The pipeline is also available on WorkflowHub: https://workflowhub.eu/workflows/1069.

Genome evolution following an ecological shift in nectar-dwelling Acinetobacter

The bacterial genus Acinetobacter includes species found in environmental habitats like soil and water, as well as taxa adapted to be host-associated or pathogenic. High genetic diversity may allow for this habitat flexibility, but the specific genes underlying switches between habitats are poorly understood. One lineage of Acinetobacter has undergone a substantial habitat change by evolving from a presumed soil-dwelling ancestral state to thrive in floral nectar. Here, we compared the genomes of floral-dwelling and pollinator-associated Acinetobacter, including newly described species, with genomes from relatives found in other environments to determine the genomic changes associated with this ecological shift. Following one evolutionary origin of floral nectar adaptation, nectar-dwelling Acinetobacter taxa have undergone reduction in genome size compared with relatives and have experienced dynamic gene gains and losses as they diversified. Gene content changes suggest a shift to metabolism of monosaccharides rather than diverse carbohydrates, and scavenging of nitrogen sources, which we predict to be beneficial in nectar environments. Gene gains appear to result from duplication events, evolutionary divergence, and horizontal gene transfer. Most notably, nectar-dwelling Acinetobacter acquired the ability to degrade pectin from plant pathogens, and the genes underlying this ability have duplicated and are under selection within the clade. We hypothesize that this ability was a key trait for adaptation to floral nectar, as it could improve access to nutrients in the nutritionally unbalanced habitat of nectar. These results identify the genomic changes and traits coinciding with a dramatic habitat switch from soil to floral nectar.

IMPORTANCE

Many bacteria, including the genus Acinetobacter, commonly evolve to exploit new habitats. However, the genetic changes that underlie habitat switches are often unknown. Floral nectar is home to specialized microbes that can grow in this nutritionally unbalanced habitat. Several specialized Acinetobacter species evolved from soil-dwelling relatives to become common and abundant in floral nectar. Here, we investigate the genomic adaptations required to successfully colonize a novel habitat like floral nectar. We performed comparative genomics analyses between nectar-dwelling Acinetobacter and Acinetobacter species from other environments, like soil and water. We find that although gene loss coincided with the switch to living in nectar, gains of specific genes from other bacteria may have been particularly important for this ecological change. Acinetobacter living in nectar gained genes for degrading pectin, a plant polysaccharide, which may improve access to nutrients in their environment. These findings shed light on how evolutionary novelty evolves in bacteria.

Somy evolution in the honey bee infecting trypanosomatid parasite Lotmaria passim

Lotmaria passim is a ubiquitous trypanosomatid parasite of honey bees nestled within the medically important subfamily Leishmaniinae. Although this parasite is associated with honey bee colony losses, the original draft genome-which was completed before its differentiation from the closely related Crithidia mellificae-has remained the reference for this species despite lacking improvements from newer methodologies. Here, we report the updated sequencing, assembly, and annotation of the BRL-type (Bee Research Laboratory) strain (ATCC PRA-422) of Lotmaria passim. The nuclear genome assembly has been resolved into 31 complete chromosomes and is paired with an assembled kinetoplast genome consisting of a maxicircle and 30 minicircle sequences. The assembly spans 33.7 Mb and contains very little repetitive content, from which our annotation of both the nuclear assembly and kinetoplast predicted 10,288 protein-coding genes. Analyses of the assembly revealed evidence of a recent chromosomal duplication event within chromosomes 5 and 6 and provided evidence for a high level of aneuploidy in this species, mirroring the genomic flexibility employed by other trypanosomatids as a means of adaptation to different environments. This high-quality reference can therefore provide insights into adaptations of trypanosomatids to the thermally regulated, acidic, and phytochemically rich honey bee hindgut niche, which offers parallels to the challenges faced by other Leishmaniinae during the challenges they undergo within insect vectors, during infection of mammals, and exposure to antiparasitic drugs throughout their multi-host life cycles. This reference will also facilitate investigations of strain-specific genomic polymorphisms, their role in pathogenicity, and the development of treatments for pollinator infection.

Shrinking in the dark: Parallel endosymbiont genome erosions are associated with repeated host transitions to an underground life

Microbial symbioses have had profound impacts on the evolution of animals. Conversely, changes in host biology may impact the evolutionary trajectory of symbionts themselves. Blattabacterium cuenoti is present in almost all cockroach species and enables hosts to subsist on a nutrient-poor diet. To investigate if host biology has impacted Blattabacterium at the genomic level, we sequenced and analyzed 25 genomes from Australian soil-burrowing cockroaches (Blaberidae: Panesthiinae), which have undergone at least seven separate subterranean, subsocial transitions from above-ground, wood-feeding ancestors. We find at least three independent instances of genome erosion have occurred in Blattabacterium strains exclusive to Australian soil-burrowing cockroaches. These shrinkages have involved the repeated inactivation of genes involved in amino acid biosynthesis and nitrogen recycling, the core role of Blattabacterium in the host-symbiont relationship. The most drastic of these erosions have occurred in hosts thought to have transitioned underground the earliest relative to other lineages, further suggestive of a link between gene loss in Blattabacterium and the burrowing behavior of hosts. As Blattabacterium is unable to fulfill its core function in certain host lineages, these findings suggest soil-burrowing cockroaches must acquire these nutrients from novel sources. Our study represents one of the first cases, to our knowledge, of parallel host adaptations leading to concomitant parallelism in their mutualistic symbionts, further underscoring the intimate relationship between these two partners.

A unique symbiosome in an anaerobic single-celled eukaryote

Symbiotic relationships between eukaryotes and prokaryotes played pivotal roles in the evolution of life and drove the emergence of specialized symbiotic structures in animals, plants and fungi. The host-evolved symbiotic structures of microbial eukaryotes - the vast majority of such hosts in nature - remain largely unstudied. Here we describe highly structured symbiosomes within three free-living anaerobic protists (Anaeramoeba spp.). We dissect this symbiosis using complete genome sequencing and transcriptomics of host and symbiont cells coupled with fluorescence in situ hybridization, and 3D reconstruction using focused-ion-beam scanning electron microscopy. The emergence of the symbiosome is underpinned by expansion of gene families encoding regulators of membrane trafficking and phagosomal maturation and extensive bacteria-to-eukaryote lateral transfer. The symbionts reside deep within a symbiosomal membrane network that enables metabolic syntrophy by precisely positioning sulfate-reducing bacteria alongside host hydrogenosomes. Importantly, the symbionts maintain connections to the Anaeramoeba plasma membrane, blurring traditional boundaries between ecto- and endosymbiosis.

Genetic potential for aerobic respiration and denitrification in globally distributed respiratory endosymbionts

The endosymbiont Candidatus Azoamicus ciliaticola was proposed to generate ATP for its eukaryotic host, an anaerobic ciliate of the Plagiopylea class, fulfilling a function analogous to mitochondria in other eukaryotic cells. The discovery of this respiratory endosymbiosis has major implications for both evolutionary history and ecology of microbial eukaryotes. However, with only a single species described, knowledge of its environmental distribution and diversity is limited. Here we report four complete, circular metagenome assembled genomes (cMAGs) representing respiratory endosymbionts inhabiting groundwater in California, Ohio, and Germany. These cMAGs form two lineages comprising a monophyletic clade within the uncharacterized gammaproteobacterial order UBA6186, enabling evolutionary analysis of their key protein complexes. Strikingly, all four cMAGs encode a cytochrome cbb3 oxidase, which indicates that these endosymbionts have the capacity for aerobic respiration. Accordingly, we detect these respiratory endosymbionts in diverse habitats worldwide, thus further expanding the ecological scope of this respiratory symbiosis.

Revisiting the functional annotation of TriTryp using sequence similarity tools

Trypanosomatids are the causative agents of deadly diseases in humans and livestock. Given the high phylogenetic distance of trypanosomatids from model organisms, these organisms have ample unannotated genes. Manual functional annotation is time-consuming, highlighting the importance of automated functional annotation tools. The development of automated functional tools is a hot research topic, and multiple tools have been developed for the task. PANNZER2 is an automated functional annotation tool that merely relies on the sequence similarity of the query to the annotated proteins. We tried PANNZER2 on Trypanosoma brucei, the most studied organism among trypanosomatids, to see if it could improve our knowledge of the functions of the genes. Even with the availability of automated annotation tools like InterPro2GO in databases such as TriTrypDB, PANNZER2 has made surprisingly confident predictions for some hypothetical proteins in T. brucei. In this study, we identify gaps in such annotations because of not employing pairwise sequence alignment tools in TriTrypDB's automated annotation process. Our findings demonstrate that even the use of stringent cutoffs can successfully annotate a significant number of proteins. Additionally, we discovered that adjusting the open reading frames in certain genes leads to sequences with increased sequence signature coverage-characterized by the length covered by at least one sequence signature-compared to the original sequences. This enhanced sequence signature coverage suggests these genomic fragments could be pseudogenes. To facilitate further exploration, we developed a script to help identify potential pseudogenes within an organism's genome, offering researchers a new tool for genomic analysis and understanding. We extended all our analysis to Trypanosoma cruzi and Leishmania major to assess the impact of this approach across different species. Our study demonstrates that by utilizing pairwise sequence similarity alignment, even with stringent cutoffs, we can attribute 2986, 3953, and 3798 new GO terms to the genomes of T. brucei, T. cruzi, and L. major. Additionally, we found that 210, 239, and 29 genes exhibit increased sequence signature coverage following frame correction, suggesting the presence of pseudogenes.

Discovery of a novel Wolbachia in Heterodera expands nematode host distribution

Bioinformatics sequence data mining can reveal hidden microbial symbionts that might normally be filtered and removed as contaminants. Data mining can be helpful to detect Wolbachia, a widespread bacterial endosymbiont in insects and filarial nematodes whose distribution in plant-parasitic nematodes (PPNs) remains underexplored. To date, Wolbachia has only been reported a few PPNs, yet nematode-infecting Wolbachia may have been widespread in the evolutionary history of the phylum based on evidence of horizontal gene transfers, suggesting there may be undiscovered Wolbachia infections in PPNs. The goal of this study was to more broadly sample PPN Wolbachia strains in tylenchid nematodes to enable further comparative genomic analyses that may reveal Wolbachia's role and identify targets for biocontrol. Published whole-genome shotgun assemblies and their raw sequence data from 33 Meloidogyne spp. assemblies, seven Globodera spp. assemblies, and seven Heterodera spp. assemblies were analyzed to look for Wolbachia. No Wolbachia was found in Meloidogyne spp. and Globodera spp., but among seven genome assemblies for Heterodera spp., an H. schachtii assembly from the Netherlands was found to have a large Wolbachia-like sequence that, when re-assembled from reads, formed a complete, circular genome. Detailed analyses comparing read coverage, GC content, pseudogenes, and phylogenomic patterns clearly demonstrated that the H. schachtii Wolbachia represented a novel strain (hereafter, denoted wHet). Phylogenomic tree construction with PhyloBayes showed wHet was most closely related to another PPN Wolbachia, wTex, while 16S rRNA gene analysis showed it clustered with other Heterodera Wolbachia assembled from sequence databases. Pseudogenes in wHet suggested relatedness to the PPN clade, as did the lack of significantly enriched GO terms compared to PPN Wolbachia strains. It remains unclear whether the lack of Wolbachia in other published H. schachtii isolates represents the true absence of the endosymbiont from some hosts.

Rapid Loss of Nutritional Symbionts in an Endemic Hawaiian Herbivore Radiation Is Associated with Plant Galling Habit

Insect herbivores frequently cospeciate with symbionts that enable them to survive on nutritionally unbalanced diets. While ancient symbiont gain and loss events have been pivotal for insect diversification and feeding niche specialization, evidence of recent events is scarce. We examine the recent loss of nutritional symbionts (in as little as 1 MY) in sap-feeding Pariaconus, an endemic Hawaiian insect genus that has undergone adaptive radiation, evolving various galling and free-living ecologies on a single host-plant species, Metrosideros polymorpha within the last ∼5 MY. Using 16S rRNA sequencing, we investigated the bacterial microbiomes of 19 Pariaconus species and identified distinct symbiont profiles associated with specific host-plant ecologies. Phylogenetic analyses and metagenomic reconstructions revealed significant differences in microbial diversity and functions among psyllids with different host-plant ecologies. Within a few millions of years, Pariaconus species convergently evolved the closed-gall habit twice. This shift to enclosed galls coincided with the loss of the Morganella-like symbiont that provides the essential amino acid arginine to free-living and open-gall sister species. After the Pariaconus lineage left Kauai and colonized younger islands, both open- and closed-gall species lost the Dickeya-like symbiont. This symbiont is crucial for synthesizing essential amino acids (phenylalanine, tyrosine, and lysine) as well as B vitamins in free-living species. The recurrent loss of these symbionts in galling species reinforces evidence that galls are nutrient sinks and, combined with the rapidity of the evolutionary timeline, highlights the dynamic role of insect-symbiont relationships during the diversification of feeding ecologies. We propose new Candidatus names for the novel Morganella-like and Dickeya-like symbionts.

Franco D, Szklarzewicz T, Stroiński A, Łukasik P. Facultatively intrabacterial localization of a planthopper endosymbiont as an adaptation to its vertical transmission

Transovarial transmission is the most reliable way of passing on essential nutrient-providing endosymbionts from mothers to offspring. However, not all endosymbiotic microbes follow the complex path through the female host tissues to oocytes on their own. Here, we demonstrate an unusual transmission strategy adopted by one of the endosymbionts of the planthopper Trypetimorpha occidentalis (Hemiptera: Tropiduchidae) from Bulgaria. In this species, an Acetobacteraceae endosymbiont is transmitted transovarially within deep invaginations of cellular membranes of an ancient endosymbiont Sulcia-strikingly resembling recently described plant virus transmission. However, in males, Acetobacteraceae colonizes the same bacteriocytes as Sulcia but remains unenveloped. Then, the unusual endobacterial localization of Acetobacteraceae observed in females appears to be a unique adaptation to maternal transmission. Further, the symbiont's genomic features, including encoding essential amino acid biosynthetic pathways and its similarity to a recently described psyllid symbiont, suggest a unique combination of the ability to horizontally transmit among species and confer nutritional benefits. The close association with Acetobacteraceae symbiont correlates with the so-far-unreported level of genomic erosion of ancient nutritional symbionts of this planthopper. In Sulcia, this is reflected in substantial changes in genomic organization, reported for the first time in the symbiont renowned for its genomic stability. In Vidania, substantial gene loss resulted in one of the smallest genomes known, at 108.6 kb. Thus, the symbionts of T. occidentalis display a combination of unusual adaptations and genomic features that expand our understanding of how insect-microbe symbioses may transmit and evolve.IMPORTANCEReliable transmission across host generations is a major challenge for bacteria that associate with insects, and independently established symbionts have addressed this challenge in different ways. The facultatively endobacterial localization of Acetobacteraceae symbiont, enveloped by cells of ancient nutritional endosymbiont Sulcia in females but not males of the planthopper Trypetimorpha occidentalis, appears to be a unique adaptation to maternal transmission. Acetobacteraceae's genomic features indicate its unusual evolutionary history, and the genomic erosion experienced by ancient nutritional symbionts demonstrates the apparent consequences of such close association. Combined, this multi-partite symbiosis expands our understanding of the diversity of strategies that insect symbioses form and some of their evolutionary consequences.

Draft genome of the marine bacterium Alteromonas gracilis strain J4 isolated from the green coenocytic alga Caulerpa prolifera

Here, we present the draft genome sequence of Alteromonas gracilis strain J4, isolated from the green macroalga Caulerpa prolifera. The draft genome is 4,492,914 bp in size and contains 4,719 coding DNA sequences, 67 tRNAs, and 16 rRNA-coding genes. Strain J4 may exhibit host growth-promoting properties.

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BACKGROUND

To better understand the influence of habitat on the genetic content of bacteria, with a focus on members of Candidate Phyla Radiation (CPR) bacteria, we studied the effects of transitioning from soil via seepage waters to groundwater on genomic composition of ultra-small Parcubacteria, the dominating CPR class in seepage waters, using genome resolved metagenomics.

RESULTS

Bacterial metagenome-assembled genomes (MAGs), (318 total, 32 of Parcubacteria) were generated from seepage waters and compared directly to groundwater counterparts. The estimated average genome sizes of members of major phyla Proteobacteria, Bacteroidota and Cand. Patescibacteria (Candidate Phyla Radiation - CPR bacteria) were significantly higher in soil-seepage water as compared to their groundwater counterparts. Seepage water Parcubacteria (Paceibacteria) exhibited 1.18-fold greater mean genome size and 2-fold lower mean proportion of pseudogenes than those in groundwater. Bacteroidota and Proteobacteria also showed a similar trend of reduced genomes in groundwater compared to seepage. While exploring gene loss and adaptive gains in closely related CPR lineages in groundwater, we identified a membrane protein, and a lipoglycopeptide resistance gene unique to a seepage Parcubacterium genome. A nitrite reductase gene was also identified and was unique to the groundwater Parcubacteria genomes, likely acquired from other planktonic microbes via horizontal gene transfer.

CONCLUSIONS

Overall, our data suggest that bacteria in seepage waters, including ultra-small Parcubacteria, have significantly larger genomes and higher metabolic enrichment than their groundwater counterparts, highlighting possible genome streamlining of the latter in response to habitat selection in an oligotrophic environment.

In vitro, in planta, and comparative genomic analyses of Pseudomonas syringae pv. syringae strains of pepper (Capsicum annuum var. annuum)

Pseudomonas syringae pv. syringae (Pss) is an emerging phytopathogen that causes Pseudomonas leaf spot (PLS) disease in pepper plants. Pss can cause serious economic damage to pepper production, yet very little is known about the virulence factors carried by Pss that cause disease in pepper seedlings. In this study, Pss strains isolated from pepper plants showing PLS symptoms in Ohio between 2013 and 2021 (n = 16) showed varying degrees of virulence (Pss populations and disease symptoms on leaves) on 6-week-old pepper seedlings. In vitro studies assessing growth in nutrient-limited conditions, biofilm production, and motility also showed varying degrees of virulence, but in vitro and in planta variation in virulence between Pss strains did not correlate. Comparative whole-genome sequencing studies identified notable virulence genes including 30 biofilm genes, 87 motility genes, and 106 secretion system genes. Additionally, a total of 27 antimicrobial resistance genes were found. A multivariate correlation analysis and Scoary analysis based on variation in gene content (n = 812 variable genes) and single nucleotide polymorphisms within virulence genes identified no significant correlations with disease severity, likely due to our limited sample size. In summary, our study explored the virulence and antimicrobial gene content of Pss in pepper seedlings as a first step toward understanding the virulence and pathogenicity of Pss in pepper seedlings. Further studies with additional pepper Pss strains will facilitate defining genes in Pss that correlate with its virulence in pepper seedlings, which can facilitate the development of effective measures to control Pss in pepper and other related P. syringae pathovars.

IMPORTANCE

Pseudomonas leaf spot (PLS) caused by Pseudomonas syringae pv. syringae (Pss) causes significant losses to the pepper industry. Highly virulent Pss strains under optimal environmental conditions (cool-moderate temperatures, high moisture) can cause severe necrotic lesions on pepper leaves that consequently can decrease pepper yield if the disease persists. Hence, it is important to understand the virulence mechanisms of Pss to be able to effectively control PLS in peppers. In our study, in vitro, in planta, and whole-genome sequence analyses were conducted to better understand the virulence and pathogenicity characteristics of Pss strains in peppers. Our findings fill a knowledge gap regarding potential virulence and pathogenicity characteristics of Pss in peppers, including virulence and antimicrobial gene content. Our study helps pave a path to further identify the role of specific virulence genes in causing disease in peppers, which can have implications in developing strategies to effectively control PLS in peppers.

Facultative symbiont virulence determines horizontal transmission rate without host specificity in Dictyostelium discoideum social amoebas

In facultative symbioses, only a fraction of hosts are associated with symbionts. Specific host and symbiont pairings may be the result of host-symbiont coevolution driven by reciprocal selection or priority effects pertaining to which potential symbiont is associated with a host first. Distinguishing between these possibilities is important for understanding the evolutionary forces that affect facultative symbioses. We used the social amoeba, Dictyostelium discoideum, and its symbiont, Paraburkholderia bonniea, to determine whether ongoing coevolution affects which host-symbiont strain pairs naturally cooccur within a facultative symbiosis. Relative to other Paraburkholderia, including another symbiont of D. discoideum, P. bonniea features a reduced genome size that indicates a significant history of coevolution with its host. We hypothesized that ongoing host-symbiont coevolution would lead to higher fitness for naturally cooccurring (native) host and symbiont pairings compared to novel pairings. We show for the first time that P. bonniea symbionts can horizontally transmit to new amoeba hosts when hosts aggregate together during the social stage of their life cycle. Here we find evidence for a virulence-transmission trade-off without host specificity. Although symbiont strains were significantly variable in virulence and horizontal transmission rate, hosts and symbionts responded similarly to associations in native and novel pairings. We go on to identify candidate virulence factors in the genomes of P. bonniea strains that may contribute to variation in virulence. We conclude that ongoing coevolution is unlikely for D. discoideum and P. bonniea. The system instead appears to represent a stable facultative symbiosis in which naturally cooccurring P. bonniea host and symbiont pairings are the result of priority effects.

Detection of Anaplasma and Ehrlichia bacteria in humans, wildlife, and ticks in the Amazon rainforest

Tick-borne bacteria of the genera Ehrlichia and Anaplasma cause several emerging human infectious diseases worldwide. In this study, we conduct an extensive survey for Ehrlichia and Anaplasma infections in the rainforests of the Amazon biome of French Guiana. Through molecular genetics and metagenomics reconstruction, we observe a high indigenous biodiversity of infections circulating among humans, wildlife, and ticks inhabiting these ecosystems. Molecular typing identifies these infections as highly endemic, with a majority of new strains and putative species specific to French Guiana. They are detected in unusual rainforest wild animals, suggesting they have distinctive sylvatic transmission cycles. They also present potential health hazards, as revealed by the detection of Candidatus Anaplasma sparouinense in human red blood cells and that of a new close relative of the human pathogen Ehrlichia ewingii, Candidatus Ehrlichia cajennense, in the tick species that most frequently bite humans in South America. The genome assembly of three new putative species obtained from human, sloth, and tick metagenomes further reveals the presence of major homologs of Ehrlichia and Anaplasma virulence factors. These observations converge to classify health hazards associated with Ehrlichia and Anaplasma infections in the Amazon biome as distinct from those in the Northern Hemisphere.

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.

Paleocene origin of a streamlined digestive symbiosis in leaf beetles

Timing the acquisition of a beneficial microbe relative to the evolutionary history of its host can shed light on the adaptive impact of a partnership. Here, we investigated the onset and molecular evolution of an obligate symbiosis between Cassidinae leaf beetles and Candidatus Stammera capleta, a γ-proteobacterium. Residing extracellularly within foregut symbiotic organs, Stammera upgrades the digestive physiology of its host by supplementing plant cell wall-degrading enzymes. We observe that Stammera is a shared symbiont across tortoise and hispine beetles that collectively comprise the Cassidinae subfamily, despite differences in their folivorous habits. In contrast to its transcriptional profile during vertical transmission, Stammera elevates the expression of genes encoding digestive enzymes while in the foregut symbiotic organs, matching the nutritional requirements of its host. Despite the widespread distribution of Stammera across Cassidinae beetles, symbiont acquisition during the Paleocene (∼62 mya) did not coincide with the origin of the subfamily. Early diverging lineages lack the symbiont and the specialized organs that house it. Reconstructing the ancestral state of host-beneficial factors revealed that Stammera encoded three digestive enzymes at the onset of symbiosis, including polygalacturonase-a pectinase that is universally shared. Although non-symbiotic cassidines encode polygalacturonase endogenously, their repertoire of plant cell wall-degrading enzymes is more limited compared with symbiotic beetles supplemented with digestive enzymes from Stammera. Highlighting the potential impact of a symbiotic condition and an upgraded metabolic potential, Stammera-harboring beetles exploit a greater variety of plants and are more speciose compared with non-symbiotic members of the Cassidinae.

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.

Evolutionary insights into the emergence of virulent Leptospira spirochetes

Pathogenic Leptospira are spirochete bacteria which cause leptospirosis, a re-emerging zoonotic disease of global importance. Here, we use a recently described lineage of environmental-adapted leptospires, which are evolutionarily the closest relatives of the highly virulent Leptospira species, to explore the key phenotypic traits and genetic determinants of Leptospira virulence. Through a comprehensive approach integrating phylogenomic comparisons with in vitro and in vivo phenotyping studies, we show that the evolution towards pathogenicity is associated with both a decrease of the ability to survive in the environment and the acquisition of strategies that enable successful host colonization. This includes the evasion of the human complement system and the adaptations to avoid activation of the innate immune cells. Moreover, our analysis reveals specific genetic determinants that have undergone positive selection during the course of evolution in Leptospira, contributing directly to virulence and host adaptation as demonstrated by gain-of-function and knock-down studies. Taken together, our findings define a new vision on Leptospira pathogenicity, identifying virulence attributes associated with clinically relevant species, and provide insights into the evolution and emergence of these life-threatening pathogens.

Retention of an Endosymbiont for the Production of a Single Molecule

Sap-feeding insects often maintain two or more nutritional endosymbionts that act in concert to produce compounds essential for insect survival. Many mealybugs have endosymbionts in a nested configuration: one or two bacterial species reside within the cytoplasm of another bacterium, and together, these bacteria have genomes that encode interdependent sets of genes needed to produce key nutritional molecules. Here, we show that the mealybug Pseudococcus viburni has three endosymbionts, one of which contributes only two unique genes that produce the host nutrition-related molecule chorismate. All three bacterial endosymbionts have tiny genomes, suggesting that they have been coevolving inside their insect host for millions of years.

Single-cell genomics revealed Candidatus Grellia alia sp. nov. as an endosymbiont of Eutreptiella sp. (Euglenophyceae)

Though endosymbioses between protists and prokaryotes are widespread, certain host lineages have received disproportionate attention what may indicate either a predisposition to such interactions or limited studies on certain protist groups due to lack of cultures. The euglenids represent one such group in spite of microscopic observations showing intracellular bacteria in some strains. Here, we perform a comprehensive molecular analysis of a previously identified endosymbiont in the Eutreptiella sp. CCMP3347 using a single cell approach and bulk culture sequencing. The genome reconstruction of this endosymbiont allowed the description of a new endosymbiont Candidatus Grellia alia sp. nov. from the family Midichloriaceae. Comparative genomics revealed a remarkably complete conjugative type IV secretion system present in three copies on the plasmid sequences of the studied endosymbiont, a feature missing in the closely related Grellia incantans. This study addresses the challenge of limited host cultures with endosymbionts by showing that the genomes of endosymbionts reconstructed from single host cells have the completeness and contiguity that matches or exceeds those coming from bulk cultures. This paves the way for further studies of endosymbionts in euglenids and other protist groups. The research also provides the opportunity to study the diversity of endosymbionts in natural populations.

Speciation Features of Ferdinandcohnia quinoae sp. nov to Adapt to the Plant Host

The bacterial strain SECRCQ15T was isolated from seeds of Chenopodium quinoa in Spain. Phylogenetic, chemotaxonomic, and phenotypic analyses, as well as genome similarity indices, support the classification of the strain into a novel species of the genus Ferdinandcohnia, for which we propose the name Ferdinandcohnia quinoae sp. nov. To dig deep into the speciation features of the strain SECRCQ15T, we performed a comparative genomic analysis of the genome of this strain and those of the type strains of species from the genus Ferdinandcohnia. We found several genes related with plant growth-promoting mechanisms within the SECRCQ15T genome. We also found that singletons of F. quinoae SECRCQ15T are mainly related to the use of carbohydrates, which is a common trait of plant-associated bacteria. To further reveal speciation events in this strain, we revealed genes undergoing diversifying selection (e.g., genes encoding ribosomal proteins) and functions likely lost due to pseudogenization. Also, we found that this novel species contains 138 plant-associated gene-cluster functions that are unique within the genus Ferdinandcohnia. These features may explain both the ecological and taxonomical differentiation of this new taxon.

Comparative genomics of a novel Erwinia species associated with the Highland midge (Culicoides impunctatus)

Erwinia (Enterobacterales: Erwiniaceae) are a group of cosmopolitan bacteria best known as the causative agents of various plant diseases. However, other species in this genus have been found to play important roles as insect endosymbionts supplementing the diet of their hosts. Here, I describe Candidatus Erwinia impunctatus (Erwimp) associated with the Highland midge Culicoides impunctatus (Diptera: Ceratopogonidae), an abundant biting pest in the Scottish Highlands. The genome of this new Erwinia species was assembled using hybrid long and short read techniques, and a comparative analysis was undertaken with other members of the genus to understand its potential ecological niche and impact. Genome composition analysis revealed that Erwimp is similar to other endophytic and ectophytic species in the genus and is unlikely to be restricted to its insect host. Evidence for an additional plant host includes the presence of a carotenoid synthesis operon implicated as a virulence factor in plant-associated members in the sister genus Pantoea. Unique features of Erwimp include several copies of intimin-like proteins which, along with signs of genome pseudogenization and a loss of certain metabolic pathways, suggests an element of host restriction seen elsewhere in the genus. Furthermore, a screening of individuals over two field seasons revealed the absence of the bacteria in Culicoides impunctatus during the second year indicating this microbe-insect interaction is likely to be transient. These data suggest that Culicoides impunctatus may have an important role to play beyond a biting nuisance, as an insect vector transmitting Erwimp alongside any conferred impacts to surrounding biota.

Emergence of Poultry-Associated Human Salmonella enterica Serovar Abortusovis Infections, New South Wales, Australia

Salmonella enterica serovar Abortusovis is a ovine-adapted pathogen that causes spontaneous abortion. Salmonella Abortusovis was reported in poultry in 2009 and has since been reported in human infections in New South Wales, Australia. Phylogenomic analysis revealed a clade of 51 closely related isolates from Australia originating in 2004. That clade was genetically distinct from ovine-associated isolates. The clade was widespread in New South Wales poultry production facilities but was only responsible for sporadic human infections. Some known virulence factors associated with human infections were only found in the poultry-associated clade, some of which were acquired through prophages and plasmids. Furthermore, the ovine-associated clade showed signs of genome decay, but the poultry-associated clade did not. Those genomic changes most likely led to differences in host range and disease type. Surveillance using the newly identified genetic markers will be vital for tracking Salmonella Abortusovis transmission in animals and to humans and preventing future outbreaks.

Draft genome sequence of Wigglesworthia glossinidia "palpalis gambiensis" isolate

The 0.719 Mb genome of the tsetse endosymbiont, Wigglesworthia glossinidia, from Glossina palpalis gambiensis is presented. This Wigglesworthia genome retains 611 protein-coding sequences and a 25.3% GC content. A cryptic plasmid is conserved, between Wigglesworthia isolates, suggesting functional significance. This genome adds a further dimension to characterize Wigglesworthia lineage-based differences.

Pseudogenes act as a neutral reference for detecting selection in prokaryotic pangenomes

A long-standing question is to what degree genetic drift and selection drive the divergence in rare accessory gene content between closely related bacteria. Rare genes, including singletons, make up a large proportion of pangenomes (all genes in a set of genomes), but it remains unclear how many such genes are adaptive, deleterious or neutral to their host genome. Estimates of species' effective population sizes (Ne) are positively associated with pangenome size and fluidity, which has independently been interpreted as evidence for both neutral and adaptive pangenome models. We hypothesized that pseudogenes, used as a neutral reference, could be used to distinguish these models. We find that most functional categories are depleted for rare pseudogenes when a genome encodes only a single intact copy of a gene family. In contrast, transposons are enriched in pseudogenes, suggesting they are mostly neutral or deleterious to the host genome. Thus, even if individual rare accessory genes vary in their effects on host fitness, we can confidently reject a model of entirely neutral or deleterious rare genes. We also define the ratio of singleton intact genes to singleton pseudogenes (si/sp) within a pangenome, compare this measure across 668 prokaryotic species and detect a signal consistent with the adaptive value of many rare accessory genes. Taken together, our work demonstrates that comparing with pseudogenes can improve inferences of the evolutionary forces driving pangenome variation.

Pseudogenomic insights into the evolution of Mycobacterium ulcerans

BACKGROUND

Buruli ulcer (BU) disease, caused by Mycobacterium ulcerans (MU), and characterized by necrotic ulcers is still a health problem in Africa and Australia. The genome of the bacterium has several pseudogenes due to recent evolutionary events and environmental pressures. Pseudogenes are genetic elements regarded as nonessential in bacteria, however, they are less studied due to limited available tools to provide understanding of their evolution and roles in MU pathogenicity.

RESULTS

This study developed a bioinformatic pipeline to profile the pseudogenomes of sequenced MU clinical isolates from different countries. One hundred and seventy-two MU genomes analyzed revealed that pseudogenomes of African strains corresponded to the two African lineages 1 and 2. Pseudogenomes were lineage and location specific and African lineage 1 was further divided into A and B. Lineage 2 had less relaxation in positive selection than lineage 1 which may signify different evolutionary points. Based on the Gil-Latorre model, African MU strains may be in the latter stages of evolutionary adaption and are adapting to an environment rich in metabolic resources with a lower temperature and decreased UV radiation. The environment fosters oxidative metabolism and MU may be less reliant on some secondary metabolites. In-house pseudogenomes from Ghana and Cote d'Ivoire were different from other African strains, however, they were identified as African strains.

CONCLUSION

Our bioinformatic pipeline provides pseudogenomic insights to complement other whole genome analyses, providing a better view of the evolution of the genome of MU and suggest an adaptation model which is important in understanding transmission. MU pseudogene profiles vary based on lineage and country, and an apparent reduction in insertion sequences used for the detection of MU which may adversely affect the sensitivity of diagnosis.

Gene content, phage cycle regulation model and prophage inactivation disclosed by prophage genomics in the Helicobacter pylori Genome Project

Prophages can have major clinical implications through their ability to change pathogenic bacterial traits. There is limited understanding of the prophage role in ecological, evolutionary, adaptive processes and pathogenicity of Helicobacter pylori, a widespread bacterium causally associated with gastric cancer. Inferring the exact prophage genomic location and completeness requires complete genomes. The international Helicobacter pylori Genome Project (HpGP) dataset comprises 1011 H. pylori complete clinical genomes enriched with epigenetic data. We thoroughly evaluated the H. pylori prophage genomic content in the HpGP dataset. We investigated population evolutionary dynamics through phylogenetic and pangenome analyses. Additionally, we identified genome rearrangements and assessed the impact of prophage presence on bacterial gene disruption and methylome. We found that 29.5% (298) of the HpGP genomes contain prophages, of which only 32.2% (96) were complete, minimizing the burden of prophage carriage. The prevalence of H. pylori prophage sequences was variable by geography and ancestry, but not by disease status of the human host. Prophage insertion occasionally results in gene disruption that can change the global bacterial epigenome. Gene function prediction allowed the development of the first model for lysogenic-lytic cycle regulation in H. pylori. We have disclosed new prophage inactivation mechanisms that appear to occur by genome rearrangement, merger with other mobile elements, and pseudogene accumulation. Our analysis provides a comprehensive framework for H. pylori prophage biological and genomics, offering insights into lysogeny regulation and bacterial adaptation to prophages.

Facultative endosymbiosis between cellulolytic protists and methanogenic archaea in the gut of the Formosan termite Coptotermes formosanus

Anaerobic protists frequently harbour methanogenic archaea, which apparently contribute to the hosts' fermentative metabolism by consuming excess H2. However, the ecological properties of endosymbiotic methanogens remain elusive in many cases. Here we investigated the ecology and genome of the endosymbiotic methanogen of the Cononympha protists in the hindgut of the termite Coptotermes formosanus. Microscopic and 16S rRNA amplicon sequencing analyses revealed that a single species, designated here "Candidatus Methanobrevibacter cononymphae", is associated with both Cononympha leidyi and Cononympha koidzumii and that its infection rate in Cononympha cells varied from 0.0% to 99.8% among termite colonies. Fine-scale network analysis indicated that multiple 16S rRNA sequence variants coexisted within a single host cell and that identical variants were present in both Cononympha species and also on the gut wall. Thus, "Ca. Methanobrevibacter cononymphae" is a facultative endosymbiont, transmitted vertically with frequent exchanges with the gut environment. Indeed, transmission electron microscopy showed escape or uptake of methanogens from/by a Cononympha cell. The genome of "Ca. Methanobrevibacter cononymphae" showed features consistent with its facultative lifestyle: i.e., the genome size (2.7 Mbp) comparable to those of free-living relatives; the pseudogenization of the formate dehydrogenase gene fdhA, unnecessary within the non-formate-producing host cell; the dependence on abundant acetate in the host cell as an essential carbon source; and the presence of a catalase gene, required for colonization on the microoxic gut wall. Our study revealed a versatile endosymbiosis between the methanogen and protists, which may be a strategy responding to changing conditions in the termite gut.

Multiple independent losses of sporulation and peptidoglycan in the Mycoplasmatales and related orders of the class Bacilli

Many peptidoglycan-deficient bacteria such as the Mycoplasmatales are known host-associated lineages, lacking the environmental resistance mechanisms and metabolic capabilities necessary for a free-living lifestyle. Several peptidoglycan-deficient and non-sporulating orders of interest are thought to be descended from Gram-positive sporulating Bacilli through reductive evolution. Here we annotate 2650 genomes belonging to the class Bacilli, according to the Genome Taxonomy Database, to predict the peptidoglycan and sporulation phenotypes of three novel orders, RFN20, RF39 and ML615J-28, known only through environmental sequence surveys. These lineages are interspersed between peptidoglycan-deficient non-sporulating orders including the Mycoplasmatales and Acholeplasmatales, and more typical Gram-positive orders such as the Erysipelotrichales and Staphylococcales. We use the extant genotypes to perform ancestral state reconstructions. The novel orders are predicted to have small genomes with minimal metabolic capabilities and to comprise a mix of peptidoglycan-deficient and/or non-sporulating species. In contrast to expectations based on cultured representatives, the order Erysipelotrichales lacks many of the genes involved in peptidoglycan and endospore formation. The reconstructed evolutionary history of these traits suggests multiple independent whole-genome reductions and loss of phenotype via intermediate transition states that continue into the present. We suggest that the evolutionary history of the reduced-genome lineages within the class Bacilli is one driven by multiple independent transitions to host-associated lifestyles, with the degree of reduction in environmental resistance and metabolic capabilities correlated with degree of host association.

Pantoea bathycoeliae sp. nov and Sodalis sp. are core gut microbiome symbionts of the two-spotted stink bug

Stink bug species (Pentatomoidea superfamily) have developed an interdependence with obligate bacterial gut symbionts in specialized midgut crypts (M4 sub-region). Species of the Enterobacteriaceae family (predominantly Pantoea) are vertically transferred to their offspring and provide nutrients that cannot be obtained from plant sap food sources. However, the bacteria in the other gut compartments of stink bugs have rarely been investigated. The two-spotted stink bug, Bathycoelia distincta, is a serious pest of macadamias in South Africa. Nothing is currently known regarding its gut microbiome or how symbionts are transferred between insect generations. In this study, the consistency of B. distincta gut bacteria across geographic locations and life stages was determined with 16S rRNA metabarcoding, considering both the M4 and other gut compartments. A novel Pantoea species was found to be the primary M4 gut symbiont and is vertically transferred to the offspring. The other gut compartments had a low bacterial diversity and genera varied between stink bug populations but a Sodalis species was prominent in all populations. Sequence data of the M4 compartment were used to produce high-quality metagenome-assembled genomes (MAGs) for the Pantoea and Sodalis species. Functional analyses suggested a similar role in nutrient provision for the host, yet also unique metabolites produced by each species. The Sodalis sp. also had additional traits, such as secretion systems, that likely allowed it to establish itself in the host. The Pantoea species was described as Pantoea bathycoeliae sp. nov based on the rules of the SeqCode.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

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.

Comparative genomics of Mollicutes-related endobacteria supports a late invasion into Mucoromycota fungi

Diverse members of early-diverging Mucoromycota, including mycorrhizal taxa and soil-associated Mortierellaceae, are known to harbor Mollicutes-related endobacteria (MRE). It has been hypothesized that MRE were acquired by a common ancestor and transmitted vertically. Alternatively, MRE endosymbionts could have invaded after the divergence of Mucoromycota lineages and subsequently spread to new hosts horizontally. To better understand the evolutionary history of MRE symbionts, we generated and analyzed four complete MRE genomes from two Mortierellaceae genera: Linnemannia (MRE-L) and Benniella (MRE-B). These genomes include the smallest known of fungal endosymbionts and showed signals of a tight relationship with hosts including a reduced functional capacity and genes transferred from fungal hosts to MRE. Phylogenetic reconstruction including nine MRE from mycorrhizal fungi revealed that MRE-B genomes are more closely related to MRE from Glomeromycotina than MRE-L from the same host family. We posit that reductions in genome size, GC content, pseudogene content, and repeat content in MRE-L may reflect a longer-term relationship with their fungal hosts. These data indicate Linnemannia and Benniella MRE were likely acquired independently after their fungal hosts diverged from a common ancestor. This work expands upon foundational knowledge on minimal genomes and provides insights into the evolution of bacterial endosymbionts.

Genomic Characteristics of Emerging Intraerythrocytic Anaplasma capra and High Prevalence in Goats, China

Anaplasma capra is an emerging tickborne human pathogen initially recognized in China in 2015; it has been reported in ticks and in a wide range of domestic and wild animals worldwide. We describe whole-genome sequences of 2 A. capra strains from metagenomic sequencing of purified erythrocytes from infected goats in China. The genome of A. capra was the smallest among members of the genus Anaplasma. The genomes of the 2 A. capra strains contained comparable G+C content and numbers of pseudogenes with intraerythrocytic Anaplasma species. The 2 A. capra strains had 54 unique genes. The prevalence of A. capra was high among goats in the 2 endemic areas. Phylogenetic analyses revealed that the A. capra strains detected in this study were basically classified into 2 subclusters with those previously detected in Asia. Our findings clarify details of the genomic characteristics of A. capra and shed light on its genetic diversity.

A single cryptomonad cell harbors a complex community of organelles, bacteria, a phage, and selfish elements

Symbiosis between prokaryotes and microbial eukaryotes (protists) has broadly impacted both evolution and ecology. Endosymbiosis led to mitochondria and plastids, the latter spreading across the tree of eukaryotes by subsequent rounds of endosymbiosis. Present-day endosymbionts in protists remain both common and diverse, although what function they serve is often unknown. Here, we describe a highly complex community of endosymbionts and a bacteriophage (phage) within a single cryptomonad cell. Cryptomonads are a model for organelle evolution because their secondary plastid retains a relict endosymbiont nucleus, but only one previously unidentified Cryptomonas strain (SAG 25.80) is known to harbor bacterial endosymbionts. We carried out electron microscopy and FISH imaging as well as genomic sequencing on Cryptomonas SAG 25.80, which revealed a stable, complex community even after over 50 years in continuous cultivation. We identified the host strain as Cryptomonas gyropyrenoidosa, and sequenced genomes from its mitochondria, plastid, and nucleomorph (and partially its nucleus), as well as two symbionts, Megaira polyxenophila and Grellia numerosa, and one phage (MAnkyphage) infecting M. polyxenophila. Comparing closely related endosymbionts from other hosts revealed similar metabolic and genomic features, with the exception of abundant transposons and genome plasticity in M. polyxenophila from Cryptomonas. We found an abundance of eukaryote-interacting genes as well as many toxin-antitoxin systems, including in the MAnkyphage genome that also encodes several eukaryotic-like proteins. Overall, the Cryptomonas cell is an endosymbiotic conglomeration with seven distinct evolving genomes that all show evidence of inter-lineage conflict but nevertheless remain stable, even after more than 4,000 generations in culture.

Supergroup F Wolbachia with extremely reduced genome: transition to obligate insect symbionts

BACKGROUND

Wolbachia belong to highly abundant bacteria which are frequently found in invertebrate microbiomes and manifest by a broad spectrum of lifestyles from parasitism to mutualism. Wolbachia supergroup F is a particularly interesting clade as it gave rise to symbionts of both arthropods and nematodes, and some of its members are obligate mutualists. Investigations on evolutionary transitions among the different symbiotic stages have been hampered by a lack of the known diversity and genomic data for the supergroup F members.

RESULTS

Based on amplicon screening, short- and long-read WGS approaches, and laser confocal microscopy, we characterize five new supergroup F Wolbachia strains from four chewing lice species. These strains reached different evolutionary stages and represent two remarkably different types of symbiont genomes. Three of the genomes resemble other known members of Wolbachia F supergroup, while the other two show typical signs of ongoing gene inactivation and removal (genome size, coding density, low number of pseudogenes). Particularly, wMeur1, a symbiont fixed in microbiomes of Menacanthus eurysternus across four continents, possesses a highly reduced genome of 733,850 bp. The horizontally acquired capacity for pantothenate synthesis and localization in specialized bacteriocytes suggest its obligate nutritional role.

CONCLUSIONS

The genome of wMeur1 strain, from the M. eurysternus microbiome, represents the smallest currently known Wolbachia genome and the first example of Wolbachia which has completed genomic streamlining as known from the typical obligate symbionts. This points out that despite the large amount and great diversity of the known Wolbachia strains, evolutionary potential of these bacteria still remains underexplored. The diversity of the four chewing lice microbiomes indicates that this vast parasitic group may provide suitable models for further investigations. Video Abstract.

Genomic analysis of Mycobacterium brumae sustains its nonpathogenic and immunogenic phenotype

Mycobacterium brumae is a rapid-growing, non-pathogenic Mycobacterium species, originally isolated from environmental and human samples in Barcelona, Spain. Mycobacterium brumae is not pathogenic and it's in vitro phenotype and immunogenic properties have been well characterized. However, the knowledge of its underlying genetic composition is still incomplete. In this study, we first describe the 4 Mb genome of the M. brumae type strain ATCC 51384T assembling PacBio reads, and second, we assess the low intraspecies variability by comparing the type strain with Illumina reads from three additional strains. Mycobacterium brumae genome is composed of a circular chromosome with a high GC content of 69.2% and containing 3,791 CDSs, 97 pseudogenes, one prophage and no CRISPR loci. Mycobacterium brumae has shown no pathogenic potential in in vivo experiments, and our genomic analysis confirms its phylogenetic position with other non-pathogenic and rapid growing mycobacteria. Accordingly, we determined the absence of virulence-related genes, such as ESX-1 locus and most PE/PPE genes, among others. Although the immunogenic potential of M. brumae was proved to be as high as Mycobacterium bovis BCG, the only mycobacteria licensed to treat cancer, the genomic content of M. tuberculosis T cell and B cell antigens in M. brumae genome is considerably lower than those antigens present in M. bovis BCG genome. Overall, this work provides relevant genomic data on one of the species of the mycobacterial genus with high therapeutic potential.

Complex host/symbiont integration of a multi-partner symbiotic system in the eusocial aphid Ceratovacuna japonica

Some hemipteran insects rely on multiple endosymbionts for essential nutrients. However, the evolution of multi-partner symbiotic systems is not well-established. Here, we report a co-obligate symbiosis in the eusocial aphid, Ceratovacuna japonica. 16S rRNA amplicon sequencing unveiled co-infection with a novel Arsenophonus sp. symbiont and Buchnera aphidicola, a common obligate endosymbiont in aphids. Both symbionts were housed within distinct bacteriocytes and were maternally transmitted. The Buchnera and Arsenophonus symbionts had streamlined genomes of 432,286 bp and 853,149 bp, respectively, and exhibited metabolic complementarity in riboflavin and peptidoglycan synthesis pathways. These anatomical and genomic properties were similar to those of independently evolved multi-partner symbiotic systems, such as Buchnera-Serratia in Lachninae and Periphyllus aphids, representing remarkable parallelism. Furthermore, symbiont populations and bacteriome morphology differed between reproductive and soldier castes. Our study provides the first example of co-obligate symbiosis in Hormaphidinae and gives insight into the evolutionary genetics of this complex system.

At the threshold of symbiosis: the genome of obligately endosymbiotic 'Candidatus Nebulobacter yamunensis' is almost indistinguishable from that of a cultivable strain

Comparing obligate endosymbionts with their free-living relatives is a powerful approach to investigate the evolution of symbioses, and it has led to the identification of several genomic traits consistently associated with the establishment of symbiosis. 'Candidatus Nebulobacter yamunensis' is an obligate bacterial endosymbiont of the ciliate Euplotes that seemingly depends on its host for survival. A subsequently characterized bacterial strain with an identical 16S rRNA gene sequence, named Fastidiosibacter lacustris, can instead be maintained in pure culture. We analysed the genomes of 'Candidatus Nebulobacter' and Fastidiosibacter seeking to identify key differences between their functional traits and genomic structure that might shed light on a recent transition to obligate endosymbiosis. Surprisingly, we found almost no such differences: the two genomes share a high level of sequence identity, the same overall structure, and largely overlapping sets of genes. The similarities between the genomes of the two strains are at odds with their different ecological niches, confirmed here with a parallel growth experiment. Although other pairs of closely related symbiotic/free-living bacteria have been compared in the past, 'Candidatus Nebulobacter' and Fastidiosibacter represent an extreme example proving that a small number of (unknown) factors might play a pivotal role in the earliest stages of obligate endosymbiosis establishment.

Reduced and Nonreduced Genomes in Paraburkholderia Symbionts of Social Amoebas

The social amoeba Dictyostelium discoideum is a predatory soil protist frequently used for studying host-pathogen interactions. A subset of D. discoideum strains isolated from soil persistently carry symbiotic Paraburkholderia, recently formally described as P. agricolaris, P. bonniea, and P. hayleyella. The three facultative symbiont species of D. discoideum present a unique opportunity to study a naturally occurring symbiosis in a laboratory model protist. There is a large difference in genome size between P. agricolaris (8.7 million base pairs [Mbp]) versus P. hayleyella and P. bonniea (4.1 Mbp). We took a comparative genomics approach and compared the three genomes of D. discoideum symbionts to 12 additional Paraburkholderia genomes to test for genome evolution patterns that frequently accompany host adaptation. Overall, P. agricolaris is difficult to distinguish from other Paraburkholderia based on its genome size and content, but the reduced genomes of P. bonniea and P. hayleyella display characteristics indicative of genome streamlining rather than deterioration during adaptation to their protist hosts. In addition, D. discoideum-symbiont genomes have increased secretion system and motility genes that may mediate interactions with their host. Specifically, adjacent BurBor-like type 3 and T6SS-5-like type 6 secretion system operons shared among all three D. discoideum-symbiont genomes may be important for host interaction. Horizontal transfer of these secretion system operons within the amoeba host environment may have contributed to the unique ability of these symbionts to establish and maintain a symbiotic relationship with D. discoideum. IMPORTANCE Protists are a diverse group of typically single cell eukaryotes. Bacteria and archaea that form long-term symbiotic relationships with protists may evolve in additional ways than those in relationships with multicellular eukaryotes such as plants, animals, or fungi. Social amoebas are a predatory soil protist sometimes found with symbiotic bacteria living inside their cells. They present a unique opportunity to explore a naturally occurring symbiosis in a protist frequently used for studying host-pathogen interactions. We show that one amoeba-symbiont species is similar to other related bacteria in genome size and content, while the two reduced-genome-symbiont species show characteristics of genome streamlining rather than deterioration during adaptation to their host. We also identify sets of genes present in all three amoeba-symbiont genomes that are potentially used for host-symbiont interactions. Because the amoeba symbionts are distantly related, the amoeba host environment may be where these genes were shared among symbionts.

Brucella ceti and Brucella pinnipedialis genome characterization unveils genetic features that highlight their zoonotic potential

The Gram-negative bacteria Brucella ceti and Brucella pinnipedialis circulate in marine environments primarily infecting marine mammals, where they cause an often-fatal disease named brucellosis. The increase of brucellosis among several species of cetaceans and pinnipeds, together with the report of sporadic human infections, raises concerns about the zoonotic potential of these pathogens on a large scale and may pose a threat to coastal communities worldwide. Therefore, the characterization of the B. ceti and B. pinnipedialis genetic features is a priority to better understand the pathological factors that may impact global health. Moreover, an in-depth functional analysis of the B. ceti and B. pinnipedialis genome in the context of virulence and pathogenesis was not undertaken so far. Within this picture, here we present the comparative whole-genome characterization of all B. ceti and B. pinnipedialis genomes available in public resources, uncovering a collection of genetic tools possessed by these aquatic bacterial species compared to their zoonotic terrestrial relatives. We show that B. ceti and B. pinnipedialis genomes display a wide host-range infection capability and a polyphyletic phylogeny within the genus, showing a genomic structure that fits the canonical definition of closeness. Functional genome annotation led to identifying genes related to several pathways involved in mechanisms of infection, others conferring pan-susceptibility to antimicrobials and a set of virulence genes that highlight the similarity of B. ceti and B. pinnipedialis genotypes to those of Brucella spp. displaying human-infecting phenotypes.