Gut dysbiosis in cutaneous T-cell lymphoma is characterized by shifts in relative abundances of specific bacterial taxa and decreased diversity in more advanced disease.

BACKGROUND: Cutaneous T-cell lymphoma (CTCL) patients often suffer from recurrent skin infections and profound immune dysregulation in advanced disease. The gut microbiome has been recognized to influence cancers and cutaneous conditions; however, it has not yet been studied in CTCL.
OBJECTIVES: To investigate the gut microbiome in patients with CTCL and in healthy controls.
METHODS: A case-control study was conducted between January 2019 and November 2020 at Northwestern's busy multidisciplinary CTCL clinic (Chicago, Illinois, USA) utilizing 16S ribosomal RNA gene amplicon sequencing and bioinformatics analyses to characterize the microbiota present in fecal samples of CTCL patients (n = 38) and age-matched healthy controls (n = 13) from the same geographical region.
RESULTS: Gut microbial α-diversity trended lower in patients with CTCL and was significantly lower in patients with advanced CTCL relative to controls (P = 0.015). No differences in β-diversity were identified. Specific taxa were significantly reduced in patient samples; significance was determined using adjusted P-values (q-values) that accounted for a false discovery rate threshold of 0.05. Significantly reduced taxa in patient samples included the phylum Actinobacteria (q = 0.0002), classes Coriobacteriia (q = 0.002) and Actinobacteria (q = 0.03), order Coriobacteriales (q = 0.003), and genus Anaerotruncus (q = 0.01). The families Eggerthellaceae (q = 0.0007) and Lactobacillaceae (q = 0.02) were significantly reduced in patients with high skin disease burden.
CONCLUSIONS: Gut dysbiosis can be seen in patients with CTCL compared to healthy controls and is pronounced in more advanced CTCL. The taxonomic shifts associated with CTCL are similar to those previously reported in atopic dermatitis and opposite those of psoriasis, suggesting microbial parallels to the immune profile and skin barrier differences between these conditions. These findings may suggest new microbial disease biomarkers and reveal a new angle for intervention.

Introduction

Cutaneous T‐cell lymphoma (CTCL) comprises a heterogeneous group of non‐Hodgkin’s lymphomas involving skin‐homing malignant T‐cells. While CTCL remains to be robustly understood, it is likely influenced by both host and environmental factors. Strong clinical evidence connects CTCL and the microbial world. Patients with advanced CTCL suffer from significant morbidity secondary to recurrent infections and individuals with frequent infections tend to have more advanced disease that is less responsive to CTCL therapies. , , Moreover, prolonged courses of broad‐spectrum antibiotics have been shown to reduce malignant cytokine activity and tumor burden in advanced‐stage CTCL patients. Concomitant with these phenomena is profound immune dysregulation that may be antecedent to or a consequence of pathogenic microbial activity. As such, CTCL is believed to foster global microbial dysbiosis. The microbiome is an emerging focus within this field; however, the gut microbiota of CTCL have yet to be characterized.

Gut dysbiosis has been associated with cancer and inflammatory skin disease. Studies have demonstrated the tumorigenic potential of certain bacterial taxa and their powerful immunomodulatory abilities. Immune dysfunction can also alter the gastrointestinal microcosm to allow or stimulate the growth of virulent bacteria. This finding suggests augmented gut dysbiosis may simultaneously encourage and reflect more severe immune dysfunction. Gut dysbiosis is also understood to promote systemic disease through cytokine‐induced inflammation, aberrant effector T‐cell activation, and gut epithelial barrier disruptions that result in bacterial translocation. , , Increased gastrointestinal abundance of pro‐inflammatory, immune‐sensitizing species like Ruminococcus gnavus and loss of anti‐inflammatory, protective species such as Faecalibacterium may be linked to the dysregulated cytokine signatures characteristic of atopic dermatitis (AD), , psoriasis, , , and hidradenitis suppurativa. Because gut dysbiosis has been demonstrated in other inflammatory skin conditions and cancers, we hypothesize alterations in the gut microbiome may also be associated with CTCL disease progression.

To better understand the gut microbiome of CTCL, we conducted a cross‐sectional analysis of the microbiota present in stool samples from CTCL patients and healthy controls.

