Recognizing the pseudogenes in bacterial genomes

The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

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.

Genomic insights into local-scale evolution of ocular Chlamydia trachomatis strains within and between individuals in Gambian trachoma-endemic villages

Trachoma, a neglected tropical disease caused by Chlamydia trachomatis (Ct) serovars A-C, is the leading infectious cause of blindness worldwide. Africa bears the highest burden, accounting for over 86 % of global trachoma cases. We investigated Ct serovar A (SvA) and B (SvB) whole genome sequences prior to the induction of mass antibiotic drug administration in The Gambia. Here, we explore the factors contributing to Ct strain diversification and the implications for Ct evolution within the context of ocular infection. A cohort study in 2002-2003 collected ocular swabs across nine Gambian villages during a 6 month follow-up study. To explore the genetic diversity of Ct within and between individuals, we conducted whole-genome sequencing (WGS) on a limited number (n=43) of Ct-positive samples with an omcB load ≥10 from four villages. WGS was performed using target enrichment with SureSelect and Illumina paired-end sequencing. Out of 43 WGS samples, 41 provided sufficient quality for further analysis. ompA analysis revealed that 11 samples had highest identity to ompA from strain A/HAR13 (NC_007429) and 30 had highest identity to ompA from strain B/Jali20 (NC_012686). While SvB genome sequences formed two distinct village-driven subclades, the heterogeneity of SvA sequences led to the formation of many individual branches within the Gambian SvA subclade. Comparing the Gambian SvA and SvB sequences with their reference strains, Ct A/HAR13 and Ct B/Jali20, indicated an single nucleotide polymorphism accumulation rate of 2.4×10-5 per site per year for the Gambian SvA and 1.3×10-5 per site per year for SvB variants (P<0.0001). Variant calling resulted in a total of 1371 single nucleotide variants (SNVs) with a frequency >25 % in SvA sequences, and 438 SNVs in SvB sequences. Of note, in SvA variants, highest evolutionary pressure was recorded on genes responsible for host cell modulation and intracellular survival mechanisms, whereas in SvB variants this pressure was mainly on genes essential for DNA replication/repair mechanisms and protein synthesis. A comparison of the sequences between observed separate infection events (4-20 weeks between infections) suggested that the majority of the variations accumulated in genes responsible for host-pathogen interaction such as CTA_0166 (phospholipase D-like protein), CTA_0498 (TarP) and CTA_0948 (deubiquitinase). This comparison of Ct SvA and SvB variants within a trachoma endemic population focused on their local evolutionary adaptation. We found a different variation accumulation pattern in the Gambian SvA chromosomal genes compared with SvB, hinting at the potential of Ct serovar-specific variation in diversification and evolutionary fitness. These findings may have implications for optimizing trachoma control and prevention strategies.

Comprehensive Proteomic Analysis of Brucella melitensis ATCC23457 Strain Reveals Metabolic Adaptations in Response to Nutrient Stress

In the present study, a comprehensive proteomic analysis of Brucella melitensis (B. melitensis) strain ATCC23457 was carried out to investigate proteome alterations in response to in vitro-induced nutrient stress. Our analysis resulted in the identification of 2440 proteins, including 365 hypothetical proteins and 850 potentially secretory proteins representing ~77.8% of the B. melitensis proteome. Utilizing a proteogenomics approach, we provide translational evidence for eight novel putative protein-coding genes and confirmed the coding potential of 31 putatively annotated pseudogenes, thus refining the existing genome annotation. Further, using a label-free quantitative proteomic approach, new insights into the cellular processes governed by nutrient stress, including enrichment of amino acid metabolism (E), transcription (K), energy production and conversion (C), and biogenesis (J) processes were obtained. Pathway analysis revealed the enrichment of survival and homeostasis maintenance pathways, including type IV secretion system, nitrogen metabolism, and urease pathways in response to nutrient limitation. To conclude, our analysis demonstrates the utility of in-depth proteomic analysis in enabling improved annotation of the B. melitensis genome. Further, our results indicate that B. melitensis undergoes metabolic adaptations during nutrient stress similar to other Brucella. sp, and adapts itself for long-term persistence and survival.

Uncovering Lasonolide A Biosynthesis Using Genome-Resolved Metagenomics

Invertebrates, particularly sponges, have been a dominant source of new marine natural products. For example, lasonolide A (LSA) is a potential anticancer molecule isolated from the marine sponge Forcepia sp., with nanomolar growth inhibitory activity and a unique cytotoxicity profile against the National Cancer Institute 60-cell-line screen. Here, we identified the putative biosynthetic pathway for LSA. Genomic binning of the Forcepia sponge metagenome revealed a Gram-negative bacterium belonging to the phylum Verrucomicrobia as the candidate producer of LSA. Phylogenetic analysis showed that this bacterium, here named "Candidatus Thermopylae lasonolidus," only has 88.78% 16S rRNA identity with the closest relative, Pedosphaera parvula Ellin514, indicating that it represents a new genus. The lasonolide A (las) biosynthetic gene cluster (BGC) was identified as a trans-acyltransferase (AT) polyketide synthase (PKS) pathway. Compared with its host genome, the las BGC exhibits a significantly different GC content and pentanucleotide frequency, suggesting a potential horizontal acquisition of the gene cluster. Furthermore, three copies of the putative las pathway were identified in the candidate producer genome. Differences between the three las repeats were observed, including the presence of three insertions, two single-nucleotide polymorphisms, and the absence of a stand-alone acyl carrier protein in one of the repeats. Even though the verrucomicrobial producer shows signs of genome reduction, its genome size is still fairly large (about 5 Mbp), and, compared to its closest free-living relative, it contains most of the primary metabolic pathways, suggesting that it is in the early stages of reduction. IMPORTANCE While sponges are valuable sources of bioactive natural products, a majority of these compounds are produced in small quantities by uncultured symbionts, hampering the study and clinical development of these unique compounds. Lasonolide A (LSA), isolated from marine sponge Forcepia sp., is a cytotoxic molecule active at nanomolar concentrations, which causes premature chromosome condensation, blebbing, cell contraction, and loss of cell adhesion, indicating a novel mechanism of action and making it a potential anticancer drug lead. However, its limited supply hampers progression to clinical trials. We investigated the microbiome of Forcepia sp. using culture-independent DNA sequencing, identified genes likely responsible for LSA synthesis in an uncultured bacterium, and assembled the symbiont's genome. These insights provide future opportunities for heterologous expression and cultivation efforts that may minimize LSA's supply problem.

The resistomes of Mycobacteroides abscessus complex and their possible acquisition from horizontal gene transfer

Background

Mycobacteroides abscessus complex (MABC), an emerging pathogen, causes human infections resistant to multiple antibiotics. In this study, the genome data of 1,581 MABC strains were downloaded from NCBI database for phylogenetic relatedness inference, resistance profile identification and the estimation of evolutionary pressure on resistance genes in silico.

Results

From genes associated with resistance to 28 antibiotic classes, 395 putative proteins (ARPs) were identified, based on the information in two antibiotic resistance databases (CARD and ARG-ANNOT). The ARPs most frequently identified in MABC were those associated with resistance to multiple antibiotic classes, beta-lactams and aminoglycosides. After excluding ARPs that had undergone recombination, two ARPs were predicted to be under diversifying selection and 202 under purifying selection. This wide occurrence of purifying selection suggested that the diversity of commonly shared ARPs in MABC have been reduced to achieve stability. The unequal distribution of ARPs in members of the MABC could be due to horizontal gene transfer or ARPs pseudogenization events. Most (81.5%) of the ARPs were observed in the accessory genome and 72.2% ARPs were highly homologous to proteins associated with mobile genetic elements such as plasmids, prophages and viruses. On the other hand, with TBLASTN search, only 18 of the ARPs were identified as pseudogenes.

Conclusion

Altogether, our results suggested an important role of horizontal gene transfer in shaping the resistome of MABC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08941-7.

