16S rRNA Gene Sequencing-Based Identification and Comparative Analysis of the Fecal Microbiota of Five Syntopic Lizard Species from a Low-Mountain Area in Western Bulgaria

Abstract

Studies on the gut microbiome of free-living reptiles in Europe are generally fragmentary and still missing in Bulgaria. We aimed to identify and compare the fecal microbiota profiles of five syntopic lizard species from three families: the European green lizard (Lacerta viridis), the common wall lizard (Podarcis muralis), the meadow lizard (Darevskia praticola) (Lacertidae), the European snake-eyed skink (Ablepharus kitaibelii) (Scincidae), and the European slow worm (Anguis fragilis) (Anguidae), which coinhabit a low mountainous area in the western part of the country. A high-throughput sequencing of the hypervariable V3-V4 region of the 16S rRNA gene, performed on the Illumina HiSeq2500 platform, was used. The core microbiota of lizard hosts seems to be species-specific. A dynamic phyla proportion between hosts was found. The richest alpha diversity was observed in D. praticola, and the lowest alpha diversity was observed in P. muralis and A. fragilis. Within the three lacertids, the microbiota of D. praticola and L. viridis were more closely related to each other than they were to those of P. muralis. Sharing a largely common trophic resource (all species except A. fragilis are mainly insectivorous) was not an indication of similarity in their gut microbial communities.