Methods

Ethical approval was obtained from the Northwestern University Institutional Review Board (STU00209226). Personal data and stool samples were collected from the Northwestern University Cutaneous Lymphoma specialty clinic between 2019 and 2020 in compliance with the Declaration of Helsinki. Informed consent and HIPAA Authorization for Research were obtained from all participants prior to study enrollment. Each patient had clinically‐ and biopsy‐confirmed CTCL, as reviewed by an expert dermatopathologist (JG).

Participants

Of CTCL patients (n = 38), 27 patients had been diagnosed with mycosis fungoides (MF), 5 with Sézary syndrome (SS), and 6 with non‐MF/SS CTCL (Table 1, Table S1). Patients were untreated or receiving standard‐of‐care therapies, including skin‐directed (n = 18, 82%) and select systemic treatments (n = 10, 26%) (Table S2). Patients were excluded if they used any form of antibiotics within the preceding 4 weeks. Modified Severity‐Weighted Assessment Tool (mSWAT) was assessed by the principal investigator (XAZ). The HC group (n = 13) was composed of age‐matched volunteers without CTCL or other active skin diseases from the same geographical region. Statistical analyses were completed using STATA SE (College Station, TX, USA).

Table 1

Characteristics of patients (n = 38) and healthy controls (n = 13)

	Patients	Controls	P‐value
N	38	13
Gender*
Male	27 (71.0)	7 (53.8)	0.265†
Female	11 (29.0)	6 (46.2)
Age (year)**	64.6 (17.5–83.4)	53.8 (24.4–79.1)	0.118†
Race/ethnicity*
Asian	0 (0.0)	1 (7.7)	0.235†
Black	3 (7.9)	0 (0.0)
White	30 (79.0)	10 (76.9)
White/Hispanic	4 (10.4)	1 (7.7)
Other/Hispanic	1 (2.6)	1 (7.7)
Phototype*
Light (FST I–III)	34 (89.5)	12 (92.3)	0.772†
Dark (FST IV–VI)	4 (10.5)	1 (7.7)
Comorbidities*
HTN	13 (34.2)	4 (30.8)	1.000‡
DLP	16 (42.1)	5 (38.5)	0.529‡
GERD	10 (26.3)	5 (38.5)	0.487‡
Diagnosis subtype*
MF	27 (71.0)	–
SS	5 (13.2)	–
Non‐MF/SS CTCL	6 (15.8)	–
Clinical stage*
Early (IA–IIA)	20 (52.6)	–
Mid/Late (IIB–IVB)	18 (47.4)	–
Disease duration (y)**	2.7 (0.15–29.6)	–
mSWAT**	14 (2–159)	–

CTCL, cutaneous T‐cell lymphoma; DLP, dyslipidemia; FST, Fitzpatrick skin phototype; GERD, gastroesophageal reflux; HTN, hypertension; MF, mycosis fungoides; mSWAT, modified Severity‐Weighted Assessment Tool; SS, Sézary syndrome.

*N (%); **Median (range); †Two‐tailed t‐test; ‡Fisher’s exact test.

Sample collection and DNA extraction

Patients were instructed to swab stool from toilet paper immediately after defecation using study‐provided sterile swabs and tubes. , Samples were sent by overnight mail to our facility and promptly stored at −80°C. Genomic DNA was extracted with a Maxwell® RSC Fecal Microbiome DNA Kit (Promega; Madison, WI, USA) on a Maxwell® RSC Instrument, following the manufacturer's protocol.

16S rRNA amplification and sequencing

Genomic DNA was prepared for sequencing using a two‐stage amplicon sequencing workflow, as described previously. Genomic DNA was PCR‐amplified using primers targeting the V4 region of microbial 16S rRNA genes. The primers, 515F modified and 806R modified, contained 5′ linker sequences compatible with Access Array primers for Illumina sequencers (Fluidigm; South San Francisco, CA). PCRs were performed in a total volume of 10 μL using MyTaq™ HS 2X Mix (Bioline), primers at 500 nmol/L concentration, and approximately 1000 copies per reaction of a synthetic double‐stranded DNA template (described below). Thermocycling conditions were 95°C for 5′ (initial denaturation), followed by 28 cycles of 95°C for 30 s, 55°C for 45 s, and 72°C for 30 s. The second‐stage PCR reaction contained 1 μL of PCR product from each reaction and a unique primer pair of Access Array primers. Thermocycling conditions were 95°C for 5′ (initial denaturation), followed by 8 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s. Libraries were pooled and sequenced on an Illumina MiniSeq sequencer (Illumina; San Diego, CA, USA) with 15% phiX spike‐in and paired‐end 2 × 153 base sequencing reads.

