An effect of 16S rRNA intercistronic variability on coevolutionary analysis in symbiotic bacteria: Molecular phylogeny of Arsenophonus triatominarum

The genes of ribosomal RNA are the most popular and frequently used markers for bacterial phylogeny and reconstruction of insect–symbiont coevolution. In primary symbionts, such as Buchnera and Wigglesworthia, genome economization leads to the establishment of a single copy of these sequences. In phylogenetic studies, they provide sufficient information and yield phylogenetic trees congruent with host evolution. In contrast, other symbiotic lineages (e.g., the genus Arsenophonus) carry a higher number of rRNA copies in their genomes, which may have serious consequences for phylogenetic inference. In this study, we show that in Arsenophonus triatominarum the degree of heterogeneity can affect reconstruction of phylogenetic relationships and mask possible coevolution between the symbiont and its host. Phylogenetic arrangement of individual rRNA copies was used, together with a calculation of their divergence time, to demonstrate that the incongruent 16S rDNA trees and low nucleotide diversity in the secondary symbiont could be reconciled with the coevolutionary scenario.