Characterization of the “gut microbiota-immunity axis” and microbial lipid metabolites in atrophic and potential celiac disease

Introduction Potential celiac disease (pCD) is characterized by genetic predisposition, positive anti-endomysial and anti-tissue transglutaminase antibodies, but a normal or almost normal jejunal mucosa (e.g., minor histological abnormalities without villous atrophy). To gain further insights into basic mechanisms involved in the development of intestinal villous atrophy, we evaluated and compared the microbial, lipid, and immunological signatures of pCD and atrophic CD (aCD). Materials and methods This study included 17 aCD patients, 10 pCD patients, and 12 healthy controls (HC). Serum samples from all participants were collected to analyze free fatty acids (FFAs). Duodenal mucosa samples of aCD and pCD patients were taken to evaluate histology, tissue microbiota composition, and mucosal immune response. Results We found no significant differences in the mucosa-associated microbiota composition of pCD and aCD patients. On the other hand, in pCD patients, the overall abundance of serum FFAs showed relevant and significant differences in comparison with aCD patients and HC. In detail, compared to HC, pCD patients displayed higher levels of propionic, butyric, valeric, 2-ethylhexanoic, tetradecanoic, hexadecanoic, and octadecanoic acids. Instead, aCD patients showed increased levels of propionic, isohexanoic, and 2-ethylhexanoic acids, and a lower abundance of isovaleric and 2-methylbutyricacids when compared to HC. In addition, compared to aCD patients, pCD patients showed a higher abundance of isobutyric and octadecanoic acid. Finally, the immunological analysis of duodenal biopsy revealed a lower percentage of CD4+ T lymphocytes in pCD infiltrate compared to that observed in aCD patients. The functional characterization of T cells documented a pro-inflammatory immune response in both aCD and pCD patients, but the pCD patients showed a higher percentage of Th0/Th17 and a lower percentage of Th1/Th17. Conclusion The results of the present study show, for the first time, that the duodenal microbiota of patients with pCD does not differ substantially from that of aCD; however, serum FFAs and local T cells displayed a distinctive profile between pCD, aCD, and HC. In conclusion, our result may help to shed new light on the “gut microbiota-immunity axis,” lipid metabolites, and duodenal immune response in overt CD and pCD patients, opening new paradigms in understanding the pathogenesis behind CD progression.


