Massive expansion of Ubiquitination-related gene families within the Chlamydiae

Amoebae as training grounds for microbial pathogens

Grazing of amoebae on microorganisms represents one of the oldest predator-prey dynamic relationships in nature. It represents a genetic "melting pot" for an ancient and continuous multi-directional inter- and intra-kingdom horizontal gene transfer between amoebae and its preys, intracellular microbial residents, endosymbionts, and giant viruses, which has shaped the evolution, selection, and adaptation of microbes that evade degradation by predatory amoeba. Unicellular phagocytic amoebae are thought to be the ancient ancestors of macrophages with highly conserved eukaryotic processes. Selection and evolution of microbes within amoeba through their evolution to target highly conserved eukaryotic processes have facilitated the expansion of their host range to mammals, causing various infectious diseases. Legionella and environmental Chlamydia harbor an immense number of eukaryotic-like proteins that are involved in ubiquitin-related processes or are tandem repeats-containing proteins involved in protein-protein and protein-chromatin interactions. Some of these eukaryotic-like proteins exhibit novel domain architecture and novel enzymatic functions absent in mammalian cells, such as ubiquitin ligases, likely acquired from amoebae. Mammalian cells and amoebae may respond similarly to microbial factors that target highly conserved eukaryotic processes, but mammalian cells may undergo an accidental response to amoeba-adapted microbial factors. We discuss specific examples of microbes that have evolved to evade amoeba predation, including the bacterial pathogens- Legionella, Chlamydia, Coxiella, Rickettssia, Francisella, Mycobacteria, Salmonella, Bartonella, Rhodococcus, Pseudomonas, Vibrio, Helicobacter, Campylobacter, and Aliarcobacter. We also discuss the fungi Cryptococcus, and Asperigillus, as well as amoebae mimiviruses/giant viruses. We propose that amoeba-microbe interactions will continue to be a major "training ground" for the evolution, selection, adaptation, and emergence of microbial pathogens equipped with unique pathogenic tools to infect mammalian hosts. However, our progress will continue to be highly dependent on additional genomic, biochemical, and cellular data of unicellular eukaryotes.

Bacterial usurpation of the OTU deubiquitinase fold

The extensive cellular signalling events controlled by posttranslational ubiquitination are tightly regulated through the action of specialized proteases termed deubiquitinases. Among them, the OTU family of deubiquitinases can play very specialized roles in the regulation of discrete subtypes of ubiquitin signals that control specific cellular functions. To exert control over host cellular functions, some pathogenic bacteria have usurped the OTU deubiquitinase fold as a secreted virulence factor that interferes with ubiquitination inside infected cells. Herein, we provide a review of the function of bacterial OTU deubiquitinases during infection, the structural basis for their deubiquitinase activities and the bioinformatic approaches leading to their identification. Understanding bacterial OTU deubiquitinases holds the potential for discoveries not only in bacterial pathogenesis but in eukaryotic biology as well.

Identification of homologs of the Chlamydia trachomatis effector CteG reveals a family of Chlamydiaceae type III secreted proteins that can be delivered into host cells

Chlamydiae are a large group of obligate endosymbionts of eukaryotes that includes the Chlamydiaceae family, comprising several animal pathogens. Among Chlamydiaceae, Chlamydia trachomatis causes widespread ocular and urogenital infections in humans. Like many bacterial pathogens, all Chlamydiae manipulate host cells by injecting them with type III secretion effector proteins. We previously characterized the C. trachomatis effector CteG, which localizes at the host cell Golgi and plasma membrane during distinct phases of the chlamydial infectious cycle. Here, we show that CteG is a Chlamydiaceae-specific effector with over 60 homologs phylogenetically categorized into two distinct clades (CteG I and CteG II) and exhibiting several inparalogs and outparalogs. Notably, cteG I homologs are syntenic to C. trachomatis cteG, whereas cteG II homologs are syntenic among themselves but not with C. trachomatis cteG. This indicates a complex evolution of cteG homologs, which is unique among C. trachomatis effectors, marked by numerous events of gene duplication and loss. Despite relatively modest sequence conservation, nearly all tested CteG I and CteG II proteins were identified as type III secretion substrates using Yersinia as a heterologous bacterial host. Moreover, most of the type III secreted CteG I and CteG II homologs were delivered by C. trachomatis into host cells, where they localized at the Golgi region and cell periphery. Overall, this provided insights into the evolution of bacterial effectors and revealed a Chlamydiaceae family of type III secreted proteins that underwent substantial divergence during evolution while conserving the capacity to localize at specific host cell compartments.

The SET and ankyrin domains of the secreted Legionella pneumophila histone methyltransferase work together to modify host chromatin

IMPORTANCE

Legionella pneumophila is an intracellular bacterium responsible of Legionnaires' disease, a severe pneumonia that is often fatal when not treated promptly. The pathogen's ability to efficiently colonize the host resides in its ability to replicate intracellularly. Essential for intracellular replication is translocation of many different protein effectors via a specialized secretion system. One of them, called RomA, binds and directly modifies the host chromatin at a unique site (tri-methylation of lysine 14 of histone H3 [H3K14me]). However, the molecular mechanisms of binding are not known. Here, we resolve this question through structural characterization of RomA together with the H3 peptide. We specifically reveal an active role of the ankyrin repeats located in its C-terminal in the interaction with the histone H3 tail. Indeed, without the ankyrin domains, RomA loses its ability to act as histone methyltransferase. These results discover the molecular mechanisms by which a bacterial histone methyltransferase that is conserved in L. pneumophila strains acts to modify chromatin.

