Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the “Genomes from Earth’s Microbiomes” initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

Highlights d Chlamydial plasmids coevolved with their bacterial hosts over a billion years d Recombination with extrachromosomal elements and viruses shaped plasmid gene content d Plasmid-mediated chromosomal gene mobilization and transfer drove genome evolution d Plasmids contributed to adaptation of chlamydiae to diverse eukaryotic hosts

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.

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 16S 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.

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well-known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analyzed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, ‘Candidatus Tisiphia’ (formerly Torix group Rickettsia), Cardinium, and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, ‘Ca. Tisiphia’ and Cardinium endosymbionts. We provide insights in the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, ‘Ca. Tisiphia’ and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on earth and show the usefulness of transposable elements as an evolutionary “contact-tracing” tool.Data summaryAll supporting data, code and protocols have been provided within the article or through supplementary data files. Seven supplementary figures and seven supplementary tables are available with the online version of this article. Sequencing data used in this study was generated and previously published by Hendrickx et al., 2021. Genome assemblies generated in this study have been deposited under the project PRJEB52003 at DDBJ/ENA/GenBank. The MAG of R. oedothoracis OV001 was deposited at DDBJ/ENA/GenBank under the sample SAMN28026840. The genome of ‘Candidatus Rhabdochlamydia oedothoracis W744×776’ was previously published by Halter et al., 2022 and is available at DDBJ/ENA/GenBank (accession: CP075587-CP075588). The collection of genomes and proteomes, all files for phylogenetic analyses including gene alignments, concatenated alignments, and tree files, and original output files of the HGT and SNP predictions used in this study are available at zenodo (https://doi.org/10.5281/zenodo.6362846).