A synthetic double‐stranded DNA spike‐in was created as a gBLOCK by Integrated DNA Technologies (IDT; Coralville, IA). The design basis was a 999 base pair (bp) region of the 16S rRNA gene of Rhodanobacter denitrificans strain 2APBS1T (NC_020541). Portions of V1, V2, and V4 variable regions were replaced by eukaryotic mRNA sequences (Apostichopus japonicus glyceraldehyde‐3‐phosphate dehydrogenase mRNA, HQ292612 and Strongylocentrotus intermedius glyceraldehyde‐3‐phosphate dehydrogenase mRNA, KC775387). Primer sites were preserved, and the overall length of the synthetic DNA did not differ from the equivalent R. denitrificans fragment. PCR amplicons generated from this synthetic DNA do not differ in size from bacterial amplicons and can only be identified and removed through postsequencing bioinformatics analysis. The sequence can be accessed via GenBank (OK324963).

Basic processing

Forward and reverse reads were merged using PEAR v.0.9.6. Merged reads were trimmed using cutadapt v1.18 to remove ambiguous nucleotides and primer sequences, and trimmed based on the quality threshold of P = 0.01. Reads lacking the primer sequence and/or sequences less than 225 bp following merging and quality trimming were discarded. Chimeric sequences were identified and removed using the USEARCH algorithm with a comparison to Silva v132 reference sequence. , Amplicon sequence variants (ASVs) were identified using DADA2 v1.18 and annotated taxonomically using the Naïve Bayesian classifier included in DADA2 with the Silva v132 training set. Synthetic spike‐in sequences were removed before proceeding with downstream bioinformatics analyses.

Alpha diversity analyses

Shannon indices were calculated with default parameters (i.e., base = e) in R using the vegan library v2.5‐6. The data were rarefied to a depth of 5000 counts/sample after removal of spike‐in sequences. A generalized linear model assuming Gaussian distribution was utilized for index modeling and significance (ANOVA) was tested using the F‐test. Post‐hoc, pairwise analyses were performed using the Mann–Whitney test. Plots were generated using GraphPad Prism v9.2 (GraphPad; San Diego, CA, USA).

Beta diversity analyses

The normalized data were square‐root transformed and Bray–Curtis indices were calculated without autotransformation in R using the metaMDSdist function in the vegan library v2.5‐6. The dissimilarity indices were modelled and tested for significance with the sample covariates (PERMANOVA). Additional comparisons of the individual covariates were performed using ANOSIM. Plots were generated in R using the ggplot2 library.

Taxonomic differential analysis

Differential analyses of taxa as compared to experimental covariates were performed using edgeR v3.28.1 on raw sequence counts. The data were filtered to remove sequences of chloroplast or mitochondrial origin and taxa accounting for less than 0.1% of the total sequence count. Data were normalized as counts per million and fit using a negative binomial generalized linear model using experimental covariates. Statistical tests were performed using a likelihood ratio test. Adjusted P‐values (q‐values) were calculated using the Benjamini–Hochberg false discovery rate (FDR) correction. Significant taxa were determined based on an FDR threshold of 5%. Plots were generated using GraphPad Prism v9.2.

Data accession

The raw 16S rRNA sequences reported here are accessible on the NCBI Short Read Archive (PRJNA767860).

Results

Clinical characteristics of CTCL patients and healthy controls

Thirty‐eight unique CTCL patients and 13 age‐matched HC were enrolled in this study. All subjects were from the same geographical region (Chicago metropolitan area, United States) to reduce the impact of environmental variation on the microbiota. Four HC‐CTCL pairs sharing a home were selected for even closer matching. To avoid bias in sample collection, manipulation, and analysis, we concurrently enrolled patients and controls rather than rely on publicly available human microbiome data.

