Results of electronic database search and study selection
We identified 3,874 potentially relevant publications by searching three electronic databases (Figure 1, Methods). We screened the title and abstract of all publications identified in the initial database search. Primary screening identified 238 studies for full-text retrieval and eligibility checking. A total of 82 case-control studies were included (Supplementary Table 3-4).
Characteristics and quality assessment of the included studies
We included 37 autoimmune and 45 cancer studies, amassing 4,208 healthy human controls and 5,957 disease cases from 27 countries (Table 1, Supplementary Table 5). We captured nine types of autoimmune diseases: type 1 diabetes (n=13) [23–35], arthritis (n=9) [36–44], multiple sclerosis (n=7) [36,45–50], systemic lupus erythematosus (n=4) [51–54], Graves' disease (n=2) [55,56], primary Sjögren's syndrome (n=1) [53], pemphigus vulgaris (n=1) [57], anti‐NMDAR encephalitis (n=1) [58], and autoimmune hepatitis (n=1) [59]. The cancer studies constituted ten types of cancer including colorectal cancer and adenoma (n=27) [60–86], gastric cancer and adenocarcinoma (n=6) [87–92], pancreatic cancer and ductal adenocarcinoma (n=3) [93–95], prostate cancer (n=2) [96,97], cervical cancer (n=2) [98,99], lung cancer (n=1) [100], thyroid papillary cancer (n=1) [101], breast cancer (n=1) [102], hepatocellular carcinoma (n=1) [103], and GIT neoplasia (n=1) [104]. Overall, the majority of the studies (79%) have low risk of bias (based on the nine criteria of the Newcastle-Ottawa Scale) (Figure 2). However, there was a high risk of bias in some indicators such as bias due to lack of representative cases (49%) and improper selection of controls (44%). The lack of representative case bias was mainly due to the absence of consecutive or obviously representative series of cases, and the improper selection of controls bias was due to inclusion of controls from hospital or no description about the selection of study control.
Significant associations of alpha-diversity in cancer and autoimmune studies have opposite trends
The majority of studies (65/82) reported alpha diversity as a measure of microbiome diversity within cases or controls (Figure 3a). The most frequently used alpha-diversity indices were Shannon, Chao and Simpson, reported in 56, 29 and 22 studies, respectively. Overall, most indices showed no significant difference in alpha-diversities between cases and controls in both cancer (62%) and autoimmune (67%) studies. However, an increasing trend of alpha-diversity in cancer studies and an opposite decreasing trend in autoimmune studies was observed. The overall alpha-diversity was significantly increased among cases in 25% of the cancer studies and decreased in 13%. Conversely, the alpha-diversity was significantly decreased among cases in 30% of the autoimmune studies and increased in 3% (Figure 3a).
Cancer and autoimmunity contain inverted and distinct microbiome signatures
A microbiome feature can be simultaneously found to be increased in cancer and decreased in autoimmunity or vice versa. We called this type of association the "Janus-like" effect, to allude to their bidirectionality between cancer and autoimmune studies. We identified features that fall into two categories: i) cancer-increased and ii) autoimmune-increased. The 'cancer-increased' category represents the microbiome features that were found to be increased among the cancer cases but decreased or not significantly associated among the autoimmune cases. On the other hand, the 'autoimmune-increased' category represents the microbiome features that were found to be increased in autoimmune cases but decreased or not significantly associated in the cancer cases.
Genus
We identified 214 distinct genus-level associations with either cancer or autoimmune diseases across 73/82 (89%) studies. Of these, 83 genera were reported in more than one cancer or autoimmune study. We identified 30 genera in the cancer-increased category and 10 genera in the autoimmune-increased category (Figure 3b). In the cancer-increased category, six genera (Enterococcus, Parabacteroides, Odoribacter, Paraprevotella, Desulfovibrio, and Oxalobacter) show true opposite directionality (i.e. increased in cancer but decreased in autoimmune diseases) between cancer and autoimmune diseases, whereas 24 genera in the cancer-increased category were increased in cancer but absent in autoimmunity. Of them, Fusobacterium, Peptostreptococcus and Porphyromonas were the most frequently increased genera among the cancer cases which were found to be associated in 13, 11 and eight cancer studies, respectively. In the autoimmune-increased category, Bacteroides and Parasutterella were inversely associated between cancer and autoimmune cases. The opposite directionality of Bacteroides between cancer and autoimmune diseases is the most prominent, which was found to be increased in 12 autoimmune studies and decreased in six cancer studies. The remaining eight genera in the autoimmune-increased category were always increased in autoimmunity but absent in cancer.
We next explored whether these associations are consistent with the relative abundance of the genus between cases and healthy controls. We extracted the relative abundance data of these significantly associated genera, where available, and conducted a random-effects meta-analysis. We conducted the meta-analysis only on genera for which we found the relative abundance data in at least four studies. By this criteria, six genera had sufficient abundance data available in at least four studies. In the autoimmune-increased genus category, meta-analysis was performed on Bacteroides and Eggerthella. Meta-analysis of Bacteroides involved 10 studies including six autoimmune and four cancer studies. The genus Bacteroides was increased among the autoimmune cases by an overall +0.24 log fold change (95% CI, 0.12 to 0.36) compared with healthy controls, whereas it was decreased among the cancer cases by an overall −0.32 log fold change (95% CI, −0.55 to −0.09) (Figure 3c). Between-study heterogeneity among the autoimmune and cancer studies were 48% and 66%, respectively. Four autoimmune studies were included in the meta-analysis of Eggerthella, which showed an increased abundance in autoimmune cases by an overall +1.71 log fold change (95% CI, 1.28 to 2.14) with a between-study heterogeneity of 0% (Figure 3d). In the cancer-increased genus category, meta-analysis was performed on Enterococcus, Parabacteroides, Peptostreptococcus and Fusobacterium. A total of eight studies were included in the meta-analysis of Enterococcus, which showed an overall increase of +1.89 log fold (95% CI, 0.72 to 3.05) among the cancer cases and decrease by −0.70 log fold (95% CI, −1.97 to 0.58) (Figure 3e). The between-study heterogeneity was 76% and 93% among the cancer and autoimmune studies, respectively. We included six studies in the meta-analysis of Parabacteroides (Figure 3f). It accounted for an overall +0.93 log fold increase (95% CI, −1.35 to 3.21) in cancer cases and −0.70 log fold decrease (95% CI, −1.02 to −0.38) in autoimmune cases with 82% and 19% between-study heterogeneity, respectively. Meta-analysis of Peptostreptococcus includes five cancer studies which shows an overall +3.75 log fold increase (95% CI, 0.90 to 6.60) among the cancer cases where the between-study heterogeneity was 85% (Figure 3g). Five cancer studies were included in the meta-analysis of Fusobacterium (Figure 3h). The results showed an overall +0.54 log fold increase (95% CI, 0.17 to 0.92) of Fusobacterium among the cancer cases with a between study-heterogeneity of 64%.
Species
Twenty-eight studies (34%) reported significant species-level taxonomic associations in cancer and autoimmune diseases. A total of 159 species were significantly associated in at least one study, of which 28 species were associated in more than one cancer or autoimmune study. Finally, 20 species were identified in the cancer-increased category and seven species in the autoimmune-increased category (Figure 4a). Fusobacterium nucleatum was the most common species in the cancer-increased category, appearing in eight cancer studies. However, seven of these eight cancer studies were on colorectal cancer and one was on breast cancer. This might actually represent a skew, because colorectal cancer was the most represented cancer among the cancer studies (27/45). Therefore, it is reasonable that F. nucleatum would be the most common cancer-increased species and it may not indicate that F. nucleatum was increased in cancer generally. None of the cancer-increased species was found to be significantly associated in the autoimmune studies. Among the autoimmune-increased species, Bifidobacterium longum and Streptococcus salivarius were associated inversely between cancer and autoimmune diseases. Interestingly, four members of the Eubacterium genus, E. eligens, E. hallii, E. rectale and E. ventriosum were decreased in cancer cases in multiple studies. However, none of these Eubacterium species were significantly associated with autoimmune diseases.
Metabolic pathways
Associations between predicted microbial metabolic pathways and diseases were only reported in 20 studies (24%). A total of 405 predicted metabolic pathways were shown to be significantly associated in at least one cancer or autoimmune study. Most of the pathways were found in only one study, whereas 48 were associated with cancer and autoimmunity in multiple studies. None of the pathways exhibited an inverse association between cancer and autoimmunity, as observed in several genera and species. However, some pathways were increased or decreased in multiple cancer and autoimmune studies such as lipopolysaccharide biosynthesis, which was found to be increased in three autoimmune studies (Figure 4b).
Disease-associated genera and species are taxonomically and functionally diverse
Taxonomic diversity
To understand the phylogenetic and taxonomic relatedness among the associated genera and species, we used the NCBI common tree taxonomy browser tool for the construction of taxonomic trees. All 40 cancer-increased and autoimmune-increased genera belonged to six bacterial phyla. The majority of the genera (n=18) belonged to the phylum Firmicutes followed by Proteobacteria (n=10) (Figure 5a). All 27 species were represented by the same phylum except Verrucomicrobia (Figure 5b). The Firmicutes phylum represents the majority of the species (n=13) followed by Bacteroidetes (n=6). We observed some interesting trends across phyla. The members (both genus and species) of Fusobacteria were always found to be increased in cancers. All the genera and species in Actinobacteria showed a positive association both in cancer and autoimmunity. Interestingly, most of the genera and species of the phylum Bacteroidetes were found to be increasingly associated with cancer and decreasingly associated with autoimmunity. The only exception is the Bacteroides genus which showed an opposite directionality compared with other members of this phylum (Figure 5b). Similarly, the majority of the members of the Proteobacteria showed a positive association with cancer or a negative association with autoimmunity, with the exception of the genera Parasutterella and Sutterella, which showed an opposite directionality. Most of the genera of the Firmicutes phylum were found to be increased both in cancer and autoimmunity. However, the species of the Firmicutes phylum could be divided into two groups considering opposite directionality with cancer and autoimmunity. One group of species were found to be increased in cancer (n=8 species) and another group was found to be decreased in cancer (n=5 species) (Figure 5b).
Functional diversity
To predict the shared or unique metabolic functions among the cancer-increased or autoimmune-increased genera and species, we identified metabolic pathways using MetAboliC pAthways DAtabase for Microbial taxonomic groups (MACADAM) [105], which employs pathway tools based on the MetaCyc database [106] that includes metabolic pathways as well as associated metabolites, reactions and enzymes. We first extracted all the metabolic pathways that are present among cancer-increased and autoimmune-increased genera and species. We then identified the pathways that were exclusively present either in cancer-increased or autoimmune-increased groups. Of the 936 metabolic pathways that were detected together in two groups, 303 pathways were found exclusively in cancer-increased taxa. On the other hand, only 17 pathways were present exclusively in autoimmune-increased taxa (Figure 5c, Supplementary file 1). Among the 303 cancer-increased pathways, 119 had a complete pathway score of 1, meaning that some members of the cancer-increased taxa contain all enzymes to complete the metabolic pathway. The most frequent pathways in the cancer-increased groups belonged to aromatic compound degradation (38/303) followed by non-carbon nutrient metabolism (23/303), amino-acid degradation (21/303) and cofactor biosynthesis (16/303). Interestingly, 17 of the 23 non-carbon nutrient metabolism pathways were related to sulfur-metabolism. Of the autoimmune-increased pathways (n=17), four complete metabolic pathways were contained in some members of the autoimmune-increased taxa. (Figure 5d, Supplementary file 1).
Because we observed taxonomically-opposing directionality within the members of Bacteroidetes, Proteobacteria and Firmicutes (Figure 5a-b), we further explored the functional similarities or dissimilarities within the members of these phyla. In the Bacteroidetes phylum, four genera and six species were associated in the cancer-increased group (Figure 5a-b), whereas only the Bacteroides genus belonged to the autoimmune-increased group. By comparing the potential metabolic functions between the cancer-increased taxa and the Bacteroides genus (autoimmune-increased) under the Bacteroidetesphylum, we detected 53 pathways exclusively present among the cancer-increased taxa and 39 pathways that were exclusively present in the genus Bacteroides (Supplementary file 1). In the Firmicutes phylum, we compared the metabolic functions between cancer-increased (i.e. Anaerococcus, Christensenella, Dialister pneumosintes, Enterococcus, Filifactor, Gemella morbillorum, Lactobacillus, Lactococcus, Megasphaera, Mogibacterium, Parvimonas, Parvimonas micra, Selenomonas, Selenomonas sputigena, and Turicibacter) and autoimmune-increased (i.e. Acidaminococcus intestini, Blautia, Eubacterium eligens, Eubacterium hallii, Eubacterium rectale, Flavonifractor, Gemella, Lachnoclostridium, and Streptococcus salivarius) species. A total of 162 pathways were exclusively present in the cancer-increased group, and 21 pathways were found exclusively in the autoimmune-increased group under the phylum Firmicutes (Supplementary file 1). Similarly, in the Proteobacteria phylum, we compared the metabolic functions in the cancer-increased group (i.e. Aquabacterium, Bilophila wadsworthia, Delftia, Desulfovibrio, Escherichia coli, Klebsiella, Morganella, Oxalobacter and Proteus) and autoimmune-increased group (i.e. Sutterella). A total of 700 pathways were exclusively present in the cancer-increased group and no pathways were found exclusively in the autoimmune-increased taxa in this phylum (Supplementary file 1).