The nature and dynamics of bacterial genomes

Gene pseudogenization in fertility-associated genes in cheetah (Acinonyx jubatus), a species with long-term low effective population size

We are witnessing an ongoing global biodiversity crisis, and an increasing number of mammalian populations are at risk of decline. Species that have survived severe historic bottlenecks, such as the cheetah (Acinonyx jubatus) exhibit symptoms of inbreeding depression including reproductive and developmental defects. Although it has long been suggested that such defects stem from an accumulation of weakly deleterious mutations, the implications of such mutations leading to pseudogenization has not been assessed. Here, we use comparative analysis of eight felid genomes to better understand the impacts of deleterious mutations in the cheetah. We find novel pseudogenization events specific to the cheetah. Through careful curation, we identify 65 genes with previously unreported premature termination codons (PTCs) that likely affect gene function. With the addition of population data (n = 6), we find 22 of these PTCs in at least one resequenced individual, four of which (DEFB116, ARL13A, CFAP119, and NT5DC4) are also found in a more recent reference genome. Mutations within three of these genes are linked with sterility, including azoospermia, which is common in cheetahs. Our results highlight the power of comparative genomic approaches for the discovery of novel causative variants in declining species.

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.

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.

Live tracking of a plant pathogen outbreak reveals rapid and successive, multidecade plasmid reduction

Quickly understanding the genomic changes that lead to pathogen emergence is necessary to launch mitigation efforts and reduce harm. In this study, we tracked in real time a 2022 bacterial plant disease outbreak in U.S. geraniums (Pelargonium × hortorum) caused by Xhp2022, a novel lineage of Xanthomonas hortorum. Genomes from 31 Xhp2022 isolates from seven states showed limited chromosomal variation and all contained a single plasmid (p93). Time tree and single nucleotide polymorphism whole-genome analysis estimated that Xhp2022 emerged within the last decade. The phylogenomic analysis determined that p93 resulted from the cointegration of three plasmids (p31, p45, and p66) found sporadically across isolates from previous outbreaks. Although p93 had a 49 kb nucleotide reduction, it retained putative fitness genes, which became predominant in the 2022 outbreak. Overall, we demonstrated, through rapid whole-genome sequencing and analysis, a recent, traceable event of genome reduction for niche adaptation typically observed over millennia in obligate and fastidious pathogens.IMPORTANCEThe geranium industry, valued at $4 million annually, faces an ongoing Xanthomonas hortorum pv. pelargonii (Xhp) pathogen outbreak. To track and describe the outbreak, we compared the genome structure across historical and globally distributed isolates. Our research revealed Xhp population has not had chromosome rearrangements since 1974 and has three distinct plasmids. In 2012, we found all three plasmids in individual Xhp isolates. However, in 2022, the three plasmids co-integrated into one plasmid named p93. p93 retained putative fitness genes but lost extraneous genomic material. Our findings show that the 2022 strain group of the bacterial plant pathogen Xanthomonas hortorum underwent a plasmid reduction. We also observed several Xanthomonas species from different years, hosts, and continents have similar plasmids to p93, possibly due to shared agricultural settings. We noticed parallels between genome efficiency and reduction that we see across millennia with obligate parasites with increased niche specificity.

Contraction and expansion dynamics: deciphering genomic underpinnings of growth rate and pathogenicity in Mycobacterium

BACKGROUND

Mycobacterium bacteria, encompassing both slow growth (SGM) and rapid growth mycobacteria (RGM), along with true pathogenic (TP), opportunistic pathogenic (OP), and non-pathogenic (NP) types, exhibit diverse phenotypes. Yet, the genetic underpinnings of these variations remain elusive.

METHODS

Here, We conducted a comprehensive comparative genomics study involving 53 Mycobacterium species to unveil the genomic drivers behind growth rate and pathogenicity disparities.

RESULTS

Our core/pan-genome analysis highlighted 1,307 shared gene families, revealing an open pan-genome structure. A phylogenetic tree highlighted clear boundaries between SGM and RGM, as well as TP and other species. Gene family contraction emerged as the primary alteration associated with growth and pathogenicity transitions. Specifically, ABC transporters for amino acids and inorganic ions, along with quorum sensing genes, exhibited significant contractions in SGM species, potentially influencing their distinct traits. Conversely, TP strains displayed contraction in lipid and secondary metabolite biosynthesis and metabolism-related genes. Across the 53 species, we identified 26 core and 64 accessory virulence factors. Remarkably, TP and OP strains stood out for their expanded mycobactin biosynthesis and type VII secretion system gene families, pivotal for their pathogenicity.

CONCLUSION

Our findings underscore the importance of gene family contraction in nucleic acids, ions, and substance metabolism for host adaptation, while emphasizing the significance of virulence gene family expansion, including type VII secretion systems and mycobactin biosynthesis, in driving mycobacterial pathogenicity.

CPGminer: An Interactive Dashboard to Explore the Genomic Features and Taxonomy of Complete Prokaryotic Genomes

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

Pangenome analysis provides insights into the genetic diversity, metabolic versatility, and evolution of the genus Flavobacterium

Members of the genus Flavobacterium are widely distributed and produce various polysaccharide-degrading enzymes. Many species in the genus have been isolated and characterized. However, few studies have focused on marine isolates or fish pathogens, and in-depth genomic analyses, particularly comparative analyses of isolates from different habitat types, are lacking. Here, we isolated 20 strains of the genus from various environments in South Korea and sequenced their full-length genomes. Combined with published sequence data, we examined genomic traits, evolution, environmental adaptation, and putative metabolic functions in total 187 genomes of isolated species in Flavobacterium categorized as marine, host-associated, and terrestrial including freshwater. A pangenome analysis revealed a correlation between genome size and coding or noncoding density. Flavobacterium spp. had high levels of diversity, allowing for novel gene repertories via recombination events. Defense-related genes only accounted for approximately 3% of predicted genes in all Flavobacterium genomes. While genes involved in metabolic pathways did not differ with respect to isolation source, there was substantial variation in genomic traits; in particular, the abundances of tRNAs and rRNAs were higher in the host-associdated group than in other groups. One genome in the host-associated group contained a Microviridae prophage closely related to an enterobacteria phage. The proteorhodopsin gene was only identified in four terrestrial strains isolated for this study. Furthermore, recombination events clearly influenced genomic diversity and may contribute to the response to environmental stress. These findings shed light on the high genetic variation in Flavobacterium and functional roles in diverse ecosystems as a result of their metabolic versatility. IMPORTANCE The genus Flavobacterium is a diverse group of bacteria that are found in a variety of environments. While most species of this genus are harmless and utilize organic substrates such as proteins and polysaccharides, some members may play a significant role in the cycling for organic substances within their environments. Nevertheless, little is known about the genomic dynamics and/or metabolic capacity of Flavobacterium. Here, we found that Flavobacterium species may have an open pangenome, containing a variety of diverse and novel gene repertoires. Intriguingly, we discovered that one genome (classified into host-associated group) contained a Microviridae prophage closely related to that of enterobacteria. Proteorhodopsin may be expressed under conditions of light or oxygen pressure in some strains isolated for this study. Our findings significantly contribute to the understanding of the members of the genus Flavobacterium diversity exploration and will provide a framework for the way for future ecological characterizations.

Genomic distribution of the insertion sequence IS711 reveal a potential role in Brucella genome plasticity and host preference