Introduction
Celiac disease (CD) is an immune-mediated systemic disorder elicited, in genetically susceptible individuals, by the ingestion of gluten and related prolamines present in barley and rye. CD is a common inflammatory disease that may affect not only the small intestine but also many extra-intestinal sites, with an estimated prevalence of about 1% in the European and North American populations (Karell et al., 2003;Ciccocioppo et al., 2005). Gluten acts in concert with the HLA-DQ2 or -DQ8 genes of predisposition, the mechanistically involved autoantigen tissue transglutaminase (TG2), and the activation of adaptive immune response (CD4 + T and B cells) leading to enterocyte damage and the development of villous atrophy, the histological hallmark of CD. In particular, in CD patients, incompletely digested gliadin peptides translocate, largely via trans-epithelial transport, from the intestinal lumen to the lamina propria, where they stimulate the release of TG2. Once secreted, the enzyme catalyzes the deamidation of specific glutamine in the gliadin peptides into glutamic acid, thus creating epitopes with higher affinity for the molecules DQ2/DQ8 of major histocompatibility complex class II (MHCII). Deamidated gluten peptides bind to HLA-DQ2 or -DQ8 expressed on antigen-presenting cells (APCs) and induce the activation of CD4 + T cells, which are responsible for intraepithelial lymphocyte (IEL) infiltration, crypt hyperplasia, and villous atrophy, and the subsequent production of anti-endomysial (EMA), anti-tissue transglutaminase (tTG2), and anti-gliadin (AGA) antibodies (Jabri and Sollid, 2009;Monteleone et al., 2010;Husby et al., 2012;Rubio-Tapia et al., 2012;van Bergen et al., 2015). Gluten-specific CD4 + T cells have also been involved as the primary driver of acute cytokine release and the onset of digestive symptoms after gluten ingestion. In addition, increasing evidence indicates that the induction of a glutenspecific adaptive CD4 + T-cell response must be preceded by the activation of the innate immune system and that mast cells, key players of the innate immune response, contribute to the pathogenesis of CD [as reviewed in Anderson (2020) and Frossi et al. (2019)].
The innate and adaptive immune responses, working in synergism, lead to villous atrophy and a variable combination of clinical manifestations, from the classical syndrome of malabsorption to light or even the absent of symptoms (Ludvigsson et al., 2013).
On the other hand, potential CD (pCD) is a condition characterized by the presence of positive CD serology and genetic susceptibility, but a normal (Marsh 0) or almost normal (Marsh 1-2) jejunal mucosa without villous atrophy (Troncone et al., 1996). The expression pCD was first proposed by Ferguson (Ferguson et al., 1993) for patients who do not have and never have had a jejunal biopsy consistent with overt CD and yet have immunological abnormalities similar to those found in CD. In detail, features that are good candidate markers of pCD include the presence of serum endomysial antibodies, a high count of villous intraepithelial lymphocytes, an increased density of IEL expressing γδ T-cell receptor, and signs of activated mucosal cell-mediated immunity, such as expression of CD25 and B7 by lamina propria mononuclear cells. Moreover, pCD is also characterized by enhanced expression of MHC class II molecules on the epithelium and on adhesion molecules in the lamina propria, all of which are enhanced by in vitro gluten challenge (Maiuri et al., 1996), a coeliac-like intestinal antibody pattern, and a positive rectal gluten challenge (Troncone et al., 1996).
This particular condition is usually considered a clinical challenge because, though pCD represents the early CD stage and its diagnostic criteria are clear, many questions are still unsettled and the results of the studies conducted so far are still conflicting (Kurppa et al., 2009;Tosco et al., 2011;Lionetti et al., 2012;Biagi et al., 2013;Zanini et al., 2013;Auricchio et al., 2014).
In addition to the involvement of the immune system, recent evidence reported several modifications in the intestinal tract microenvironment of CD patients, suggesting a role for gut microbiota in CD onset and persistence (Abdukhakimova et al., 2021). It is known that the gut bacteria play many fundamental roles for the host, such as the synthesis of many nutrients and metabolites (Gibiino et al., 2017;Thursby and Juge, 2017;Valdes et al., 2018), maintenance of the intestinal epithelial integrity (Khosravi and Mazmanian, 2013), and notably, the modulation of immune responses (Brestoff and Artis, 2013;Woo and Alenghat, 2017).
Furthermore, the production of short-chain fatty acids (SCFAs), the major bacterial fermentative end products, reflects the intestinal microbiota composition and activity. In particular, given that SCFAs are crucial to maintaining the host's normal gut physiology and metabolic functions and that a part of them enters the systemic circulation (den Besten et al., 2013), it is tempting to speculate that the gut microbiota could be also related to free fatty acid (FFA) levels. FFAs, classified as short-chain fatty acids (SCFAs), medium-chain fatty acids (MCFAs), and long-chain fatty acids (LCFAs), derive from microbial and host metabolism and, especially MCFAs and LCFAs, are introduced into the diet (i.e., milk and dairy products). Since FFAs modulate the production of chemokines and cytokines (Frommer et al., 2015;Honda et al., 2015;Hung et al., 2015), an altered FFA profile has been associated with the risk of developing a range of disorders in which the immune system is involved, including CD (Aghdassi et al., 2011;Nishi et al., 2014;van Hees et al., 2014;Dai et al., 2015;Rodriguez-Carrio et al., 2016). Therefore, given that pCD is a valuable biological model of the pathway leading to small intestinal mucosal damage in genetically predisposed subjects, our present study aimed to evaluate whether microbial, lipid, and immunological signatures could better characterize pCD from the overt CD condition. In addition, we also evaluated whether and how immunological peculiarity affects the composition of the duodenal microbiota and its metabolic activity, opening new paradigms in understanding the basic mechanisms involved in the development of small intestinal villous atrophy.

