I am doing some metagenomic analysis of samples which can only contain microbes due to the experimental setup. For this, it would be useful to be able to distinguish whether a specific organism is a microbe or not, as it allows:

  • to detect potentially false findings,
  • to avoid an unnecessary inflation of the metagenomic database.

For animals and plants, this is straightforward, as they are never microbes. However, for fungi (and some other clades) it gets complicated. At the very least, I failed to find a simple list of microbial or non-microbial clades.

Ideally, I would like to have a function that maps a taxonomic ID to a boolean stating whether this organism is microbial. I acknowledge that this is utopic since microbe is not that well defined, there are facultatively multicellular organisms, and for some cases it may simply not be known. I am therefore happy with any reasonable approximation such as unicellularity. In particular, given my application, a small chance of falsely identifying an organism as microbial is not an issue.

I am currently working with the NCBI database accessed via the ETE 3 module for Python, but I am open to use other databases and tools.

  • $\begingroup$ Hi, you may eventually try to find marker genes related to intercellular communication? I do not know whether there are proteins specific to algae synapsis or cell junctions in fungi for example. Howver, it may be quite challenging to identify such markers! $\endgroup$ Jun 4 '20 at 9:31
  • $\begingroup$ @thomasdugedebernonville: That multicellularity has evolved roughly thirty times and lost a few times, I doubt that anything universal exists. Cataloguing the tree of life by marking the nodes where multicellularity emerged and ceased certainly beauts that. $\endgroup$
    – Wrzlprmft
    Jun 4 '20 at 11:00
  • $\begingroup$ Yes but convergent evolution probably drives the appearance of conserved functions at cell boundaries. However, proteins with similar roles will probably too divergent for your purposes. By extending your statement " anything universal exists", it is gonna be pretty tricky to get a reasonable approximation for your classification. Hope you will be able too, however! $\endgroup$ Jun 4 '20 at 11:47
  • $\begingroup$ @thomasdugedebernonville: Mind that my point is not that I want something universal, but since there is no reason to expect something universal, I would have to check all cases of multicellularity anyway, in which case I might as well build a database of those myself. $\endgroup$
    – Wrzlprmft
    Jun 4 '20 at 13:26
  • 1
    $\begingroup$ you may have a look at this publication maybe? academic.oup.com/gbe/article/8/4/1279/2574110 $\endgroup$ Jun 4 '20 at 14:19

I would start by saying that your statement that no animals/plants are microbes is false. There are indeed microbial algae (plants) and microbial animals.

So it is indeed difficult. (I'll briefly note that unless your protocol specifically gets rid of environmental DNA, you will probably have some contaminating non-microbial DNA no matter what.)

However, we can use a different version of your proposition to make a simple (over-simple, most would argue) approximation: all prokaryotic organisms, e.g. bacteria and archaea, are microbial, and all eukaryotic organisms are not. This is not hard to do if you are already using NCBI taxonomy, as you can use e.g. this mapping. It has the advantage that all bacteria and archaea are definitely microbial.

Of course it's quite wrong, as you already know. In principle you could set up a complex decision tree sorting this out, though I imagine that's the thing you wish to avoid.

In terms of other genome databases, you might consider looking into GOLD, which has extensive metadata in a slightly less arcane setup than NCBI, which might be helpful. But there is a high rate of missingness there, and also there is still not necessarily a great "microbe"/"not microbe" mapping.

There are also existing databases of microorganisms, which may or may not satisfy your requirements; I have little familiarity with most of those so caveat emptor.

I believe (hope) that NCBI taxid is a good way to do lookups in all such databases, but it may not be true.

  • $\begingroup$ Thank you for your answer, but your simple approximation is indeed too simple as the vast majority of eukaryotes in my samples (and probably others with similar issue) are unicellular microbes, mostly fungi. — In principle you could set up a complex decision tree sorting this out, though I imagine that's the thing you wish to avoid. – Well, if I go by unicellularity, I would “only” need a list of the roughly thirty clades that developed multicellularity and the few cases which lost it. So far, this seems like the most simple solution to me. $\endgroup$
    – Wrzlprmft
    Jun 4 '20 at 9:08
  • $\begingroup$ @Wrzlprmft if you feel confident to implement this decision tree, then I think that might be the best approach. My concern in that case is that basal eukaryotes are relatively quite poorly described and that might affect your results, but I don't have a better alternative. If someone hasn't already implemented a procedure like this, it would honestly be a bit of a public service to the metagenomics community for such an approach to be publicly available. $\endgroup$ Jun 4 '20 at 16:51

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.