Domesticated equine species and their derived hybrids differ in their fecal microbiota
Animal Microbiome volume 2, Article number: 8 (2020)
Compared to horses and ponies, donkeys have increased degradation of dietary fiber. The longer total mean retention time of feed in the donkey gut has been proposed to be the basis of this, because of the increased time available for feed to be acted upon by enzymes and the gut microbiota. However, differences in terms of microbial concentrations and/or community composition in the hindgut may also underpin the increased degradation of fiber in donkeys. Therefore, a study was conducted to assess if differences existed between the fecal microbiota of pony, donkey and hybrids derived from them (i.e. pony × donkey) when fed the same forage diet.
Fecal community composition of prokaryotes and anaerobic fungi significantly differed between equine types. The relative abundance of two bacterial genera was significantly higher in donkey compared to both pony and pony x donkey: Lachnoclostridium 10 and ‘probable genus 10’ from the Lachnospiraceae family. The relative abundance of Piromyces was significantly lower in donkey compared to pony × donkey, with pony not significantly differing from either of the other equine types. In contrast, the uncultivated genus SK3 was only found in donkey (4 of the 8 animals). The number of anaerobic fungal OTUs was also significantly higher in donkey than in the other two equine types, with no significant differences found between pony and pony × donkey. Equine types did not significantly differ with respect to prokaryotic alpha diversity, fecal dry matter content or fecal concentrations of bacteria, archaea and anaerobic fungi.
Donkey fecal microbiota differed from that of both pony and pony × donkey. These differences related to a higher relative abundance and diversity of taxa with known, or speculated, roles in plant material degradation. These findings are consistent with the previously reported increased fiber degradation in donkeys compared to ponies, and suggest that the hindgut microbiota plays a role. This offers novel opportunities for pony and pony × donkey to extract more energy from dietary fiber via microbial mediated strategies. This could potentially decrease the need for energy dense feeds which are a risk factor for gut-mediated disease.
Hindgut microbial fermentation of plant material in equines results in the generation of volatile fatty acids, which are a major source of energy for the equine host [1,2,3]. The improvement of fiber utilization would decrease the need for energy dense concentrate feeds, which have been linked to increased risk of colic and adverse effects on the hindgut microbiota [4, 5]. Furthermore, in regions where only limited feed sources and poor quality forage are available, the efficient use of available fibrous feed sources by the animal is important to be able to meet its energy requirements. Optimizing fiber utilization can, therefore, have a beneficial impact on the health and welfare of equines.
The world’s 116 million domesticated equines includes horses/ponies (52%), donkeys (40%) and mules/hinnies (8%) . Domesticated donkeys have evolved from ancestors that inhabited semi-arid and often mountainous environments which had sparse vegetation . This is in contrast to the grassland plains where horses evolved. As such, it is not surprising that as part of their evolution donkeys developed an improved ability to degrade fiber which enables them to survive better on lignin rich, low energy, fibrous plants .
Compared to horses/ponies, donkeys have increased dietary fiber digestibility and decreased feed intake [8, 9]. This is true regardless of dietary fiber content , although differences in dry matter (DM) digestion between species become more pronounced when diet quality decreases . The increased mean retention time (MRT) of feed particles in the donkey gut, compared to horses/ponies [8, 9], has been suggested to be the cause of this increased fiber digestibility. It has been hypothesized that this is due to the increased time available for the dietary material to be acted upon by enzymes and the gut microbiota . However, it is also possible that different host species select for/harbor different gut microbial communities, as has been previously shown in a study with two species of deer and a hybrid derived from them .
The equine hindgut microbiota is comprised of bacteria, anaerobic fungi, archaea, protozoa, and viruses. Of these, mainly bacteria have been studied to date . Anaerobic fungi and certain bacteria are the only hindgut microbes directly involved in fiber degradation. Archaea and viruses are unable to degrade dietary material, and protozoa make a limited contribution to fiber degradation in the horse hindgut . It has been previously reported that the total number of anaerobic bacteria and cellulolytic bacteria in the equine caecum does not significantly differ between ponies and donkeys fed the same diet . However, a donkey anaerobic fungal isolate of Piromyces citronii was found to be superior to that of a pony isolate of the same species in terms of degrading cellulose . Anaerobic fungi are powerful fiber degraders due to their combined invasive growth and broad array of highly effective plant degrading enzymes [16, 17]. In ruminants fed poor quality herbage or straw anaerobic fungi are known to increase feed dry matter digestibility by 7–9% , however, their contribution to feed degradation in the equine hindgut remains to be determined.
Whilst numerous studies have explored the composition of the bacterial community in the hindgut of horses and ponies [5, 12, 19,20,21,22,23,24,25,26,27,28,29,30,31], only limited studies have been performed on other equine species to date. Liu et al.  described the fecal microbiota of domesticated donkeys, and noted differences in the relative abundances of bacterial phyla compared to those previously reported in horses and ponies. The fecal microbiota of the Tibetan wild ass has also been studied , which is an equine species that is more closely related to the domesticated donkey than horse . Captive Tibetan wild ass was reported to have decreased prokaryotic alpha diversity and a different community composition relative to wild Tibetan wild ass, and authors speculated that this difference was mainly diet related . No studies to date have been performed on mules or hinnies despite their global population of 9.6 million , the majority of which are purpose bred working animals.
In this study, fecal samples from domesticated ponies (Equus caballus), donkeys (Equus africanus asinus) and their derived hybrids (i.e. pony × donkey) fed the same forage diet were characterized in order to assess if the microbiota present in the equine hindgut differs between domesticated equine species and their derived hybrids. The fecal prokaryotic and anaerobic fungal community composition was determined using barcoded amplicon sequencing, and microbial concentrations measured using quantitative PCR. It was hypothesized that the fecal microbiota of donkey differed from that of pony in relation to taxa involved in fiber degradation, with the hybrid being intermediate relative to donkey and pony.
The animals used in this study were all healthy adults that had no known history of any gut-mediated disease. The animals included ponies (n = 8), donkeys (n = 8) and mules/hinnies (n = 8), and each equine type was composed of three females and five males. As the parentage of the mules/hinnies was not known, they are referred to as pony x donkey in this study. Further details of the individual animals are given in Additional file 1: Table S1.
Fecal dry Matter Content & Microbial Concentrations
The percentage fecal DM content (average ± SD) of pony (18.2 ± 2.0), pony × donkey (19.0 ± 1.7) and donkey (20.3 ± 1.5) did not significantly differ between equine types (P = 0.07). Equine types were also not significantly different with respect to fecal concentrations of bacteria, archaea or anaerobic fungi when expressed on a fresh weight basis (Table 1). When analyzed on a DM basis, equine types also did not differ in fecal concentrations of bacteria (P = 0.54), archaea (P = 0.70) and anaerobic fungi (P = 0.14) (data not shown).
Prokaryotic community composition
Bacteria (97.1 ± 1.14% of the 16S rRNA gene sequences) were represented by 1203 different operational taxonomic units (OTUs). Of these, 319 OTUs were detected at least once in all three equine types, and 18 OTUs were detected in all animals. The 1203 bacterial OTUs could be summarized to 172 different genus-level phylogenetic groupings. The archaea (2.9 ± 1.14% of the 16S rRNA gene sequences) were represented by ten OTUs. Of these, four OTUs were detected at least once in all three equine types and one OTU was detected in all animals. The ten OTUs could be summarized to two different genus-level phylogenetic groupings. Of the 17 phyla detected in total, the following six were predominant (> 1%): Firmicutes, Bacteroidetes, Verrucomicrobia, Fibrobacteres, Spirochaetes and Euryarchaeota (Additional file 2: Figure S1).
We did not observe a significant difference between equine types (P = 0.43) with respect to the number of OTUs observed in pony (average ± SD; 265 ± 9.9), donkey (258 ± 12.0) and pony × donkey (267 ± 18.7). Similarly, the ‘Phylogenetic Diversity’ metric used to assess alpha diversity was also not significantly different between equine types (P = 0.30): pony (average ± SD; 17.5 ± 0.44), donkey (17.8 ± 0.30) and pony × donkey (17.5 ± 0.63).
In terms of beta diversity at the OTU level, the prokaryotic community composition in donkey separated from both pony and pony × donkey along the first axis in the unweighted UniFrac principal co-ordinate analysis (PCoA), whilst pony and pony × donkey did not differ (Fig. 1a). Using weighted UniFrac PCoA, the separation of donkey from pony and pony × donkey did not occur along one axis, instead separating along a diagonal line (Fig. 1b).
Redundancy analysis (RDA) using genus level groupings of the OTUs indicated that equine type significantly contributed to explaining the observed variation in fecal prokaryotic community composition (P < 0.01), accounting for 22.8% of the overall variation in the dataset. Donkey separated from both pony and pony × donkey along the horizontal canonical axis, which represented 17.6% of the variation in prokaryotic community composition (Fig. 2). Pony and pony × donkey separated from each other along the second canonical axis, which represented 5.3% of the variation in the prokaryotic composition data. Five genus level phylogenetic groupings were strongly positively associated with donkey. Three of these could be annotated to the genus level (Sarcina, Ruminococcaceae NK4A214 group and Lachnoclostridium 10), whilst the other two could only be annotated to the family (Lachnospiraceae) or order (‘uncultured rumen bacterium’ in the Pla4 lineage of the phylum Planctomycetes) level. One genus level phylogenetic grouping was clearly positively associated with pony, and was annotated as an uncultured bacterium belonging to the family ‘gir-aah93h0’ within the order Bacteroidales. No genus level phylogenetic grouping was positively associated with pony × donkey.
A Kruskal-Wallis test indicated that the relative abundance of two of the genus level groups significantly different between equine types: the genus Lachnoclostridium 10 (P = 0.03) and an uncultured genus called ‘probable genus 10’ from the Lachnospiraceae family (P = 0.01). Dunn’s Sidak post-hoc analysis showed that the relative abundance of Lachnoclostridium 10 was significantly higher in donkey (average percentage relative abundance ± SD, 1.43 ± 0.91) compared to both pony × donkey (0.05 ± 0.11, P < 0.001) and pony (0.09 ± 0.14, P < 0.01). We did not observe a significant difference between pony × donkey and pony (P = 0.98). Lachnoclostridium 10 was detected in all eight of the donkeys sampled, whereas for pony × donkey and pony it was detected in only two and three animals, respectively, at lower relative abundances compared to donkey. With ‘probable genus 10’, Dunn’s Sidak post-hoc analysis also showed that the relative abundance of this genus was significantly higher in donkey (0.58 ± 0.19) compared to pony × donkey (0.07 ± 0.12, P < 0.01) and pony (not detected, P < 0.0001). We did not observe a significant difference between pony × donkey and pony (P = 0.73). This genus was detected in all eight donkeys, whereas in the pony × donkey it was detected in only three of the animals at lower relative abundances compared to donkey.
Anaerobic fungal community composition
In the anaerobic fungal sequence data 72 OTUs were detected. Of these, 13 OTUs were detected at least once in all three equine types, and three OTUs were detected in all animals. The 72 OTUs could be summarized to five different genus level groups. The genus Caecomyces was most predominant in the anaerobic fungal community of all of the animals (Fig. 3). The number of OTUs observed was significantly different between equine types (P < 0.01), with donkeys (average ± SD; 17 ± 5.2) having more OTUs than both pony (10 ± 4.6) and pony × donkey (10 ± 2.97). However, the ‘Phylogenetic Diversity’ metric used to assess alpha diversity did not significantly differ between equine types (P = 0.27): donkey (average ± SD; 0.40 ± 0.186), pony (0.28 ± 0.177) and pony × donkey (0.31 ± 0.072).
In terms of beta diversity at the OTU level, weighted and unweighted UniFrac PCoA showed that pony × donkey anaerobic fungal community composition varied less between individuals than that of donkey and pony (Fig. 4). Separation of pony × donkey from donkey occurred along the first axis of the unweighted PCoA, but no separation of these two equine types from pony was seen (Fig. 4a). No obvious separation of the samples by equine type was observed in the weighted PCoA (Fig. 4b), presumably due to the predominance of Caecomyces in all animals.
RDA using genus level groupings of the OTUs showed that equine type significantly contributed to explaining the observed variation in the fecal anaerobic fungal community composition (P < 0.01), and accounted for 23.6% of the overall variation in the dataset. Caecomyces was positively associated with pony and pony × donkey. Piromyces and the uncultivated genus SK3 were positively associated with pony × donkey and donkey, respectively (Fig. 5). The relative abundance of genera that could not be annotated also seemed to be higher in donkey compared to pony and pony × donkey.
A Kruskal Wallis test showed that only the genus level groups SK3 (P = 0.03) and Piromyces (P = 0.05) were significantly different between equine types. SK3 was detected in only four of the eight donkeys sampled, and none of the animals from the equine types pony × donkey and pony. A Dunn’s Sidak post-hoc analysis indicated that the relative abundance of SK3 was significantly higher in donkey (average percentage relative abundance ± SD; 2.5 ± 3.45) compared to pony × donkey (P = 0.03) and pony (P = 0.03). In contrast, Piromyces was detected in all three equine types. Dunn’s Sidak post-hoc analysis showed that the relative abundance of Piromyces was significantly decreased in donkey (0.54 ± 0.88) compared to pony × donkey (3.42 ± 2.68) (P = 0.01), with pony (2.13 ± 2.99) not significantly differing from either donkey (P = 0.52) or pony × donkey (P = 0.26). In pony × donkey, Piromyces was detected in all eight animals, whereas in donkey and pony, it was detected in only three and four animals, respectively.
As fecal microbiota is mainly representative of the distal region of the equine hindgut , the findings of this study need to be interpreted with this in mind. Fecal microbial concentrations or fecal DM content did not differ between equine types. This new finding suggests that the reported increased fiber degradation in donkeys, compared to ponies, is unlikely to be due to higher concentrations of microbes in the hindgut.
The predominance of the phylum Firmicutes in the fecal microbiota of all equines observed in this study is in line with previous reports for horse/pony [23, 26,27,28, 31, 35,36,37], donkey  and captive Tibetan wild ass . In this study, the fecal prokaryotic community composition was shown to significantly differ between equine types, with donkey being most distinct from both pony and pony × donkey. This novel finding appears to be underpinned by two bacterial genera from the Lachnospiraceae family that had a higher relative abundance and prevalence in donkey compared to both pony and pony × donkey: Lachnoclostridium 10 and ‘probable genus 10’.
Within the SILVA database version 128, the genus Lachnoclostridium 10 contains only one characterized species, Lachnoclostridium phytofermentans. The type strain of this species, initially published as Clostridium phytofermentans , has been effectively (but not validly) published as the type species for the new genus Lachnoclostridium, which included a total of 30 different validly described species at the time of its definition in 2013 . L. phytofermentans is an obligately anaerobic fibrolytic bacterium that can ferment a wide range of plant polysaccharides . This ability appears to be due to the numerous and diverse range of glycosyl hydrolases encoded within its genome, many of which have been acquired by horizontal gene transfer . Whilst the L. phytofermentans type strain was isolated from forest soil, its optimum temperature for growth was reported to be 37 °C with growth observed at pH 6.0–9.0 . One of the other few cultivated members of the genus Lachnoclostridium 10 listed in the SILVA database is the rumen bacterium FE2016, which was isolated as part of a study that cultivated plant-attached rumen bacteria (GenBank Accession No. KF698008). As part of the rumen Hungate 1000 project, several ruminal isolates described as belonging to the genus Lachnoclostridium have also been genome-sequenced .
No characterized species exist within ‘probable genus 10’. However, within the SILVA database this genus does contain 14 cultured bacterial isolates. All of these isolates are of ruminal origin. Three of the isolates were obtained from studies that used modified culturing conditions to isolate novel rumen bacteria: CA43 , P18  and NK4A212 . Isolate NK4A212 was proposed as a new genus , whilst isolate CA43 was shown to have carboxymethylcellulase and xylanase activity . The remaining 11 isolates are all unpublished, however, ten of them were isolated from the same ‘plant-attached rumen bacteria’ study as the rumen bacterium FE2016.
The isolation of plant-attached rumen bacteria belonging to Lachnoclostridium 10 and ‘probable genus 10’, combined with the demonstrated activities of these genera in terms of their ability to degrade polymers present in plant fibers, indicates they are likely to play a role in ruminal fiber degradation. This is likely to also be true of fiber degradation occurring in the equine hindgut. If this is the case, then the higher relative abundances and prevalence of these two genera in donkey compared to pony may contribute to the previously reported increased ability of donkeys, relative to horses/ponies, to digest fiber [8, 9].
Like the prokaryotic community, another novel finding in this study was that the anaerobic fungal community composition was also different between equine types. A higher number of anaerobic fungal OTUs was found in donkey compared to both pony and pony × donkey. As anaerobic fungi can vary in terms of their growth rate, substrate preferences and fiber degrading activity [16, 45], a larger diversity of anaerobic fungi may enable donkeys to utilize fibrous plant material more effectively compared to equines with less diverse anaerobic fungal populations. As anaerobic fungal pure cultures have been reported to have 3–29 different OTUs , it is likely that this increased number of OTUs in donkey is representative of one or two extra anaerobic fungal species being present.
Only donkey contained the uncultivated genus SK3 [47, 48], although the genus was only present in half of the animals sampled. Anaerobic fungal sequences belonging to SK3 have been previously obtained from cow manure , as well as the rumen of sheep, cattle and deer . SK3 was not previously found in cultivation independent studies of the hindgut digesta of a pony  or the fecal material of five different types of domesticated and non-domesticated equines . As such, this appears to be the first report of this genus in equines, and indicates that it is not ruminant specific. SK3 is a sister group to the two bulbous genera Caecomyces and Cyllamyces , suggesting that SK3 may be a third genus that has this type of rhizomycelium.
Caecomyces predominated all of the animals in this study, and has previously been reported to occur in the pony caecum and equine feces [51, 52]. The uncultivated genera AL1 (=NG1) and AL3 (=NG3) were previously found to be predominant in the feces of different equine species , and AL1 has also been detected along the equine hindgut . However, neither of these genera were detected in this study. Cultivation based studies of anaerobic fungi from domesticated equines have most commonly resulted in the isolation of Piromyces [16, 53,54,55], although other genera have also been reported [52, 56]. In this study, Piromyces was detected in all three equine types as a minor genus, but was not detected in all animals. It has been previously shown that a Piromyces citronii isolate from a donkey degraded cellulose more rapidly and to a greater extent than a P. citronii isolate from a pony . Therefore, whilst no differences in Piromyces existed between pony and donkey in this study, it cannot be concluded that the fibrolytic activity of the Piromyces in donkey was comparable to that of pony. To this end, studies that assess anaerobic fungal community composition and function at strain-level resolution are warranted.
From this study, the findings clearly indicate that fecal microbial community composition differs between different equine types. When comparing donkey and pony, the observed differences all related to an increase in the relative abundance or diversity of taxa with known, or potential, roles in plant material degradation. These findings are consistent with previous reports of donkeys having an increased ability to digest fiber relative to horses/ponies [8, 9]. As such, this information may provide promising avenues to enhance fiber degradation in the equine hindgut, for example by assessing the probiotic potential of Lachnoclostridium phytofermentans or anaerobic fungi to enhance fiber degradation in the equine hindgut.
The driver(s) for the observed difference in fecal microbiota between donkey and pony in this study remains to be determined. For example, the observed differences in the fecal microbiota may be the cause or an effect of the previously reported  decreased feed intake and increased MRT in donkeys compared to ponies. Further investigations are therefore now needed to expand on the novel findings reported here, particularly as measurements of intake, digestibility, MRT and fermentation metabolites were not performed in this study alongside the microbiota analysis.
Additional insight may also be gained by studying mules and hinnies separately in future studies. Differences may exist between mules and hinnies due to microbiota acquired from the mother during birth and early life, particularly due to coprophagy . Coprophagic behavior has been speculated to imprint feed selective behavior , which also has important implications for the host due to diet having a major impact on the equine hindgut microbiome .
When fed the same forage diet, fecal prokaryotic and anaerobic fungal community composition significantly differed between equine types, but not prokaryotic alpha diversity, fecal microbial concentrations or fecal DM content. Donkeys, compared to both pony and pony × donkey, had higher relative abundances of two bacterial genera, Lachnoclostridium 10 and ‘probable genus 10’, that have known, or a potential, role in plant fiber degradation. Furthermore, donkeys also had an increased number of anaerobic fungal OTUs and a higher relative abundance of the uncultivated anaerobic fungal genus SK3 compared with pony and pony × donkey, but a lower relative abundance of Piromyces compared with pony × donkey. These findings are consistent with the previously reported increased fiber degradation in donkeys compared to ponies, and suggests that the hindgut microbiota play a role. This offers novel opportunities to generate more energy from dietary fiber through management of the microbiota, decreasing the need for energy dense feeds which are a risk factor for gut-mediated disease.
Animals, diet and management
This study was conducted to assess differences in fecal microbiota between two different equine species (pony, donkey) and their derived hybrids (i.e. pony × donkey). The adult animals (n = 8 for each type) used in the experiment were selected for this study based on the following criteria: (i) clinically healthy with no reported medication/illness in the 6 months prior to the study, (ii) no known history of gut-related disorders and (iii) history (health, diet and management) information available for a minimum of 12 months prior to the start of the experiment. Furthermore, the animals involved in this study were all clinically healthy in terms of their parasite profiles, and none of the animals had a history of endoparasite related disease. Further details of the individual animals in terms of their age, sex and weight are given in Additional file 1: Table S1.
During the study, all the animals were group housed in open barns on the same farm, with free access to an outside yard area, water and a mineral lick. For a 4 week period, all animals were provided straw ad libitum and haylage. Haylage was supplied to the group housed animals at amounts known to result in maintenance of consistent bodyweight. All animals were fed exactly the same batches of forage for the last 2 weeks of the 4 week period.
Fecal sample collection and determination of dry matter content
On the last day of the 4 week period, for each animal the first feces produced after 9 A.M. was collected from the ground immediately following defecation. Parts of the feces that were visibly free of dirt, bedding etc. were placed into a clean bucket and then a pre-weighed tube was filled (approx. 20–30 g wet weight). The filled tubes were then weighed before being placed on wet ice. Samples were kept on wet ice for a maximum of 1 hour before being stored at − 20 °C. Fecal samples were then freeze-dried to a constant weight. For each sample, the percentage fecal dry matter content was then calculated using the original wet weight and the final freeze-dried weight.
The freeze dried fecal material was broken up by hand, and any large fibrous particles cut into smaller pieces using a sterile scalpel. The material was then placed into a mortar and manually ground with a pestle. Total DNA was extracted from 25 mg of the freeze-dried and ground fecal samples using the MoBio PowerSoil DNA isolation kit (QIAGEN Benelux BV, Venlo, Netherlands). The manufacturer’s protocol was followed except that after the addition of buffer C1, the samples in the PowerBead tubes were processed in a bead beater (Precellys 24, Bertin technologies, Montigny-le-Bretonneux, France) for 3 × 1 min at 5.5 m/s. DNA extracts were then further purified using the Zymo Research OneStep PCR inhibitor removal kit (BaseClear Lab Products, Leiden, Netherlands) following manufacturer’s instructions. The purity of the resulting DNA extract was assessed using a NanoDrop ND-1000 spectrophotometer (NanoDrop® Technologies, Wilmington, DE, USA), and the quantity determined using a Qubit dsDNA BR assay (Thermo Scientific, Breda, Netherlands).
Determination of microbial concentrations
For absolute quantification of bacteria and archaea, SYBR green qPCR assays were performed with sample DNA extracts using a CFX384 Touch™ Real-Time PCR Detection System (Bio-Rad Laboratories BV, Veenendaal, Netherlands) as previously described . All qPCR analyses were carried out in triplicate with a reaction volume of 10 μL and 2 ng of sample DNA extracts. Equine specific standard curves (108 to 102 amplicon copies/μL) for the assays were prepared using purified PCR amplicons generated from an equine fecal DNA extract using the primers and cycling conditions previously described for the preparation of qPCR standards . The bacterial PCR amplicons were generated using the primers 27F (5′- AGAGTTTGATCCTGGCTCAG-3′  and PROK1492R (5′- GGWTACCTTGTTACGACTT-3′ ). The archaeal PCR amplicons were generated using the primers 25F (5′- CYGGTTGATCCTGCCRG-3′  and PROK1492R (5′- GGWTACCTTGTTACGACTT-3′ ).
For absolute quantification of anaerobic fungi, a Taqman probe based method was used as previously described  with the exception that a CFX384 Touch™ Real-Time PCR Detection System (Bio-Rad Laboratories BV) was used. All qPCR analyses were carried out in triplicate with a reaction volume of 10 μL, and 20 ng of sample DNA extracts were used. Standard curves (108 to 101 amplicon copies/μL) for the assays were prepared using purified PCR amplicons generated from Neocallimastix frontalis strain RE1 DNA (kindly provided by Dr. Tony Callaghan, Bavarian State Research Center for Agriculture, Freising, Germany). The PCR amplicon was generated using the primers Neo18SF (5′-AATCCTTCGGATTGGCT-3′  and AF LSU reverse (5′-CTTGTTAAMYRAAAAGTGCATT-3′ ).
Prokaryotic community composition analysis
For 16S rRNA gene based prokaryotic community composition profiling, barcoded amplicons from the V4 region of 16S rRNA genes were generated from the DNA extracts. Primers for the V4 region and individual sample-specific barcoding strategy were as previously described . PCR was performed in a total volume of 50 μL containing 1× HF buffer (Finnzymes, Vantaa, Finland), 1 μL dNTP Mix (10 mM; Promega Benelux, Leiden, Netherlands), 1 U of Phusion® Hot Start II High-Fidelity DNA polymerase (Finnzymes), 500 nM of each sample-specific barcoded primer and 2 ng of sample DNA. The cycling conditions consisted of an initial denaturation at 98 °C for 30 s followed by 25 cycles of 98 °C for 10 s, 56 °C for 10 s and 72 °C for 10 s, and then a final extension at 72 °C for 7 min. Triplicate sample-specific barcoded PCR reactions were prepared for each sample, along with a non-template control (NTC) reaction. The presence of the sample-specific barcoded PCR products was assessed by agarose gel electrophoresis on a 2% (w/v) agarose gel containing 1× SYBR® Safe (Thermo Scientific), and the NTC reactions were confirmed to be negative. Pooled triplicate sample-specific barcoded reactions were then purified using HighPrep™ (MagBio Europe Ltd., Kent, United Kingdom) and quantified using a Qubit dsDNA BR Assay Kit (Thermo Scientific). Purified sample-specific barcoded PCR products were mixed in equimolar amounts into pools together with defined synthetic mock communities which allow assessment of potential technical biases . Pools then underwent adaptor ligation followed by sequencing on the Illumina HiSeq4000 using 150 paired end (PE) sequencing (GATC-Biotech, Konstanz, Germany, now part of Eurofins Genomics Germany GmbH).
The 16S rRNA gene sequencing data was analyzed using NG-Tax 2.0 , which executes four major tasks: demultiplexing and amplicon read cleaning, OTU-picking, denoising, and taxonomic assignment. NG-Tax 2.0 defines OTUs using an open reference approach, and OTUs are defined as unique sequences that are above a user-defined minimum abundance threshold. NG-Tax 2.0 was run with the following default settings: 70 nt read length (i.e. 140 nt in total due to being paired-end data), ratio OTU abundance 2.0, classify ratio 0.8, minimum percentage threshold 0.1%, identity level 100% and error correction of one mismatch (98.5%). Paired-end libraries were filtered to contain only read pairs with perfectly matching barcodes, and those barcodes were used to demultiplex reads by sample. The chimera detection process in NG-Tax uses the following condition: if the forward and reverse read of the OTU are identical to two different OTUs in the same sample and the abundance of the matched OTUs are at least twice of the abundance, then the OTU is marked as chimeric. Taxonomy was assigned to OTUs in NG-Tax 2.0 as previously described  using the 128 version of the SILVA 16S rRNA gene reference database .
Anaerobic fungal community composition analysis
For anaerobic fungal community composition profiling, barcoded amplicons comprising the partial 18S rRNA gene (~ 130 bp), full ITS1 region and partial 5.8S rRNA gene (~ 31 bp) were generated using a 2-step PCR strategy with a SensoQuest Labcycler as previously described . The first PCR step was performed using previously published ARISA primers  with the addition of UniTag adapters (underlined): Neo 18S For 5′-GAGCCGTAGCCAGTCTGCAATCCTTCGGATTGGCT-3′ and Neo 5.8S Rev. 5′-GCCGTGACCGTGACATCGCGAGAACCAAGAGATCCA-3′. PCR was performed in a total volume of 25 μL containing 1× HF buffer, 1 μL dNTP Mix (10 mM), 1 U of Phusion® Hot Start II High-Fidelity DNA polymerase, 500 nM of each primer and 2 ng of sample DNA. The cycling conditions consisted of an initial denaturation at 98 °C for 3 min followed by 40 cycles of 98 °C for 10 s, 58 °C for 30 s and 72 °C for 30 s, and then a final extension at 72 °C for 6 min. Triplicate PCR reactions were prepared for each sample, along with a non-template control (NTC) reaction. The presence of the PCR products was assessed by agarose gel electrophoresis on a 2% (w/v) agarose gel containing 1× SYBR® Safe. Pooled triplicate reactions, as well as the negative NTC reaction (due to the high number of PCR cycles), were then purified using HighPrep™.
The second PCR step was then employed to add an 8 nucleotide sample specific barcode to the 5′- and 3′- end of the PCR products as previously described . Each PCR reaction, with a final volume of 100 μL, contained 5 μL of the purified first step PCR product, 5 μL each of barcoded forward and reverse primers (10 μM), 2 μL dNTP Mix (10 mM), 2 U of Phusion® Hot Start II High-Fidelity DNA polymerase and 1× HF buffer. Amplification consisted of an initial denaturation at 98 °C for 30 s followed by 5 cycles of 98 °C for 10 s, 52 °C for 20 s and 72 °C for 20 s, and then a final extension at 72 °C for 10 min. Barcoded PCR products were then purified using the HighPrep™ and quantified using a Qubit dsDNA BR Assay Kit. Purified barcoded PCR products were then pooled in equimolar amounts along with defined synthetic mock communities . Pools were then sequenced on the Illumina HiSeq 2500 using the Rapid Run 300 bp PE sequencing mode (GATC-Biotech).
The anaerobic fungal sequence data was then analyzed using NG-Tax 2.0 as previously described . NG-Tax 2.0 was run using the default parameters (as specified earlier) except for the following: 150 nt read length (i.e. 300 nt in total due to being paired-end data), minimum percentage threshold 0.6% and error correction of one mismatch (99.3%). As the barcoded amplicon primers used were not within the AF-ITS1 database used for OTU annotation (which is a requirement for annotation by NG-Tax), an empty database file (emptydb.fasta.gz which is available at http://download.systemsbiology.nl/ngtax/databases/) was used and the OTUs then subsequently annotated manually.
Fasta files of the OTUs from the NG-Tax generated biom file were extracted using the script otuseq_export.py (which is available at https://gitlab.com/wurssb/gen_fake_mocks/tree/master/paper_data). The OTUs were annotated using BLASTN searches against the AF-ITS1 database  (version 3.3, available from www.anaerobicfungi.org) using default settings with “-num_alignments 10” (BLAST version 2.4.0). For OTUs that could not be annotated by the AF-ITS1 database, BLASTN searches were performed against the NCBI database. Cut-off levels for OTU annotations were determined based on the mean percentage similarities of full-length sequences in the AF-ITS1 database within clade and within genus. These cut-off levels were > 98% for clade and > 95% for genus. The NG-Tax generated biom file was converted to a tab delimited table to enable OTU annotations to be added. The OTUs that were clearly associated with the NTC sample were also manually removed from the tab delimited table at this stage, along with any OTUs that were not anaerobic fungal in origin. The resulting tab delimited table was then converted back to a biom file.
Microbial composition summary box plots and UniFrac based Principal Coordinate Analysis (PCoA) were generated within R (version 3.4.1)  using the following libraries and packages: microbiome (https://microbiome.github.io/tutorials/), microbiomeutilities (https://github.com/microsud/microbiomeutilities), RColorBrewer , magrittr , phyloseq , picante , nlme , vegan , lattice , permute , ape , ggplot2 , and ggpubr . QIIME 1  was used to generate taxa genus level biom tables using the script “summarize_taxa.py”. Redundancy analysis (RDA) was performed using Canoco 5  to assess the relationship between genus-level phylogenetic groupings of the OTUs and equine type. The QIIME 1 script “group.signficance.py” was used to test differences in relative abundance of individual genera between equine types using Kruskal Wallis with Bonferroni correction of P values. A Dunn’s-Sidak post-hoc test with Bonferroni correction was then performed on genera that were significantly different between equine types (P < 0.05), in order to determine which equine types significantly differed from each other (MATLAB). The QIIME 1 script “alpha.diversity.py” was used to determine per sample the number of observed OTUs and the Phylogenetic Diversity (i.e. PD_whole_tree) value. Number of observed OTUs, ‘Phylogenetic Diversity’ and fecal dry matter content were analyzed using a one-way ANOVA with equine type as a single independent factor (Genstat 18th Edition, VSN International Ltd.). All qPCR data was analyzed in the same manner after a Log10 transformation. For all statistical tests the significance threshold was alpha = 0.05. P values for multivariate data were all Bonferroni corrected (as indicated above).
Availability of data and materials
The datasets and material supporting the conclusions of this article are provided as follows. Additional information is provided in Additional files 1 and 2. The barcoded amplicon sequence data is deposited in the European Nucleotide Archive under the study accession number PRJEB32772. All the sample barcodes, R codes and data used in the analysis are available at https://github.com/mibwurrepo/EdwardsJ_2019_Equine_Type_Comparison, unless indicated otherwise.
Anaerobic fungal ITS1
Internal Transcribed Spacer 1
Mean retention time
Operational Taxonomic Unit
Principal co-ordinate analysis
Polymerase chain reaction
Argenzio RA, Southworth M, Stevens CE. Sites of organic acid production and absorption in the equine gastrointestinal tract. Am J Phys. 1974;22:1043–50.
Glinsky MJ, Smith RM, Spires HR, Davis CL. Measurement of volatile fatty acid production rates in the cecum of the pony. J Anim Sci. 1976;42:1465–70.
Bergman EN. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev. 1990;70:567–90.
Durham AE. Intestinal disease. In: Geor RJ, Harris PA, Coenen M, editors. Equine applied and clinical nutrition: health, welfare and performance. Philadelphia: Saunders Elsevier Ltd; 2013. p. 568–81.
Julliand V, Grimm P. The impact of diet on the hindgut microbiome. J Equine Vet Sci. 2017;52:23–8.
FAOSTAT. Food and Agriculture Organization for the United Nations, Rome. 2017. http://faostat.fao.org. Accessed 11 April 2019.
Burden F, Thiemann A. Donkeys are different. J Equine Vet Sci. 2015;35:376–82.
Pearson RA, Merritt JB. Intake, digestion and gastrointestinal transit time in resting donkeys and ponies and exercised donkeys given ad libitum hay and straw diets. Equine Vet J. 1991;23:339–43.
Cuddeford D, Pearson RA, Archibald RF, Muirhead R. Digestibility and gastro-intestinal transit time of diets containing different proportions of alfalfa and oat straw given to thoroughbreds, Shetland ponies, Highland ponies and donkeys. Anim Sci. 1995;61:407–17.
Smith DG, Pearson RA. A review of the factors affecting the survival of donkeys in semi-arid regions of sub-Saharan Africa. Trop Anim Health Prod. 2005;37:S1–19.
Li Z, Wright ADG, Si H, Wang X, Qian W, Zhang Z, et al. Changes in the rumen microbiome and metabolites reveal the effect of host genetics on hybrid crosses. Environ Microbiol Rep. 2016;8:1016–23.
Julliand V, Grimm P. Horse species symposium: the microbiome of the horse hindgut: history and current knowledge. J Anim Sci. 2016;94:2262–74.
Moore BE, Dehority BA. Effects of diet and hindgut defaunation on diet digestibility and microbial concentrations in the cecum and colon of the horse. J Anim Sci. 1993;71:3350–8.
Julliand V, De Vaux A, Millet L, Fonty G. Identification of Ruminococcus flavefaciens as the predominant cellulolytic bacterial species of the equine cecum. Appl Environ Microbiol. 1999;65:3738–41.
Julliand V, Riondet C, De Vaux A, Alcaraz G, Fonty G. Comparison of metabolic activities between Piromyces citronii, an equine fungal species, and Piromyces communis, a ruminal species. Anim Feed Sci Technol. 1998;70:161–8.
Solomon KV, Haitjema CH, Henske JK, Gilmore SP, Borges-Rivera D, Lipzen A, et al. Early-branching gut fungi possess large, comprehensive array of biomass-degrading enzymes. Science. 2016;351:1192–5.
Gruninger RJ, Puniya AK, Callaghan TM, Edwards JE, Youssef N, Dagar SS, et al. Anaerobic fungi (phylum Neocallimastigomycota): advances in understanding their taxonomy, life cycle, ecology, role and biotechnological potential. FEMS Microbiol Ecol. 2014;90:1–17.
Gordon GLR, Phillips MW. The role of anaerobic gut fungi in ruminants. Nutr Res Rev. 2005;11:133–68.
Blackmore TM, Dugdale A, Argo CM, Curtis G, Pinloche E, Harris PA, et al. Strong stability and host specific bacterial Community in Faeces of ponies. PLoS One. 2013;8:e75079.
Morrison PK, Newbold CJ, Jones E, Worgan HJ, Grove-White DH, Dugdale AH, et al. The equine gastrointestinal microbiome: impacts of age and obesity. Front Microbiol. 2018;9:3017.
Grimm P, Philippeau C, Julliand V. Faecal parameters as biomarkers of the equine hindgut microbial ecosystem under dietary change. Animal. 2017;11:1136–45.
Daly K, Stewart CS, Flint HJ, Shirazi-Beechey SP. Bacterial diversity within the equine large intestine as revealed by molecular analysis of cloned 16S rRNA genes. FEMS Microbiol Ecol. 2001;38:141–51.
Shepherd ML, Swecker WS, Jensen RV, Ponder MA. Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons. FEMS Microbiol Lett. 2012;326:62–8.
Costa MC, Stämpfli HR, Allen-Vercoe E, Weese JS. Development of the faecal microbiota in foals. Equine Vet J. 2016;48:681–8.
Costa MC, Stämpfli HR, Arroyo LG, Allen-Vercoe E, Gomes RG, Weese J. Changes in the equine fecal microbiota associated with the use of systemic antimicrobial drugs. BMC Vet Res. 2015;11:19.
Rodriguez C, Taminiau B, Brévers B, Avesani V, Van Broeck J, Leroux A, et al. Faecal microbiota characterisation of horses using 16 rdna barcoded pyrosequencing, and carriage rate of Clostridium difficile at hospital admission. BMC Microbiol. 2015;15:181.
Fernandes KA, Kittelmann S, Rogers CW, Gee EK, Bolwell CF, Bermingham EN, et al. Faecal microbiota of forage-fed horses in New Zealand and the population dynamics of microbial communities following dietary change. PLoS One. 2014;9:e112846.
Steelman SM, Chowdhary BP, Dowd S, Suchodolski J, Janečka JE. Pyrosequencing of 16S rRNA genes in fecal samples reveals high diversity of hindgut microflora in horses and potential links to chronic laminitis. BMC Vet Res. 2012;8:231.
Dougal K, Harris PA, Girdwood SE, Creevey CJ, Curtis GC, Barfoot CF, et al. Changes in the total fecal bacterial population in individual horses maintained on a restricted diet over 6 weeks. Front Microbiol. 2017;8:1502.
Dougal K, Harris PA, Edwards A, Pachebat JA, Blackmore TM, Worgan HJ, et al. A comparison of the microbiome and the metabolome of different regions of the equine hindgut. FEMS Microbiol Ecol. 2012;82:642–52.
Dougal K, de la Fuente G, Harris PA, Girdwood SE, Pinloche E, Newbold CJ. Identification of a Core bacterial community within the large intestine of the horse. PLoS One. 2013;8:e77660.
Liu X, Fan H, Ding X, Hong Z, Nei Y, Liu Z, et al. Analysis of the gut microbiota by high-throughput sequencing of the v5-v6 regions of the 16s rRNA gene in donkey. Curr Microbiol. 2014;68:657–62.
Gao H, Chi X, Qin W, Wang L, Song P, Cai Z, et al. Comparison of the gut microbiota composition between the wild and captive Tibetan wild ass (Equus kiang). J Appl Microbiol. 2019;126:1869–78.
Jónsson H, Schubert M, Seguin-Orlando A, Ginolhac A, Petersen L, Fumagalli M, et al. Speciation with gene flow in equids despite extensive chromosomal plasticity. Proc Natl Acad Sci. 2014;111:18655–60.
Costa MC, Arroyo LG, Allen-Vercoe E, Stämpfli HR, Kim PT, Sturgeon A, et al. Comparison of the fecal microbiota of healthy horses and horses with colitis by high throughput sequencing of the V3-V5 region of the 16S rRNA gene. PLoS One. 2012;7:e41484.
Dougal K, de la Fuente G, Harris PA, Girdwood SE, Pinloche E, Geor RJ, et al. Characterisation of the faecal bacterial community in adult and elderly horses fed a high fibre, high oil or high starch diet using 454 pyrosequencing. PLoS One. 2014;9:e87424.
Costa MC, Silva G, Ramos RV, Staempfli HR, Arroyo LG, Kim P, et al. Characterization and comparison of the bacterial microbiota in different gastrointestinal tract compartments in horses. Vet J. 2015;205:74–80.
Warnick TA, Methé BA, Leschine SB. Clostridium phytofermentans sp. nov., a cellulolytic mesophile from forest soil. Int J Syst Evol Microbiol. 2002;52:1155–60.
Yutin N, Galperin MY. A genomic update on clostridial phylogeny: gram-negative spore formers and other misplaced clostridia. Environ Microbiol. 2013;15:2631–41.
Petit E, Coppi MV, Hayes JC, Tolonen AC, Warnick T, Latouf WG, et al. Genome and transcriptome of Clostridium phytofermentans, catalyst for the direct conversion of plant feedstocks to fuels. PLoS One. 2015;10:e0118285.
Seshadri R, Leahy SC, Attwood GT, Teh KH, Lambie SC, Cookson AL, et al. Cultivation and sequencing of rumen microbiome members from the Hungate1000 collection. Nat Biotechnol. 2018;36:359–36.
Nyonyo T, Shinkai T, Mitsumori M. Improved culturability of cellulolytic rumen bacteria and phylogenetic diversity of culturable cellulolytic and xylanolytic bacteria newly isolated from the bovine rumen. FEMS Microbiol Ecol. 2014;88:528–37.
Nyonyo T, Shinkai T, Tajima A, Mitsumori M. Effect of media composition, including gelling agents, on isolation of previously uncultured rumen bacteria. Lett Appl Microbiol. 2013;56:63–70.
Kenters N, Henderson G, Jeyanathan J, Kittelmann S, Janssen PH. Isolation of previously uncultured rumen bacteria by dilution to extinction using a new liquid culture medium. J Microbiol Methods. 2011;84:52–60.
Hanafy RA, Elshahed MS, Youssef NH. Feramyces austinii, gen. Nov., sp. nov., an anaerobic gut fungus from rumen and fecal samples of wild barbary sheep and fallow deer. Mycologia. 2018;110:513–25.
Edwards JE, Hermes GDA, Kittelmann S, Nijsse B, Smidt H. Assessment of the accuracy of Neocallimastigomycota ITS1 based barcoded amplicon sequencing for community composition analysis. Front Microbiol. 2019;10:2370.
Koetschan C, Kittelmann S, Lu J, Al-Halbouni D, Jarvis GN, Müller T, et al. Internal transcribed spacer 1 secondary structure analysis reveals a common core throughout the anaerobic fungi (Neocallimastigomycota). PLoS One. 2014;9:e91928.
Kittelmann S, Naylor GE, Koolaard JP, Janssen PH. A proposed taxonomy of anaerobic fungi (class Neocallimastigomycetes) suitable for large-scale sequence-based community structure analysis. PLoS One. 2012;7:e36866.
Fliegerová K, Mrázek J, Hoffmann K, Zábranská J, Voigt K. Diversity of anaerobic fungi within cow manure determined by ITS1 analysis. Folia Microbiol. 2010;55:319–25.
Mura E, Edwards J, Kittelmann S, Kaerger K, Voigt K, Mrázek J, et al. Anaerobic fungal communities differ along the horse digestive tract. Fungal Biol. 2019;123:240–6.
Liggenstoffer AS, Youssef NH, Couger MB, Elshahed MS. Phylogenetic diversity and community structure of anaerobic gut fungi (phylum Neocallimastigomycota) in ruminant and non-ruminant herbivores. ISME J. 2010;4:1225–35.
Gold JJ, Brent Heath I, Bauchop T. Ultrastructural description of a new chytrid genus of caecum anaerobe, Caecomyces equi gen. nov., sp. nov., assigned to the Neocallimastigaceae. BioSystems. 1988;21:403–15.
Orpin CG. Isolation of cellulolytic Phycomycete Fungi from the Caecum of the horse. Microbiology. 1981;123:287–96.
Li J, Heath IB, Bauchop T. Piromyces mae and Piromyces dumbonica , two new species of uniflagellate anaerobic chytridiomycete fungi from the hindgut of the horse and elephant. Can J Bot. 1990;68:1021–33.
Gaillard-Martinie B, Breton A, Dusser M, Julliand V. Piromyces citronii sp. nov., a strictly anaerobic fungus from the equine caecum: a morphological, metabolic, and ultrastructural study. FEMS Microbiol Lett. 1995;130:321–6.
Callaghan TM, Podmirseg SM, Hohlweck D, Edwards JE, Puniya AK, Dagar SS, et al. Buwchfawromyces eastonii gen. nov., sp. nov.: a new anaerobic fungus (Neocallimastigomycota) isolated from buffalo faeces. MycoKeys. 2015;9:11–28.
Francis-Smith K, Wood-Gush DGM. Coprophagia as seen in thoroughbred foals. Equine Vet J. 1977;9:155–7.
Marinier SL, Alexander AJ. Coprophagy as an avenue for foals of the domestic horse to learn food preferences from their dams. J Theor Biol. 1995;173:121–4.
van Lingen HJ, Edwards JE, Vaidya JD, van Gastelen S, Saccenti E, van den Bogert B, et al. Diurnal dynamics of gaseous and dissolved metabolites and microbiota composition in the bovine rumen. Front Microbiol. 2017;8:425.
Lane D. 16S/23S rRNA sequencing. In: Stackebrandt E, Goodfellow M, editors. Nucleic acid techniques in bacterial systematics. Wiley: New York; 1991. p. 115–48.
Suzuki MT, Taylor LT, DeLong EF. Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5′-nuclease assays. Appl Environ Microbiol. 2000;66:4605–14.
Dojka MA, Hugenholtz P, Haack SK, Pace NR. Microbial diversity in a hydrocarbon- and chlorinated-solvent contaminated aquifer undergoing intrinsic bioremediation. Appl Environ Microbiol. 1998;64:3869–77.
Edwards JE, Kingston-Smith AH, Jimenez HR, Huws SA, Skøt KP, Griffith GW, et al. Dynamics of initial colonization of nonconserved perennial ryegrass by anaerobic fungi in the bovine rumen. FEMS Microbiol Ecol. 2008;66:537–45.
Dollhofer V, Callaghan TM, Griffith GW, Lebuhn M, Bauer J. Presence and transcriptional activity of anaerobic fungi in agricultural biogas plants. Bioresour Technol. 2017;235:131–9.
Ramiro-Garcia J, Hermes GDA, Giatsis C, Sipkema D, Zoetendal EG, Schaap PJ, et al. NG-Tax, a highly accurate and validated pipeline for analysis of 16S rRNA amplicons from complex biomes. F1000Research. 2018;5:1791.
Poncheewin W, Hermes GDA, van Dam JCJ, Koehorst JJ, Smidt HSPJ. NG-tax 2.0: a semantic framework for high-throughput amplicon analysis. Front Genet. 2019. https://doi.org/10.3389/fgene.2019.01366.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. URL http://www.R-project.org/.
Neuwirth E. RColorBrewer: ColorBrewer palettes. R package version 1.1–2. 2014. https://cran.r-project.org/package=RColorBrewer.
Milton Bache S, Wickham H. magrittr: a forward-pipe operator for R. R package version 1.5. 2014. https://cran.r-project.org/package=magrittr.
McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.
Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. _nlme: linear and nonlinear mixed effects models_. R package version 3.1–131. 2017. https://cran.r-project.org/package=nlme.
Oksanen J, Guillaume Blanchet F, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. R package version 2.5–3. 2018. https://cran.r-project.org/package=vegan.
Deepayan S. Lattice: multivariate data visualization with R. New York: Springer; 2008.
Simpson GL. permute: functions for generating restricted permutations of data. 2016. https://cran.r-project.org/package=permute.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer-Verlag; 2016.
Kassambara A. ggpubr: “ggplot2” based publication ready plots. R package version 0.2. https://cran.r-project.org/package=ggpubr.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Šmilauer P, Leps J. Multivariate analysis of ecological data using Canoco 5. 2nd ed. New York: Cambridge University Press; 2014.
We would like to acknowledge all of the members of staff from The Donkey Sanctuary that assisted with this study.
This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement number 706899.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Table S1. Details of the animals used in this study.
: Figure S1. Boxplot showing the six main bacterial and archaeal phyla detected in the different equine types. The minor phyla (< 1%) are grouped as ‘Other’. Boxes show the 25th and 75th percentiles with the median represented by a horizontal line. Whiskers show the data range with the exception of any outliers, which are indicated as data points.
About this article
Cite this article
Edwards, J.E., Schennink, A., Burden, F. et al. Domesticated equine species and their derived hybrids differ in their fecal microbiota. anim microbiome 2, 8 (2020). https://doi.org/10.1186/s42523-020-00027-7
- Barcoded amplicon sequencing
- Anaerobic fungi