Twenty‐seven (71.0%) CTCL patients were male; median patient age was 64.6 (range 17.5–83.4) years. Seven (53.8%) HC were male; the median age was 53.8 (range 24.4–79.1) years. There was no significant difference in gender, age, race/ethnicity, or phototype between the groups (Table 1). Twenty patients had early‐stage disease (stage IA‐IIA; 52.6%), while 18 had mid/late‐stage disease (stage IIB–IVB; 47.4%); stage IB was the most common overall (n = 12, 31.6%). The median disease duration from diagnosis to sample collection was 5.4 (range 0.15–29.6) years and the median mSWAT was 14 (range 2–159). The great majority (73.7%) of patients were not on any systemic CTCL therapies. Six (15.8%) patients had a remote history of non‐CTCL cancer, but all were in remission at the time of sample collection.

The study groups shared similar comorbidity profiles, the most common comorbidities being hypertension (CTCL: n = 13, 34.2% vs. HC: n = 4, 30.8%; Fisher’s exact test P = 1.00), dyslipidemia (CTCL: n = 16, 42.1% vs. HC: n = 5, 38.5%; P = 0.529), and gastroesophageal reflux (CTCL: n = 10, 26.3% vs. HC: n = 5, 38.5%; P = 0.487). There were no differences in surveyed dietary patterns, including processed food intake, dairy intake, organic/hormone‐free meat intake, pre/probiotic use, and alcohol consumption (Table S4).

Sequencing and taxonomic characteristics of sample cohort

Bacterial DNA extracted from stool samples was used for 16S rRNA gene amplicon sequencing. A total of 4 914 148 paired‐end reads were generated with an average of 96 355 reads per sample. Quality filtering of the reads produced a total of 4 194 176 reads with an average of 82 239 reads per sample. Swab, reagent, and PCR controls were negative for any significant contamination. In total, 472 genera, 152 families, 73 orders, 33 classes, and 18 phyla of microorganisms were detected in all samples. The 4 most abundant phyla in both subject groups were Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria – all normally found in the human gut (Fig. 1).

Figure 1

Relative sequence abundance of bacterial taxa in fecal samples at the phylum, class, order, family, and genus levels. Relative sequence abundances (%) were calculated for the cutaneous T‐cell lymphoma (CTCL) patient and healthy control (HC) cohorts at each taxonomic level. Phylum and class were filtered to highlight all taxa with greater than 1.0% relative abundance, order was filtered to greater than 5.0% relative abundance, and family and genus were filtered to greater than 10.0% relative abundance. Each taxonomic level is visualized by the individual subject (left) and the mean relative abundance of each bacterial taxa (right). The mean relative abundances are also delineated for each level (right).

There were no significant differences in α‐diversity at any taxonomic level, but α‐diversity trended lower in CTCL patients vs. controls at the genus level (Figs 2a and S1; Kruskal–Wallis P = 0.17). The β‐diversity analysis demonstrated no global differences in microbial community structure between patients and controls (Fig. 2c; PERMANOVA R 2 = 0.0193, P = 0.488; ANOSIM R = −0.013, P = 0.627).

Figure 2

α‐ and β‐diversity of the gut microbiota of CTCL patient and HC cohorts. (a) α‐diversity trended lower among all CTCL patients compared to HC but was not statistically significant (P = 0.17), as represented by the Shannon diversity score. Dots are colour‐coded for mycosis fungoides (MF)/Sézary syndrome (SS) clinical stage (left) and modified Severity‐Weighted Assessment Tool (mSWAT) (right) divisions. Group medians are denoted by coloured horizontal bars. (b) Among advanced CTCL patients, α‐diversity was significantly lower compared to HC (P = 0.015), as represented by the Shannon diversity score. Dots are colour‐coded for mSWAT divisions and group medians are denoted by coloured horizontal bars. (c) Multidimensional scaling (MDS) plots of gut microbial communities based on Bray–Curtis dissimilarity analysis performed at the taxonomic level of genus shows no global differences in gut microbial community structure between CTCL patient and HC samples (PERMANOVA R 2 = 0.019, P = 0.49) or between advanced CTCL patient and HC samples (R 2 = 0.038, P = 0.15).