The Insertion Sequence 711 (IS711) is linked to the Brucella genus. Mapping the genomic distribution of IS711 can help understand this insertion element's biological and evolutionary role. This work aimed to delineate the genomic distribution of the IS711 element and to study its association with Brucella evolution. A total of 124 genomes representing 9 Brucella species were searched using BLASTn sequence alignment tool to identify complete and truncated copies of IS711. Based on the genomic context, each IS711 locus was assigned a code using the initial letters of its neighboring genes. Various tools were used to annotate the neighboring genes and determine the shared synteny around orthologous IS711 loci. The tool Islandviewer 4 was used to scan for genomic islands. The Codon Tree method was used to build phylogenetic trees of B. melitensis, B. abortus, and B. suis genomes. The phylogenetic trees of the three species were analyzed, taking into account the genomic distribution patterns of IS711. The result of IS711 frequency analysis showed a relatively conserved number of copies/genome for the different species and for some biovars. The analysis showed that Brucella species with a relatively low IS711 copy number (4-8 copies/genome) are linked to domestic animals as primary hosts and have potential for zoonotic transmission. However, species with a relatively higher copy number (12-30 copies/genome) are less zoonotic and tend to be linked with wild animals as primary hosts. Analyzing the genomic distribution map of IS711 loci showed several unique patterns of IS711 distribution that are correlated with the evolution of Brucella species and biovars. The results also showed that 46.2% of the conserved IS711 elements are located within genomic islands. Based on our results and previous data, we postulate a model explaining the IS711 role in Brucella evolution. We assume that during the transition from a free-living to an intracellular lifestyle, a descendant of the Brucella genus had acquired a progenitor sequence of the IS711. Subsequently, a burst in IS711 transposition occurred. This parasitic expansion can be deleterious and has to be counteracted by evolutionary forces to prevent lineage extension and to promote adaptation to host. Similar to other plasmid-free pathogenic α-Proteobacteria bacteria, the balance of expansion and reduction of insertion elements could be one of the mechanisms to control genome reduction and streamlining. We hypothesize that the IS711-mediated genomic changes and other small sequence nucleotide changes in specific orthologous genes could significantly contribute to Brucella's evolution and adaptation to different animal hosts.

Genomic and metabolic features of Bacillus cereus, inhibiting the growth of Sclerotinia sclerotiorum by synthesizing secondary metabolites

We investigated the biocontrol mechanism of Bacillus cereus CF4-51 to find powerful microbes that effectively control Sclerotinia sclerotiorum. To assess its inhibitory effect on fungal growth, the plant pathogen (S. sclerotiorum) was co-cultured with Bacillus cereus. Scanning electron microscope (SEM) was used to study the morphology of S. sclerotiorum treated with CF4-51 biofumigant. The expression of sclerotium formation-related genes was analyzed by qRT-PCR. We performed whole genome sequencing of CF4-51 by PacBio Sequel platform. Lipopeptides were extracted from strain CF4-51 according to the method of hydrochloric acid precipitation and methanol dissolution. The volatiles CF4-51 were identified using gas chromatography-mass spectrometry (GC-MS). We found that the volatile organic compounds (VOCs) released by CF4-51 damaged the S. sclerotiorum hyphae and inhibited the formation of sclerotia. The qRT-PCR data revealed the down-regulated expression of the genes involved in sclerotial formation. Moreover, we analyzed the B. cereus CF4-51 genome and metabolites. The genome consisted of 5.35 Mb, with a GC content of 35.74%. An examination of the genome revealed the presence of several gene clusters for the biosynthesis of antibiotics, siderophores, and various other bioactive compounds, including those belonging to the NRPS-like, LAP, RIPP-like, NRPS, betalactone, CDPS, terpene, ladderane, ranthipeptide, and lanthipeptide (class II) categories. A gas chromatography-tandem mass spectrometry analysis identified 45 VOCs produced by strain CF4-51. Among these, technical grade formulations of five were chosen for further study: 2-Pentadecanone, 6,10,14-trimethyl-,1,2-Benzenedicarboxylic acid, bis(2-methylpropyl) ester, Dibutyl phthalate, Cyclododecane, Heptadecane. the five major constituents play important roles in the antifungal activity of the VOCs CF4-51 on the growth of S. sclerotiorum. The secondary metabolites produced by strain CF4-51are critical for the inhibition of S. sclerotiorum hyphal growth and sclerotial formation.

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.

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.

Description of three novel species of Scandinavium: Scandinavium hiltneri sp. nov., Scandinavium manionii sp. nov. and Scandinavium tedordense sp. nov., isolated from the oak rhizosphere and bleeding cankers of broadleaf hosts

While investigating the bacterial populations of environmental samples taken from a mix of healthy and Acute Oak Decline afflicted Quercus robur (pedunculate or English oak) rhizosphere soil samples and swabs of bleeding lesions on Tilia spp. (lime) and Quercus rubra (red oak) trees, several strains belonging to the order Enterobacterales were isolated using selective media and enrichment broth. Seven strains from the Q. robur rhizosphere, three strains from Tilia spp. and one from Q. rubra were investigated, with their taxonomic status determined via a polyphasic taxonomic approach. Initially stains were identified as potential members of the recently described genus Scandinavium, based on the partial sequencing of three housekeeping genes. Further analysis of phenotypic traits, including fatty acid profiles, coupled with 16S rRNA gene and phylogenomic analysis of whole genome sequences were applied to a subset of the strains. Phylogenetic and phylogenomic analysis repeatedly placed the isolates in a monophyletic clade within Scandinavium, with four distinct clusters observed, one of which corresponded to Scandinavium goeteborgense, the type species of the genus. The remaining three clusters could be phenotypically and genotypically differentiated from each other and S. goeteborgense. As such, we describe three novel species of the genus, for which we propose the names Scandinavium hiltneri sp. nov. (type strain H11S7^T = LMG 32612^T = CCUG 76179^T), Scandinavium manionii sp. nov. (type strain H17S15^T = LMG 32613^T = CCUG 76183^T) and Scandinavium tedordense sp. nov. (type strain TWS1a^T = LMG 32614^T = CCUG 76188^T). Additionally, the descriptions of the genus Scandinavium and the type species, S. goeteborgense, are emended.

Bacterioplankton communities reveal horizontal and vertical influence of an Island Mass Effect

Coral reefs are highly productive ecosystems with distinct biogeochemistry and biology nestled within unproductive oligotrophic gyres. Coral reef islands have often been associated with a nearshore enhancement in phytoplankton, a phenomenon known as the Island Mass Effect (IME). Despite being documented more than 60 years ago, much remains unknown about the extent and drivers of IMEs. Here we utilized 16S rRNA gene metabarcoding as a biological tracer to elucidate horizontal and vertical influence of an IME around the islands of Mo'orea and Tahiti, French Polynesia. We show that those nearshore oceanic stations with elevated chlorophyll a included bacterioplankton found in high abundance in the reef environment, suggesting advection of reef water is the source of altered nearshore biogeochemistry. We also observed communities in the nearshore deep chlorophyll maximum (DCM) with enhanced abundances of upper euphotic bacterioplankton that correlated with intrusions of low-density, O₂ rich water, suggesting island influence extends into the DCM.

Non-synonymous to synonymous substitutions suggest that orthologs tend to keep their functions, while paralogs are a source of functional novelty

Orthologs separate after lineages split from each other and paralogs after gene duplications. Thus, orthologs are expected to remain more functionally coherent across lineages, while paralogs have been proposed as a source of new functions. Because protein functional divergence follows from non-synonymous substitutions, we performed an analysis based on the ratio of non-synonymous to synonymous substitutions (dN/dS), as proxy for functional divergence. We used five working definitions of orthology, including reciprocal best hits (RBH), among other definitions based on network analyses and clustering. The results showed that orthologs, by all definitions tested, had values of dN/dS noticeably lower than those of paralogs, suggesting that orthologs generally tend to be more functionally stable than paralogs. The differences in dN/dS ratios remained suggesting the functional stability of orthologs after eliminating gene comparisons with potential problems, such as genes with high codon usage biases, low coverage of either of the aligned sequences, or sequences with very high similarities. Separation by percent identity of the encoded proteins showed that the differences between the dN/dS ratios of orthologs and paralogs were more evident at high sequence identity, less so as identity dropped. The last results suggest that the differences between dN/dS ratios were partially related to differences in protein identity. However, they also suggested that paralogs undergo functional divergence relatively early after duplication. Our analyses indicate that choosing orthologs as probably functionally coherent remains the right approach in comparative genomics.