Comparative analysis of Vibrio cholerae isolates from Ghana reveals variations in genome architecture and adaptation of outbreak and environmental strains

Recurrent epidemics of cholera denote robust adaptive mechanisms of Vibrio cholerae for ecological shifting and persistence despite variable stress conditions. Tracking the evolution of pathobiological traits requires comparative genomic studies of isolates from endemic areas. Here, we investigated the genetic differentiation among V. cholerae clinical and environmental isolates by highlighting the genomic divergence associated with gene decay, genome plasticity, and the acquisition of virulence and adaptive traits. The clinical isolates showed high phylogenetic relatedness due to a higher frequency of shared orthologs and fewer gene variants in contrast to the evolutionarily divergent environmental strains. Divergence of the environmental isolates is linked to extensive genomic rearrangements in regions containing mobile genetic elements resulting in numerous breakpoints, relocations, and insertions coupled with the loss of virulence determinants acf, zot, tcp, and ctx in the genomic islands. Also, four isolates possessed the CRISPR-Cas systems with spacers specific for Vibrio phages and plasmids. Genome synteny and homology analysis of the CRISPR-Cas systems suggest horizontal acquisition. The marked differences in the distribution of other phage and plasmid defense systems such as Zorya, DdmABC, DdmDE, and type-I Restriction Modification systems among the isolates indicated a higher propensity for plasmid or phage disseminated traits in the environmental isolates. Our results reveal that V. cholerae strains undergo extensive genomic rearrangements coupled with gene acquisition, reflecting their adaptation during ecological shifts and pathogenicity.

The rate and role of pseudogenes of the Mycobacterium tuberculosis complex

Whole-genome sequence analyses have significantly contributed to the understanding of virulence and evolution of the Mycobacterium tuberculosis complex (MTBC), the causative pathogens of tuberculosis. Most MTBC evolutionary studies are focused on single nucleotide polymorphisms and deletions, but rare studies have evaluated gene content, whereas none has comprehensively evaluated pseudogenes. Accordingly, we describe an extensive study focused on quantifying and predicting possible functions of MTBC and Mycobacterium canettii pseudogenes. Using NCBI's PGAP-detected pseudogenes, we analysed 25 837 pseudogenes from 158 MTBC and M. canetii strains and combined transcriptomics and proteomics of M. tuberculosis H37Rv to gain insights about pseudogenes' expression. Our results indicate significant variability concerning rate and conservancy of in silico predicted pseudogenes among different ecotypes and lineages of tuberculous mycobacteria and pseudogenization of important virulence factors and genes of the metabolism and antimicrobial resistance/tolerance. We show that in silico predicted pseudogenes contribute considerably to MTBC genetic diversity at the population level. Moreover, the transcription machinery of M. tuberculosis can fully transcribe most pseudogenes, indicating intact promoters and recent pseudogene evolutionary emergence. Proteomics of M. tuberculosis and close evaluation of mutational lesions driving pseudogenization suggest that few in silico predicted pseudogenes are likely capable of neofunctionalization, nonsense mutation reversal, or phase variation, contradicting the classical definition of pseudogenes. Such findings indicate that genome annotation should be accompanied by proteomics and protein function assays to improve its accuracy. While indels and insertion sequences are the main drivers of the observed mutational lesions in these species, population bottlenecks and genetic drift are likely the evolutionary processes acting on pseudogenes' emergence over time. Our findings unveil a new perspective on MTBC's evolution and genetic diversity.

Unraveling the metabolism of Mycobacterium caprae using comparative genomics

In our laboratory, Mycobacterium caprae has poor growth in standard medium (SM) 7H9-OADC supplemented with pyruvate and Tween-80. Our objectives were to identify mutations affecting M. caprae metabolism and use this information to design a culture medium to improve its growth. We selected 77 M. caprae genomes and sequenced M. caprae NLA000201913 used in our experiments. Mutations present in >95% of the strains compared to Mycobacterium tuberculosis H37Rv were analyzed in silico for their deleterious effects on proteins of metabolic pathways. Apart from the known defect in the pyruvate kinase, M. caprae has important lesions in enzymes of the TCA cycle, methylmalonyl cycle, B₁₂ metabolism, and electron-transport chain. We provide evidence of enzymatic redundancy elimination and epistatic mutations, and possible production of toxic metabolites hindering M. caprae growth in vitro. A newly designed SM supplemented with l-glutamate allowed faster growth and increased final microbial mass of M. caprae. However, possible accumulation of metabolic waste-products and/or nutritional limitations halted M. caprae growth prior to a M. tuberculosis-like stationary phase. Our findings suggest that M. caprae relies on GABA and/or glyoxylate shunts for in vitro growth in routine media. The newly developed medium will improve experiments with this bacterium by allowing faster growth in vitro.

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.

The Di-Symbiotic Systems in the Aphids Sipha maydis and Periphyllus lyropictus Provide a Contrasting Picture of Recent Co-Obligate Nutritional Endosymbiosis in Aphids

Dependence on multiple nutritional bacterial symbionts forming a metabolic unit has repeatedly evolved in many insect species that feed on nutritionally unbalanced diets such as plant sap. This is the case for aphids of the subfamilies Lachninae and Chaitophorinae, which have evolved di-symbiotic systems in which the ancient obligate nutritional symbiont Buchnera aphidicola is metabolically complemented by an additional nutritional symbiont acquired more recently. Deciphering how different symbionts integrate both metabolically and anatomically in such systems is crucial to understanding how complex nutritional symbiotic systems function and evolve. In this study, we sequenced and analyzed the genomes of the symbionts B. aphidicola and Serratia symbiotica associated with the Chaitophorinae aphids Sipha maydis and Periphyllus lyropictus. Our results show that, in these two species, B. aphidicola and S. symbiotica complement each other metabolically (and their hosts) for the biosynthesis of essential amino acids and vitamins, but with distinct metabolic reactions supported by each symbiont depending on the host species. Furthermore, the S. symbiotica symbiont associated with S. maydis appears to be strictly compartmentalized into the specialized host cells housing symbionts in aphids, the bacteriocytes, whereas the S. symbiotica symbiont associated with P. lyropictus exhibits a highly invasive phenotype, presumably because it is capable of expressing a larger set of virulence factors, including a complete flagellum for bacterial motility. Such contrasting levels of metabolic and anatomical integration for two S. symbiotica symbionts that were recently acquired as nutritional co-obligate partners reflect distinct coevolutionary processes specific to each association.

Whokaryote: distinguishing eukaryotic and prokaryotic contigs in metagenomes based on gene structure

Metagenomics has become a prominent technology to study the functional potential of all organisms in a microbial community. Most studies focus on the bacterial content of these communities, while ignoring eukaryotic microbes. Indeed, many metagenomics analysis pipelines silently assume that all contigs in a metagenome are prokaryotic, likely resulting in less accurate annotation of eukaryotes in metagenomes. Early detection of eukaryotic contigs allows for eukaryote-specific gene prediction and functional annotation. Here, we developed a classifier that distinguishes eukaryotic from prokaryotic contigs based on foundational differences between these taxa in terms of gene structure. We first developed Whokaryote, a random forest classifier that uses intergenic distance, gene density and gene length as the most important features. We show that, with an estimated recall, precision and accuracy of 94, 96 and 95 %, respectively, this classifier with features grounded in biology can perform almost as well as the classifiers EukRep and Tiara, which use k-mer frequencies as features. By retraining our classifier with Tiara predictions as an additional feature, the weaknesses of both types of classifiers are compensated; the result is Whokaryote+Tiara, an enhanced classifier that outperforms all individual classifiers, with an F1 score of 0.99 for both eukaryotes and prokaryotes, while still being fast. In a reanalysis of metagenome data from a disease-suppressive plant endospheric microbial community, we show how using Whokaryote+Tiara to select contigs for eukaryotic gene prediction facilitates the discovery of several biosynthetic gene clusters that were missed in the original study. Whokaryote (+Tiara) is wrapped in an easily installable package and is freely available from https://github.com/LottePronk/whokaryote.

