Molecular genetic methods can distinguish divergent evolutionary lineages in what previously appeared to be single species, but it is not always clear what functional differences exist between such cryptic species. We used a metabolomic approach to profile biochemical phenotype (metabotype) differences between two putative cryptic species of the earthworm Lumbricus rubellus. There were no straightforward metabolite biomarkers of lineage, i.e. no metabolites that were always at higher concentration in one lineage. Multivariate methods, however, identified a small number of metabolites that together helped distinguish the lineages, including uncommon metabolites such as Nε-trimethyllysine, which is not usually found at high concentrations. This approach could be useful for characterizing functional trait differences, especially as it is applicable to essentially any species group, irrespective of its genome sequencing status.