Abstract
The association between anaerobic ciliates and methanogenic archaea has been recognized for over a century. Nevertheless, knowledge of these associations is limited to a few ciliate species, and so the identification of patterns of host-symbiont specificity has been largely speculative. In this study, we integrated microscopy and genetic identification to survey the methanogenic symbionts of 32 free-living anaerobic ciliate species, mainly from the order Metopida. Based on Sanger and Illumina sequencing of the 16S rRNA gene, our results show that a single methanogenic symbiont population, belonging to Methanobacterium, Methanoregula, or Methanocorpusculum, is dominant in each host strain. Moreover, the host's taxonomy (genus and above) and environment (i.e. endobiotic, marine/brackish, or freshwater) are linked with the methanogen identity at the genus level, demonstrating a strong specificity and fidelity in the association. We also established cultures containing artificially co-occurring anaerobic ciliate species harboring different methanogenic symbionts. This revealed that the host-methanogen relationship is stable over short timescales in cultures without evidence of methanogenic symbiont exchanges, although our intraspecific survey indicated that metopids also tend to replace their methanogens over longer evolutionary timescales. Therefore, anaerobic ciliates have adapted a mixed transmission mode to maintain and replace their methanogenic symbionts, allowing them to thrive in oxygen-depleted environments.