"Pseudo-pseudogenes" in bacterial genomes: Proteogenomics reveals a wide but low protein expression of pseudogenes in Salmonella enterica

Pseudogenes (genes disrupted by frameshift or in-frame stop codons) are ubiquitously present in the bacterial genome and considered as nonfunctional fossil. Here, we used RNA-seq and mass-spectrometry technologies to measure the transcriptomes and proteomes of Salmonella enterica serovars Paratyphi A and Typhi. All pseudogenes’ mRNA sequences remained disrupted, and were present at comparable levels to their intact homologs. At the protein level, however, 101 out of 161 pseudogenes suggested successful translation, with their low expression regardless of growth conditions, genetic background and pseudogenization causes. The majority of frameshifting detected was compensatory for -1 frameshift mutations. Readthrough of in-frame stop codons primarily involved UAG; and cytosine was the most frequent base adjacent to the codon. Using a fluorescence reporter system, fifteen pseudogenes were confirmed to express successfully in vivo in Escherichia coli. Expression of the intact copy of the fifteen pseudogenes in S. Typhi affected bacterial pathogenesis as revealed in human macrophage and epithelial cell infection models. The above findings suggest the need to revisit the nonstandard translation mechanism as well as the biological role of pseudogenes in the bacterial genome.

Adaptation and genomic erosion in fragmented Pseudomonas aeruginosa populations in the sinuses of people with cystic fibrosis

Pseudomonas aeruginosa notoriously adapts to the airways of people with cystic fibrosis (CF), yet how infection-site biogeography and associated evolutionary processes vary as lifelong infections progress remains unclear. Here we test the hypothesis that early adaptations promoting aggregation influence evolutionary-genetic trajectories by examining longitudinal P. aeruginosa from the sinuses of six adults with CF. Highly host-adapted lineages harbored mutator genotypes displaying signatures of early genome degradation associated with recent host restriction. Using an advanced imaging technique (MiPACT-HCR [microbial identification after passive clarity technique]), we find population structure tracks with genome degradation, with the most host-adapted, genome-degraded P. aeruginosa (the mutators) residing in small, sparse aggregates. We propose that following initial adaptive evolution in larger populations under strong selection for aggregation, P. aeruginosa persists in small, fragmented populations that experience stronger effects of genetic drift. These conditions enrich for mutators and promote degenerative genome evolution. Our findings underscore the importance of infection-site biogeography to pathogen evolution.

Escherichia coli Small Proteome

Escherichia coli was one of the first species to have its genome sequenced and remains one of the best-characterized model organisms. Thus, it is perhaps surprising that recent studies have shown that a substantial number of genes have been overlooked. Genes encoding more than 140 small proteins, defined as those containing 50 or fewer amino acids, have been identified in E. coli in the past 10 years, and there is substantial evidence indicating that many more remain to be discovered. This review covers the methods that have been successful in identifying small proteins and the short open reading frames that encode them. The small proteins that have been functionally characterized to date in this model organism are also discussed. It is hoped that the review, along with the associated databases of known as well as predicted but undetected small proteins, will aid in and provide a roadmap for the continued identification and characterization of these proteins in E. coli as well as other bacteria.

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.

At the Gate of Mutualism: Identification of Genomic Traits Predisposing to Insect-Bacterial Symbiosis in Pathogenic Strains of the Aphid Symbiont Serratia symbiotica

Mutualistic associations between insects and heritable bacterial symbionts are ubiquitous in nature. The aphid symbiont Serratia symbiotica is a valuable candidate for studying the evolution of bacterial symbiosis in insects because it includes a wide diversity of strains that reflect the diverse relationships in which bacteria can be engaged with insects, from pathogenic interactions to obligate intracellular mutualism. The recent discovery of culturable strains, which are hypothesized to resemble the ancestors of intracellular strains, provide an opportunity to study the mechanisms underlying bacterial symbiosis in its early stages. In this study, we analyzed the genomes of three of these culturable strains that are pathogenic to aphid hosts, and performed comparative genomic analyses including mutualistic host-dependent strains. All three genomes are larger than those of the host-restricted S. symbiotica strains described so far, and show significant enrichment in pseudogenes and mobile elements, suggesting that these three pathogenic strains are in the early stages of the adaptation to their host. Compared to their intracellular mutualistic relatives, the three strains harbor a greater diversity of genes coding for virulence factors and metabolic pathways, suggesting that they are likely adapted to infect new hosts and are a potential source of metabolic innovation for insects. The presence in their genomes of secondary metabolism gene clusters associated with the production of antimicrobial compounds and phytotoxins supports the hypothesis that S. symbiotia symbionts evolved from plant-associated strains and that plants may serve as intermediate hosts. Mutualistic associations between insects and bacteria are the result of independent transitions to endosymbiosis initiated by the acquisition of environmental progenitors. In this context, the genomes of free-living S. symbiotica strains provide a rare opportunity to study the inventory of genes held by bacterial associates of insects that are at the gateway to a host-dependent lifestyle.

Analysis of Stop Codons within Prokaryotic Protein-Coding Genes Suggests Frequent Readthrough Events

Nonsense mutations turn a coding (sense) codon into an in-frame stop codon that is assumed to result in a truncated protein product. Thus, nonsense substitutions are the hallmark of pseudogenes and are used to identify them. Here we show that in-frame stop codons within bacterial protein-coding genes are widespread. Their evolutionary conservation suggests that many of them are not pseudogenes, since they maintain dN/dS values (ratios of substitution rates at non-synonymous and synonymous sites) significantly lower than 1 (this is a signature of purifying selection in protein-coding regions). We also found that double substitutions in codons—where an intermediate step is a nonsense substitution—show a higher rate of evolution compared to null models, indicating that a stop codon was introduced and then changed back to sense via positive selection. This further supports the notion that nonsense substitutions in bacteria are relatively common and do not necessarily cause pseudogenization. In-frame stop codons may be an important mechanism of regulation: Such codons are likely to cause a substantial decrease of protein expression levels.

Computational Methods for Pseudogene Annotation Based on Sequence Homology

The number of complete genome sequences explodes more and more with each passing year. Thus, methods for genome annotation need to be honed constantly to handle the deluge of information. Annotation of pseudogenes (i.e., gene copies that appear not to make a functional protein) in genomes is a persistent problem; here, we overview pseudogene annotation methods that are based on the detection of sequence homology in genomic DNA.

An assessment of genome annotation coverage across the bacterial tree of life

Although gene-finding in bacterial genomes is relatively straightforward, the automated assignment of gene function is still challenging, resulting in a vast quantity of hypothetical sequences of unknown function. But how prevalent are hypothetical sequences across bacteria, what proportion of genes in different bacterial genomes remain unannotated, and what factors affect annotation completeness? To address these questions, we surveyed over 27 000 bacterial genomes from the Genome Taxonomy Database, and measured genome annotation completeness as a function of annotation method, taxonomy, genome size, 'research bias' and publication date. Our analysis revealed that 52 and 79 % of the average bacterial proteome could be functionally annotated based on protein and domain-based homology searches, respectively. Annotation coverage using protein homology search varied significantly from as low as 14 % in some species to as high as 98 % in others. We found that taxonomy is a major factor influencing annotation completeness, with distinct trends observed across the microbial tree (e.g. the lowest level of completeness was found in the Patescibacteria lineage). Most lineages showed a significant association between genome size and annotation incompleteness, likely reflecting a greater degree of uncharacterized sequences in 'accessory' proteomes than in 'core' proteomes. Finally, research bias, as measured by publication volume, was also an important factor influencing genome annotation completeness, with early model organisms showing high completeness levels relative to other genomes in their own taxonomic lineages. Our work highlights the disparity in annotation coverage across the bacterial tree of life and emphasizes a need for more experimental characterization of accessory proteomes as well as understudied lineages.