Microbial genomic trait evolution is dominated by frequent and rare pulsed evolution

On the macroevolutionary time scale, does trait evolution proceed gradually or by rapid bursts (pulses) separated by prolonged periods of stasis or slow evolution? Although studies have shown that pulsed evolution is prevalent in animals, our knowledge about the tempo and mode of evolution across the tree of life is very limited. This long-standing debate calls for a test in bacteria and archaea, the most ancient and diverse forms of life with unique population genetic properties. Using a likelihood-based framework, we show that pulsed evolution is not only present but also prevalent and predominant in microbial genomic trait evolution. We detected two distinct types of pulsed evolution (small frequent and large rare jumps) that are predicted by the punctuated equilibrium and quantum evolution theories. Our findings suggest that major bacterial lineages could have originated in quick bursts and that pulsed evolution is a common theme across the tree of life.

Pseudofinder: detection of pseudogenes in prokaryotic genomes

Prokaryotic genomes are usually densely packed with intact and functional genes. However, in certain contexts, such as after recent ecological shifts or extreme population bottlenecks, broken and nonfunctional gene fragments can quickly accumulate and form a substantial fraction of the genome. Identification of these broken genes, called pseudogenes, is a critical step for understanding the evolutionary forces acting upon, and the functional potential encoded within, prokaryotic genomes. Here, we present Pseudofinder, an open-source software dedicated to pseudogene identification and analysis in bacterial and archaeal genomes. We demonstrate that Pseudofinder’s multi-pronged, reference-based approach can detect a wide variety of pseudogenes, including those that are highly degraded and typically missed by gene-calling pipelines, as well newly formed pseudogenes containing only one or a few inactivating mutations. Additionally, Pseudofinder can detect genes that lack inactivating substitutions but experiencing relaxed selection. Implementation of Pseudofinder in annotation pipelines will allow more precise estimations of the functional potential of sequenced microbes, while also generating new hypotheses related to the evolutionary dynamics of bacterial and archaeal genomes.

Distinct Potentially Adaptive Accumulation of Truncation Mutations in Salmonella enterica serovar Typhi and Salmonella enterica serovar Paratyphi A

Gene inactivation through the accumulation of truncation (or premature stop codon) mutations is a common mode of evolution in bacteria. It is frequently believed to result from reductive evolutionary processes allowing purging of superfluous traits. However, several works have demonstrated that, similar to the occurrences of inactivating nonsynonymous (i.e., amino acid replacement) mutations under positive selection pressures, truncation mutations can also be adaptive where specific traits deleterious in particular environmental conditions need to be inactivated for survival. Here, we performed a comparative analysis of genome-wide accumulation of truncation mutations in Salmonella enterica serovar Typhi and Salmonella enterica serovar Paratyphi A. Considering the known convergent evolutionary trajectories in these two serovars, we expected a strong overlap of truncated genes in S. Typhi and S. Paratyphi A, emerging through either reductive or adaptive dynamics. However, we detected a distinct set of core truncated genes encoding different overrepresented functional clusters in each serovar. In 54% and 28% truncated genes in S. Typhi and S. Paratyphi A, respectively, inactivating mutations were acquired by only different subsets of isolates, instead of all isolates analyzed for that serovar. Importantly, 62% truncated genes (P < 0.001) in S. Typhi and S. Paratyphi A were also targeted by convergent amino acid mutations in different serovars, suggesting those genes to be under selection pressures. Our findings indicate significant presence of potentially adaptive truncation mutations in conjunction with the ones emerging due to reductive evolution. Further experimental and large-scale bioinformatic studies are necessary to better explore the impact of such adaptive footprints of truncation mutations in the evolution of bacterial virulence. IMPORTANCE Detecting the adaptive mutations leading to gene inactivation or loss of function is crucial for understanding their contribution in the evolution of bacterial virulence and antibiotic resistance. Such inactivating mutations, apart from being of nonsynonymous (i.e., amino acid replacement) nature, can also be truncation mutations, abruptly trimming the length of encoded proteins. Importantly, the notion of reductive evolutionary dynamics is primarily accepted toward the accumulation of truncation mutations. However, our case study on S. Typhi and S. Paratyphi A, two human-restricted systemically invasive pathogens exerting similar clinical manifestations, indicated that a significant proportion of truncation mutations emerge from positive selection pressures. The candidate genes from our study will enable directed functional assays for deciphering the adaptive role of truncation mutations in S. Typhi and S. Paratyphi A pathogenesis. Also, our genome-level analytical approach will pave the way to understand the contribution of truncation mutations in the adaptive evolution of other bacterial pathogens.

Characterisation of the symbionts in the Mediterranean fruit fly gut

Symbioses between bacteria and their insect hosts can range from loose associations through to obligate interdependence. While fundamental evolutionary insights have been gained from the in-depth study of obligate mutualisms, there is increasing interest in the evolutionary potential of flexible symbiotic associations between hosts and their gut microbiomes. Understanding relationships between microbes and hosts also offers the potential for exploitation for insect control. Here, we investigate the gut microbiome of a global agricultural pest, the Mediterranean fruit fly (Ceratitis capitata). We used 16S rRNA profiling to compare the gut microbiomes of laboratory and wild strains raised on different diets and from flies collected from various natural plant hosts. The results showed that medfly guts harbour a simple microbiome that is primarily determined by the larval diet. However, regardless of the laboratory diet or natural plant host on which flies were raised, Klebsiella spp. dominated medfly microbiomes and were resistant to removal by antibiotic treatment. We sequenced the genome of the dominant putative Klebsiella spp. (‘Medkleb’) isolated from the gut of the Toliman wild-type strain. Genome-wide ANI analysis placed Medkleb within the K. oxytoca / michiganensis group. Species level taxonomy for Medkleb was resolved using a mutli-locus phylogenetic approach - and molecular, sequence and phenotypic analyses all supported its identity as K. michiganensis. Medkleb has a genome size (5825435 bp) which is 1.6 standard deviations smaller than the mean genome size of free-living Klebsiella spp. Medkleb also lacks some genes involved in environmental sensing. Moreover, the Medkleb genome contains at least two recently acquired unique genomic islands as well as genes that encode pectinolytic enzymes capable of degrading plant cell walls. This may be advantageous given that the medfly diet includes unripe fruits containing high proportions of pectin. The results suggest that the medfly harbours a commensal gut bacterium that may have developed a mutualistic association with its host and provide nutritional benefits.

A high-quality reference genome for the fish pathogen Streptococcus iniae

Fish mortality caused by Streptococcus iniae is a major economic problem in aquaculture in warm and temperate regions globally. There is also risk of zoonotic infection by S. iniae through handling of contaminated fish. In this study, we present the complete genome sequence of S. iniae strain QMA0248, isolated from farmed barramundi in South Australia. The 2.12 Mb genome of S. iniae QMA0248 carries a 32 kb prophage, a 12 kb genomic island and 92 discrete insertion sequence (IS) elements. These include nine novel IS types that belong mostly to the IS3 family. Comparative and phylogenetic analysis between S. iniae QMA0248 and publicly available complete S. iniae genomes revealed discrepancies that are probably due to misassembly in the genomes of isolates ISET0901 and ISNO. Long-range PCR confirmed five rRNA loci in the PacBio assembly of QMA0248, and, unlike S. iniae 89353, no tandemly repeated rRNA loci in the consensus genome. However, we found sequence read evidence that the tandem rRNA repeat existed within a subpopulation of the original QMA0248 culture. Subsequent nanopore sequencing revealed that the tandem rRNA repeat was the most prevalent genotype, suggesting that there is selective pressure to maintain fewer rRNA copies under uncertain laboratory conditions. Our study not only highlights assembly problems in existing genomes, but provides a high-quality reference genome for S. iniae QMA0248, including manually curated mobile genetic elements, that will assist future S. iniae comparative genomic and evolutionary studies.