The Fish Pathogen "Candidatus Clavichlamydia salmonicola"-A Missing Link in the Evolution of Chlamydial Pathogens of Humans

Chlamydiae like Chlamydia trachomatis and Chlamydia psittaci are well-known human and animal pathogens. Yet, the chlamydiae are a much larger group of evolutionary ancient obligate intracellular bacteria that includes predominantly symbionts of protists and diverse animals. This makes them ideal model organisms to study evolutionary transitions from symbionts in microbial eukaryotes to pathogens of humans. To this end, comparative genome analysis has served as an important tool. Genome sequence data for many chlamydial lineages are, however, still lacking, hampering our understanding of their evolutionary history. Here, we determined the first high-quality draft genome sequence of the fish pathogen "Candidatus Clavichlamydia salmonicola", representing a separate genus within the human and animal pathogenic Chlamydiaceae. The "Ca. Clavichlamydia salmonicola" genome harbors genes that so far have been exclusively found in Chlamydia species suggesting that basic mechanisms important for the interaction with chordate hosts have evolved stepwise in the history of chlamydiae. Thus, the genome sequence of "Ca. Clavichlamydia salmonicola" allows to constrain candidate genes to further understand the evolution of chlamydial virulence mechanisms required to infect mammals.

Ecology and evolution of chlamydial symbionts of arthropods

The phylum Chlamydiae consists of obligate intracellular bacteria including major human pathogens and diverse environmental representatives. Here we investigated the Rhabdochlamydiaceae, which is predicted to be the largest and most diverse chlamydial family, with the few described members known to infect arthropod hosts. Using published 16 S rRNA gene sequence data we identified at least 388 genus-level lineages containing about 14 051 putative species within this family. We show that rhabdochlamydiae are mainly found in freshwater and soil environments, suggesting the existence of diverse, yet unknown hosts. Next, we used a comprehensive genome dataset including metagenome assembled genomes classified as members of the family Rhabdochlamydiaceae, and we added novel complete genome sequences of Rhabdochlamydia porcellionis infecting the woodlouse Porcellio scaber, and of 'Candidatus R. oedothoracis' associated with the linyphiid dwarf spider Oedothorax gibbosus. Comparative analysis of basic genome features and gene content with reference genomes of well-studied chlamydial families with known host ranges, namely Parachlamydiaceae (protist hosts) and Chlamydiaceae (human and other vertebrate hosts) suggested distinct niches for members of the Rhabdochlamydiaceae. We propose that members of the family represent intermediate stages of adaptation of chlamydiae from protists to vertebrate hosts. Within the genus Rhabdochlamydia, pronounced genome size reduction could be observed (1.49-1.93 Mb). The abundance and genomic distribution of transposases suggests transposable element expansion and subsequent gene inactivation as a mechanism of genome streamlining during adaptation to new hosts. This type of genome reduction has never been described before for any member of the phylum Chlamydiae. This study provides new insights into the molecular ecology, genomic diversity, and evolution of representatives of one of the most divergent chlamydial families.

Annotation and Analysis of 3902 Odorant Receptor Protein Sequences from 21 Insect Species Provide Insights into the Evolution of Odorant Receptor Gene Families in Solitary and Social Insects

The gene family of insect olfactory receptors (ORs) has expanded greatly over the course of evolution. ORs enable insects to detect volatile chemicals and therefore play an important role in social interactions, enemy and prey recognition, and foraging. The sequences of several thousand ORs are known, but their specific function or their ligands have only been identified for very few of them. To advance the functional characterization of ORs, we have assembled, curated, and aligned the sequences of 3902 ORs from 21 insect species, which we provide as an annotated online resource. Using functionally characterized proteins from the fly Drosophila melanogaster, the mosquito Anopheles gambiae and the ant Harpegnathos saltator, we identified amino acid positions that best predict response to ligands. We examined the conservation of these predicted relevant residues in all OR subfamilies; the results showed that the subfamilies that expanded strongly in social insects had a high degree of conservation in their binding sites. This suggests that the ORs of social insect families are typically finely tuned and exhibit sensitivity to very similar odorants. Our novel approach provides a powerful tool to exploit functional information from a limited number of genes to study the functional evolution of large gene families.

The evolution and role of eukaryotic-like domains in environmental intracellular bacteria: the battle with a eukaryotic cell

Intracellular pathogens that are able to thrive in different environments, such as Legionella spp. which preferentially live in protozoa in aquatic environments or environmental Chlamydiae which replicate either within protozoa or a range of animals, possess a plethora of cellular biology tools to influence their eukaryotic host. The host manipulation tools that evolved in the interaction with protozoa, confer these bacteria the capacity to also infect phylogenetically distinct eukaryotic cells, such as macrophages and thus they can also be human pathogens. To manipulate the host cell, bacteria use protein secretion systems and molecular effectors. Although these molecular effectors are encoded in bacteria, they are expressed and function in a eukaryotic context often mimicking or inhibiting eukaryotic proteins. Indeed, many of these effectors have eukaryotic-like domains. In this review we propose that the main pathways environmental intracellular bacteria need to subvert in order to establish the host eukaryotic cell as a replication niche are chromatin remodelling, ubiquitination signalling, and modulation of protein-protein interactions via tandem repeat domains. We then provide mechanistic insight into how these proteins might have evolved as molecular weapons. Finally, we highlight that in environmental intracellular bacteria the number of eukaryotic-like domains and proteins is considerably higher than in intracellular bacteria specialised to an isolated niche, such as obligate intracellular human pathogens. As mimics of eukaryotic proteins are critical components of host pathogen interactions, this distribution of eukaryotic-like domains suggests that the environment has selected them.