Accurate reconstruction of bacterial pan- and core genomes with PEPPAN

Bacterial genomes can contain traces of a complex evolutionary history, including extensive homologous recombination, gene loss, gene duplications, and horizontal gene transfer. To reconstruct the phylogenetic and population history of a set of multiple bacteria, it is necessary to examine their pangenome, the composite of all the genes in the set. Here we introduce PEPPAN, a novel pipeline that can reliably construct pangenomes from thousands of genetically diverse bacterial genomes that represent the diversity of an entire genus. PEPPAN outperforms existing pangenome methods by providing consistent gene and pseudogene annotations extended by similarity-based gene predictions, and identifying and excluding paralogs by combining tree- and synteny-based approaches. The PEPPAN package additionally includes PEPPAN_parser, which implements additional downstream analyses, including the calculation of trees based on accessory gene content or allelic differences between core genes. To test the accuracy of PEPPAN, we implemented SimPan, a novel pipeline for simulating the evolution of bacterial pangenomes. We compared the accuracy and speed of PEPPAN with four state-of-the-art pangenome pipelines using both empirical and simulated data sets. PEPPAN was more accurate and more specific than any of the other pipelines and was almost as fast as any of them. As a case study, we used PEPPAN to construct a pangenome of approximately 40,000 genes from 3052 representative genomes spanning at least 80 species of Streptococcus The resulting gene and allelic trees provide an unprecedented overview of the genomic diversity of the entire Streptococcus genus.

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.

Silence as a way of niche adaptation: mecC-MRSA with variations in the accessory gene regulator (agr) functionality express kaleidoscopic phenotypes

Functionality of the accessory gene regulator (agr) quorum sensing system is an important factor promoting either acute or chronic infections by the notorious opportunistic human and veterinary pathogen Staphylococcus aureus. Spontaneous alterations of the agr system are known to frequently occur in human healthcare-associated S. aureus lineages. However, data on agr integrity and function are sparse regarding other major clonal lineages. Here we report on the agr system functionality and activity level in mecC-carrying methicillin resistant S. aureus (MRSA) of various animal origins (n = 33) obtained in Europe as well as in closely related human isolates (n = 12). Whole genome analysis assigned all isolates to four clonal complexes (CC) with distinct agr types (CC599 agr I, CC49 agr II, CC130 agr III and CC1943 agr IV). Agr functionality was assessed by a combination of phenotypic assays and proteome analysis. In each CC, isolates with varying agr activity levels were detected, including the presence of completely non-functional variants. Genomic comparison of the agr I-IV encoding regions associated these phenotypic differences with variations in the agrA and agrC genes. The genomic changes were detected independently in divergent lineages, suggesting that agr variation might foster viability and adaptation of emerging MRSA lineages to distinct ecological niches.

Comparative genomics reveals different population structures associated with host and geographic origin in antimicrobial-resistant Salmonella enterica

Genetic variation in a pathogen, including the causative agent of salmonellosis, Salmonella enterica, can occur as a result of eco-evolutionary forces triggered by dissimilarities of ecological niches. Here, we applied comparative genomics to study 90 antimicrobial resistant (AMR) S. enterica isolates from bovine and human hosts in New York and Washington states to understand host- and geographic-associated population structure. Results revealed distinct presence/absence profiles of functional genes and pseudogenes (e.g., virulence genes) associated with bovine and human isolates. Notably, bovine isolates contained significantly more transposase genes but fewer transposase pseudogenes than human isolates, suggesting the occurrence of large-scale transposition in genomes of bovine and human isolates at different times. The high correlation between transposase genes and AMR genes, as well as plasmid replicons, highlights the potential role of horizontally transferred transposons in promoting adaptation to antibiotics. By contrast, a number of potentially geographic-associated single-nucleotide polymorphisms (SNPs), rather than geographic-associated genes, were identified. Interestingly, 38% of these SNPs were in genes annotated as cell surface protein-encoding genes, including some essential for antibiotic resistance and host colonization. Overall, different evolutionary forces and limited recent inter-population transmission appear to shape AMR S. enterica population structure in different hosts and geographic origins.

Horizontal Gene Transfer to a Defensive Symbiont with a Reduced Genome in a Multipartite Beetle Microbiome

Symbiotic mutualisms of bacteria and animals are ubiquitous in nature, running a continuum from facultative to obligate from the perspectives of both partners. The loss of functions required for living independently but not within a host gives rise to reduced genomes in many symbionts. Although the phenomenon of genome reduction can be explained by existing evolutionary models, the initiation of the process is not well understood. Here, we describe the microbiome associated with the eggs of the beetle Lagria villosa, consisting of multiple bacterial symbionts related to Burkholderia gladioli, including a reduced-genome symbiont thought to be the exclusive producer of the defensive compound lagriamide. We show that the putative lagriamide-producing symbiont is the only member of the microbiome undergoing genome reduction and that it has already lost the majority of its primary metabolism and DNA repair pathways. The key step preceding genome reduction in the symbiont was likely the horizontal acquisition of the putative lagriamide lga biosynthetic gene cluster. Unexpectedly, we uncovered evidence of additional horizontal transfers to the symbiont's genome while genome reduction was occurring and despite a current lack of genes needed for homologous recombination. These gene gains may have given the genome-reduced symbiont a selective advantage in the microbiome, especially given the maintenance of the large lga gene cluster despite ongoing genome reduction.IMPORTANCE Associations between microorganisms and an animal, plant, or fungal host can result in increased dependence over time. This process is due partly to the bacterium not needing to produce nutrients that the host provides, leading to loss of genes that it would need to live independently and to a consequent reduction in genome size. It is often thought that genome reduction is aided by genetic isolation-bacteria that live in monocultures in special host organs, or inside host cells, have less access to other bacterial species from which they can obtain genes. Here, we describe exposure of a genome-reduced beetle symbiont to a community of related bacteria with nonreduced genomes. We show that the symbiont has acquired genes from other bacteria despite going through genome reduction, suggesting that isolation has not yet played a major role in this case of genome reduction, with horizontal gene gains still offering a potential route for adaptation.

Comparative Genomics of Streptococcus thermophilus Support Important Traits Concerning the Evolution, Biology and Technological Properties of the Species

Streptococcus thermophilus is a major starter for the dairy industry with great economic importance. In this study we analyzed 23 fully sequenced genomes of S. thermophilus to highlight novel aspects of the evolution, biology and technological properties of this species. Pan/core genome analysis revealed that the species has an important number of conserved genes and that the pan genome is probably going to be closed soon. According to whole genome phylogeny and average nucleotide identity (ANI) analysis, most S. thermophilus strains were grouped in two major clusters (i.e., clusters A and B). More specifically, cluster A includes strains with chromosomes above 1.83 Mbp, while cluster B includes chromosomes below this threshold. This observation suggests that strains belonging to the two clusters may be differentiated by gene gain or gene loss events. Furthermore, certain strains of cluster A could be further subdivided in subgroups, i.e., subgroup I (ASCC 1275, DGCC 7710, KLDS SM, MN-BM-A02, and ND07), II (MN-BM-A01 and MN-ZLW-002), III (LMD-9 and SMQ-301), and IV (APC151 and ND03). In cluster B certain strains formed one distinct subgroup, i.e., subgroup I (CNRZ1066, CS8, EPS, and S9). Clusters and subgroups observed for S. thermophilus indicate the existence of lineages within the species, an observation which was further supported to a variable degree by the distribution and/or the architecture of several genomic traits. These would include exopolysaccharide (EPS) gene clusters, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)-CRISPR associated (Cas) systems, as well as restriction-modification (R-M) systems and genomic islands (GIs). Of note, the histidine biosynthetic cluster was found present in all cluster A strains (plus strain NCTC12958^T) but was absent from all strains in cluster B. Other loci related to lactose/galactose catabolism and urea metabolism, aminopeptidases, the majority of amino acid and peptide transporters, as well as amino acid biosynthetic pathways were found to be conserved in all strains suggesting their central role for the species. Our study highlights the necessity of sequencing and analyzing more S. thermophilus complete genomes to further elucidate important aspects of strain diversity within this starter culture that may be related to its application in the dairy industry.