Gut microbiome of CTCL patients shows shifts of specific bacterial taxa

Taxon‐by‐taxon analysis revealed that the abundance of certain bacterial populations in CTCL patients is significantly different than that of HC (Fig. 3; Table S2A). When compared to controls, the CTCL group showed significant decreases in the relative abundance of bacteria from the phylum Actinobacteria (q < 0.001), classes Coriobacteriia (q < 0.01) and Actinobacteria (q = 0.03), and order Coriobacteriales (q < 0.01) (Fig. 3a–c). Reduced mean relative abundance of Actinobacteria at the phylum and class levels correlated with increased stage and mSWAT, but this trend was less compelling within the Coriobacteriia and Coriobacteriales data.

Figure 3

Specific gut bacterial taxa differ between CTCL patients and controls. Dot plots illustrate the relative sequence abundance (%) of taxa that were significantly different in CTCL patients organized by MF/SS clinical stage (left) and mSWAT (right) vs. HC at the (a) phylum, (b) class, (c) order, (d) family, and (e) genus levels. Data are shown on a log scale. Group means are denoted by coloured horizontal bars. †Advanced CTCL only.

At the genus level, Anaerotruncus (q = 0.01) was significantly less abundant in the gut microbiota of CTCL patients vs. HC (Fig. 3e) and mean relative abundance was directly related to clinical stage and mSWAT. Unclassified Eggerthellaceae (q < 0.01) and unclassified Enterobacteriaceae (q < 0.01) were also less abundant, but they did not fit into defined clades and were not graphed. Other genera that trended towards significance (P < 0.05 but q > 0.05) included Bifidobacterium, Collinsella, unclassified Clostridiales Family XIII, Romboutsia, Angelakisella, and unclassified Erysipelotrichaceae, which were more abundant in controls vs. patients; and Prevotella, Erysipelatoclostridium, Faecalitalea, unclassified Burkholderiaceae, Solobacterium, Lawsonella, and Dielma, which were more abundant in patients vs. controls.

CTCL patients with greater active disease burden are associated with distinct microbial communities and reduced microbial richness compared to healthy controls

We separately examined patients with greater active disease burden – defined as CTCL stage IB or higher and substantial skin involvement (mSWAT > 10) at the time of sample collection (designated ‘advanced CTCL’) – because they are more likely to have systemic rather than skin‐only immune dysregulation and dysbiosis. In this analysis, α‐diversity at the genus level was significantly lower in patients compared to controls (Fig. 2b; P = 0.015). β‐diversity analysis revealed significant dissimilarity between groups at the phylum (PERMANOVA R 2 = 0.115, P = 0.016), class (R 2 = 0.091, P = 0.014), order (R 2 = 0.072, P = 0.023), and family (R 2 = 0.054, P = 0.049) levels, but not at the genus level (R 2 = 0.038, P = 0.149) (Fig. 2c). Differential taxonomic analysis of this more advanced cohort (Table S2B) identified a significantly decreased relative abundance of Eggerthellaceae (q = 0.01) and Lactobacillaceae (q = 0.01) at the family level, but not at the genus level (Fig. 3d). Reduced mean relative abundance of Eggerthellaceae was associated with lower clinical stage and mSWAT, but the inverse was observed for Lactobacillaceae. Those genera approaching significance in this analysis mirrored those in the full dataset. Additional genera with relatively lower abundance in advanced CTCL vs. controls included Lactobacillus, Lachnospiraceae ND3007 group, and Oxalobacter. Sellimonas and unclassified Christensenellaceae were relatively more abundant in the gut of patients with advanced CTCL compared to controls.

Discussion

Using 16S rRNA gene amplicon sequencing, we elucidated the gut microbial profiles of 38 CTCL patients and 13 healthy, age‐matched individuals. Our results suggest bacterial dysbiosis is more pronounced in CTCL patients with more advanced disease. To the best of our knowledge, this is the first study to characterize the gut microbiome in CTCL and describe the alterations therein. It is also one of the largest microbiome sample sets for this disease to date. The rarity of CTCL underscores the importance of this carefully curated dataset, despite its relatively smaller sample size compared to studies examining substantially more common cancers.