Patient recruitment
Twenty-seven CD patients with positive IgA antiendomysial (EMA) and anti-tissue transglutaminase (tTGA) antibodies consecutively observed at the Tuscany Regional Referral Center for adult CD were enrolled in this observational study from January 2018 to December 2019. Twelve healthy controls, with negative EMA and tTGA antibodies, and without any other clinical problems, were enrolled among the internal staff participating in the study.
Inclusion criteria were as follows: Age between 18 and 70 years, absence of any form of immunodeficiency, in particular selective IgA deficiency, and patients resident in Tuscany who, in the past 5 years, have not made any trips to countries outside Europe. This last criterion is very relevant for the study of the composition of intestinal microbiota. All patients recruited received adequate information on the study, the objectives, and how they will have to provide the organic samples for analysis. Each patient signed and dated the informed consent.
Exclusion criteria were as follows: Treatment with antibiotics or probiotics during the previous 2 months, acute gastrointestinal infections 1 month before the enrolment, pregnant and breastfeeding women (ongoing or scheduled for the next 48 weeks), the concomitant presence of established malignant neoplasms or chronic inflammatory bowel diseases (Crohn's disease and ulcerative colitis), and patients who have used immunosuppressive drugs in the previous 3 months.

Patient classification
According to mucosal histology, patients were divided into two groups: aCD (5 men, 12 women, mean age 35.6 years, age range 17-56 years) and pCD (two men, 8 women, mean age 38.8 years, age range 21-50 years). In agreement with the definition of Troncone (Troncone et al., 1996), we decided to include cases with Marsh 2 in the pCD group, given the absence of villous atrophy in these patients. Moreover, 12 healthy controls (five men, seven women, mean age 33.6 years, age range 23-50 years), with negative EMA and tTGA antibodies, served as a control group (HC). For ethical reasons, esophagogastroduodenoscopy with small bowel biopsies was not performed in this group. All the study participants were on a gluten-containing diet and did not consume any drug at the time of examination.

Sample collection
After an overnight fast, venous blood samples were collected, and three serum aliquots per patient were immediately frozen and stored at -80 • C until further use. All patients underwent an esophagogastroduodenoscopy, and 4-6 small bowel biopsies were taken from the distal part of the duodenum. Small intestinal mucosal damage was graduated according to the classification of Marsh modified by Oberhuber et al. (1999).

Antibody testing
The tTGA levels were measured by a commercially available enzyme-linked immunosorbent assay kit (EutTG, Eurospital, Trieste, Italy) that employs human recombinant tTG as antigen. EMA antibodies were determined by indirect immunofluorescence, using tissue sections of monkey esophagus (Eurospital) as previously reported (Vignoli et al., 2019).
Deoxyribonucleic acid extraction and bioinformatics analysis of 16S ribosomal ribonucleic acid DNA extraction and bioinformatics analysis of 16S rRNA were performed as previously described , Briefly, total DNA was extracted from duodenal mucosa samples of aCD and pCD patients using the DNeasyPowerLyzerPowerSoil Kit (Qiagen, Hilden, Germany) from frozen samples (-80 • C) according to the manufacturer's instructions. The quality and quantity of extracted DNA were assessed using the Qubit Fluorometer (Thermo Fisher Scientific) and then genomic DNA was frozen at -20 • C.
Extracted DNA samples were sent to NEXT Genomics (Sesto Fiorentino, Italy) where amplicons of the variable V3-V4 region of the bacterial 16s rRNA gene were sequenced in paired-end mode (2 × 300 cycles) on the Illumina MiSeq platform, according to the Illumina 16S Metagenomic Sequencing Library Preparation protocol (Pagliai et al., 2020;Russo et al., 2021). Raw sequences were processed following the software pipeline MICCA (Albanese et al., 2015). Pairedend reads were assembled using the "mergepairs" command, maintaining a minimum overlap of 100 bp and an edit distance in the maximum overlap of 32 bp. Subsequently, the sequences were cut with the "trim" command in order to remove the primers and eventually eliminate the reads with imperfections in primer sequences. All the reads with a length lower than 350 bp and with an error rate higher than or equal to 0.5 were removed with the "filter" command. Cleaned reads were eventually merged into a single file with the "merge" command and transformed into a fasta file. The OTUs were generated using the "otu" command in "denovo_greedy" mode, setting a 97% identity and performing an automatic removal of chimeras with the "-c" option. The longest sequence of each OTU was used for the taxonomic assignment using the "classify" command in "rdp" mode, i.e., using the RDP Bayesian classifier (10.1093/nar/gki038) that is able to obtain classification and confidence for taxonomic ranks up to genus level.

Microbial community analysis
Statistical analyses on the bacterial community were performed in R (R Development Core Team, 2014) with the help of the packages phyloseq 1.26.1 (McMurdie and Holmes, 2013), DESeq2 1.22.2 (Love et al., 2014), and other packages satisfying their dependencies, in particular vegan 2.5-5 (Willis and Bunge, 2015). Rarefaction analysis on OTUs was performed using the function rarecurve (step 50 reads), further processed to highlight saturated samples (arbitrarily defined as saturated samples with a final slope in the rarefaction curve with an increment in OTU number per reads < 1e-2). For the cluster analysis (complete clustering on Euclidean distance) of the entire community, the OTU table was first normalized using the total OTU counts of each sample and then adjusted using square root transformation. The coverage was calculated by Good's estimator using the formula: (1 -n/N) × 100, where n is the number of sequences found once in a sample (singletons), and N is the total number of sequences in that sample.
Richness, Shannon, Chao 1, and evenness indices were used to characterize the ecological properties of each sample using the function estimate_richness from phyloseq (McMurdie and Holmes, 2013). The evenness index was calculated using the formula E = S/log(R), where S is the Shannon diversity index and R is the number of OTUs in the sample. Differences in all indices between overt CD and pCD were tested using a paired Wilcoxon signed-rank test. The differential analysis of abundance at the OTUs, as well as at the different taxonomic ranks (created using the tax_glom function in phyloseq), was performed with DESeq2 (Love et al., 2014).

Evaluation of free fatty acids by gas chromatography-mass spectrometry analysis
The analysis of FFAs was performed by Agilent gas chromatography-mass spectrometry (GC-MS) system composed of 5971 single quadrupole mass spectrometer, 5890 gas-chromatograph, and 7673 autosampler, with a dedicated previously described protocol (Baldi et al., 2021). The chemicals, GC-MS conditions, and calibration parameters were reported in supporting information SX. Just before the analysis, each sample was thawed. The FFAs were extracted as follows: an aliquot of 300 µl of serum sample was added to 10 µl of ISTD mixture, 100 µl of tert-butyl methyl ether, and 20 µl of 6 M HCl + 0.5 M NaCl solution in 0.5 ml centrifuge tube. Afterward, each tube was stirred in a vortex for 2 min, Frontiers in Microbiology 05 frontiersin.org centrifuged at 10,000 rpm for 5 min, and finally, the solvent layer was transferred to a vial with a microvolume insert and analyzed.

Statistical analysis
The statistical analysis of immunological and FFAs data was performed using GraphPad Prism (v. 5) software, and the data were expressed as the mean ± standard deviation (SD). The comparisons between aCD, pCD patients, and HC were evaluated with the Kruskal-Wallis test, while the comparisons between the Tcc groups were evaluated with the two-tailed Chisquare test or Fisher's exact test, when appropriate (in detail, when ≤ 20% of cell counts were less than 5, we used the Chisquare test; if > 20% of expected cell counts were less than 5, we used Fisher's exact test). P-values < 0.05 were considered statistically significant.