Adaptation of host transmission cycle during Clostridium difficile speciation

Bacterial speciation is a fundamental evolutionary process characterized by diverging genotypic and phenotypic properties. However, the selective forces that affect genetic adaptations and how they relate to the biological changes that underpin the formation of a new bacterial species remain poorly understood. Here, we show that the spore-forming, healthcare-associated enteropathogen Clostridium difficile is actively undergoing speciation. Through large-scale genomic analysis of 906 strains, we demonstrate that the ongoing speciation process is linked to positive selection on core genes in the newly forming species that are involved in sporulation and the metabolism of simple dietary sugars. Functional validation shows that the new C. difficile produces spores that are more resistant and have increased sporulation and host colonization capacity when glucose or fructose is available for metabolism. Thus, we report the formation of an emerging C. difficile species, selected for metabolizing simple dietary sugars and producing high levels of resistant spores, that is adapted for healthcare-mediated transmission.

Evolution of a major bovine mastitic genotype (rpoB sequence type 10-2) of Staphylococcus aureus in cows

Staphylococcus aureus is the major pathogen leading to bovine mastitis globally while livestock-associated methicillin resistant S. aureus (LA-MRSA) has become a potential threat to public health. MRSA from bovine mastitis is not common but a methicillin susceptible S. aureus (MSSA) genotype, rpoB sequence type (RST)10-2 (RST10-2), is prevalent in Korea. To date, many genomic sequences from S. aureus have been elucidated, but the complete genome sequences of RST10-2 MSSA from bovine mastitis has never been reported. In this study, we determined the complete genome sequence of two RST10-2 MSSA that differ from each other in staphylococcal protein A and molecular prophage types [PMB64-1 (t2489/ mPPT0) and PMB81-4 (t127/mPPT1-2-3)] and conducted a comparative genomics study. The genomic sequences of PMB64-1 and PMB81-4 were more homologous to the representative human RST10-2 strains (MSSA476, MW2 etc.) compared to other RSTs. Most of them shared five common pseudogenes, along with high amino acid identity of four variable virulence genes that were identified in this study. However, PMB64-1 and PMB81-4 acquired different strainspecific pseudogenes and mobile genetic elements than the human strains. The unique pseudogene profile and high identity of the virulence genes were verified in RST10-2 field strains from bovine mastitis. Thus, bovine mastitic RST10-2 MSSA may have an evolutionary relationship with the human RST10-2 community-associated (CA) MSSA and CA-MRSA strains but may have adapted to cows.

Nanopore ultra-long read sequencing technology for antimicrobial resistance detection in Mannheimia haemolytica

Disruptive innovations in long-range, cost-effective direct template nucleic acid sequencing are transforming clinical and diagnostic medicine. A multidrug resistant strain and a pan-susceptible strain of Mannheimia haemolytica, isolated from pneumonic bovine lung samples, were sequenced at 146× and 111× coverage, respectively with Oxford Nanopore Technologies MinION. De novo assembly produced a complete genome for the non-resistant strain and a nearly complete assembly for the drug resistant strain. Functional annotation using RAST (Rapid Annotations using Subsystems Technology), CARD (Comprehensive Antibiotic Resistance Database) and ResFinder databases identified genes conferring resistance to different classes of antibiotics including β-lactams, tetracyclines, lincosamides, phenicols, aminoglycosides, sulfonamides and macrolides. Resistance phenotypes of the M. haemolytica strains were determined by minimum inhibitory concentration (MIC) of the antibiotics. Sequencing with a highly portable MinION device corresponded to MIC assays with most of the antimicrobial resistant determinants being identified with as few as 5437 reads, except for the genes responsible for resistance to Fluoroquinolones. The resulting quality assemblies and AMR gene annotation highlight the efficiency of ultra-long read, whole-genome sequencing (WGS) as a valuable tool in diagnostic veterinary medicine.

Factors driving effective population size and pan-genome evolution in bacteria

Background

Knowledge of population-level processes is essential to understanding the efficacy of selection operating within a species. However, attempts at estimating effective population sizes (Ne) are particularly challenging in bacteria due to their extremely large census populations sizes, varying rates of recombination and arbitrary species boundaries.

Results

In this study, we estimated Ne for 153 species (152 bacteria and one archaeon) defined under a common framework and found that ecological lifestyle and growth rate were major predictors of Ne; and that contrary to theoretical expectations, Ne was unaffected by recombination rate. Additionally, we found that Ne shapes the evolution and diversity of total gene repertoires of prokaryotic species.

Conclusion

Together, these results point to a new model of genome architecture evolution in prokaryotes, in which pan-genome sizes, not individual genome sizes, are governed by drift-barrier evolution.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1272-4) contains supplementary material, which is available to authorized users.

Complete Genome Sequence of α-1,3-Glucanase-Producing Strain Paracoccus mutanolyticus RSP-02

A mutanase (α-1,3 glucanase)-producing bacterial strain of Paracoccus mutanolyticus was isolated from soil samples rich in cellulosic waste. Here, we report the whole-genome sequencing and annotation of P. mutanolyticus, which has a genome size of around 3.5 Mb and the potential to degrade water-insoluble α-1,3 glucans with an overall G+C content of 67.4%.

A mutanase (α-1,3 glucanase)-producing bacterial strain of Paracoccus mutanolyticus was isolated from soil samples rich in cellulosic waste. Here, we report the whole-genome sequencing and annotation of P. mutanolyticus, which has a genome size of around 3.5 Mb and the potential to degrade water-insoluble α-1,3 glucans with an overall G+C content of 67.4%.

Microevolution of Streptococcus agalactiae ST-261 from Australia Indicates Dissemination via Imported Tilapia and Ongoing Adaptation to Marine Hosts or Environment

Streptococcus agalactiae (group B Streptococcus [GBS]) causes disease in a wide range of animals. The serotype Ib lineage is highly adapted to aquatic hosts, exhibiting substantial genome reduction compared with terrestrial conspecifics. Here, we sequence genomes from 40 GBS isolates, including 25 isolates from wild fish and captive stingrays in Australia, six local veterinary or human clinical isolates, and nine isolates from farmed tilapia in Honduras, and compared them with 42 genomes from public databases. Phylogenetic analysis based on nonrecombinant core-genome single nucleotide polymorphisms (SNPs) indicated that aquatic serotype Ib isolates from Queensland were distantly related to local veterinary and human clinical isolates. In contrast, Australian aquatic isolates are most closely related to a tilapia isolate from Israel, differing by only 63 core-genome SNPs. A consensus minimum spanning tree based on core-genome SNPs indicates the dissemination of sequence type 261 (ST-261) from an ancestral tilapia strain, which is congruent with several introductions of tilapia into Australia from Israel during the 1970s and 1980s. Pangenome analysis identified 1,440 genes as core, with the majority being dispensable or strain specific, with non-protein-coding intergenic regions (IGRs) divided among core and strain-specific genes. Aquatic serotype Ib strains have lost many virulence factors during adaptation, but six adhesins were well conserved across the aquatic isolates and might be critical for virulence in fish and for targets in vaccine development. The close relationship among recent ST-261 isolates from Ghana, the United States, and China with the Israeli tilapia isolate from 1988 implicates the global trade in tilapia seed for aquaculture in the widespread dissemination of serotype Ib fish-adapted GBS.IMPORTANCEStreptococcus agalactiae (GBS) is a significant pathogen of humans and animals. Some lineages have become adapted to particular hosts, and serotype Ib is highly specialized to fish. Here, we show that this lineage is likely to have been distributed widely by the global trade in tilapia for aquaculture, with probable introduction into Australia in the 1970s and subsequent dissemination in wild fish populations. We report here the variability in the polysaccharide capsule among this lineage but identify a cohort of common surface proteins that may be a focus of future vaccine development to reduce the biosecurity risk in international fish trade.