Comparative Genome Analysis of ‘Candidatus Phytoplasma luffae’ Reveals the Influential Roles of Potential Mobile Units in Phytoplasma Evolution

Phytoplasmas are insect-transmitted plant pathogens that cause substantial losses in agriculture. In addition to economic impact, phytoplasmas induce distinct disease symptoms in infected plants, thus attracting attention for research on molecular plant-microbe interactions and plant developmental processes. Due to the difficulty of establishing an axenic culture of these bacteria, culture-independent genome characterization is a crucial tool for phytoplasma research. However, phytoplasma genomes have strong nucleotide composition biases and are repetitive, which make it challenging to produce complete assemblies. In this study, we utilized Illumina and Oxford Nanopore sequencing technologies to obtain the complete genome sequence of ‘Candidatus Phytoplasma luffae’ strain NCHU2019 that is associated with witches’ broom disease of loofah (Luffa aegyptiaca) in Taiwan. The fully assembled circular chromosome is 769 kb in size and is the first representative genome sequence of group 16SrVIII phytoplasmas. Comparative analysis with other phytoplasmas revealed that NCHU2019 has a remarkably repetitive genome, possessing a pair of 75 kb repeats and at least 13 potential mobile units (PMUs) that account for ∼25% of its chromosome. This level of genome repetitiveness is exceptional for bacteria, particularly among obligate pathogens with reduced genomes. Our genus-level analysis of PMUs demonstrated that these phytoplasma-specific mobile genetic elements can be classified into three major types that differ in gene organization and phylogenetic distribution. Notably, PMU abundance explains nearly 80% of the variance in phytoplasma genome sizes, a finding that provides a quantitative estimate for the importance of PMUs in phytoplasma genome variability. Finally, our investigation found that in addition to horizontal gene transfer, PMUs also contribute to intra-genomic duplications of effector genes, which may provide redundancy for subfunctionalization or neofunctionalization. Taken together, this work improves the taxon sampling for phytoplasma genome research and provides novel information regarding the roles of mobile genetic elements in phytoplasma evolution.

A Peroxide-Responding sRNA Evolved from a Peroxidase mRNA

Small RNAs (sRNAs) are important gene regulators in bacteria, but it is unclear how new sRNAs originate and become part of regulatory networks that coordinate bacterial response to environmental stimuli. Using a covariance modeling-based approach, we analyzed the presence of hundreds of sRNAs in more than a thousand genomes across Enterobacterales, a bacterial order with a confluence of factors that allows robust genome-scale sRNA analyses: several well-studied organisms with fairly conserved genome structures, an established phylogeny, and substantial nucleotide diversity within a narrow evolutionary space. We discovered that a majority of sRNAs arose recently, and uncovered protein-coding genes as a potential source from which new sRNAs arise. A detailed investigation of the emergence of OxyS, a peroxide-responding sRNA, revealed that it evolved from a fragment of a peroxidase messenger RNA. Importantly, although it replaced the ancestral peroxidase, OxyS continues to be part of the ancestral peroxide-response regulon, indicating that an sRNA that arises from a protein-coding gene would inherently be part of the parental protein's regulatory network. This new insight provides a fresh framework for understanding sRNA origin and regulatory integration in bacteria.

Brenneria goodwinii growth in vitro is improved by competitive interactions with other bacterial species associated with Acute Oak Decline

•

Mutually competitive interactions prevail in pairwise cultures of AOD bacteria.

•

In vitro growth of brenneria goodwinii is improved by bacterial competition.

•

Co-culturing of AOD bacteria indicates evolving improved fitness of B. goodwinii.

•

B. goodwinii and R. victoriana can collaborate in vitro to outcompete G. quercinecans.

Brenneria goodwinii, Rahnella victoriana and Gibbsiella quercinecans are three bacterial species frequently isolated together from oak displaying symptoms of Acute Oak Decline (AOD), which include weeping patches on trunks. All three bacterial species play a role in lesion formation in the current episode of AOD in Britain, although B. goodwinii is the most dominant. The ongoing research into stem lesion formation characteristic of this polybacterial syndrome has been focussed primarily on the pathogenicity, identification and taxonomy of these bacteria. As all three species were newly classified within the past ten years, there are many unanswered questions regarding their ecology and interactions with each other. To determine the effect of bacterial interactions on fitness in vitro, we examined pairwise (diculture) and multispecies (triculture) interactions between B. goodwinii, R. victoriana and G. quercinecans in oak leaf media microcosms. Additionally, the effect of co-culturing on the evolution of these species was determined and the evolved B. goodwinii strains were examined further by whole genome sequencing. Our results indicate that B. goodwinii thrived in monoculture with significantly higher viable cell counts than the other two species. Additionally, B. goodwinii performed well in pairwise culture with mutually competitive interactions observed between B. goodwinii and R. victoriana, and between B. goodwinii and G. quercinecans. In the multispecies triculture, B. goodwinii and R. victoriana appeared to exhibit co-ordinated behaviour to outcompete G. quercinecans. After four weeks B. goodwinii grown in co-culture with the other two species developed greater evolved fitness than the strain grown in monoculture as reflected by the increased viable cell counts. The competitive interactions taking place between the threes species indicated evolving improved fitness of B. goodwinii in vitro, that gave it a growth advantage over both R. victoriana and G. quercinecans which showed no significant changes in fitness. Overall, B. goodwinii gains greater benefit in terms of fitness from in vitro competitive interaction with the other two species.

In Vitro Probiotic Characteristics and Whole Genome Sequence Analysis of Lactobacillus Strains Isolated from Cattle-Yak Milk

Cattle-yak milk is an important raw material and an indispensable source of high-quality food for local farmers and herdsmen to produce ghee, milk residue, yogurt, and other dairy products. In this study, Lactobacillus strains were isolated from cattle-yak milk for potential probiotic candidates using a series of in vitro tests, including probiotic characterization and safety evaluation (antibiotic susceptibility and hemolytic ability). The results found that the Lactobacillus rhamnosus CY12 strain showed a high survival rate in bile salts, under acid conditions, and in the gastrointestinal juice environment, as well as showing high antimicrobial activity and adhesive potential. The safety evaluation showed that all strains were considered non-hemolytic. In addition, the whole-genome sequencing indicated that the strain CY12 spanned 2,506,167 bp, with an average length of 881 bp; the GC content in the gene region (%) was 47.35, contained 1347 protein-coding sequences, and accounted for 85.72% of the genome. The genome annotation showed that genes mainly focused on the immune system process, metabolic process, carbohydrate utilization, carbon metabolism, galactose metabolism, and biological adhesion, etc. This study revealed that the Lactobacillus rhamnosus CY12 strain might be an excellent potential probiotic in the development of feed additives for animals and has the ability to promote health.

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.

The Evolution of Interdependence in a Four-Way Mealybug Symbiosis

Mealybugs are insects that maintain intracellular bacterial symbionts to supplement their nutrient-poor plant sap diets. Some mealybugs have a single betaproteobacterial endosymbiont, a Candidatus Tremblaya species (hereafter Tremblaya) that alone provides the insect with its required nutrients. Other mealybugs have two nutritional endosymbionts that together provision these same nutrients, where Tremblaya has gained a gammaproteobacterial partner that resides in its cytoplasm. Previous work had established that Pseudococcus longispinus mealybugs maintain not one but two species of gammaproteobacterial endosymbionts along with Tremblaya. Preliminary genomic analyses suggested that these two gammaproteobacterial endosymbionts have large genomes with features consistent with a relatively recent origin as insect endosymbionts, but the patterns of genomic complementarity between members of the symbiosis and their relative cellular locations were unknown. Here, using long-read sequencing and various types of microscopy, we show that the two gammaproteobacterial symbionts of P. longispinus are mixed together within Tremblaya cells, and that their genomes are somewhat reduced in size compared with their closest nonendosymbiotic relatives. Both gammaproteobacterial genomes contain thousands of pseudogenes, consistent with a relatively recent shift from a free-living to an endosymbiotic lifestyle. Biosynthetic pathways of key metabolites are partitioned in complex interdependent patterns among the two gammaproteobacterial genomes, the Tremblaya genome, and horizontally acquired bacterial genes that are encoded on the mealybug nuclear genome. Although these two gammaproteobacterial endosymbionts have been acquired recently in evolutionary time, they have already evolved codependencies with each other, Tremblaya, and their insect host.