The life cycle-dependent transcriptional profile of the obligate intracellular amoeba symbiont Amoebophilus asiaticus

Free-living amoebae often harbor obligate intracellular bacterial symbionts. Amoebophilus (A.) asiaticus is a representative of a lineage of amoeba symbionts in the phylum Bacteroidota. Here, we analyse the transcriptome of A. asiaticus strain 5a2 at four time points during its infection cycle and replication within the Acanthamoeba host using RNA sequencing. Our results reveal a dynamic transcriptional landscape throughout different A. asiaticus life cycle stages. Many intracellular bacteria and pathogens utilize eukaryotic-like proteins (ELPs) for host cell interaction and the A. asiaticus 5a2 genome shows a particularly high abundance of ELPs. We show the expression of all genes encoding ELPs and found many ELPs to be differentially expressed. At the replicative stage of A. asiaticus, ankyrin repeat proteins and tetratricopeptide/Sel1-like repeat proteins were upregulated. At the later time points, high expression levels of a type 6 secretion system that likely prepares for a new infection cycle after lysing its host, were found. This study reveals comprehensive insights into the intracellular lifestyle of A. asiaticus and highlights candidate genes for host cell interaction. The results from this study have implications for other intracellular bacteria such as other amoeba-associated bacteria and the arthropod symbionts Cardinium forming the sister lineage of A. asiaticus.

Conservation lost: host-pathogen battles drive diversification and expansion of gene families

One of the strongest drivers in evolution is the struggle to survive a host-pathogen battle. This pressure selects for diversity among the factors directly involved in this battle, including virulence factors deployed by pathogens, their corresponding host targets, and host immune factors. A logical outcome of this diversification is that over time, the sequence of many immune factors will not be evolutionarily conserved across a broad range of species. Thus, while universal sequence conservation is often hailed as the hallmark of the importance of a particular gene, the immune system does not necessarily play by these rules when defending against co-evolving pathogens. This loss of sequence conservation is in contrast to many signaling pathways in development and basic cell biology that are not targeted by pathogens. In addition to diversification, another consequence of host-pathogen battles can be an amplification in gene number, thus leading to large gene families that have sequence relatively specific to a particular strain, species, or clade. Here we highlight this general theme across a variety of pathogen virulence factors and host immune factors. We summarize the wide range and number across species of these expanded, lineage-specific host-pathogen factors including ubiquitin ligases, nucleotide-binding leucine-rich repeat receptors, GTPases, and proteins without obvious biochemical function but that nonetheless play key roles in immunity.

Chlamydiae in the Environment

Chlamydiae have been known for more than a century as major pathogens of humans. Yet they are also found ubiquitously in the environment where they thrive within protists and in an unmatched wide range of animals. This review summarizes recent advances in understanding chlamydial diversity and distribution in nature. Studying these environmental chlamydiae provides a novel perspective on basic chlamydial biology and evolution. A picture is beginning to emerge with chlamydiae representing one of the evolutionarily most ancient and successful groups of obligate intracellular bacteria.

Draft Genome Sequences of Chlamydiales Bacterium STE3 and Neochlamydia sp. Strain AcF84, Endosymbionts of Acanthamoeba spp

Chlamydiales bacterium STE3 and Neochlamydia sp. strain AcF84 are obligate intracellular symbionts of Acanthamoeba spp. isolated from the biofilm of a littoral cave wall and gills from striped tiger leaf fish, respectively. We report the draft genome sequences of these two environmental chlamydiae affiliated with the family Parachlamydiaceae.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

The Hidden Genomics of Chlamydia trachomatis

The application of whole-genome sequencing has moved us on from sequencing single genomes to defining unravelling population structures in different niches, and at the -species, -serotype or even -genus level, and in local, national and global settings. This has been instrumental in cataloguing and revealing a huge a range of diversity in this bacterium, when at first we thought there was little. Genomics has challenged assumptions, added insight, as well as confusion and glimpses of truths. What is clear is that at a time when we start to realise the extent and nature of the diversity contained within a genus or a species like this, the huge depth of knowledge communities have developed, through cell biology, as well as the new found molecular approaches will be more precious than ever to link genotype to phenotype. Here we detail the technological developments and insights we have seen during the relatively short time since we began to see the hidden genome of Chlamydia trachomatis.

Symbiont-Mediated Defense against Legionella pneumophila in Amoebae

Bacterial pathogens are generally investigated in the context of disease. To prevent outbreaks, it is essential to understand their lifestyle and interactions with other microbes in their natural environment. Legionella pneumophila is an important human respiratory pathogen that survives and multiplies in biofilms or intracellularly within protists, such as amoebae. Importantly, transmission to humans occurs from these environmental sources. Legionella infection generally leads to rapid host cell lysis. It was therefore surprising to observe that amoebae, including fresh environmental isolates, were well protected during Legionella infection when the bacterial symbiont Protochlamydia amoebophila was also present. Legionella was not prevented from invading amoebae but was impeded in its ability to develop fully virulent progeny and were ultimately cleared in the presence of the symbiont. This study highlights how ecology and virulence of an important human pathogen is affected by a defensive amoeba symbiont, with possibly major consequences for public health.

