Our most important findings are the higher Blastocystis prevalence in the UNEME cohort than in the FACSA cohort (Table 2). The faecal-oral transmission of Blastocystis is due to poor hygiene practices, exposure to animals infected with the parasite, and intake of contaminated water or food (17). This transmission route is probably most common since the most UNEME cohort subjects live in rural areas (Table 1).
The F/B ratio comparison between cohorts showed dysbiosis in UNEME cohort subjects (Table 2) according to a high F/B ratio (Median of 1.40, IR: 0.32-4.03) compared to the subjects of the FACSA cohort (Median of 0.830, IR: 0.07-3.50, p <0.05). This dysbiosis agrees with previous reports, in which a high F/B ratio is found in metabolically ill subjects (4), this relationship implies a predisposition to disease states (29). Likewise, a low Firmicutes/Bacteroidetes ratio is related to weight loss (30), which corresponds with FACSA cohort individuals with lower BMI in this study (Table 1).
Obese and diabetic individuals, compared to healthy individuals, have a higher relative abundance of Firmicutes and a reduced abundance of Bacteroidetes (10), as well as a low microbial gene count and a dominance in the genera Bacteroides and Ruminococcus, all this is associated with a more remarkable ability to obtain energy from the diet, systemic inflammation, adiposity, insulin resistance, and dyslipidemia (31,32). Bacteroidetes are known to produce mainly acetate and propionate, while Firmicutes produce more butyrate, attributed to anti-inflammatory activities, regulation of energy metabolism, and increases leptin.
The increase in Firmicutes in metabolically ill patients could cause an increase in the butyrate production leading to antiobesogenic effects, which is contradictory. It has been speculated that in obese subjects, the butyrate-producing bacteria decrease and are replaced by other bacteria belonging to the same phylum (33).
Blastocystis is a widely distributed organism with great adaptability that could colonize healthy and diseased subjects. We analyzed the Firmicutes/Bacteroidetes relationship and the presence of Blastocystis in these two cohorts, and the results showed an association with a low F/B ratio (Table 3). Some studies have suggested that the intestine's dysbiosis prevents colonization by anaerobic eukaryote Blastocystis. The decreases in butyrate available for the colonocytes' oxidative metabolism allow the increase of the luminal oxygen concentrations and, therefore, the proliferation of pathogenic bacteria. Thus, this environment would cause oxidative stress on Blastocystis, and therefore its survival would be affected (34,35).
Surprisingly our findings suggest a possible adaptation of Blastocystis to oxidative stress and low-grade inflammation that has been observed in metabolically ill patients, this low-grade inflammation was found in the UNEME cohort but not in the FACSA cohort (Table 1). Although we do not have data of interleukin expression, the data on the white formula showed higher levels of leukocytes, lymphocytes, neutrophils, and platelets in UNEME subjects, which suggest low-grade inflammation (36,37). Besides high Blastocystis prevalence in both cohorts, these results are unlikely caused by the protozoa. The adaptation of Blastocystis may have been due to its alternative oxidase mechanism (OXA) (38) and the five families of genes acquired by lateral gene transfer involved in response to oxidative stress in Blastocystis (39). Therefore, an essential role for Blastocystis is suggested in terms of this low F/B ratio, which may have a beneficial effect on obese individuals in the UNEME Cohort.
The relationship of Blastocystis with the gut microbiota is a subject of debate, as it has been linked to a low F/B ratio and irritable bowel disease. However, Audebert et al. (2016) found a greater abundance of Clostridiales at the class level and a greater abundance of Rumminococcaceae and Prevotellaaceae at the family level in subjects with Blastocystis, while Enterobacteriaceae increased in patients without Blastocystis. It has been suggested that Blastocystis is not associated with dysbiosis observed in intestinal, metabolic diseases or infections commonly associated with inflammation of the lower gastrointestinal tract, instead, colonization by this parasite could be associated with a healthy intestinal microbiota (18). Although still controversial, Defaye et al. (2020) reported the decrease in the Firmicutes/Bacteroidetes ratio (p = 0.06) (40) in a model of IBS (Irritable Bowel Syndrome) in rats infected with Blastocystis ST4 from healthy humans. This partially agrees with our findings since in the FACSA cohort ST2 had an OR = 3.12 (0.61-15.88 95% CI), ST3 OR = 1.97 (0.86-4.52 95% CI), and ST4 OR = 2.34 (0.78-7.02 95% CI), p> 0.05 (Table 4) are not associated with Low F/B ratio. However, there was an association with low F/B at the gender level OR = 3.78 (2.10-6.81 CI 95%), for the subtypes ST1 OR = 3.99 (1.07-14.79 CI 95%) and ST7 OR = 5.44 (1.16- 25.52 95% CI) p <0.05, (Table 4). Regarding ST7, Yason et al. (2019) reported that the presence of ST7 decreases the Bifidobacterium and Lactobacillus populations while increasing the Escherichia populations (41).
Regarding ST1, the presence of an alternative oxidase provides a partially dependent metabolism of molecular oxygen to resist the stress that the high oxygen concentration entails, which may occur in the intestine of the FACSA subjects due to the low F/B ratio (39). Also, in the UNEME cohort, a low F/B ratio was found concerning gender OR = 4.24 p <0.05 (Table 5) but not for subtypes.
The low F/B ratio found in this study might imply that the subjects of the two cohorts infected by Blastocystis present gastrointestinal symptoms, however, the analysis of the association of the Blastocystis prevalence with gastrointestinal symptoms showed an inverse association between abdominal pain and ST1. At the same time, ST4 was inversely associated with abdominal distension in the FACSA cohort in previous data published (42) (Table S2 and S3 of the supplement), however, no association of Blastocystis with gastrointestinal symptoms was found in the UNEME cohort (Table S4 and S5 of the supplement). Contrary to the results obtained, a study found that Blastocystis negatively correlates with Bacteroidetes (18), while the phylum Firmicutes presents a positive correlation in Blastocystis positive samples (21). These results could be biased in type 2 diabetes subjects due to the low carbohydrate consumption. Since Firmicutes have more genes for the enzymes involved in their metabolism, a decrease in carbohydrate consumption could leads to Firmicutes decrease and Bacteroidetes increase (10,43), however, we did not analyse diet in this study.
Possibly, dysbiosis does not always lead to inflammation or disease, but certain conditions that lead to an inflammatory state would have to be present to affect the individual. For example, in obese individuals, the proportion of Firmicutes and Proteobacteria increases, compared to Bacteroidetes, in this case, an inflammatory environment is present (44). In the UNEME cohort, we observed dysbiosis since the Firmicutes/Bacteroidetes ratio increased in the UAR of Bacteroidetes, but this does not necessarily imply inflammation. Additionally, there is evidence that Blastocystis modulates the immune system through IL-22 release that stimulates the mucus production, alleviates colitis symptoms (45) and induces an immune response with a predominance of the Th2 cell response, favouring an anti-inflammatory environment (46).
Regarding ST4, although the prevalence was not high, it was more frequent in FACSA cohort (clinically healthy subjects) subjects with lower BMI (Table 1 and 2) and subjects with lower BMI in the total sample (Table S6 supplement). These findings agree with Beghini et al. (2017) that found a strong negative correlation between BMI and Blastocystis prevalence (21). Also, consistent with findings from the Danish subjects study (20), the difference in Blastocystis prevalence between average weight and obese subjects (p = 5E-03), average weight and overweight (p = 0.01), and between non-overweight and overweight (p = 0.02) was significant. Between specific subtypes, only ST4 reached statistical significance (p = 0.03 between average weight and obese). Besides, Tito et al. (2019) found a positive and significant correlation (R = 0.26 p = 0.00028) between ST4 and Akkermansia and Methanobrevibacter (19), the first is an abundant bacterium in healthy people that degrades intestinal mucin, which is associated with weight loss, the second is a methanogenic archaeon that plays an essential role in carbohydrate digestion and may protect against weight gain (31).
ST3 was the most prevalent in the FACSA cohort, 29.79%, and the second most prevalent, 25.93%, in the UNEME cohort (Table 2). This subtype was not associated with intestinal dysbiosis in either of the two cohorts (Table 4 and Table 5), in agreement with previous reports (48,56,57). It has been reported a higher bacterial diversity in ST3-Blastocystis-carriers (high abundance of Prevotella, Methanobrevibacter, and Ruminococcus), while a high percentage of Bacteroides found in Blastocystis-free subjects (57). Asnicar F. et al. (2021) reported interesting findings of the presence of Prevotella copri and Blastocystis spp as markers of improved postprandial glucose response, both were strongly linked with favourable glucose homeostasis and a decrease of the estimated visceral adipose tissue mass (58).
In future studies, we will analyse the dietary habits and the composition of their microbiota, including Blastocystis, in both cohorts. One of the limitations of our work was qualitative PCR, which only identifies the presence or absence of ST. The implementation of a more sensitive molecular technique, Next generation amplicon sequencing like the one used by Maloney (2019), could give us additional information, such as the most predominant ST or the existence of more than two subtypes in a single sample (59). It could also help with the identification of genotypes that were not detected with the primers used. A more accurate assessment of Blastocystis diversity is the key to understand the transmission mechanism and its pathogenicity in our population. Another limitation was that the majority of the subjects in the UNEME cohort were obese, and the analysis of the comparison with thin individuals between the two cohorts could not be carried out. Therefore, the objective is to increase the sample size in this cohort concerning this group.