Machine learning identifies signatures of host adaptation in the bacterial pathogen Salmonella enterica

Emerging pathogens are a major threat to public health, however understanding how pathogens adapt to new niches remains a challenge. New methods are urgently required to provide functional insights into pathogens from the massive genomic data sets now being generated from routine pathogen surveillance for epidemiological purposes. Here, we measure the burden of atypical mutations in protein coding genes across independently evolved Salmonella enterica lineages, and use these as input to train a random forest classifier to identify strains associated with extraintestinal disease. Members of the species fall along a continuum, from pathovars which cause gastrointestinal infection and low mortality, associated with a broad host-range, to those that cause invasive infection and high mortality, associated with a narrowed host range. Our random forest classifier learned to perfectly discriminate long-established gastrointestinal and invasive serovars of Salmonella. Additionally, it was able to discriminate recently emerged Salmonella Enteritidis and Typhimurium lineages associated with invasive disease in immunocompromised populations in sub-Saharan Africa, and within-host adaptation to invasive infection. We dissect the architecture of the model to identify the genes that were most informative of phenotype, revealing a common theme of degradation of metabolic pathways in extraintestinal lineages. This approach accurately identifies patterns of gene degradation and diversifying selection specific to invasive serovars that have been captured by more labour-intensive investigations, but can be readily scaled to larger analyses.

REPARATION: ribosome profiling assisted (re-)annotation of bacterial genomes

