Despite recent progress in our understanding of the association between the gut microbiome and colorectal cancer (CRC), multi-kingdom gut microbiome dysbiosis in CRC across cohorts is unexplored. We investigated four-kingdom microbiota alterations using CRC metagenomic datasets of 1,368 samples from 8 distinct geographical cohorts. Integrated analysis identified 20 archaeal, 27 bacterial, 20 fungal and 21 viral species for each single-kingdom diagnostic model. However, our data revealed superior diagnostic accuracy for models constructed with multi-kingdom markers, in particular the addition of fungal species. Specifically, 16 multi-kingdom markers including 11 bacterial, 4 fungal and 1 archaeal feature, achieved good performance in diagnosing patients with CRC (area under the receiver operating characteristic curve (AUROC) = 0.83) and maintained accuracy across 3 independent cohorts. Coabundance analysis of the ecological network revealed associations between bacterial and fungal species, such as Talaromyces islandicus and Clostridium saccharobutylicum. Using metagenome shotgun sequencing data, the predictive power of the microbial functional potential was explored and elevated D-amino acid metabolism and butanoate metabolism were observed in CRC. Interestingly, the diagnostic model based on functional EggNOG genes achieved high accuracy (AUROC = 0.86). Collectively, our findings uncovered CRC-associated microbiota common across cohorts and demonstrate the applicability of multi-kingdom and functional markers as CRC diagnostic tools and, potentially, as therapeutic targets for the treatment of CRC.

Microbial signatures show remarkable potentials in predicting colorectal cancer (CRC). This study aimed to evaluate the diagnostic powers of multimodal microbial signatures, multi-kingdom species, genes, and single-nucleotide variants (SNVs) for detecting precancerous adenomas. We performed cross-cohort analyses on whole metagenome sequencing data of 750 samples via xMarkerFinder to identify adenoma-associated microbial multimodal signatures. Our data revealed that fungal species outperformed species from other kingdoms with an area under the ROC curve (AUC) of 0.71 in distinguishing adenomas from controls. The microbial SNVs, including dark SNVs with synonymous mutations, displayed the strongest diagnostic capability with an AUC value of 0.89, sensitivity of 0.79, specificity of 0.85, and Matthews correlation coefficient (MCC) of 0.74. SNV biomarkers also exhibited outstanding performances in three independent validation cohorts (AUCs = 0.83, 0.82, 0.76; sensitivity = 1.0, 0.72, 0.93; specificity = 0.67, 0.81, 0.67, MCCs = 0.69, 0.83, 0.72) with high disease specificity for adenoma. In further support of the above results, functional analyses revealed more frequent inter-kingdom associations between bacteria and fungi, and abnormalities in quorum sensing, purine and butanoate metabolism in adenoma, which were validated in a newly recruited cohort via qRT-PCR. Therefore, these data extend our understanding of adenoma-associated multimodal alterations in the gut microbiome and provide a rationale of microbial SNVs for the early detection of CRC.

Excess body weight (EBW) increases the risk of colorectal cancer (CRC) and is linked to lower colonoscopy compliance. Here, we extensively analyzed 981 metagenome samples from multiple cohorts to pinpoint the specific microbial signatures and their potential capability distinguishing EBW patients with CRC. The gut microbiome displayed considerable variations between EBW and lean CRC. We identify 44 and 37 distinct multi-kingdom microbial species differentiating CRC and controls in EBW and lean populations, respectively. Unique bacterial-fungal associations are also observed between EBW-CRC and lean-CRC. Our analysis revealed specific microbial functions in EBW-CRC, including D-Arginine and D-ornithine metabolism, and lipopolysaccharide biosynthesis. The best-performing classifier for EBW-CRC, comprising 12 bacterial and three fungal species, achieved an AUROC of 0.90, which was robustly validated across three independent cohorts (AUROC = 0.96, 0.94, and 0.80). Pathogenic microbial species, Anaerobutyricum hallii, Clostridioides difficile and Fusobacterium nucleatum, are EBW-CRC specific signatures. This work unearths the specific multi-kingdom microbial signatures for EBW-CRC and lean CRC, which may contribute to precision diagnosis and treatment of CRC.

Microbial metabolism of bile acids (BAs) is crucial for maintaining homeostasis in vertebrate hosts and environments. Although certain organisms involved in bile acid metabolism have been identified, a global, comprehensive elucidation of the microbes, metabolic enzymes, and bile acid remains incomplete. To bridge this gap, we employed hidden Markov models to systematically search in a large-scale and high-quality search library comprising 28,813 RefSeq multi-kingdom microbial complete genomes, enabling us to construct a secondary bile acid production gene catalog. This catalog greatly expanded the distribution of secondary bile acid production genes across 11 phyla, encompassing bacteria, archaea, and fungi, and extended to 14 habitats spanning hosts and environmental contexts. Furthermore, we highlighted the associations between secondary bile acids (SBAs) and gastrointestinal and hepatic disorders, including inflammatory bowel disease (IBD), colorectal cancer (CRC), and nonalcoholic fatty liver disease (NAFLD), further elucidating disease-specific alterations in secondary bile acid production genes. Additionally, we proposed the pig as a particularly suitable animal model for investigating secondary bile acid production in humans, given its closely aligned secondary bile acid production gene composition. This gene catalog provides a comprehensive and reliable foundation for future studies on microbial bile acid metabolism, offering new insights into the microbial contributions to health and disease. IMPORTANCE: Bile acid metabolism is an important function in both host and environmental microorganisms. The existing functional annotations from single source pose limitations on cross-habitat analysis. Our construction of a systematic secondary bile acid production gene catalog encompassing numerous high-quality reference sequences propelled research on bile acid metabolism in the global microbiome, holding significance for the concept of One Health. We further highlighted the potential of the microbiota-secondary bile acid axis as a target for the treatment of hepatic and intestinal diseases, as well as the varying feasibility of using animal models for studying human bile acid metabolism. This gene catalog offers a solid groundwork for investigating microbial bile acid metabolism across different compartments, including humans, animals, plants, and environments, shedding light on the contributions of microorganisms to One Health.