The Genomics and Cell Biology of Host-Beneficial Intracellular Infections

Microbes gain access to eukaryotic cells as food for bacteria-grazing protists, for host protection by microbe-killing immune cells, or for microbial benefit when pathogens enter host cells to replicate. But microbes can also gain access to a host cell and become an important-often required-beneficial partner. The oldest beneficial microbial infections are the ancient eukaryotic organelles now called the mitochondrion and plastid. But numerous other host-beneficial intracellular infections occur throughout eukaryotes. Here I review the genomics and cell biology of these interactions with a focus on intracellular bacteria. The genomes of host-beneficial intracellular bacteria have features that span a previously unfilled gap between pathogens and organelles. Host cell adaptations to allow the intracellular persistence of beneficial bacteria are found along with evidence for the microbial manipulation of host cells, but the cellular mechanisms of beneficial bacterial infections are not well understood. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Comparative genomics reveals environmental adaptation differences between Cronobacter species

The genus Cronobacter is an opportunistic food-borne pathogen which is able to adapt to diverse environments and shows considerable genetic diversity. Genomic analysis can be used to reveal the variation across the genus with respect to virulence, drug resistance and factors involved in horizontal gene transfer mechanisms, such as integrons, conjugative plasmids, and recombinases. In this study, whole-genome comparative analysis of 27 Cronobacter genomes (12 existing and 15 newly assembled genomes) was performed. A total of 110,010 protein-coding genes were grouped into 11,262 clusters, including 2637 core genes, 4814 strain-specific genes and 3811 dispensable genes. Clusters of Orthologous Groups (COG) analysis indicated that 97.35% of the core genes were for substrate transport and metabolism, and the antibiotic resistance genetic determinants were classified into 136 antibiotic resistance ontologies (AROs). A total of 80 genomic islands (GIs) were identified which contained several type VI secretion system gene clusters, and these were likely to have been acquired by horizontal gene transfer. This study has generated a comprehensive characterization of the genus Cronobacter, leading to a better understanding of the mechanisms and ecological functions among the genome features, speciation, and environmental adaptation strategies.

Genome and Pangenome Analysis of Lactobacillus hilgardii FLUB-A New Strain Isolated from Mead

The production of mead holds great value for the Polish liquor industry, which is why the bacterium that spoils mead has become an object of concern and scientific interest. This article describes, for the first time, Lactobacillus hilgardii FLUB newly isolated from mead, as a mead spoilage bacteria. Whole genome sequencing of L. hilgardii FLUB revealed a 3 Mbp chromosome and five plasmids, which is the largest reported genome of this species. An extensive phylogenetic analysis and digital DNA-DNA hybridization confirmed the membership of the strain in the L. hilgardii species. The genome of L. hilgardii FLUB encodes 3043 genes, 2871 of which are protein coding sequences, 79 code for RNA, and 93 are pseudogenes. L. hilgardii FLUB possesses three clustered regularly interspaced short palindromic repeats (CRISPR), eight genomic islands (44,155 bp to 6345 bp), and three (two intact and one incomplete) prophage regions. For the first time, the characteristics of the genome of this species were described and a pangenomic analysis was performed. The concept of the pangenome was used not only to establish the genetic repertoire of this species, but primarily to highlight the unique characteristics of L. hilgardii FLUB. The core of the genome of L. hilgardii is centered around genes related to the storage and processing of genetic information, as well as to carbohydrate and amino acid metabolism. Strains with such a genetic constitution can effectively adapt to environmental changes. L. hilgardii FLUB is distinguished by an extensive cluster of metabolic genes, arsenic detoxification genes, and unique surface layer proteins. Variants of MRS broth with ethanol (10-20%), glucose (2-25%), and fructose (2-24%) were prepared to test the strain's growth preferences using Bioscreen C and the PYTHON script. L. hilgardii FLUB was found to be more resistant than a reference strain to high concentrations of alcohol (18%) and sugars (25%). It exhibited greater preference for fructose than glucose, which suggests it has a fructophilic nature. Comparative genomic analysis supported by experimental research imitating the conditions of alcoholic beverages confirmed the niche specialization of L. hilgardii FLUB to the mead environment.

Gene loss through pseudogenization contributes to the ecological diversification of a generalist Roseobacter lineage

Ecologically relevant genes generally show patchy distributions among related bacterial genomes. This is commonly attributed to lateral gene transfer, whereas the opposite mechanism-gene loss-has rarely been explored. Pseudogenization is a major mechanism underlying gene loss, and pseudogenes are best characterized by comparing closely related genomes because of their short life spans. To explore the role of pseudogenization in microbial ecological diversification, we apply rigorous methods to characterize pseudogenes in the 279 newly sequenced Ruegeria isolates of the globally abundant Roseobacter group collected from two typical coastal habitats in Hong Kong, the coral Platygyra acuta and the macroalga Sargassum hemiphyllum. Pseudogenes contribute to ~16% of the accessory genomes of these strains. Ancestral state reconstruction reveals that many pseudogenization events are correlated with ancestral niche shifts. Specifically, genes related to resource scavenging and energy acquisition were often pseudogenized when roseobacters inhabiting carbon-limited and energy-poor coral skeleton switched to other resource-richer niches. For roseobacters inhabiting the macroalgal niches, genes for nitrogen regulation and carbohydrate utilization were important but became dispensable upon shift to coral skeleton where nitrate is abundant but carbohydrates are less available. Whereas low-energy-demanding secondary transporters are more favorable in coral skeleton, ATP-driven primary transporters are preferentially kept in the energy-replete macroalgal niches. Moreover, a large proportion of these families mediate organismal interactions, suggesting their rapid losses by pseudogenization as a potential response to host and niche shift. These findings illustrate an important role of pseudogenization in shaping genome content and driving ecological diversification of marine roseobacters.

Molecular insights into the genome dynamics and interactions between core and acquired genomes of Vibrio cholerae

Bacterial species are hosts to horizontally acquired mobile genetic elements (MGEs), which encode virulence, toxin, antimicrobial resistance, and other metabolic functions. The bipartite genome of Vibrio cholerae harbors sporadic and conserved MGEs that contribute in the disease development and survival of the pathogens. For a comprehensive understanding of dynamics of MGEs in the bacterial genome, we engineered the genome of V. cholerae and examined in vitro and in vivo stability of genomic islands (GIs), integrative conjugative elements (ICEs), and prophages. Recombinant vectors carrying the integration module of these GIs, ICE and CTXΦ, helped us to understand the efficiency of integrations of MGEs in the V. cholerae chromosome. We have deleted more than 250 acquired genes from 6 different loci in the V. cholerae chromosome and showed contribution of CTX prophage in the essentiality of SOS response master regulator LexA, which is otherwise not essential for viability in other bacteria, including Escherichia coli In addition, we observed that the core genome-encoded RecA helps CTXΦ to bypass V. cholerae immunity and allow it to replicate in the host bacterium in the presence of similar prophage in the chromosome. Finally, our proteomics analysis reveals the importance of MGEs in modulating the levels of cellular proteome. This study engineered the genome of V. cholerae to remove all of the GIs, ICEs, and prophages and revealed important interactions between core and acquired genomes.

Selection-Driven Gene Inactivation in Salmonella