Legionella pneumophila is an important opportunistic pathogen for which environmental reservoirs are crucial for the infection of humans. In the environment, free-living amoebae represent key hosts providing nutrients and shelter for highly efficient intracellular proliferation of L. pneumophila, which eventually leads to lysis of the protist. However, the significance of other bacterial players for L. pneumophila ecology is poorly understood. In this study, we used a ubiquitous amoeba and bacterial endosymbiont to investigate the impact of this common association on L. pneumophila infection. We demonstrate that L. pneumophila proliferation was severely suppressed in Acanthamoeba castellanii harboring the chlamydial symbiont Protochlamydia amoebophila. The amoebae survived the infection and were able to resume growth. Different environmental amoeba isolates containing the symbiont were equally well protected as different L. pneumophila isolates were diminished, suggesting ecological relevance of this symbiont-mediated defense. Furthermore, protection was not mediated by impaired L. pneumophila uptake. Instead, we observed reduced virulence of L. pneumophila released from symbiont-containing amoebae. Pronounced gene expression changes in the presence of the symbiont indicate that interference with the transition to the transmissive phase impedes the L. pneumophila infection. Finally, our data show that the defensive response of amoebae harboring P. amoebophila leaves the amoebae with superior fitness reminiscent of immunological memory. Given that mutualistic associations between bacteria and amoebae are widely distributed, P. amoebophila and potentially other amoeba endosymbionts could be key in shaping environmental survival, abundance, and virulence of this important pathogen, thereby affecting the frequency of human infection.

The draft genome of strain cCpun from biting midges confirms insect Cardinium are not a monophyletic group and reveals a novel gene family expansion in a symbiont

BACKGROUND

It is estimated that 13% of arthropod species carry the heritable symbiont Cardinium hertigii. 16S rRNA and gyrB sequence divides this species into at least four groups (A-D), with the A group infecting a range of arthropods, the B group infecting nematode worms, the C group infecting Culicoides biting midges, and the D group associated with the marine copepod Nitocra spinipes. To date, genome sequence has only been available for strains from groups A and B, impeding general understanding of the evolutionary history of the radiation. We present a draft genome sequence for a C group Cardinium, motivated both by the paucity of genomic information outside of the A and B group, and the importance of Culicoides biting midge hosts as arbovirus vectors.

METHODS

We reconstructed the genome of cCpun, a Cardinium strain from group C that naturally infects Culicoides punctatus, through Illumina sequencing of infected host specimens.

RESULTS

The draft genome presented has high completeness, with BUSCO scores comparable to closed group A Cardinium genomes. Phylogenomic analysis based on concatenated single copy core proteins do not support Cardinium from arthropod hosts as a monophyletic group, with nematode Cardinium strains nested within the two groups infecting arthropod hosts. Analysis of the genome of cCpun revealed expansion of a variety of gene families classically considered important in symbiosis (e.g., ankyrin domain containing genes), and one set-characterized by DUF1703 domains-not previously associated with symbiotic lifestyle. This protein group encodes putative secreted nucleases, and the cCpun genome carried at least 25 widely divergent paralogs, 24 of which shared a common ancestor in the C group. The genome revealed no evidence in support of B vitamin provisioning to its haematophagous host, and indeed suggests Cardinium may be a net importer of biotin.

DISCUSSION

These data indicate strains of Cardinium within nematodes cluster within Cardinium strains found in insects. The draft genome of cCpun further produces new hypotheses as to the interaction of the symbiont with the midge host, in particular the biological role of DUF1703 nuclease proteins that are predicted as being secreted by cCpun. In contrast, the coding content of this genome provides no support for a role for the symbiont in provisioning the host with B vitamins.

Metagenomic Analysis of Fish-Associated Ca. Parilichlamydiaceae Reveals Striking Metabolic Similarities to the Terrestrial Chlamydiaceae

Chlamydiae are an example of obligate intracellular bacteria that possess highly reduced, compact genomes (1.0-3.5 Mbp), reflective of their abilities to sequester many essential nutrients from the host that they no longer need to synthesize themselves. The Chlamydiae is a phylum with a very wide host range spanning mammals, birds, fish, invertebrates, and unicellular protists. This ecological and phylogenetic diversity offers ongoing opportunities to study intracellular survival and metabolic pathways and adaptations. Of particular evolutionary significance are Chlamydiae from the recently proposed Ca. Parilichlamydiaceae, the earliest diverging clade in this phylum, species of which are found only in aquatic vertebrates. Gill extracts from three Chlamydiales-positive Australian aquaculture species (Yellowtail kingfish, Striped trumpeter, and Barramundi) were subject to DNA preparation to deplete host DNA and enrich microbial DNA, prior to metagenome sequencing. We assembled chlamydial genomes corresponding to three Ca. Parilichlamydiaceae species from gill metagenomes, and conducted functional genomics comparisons with diverse members of the phylum. This revealed highly reduced genomes more similar in size to the terrestrial Chlamydiaceae, standing in contrast to members of the Chlamydiae with a demonstrated cosmopolitan host range. We describe a reduction in genes encoding synthesis of nucleotides and amino acids, among other nutrients, and an enrichment of predicted transport proteins. Ca. Parilichlamydiaceae share 342 orthologs with other chlamydial families. We hypothesize that the genome reduction exhibited by Ca. Parilichlamydiaceae and Chlamydiaceae is an example of within-phylum convergent evolution. The factors driving these events remain to be elucidated.

