Metagenomics uncovers dietary adaptations for chitin digestion in the gut microbiota of convergent myrmecophagous mammals

<strong>Metagenomics uncovers dietary adaptations for chitin digestion in the gut microbiota of convergent myrmecophagous mammals</strong> Sophie Teullet^a,#, Marie-Ka Tilak^a, Amandine Magdeleine^a, Roxane Schaub^b,c, Nora M. Weyer^d, Wendy Panaino^d,e, Andrea Fuller^d, W. J. Loughry^f, Nico L. Avenant^g, Benoit de Thoisy^h,i, Guillaume Borrel^j and Frédéric Delsuc^a,# ^aInstitut des Sciences de l’Evolution de Montpellier (ISEM), Univ Montpellier, CNRS, IRD, Montpellier, France ^bCIC AG/Inserm 1424, Centre Hospitalier de Cayenne Andrée Rosemon, Cayenne, French Guiana ^cTropical Biome and immunopathology, Université de Guyane, Labex CEBA, DFR Santé, Cayenne, French Guiana ^dBrain Function Research Group, School of Physiology, University of the Witwatersrand, Johannesburg, South Africa ^eCentre for African Ecology, School of Animals, Plant, and Environmental Sciences, University of the Witwatersrand, Johannesburg, South Africa ^fDepartment of Biology, Valdosta State University, Valdosta, GA, USA ^gNational Museum and Centre for Environmental Management, University of the Free State, Bloemfontein, South Africa ^hInstitut Pasteur de la Guyane, Cayenne, French Guiana, France ⁱKwata NGO, Cayenne, French Guiana, France ^jInstitut Pasteur, Université Paris Cité, UMR CNRS 6047, Evolutionary Biology of the Microbial Cell, Paris, France ^#Corresponding authors: [e-mail]; [e-mail] <em><strong>Abstract</strong></em> In mammals, myrmecophagy (ant and termite consumption) represents a striking example of dietary convergence. This trait evolved independently at least five times in placentals with myrmecophagous species comprising aardvarks, anteaters, some armadillos, pangolins, and aardwolves. The gut microbiome plays an important role in dietary adaptation, and previous analyses of 16S rRNA metabarcoding data have revealed convergence in the composition of the gut microbiota among some myrmecophagous species. However, the functions performed by these gut bacterial symbionts and their potential role in the digestion of prey chitinous exoskeletons remain open questions. Using long- and short-read sequencing of fecal samples, we generated 29 gut metagenomes from nine myrmecophagous and closely related insectivorous species sampled in French Guiana, South Africa, and the USA. From these, we reconstructed 314 high-quality bacterial genome bins of which 132 carried chitinase genes, highlighting their potential role in insect prey digestion. These chitinolytic bacteria belonged mainly to the family Lachnospiraceae, and some were likely convergently recruited in the different myrmecophagous species as they were detected in several host orders (i.e., <em>Enterococcus faecalis</em>, <em>Blautia</em> sp), suggesting that they could be directly involved in the adaptation to myrmecophagy. Others were found to be more host-specific, possibly reflecting phylogenetic constraints and environmental influences. Overall, our results highlight the potential role of the gut microbiome in chitin digestion in myrmecophagous mammals and provide the basis for future comparative studies performed at the mammalian scale to further unravel the mechanisms underlying the convergent adaptation to myrmecophagy. <em><strong>Figures and Tables </strong></em> <strong>Main_figures.zip</strong> <strong>Figure. 1. Phylogenetic position of the 314 high-quality selected bins reconstructed from 29 gut metagenomes of the nine focal myrmecophagous species within a reference prokaryotic phylogeny.</strong> A: Phylogeny of the 314 selected bins (red branches) with 2496 prokaryote reference genomes. Circles respectively indicate (from inner to outer circles): the bacterial phyla and kingdom to which these genome bins were assigned based on the Genome Taxonomy Database (Chaumeil <em>et al</em>, 2020). Clades, where a subtree was defined, are highlighted in blue for the Firmicutes (Fig 1B), green for the Bacteroidetes, and pink for the Proteobacteria (Figs S1 A and B respectively). B: Subtree within Firmicutes showing myrmecophagous-specific clades (blue highlights; dark blue corresponds to the three clades mentioned in the results, light blue to the other clades). The outer circle indicates the bacterial family to which these genome bins were assigned based on the Genome Taxonomy Database. <strong>Figure. 2. Phylogeny of the 394 GH18 sequences identified in 132 high-quality selected bins reconstructed from 29 gut metagenomes of the nine focal myrmecophagous species and relatives. </strong>Red branches indicate the 237 sequences having an active chitinolytic site (DXXDXDXE). Circles respectively indicate (from inner to outer circles): the bacterial family and phyla of the bin the sequence was retrieved from. Colored sequence names indicate the host species. Colored circles at certain nodes indicate enzymes to which sequences are similar when blasting them against the NCBI non-redundant protein database. <strong>Figure 3. Detection of the 314 high-quality selected bins (lines) in the 29 gut metagenomes (columns) of the nine focal species.</strong> Each square indicates the detection of a bin in a sample as estimated by anvi’o v7 (Eren <em>et al</em>, 2021). Names of bins are indicated on the left with red indicating bins detected in at least one soil sample (detection &gt; 0.25) (Fig S4 and detection table available via Zenodo). Phylogenetic relationships of host species, distinguished by different color strips, are represented at the bottom of the graph. Columns on the right indicate (from left to right): the number of GH18 sequences identified in each bin (from 0 to 17), the bin’s taxonomic phylum, class, order, and family. The phylogeny of the 314 selected bins inferred with PhyloPhlAn v3.0.58 (Asnicar <em>et al</em>, 2020) is also represented on the right of the graph (see Fig S2). Silhouettes were downloaded from phylopic.org. <strong>Figure 4. Distribution of chitinolytic selected bins (red links) among the nine focal myrmecophagous species and relatives. </strong>Phylogenies of the 314 high-quality selected bins (Fig S2) and of the nine host species (downloaded from timetree.org) are represented respectively on the left and the right of the graph. Links illustrate, for each bin, in which host species the bin was detected (detection threshold &gt; 0.25). Red links indicate bins in which at least one GH18 sequence with an active chitinolytic site (DXXDXDXE) was found (chitinolytic bins). The size of the circles at the tips of the host phylogeny is proportional to the number of samples (n = 1 for <em>D. kap</em>; n = 2 for <em>D. nov</em>, <em>C. uni </em>and <em>M. tri</em>; n = 3 for <em>T. tet</em> and <em>O. afe</em>; n = 4 for <em>D. sp. nov FG</em>; n = 6 for <em>P. cri</em> and <em>S. tem</em>). Bins’ names are indicated at the tip of the bins’ phylogeny and main bacterial phyla are indicated by colored vertical bars. This graph was done with the cophylo R package within the phytools suite (Revell, 2012). Silhouettes downloaded from phylopic.org. <strong>Table 1. Detailed sample information for the 33 fecal samples collected.</strong> <em><strong>Supplementary Materials</strong></em> <strong>Supplementary_material_Teullet_etal_2023.pdf </strong>contains supplementary figures (S1-4) and tables (S2-3). <strong>Supplementary_results_Teullet_etal_2023.pdf</strong> includes a comparison of genome statistics of the selected bins reconstructed from the long-read vs the short-read datasets, a phylogeny of the set of selected bins before dereplication (n = 407) and a comparison of the distribution of shared and specific genome bins carrying GH18 among host orders. <strong>Supplementary_material_files.zip </strong>contains a file for each of the supplementary material figures (S1-4) and tables (S2-3). <strong>Supplementary_results_files.zip</strong> contains a file for each of the three figures presented in the supplementary results. <strong>Table S1. Raw results of the different analyses conducted on each gut metagenome to reconstruct high-quality genome bins from raw metagenomic data for each dataset</strong> (long- and short-reads). <em><strong>Zenodo supplementary files</strong></em> <em>Assemblies</em> <strong>Long-read_metagenomic_assemblies_polished.zip</strong> contains 31 long-read metagenomes assembled with metaFlye strain v2.9 and polished with short reads using Pilon v1.4, which were used for binning. <strong>Long-read_metagenomic_assemblies_not_polished.zip</strong> contains 33 long-read metagenomes assembled with metaFlye strain v2.9 before polishing. <strong>Short-read_metagenomic_assemblies.zip</strong> contains 31 short-read metagenomes assembled with metaSPAdes and MEGAHIT. <em>N.B</em>: Two samples were not sequenced using Illumina short reads (DASY_M1746 and DASY_VLD168), only long reads were generated and assembled for these two samples and are made available here. As these assemblies could not be polished, these samples were not included in downstream analyses. Two samples were highly contaminated by host reads (CAB_M3141 and MYR_M5295) and not used for downstream analyses. As they were still assembled with the other samples, the corresponding metagenomes are made available here. <em>Binning: genome bins and dereplication results</em> <strong>High-quality_selected_bins_dereplicated.zip</strong> contains the 314 high-quality selected bins (&gt;90% completion, &lt;5% redundancy) reconstructed from long- and short-read metagenomes with metaBAT2 and dereplicated with dRep at 98% ANI. <strong>metaBAT2_short-read_assemblies_bins.zip </strong>contains all bins reconstructed from the short-read assemblies with metaBAT2 (i.e, output of metaBAT2). <strong>metaBAT2_long-read_assemblies_bins.zip</strong> contains all bins reconstructed from the long-read polished assemblies with metaBAT2 (i.e, output of metaBAT2). <strong>Output_dRep_98ANI_407_bins_long-short-reads.zip</strong> contains the output of the dereplication analysis done on the set of 407 high-quality selected genome bins reconstructed from long- (n = 201) and short-read (n = 206; labeled "spad") metagenomes. It was done with dRep using default parameters. After this step, the final dataset included 314 high-quality non-redundant genome bins. This folder includes: The primary clustering of selected genome bins using the Mash algorithm with an ANI threshold of 90%. The secondary clustering of selected genome bins using the fastANI algorithm with an ANI threshold of 98%. The clustering score attributed to each genome bin during dereplication. * symbol indicates genomes chosen to be the representative genomes of their cluster. <em>Phylogenetic datasets</em> <strong>Fig1_supplementary_files_alignment_tree.zip</strong> contains the alignment of the concatenated markers and the final tree reconstructed by PhyloPhlAn v3.0.58 (Asnicar <em>et al</em>, 2020) for the 314 high-quality selected genome bins and the 2496 prokaryote reference genomes. <strong>Fig2_supplementary_files_alignment_tree.zip</strong> contains the alignment of the 394 GH18 sequences identified in 132 high-quality selected bins done with MAFFT v7.450 and the corresponding tree inferred with RAxML v8.2.11 within Geneious Prime 2022.0.2. <strong>Phylogeny_314_selected_bins_dRep.zip</strong> contains the alignment of the concatenated markers and the final tree reconstructed by PhyloPhlAn v3.0.58 (Asnicar <em>et al</em>, 2020) for the 314 high-quality selected and dereplicated genome bins. <strong>Phylogeny_407_selected_bins_nodRep.zip</strong> contains the alignment of the concatenated markers and the final tree reconstructed by PhyloPhlAn v3.0.58 (Asnicar <em>et al</em>, 2020) for the 407 high-quality selected genome bins before dereplication. <em>Bins across samples</em> <strong>detection_bins_across_gut_metagenomes.txt </strong>is a tab-delimited file containing the detection values inferred by anvi'o v7 (Eren <em>et al</em>, 2021) for the 314 high-quality selected bins across the 29 gut metagenomes from the nine focal myrmecophagous species. <strong>abundance_bins_gut_metagenomes.txt </strong>is a tab-delimited file containing the absolute abundance values inferred by anvi'o v7 (Eren <em>et al</em>, 2021) for the 314 high-quality selected bins across the 29 gut metagenomes from the nine focal myrmecophagous species. <strong>detection_bins_across_soil_samples.txt</strong> is a tab-delimited file containing the detection values inferred by anvi'o v7 (Eren <em>et al</em>, 2021) for the 140 high-quality selected bins reconstructed from the aardvark, ground pangolin and southern aardwolf gut metagenomes across the eight soil samples collected on sample sites in South Africa.