Bacterial genes are sometimes found to be inactivated by mutation. This inactivation may be observable simply because selection for function is intermittent or too weak to eliminate inactive alleles quickly. Here, I investigate cases in Salmonella enterica where inactivation is instead positively selected. These are identified by a rate of introduction of premature stop codons to a gene that is higher than expected under selective neutrality, as assessed by comparison to the rate of synonymous changes. I identify 84 genes that meet this criterion at a 10% false discovery rate. Many of these genes are involved in virulence, motility and chemotaxis, biofilm formation, and resistance to antibiotics or other toxic substances. It is hypothesized that most of these genes are subject to an ongoing process in which inactivation is favored under rare conditions, but the inactivated allele is deleterious under most other conditions and is subsequently driven to extinction by purifying selection.

Emerging Priorities for Microbiome Research

Microbiome research has increased dramatically in recent years, driven by advances in technology and significant reductions in the cost of analysis. Such research has unlocked a wealth of data, which has yielded tremendous insight into the nature of the microbial communities, including their interactions and effects, both within a host and in an external environment as part of an ecological community. Understanding the role of microbiota, including their dynamic interactions with their hosts and other microbes, can enable the engineering of new diagnostic techniques and interventional strategies that can be used in a diverse spectrum of fields, spanning from ecology and agriculture to medicine and from forensics to exobiology. From June 19-23 in 2017, the NIH and NSF jointly held an Innovation Lab on Quantitative Approaches to Biomedical Data Science Challenges in our Understanding of the Microbiome. This review is inspired by some of the topics that arose as priority areas from this unique, interactive workshop. The goal of this review is to summarize the Innovation Lab's findings by introducing the reader to emerging challenges, exciting potential, and current directions in microbiome research. The review is broken into five key topic areas: (1) interactions between microbes and the human body, (2) evolution and ecology of microbes, including the role played by the environment and microbe-microbe interactions, (3) analytical and mathematical methods currently used in microbiome research, (4) leveraging knowledge of microbial composition and interactions to develop engineering solutions, and (5) interventional approaches and engineered microbiota that may be enabled by selectively altering microbial composition. As such, this review seeks to arm the reader with a broad understanding of the priorities and challenges in microbiome research today and provide inspiration for future investigation and multi-disciplinary collaboration.

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

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

Genomic analysis of bacteria in the Acute Oak Decline pathobiome

The UK’s native oak is under serious threat from Acute Oak Decline (AOD). Stem tissue necrosis is a primary symptom of AOD and several bacteria are associated with necrotic lesions. Two members of the lesion pathobiome, Brenneria goodwinii and Gibbsiella quercinecans, have been identified as causative agents of tissue necrosis. However, additional bacteria including Lonsdalea britannica and Rahnella species have been detected in the lesion microbiome, but their role in tissue degradation is unclear. Consequently, information on potential genome-encoded mechanisms for tissue necrosis is critical to understand the role and mechanisms used by bacterial members of the lesion pathobiome in the aetiology of AOD. Here, the whole genomes of bacteria isolated from AOD-affected trees were sequenced, annotated and compared against canonical bacterial phytopathogens and non-pathogenic symbionts. Using orthologous gene inference methods, shared virulence genes that retain the same function were identified. Furthermore, functional annotation of phytopathogenic virulence genes demonstrated that all studied members of the AOD lesion microbiota possessed genes associated with phytopathogens. However, the genome of B. goodwinii was the most characteristic of a necrogenic phytopathogen, corroborating previous pathological and metatranscriptomic studies that implicate it as the key causal agent of AOD lesions. Furthermore, we investigated the genome sequences of other AOD lesion microbiota to understand the potential ability of microbes to cause disease or contribute to pathogenic potential of organisms isolated from this complex pathobiome. The role of these members remains uncertain but some such as G. quercinecans may contribute to tissue necrosis through the release of necrotizing enzymes and may help more dangerous pathogens activate and realize their pathogenic potential or they may contribute as secondary/opportunistic pathogens with the potential to act as accessory species for B. goodwinii. We demonstrate that in combination with ecological data, whole genome sequencing provides key insights into the pathogenic potential of bacterial species whether they be phytopathogens, part-contributors or stimulators of the pathobiome.

Changes detected in the genome sequences of Escherichia coli, Listeria monocytogenes, Vibrio parahaemolyticus, and Salmonella enterica after serial subculturing

Whole genome sequencing (WGS) is rapidly replacing other molecular techniques for identifying and subtyping bacterial isolates. The resolution or discrimination offered by WGS is significantly higher than that offered by other molecular techniques, and WGS readily allows infrequent differences that occur between 2 closely related strains to be found. In this investigation, WGS was used to identify the changes that occurred in the genomes of 13 strains of bacterial foodborne pathogens after 100 serial subcultures. Pure cultures of Shiga-toxin-producing Escherichia coli, Salmonella enterica, Listeria monocytogenes, and Vibrio parahaemolyticus were subcultured daily for 100 successive days. The 1st and 100th subcultures were whole-genome sequenced using short-read sequencing. Single nucleotide polymorphisms (SNPs) were identified between the 1st and final culture using 2 different approaches, and multilocus sequence typing of the whole genome was also performed to detect any changes at the allelic level. The number of observed genomic changes varied by strain, species, and the SNP caller used. This study provides insight into the genomic variation that can be detected using next-generation sequencing and analysis methods after repeated subculturing of 4 important bacterial pathogens.

Seasons of change: Mechanisms of genome evolution in human fungal pathogens

Fungi are a diverse kingdom of organisms capable of thriving in various niches across the world including those in close association with multicellular eukaryotes. Fungal pathogens that contribute to human disease reside both within the host as commensal organisms of the microbiota and the environment. Their niche of origin dictates how infection initiates but also places specific selective pressures on the fungal pathogen that contributes to its genome organization and genetic repertoire. Recent efforts to catalogue genomic variation among major human fungal pathogens have unveiled evolutionary themes that shape the fungal genome. Mechanisms ranging from large scale changes such as aneuploidy and ploidy cycling as well as more targeted mutations like base substitutions and gene copy number variations contribute to the evolution of these species, which are often under multiple competing selective pressures with their host, environment, and other microbes. Here, we provide an overview of the major selective pressures and mechanisms acting to evolve the genome of clinically important fungal pathogens of humans.

Evolution of host support for two ancient bacterial symbionts with differentially degraded genomes in a leafhopper host

Plant sap-feeding insects (Hemiptera) rely on bacterial symbionts for nutrition absent in their diets. These bacteria experience extreme genome reduction and require genetic resources from their hosts, particularly for basic cellular processes other than nutrition synthesis. The host-derived mechanisms that complete these processes have remained poorly understood. It is also unclear how hosts meet the distinct needs of multiple bacterial partners with differentially degraded genomes. To address these questions, we investigated the cell-specific gene-expression patterns in the symbiotic organs of the aster leafhopper (ALF), Macrosteles quadrilineatus (Cicadellidae). ALF harbors two intracellular symbionts that have two of the smallest known bacterial genomes: Nasuia (112 kb) and Sulcia (190 kb). Symbionts are segregated into distinct host cell types (bacteriocytes) and vary widely in their basic cellular capabilities. ALF differentially expresses thousands of genes between the bacteriocyte types to meet the functional needs of each symbiont, including the provisioning of metabolites and support of cellular processes. For example, the host highly expresses genes in the bacteriocytes that likely complement gene losses in nucleic acid synthesis, DNA repair mechanisms, transcription, and translation. Such genes are required to function in the bacterial cytosol. Many host genes comprising these support mechanisms are derived from the evolution of novel functional traits via horizontally transferred genes, reassigned mitochondrial support genes, and gene duplications with bacteriocyte-specific expression. Comparison across other hemipteran lineages reveals that hosts generally support the incomplete symbiont cellular processes, but the origins of these support mechanisms are generally specific to the host-symbiont system.

Horizontal Gene Acquisitions Contributed to Genome Expansion in Insect-Symbiotic Spiroplasma clarkii