Environmental Metagenomic Assemblies Reveal Seven New Highly Divergent Chlamydial Lineages and Hallmarks of a Conserved Intracellular Lifestyle

The Chlamydiae phylum exclusively encompasses bacteria sharing a similar obligate intracellular life cycle. Existing 16S rDNA data support a high diversity within the phylum, however genomic data remain scarce owing to the difficulty in isolating strains using culture systems with eukaryotic cells. Yet, Chlamydiae genome data extracted from large scale metagenomic studies might help fill this gap. This work compares 33 cultured and 27 environmental, uncultured chlamydial genomes, in order to clarify the phylogenetic relatedness of the new chlamydial clades and to investigate the genetic diversity of the Chlamydiae phylum. The analysis of published chlamydial genomes from metagenomics bins and single cell sequencing allowed the identification of seven new deeply branching chlamydial clades sharing genetic hallmarks of parasitic Chlamydiae. Comparative genomics suggests important biological differences between those clades, including loss of many proteins involved in cell division in the genus Similichlamydia, and loss of respiratory chain and tricarboxylic acid cycle in several species. Comparative analyses of chlamydial genomes with two proteobacterial orders, the Rhizobiales and the Rickettsiales showed that genomes of different Rhizobiales families are much more similar than genomes of different Rickettsiales families. On the other hand, the chlamydial 16S rRNAs exhibit a higher sequence conservation than their Rickettsiales counterparts, while chlamydial proteins exhibit increased sequence divergence. Studying the diversity and genome plasticity of the entire Chlamydiae phylum is of major interest to better understand the emergence and evolution of this ubiquitous and ancient clade of obligate intracellular bacteria.

Transcriptome Sequencing Reveals Novel Candidate Genes for Cardinium hertigii-Caused Cytoplasmic Incompatibility and Host-Cell Interaction

The majority of insects carry maternally inherited intracellular bacteria that are important in their hosts’ biology, ecology, and evolution. Some of these bacterial symbionts cause a reproductive failure known as cytoplasmic incompatibility (CI). In CI, the mating of symbiont-infected males and uninfected females produces few or no daughters. The CI symbiont then spreads and can have a significant impact on the insect host population. Cardinium, a bacterial endosymbiont of the parasitoid wasp Encarsia in the Bacteroidetes, is the only bacterial lineage known to cause CI outside the Alphaproteobacteria, where Wolbachia and another recently discovered CI symbiont reside. Here, we sought insight into the gene expression of a CI-inducing Cardinium strain in its natural host, Encarsia suzannae. Our study provides the first insights into the Cardinium transcriptome and provides support for the hypothesis that Wolbachia and Cardinium target similar host pathways with distinct and largely unrelated sets of genes.

Cytoplasmic incompatibility (CI) is an intriguing, widespread, symbiont-induced reproductive failure that decreases offspring production of arthropods through crossing incompatibility of infected males with uninfected females or with females infected with a distinct symbiont genotype. For years, the molecular mechanism of CI remained unknown. Recent genomic, proteomic, biochemical, and cell biological studies have contributed to understanding of CI in the alphaproteobacterium Wolbachia and implicate genes associated with the WO prophage. Besides a recently discovered additional lineage of alphaproteobacterial symbionts only moderately related to Wolbachia, Cardinium (Bacteroidetes) is the only other symbiont known to cause CI, and genomic evidence suggests that it has very little homology with Wolbachia and evolved this phenotype independently. Here, we present the first transcriptomic study of the CI Cardinium strain cEper1, in its natural host, Encarsia suzannae, to detect important CI candidates and genes involved in the insect-Cardinium symbiosis. Highly expressed transcripts included genes involved in manipulating ubiquitination, apoptosis, and host DNA. Female-biased genes encoding ribosomal proteins suggest an increase in general translational activity of Cardinium in female wasps. The results confirm previous genomic analyses that indicated that Wolbachia and Cardinium utilize different genes to induce CI, and transcriptome patterns further highlight expression of some common pathways that these bacteria use to interact with the host and potentially cause this enigmatic and fundamental manipulation of host reproduction.

IMPORTANCE The majority of insects carry maternally inherited intracellular bacteria that are important in their hosts’ biology, ecology, and evolution. Some of these bacterial symbionts cause a reproductive failure known as cytoplasmic incompatibility (CI). In CI, the mating of symbiont-infected males and uninfected females produces few or no daughters. The CI symbiont then spreads and can have a significant impact on the insect host population. Cardinium, a bacterial endosymbiont of the parasitoid wasp Encarsia in the Bacteroidetes, is the only bacterial lineage known to cause CI outside the Alphaproteobacteria, where Wolbachia and another recently discovered CI symbiont reside. Here, we sought insight into the gene expression of a CI-inducing Cardinium strain in its natural host, Encarsia suzannae. Our study provides the first insights into the Cardinium transcriptome and provides support for the hypothesis that Wolbachia and Cardinium target similar host pathways with distinct and largely unrelated sets of genes.

Unexpected genomic features in widespread intracellular bacteria: evidence for motility of marine chlamydiae