Prokaryotic genome annotation is highly dependent on automated methods, as manual curation cannot keep up with the exponential growth of sequenced genomes. Current automated methods depend heavily on sequence composition and often underestimate the complexity of the proteome. We developed RibosomeE Profiling Assisted (re-)AnnotaTION (REPARATION), a de novo machine learning algorithm that takes advantage of experimental protein synthesis evidence from ribosome profiling (Ribo-seq) to delineate translated open reading frames (ORFs) in bacteria, independent of genome annotation (https://github.com/Biobix/REPARATION). REPARATION evaluates all possible ORFs in the genome and estimates minimum thresholds based on a growth curve model to screen for spurious ORFs. We applied REPARATION to three annotated bacterial species to obtain a more comprehensive mapping of their translation landscape in support of experimental data. In all cases, we identified hundreds of novel (small) ORFs including variants of previously annotated ORFs and >70% of all (variants of) annotated protein coding ORFs were predicted by REPARATION to be translated. Our predictions are supported by matching mass spectrometry proteomics data, sequence composition and conservation analysis. REPARATION is unique in that it makes use of experimental translation evidence to intrinsically perform a de novo ORF delineation in bacterial genomes irrespective of the sequence features linked to open reading frames.

Increased Biosynthetic Gene Dosage in a Genome-Reduced Defensive Bacterial Symbiont

Secondary metabolites, which are small-molecule organic compounds produced by living organisms, provide or inspire drugs for many different diseases. These natural products have evolved over millions of years to provide a survival benefit to the producing organism and often display potent biological activity with important therapeutic applications. For instance, defensive compounds in the environment may be cytotoxic to eukaryotic cells, a property exploitable for cancer treatment. Here, we describe the genome of an uncultured symbiotic bacterium that makes such a cytotoxic metabolite. This symbiont is losing genes that do not endow a selective advantage in a hospitable host environment. Secondary metabolism genes, however, are repeated multiple times in the genome, directly demonstrating their selective advantage. This finding shows the strength of selective forces in symbiotic relationships and suggests that uncultured bacteria in such relationships should be targeted for drug discovery efforts.

A symbiotic lifestyle frequently results in genome reduction in bacteria; the isolation of small populations promotes genetic drift and the fixation of deletions and deleterious mutations over time. Transitions in lifestyle, including host restriction or adaptation to an intracellular habitat, are thought to precipitate a wave of sequence degradation events and consequent proliferation of pseudogenes. We describe here a verrucomicrobial symbiont of the tunicate Lissoclinum sp. that appears to be undergoing such a transition, with low coding density and many identifiable pseudogenes. However, despite the overall drive toward genome reduction, this symbiont maintains seven copies of a large polyketide synthase (PKS) pathway for the mandelalides (mnd), cytotoxic compounds that likely constitute a chemical defense for the host. There is evidence of ongoing degradation in a small number of these repeats—including variable borders, internal deletions, and single nucleotide polymorphisms (SNPs). However, the gene dosage of most of the pathway is increased at least 5-fold. Correspondingly, this single pathway accounts for 19% of the genome by length and 25.8% of the coding capacity. This increased gene dosage in the face of generalized sequence degradation and genome reduction suggests that mnd genes are under strong purifying selection and are important to the symbiotic relationship.

IMPORTANCE Secondary metabolites, which are small-molecule organic compounds produced by living organisms, provide or inspire drugs for many different diseases. These natural products have evolved over millions of years to provide a survival benefit to the producing organism and often display potent biological activity with important therapeutic applications. For instance, defensive compounds in the environment may be cytotoxic to eukaryotic cells, a property exploitable for cancer treatment. Here, we describe the genome of an uncultured symbiotic bacterium that makes such a cytotoxic metabolite. This symbiont is losing genes that do not endow a selective advantage in a hospitable host environment. Secondary metabolism genes, however, are repeated multiple times in the genome, directly demonstrating their selective advantage. This finding shows the strength of selective forces in symbiotic relationships and suggests that uncultured bacteria in such relationships should be targeted for drug discovery efforts.

Author Video: An author video summary of this article is available.

Novel coordination of lipopolysaccharide modifications in Vibrio cholerae promotes CAMP resistance

In the environment and during infection, the human intestinal pathogen Vibrio cholerae must overcome noxious compounds that damage the bacterial outer membrane. The El Tor and classical biotypes of O1 V. cholerae show striking differences in their resistance to membrane disrupting cationic antimicrobial peptides (CAMPs), such as polymyxins. The classical biotype is susceptible to CAMPs, but current pandemic El Tor biotype isolates gain CAMP resistance by altering the net charge of their cell surface through glycine modification of lipid A. Here we report a second lipid A modification mechanism that only functions in the V. cholerae El Tor biotype. We identify a functional EptA ortholog responsible for the transfer of the amino-residue phosphoethanolamine (pEtN) to the lipid A of V. cholerae El Tor that is not functional in the classical biotype. We previously reported that mildly acidic growth conditions (pH 5.8) downregulate expression of genes encoding the glycine modification machinery. In this report, growth at pH 5.8 increases expression of eptA with concomitant pEtN modification suggesting coordinated regulation of these LPS modification systems. Similarly, efficient pEtN lipid A substitution is seen in the absence of lipid A glycinylation. We further demonstrate EptA orthologs from non-cholerae Vibrio species are functional.

The Evolution of Bacterial Genome Architecture

The genome architecture of bacteria and eukaryotes evolves in opposite directions when subject to genetic drift, a difference that can be ascribed to the fact that bacteria exhibit a mutational bias that deletes superfluous sequences, whereas eukaryotes are biased toward large insertions. Expansion of eukaryotic genomes occurs through the addition of non-functional sequences, such as repetitive sequences and transposable elements, whereas variation in bacterial genome size is largely due to the acquisition and loss of functional accessory genes. These properties create the situation in which eukaryotes with very similar numbers of genes can have vastly different genome sizes, while in bacteria, gene number scales linearly with genome size. Some bacterial genomes, however, particularly those of species that undergo bottlenecks due to recent association with hosts, accumulate pseudogenes and mobile elements, conferring them a low gene content relative to their genome size. These non-functional sequences are gradually eroded and eliminated after long-term association with hosts, with the result that obligate symbionts have the smallest genomes of any cellular organism. The architecture of bacterial genomes is shaped by complex and diverse processes, but for most bacterial species, genome size is governed by a non-adaptive process, i.e., genetic drift coupled with a mutational bias toward deletions. Thus, bacteria with small effective population sizes typically have the smallest genomes. Some marine bacteria counter this near-universal trend: despite having immense population sizes, selection, not drift, acts to reduce genome size in response to metabolic constraints in their nutrient-limited environment.

Natural coinfection by Streptococcus agalactiae and Francisella noatunensis subsp. orientalis in farmed Nile tilapia (Oreochromis niloticus L.)

Streptococcus agalactiae and Francisella noatunensis subsp. orientalis (Fno) are important pathogens for farm-raised tilapia worldwide. There are no reports of coinfection caused by S. agalactiae and Fno in fish. This study aimed to determine the aetiology of atypical mortalities in a cage farm of Nile tilapia and to characterize the genetic diversity of the isolates. Fifty-two fish were sampled and subjected to parasitological and bacteriological examination. The S. agalactiae and Fno isolates were genotyped using MLST and REP-PCR, respectively. Whole-genome sequencing was performed to confirm the MLST results. Seven fish were shown coinfected by S. agalactiae and Fno. Chronic hypoxia and a reduction in the water temperature were determined as risk factors for coinfection. Fno isolates were shown clonally related in REP-PCR. The MLST analysis revealed that the S. agalactiae isolates from seven coinfected fish were negative for the glcK gene; however, these were determined to be members of clonal complex CC-552. This is the first description of coinfection by S. agalactiae and Fno in farm-raised Nile tilapia. The coinfection was predisposed by chronic hypoxia and was caused by the main genotypes of S. agalactiae and Fno reported in Brazil. Finally, a new S. agalactiae genotype with glcK gene partially deleted was described.

Comparative Genomics of H. pylori and Non-Pylori Helicobacter Species to Identify New Regions Associated with Its Pathogenicity and Adaptability

The genus Helicobacter is a group of Gram-negative, helical-shaped pathogens consisting of at least 36 bacterial species. Helicobacter pylori (H. pylori), infecting more than 50% of the human population, is considered as the major cause of gastritis, peptic ulcer, and gastric cancer. However, the genetic underpinnings of H. pylori that are responsible for its large scale epidemic and gastrointestinal environment adaption within human beings remain unclear. Core-pan genome analysis was performed among 75 representative H. pylori and 24 non-pylori Helicobacter genomes. There were 1173 conserved protein families of H. pylori and 673 of all 99 Helicobacter genus strains. We found 79 genome unique regions, a total of 202,359bp, shared by at least 80% of the H. pylori but lacked in non-pylori Helicobacter species. The operons, genes, and sRNAs within the H. pylori unique regions were considered as potential ones associated with its pathogenicity and adaptability, and the relativity among them has been partially confirmed by functional annotation analysis. However, functions of at least 54 genes and 10 sRNAs were still unclear. Our analysis of protein-protein interaction showed that 30 genes within them may have the cooperation relationship.

Lack of Overt Genome Reduction in the Bryostatin-Producing Bryozoan Symbiont "Candidatus Endobugula sertula"

The uncultured bacterial symbiont "Candidatus Endobugula sertula" is known to produce cytotoxic compounds called bryostatins, which protect the larvae of its host, Bugula neritina The symbiont has never been successfully cultured, and it was thought that its genome might be significantly reduced. Here, we took a shotgun metagenomics and metatranscriptomics approach to assemble and characterize the genome of "Ca Endobugula sertula." We found that it had specific metabolic deficiencies in the biosynthesis of certain amino acids but few other signs of genome degradation, such as small size, abundant pseudogenes, and low coding density. We also identified homologs to genes associated with insect pathogenesis in other gammaproteobacteria, and these genes may be involved in host-symbiont interactions and vertical transmission. Metatranscriptomics revealed that these genes were highly expressed in a reproductive host, along with bry genes for the biosynthesis of bryostatins. We identified two new putative bry genes fragmented from the main bry operon, accounting for previously missing enzymatic functions in the pathway. We also determined that a gene previously assigned to the pathway, bryS, is not expressed in reproductive tissue, suggesting that it is not involved in the production of bryostatins. Our findings suggest that "Ca Endobugula sertula" may be able to live outside the host if its metabolic deficiencies are alleviated by medium components, which is consistent with recent findings that it may be possible for "Ca Endobugula sertula" to be transmitted horizontally.

IMPORTANCE

The bryostatins are potent protein kinase C activators that have been evaluated in clinical trials for a number of indications, including cancer and Alzheimer's disease. There is, therefore, considerable interest in securing a renewable supply of these compounds, which is currently only possible through aquaculture of Bugula neritina and total chemical synthesis. However, these approaches are labor-intensive and low-yielding and thus preclude the use of bryostatins as a viable therapeutic agent. Our genome assembly and transcriptome analysis for "Ca Endobugula sertula" shed light on the metabolism of this symbiont, potentially aiding isolation and culturing efforts. Our identification of additional bry genes may also facilitate efforts to express the complete pathway heterologously.

An emerging zoonotic clone in the Netherlands provides clues to virulence and zoonotic potential of Streptococcus suis

Streptococcus suis is a zoonotic swine pathogen and a major public health concern in Asia, where it emerged as an important cause of bacterial meningitis in adults. While associated with food-borne transmission in Asia, zoonotic S. suis infections are mainly occupational hazards elsewhere. To identify genomic differences that can explain zoonotic potential, we compared whole genomes of 98 S. suis isolates from human patients and pigs with invasive disease in the Netherlands, and validated our observations with 18 complete and publicly available sequences. Zoonotic isolates have smaller genomes than non-zoonotic isolates, but contain more virulence factors. We identified a zoonotic S. suis clone that diverged from a non-zoonotic clone by means of gene loss, a capsule switch, and acquisition of a two-component signalling system in the late 19th century, when foreign pig breeds were introduced. Our results indicate that zoonotic potential of S. suis results from gene loss, recombination and horizontal gene transfer events.

Reannotation of Yersinia pestis Strain 91001 Based on Omics Data

Yersinia pestis is among the most dangerous human pathogens, and systematic research of this pathogen is important in bacterial pathogenomics research. To fully interpret the biological functions, physiological characteristics, and pathogenesis of Y. pestis, a comprehensive annotation of its entire genome is necessary. The emergence of omics-based research has brought new opportunities to better annotate the genome of this pathogen. Here, the complete genome of Y. pestis strain 91001 was reannotated using genomics and proteogenomics data. One hundred and thirty-seven unreliable coding sequences were removed, and 41 homologous genes were relocated with their translational initiation sites, while the functions of seven pseudogenes and 392 hypothetical genes were revised. Moreover, annotations of noncoding RNAs, repeat sequences, and transposable elements have also been incorporated. The reannotated results are freely available at http://tody.bmi.ac.cn.

Ecological Opportunity, Evolution, and the Emergence of Flea-Borne Plague

The plague bacillus Yersinia pestis is unique among the pathogenic Enterobacteriaceae in utilizing an arthropod-borne transmission route. Transmission by fleabite is a recent evolutionary adaptation that followed the divergence of Y. pestis from the closely related food- and waterborne enteric pathogen Yersinia pseudotuberculosis A combination of population genetics, comparative genomics, and investigations of Yersinia-flea interactions have disclosed the important steps in the evolution and emergence of Y. pestis as a flea-borne pathogen. Only a few genetic changes, representing both gene gain by lateral transfer and gene loss by loss-of-function mutation (pseudogenization), were fundamental to this process. The emergence of Y. pestis fits evolutionary theories that emphasize ecological opportunity in adaptive diversification and rapid emergence of new species.

Evolutionary direction of processed pseudogenes

While some pseudogenes have been reported to play important roles in gene regulation, little is known about the possible relationship between pseudogene functions and evolutionary process of pseudogenes, or about the forces responsible for the pseudogene evolution. In this study, we characterized human processed pseudogenes in terms of evolutionary dynamics. Our results show that pseudogenes tend to evolve toward: lower GC content, strong dinucleotide bias, reduced abundance of transcription factor binding motifs and short palindromes, and decreased ability to form nucleosomes. We explored possible evolutionary forces that shaped the evolution pattern of pseudogenes, and concluded that mutations in pseudogenes are likely determined, at least partially, by neighbor-dependent mutational bias and recombination-associated selection.

Functional profiling of cyanobacterial genomes and its role in ecological adaptations

With the availability of complete genome sequences of many cyanobacterial species, it is becoming feasible to study the broad prospective of the environmental adaptation and the overall changes at transcriptional and translational level in these organisms. In the evolutionary phase, niche-specific competitive forces have resulted in specific features of the cyanobacterial genomes. In this study, functional composition of the 84 different cyanobacterial genomes and their adaptations to different environments was examined by identifying the genomic composition for specific cellular processes, which reflect their genomic functional profile and ecological adaptation. It was identified that among cyanobacterial genomes, metabolic genes have major share over other categories and differentiation of genomic functional profile was observed for the species inhabiting different habitats. The cyanobacteria of freshwater and other habitats accumulate large number of poorly characterized genes. Strain specific functions were also reported in many cyanobacterial members, of which an important feature was the occurrence of phage-related sequences. From this study, it can be speculated that habitat is one of the major factors in giving the shape of functional composition of cyanobacterial genomes towards their ecological adaptations.

Realizing Microbial Evolution

Genome sequences have become the new phenotype for microbial evolutionists. The patterns of diversity revealed in the first 100 bacterial genomes fostered development of a comprehensive framework that can explain their contents, organization, and evolution.

Discovering novel enzymes by functional screening of plurigenomic libraries from alga-associated Flavobacteriia and Gammaproteobacteria

Alga-associated microorganisms, in the context of their numerous interactions with the host and the complexity of the marine environment, are known to produce diverse hydrolytic enzymes with original biochemistry. We recently isolated several macroalgal-polysaccharide-degrading bacteria from the surface of the brown alga Ascophyllum nodosum. These active isolates belong to two classes: the Flavobacteriia and the Gammaproteobacteria. In the present study, we constructed two "plurigenomic" (with multiple bacterial genomes) libraries with the 5 most interesting isolates (regarding their phylogeny and their enzymatic activities) of each class (Fv and Gm libraries). Both libraries were screened for diverse hydrolytic activities. Five activities, out of the 48 previously identified in the natural polysaccharolytic isolates, were recovered by functional screening: a xylanase (GmXyl7), a beta-glucosidase (GmBg1), an esterase (GmEst7) and two iota-carrageenases (Fvi2.5 and Gmi1.3). We discuss here the potential role of the used host-cell, the average DNA insert-sizes and the used restriction enzymes on the divergent screening yields obtained for both libraries and get deeper inside the "great screen anomaly". Interestingly, the discovered esterase probably stands for a novel family of homoserine o-acetyltransferase-like-esterases, while the two iota-carrageenases represent new members of the poorly known GH82 family (containing only 19 proteins since its description in 2000). These original results demonstrate the efficiency of our uncommon "plurigenomic" library approach and the underexplored potential of alga-associated cultivable microbiota for the identification of novel and algal-specific enzymes.

Quantitative Proteomics of the Infectious and Replicative Forms of Chlamydia trachomatis

The obligate intracellular developmental cycle of Chlamydia trachomatis presents significant challenges in defining its proteome. In this study we have applied quantitative proteomics to both the intracellular reticulate body (RB) and the extracellular elementary body (EB) from C. trachomatis. We used C. trachomatis L2 as a model chlamydial isolate for our study since it has a high infectivity:particle ratio and there is an excellent quality genome sequence. EBs and RBs (>99% pure) were quantified by chromosomal and plasmid copy number using PCR, from which the concentrations of chlamydial proteins per bacterial cell/genome were determined. RBs harvested at 15h post infection (PI) were purified by three successive rounds of gradient centrifugation. This is the earliest possible time to obtain purified RBs, free from host cell components in quantity, within the constraints of the technology. EBs were purified at 48h PI. We then used two-dimensional reverse phase UPLC to fractionate RB or EB peptides before mass spectroscopic analysis, providing absolute amount estimates of chlamydial proteins. The ability to express the data as molecules per cell gave ranking in both abundance and energy requirements for synthesis, allowing meaningful identification of rate-limiting components. The study assigned 562 proteins with high confidence and provided absolute estimates of protein concentration for 489 proteins. Interestingly, the data showed an increase in TTS capacity at 15h PI. Most of the enzymes involved in peptidoglycan biosynthesis were detected along with high levels of muramidase (in EBs) suggesting breakdown of peptidoglycan occurs in the non-dividing form of the microorganism. All the genome-encoded enzymes for glycolysis, pentose phosphate pathway and tricarboxylic acid cycle were identified and quantified; these data supported the observation that the EB is metabolically active. The availability of detailed, accurate quantitative proteomic data will be invaluable for investigations into gene regulation and function.

A zero-inflated Poisson model for insertion tolerance analysis of genes based on Tn-seq data

MOTIVATION

Transposon insertion sequencing (Tn-seq) is an emerging technology that combines transposon mutagenesis with next-generation sequencing technologies for the identification of genes related to bacterial survival. The resulting data from Tn-seq experiments consist of sequence reads mapped to millions of potential transposon insertion sites and a large portion of insertion sites have zero mapped reads. Novel statistical method for Tn-seq data analysis is needed to infer functions of genes on bacterial growth.

RESULTS

In this article, we propose a zero-inflated Poisson model for analyzing the Tn-seq data that are high-dimensional and with an excess of zeros. Maximum likelihood estimates of model parameters are obtained using an expectation-maximization (EM) algorithm, and pseudogenes are utilized to construct appropriate statistical tests for the transposon insertion tolerance of normal genes of interest. We propose a multiple testing procedure that categorizes genes into each of the three states, hypo-tolerant, tolerant and hyper-tolerant, while controlling the false discovery rate. We evaluate the proposed method with simulation studies and apply the proposed method to a real Tn-seq data from an experiment that studied the bacterial pathogen, Campylobacter jejuniAvailability and implementation: We provide R code for implementing our proposed method at http://github.com/ffliu/TnSeq A user's guide with example data analysis is also available there.

CONTACT

pliu@iastate.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic dissection of Australian Bordetella pertussis isolates from the 2008-2012 epidemic

OBJECTIVES

Despite high pertussis vaccination coverage, Australia experienced a prolonged epidemic in 2008-2012. The predominant Bordetella pertussis genotype harboured pertussis toxin promoter allele, ptxP3, and pertactin gene allele, prn2. The emergence and expansion of prn non-expressing isolates (Prn negative), were also observed. We aimed to investigate the microevolution and genomic diversity of epidemic B. pertussis isolates.

METHODS

We sequenced 22 B. pertussis isolates collected in 2008-2012 from two states of Australia which are geographically widely separated. Ten of the 22 were Prn negative isolates with three different modes of silencing of prn (prn::IS481F, prn::IS481R and prn::IS1002). Five pre-epidemic isolates were also sequenced for comparison.

RESULTS

Five single nucleotide polymorphisms were common in the epidemic isolates and differentiated them from pre-epidemic isolates. The Australian epidemic isolates can be divided into five lineages (EL1-EL5) with EL1 containing only Prn negative isolates. Comparison with global isolates showed that three lineages remained geographically and temporally distinct whereas two lineages mixed with isolates from 2012 UK outbreak.

CONCLUSION

Our results suggest significant diversification and the microevolution of B. pertussis within the 2008-2012 Australian epidemic.