Characteristics of the study population
Patients were divided into two groups: aCD (5 men, 12 women, mean age 35.6 years, age range 17-56 years) and pCD (2 men, 8 women, mean age 38.8 years, age range 21-50 years), in accordance with mucosal histology. In more detail, nine aCD patients had partial villous atrophy (Marsh 3 A), six subtotal villous atrophy (Marsh 3B), and two total villous atrophy (Marsh 3C); among pCD patients, histological examination revealed an apparently normal mucosa (Marsh 0) in three subjects, just an increase in intraepithelial lymphocytes (Marsh 1) in five, and increased intraepithelial lymphocytes coupled with crypt hyperplasia (Marsh 2) in two patients ( Table 1). Twelve healthy controls (5 men, 7 women, mean age 33.6 years, age range 23-50 years), with negative EMA and tTGA antibodies, served as a control group (HC).
Analysis of microbiota composition and free fatty acids of atrophic CD and potential celiac disease groups Comparison of tissue microbiota composition between atrophic CD and potential celiac disease patients We first analyzed the microbiota composition of duodenal mucosa samples obtained for both CD and pCD patients. We obtained a total of 737,729 reads, and after all preprocessing steps (pair merging, trimming, quality filtering, and chimeras removal), a total of 523,825 reads were available for further analysis (  Taxonomic composition of aCD and pCD patients' microbiota Coordinate plot showing the relative abundance of the five most abundant phyla in each aCD and pCD duodenal biopsy. As reported in Figure 1, the analysis of the taxonomic composition revealed that more than 99% of the sequences collected could be ascribed to the five most representative phyla: Proteobacteria (72.98%), Firmicutes (11.92%), Bacteroidetes (10.36%), Actinobacteria (2.33%), and Fusobacteria (2.16%). As shown in Figure 2, no significant difference was reported for richness and Shannon and Chao1 alpha diversity indices. Moreover, in order to investigate the similarity of patients in terms of taxonomy abundance profiles, we evaluated PCoA and NMDS and completed the hierarchical clustering using several data transformations and the Bray-Curtis dissimilarity as a distance metric. Surprisingly and despite our efforts, no evident bacterial groupings were observable between aCD and pCD patients. As evidence of our fine analysis, in Figure 3, we reported several multivariate plots built upon square root transformed percent abundances, namely, a) hierarchical clustering of Euclidean distance, b) hierarchical clustering of the top five most abundant taxa only, and c) a PCoA of Bray-Curtis distance.

Evaluation of serum free fatty acids profiles
We evaluated the metabolic profile of the patients performing the qualitative and quantitative analysis of serum FFAs, namely, linear SCFAs (acetic, propionic, butyric, and Boxplots reporting alpha diversity indices (respectively, Richness, Shannon index, and Chao1 index) in aCD and pCD samples. Alpha diversity indexes are composite indexes reflecting abundance and consistency. Boxes represent the interquartile range (IQR) between the first and third quartiles (25th and 75th percentiles, respectively), and the horizontal line inside the box defines the median. Whiskers represent the lowest and highest values within 1.5 times the IQR from the first and third quartiles, respectively. P-values less than 0.05 were considered statistically significant.

Characterization of mucosal infiltrating T cells in atrophic CD and potential celiac disease patients
We evaluated differences in the T-cell infiltration using biopsies of the duodenal mucosa of aCD and pCD patients; for ethical reasons, we were not able to obtain duodenal biopsies from HC to be used to compare the results.
We cloned and expanded in vivo-activated TILs, obtaining a total number of 267 T-cell clones: 206 Tcc from aCD patients and 61 Tcc from pCD patients.

Discussion
In this original study, we show for the first time that the duodenal microbiota of patients with pCD does not differ  Boxplots representing the overall abundance of FFAs in aCD patients, pCD patients, and healthy controls (µmol/L). P-values of the intergroup comparisons were assessed with Kruskal-Wallis test. P-values less than 0.05 were considered statistically significant. The asterisks * represent p-values, *p < 0.05, **p < 0.01, ***p < 0.001. substantially from that of aCD. On the other hand, we observed differences in the local immune T-cell population (increased Th0/Th17 and reduced Th1/Th17 in pCD), while serum FFAs displayed a distinctive profile between pCD, aCD, and HC.
Although great advances have been made in understanding the role of CD adaptive immune mechanisms in response to gluten peptides, and their crosstalk with the intestinal microbiota, the complete set of pathogenic events causing the development of tissue lesions remains unclear. In our study, therefore, we explored for the first time the complex crosstalk between the microbiota, lipid metabolites, and the T-cell profile in both aCD and pCD patients, to better characterize aspects still unclear in the pCD condition and assess further significant differences between the two faces of celiac disease.
First, we compared the mucosal microbiota composition in the duodenal biopsies of aCD and pCD patients. Since most studies have been conducted using both different sequencing methods and inhomogeneous sampling, the findings on the gut microbiota composition of aCD patients are highly heterogeneous and sometimes contrasting.
were Firmicutes, Bacteroidetes, and Actinobacteria, whereas CD children displayed higher abundances of Proteobacteria, Bacteroidetes, and Actinobacteria. In accordance with these findings, we observed an abundance of Firmicutes, Bacteroidetes, and Actinobacteria in both CD and pCD patients; however, in contrast with the previous results, we observed high percentages of Proteobacteria (72.98%) and Fusobacteria in adult CD and pCD patients. The high abundance of Gram-negative Proteobacteria may result in altered exposure of intestinal epithelial cells to bacterial lipopolysaccharides (LPS) and some metabolites, such as SCFAs (Akobeng et al., 2020). Intriguingly, no statistically different microbial taxa or bacterial groupings were displayed between the two patients' cohorts, suggesting the existence of a similar microbial consortium in both conditions. Even though we demonstrated the presence of a stable microbiota composition between overt CD and pCD patients, based on our previous metabolic results (Bertini et al., 2009;Bernini et al., 2011), we hypothesized that the functional metabolic activity of gut microbiota could be differently regulated between the two pathologic conditions. Indeed, in one of our previous studies, by employing a non-invasive metabonomic approach, we have first shown that CD has a welldefined metabonomic signature (Bertini et al., 2009). Notably, sera of CD patients were characterized by lower levels of several metabolites, such as amino acids, lipids, pyruvate, and choline, and by higher levels of glucose and 3-hydroxybutyric acid, while urine samples showed altered levels of, among others, indoxyl sulfate, meta- [hydroxyphenyl]propionic acid, and phenylacetylglycine. Moreover, in another study addressed to pCD and CD patients as well as HC, we observed, examining serum and urine samples, that pCD largely shares the metabonomic signature of aCD (Bernini et al., 2011). Most metabolites found to be significantly different between control and CD subjects were also altered in pCD (Bernini et al., 2011).
Regarding the free fatty acids, altered circulating levels of SCFAs, MCFAs, and LCFAs and LCFAs, have been related to CD (Niccolai et al., 2017;Baldi et al., 2021). In our previous study, we compared the serum FFA profile of patients with CD, adenomatous polyposis (AP), and colorectal cancer (CRC) to HC (Baldi et al., 2021). We observed that HC showed a different composition of FFAs as compared to CRC, AP, and CD patients. Moreover, we demonstrated a strong positive association between CD and butyric acid, suggesting for the first time that it could represent a potential biomarker for CD screening. However, no information addressing the serum SCFA composition of pCD patients is currently available. So, as they provide information on the general health status of both the host and its intestinal microbial ecosystem, we evaluated the abundance of each circulating FFA. In detail, compared to HC, pCD patients displayed higher levels of propionic, butyric, valeric, 2-ethylhexanoic, tetradecanoic, hexadecanoic, and octadecanoic acids. Propionic, butyric, and valeric acids are the most common microbiota-produced SCFAs that, as well established, exert positive metabolic effects on enterocytes, act as tumor suppressor agents, and regulate host gene expression via histone deacetylase inhibition (HDAC; Nankova et al., 2014;Yuille et al., 2018).
In contrast to SCFAs, MCFAs and LCFAs arise mostly from dietary triglycerides and are important regulators of energy metabolism, gene expression, and immune processes (Schonfeld and Wojtczak, 2016;Pinkosky et al., 2020). In particular, they can antagonize the anti-inflammatory activities of SCFAs, as they support the Th1 and Th17 differentiation, so both MCFAs and LCFAs could be involved in the onset and persistence of duodenal inflammation in pCD patients, compared to HC.
In aCD patients, the results of our study showed increased levels of propionic, isohexanoic, and 2-ethylhexanoic acids, and a lower abundance of isovaleric and 2-methylbutyric acids as compared to HC. In addition, we found that, compared to aCD patients, pCD patients showed higher abundances of isobutyric and octadecanoic acid. Isobutyric acid is a branched SCFA produced by valine and leucine Percentage of CD4 + and CD8 + T-cell clones, respectively, obtained from the duodenal mucosa of patients with atrophic CD (A) and potential CD (B). (A) The cytokine phenotype percentage of CD4 + T-cell clones obtained from the duodenal mucosa of patients with atrophic CD (orange) and potential CD (blue). P-values of the intergroup comparisons were assessed with Fisher's exact test. P-values less than 0.05 were considered statistically significant. The asterisks * represent p-values, *p < 0.05. (B) The cytokine phenotype percentage distribution of CD8 + T-cell clones obtained from the duodenal mucosa of patients with atrophic CD (orange) and potential CD (blue). fermentation, and a high abundance of this acid indicates an increased amino acid fermentation. So, this finding suggests that, even if the composition of the duodenal microbiota is almost similar, aCD and pCD patients displayed a different microbial functional activity, which deserves to be deepened in future studies.
As there is a documented mutual interplay between the microbiome and the immune response, to assess whether and how immunological peculiarity affects the composition of the duodenal microbiota and its metabolic activity, we also evaluated the cellular immune profile in aCD and pCD patients. In detail, the phenotypic characterization of T cells revealed, in agreement with previous studies, a lower percentage of CD4 + T cells in pCD patients compared to CD patients. It has been shown that CD4 + T cells play a key role in the inflammatory response triggered by gluten, by driving and orchestrating the gut inflammation and tissue destruction in human celiac disease (Voisine and Abadie, 2021).
Moreover, in accordance with Vitale et al., pCD patients showed a high percentage of T cells co-producing IL-17, IL-4, and IFN-γ (Th0/Th17 p < 0.00005, Tc0/Tc17; Vitale et al., 2019). Notably, IL-17 is a pro-inflammatory cytokine, but a number of studies underlined its potential involvement in the prevention of gut mucosa inflammation (Mucida and Salek-Ardakani, 2009;Cukrowska et al., 2017). In addition, as compared to atrophic CD, pCD patients showed lower percentages of T cells co-producing IFN-γ and IL-17 (Th1/Th17 p < 0.001). Of note, the Th1/Th17 subpopulation promotes damage to the enterocytes leading to the atrophy of the intestinal villi (Cukrowska et al., 2017). We hypothesize that the increase of Th0/Th17 linked with the decrease of Th1/Th17 might contrast the inflammation, contributing to preventing mucosal damage in the pCD patients. However, whether and how this immunological peculiarity can affect the composition of the duodenal microbiota and its metabolic activity remains unknown.
Therefore, although our exploratory study has some limitations, such as a limited number of patients, little information on patients' dietary pattern, or the relatively low level of taxonomic resolution of 16S rRNA gene metagenomics, we obtained innovative and interesting preliminary data, such as the same duodenal microbial architecture between aCD and pCD patients but differences in the serum FFA's profiles and in the duodenal adaptive immunity. However, we will try to overcome the mentioned limits in a future larger and more detailed investigation with aCD and pCD patients, divided into homogeneous groups considering also age, gender, nutrition, and comorbidity.
Although the precise relationship between host lipid metabolism and gut microbiota remains to be clarified, our findings suggest a possible modification in microbial functionality that needs to be further explored assessing also the contribution of the colonic microbiota, in order to detect potential new biomarkers of human pCD.

Data availability statement
The data presented in this study are deposited in the NCBI Gene Expression Omnibus repository, accession number: GSE181070 and the analysis script is available at https://github. com/matteoramazzotti/papers/tree/main/2021celiac.

Ethics statement
The studies involving human participants were reviewed and approved by Ethical Committee of the University of Florence, Italy (CE: 10443_oss, 14/02/2017). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions
ASC, AA, and ER designed the study. FR, GN, and GL collected the samples. FR, DR, SB, MP, MM, and GB performed the experiments. FR, ER, SB, and EN analyzed the data. MC and MR analyzed microbiota data. ER and FR wrote the manuscript. ER edited the manuscript. AA, DR, and ASC supervised the manuscript. ER, AA, and ASC provided for funding acquisition. All authors have approved the final draft submitted.

Funding
This work was supported by FC (Fondazione Celiachia Onlus) Grant Fellowship n • 007_FC_2016. Moreover, this project received funding from MUR under the umbrella of the European Joint Program Initiative "A Healthy Diet for a Healthy Life" (JPI-HDHL) and of the ERA-NET Cofund ERA-HDHL, ID: 1523 (GA N 696295 of the EU HORIZON 2020 Research and Innovation Programme).