Chlamydiae are obligate intracellular bacteria comprising important human pathogens and symbionts of protists. Molecular evidence indicates a tremendous diversity of chlamydiae particularly in marine environments, yet our current knowledge is based mainly on terrestrial representatives. Here we provide first insights into the biology of marine chlamydiae representing three divergent clades. Our analysis of single-cell amplified genomes revealed hallmarks of the chlamydial lifestyle, supporting the ancient origin of their characteristic developmental cycle and major virulence mechanisms. Surprisingly, these chlamydial genomes encode a complete flagellar apparatus, a previously unreported feature. We show that flagella are an ancient trait that was subject to differential gene loss among extant chlamydiae. Together with a chemotaxis system, these marine chlamydiae are likely motile, with flagella potentially playing a role during host cell infection. This study broadens our view on chlamydial biology and indicates a largely underestimated potential to adapt to different hosts and environments.

Biphasic Metabolism and Host Interaction of a Chlamydial Symbiont

Chlamydiae are obligate intracellular bacteria comprising well-known human pathogens and ubiquitous symbionts of protists, which are characterized by a unique developmental cycle. Here we comprehensively analyzed gene expression dynamics of Protochlamydia amoebophila during infection of its Acanthamoeba host by RNA sequencing. This revealed a highly dynamic transcriptional landscape, where major transcriptional shifts are conserved among chlamydial symbionts and pathogens. Our data served to propose a time-resolved model for type III protein secretion during the developmental cycle, and we provide evidence for a biphasic metabolism of P. amoebophila during infection, which involves energy parasitism and amino acids as the carbon source during initial stages and a postreplicative switch to endogenous glucose-based ATP production. This fits well with major transcriptional changes in the amoeba host, where upregulation of complex sugar breakdown precedes the P. amoebophila metabolic switch. The biphasic chlamydial metabolism represents a unique adaptation to exploit eukaryotic host cells, which likely contributed to the evolutionary success of this group of microbes. IMPORTANCE Chlamydiae are known as major bacterial pathogens of humans, causing the ancient disease trachoma, but they are also frequently found in the environment where they infect ubiquitous protists such as amoebae. All known chlamydiae require a eukaryotic host cell to thrive. Using the environmental chlamydia Protochlamydia amoebophila within its natural host, Acanthamoeba castellanii, we investigated gene expression dynamics in vivo and throughout the complete chlamydial developmental cycle for the first time. This allowed us to infer how a major virulence mechanism, the type III secretion system, is regulated and employed, and we show that the physiology of chlamydiae undergoes a complete shift regarding carbon metabolism and energy generation. This study provides comprehensive insights into the infection strategy of chlamydiae and reveals a unique adaptation to life within a eukaryotic host cell.

CRISPR System Acquisition and Evolution of an Obligate Intracellular Chlamydia-Related Bacterium

Recently, a new Chlamydia-related organism, Protochlamydia naegleriophila KNic, was discovered within a Naegleria amoeba. To decipher the mechanisms at play in the modeling of genomes from the Protochlamydia genus, we sequenced the full genome of Pr. naegleriophila, which includes a 2,885,090 bp chromosome and a 145,285 bp megaplasmid. For the first time within the Chlamydiales order, we describe the presence of a clustered regularly interspaced short palindromic repeats (CRISPR) system, the immune system of bacteria, located on the chromosome. It is composed of a small CRISPR locus comprising eight repeats and associated cas-cse genes of the subtype I-E. A CRISPR locus is also present within Chlamydia sp. Diamant, another Pr. naegleriophila strain, suggesting that the CRISPR system was acquired by a common ancestor of Pr. naegleriophila, after its divergence from Pr. amoebophila. Both nucleotide bias and comparative genomics approaches identified probable horizontal gene acquisitions within two and four genomic islands in Pr. naegleriophila KNic and Diamant genomes, respectively. The plasmid encodes an F-type conjugative system highly similar to 1) that found in the Pam100G genomic island of Pr. amoebophila UWE25 chromosome, as well as on the plasmid of Rubidus massiliensis and 2) to the three genes remaining in the chromosome of Parachlamydia acanthamoebae strains. Therefore, this conjugative system was likely acquired on an ancestral plasmid before the divergence of Parachlamydiaceae Overall, this new complete Pr. naegleriophila genome sequence enables further investigation of the dynamic processes shaping the genomes of the family Parachlamydiaceae and the genus Protochlamydia.

lolal Is an Evolutionarily New Epigenetic Regulator of dpp Transcription during Dorsal-Ventral Axis Formation

Secreted ligands in the Dpp/BMP family drive dorsal–ventral (D/V) axis formation in all Bilaterian species. However, maternal factors regulating Dpp/BMP transcription in this process are largely unknown. We identified the BTB domain protein longitudinals lacking-like (lolal) as a modifier of decapentaplegic (dpp) mutations. We show that Lolal is evolutionarily related to the Trithorax group of chromatin regulators and that lolal interacts genetically with the epigenetic factor Trithorax-like during Dpp D/V signaling. Maternally driven Lolal^HA is found in oocytes and translocates to zygotic nuclei prior to the point at which dpp transcription begins. lolal maternal and zygotic mutant embryos display significant reductions in dpp, pMad, and zerknullt expression, but they are never absent. The data suggest that lolal is required to maintain dpp transcription during D/V patterning. Phylogenetic data revealed that lolal is an evolutionarily new gene present only in insects and crustaceans. We conclude that Lolal is the first maternal protein identified with a role in dpp D/V transcriptional maintenance, that Lolal and the epigenetic protein Trithorax-like are essential for Dpp D/V signaling and that the architecture of the Dpp D/V pathway evolved in the arthropod lineage after the separation from vertebrates via the incorporation of new genes such as lolal.

