Homeostasis of the Intestinal Mucosa in Healthy Horses—Correlation between the Fecal Microbiome, Secretory Immunoglobulin A and Fecal Egg Count

Simple Summary The defense effect of the gastrointestinal mucosa in horses, as well as in other mammals, is an extremely complex process and is dependent on many parameters. Mutual communication between the host and the microbiome is essential for intestinal mucosa homeostasis. In horses, in the described system, constant exposure to nematode eggs infection is also important. The presented studies showed the existence of significant relationships between the host response and composition and metabolic activity of microbiome as well as the presence of small parasites in the digestive tract of horses. Greater microbiome diversity correlates with greater production of the neutralizing immunoglobulin A, and a key commensal bacteria phylum in horses, Firmicutes, is negatively correlated with fecal egg count. Abstract The defensive function of the intestinal mucosa depends both on the ability to secrete immunoglobulin A and communication with the mucus microbiome. In horses, the functioning of this system is also influenced by the presence of nematode eggs. Feces collected from healthy horses were examined to determine the fecal egg count, immunoglobulin A level (ELISA), microbiome composition (Next-Generation Sequencing, NGS, V3–V4 and V7–V9 hypervariable regions of the 16S rRNA gene analysis and short-chain fatty acid (SCFA) production ((high-performance liquid chromatography, HPLC). In the taxonomic analysis within the phylum, the following order of dominance was found: Firmicutes, Bacteroidota, Verrucomicrobiota and Fibrobacterota. The coefficient of phylogenetic diversity of the microbiome positively correlated with both secretory immunoglobulin A (SIgA) [μg/g of feces] (p = 0.0354, r = 0.61) and SIgA [μg/mg of fecal protein] (p = 0.0382, r = 0.6) and with the number of Cyathostomum eggs (p = 0.0023, r = 0.79). Important components of the key microbiome in horses, such as phylum Proteobacteria and species Ruminococcus flavefaciens, were positively correlated with the fecal SIgA (p < 0.05). All the obtained results indicate the existence of significant relationships between the host response (SIgA production) and composition and SCFA production in the microbiome as well as the presence of small strongyles in the digestive tract of horses.


Introduction
The defense effect of the gastrointestinal mucosa in horses, as well as in other mammals, is an extremely complex process and is dependent on many parameters. The physical barrier consists of a mucus layer and closely adjacent mucosal epithelial cells (MECs). The mucous barrier provides the first immune defense barrier for pathogens and anchoring for protective factors [1]. The outer mucous layer is colonized by the intestinal commensal microbiota. The inner layer is responsible for the protection of MECs against bacteria, including commensal ones, via the presence of secretory immunoglobulin A (SIgA) and antibacterial peptides (AMPs) such as defensins or lactoferricins [2]. Their coexistence is known as the "kill zone" [3]. A healthy gastrointestinal microbiota demonstrates a protective effect within the mucosa in three ways: (a) competition for the site of colonization and availability of nutrients; (b) within the kill zone, commensal bacteria stimulate the production of mucins, SIgA and AMPs; (c) by enhancing the host's immune response to pathogens by stimulating the cells to produce IL-22 (T cells, NK cells) and IL-1β (intestinal monocytes) [4]. The production and quantification of SIgA in the feces may be a marker of gastrointestinal-associated lymphoid tissue (GALT) function. As commensal microflora in the equine hindgut, Firmicutes, Bacteroidetes and Verrucomicrobia are amongst the predominating phyla [5][6][7][8]. An important function performed by commensal bacteria is the production of short-chain fatty acids (SCFAs). SCFAs are a source of energy, and can also have a protective and repairing effect on intestinal cells. They also promote proliferation and the de novo induction of Treg cells and suppress production of proinflammatory cytokines via macrophages [9]. When discussing the dependencies on the functioning of GALT in horses, one cannot forget about the high exposure of this area to parasites-especially small strongyles, which, despite parasite control programs that include anthelmintic, inhabit the gastrointestinal tract [10]. The interrelationships between the presence of parasites and the functioning of the mucosal "kill zone" in horses have not been fully established.
The main objective of the study was to determine the relationship between the "kill zone" components, such as SIgA and the gut microbiome, as well as the degree of invasion of small strongyles in horses. To address these questions using feces collected from healthy horses, the number of Cyathostomum eggs, SIgA level, microbiota composition and the ability to produce SCFAs were assessed. An additional goal was to try to determine the stool SIgA and SCFAs as potential markers of the functioning of the immune system within the gastrointestinal tract.

Experimental Design
The study was conducted in 60 clinically healthy adult horses of different sex, age (>2.5 y.o.) and breed, housed in three different breeding centers (detailed in Table 1). In stable No1, all test horses were fed equally-they had ad libitum access to grass pasture, high-quality hay, water and mineral licks. In stables No2 and No3, the horses had ad libitum access to a grass pasture, high-quality hay, water and mineral licks; additionally, they were fed concentrated feed according to their needs. In the groups from stables No2 and No3, it was impossible to standardize the way of feeding the horses for the purposes of the study. All horses were routinely dewormed with oral pastes containing ivermectin (18.7 mg/g) and praziquantel (140.3 mg/g) in doses of 1.07 g paste/100 kg body mass. The deworming program was defined for each stable and included deworming twice a year with a two-component preparation process. According to the schedule at the time of study, in stables No1 and No3, the horses were dewormed 6 months before sampling and in stable No2 1.5 months before sampling. According to the information from the owner, none of the horses showed symptoms of gastrointestinal diseases such as colic or diarrhea and had not been treated in the last two months. Feces samples were collected from each horse during spontaneous defecation (immediately after defecation in a volume of 100 mL for further determination and preservation) [11]. The samples were divided and protected by freezing at −80 • C. Immediately after sampling, 4 g of fresh feces was weighted and collected and stored in a refrigerator (2-8 • C) for fecal egg count (FEC) analysis. In stable No1, samples were additionally collected from the center of the stool clump in a sterile manner to avoid contamination with environmental factors: 1 g of feces was weighed and suspended in 5 mL of 70% ethanol and a sample of about 1 g of feces was placed in sterile Eppendorf tubes [12]. Both types of samples were frozen (−20 • C) immediately after sampling and then in −80 • C after transport. Additionally, in connection with the routine blood tests in stable No1, blood samples were collected from the external jugular vein into a tube with a clot activator and then centrifuged (10 min, 1500 g) to obtain serum and frozen (−80 • C). The following analyses were carried out in fecal samples from horses from the No1 stable, due to the uniform method of maintenance and feeding: FEC, fecal protein, fecal SIgA, fecal microbiome and SCFA production. In the remaining samples, due to the difficulty of systematizing conditions, the fecal microbiome was not determined. Modified McMaster method: Fecal egg count was performed to establish the number of small strongyle (Cyathostomum spp.) eggs in equine feces. FEC was performed with a modified McMaster method, as has been described before [13]. This method was chosen because it is one of the most practical ways to estimate the level of infection. Briefly, the fecal examination was performed, followed by centrifugation-an enhanced method using a sugar-salt flotation solution with a specific gravity of 1.3 g/mL [13]. The detection limit of the method was 20 eggs per gram (EPG) [13]. For the value of 0 eggs in the field of view, a value of 20 was adopted, to minimize statistical mistakes.

Fecal Protein Analysis
After thawing at room temperature, the fecal samples were prepared as described in the literature according to Marr et al., using proprietary modifications [14]. From wet feces samples, 1 g (±0.05 g) of each was weighed and suspended in 5 mL of phosphatebuffered saline (PBS) solution (pH~7.4). In order to mix them thoroughly and release the proteins, the samples were subjected to the following procedure: shaking (3 min) and resting (15 min), twice. Samples were then centrifuged for 20 min at 1600× g. The resulting centrifugation supernatant was collected, filtered through a syringe filter (syringe filter, 0.22 µm, ∅33 mm, sterile, TPP, 99722) and then placed in an Eppendorf tube containing a protease inhibitor cocktail (10 µL/1000 µL supernatant; Calbiochem Protease Inhibitor Cocktail Set I 539131). Samples were centrifuged again for 15 min at 3260× g. A total of 600 µL of supernatant was withdrawn for protein content assessment, and the excess was frozen at −80 • C for further SIgA analysis. Analysis of fecal protein of the prepared samples was performed using the commercial Pierce BCA Protein Assay Kit (ThermoFisher Scientific, Rockford, IL, USA) according to the manufacturer's instructions. After the initial determinations were made, the optimal dilution of the samples was found at the level of 1:25.

Serum IgA and Fecal SIgA Analysis
The assay was performed in the supernatant obtained as described in point 2.2 and serum samples using a commercially available ELISA kit (IgA Horse ELISA Kit, Abcam Cat. No. ab190530, Cambridge, UK) according to the manufacturer's instructions. The fecal supernatant samples and serum samples were thawed at room temperature and then diluted to 1:5 and 1:10,000, respectively. The intra-assay CV reported in these studies was CV = 1.69%, the inter-assay CV = 2.35%.

Fecal Short-Chain Fatty Acid (SCFA) Analysis
Analysis of derivatized stool extracts by high-performance liquid chromatography (HPLC) was performed as previously described [15]. Briefly, freshly picked stool samples were extracted with 70% ethanol. Debris was removed by centrifugation and 500 µL of supernatant was transferred to a new tube and mixed with 50 µL of internal standard After cooling, the mixture was shaken with 3.0 mL of phosphoric acid aqueous solution (0.5 mol/L) and 4.0 mL of ether two times for three minutes for extraction. The organic phase was extracted by shaking with 4 mL of diethyl ether and transferred to a new glass conical containing water to extract any remaining aqueous compounds. The obtained fatty acids hydrazide was dissolved in 100 µL of methanol and 20 µL was subjected to HPLC [15]. HPLC was performed using 1525 Binary HPLC Pump with a 2489 UV/Visible (UV/Vis) Detector and Phenomenex Gemini 5 µm C18 110 A (150 × 4.6 mm). The mobile phase was composed of acetonitrile-methanol-water (30:16:54), and the pH was adjusted to 4.5 with 0.1% trifluoroacetic acid. The column temperature was 50 • C; flow rate, 1.0 mL/min; and measurement wavelength, 400 nm [15]. The following SCFAs were determined: lactic acid, acetic acid (C2), propionic acid (C3), butyric acid (C4), isobutyric acid (iC4), valeric acid (C5) and isovaleric acid (iC5). The results are presented in the unit [mmol / 1g of feces]. As an important determinant of fibrolytic activity, the volatile fatty acid (VFA) ratio was counted, according to the formula: (C2 + C4)/C3 [16].

Microbiome Analysis
Microbial DNA was extracted using QIAamp PowerFecal Pro DNA Kit (Qiagen, Hilden, Germany) from frozen stool material from 12 horses. 16S rRNA gene amplicon libraries were prepared according to the standard protocol using QIAseq 16S/ITS Region Panels (Qiagen, Hilden, Germany) for V3-V4 and V7-V9 regions. Paired-end sequencing of the amplicons was performed on a MiSeq instrument, using MiSeq Reagent Kit v3, 600 cycles (Illumina, San Diego, CA, USA). To perform bioinformatic analysis of the microbiomes, QIIME2 2021.8 with the supplementary plugins was used [17]. A custom script that uses cutadapt was used to cut V3-V4 and V7-V9 primers and execute demultiplexing [18]. Using the summary method from the demux plugin, the evaluation of the quality of the reads was performed. The dada2 plugin for paired-end reads was used for fortrimming, denoising, dereplication and chimera filtering [19]. The construction of a phylogenetic tree was performed using the q2-phylogeny align-to-tree-mafft-fasttree plug procedure [20,21]. The procedure internally uses the mafft method to align multiple sequences and then masks the highly variable positions, and subsequently fasttree is used to construct a phylogenetic tree. For the V3-V4 region, rarefaction (subsampling without replacement) was conducted up to 22.039 sequences per sample and up to 18.799 for the V7-V9 region. To estimate the diversity, alpha (α) (Shannon's diversity index, observed features, Faith's phylogenetic diversity, Pielou's evenness) and beta (β) (Jaccard similarity index, Bray-Curtis dissimilarity, unweighted UniFrac, weighted UniFrac), as well as to perform principal coordinate analysis (PCoA), the q2-diversity plugin was used [21][22][23][24][25]. The naive Bayesian classifiers trained on fragments V3-V4 and V7-V9 of the 16S rRNA gene sequences derived from the SILVA 138 SSURef NR99, and the fit-classifier-naive-bayes method from the feature-classifier plugin was used to assign taxonomy to ASVs (amplicon sequence variants) [26][27][28]. The demultiplexed fastq data are publicly available on the SRA database under the accession number SAMN31431713 [26][27][28].

Statistical Analysis
Descriptive statistics are presented as mean and standard deviation. Analysis of the differences in parameters between stables was performed using the Kruskal-Wallis test with post hoc analysis using Holm correction. Analysis of correlation between parameters both within and between groups of parameters was performed using Spearman correlation coefficient. Spearman's rank correlation coefficient was performed to correlate all measurements with α-diversity metrics, and the Mantel test with 999 permutations in the case of β-diversity. False-discovery rate correction was applied to control type I error. Calculations were made using the R package for Windows (version 4.2), and for bioinformatic analysis of the microbiomes, QIIME2 2021.8 with the supplementary plugins was used [17]. Statistical results were considered significant when the p-value was below 0.05.

Fecal Egg Count, Serum IgA and sIgA Analysis Results
The results of the tests for fecal egg count, serum IgA, SigA and fecal protein content are presented in Table 2.  In stable No2, significantly lower concentrations of SIgA [µg/g of feces] and SIgA [µg/mg of fecal protein] were observed in comparison to stable No1 (p = 0.043, p = 0.013, respectively) and No3 (p = 0.043, p = 0.046, respectively). There was no correlation between the parameters (FEC, SIgA, fecal protein) and the age of the horses.
The following positive correlation between individual phyla and species and the production of SCFAs was found: C2 and Firmicutes (p= 0.0259, r = 0.65). Negative correlations were noticed between C2 and Bacteroidota (p = 0.022, r = −0.66); iC4 and Bacteroidota (p = 0.0348, r = −0.62); C4 and Spirochaetota (p = 0.0484, r = −0.59); and C4 and Escherichia Shigella p = 0.053, r = −0.57 (trend). The VFA ratio correlated with Firmicutes (p = 0.0091, r = 0.73) and Fibrobacterota (p = 0.0041, r = −0.78). The relative frequency of the genus Ruminococcus flavefaciens correlated negatively with C4 (p = 0.032, r = −0.63). The relative frequency of genus Clostridium s. stricto correlated with body weight (p = 0.0096, r = 0.71). There were negative correlations between FEC and the relative frequency of Firmicutes (p = 0.003, r = −0.78) and Desulfobacterota (p = 0.0443, r = −0.59). Relationships between some phyla and species and the level of SIgA in the feces and IgA in the serum were found; they are presented in Figure 3. Relative abundance of microbiota from horse feces at the phylum level. Taxonomy based on V3-V4 hypervariable region of 16S rRNA gene. Phyla whose relative abundance accounts for less than 0.01 are found in Others.

Discussion
In the presented study, the relationships between the production of SIgA, the microbiome and the degree of Cyathostomum infection in the gastrointestinal tract of healthy

Discussion
In the presented study, the relationships between the production of SIgA, the microbiome and the degree of Cyathostomum infection in the gastrointestinal tract of healthy horses were assessed. The mean content of SIgA in the feces of healthy horses was established at the level of 0.23 (±0.23) µg/mg of fecal protein. Equine researchers have focused only on the potential use of SIgA in the feces as a stress marker; however, the available results do not clearly confirm this relationship [14]. Studies in wild equidae in Africa have considered fecal SIgA as a parasitic nematode-associated antibody, which suggests that equids may develop an IgA response to nematode egg production when the host is in good condition [29]. Studies on the measurement of fecal immunoglobulins and their dependence on GI pathology have been carried out in dogs, and they have shown a decrease in fecal SIgA in the course of IBD [11]. The presented studies did not find any correlation between FEC and the amount of SIgA in the feces of healthy horses. However, we did not assess the presence of specific antibodies against larval or adult Cyathostomum. Significantly lower values of SIgA concentration in the feces from stable No.2 in relation to stables No.1 and No.3 (taking into account the time that had passed since the last deworming) seem to suggest a correlation between the absence of parasites in the gastrointestinal tract and a decrease in SIgA production, but confirmation of this requires further research.
The study of the microbiome was carried out on 12 horses housed in one breeding facility. This made it possible to completely unify the nutritional and climatic conditions; these horses had constant, 24 h access to a pasture and high-quality hay, water and mineral licks. They received no concentrated feed. The dominance of individual phyla was similar to that presented in the literature, usually in the following order: Firmicutes (>60%), Bacteroidetes, Proteobacteria and Verrucomicrobia [5,30,31]. However, in the presented work, the phylum Firmicutes was only about 35% in equilibrium with Bacteroidetes, and Fibrobacterota appeared as the third phylum. These differences may be due to the type of diet, which can significantly affect the composition of the microbiome [6].
In assessing the gut microbiome in specific species, it is important to know the core microbiota, which should be consistent throughout a stable in terms of bacterial components, including key microorganisms and their functions [8]. The core community at the operational taxonomic unit (OTU) level in feces is defined as "being present in all samples included in the study at least 0.1% relative abundance" [32]. According to the definition of core microbiota in our study, the phyla Firmicutes, Bacteroidota, Verrucomicrobiota, Fibrobacterota, Spirotracheota, Proteobacteria, Cyanobacteria and Desulfobacteria, and at the family level, Lachnospiraceae and Ruminococcaceae, constitute the core. However, due to high heterogeneity of the equine core microbiome, the "key microbiome" should be distinguished-species crucial for the proper functioning of the digestive tract [6]. In horses, the three most important phyla are Firmicutes, Proteobacteria and Verrucomicrobia due to having the greatest relevance in digestion, rather than due to their quantitative dominance [6]. Within the phylum Firmicutes, the Lachnospiraceae, being responsible for the production of butyrate, as well as the Ruminococceaceae and Fibrobacteraceae, with cellulolytic properties, are the most important for the proper functioning of the digestive tract of horses [6]. In the present study, the Lachnospiraceae correlates negatively with Escherichia Shigella species (p = 0.0041, r = −0.78), which may suggest the protective activity of this group. Ruminococcus flavefaciens, a representative of this phylum, shows a strong positive correlation with the amount of SIgA in the feces, which suggests their significant role in the regulation of SIgA production. Ruminococcus flavefaciens was the predominant cellulolytic bacterial species in the pony and donkey cecum [33]. The authors are aware that the evaluation of correlations at the species level with amplicon sequencing is rarely reliable; however, among the phyla, Ruminococcaceae was the dominant strain, also being relevant to the physiology of the horse's digestive tract. The second important phylum, whose role is not fully understood, is Proteobacteria, which is dominant in the equine small intestine. In the feces of the examined horses, it constituted 0.8-0.97% relative abundance [6]. Additionally, the relative prevalence of Proteobacteria moderately positively correlates with the amount of SIgA in the stool. Proteobacteria, as a group of Gram-negative bacteria often carrying LPS, is one of the possible fecal markers of intestinal homeostasis [16]. Many studies show that Proteobacteria species constitute the largest fraction of SIgA-coated microorganisms [34,35]. Additionally, intestinal Proteobacteria may influence the abundance of IgA-producing plasma cells in bone marrow and IgA in serum [35]. In the present study, however, no correlation was found between fecal Proteobacteria and serum IgA. Studies in mice have also shown that IgA specific for γ-proteobacteria can regulate the maturation of the intestinal microflora [36]. In this species, the Proteobacteria cluster will decrease with age, but in horses, the proportion of Proteobacteria (Enterobacteriaceae) increases 50 days after birth and remains unchanged in adults [7,36]. The third "key phylum" described in horses is Verrucomicrobia, which primarily colonizes the hindgut. The genus Akkermansia is described as being among those most important for the functioning of the mucosa. Its mucin-degrading action helps to maintain the integrity of the mucin layer and decrease inflammation [6]. Akkermansia mucinophila, a member of the phylum Verrucomicrobia, is highly enriched in the IgA+ fraction in humans [33]. However, in the presented study, there were a lack of correlations between Akkermansia species and SIgA in equine feces. In the V7-V9 region analysis, there was also a correlation between SIgA and Cyanobacteria (Gastranaerophilales), whose exact role in animal gut is unknown, but they are thought to have a beneficial effect for their hosts as a source of vitamins B and K [37]. However, there are no data in the literature on the significance of the other phyla, i.e., Patescibacteria and Elusimicrobia, which are also correlated with SIgA.
Diversity parameters can be used to accurately assess and analyze the variation in and composition of the gut microbiota. In the case of the presented results, the α-diversity was assessed, which made it possible to measure microbiome diversity within a single sample. To assess various aspects of community heterogeneity, the following were used: Shannon's diversity index (indirect richness), observed features (direct richness), Faith's phylogenetic diversity (diversity) and Pielou's evenness (uniformity). By analyzing V3-V4 hypervariable regions of the 16S rRNA gene, a correlation was demonstrated between indicators of evenness and richness and GALT markers. Faith's phylogenetic diversity, which takes into account the sum of the branch lengths of a phylogenetic tree connecting all species (so a higher number means more branches, more richness and a more diverse community), positively correlated with both the number of Cyathostomum eggs and SIgA.
This indicates a strong relationship between the microbiome and SIgA in the "kill zone" and the microbiome and the presence of parasites in the digestive tract of horses. Positive correlations also occurred between the amount of SIgA in the feces and the individual β-diversity metrics. β-diversity quantifies dissimilarities between communities (samples) and allowed us to determine whether the selected factors were related to the microbiome, and the obtained results suggest a clear relationship between the microbiome and GALT activity in healthy horses.
In terms of taxonomy analysis, the performance of 16S rRNA hypervariable regions varies [38]. In this study, we separately analyzed two regions of the 16S rRNA gene to estimate microbial community diversity. The resolution of these two fragments differed. On the basis of the hypervariable V3-V4 fragment, 4759 features were identified, and for the V7-V9 fragment, 3921.
The production of SCFAs by bacteria in the digestive tract of horses is closely related to their use as an energy source [30]. It is worth emphasizing that in humans, only about 5% of SCFAs produced in the large intestine is excreted in the feces, and the remaining 95% is absorbed by intestinal epithelial cells [39]. However, determination of their content in stools can provide useful information on the functioning of the microbiome. SCFAs are a kind of communicator between the microbiome and the immune system and are responsible for maintaining the balance in the anti-and pro-inflammatory reaction, including by signaling through the complex of free fatty acid receptors (GPRs) [9]. A special property of SCFAs is the induction of T-regulatory cells (Treg), which takes place through the inhibition of the enzyme histone deacetylase. The greatest inhibitory potential is seen for butyric acid, which causes proliferation and increases the functional capabilities of Treg cells [9]. Acetate is also a source of energy for regenerative and repair processes in tissues, and propionate is a precursor of gluconeogenesis. Analyzed SCFAs showed no correlation with SIgA in the feces of the tested horses. In the presented studies, the share of butyrate among all SCFAs was above 30%, which does not coincide with the results presented in the literature, where its average share does not exceed 6% [16,40]. A hay-only diet and access to pasture could be considered potential causes of these differences, but this issue requires further research. However, the VFA ratio, which is a relevant marker of fibrolytic activity [16], was about 3.473 in the presented study, which is near the normal value in the hindgut of horses fed a high-fiber diet (4.3 to 6.0) [16]. The VFA ratio correlates positively with C5 (p = 0.04, r = 0.61), which is in opposition to the results of Raspa et al., where a low level of valeric acid was associated with a fiber-rich diet [41]. However, in the presented study, the C5 value was very low in all horses. Based on data from Nadeau et al., valerate may show a pro-inflammatory effect and could be involved in gastric ulcer pathogenesis in horses; however, this cannot be compared to the current state in the gut [42]. Bacteroidota negatively correlated with iC4, which may suggest their pro-inflammatory action in horses, which was also mentioned in the work of Mach et al. They were further described to increase in number in the case of inflammatory bowel disease in humans [43,44].
It seems to be very difficult to try to answer the question of how the moderate presence of parasites influences the maintenance of homeostasis and the proper functioning of the mammalian gastrointestinal tract. There was no correlation between the degree of worming and SIgA. Studies show an increase in the α-diversity of the microbiome in the case of moderate growth parasites in both humans, non-human primates and horses (trend) [45][46][47]. In studies conducted by Peachey et al., the presence of parasites did not affect the relative frequency of phyla and the species representing the core of the microbiome [47]. A study by Walshe et al. covering the effects of anthelmintic treatment on the composition of the intestinal microbiota showed a decrease in the alpha and beta diversity of the fecal microbiota on day 7 after treatment, which reversed by day 14. These changes were accompanied by an increase in inflammatory biomarkers [48]. The use of anthelmintic drugs increases the relative abundance of Deferribacter spp. and Spirochaetes spp. and causes disturbances in carbohydrate and fiber metabolism [49]. Studies by Mach et al. indicated that Clostridiales were significantly associated with parasite egg counts, whereas Bacteroides, which is considered pro-inflammatory, was negatively associated with parasite egg counts [44]. Subsequent studies by Clark et al., 2018 showed a reduction in the phyla Ruminococcus, Clostridium and Lachnospiraceae during the egg burden in horses susceptible to infection [50]. The current research also showed a strong positive correlation between Faith's diversity and the number of strongyle eggs, and at the same time indicates a negative dependence of FEC and the frequency of Firmicutes and a negative trend between FEC and Lachnospiracea (p = 0.0712). In summary, it can be concluded that horses may have correlations between the degree of shedding (also known as the sensitivity to Cyathostomum infections) and the composition of the microflora and GI homeostasis. Additionally, the use of anthelmintic drugs affects the fecal microbiome. Considering the latest indications regarding deworming among only "high shedders" horses, the preliminary conclusions are in line with this approach [10].
The limitations of the presented paper certainly include the small number of individuals. In order to have a broader understanding of the issue, it would be appropriate to conduct a larger study, taking into account the different maintenance routines and diets in selected groups.

Conclusions
In conclusion, all the obtained results indicate the existence of significant relationships between the host response (SIgA production) and the presence of small strongyles in the digestive tract of horses. The most important factor seems to be the correlation between Faith's PD or "Key microbes" and SIgA in feces. It seems clear that the correct composition and relationship within the "kill zone" is the most important protective barrier, and its breakage may result in pathologies. In the foreseeable future, the authors plan to extend the analysis to individual parts of the gastrointestinal tract. (3) activity aimed at identifying animals; (4) capturing wild animals for biometric and systematic assessment; (5) veterinary procedures which do not cause pain, suffering, distress or permanent health impairment equal to or more invasive than the insertion of a needle. Hence, according to points 1 (blood sampling) and 5 (feces sampling) described above, the presented study did not require