Idiosyncratic genome degradation in a bacterial endosymbiont of periodical cicadas

Abstract

When a free-living bacterium transitions to a host-beneficial endosymbiotic lifestyle, it almost invariably loses a large fraction of its genome [1, 2]. The resulting small genomes often become unusually stable in size, structure, and coding capacity [3-5]. Candidatus Hodgkinia cicadicola (Hodgkinia), a bacterial endosymbiont of cicadas, sometimes exemplifies this genomic stability. The Hodgkinia genome has remained completely co-linear in some cicadas diverged by tens of millions of years [6, 7]. But in the long-lived periodical cicada Magicicada tredecim, the Hodgkinia genome has split into dozens of tiny, gene-sparse genomic circles that sometimes reside in distinct Hodgkinia cells [8]. Previous data suggested that other Magicicada species harbor similarly complex Hodgkinia populations, but the timing, number of origins, and outcomes of the splitting process were unknown. Here, by sequencing Hodgkinia metagenomes from the remaining six Magicicada species and two sister species, we show that all Magicicada species harbor Hodgkinia populations of at least twenty genomic circles each. We find little synteny among the 256 Hodgkinia circles analyzed except between the most closely related species. Individual gene phylogenies show that Hodgkinia first split in the common ancestor of Magicicada and its closest known relatives, but that most splitting has occurred within Magicicada and has given rise to highly variable Hodgkinia gene dosages between cicada species. These data show that Hodgkinia genome degradation has proceeded down different paths in different Magicicada species, and support a model of genomic degradation that is stochastic in outcome and likely nonadaptive for the host. These patterns mirror the genomic instability seen in some mitochondria.