Genome reduction is a recurring theme of symbiont evolution. The genus Spiroplasma contains species that are mostly facultative insect symbionts. The typical genome sizes of those species within the Apis clade were estimated to be ∼1.0–1.4 Mb. Intriguingly, Spiroplasma clarkii was found to have a genome size that is >30% larger than the median of other species within the same clade. To investigate the molecular evolution events that led to the genome expansion of this bacterium, we determined its complete genome sequence and inferred the evolutionary origin of each protein-coding gene based on the phylogenetic distribution of homologs. Among the 1,346 annotated protein-coding genes, 641 were originated from within the Apis clade while 233 were putatively acquired from outside of the clade (including 91 high-confidence candidates). Additionally, 472 were specific to S. clarkii without homologs in the current database (i.e., the origins remained unknown). The acquisition of protein-coding genes, rather than mobile genetic elements, appeared to be a major contributing factor of genome expansion. Notably, >50% of the high-confidence acquired genes are related to carbohydrate transport and metabolism, suggesting that these acquired genes contributed to the expansion of both genome size and metabolic capability. The findings of this work provided an interesting case against the general evolutionary trend observed among symbiotic bacteria and further demonstrated the flexibility of Spiroplasma genomes. For future studies, investigation on the functional integration of these acquired genes, as well as the inference of their contribution to fitness could improve our knowledge of symbiont evolution.

Parallel Evolution of Genome Streamlining and Cellular Bioenergetics across the Marine Radiation of a Bacterial Phylum

Understanding long-term patterns of microbial evolution is critical to advancing our knowledge of past and present role microbial life in driving global biogeochemical cycles. Historically, it has been challenging to study the evolution of environmental microbes due to difficulties in obtaining genome sequences from lineages that could not be cultivated, but recent advances in metagenomics and single-cell genomics have begun to obviate many of these hurdles. Here we present an evolutionary genomic analysis of the Marinimicrobia, a diverse bacterial group that is abundant in the global ocean. We demonstrate that distantly related Marinimicrobia species that reside in similar habitats have converged to assume similar genome architectures and cellular bioenergetics, suggesting that common factors shape the evolution of a broad array of marine lineages. These findings broaden our understanding of the evolutionary forces that have given rise to microbial life in the contemporary ocean.

Diverse bacterial and archaeal lineages drive biogeochemical cycles in the global ocean, but the evolutionary processes that have shaped their genomic properties and physiological capabilities remain obscure. Here we track the genome evolution of the globally abundant marine bacterial phylum Marinimicrobia across its diversification into modern marine environments and demonstrate that extant lineages are partitioned between epipelagic and mesopelagic habitats. Moreover, we show that these habitat preferences are associated with fundamental differences in genomic organization, cellular bioenergetics, and metabolic modalities. Multiple lineages present in epipelagic niches independently acquired genes necessary for phototrophy and environmental stress mitigation, and their genomes convergently evolved key features associated with genome streamlining. In contrast, lineages residing in mesopelagic waters independently acquired nitrate respiratory machinery and a variety of cytochromes, consistent with the use of alternative terminal electron acceptors in oxygen minimum zones (OMZs). Further, while epipelagic clades have retained an ancestral Na⁺-pumping respiratory complex, mesopelagic lineages have largely replaced this complex with canonical H⁺-pumping respiratory complex I, potentially due to the increased efficiency of the latter together with the presence of the more energy-limiting environments deep in the ocean’s interior. These parallel evolutionary trends indicate that key features of genomic streamlining and cellular bioenergetics have occurred repeatedly and congruently in disparate clades and underscore the importance of environmental conditions and nutrient dynamics in driving the evolution of diverse bacterioplankton lineages in similar ways throughout the global ocean.

Convergent Evolution among Ruminant-Pathogenic Mycoplasma Involved Extensive Gene Content Changes

Convergent evolution, a process by which organisms evolved independently to have similar traits, provides opportunities to understand adaptation. The bacterial genus Mycoplasma contains multiple species that evolved independently to become ruminant pathogens, which represents an interesting study system for investigating the process. In this work, we determined the genome sequences of 11 Entomoplasma/Mesoplasma species. This new data set, together with the other available Mollicutes genomes, provided comprehensive taxon sampling for inferring the gene content evolution that led to the emergence of Mycoplasma Mycoides cluster. Our results indicated that the most recent common ancestor (MRCA) of the Mycoides-Entomoplasmataceae clade lost ∼15% of the core genes when it diverged from the Spiroplasma Apis clade. After this initial wave of genome reduction, relatively few gene gains or losses were inferred until the emergence of the Mycoides cluster. Compared with those Entomoplasmataceae lineages that maintained the association with insects, the MRCA of the Mycoides cluster experienced a second wave of gene losses, as well as acquiring >100 novel genes through horizontal gene transfer. These gene acquisitions involved many with the Mycoplasma Hominis/Pneumoniae lineages as the putative donors, suggesting that gene exchanges among these vertebrate symbionts with distinct phylogenetic affiliations may be important in the emergence of the Mycoides cluster. These findings demonstrated that the gene content of bacterial genomes could be exceedingly dynamic, even for those symbionts with highly reduced genomes. Moreover, the emergence of novel pathogens may involve extensive remodeling of gene content, rather than acquisition of few virulence genes.

Horizontal Acquisition and Transcriptional Integration of Novel Genes in Mosquito-Associated Spiroplasma

Genetic differentiation among symbiotic bacteria is important in shaping biodiversity. The genus Spiroplasma contains species occupying diverse niches and is a model system for symbiont evolution. Previous studies have established that two mosquito-associated species have diverged extensively in their carbohydrate metabolism genes despite having a close phylogenetic relationship. Notably, although the commensal Spiroplasma diminutum lacks identifiable pathogenicity factors, the pathogenic Spiroplasma taiwanense was found to have acquired a virulence factor glpO and its associated genes through horizontal transfer. However, it is unclear if these acquired genes have been integrated into the regulatory network. In this study, we inferred the gene content evolution in these bacteria, as well as examined their transcriptomes in response to glucose availability. The results indicated that both species have many more gene acquisitions from the Mycoides-Entomoplasmataceae clade, which contains several important pathogens of ruminants, than previously thought. Moreover, several acquired genes have higher expression levels than the vertically inherited homologs, indicating possible functional replacement. Finally, the virulence factor and its functionally linked genes in S. taiwanense were up-regulated in response to glucose starvation, suggesting that these acquired genes are under expression regulation and the pathogenicity may be a stress response. In summary, although differential gene losses are a major process for symbiont divergence, gene gains are critical in counteracting genome degradation and driving diversification among facultative symbionts.

Winding paths to simplicity: genome evolution in facultative insect symbionts

Symbiosis between organisms is an important driving force in evolution. Among the diverse relationships described, extensive progress has been made in insect–bacteria symbiosis, which improved our understanding of the genome evolution in host-associated bacteria. Particularly, investigations on several obligate mutualists have pushed the limits of what we know about the minimal genomes for sustaining cellular life. To bridge the gap between those obligate symbionts with extremely reduced genomes and their non-host-restricted ancestors, this review focuses on the recent progress in genome characterization of facultative insect symbionts. Notable cases representing various types and stages of host associations, including those from multiple genera in the family Enterobacteriaceae (class Gammaproteobacteria), Wolbachia (Alphaproteobacteria) and Spiroplasma (Mollicutes), are discussed. Although several general patterns of genome reduction associated with the adoption of symbiotic relationships could be identified, extensive variation was found among these facultative symbionts. These findings are incorporated into the established conceptual frameworks to develop a more detailed evolutionary model for the discussion of possible trajectories. In summary, transitions from facultative to obligate symbiosis do not appear to be a universal one-way street; switches between hosts and lifestyles (e.g. commensalism, parasitism or mutualism) occur frequently and could be facilitated by horizontal gene transfer.

This review synthesizes the recent progress in genome characterization of insect-symbiotic bacteria, the emphases include (i) patterns of genome organization, (ii) evolutionary models and trajectories, and (iii) comparisons between facultative and obligate symbionts.

HMGA1-pseudogenes and cancer