Gut dysbiosis has been explored in myriad diseases, including cutaneous conditions, , , , , , , and hematological , , and visceral malignancies (Table S3). , , These studies have revealed certain dysbiotic signatures are associated with cytokine cascades and skewed immune activity that promote disease progression. AD and psoriasis – inflammatory skin diseases often discussed against CTCL because of their distinct skin and immunological features – are among those conditions affected by this phenomenon (Table 2).

Table 2

Immunologic, skin barrier, and skin microbiome differences between CTCL, atopic dermatitis, and psoriasis

CTCL	Atopic dermatitis	Psoriasis
Immunologic features
Th2‐predominant	Th2‐predominant	Th17/Th1‐predominant
Serum IgE levels correlate with pruritis 70 and baseline eosinophilia is a prognostic factor of poor outcomes 71	Increased IgE levels and circulating eosinophils 44	Normal IgE levels and circulating eosinophils 44
Skin barrier features
Frequent Staphylococcus aureus skin colonization	Frequent Staphylococcus aureus skin colonization	Bacterial colonization uncommon
Skin infections increasingly common with advanced disease	Skin infections common	Skin infections rare
Decreased antimicrobial peptides (S100A7, S100A8, and S100A9) in lesional skin, conferring reduced antimicrobial activity 70	Decreased S100A7 and S100A8 expression, conferring reduced antimicrobial activity 70	Enhanced S100A7 and S100A8 expression, conferring enhanced antimicrobial activity 70
Decreased filaggrin and loricrin expression in patch and plaque CTCL lesions, indicating loss of normal skin barrier function; increased expression in tumor and erythrodermic CTCL 70	Decreased filaggrin and loricrin expression, indicating loss of normal skin barrier function 70	Decreased filaggrin and loricrin expression, indicating loss of normal skin barrier function 70
Skin microbiome features
Skin bacterial shifts may correlate with disease progression; no differences are observed in diversity 72	Decrease in skin microbiome diversity correlates with increased disease severity 73	Reduced skin microbiome diversity compared to healthy individuals 74

CTCL and AD are characterized by Th2‐predominant cytokine profiles, severe pruritis, similar skin barrier defects, and frequent Staphylococcus aureus colonization and infection. , In contrast, psoriasis is Th17/Th1‐driven with opposing skin barrier features compared to CTCL and AD, and patients infrequently develop skin infections or bacterial colonization. , Consistent with these differences, we observed the direction of taxa shifts in CTCL contrasts with that in psoriasis , , , , and is more aligned with that found in AD (Table 3). , , , ,

Table 3

Gut microbiome differences between CTCL, atopic dermatitis, and psoriasis

CTCL	Atopic dermatitis	Psoriasis
Bifidobacteriaceae Bifidobacterium ↓	Bifidobacteriaceae (not seen) 10 Bifidobacterium (not seen) 10 /↓ 75	Bifidobacteriaceae ↑ 13
Coriobacteriaceae Collinsella ↓	Coriobacteriaceae ↓ 10 Collinsella (infants) ↓† 76	Coriobacteriaceae ↑ 13 Collinsella aerofaciens ↑ 11
Eggerthellaceae Unclassified Eggerthellaceae ↓	Eggerthellaceae Eggerthella ↓ 10 and (infants) ↑† 76	Eggerthellaceae ↑ 13
Prevotellaceae Prevotella ↑ Prevotella 6 ↓* Unclassified Prevotellaceae ↑	Prevotellaceae Paraprevotella (infants) ↑† 76 Prevotella (infants) ↓† 76 Prevotella stercorea (not seen) 48	Prevotellaceae ↑ 13 Prevotella ↑ 12 Prevotella copri ↓ 11
Lactobacillaceae Lactobacillus ↓*	Lactobacillaceae Lactobacillus (infants) ↓† 76	Lactobacillaceae ↓ 13
Erysipelotrichaceae Dielma ↑ Faecalitalea ↑ Erysipelatoclostridium ↑ Solobacterium ↑ Unclassified Erysipelotrichaceae ↓	Erysipelotrichaceae (not seen) 10 Bulleidia (not seen) 10 Erysipelotrichaceae incertae sedis (infants) ↓† 76	Erysipelotrichaceae ↑ 13
Clostridiales Family XIII ↓ Unclassified Family XIII ↓	Clostridiales Clostridium cluster IV ↓ 77 Clostridium cluster XI (infants) ↓† 76 Clostridium cluster XIVa (infants) ↑† 76 Clostridium cluster XIVb (infants) ↑† 76	Clostridiales Family XIII ↑ 13 Ruminococcaceae ↑ 13
Ruminococcaceae Anaerotruncus ↓ Angelakisella ↓	Ruminococcaceae ↓ Ruminococcus (infants) ↓† 76	Ruminococcaceae ↑ 13 Ruminococcus ↑ 12 , 33 Ruminococcus gnavus ↑ 11
Lachnospiraceae Lachnospiraceae ND3007 group ↓* Sellimonas ↑*	Lachnospiraceae Anaerostipes (infants) ↑† 76 Blautia ↓ 10 and ↑ 48 Coprococcus ↓ 10 and (infants) ↑† 76 Dorea (infants) ↓† 76 Roseburia (infants) ↑† 76	Lachnospiraceae ↑ 13 Coprococcus ↓ 45 Dorea formicigenerans ↑ 11
Peptostreptococcaceae Romboutsia ↓	‐	Peptostreptococcaceae ↑ 13
Burkholderiales Burkholderiaceae Unclassified Burkholderiaceae ↑	Burkholderiales Sutterellaceae Sutterella ↑ 10 and (infants) ↓† 76	Burkholderiales Burkholderiaceae ↓ 13 and ↑ 78
Oxalobacteraceae Oxalobacter ↓*