Host-Associated Genomic Features of the Novel Uncultured Intracellular Pathogen Ca. Ichthyocystis Revealed by Direct Sequencing of Epitheliocysts

Advances in single-cell and mini-metagenome sequencing have enabled important investigations into uncultured bacteria. In this study, we applied the mini-metagenome sequencing method to assemble genome drafts of the uncultured causative agents of epitheliocystis, an emerging infectious disease in the Mediterranean aquaculture species gilthead seabream. We sequenced multiple cyst samples and constructed 11 genome drafts from a novel beta-proteobacterial lineage, Candidatus Ichthyocystis. The draft genomes demonstrate features typical of pathogenic bacteria with an obligate intracellular lifestyle: a reduced genome of up to 2.6 Mb, reduced G + C content, and reduced metabolic capacity. Reconstruction of metabolic pathways reveals that Ca Ichthyocystis genomes lack all amino acid synthesis pathways, compelling them to scavenge from the fish host. All genomes encode type II, III, and IV secretion systems, a large repertoire of predicted effectors, and a type IV pilus. These are all considered to be virulence factors, required for adherence, invasion, and host manipulation. However, no evidence of lipopolysaccharide synthesis could be found. Beyond the core functions shared within the genus, alignments showed distinction into different species, characterized by alternative large gene families. These comprise up to a third of each genome, appear to have arisen through duplication and diversification, encode many effector proteins, and are seemingly critical for virulence. Thus, Ca Ichthyocystis represents a novel obligatory intracellular pathogenic beta-proteobacterial lineage. The methods used: mini-metagenome analysis and manual annotation, have generated important insights into the lifestyle and evolution of the novel, uncultured pathogens, elucidating many putative virulence factors including an unprecedented array of novel gene families.

Emerging pathogens of gilthead seabream: characterisation and genomic analysis of novel intracellular β-proteobacteria

New and emerging environmental pathogens pose some of the greatest threats to modern aquaculture, a critical source of food protein globally. As with other intensive farming practices, increasing our understanding of the biology of infections is important to improve animal welfare and husbandry. The gill infection epitheliocystis is increasingly problematic in gilthead seabream (Sparus aurata), a major Mediterranean aquaculture species. Epitheliocystis is generally associated with chlamydial bacteria, yet we were not able to localise chlamydial targets within the major gilthead seabream lesions. Two previously unidentified species within a novel β-proteobacterial genus were instead identified. These co-infecting intracellular bacteria have been characterised using high-resolution imaging and genomics, presenting the most comprehensive study on epitheliocystis agents to date. Draft genomes of the two uncultured species, Ca. Ichthyocystis hellenicum and Ca. Ichthyocystis sparus, have been de novo sequenced and annotated from preserved material. Analysis of the genomes shows a compact core indicating a metabolic dependency on the host, and an accessory genome with an unprecedented number of tandemly arrayed gene families. This study represents a critical insight into novel, emerging fish pathogens and will be used to underpin future investigations into the bacterial origins, and to develop diagnostic and treatment strategies.

Whole genome sequencing revealed host adaptation-focused genomic plasticity of pathogenic Leptospira

Leptospirosis, caused by pathogenic Leptospira spp., has recently been recognized as an emerging infectious disease worldwide. Despite its severity and global importance, knowledge about the molecular pathogenesis and virulence evolution of Leptospira spp. remains limited. Here we sequenced and analyzed 102 isolates representing global sources. A high genomic variability were observed among different Leptospira species, which was attributed to massive gene gain and loss events allowing for adaptation to specific niche conditions and changing host environments. Horizontal gene transfer and gene duplication allowed the stepwise acquisition of virulence factors in pathogenic Leptospira evolved from a recent common ancestor. More importantly, the abundant expansion of specific virulence-related protein families, such as metalloproteases-associated paralogs, were exclusively identified in pathogenic species, reflecting the importance of these protein families in the pathogenesis of leptospirosis. Our observations also indicated that positive selection played a crucial role on this bacteria adaptation to hosts. These novel findings may lead to greater understanding of the global diversity and virulence evolution of Leptospira spp.

Chlamydial seasonal dynamics and isolation of 'Candidatus Neptunochlamydia vexilliferae' from a Tyrrhenian coastal lake

The Chlamydiae are a phylum of obligate intracellular bacteria comprising important human and animal pathogens, yet their occurrence in the environment, their phylogenetic diversity and their host range has been largely underestimated. We investigated the seasonality of environmental chlamydiae in a Tyrrhenian coastal lake. By catalysed reporter deposition fluorescence in situ hybridization, we quantified the small planktonic cells and detected a peak in the abundance of environmental chlamydiae in early autumn with up to 5.9 × 10(4) cells ml(-1) . Super-resolution microscopy improved the visualization and quantification of these bacteria and enabled the detection of pleomorphic chlamydial cells in their protist host directly in an environmental sample. To isolate environmental chlamydiae together with their host, we applied a high-throughput limited dilution approach and successfully recovered a Vexillifera sp., strain harbouring chlamydiae (93% 16S rRNA sequence identity to Simkania negevensis), tentatively named 'Candidatus Neptunochlamydia vexilliferae'. Transmission electron microscopy in combination with fluorescence in situ hybridization was used to prove the intracellular location of these bacteria representing the first strain of marine chlamydiae stably maintained alongside with their host in a laboratory culture. Taken together, this study contributes to a better understanding of the distribution and diversity of environmental chlamydiae in previously neglected marine environments.