Pseudogenes are DNA sequences with high homology to the corresponding functional gene, but, because of the accumulation of various mutations, they have lost their initial functions to code for proteins. Consequently, pseudogenes have been considered until few years ago dysfunctional relatives of the corresponding ancestral genes, and then useless in the course of genome evolution. However, several studies have recently established that pseudogenes are owners of key biological functions. Indeed, some pseudogenes control the expression of functional genes by competitively binding to the miRNAs, some of them generate small interference RNAs to negatively modulate the expression of functional genes, and some of them even encode functional mutated proteins. Here, we concentrate our attention on the pseudogenes of the HMGA1 gene, that codes for the HMGA1a and HMGA1b proteins having a critical role in development and cancer progression. In this review, we analyze the family of HMGA1 pseudogenes through three aspects: classification, characterization, and their possible function and involvement in cancer.

Kullback Leibler divergence in complete bacterial and phage genomes

The amino acid content of the proteins encoded by a genome may predict the coding potential of that genome and may reflect lifestyle restrictions of the organism. Here, we calculated the Kullback–Leibler divergence from the mean amino acid content as a metric to compare the amino acid composition for a large set of bacterial and phage genome sequences. Using these data, we demonstrate that (i) there is a significant difference between amino acid utilization in different phylogenetic groups of bacteria and phages; (ii) many of the bacteria with the most skewed amino acid utilization profiles, or the bacteria that host phages with the most skewed profiles, are endosymbionts or parasites; (iii) the skews in the distribution are not restricted to certain metabolic processes but are common across all bacterial genomic subsystems; (iv) amino acid utilization profiles strongly correlate with GC content in bacterial genomes but very weakly correlate with the G+C percent in phage genomes. These findings might be exploited to distinguish coding from non-coding sequences in large data sets, such as metagenomic sequence libraries, to help in prioritizing subsequent analyses.

Genome of Ca. Pandoraea novymonadis, an Endosymbiotic Bacterium of the Trypanosomatid Novymonas esmeraldas

We have sequenced, annotated, and analyzed the genome of Ca. Pandoraea novymonadis, a recently described bacterial endosymbiont of the trypanosomatid Novymonas esmeraldas. When compared with genomes of its free-living relatives, it has all the hallmarks of the endosymbionts’ genomes, such as significantly reduced size, extensive gene loss, low GC content, numerous gene rearrangements, and low codon usage bias. In addition, Ca. P. novymonadis lacks mobile elements, has a strikingly low number of pseudogenes, and almost all genes are single copied. This suggests that it already passed the intensive period of host adaptation, which still can be observed in the genome of Polynucleobacter necessarius, a certainly recent endosymbiont. Phylogenetically, Ca. P. novymonadis is more related to P. necessarius, an intracytoplasmic bacterium of free-living ciliates, than to Ca. Kinetoplastibacterium spp., the only other known endosymbionts of trypanosomatid flagellates. As judged by the extent of the overall genome reduction and the loss of particular metabolic abilities correlating with the increasing dependence of the symbiont on its host, Ca. P. novymonadis occupies an intermediate position P. necessarius and Ca. Kinetoplastibacterium spp. We conclude that the relationships between Ca. P. novymonadis and N. esmeraldas are well-established, although not as fine-tuned as in the case of Strigomonadinae and their endosymbionts.

Microgravity validation of a novel system for RNA isolation and multiplex quantitative real time PCR analysis of gene expression on the International Space Station

The International Space Station (ISS) National Laboratory is dedicated to studying the effects of space on life and physical systems, and to developing new science and technologies for space exploration. A key aspect of achieving these goals is to operate the ISS National Lab more like an Earth-based laboratory, conducting complex end-to-end experimentation, not limited to simple microgravity exposure. Towards that end NASA developed a novel suite of molecular biology laboratory tools, reagents, and methods, named WetLab-2, uniquely designed to operate in microgravity, and to process biological samples for real-time gene expression analysis on-orbit. This includes a novel fluidic RNA Sample Preparation Module and fluid transfer devices, all-in-one lyophilized PCR assays, centrifuge, and a real-time PCR thermal cycler. Here we describe the results from the WetLab-2 validation experiments conducted in microgravity during ISS increment 47/SPX-8. Specifically, quantitative PCR was performed on a concentration series of DNA calibration standards, and Reverse Transcriptase-quantitative PCR was conducted on RNA extracted and purified on-orbit from frozen Escherichia coli and mouse liver tissue. Cycle threshold (Ct) values and PCR efficiencies obtained on-orbit from DNA standards were similar to Earth (1 g) controls. Also, on-orbit multiplex analysis of gene expression from bacterial cells and mammalian tissue RNA samples was successfully conducted in about 3 h, with data transmitted within 2 h of experiment completion. Thermal cycling in microgravity resulted in the trapping of gas bubbles inside septa cap assay tubes, causing small but measurable increases in Ct curve noise and variability. Bubble formation was successfully suppressed in a rapid follow-up on-orbit experiment using standard caps to pressurize PCR tubes and reduce gas release during heating cycles. The WetLab-2 facility now provides a novel operational on-orbit research capability for molecular biology and demonstrates the feasibility of more complex wet bench experiments in the ISS National Lab environment.

Pseudogene PDIA3P1 promotes cell proliferation, migration and invasion, and suppresses apoptosis in hepatocellular carcinoma by regulating the p53 pathway

Pseudogenes are a subclass of long non-coding (lnc) RNAs that arose from protein-coding genes, but have lost the ability to produce proteins. Pseudogenes play an important role in the pathogenesis of various tumors; however, the role of pseudogenes in hepatocellular carcinoma (HCC) is poorly understood. In this study, we investigated a novel pseudogene, PDIA3P1, which was upregulated in HCC tissues compared with paired normal adjacent tissues. The expression of PDIA3P1 was significantly correlated with tumor size, metastasis, TNM stage, and overall stage. Knockdown of PDIA3P1decreased proliferation, migration, and invasion of HCC cells. PDIA3P1 knockdown also promoted cell apoptosis and inhibited tumor growth in vivo. We performed a GeneChip assay to investigate the underlying mechanism of PDIA3P1 action on biological function, and our results suggested that PDIA3P1 may promote proliferation and inhibit apoptosis of liver cancer cells by inhibiting the p53 pathway. Thus, our data suggest that PDIA3P1 acts as an oncogene in HCC and could be a potential prognostic marker and therapeutic target for HCC.

Environmental drivers of a microbial genomic transition zone in the ocean's interior

The core properties of microbial genomes, including GC content and genome size, are known to vary widely among different bacteria and archaea ^1,2 . Several hypotheses have been proposed to explain this genomic variability, but the fundamental drivers that shape bacterial and archaeal genomic properties remain uncertain ^3-7 . Here, we report the existence of a sharp genomic transition zone below the photic zone, where bacterial and archaeal genomes and proteomes undergo a community-wide punctuated shift. Across a narrow range of increasing depth of just tens of metres, diverse microbial clades trend towards larger genome size, higher genomic GC content, and proteins with higher nitrogen but lower carbon content. These community-wide changes in genome features appear to be driven by gradients in the surrounding environmental energy and nutrient fields. Collectively, our data support hypotheses invoking nutrient limitation as a central driver in the evolution of core bacterial and archaeal genomic and proteomic properties.

LifeStyle-Specific-Islands (LiSSI): Integrated Bioinformatics Platform for Genomic Island Analysis

Distinct bacteria are able to cope with highly diverse lifestyles; for instance, they can be free living or host-associated. Thus, these organisms must possess a large and varied genomic arsenal to withstand different environmental conditions. To facilitate the identification of genomic features that might influence bacterial adaptation to a specific niche, we introduce LifeStyle-Specific-Islands (LiSSI). LiSSI combines evolutionary sequence analysis with statistical learning (Random Forest with feature selection, model tuning and robustness analysis). In summary, our strategy aims to identify conserved consecutive homology sequences (islands) in genomes and to identify the most discriminant islands for each lifestyle.