↑ enriched in patients; ↓ decreased in patients.

Advanced disease only.

Patients aged 0–12 months old.

Bifidobacteriaceae (not seen)

Bifidobacterium (not seen) /↓

Collinsella (infants) ↓†

Collinsella aerofaciens ↑

Eggerthella ↓ and (infants) ↑†

Paraprevotella (infants) ↑†

Prevotella (infants) ↓†

Prevotella stercorea (not seen)

Prevotella copri ↓

Lactobacillus (infants) ↓†

Erysipelotrichaceae (not seen)

Bulleidia (not seen)

Erysipelotrichaceae incertae sedis (infants) ↓†

Clostridium cluster IV ↓

Clostridium cluster XI (infants) ↓†

Clostridium cluster XIVa (infants) ↑†

Clostridium cluster XIVb (infants) ↑†

Ruminococcus (infants) ↓†

Ruminococcus ↑ ,

Ruminococcus gnavus ↑

Anaerostipes (infants) ↑†

Blautia ↓ and ↑

Coprococcus ↓ and (infants) ↑†

Dorea (infants) ↓†

Roseburia (infants) ↑†

Coprococcus ↓

Dorea formicigenerans ↑

Sutterella ↑ and (infants) ↓†

Burkholderiaceae ↓ and ↑

Considering these comparisons and given the bidirectional influence exercised by the gut microbiome and immune system, the question remains whether the changes in the gut microbiota are epiphenomena of disease or if dysbiosis influences disease course. Notably, Coriobacteriaceae, Lactobacillus, and Bifidobacterium, which appear reduced in CTCL, have been shown to be beneficial commensals. Bacteria from the family Coriobacteriaceae are known to strengthen gut barrier function ; species from the genus Lactobacillus are capable of cytokine‐based anti‐inflammatory activity ; and species from the genus Bifidobacterium may promote daily colonic epithelial renewal, which inhibits the overgrowth of pathogenic species. Additionally, the presence of Ruminococcaceae has been shown to inversely correlate with IL‐6 and C‐reactive protein levels. The dysbiotic signatures characterizing CTCL, AD, and psoriasis may be explained by the complex interactions constituted by the gut microbiome and immune system.