Orientia tsutsugamushi Strain Ikeda Ankyrin Repeat-Containing Proteins Recruit SCF1 Ubiquitin Ligase Machinery via Poxvirus-Like F-Box Motifs

A rising theme among intracellular microbes is the delivery of ankyrin repeat-containing effectors (Anks) that interact with target proteins to co-opt host cell functions. Orientia tsutsugamushi, an obligate intracellular bacterium and the etiologic agent of scrub typhus, encodes one of the largest Ank repertoires of any sequenced microorganism. They have been previously identified as type 1 secretion system substrates. Here, in silico and manual sequence analyses revealed that a large proportion of O. tsutsugamushi strain Ikeda Anks bear a eukaryotic/poxvirus-like F-box motif, which is known to recruit host cell SCF1 ubiquitin ligase machinery. We assessed the Anks for the ability to serve as F-box proteins. Coimmunoprecipitation assays demonstrated that F-box-containing Anks interact with overexpressed and/or endogenous SCF1 components. When coexpressed with FLAG-Ank4_01 or FLAG-Ank9, a glutathione S-transferase (GST)-tagged version of the SCF1 component SKP1 localized to subcellular sites of FLAG-Ank accumulation. The abilities of recombinant Anks to interact and colocalize with SKP1 were F-box dependent. GST-SKP1 precipitated O. tsutsugamushi-derived Ank9 from infected host cells, verifying both that the pathogen expresses Ank9 during infection and the protein's capability to bind SKP1. Aligning O. tsutsugamushi, poxviral, and eukaryotic F-box sequences delineated three F-box residues that are highly conserved and likely to be functionally important. Substitution of these residues ablated the ability of GFP-Ank9 to interact with GST-SKP1. These results demonstrate that O. tsutsugamushi strain Ikeda Anks can co-opt host cell polyubiquitination machinery, provide the first evidence that an O. tsutsugamushi Ank does so during infection, and advance overall understanding of microbial F-box proteins.

IMPORTANCE

Ankyrin repeat-containing proteins (Anks) are important virulence factors of intracellular bacteria that mediate protein-protein interactions with host cell targets. Orientia tsutsugamushi, which causes a debilitating infection called scrub typhus in one of the most densely populated regions of the world, encodes one of the largest Ank armamentariums of any sequenced bacterium. This study demonstrates that O. tsutsugamushi strain Ikeda Anks also bear F-box motifs that interact with host cell polyubiquitination machinery. By proving that an Orientia-derived Ank interacts with SKP1 in infected cells, this evidences the first bona fide Orientia effector and the first example of an endogenous F-box-containing Ank-mammalian-host ligand interaction for any intracellular bacterium. Also, importantly, this work identifies key residues that are essential for microbial F-box function.

A Rickettsiales symbiont of amoebae with ancient features

The Rickettsiae comprise intracellular bacterial symbionts and pathogens infecting diverse eukaryotes. Here, we provide a detailed characterization of 'Candidatus Jidaibacter acanthamoeba', a rickettsial symbiont of Acanthamoeba. The bacterium establishes the infection in its amoeba host within 2 h where it replicates within vacuoles. Higher bacterial loads and accelerated spread of infection at elevated temperatures were observed. The infection had a negative impact on host growth rate, although no increased levels of host cell lysis were seen. Phylogenomic analysis identified this bacterium as member of the Midichloriaceae. Its 2.4 Mb genome represents the largest among Rickettsiales and is characterized by a moderate degree of pseudogenization and a high coding density. We found an unusually large number of genes encoding proteins with eukaryotic-like domains such as ankyrins, leucine-rich repeats and tetratricopeptide repeats, which likely function in host interaction. There are a total of three divergent, independently acquired type IV secretion systems, and 35 flagellar genes representing the most complete set found in an obligate intracellular Alphaproteobacterium. The deeply branching phylogenetic position of 'Candidatus Jidaibacter acanthamoeba' together with its ancient features place it closely to the rickettsial ancestor and helps to better understand the transition from a free-living to an intracellular lifestyle.

Taxogenomics of the order Chlamydiales

Bacterial classification is a long-standing problem for taxonomists and species definition itself is constantly debated among specialists. The classification of strict intracellular bacteria such as members of the order Chlamydiales mainly relies on DNA- or protein-based phylogenetic reconstructions because these organisms exhibit few phenotypic differences and are difficult to culture. The availability of full genome sequences allows the comparison of the performance of conserved protein sequences to reconstruct Chlamydiales phylogeny. This approach permits the identification of markers that maximize the phylogenetic signal and the robustness of the inferred tree. In this study, a set of 424 core proteins was identified and concatenated to reconstruct a reference species tree. Although individual protein trees present variable topologies, we detected only few cases of incongruence with the reference species tree, which were due to horizontal gene transfers. Detailed analysis of the phylogenetic information of individual protein sequences (i) showed that phylogenies based on single randomly chosen core proteins are not reliable and (ii) led to the identification of twenty taxonomically highly reliable proteins, allowing the reconstruction of a robust tree close to the reference species tree. We recommend using these protein sequences to precisely classify newly discovered isolates at the family, genus and species levels.