Skip to main content

    Metagenomics quickstart

    Step-by-step workflows for taxonomic classification, community profiling, diversity analysis, and differential abundance testing.

    metagenomicsquickstart

    You want to know which organisms are in your samples and how their communities differ; this walks you through four common workflows — classifying metagenome reads, generating community profiles, calculating diversity metrics, and identifying differentially abundant taxa — from raw FASTQ through statistical comparison.

    Metagenomics quickstart

    This guide covers four common workflows: classifying metagenome reads, generating community profiles, calculating diversity metrics, and identifying differentially abundant taxa.

    Workflow 1: Classify metagenome reads

    Upload and classify

    Ask Ivy in chat (or submit the job directly with POST /api/v1/jobs; see the API reference):

    Classify my gut metagenome sample using Kraken2

    Upload your FASTQ file(s) when prompted, or reference a previously uploaded file. Ivy runs metagenomic classification with Kraken2 as the classifier.

    Review classification results

    FieldExample Value
    Total reads1,250,000
    Classified reads1,087,500
    Classification rate0.87
    Top taxon (taxa[0])Bacteroides fragilis, abundance 0.18

    Switch to MetaPhlAn for species-level accuracy (optional)

    Now run MetaPhlAn on the same sample for species-level profiling

    MetaPhlAn uses clade-specific marker genes and typically reports fewer but more confident species-level assignments.

    Workflow 2: Generate a community profile

    Create an abundance profile

    Generate a genus-level abundance profile from my metagenome samples

    Or configure manually:

    ParameterValue
    Methodmetaphlan (default)
    Resolutiongenus

    Review the community profile

    metagenomics_profile returns one summary per sample, not a multi-sample abundance matrix:

    FieldExample Value
    Richness187 taxa
    Evenness0.74
    Dominant taxonBacteroides
    Dominant abundance0.31
    Top phylum (phylum_abundances)Firmicutes: 0.45, Bacteroidetes: 0.38

    Repeat for each sample

    Run metagenomics_profile once per sample, then pass the resulting profiles into microbiome_diversity or microbiome_differential to compare across samples.

    Workflow 3: Calculate diversity metrics

    Compute alpha diversity

    Calculate Shannon diversity for my microbiome samples

    Or configure manually:

    ParameterValue
    Abundance table URIfrom the profiling step
    Metadata file URICSV with sample, group columns
    Metricsboth (default)

    Review alpha diversity results

    SampleShannonSimpsonChao1Observed species
    Sample_13.420.91201187
    Sample_23.280.89184172
    Sample_32.950.85169156

    Review beta diversity (optional)

    Compare community composition between all samples using Bray-Curtis dissimilarity

    Set beta_method: "bray_curtis" (the default). This returns a Bray-Curtis distance to every other sample under each sample's beta_distances field — there's no ordination (PCoA/NMDS) output.

    Sample_1Sample_2Sample_3
    Sample_10.0000.2130.387
    Sample_20.2130.0000.298
    Sample_30.3870.2980.000

    Workflow 4: Differential abundance analysis

    Define groups and run analysis

    Which taxa are differentially abundant between healthy and diseased groups?

    Or configure manually:

    ParameterValue
    Abundance table URIUploaded or generated table
    Metadata file URICSV with sample, condition columns
    Contrast"diseased_vs_healthy"
    Methodancombc (default)

    Review differentially abundant taxa

    significant_taxa is already filtered to adjusted p-value < 0.05:

    TaxonLog2 Fold ChangeAdjusted P-value
    Fusobacterium nucleatum3.210.001
    Porphyromonas gingivalis2.870.003
    Faecalibacterium prausnitzii-2.450.008
    Roseburia intestinalis-1.980.021

    A positive log2 fold change means the taxon is enriched in the second half of the contrast (here, diseased); negative means depleted.

    Interpret the results

    Log2 Fold ChangeInterpretation
    > 2.0Strong enrichment in group 2 (diseased)
    1.0 to 2.0Moderate enrichment
    -1.0 to 1.0Not biologically significant
    -2.0 to -1.0Moderate depletion
    < -2.0Strong depletion in group 2

    Next steps

    • See inputs for all parameters, database options, and method configurations
    • See outputs for result field descriptions and metric interpretation
    • Use Genomics to look up functional annotations of identified species
    • Use Literature to search for microbiome-disease association studies

    Was this page helpful?

    Loading...