Furthermore, gut microbial signatures of several malignancies share broad themes with the dysbiosis identified here. Certain bacterial subpopulations are known to influence oncogenesis through direct immune modulation and the systemic reach of bacterial metabolites. Prevotellaceae and Bacteroidaceae have been shown to promote local and distant Th17 differentiation and subsequent IL‐17 release. Multiple groups have discussed a causative link between chronic IL‐17‐mediated inflammation and malignancy. , Th17 cells may contribute to CTCL pathogenesis, but this relationship remains to be fully explored. The loss of butyrate‐producing species and enrichment of lipopolysaccharide‐secreting species have also been linked to tumor proliferation. , , , , , Butyrate, a short‐chain fatty acid, can induce cancer cell apoptosis via inhibition of histone deacetylase activity, , , a mechanism mirrored by the chemotherapy agents vorinostat and romidepsin, which are a crucial part of the treatment armamentarium for relapsed and refractory CTCL. Meanwhile, lipopolysaccharide, an endotoxin characteristic of the Proteobacteria phylum, stimulates tumorigenic cytokine cascades. Our CTCL samples were associated with the loss of butyrate producers (e.g., Bifidobacterium and Anaerotruncus in the total cohort and Lactobacillus in advanced CTCL) and Proteobacteria dysbiosis. Further research on the influence of these metabolic pathways in CTCL is warranted.

Study limitations included our small sample size and patient heterogeneity. Though most of our patient subjects were diagnosed with MF and SS, there was variation in disease subtype and staging across our cohort. Despite the wide range of disease duration amongst MF/SS patients, regression analyses demonstrated there was no significant association between disease duration and measures of gut dysbiosis. While patients with recent antibiotic use were excluded, those with advanced disease may have had a distant history of frequent antibiotic exposure with lasting microbiome consequences; however, the dysbiosis identified in these patients remains noteworthy as it still may influence immune function and disease progression. Other treatments were not considered significant confounders because the taxa shifts associated with early‐stage disease matched those of late‐stage disease. We nonetheless chose patients who were naive to or had a long and consistent history of systemic CTCL therapy given its potential to influence the gut microbiota. Among the few systemic treatments cited by patients, only methotrexate has been found to influence the microbiome, but there is no correlation between the differentially abundant taxa in our CTCL group and the taxa inhibited by this drug. Furthermore, although retinoids and immune modulators like interferon‐α may indirectly interact with the immune‐mediated microbiome, more research is needed to support this theory. The similar comorbidity profiles of the study groups helped to control for differences in non‐CTCL medication use. Moreover, our simultaneous enrollment of healthy individuals alongside patients helped control for age, geography, temporality, and data processing within this comparative analysis.

CTCL may be a disease of global dysbiosis. The link connecting CTCL and the microbial world is epitomized by the often‐reported association between disease progression and cutaneous S. aureus infections. Interestingly, Hu and colleagues demonstrated in a mouse model that fecal transplants may not only reverse gut dysbiosis but also reestablish barrier function (e.g., blood–milk barrier) and reduce S. aureus mastitis morbidity. These results further support the importance of the gut microbiota on distant host barrier sites. Indeed, this work may translate into future therapeutic clinical trials for CTCL utilizing probiotics or fecal microbial transplants, which comprise an emerging treatment focus within cancer research. , ,

This is the first study to characterize the gut dysbiosis associated with CTCL, which increases with disease severity and potentially contributes to the severe immune dysfunction concomitant with advanced disease. While multicentre, longitudinal, and treatment‐based studies will add important insights to this discussion, the microbiome and CTCL appear to be intimately connected. Future endeavours within this promising area of research may improve our understanding of CTCL pathophysiology; identify diagnostic and prognostic markers to improve advanced disease management; and possibly elucidate novel therapeutic options for this still poorly understood disease.

Table S1. Detailed demographic characteristics and comorbidities of patients (n = 38) and healthy controls (n = 13).

Table S2. Differential taxonomic analysis demonstrates unique microbial signatures at the genus level in stool samples from (A) CTCL patients vs. healthy controls, and (B) advanced CTCL patients vs. healthy controls.

Table S3. Comparison of relative taxonomic abundance in the gut microbiomes of patients vs. healthy controls as observed in CTCL and as previously reported in other cutaneous diseases and malignancies.

Table S4. Dietary characteristics of patients (n = 38) and health controls (n = 13).

Figure S1. No significant differences in α‐diversity, as represented by Shannon diversity scores, were detected in (a) CTCL patients vs. HC, and (b) advanced CTCL patients vs. HC, except for at the genus level. Group medians are represented by horizontal bars.

Click here for additional data file.