why use gap filling with carveme about metagem HOT 2 CLOSED

Hi Francisco,

Thanks a lot for the explanation, that sounds very reasonable.

And thanks for the pipeline again!

from metagem.

Hi vrmarcelino,

Good question! Indeed, we cannot be sure that the microbes will grow in isolation, and the gapfilling is meant to ensure that the microbes can grow within the community. By default, the gapfilling is done with complete media (M8 in the config file/media_db.tsv, see figure 1d here), which assumes that the microbes should be able to grow in a community/environment where all nutrients are provided. In the metaGEM paper we gapfill on M8 and then simulate on M11 to observe exchanges associated with aromatic amino acids. Also note that the gapfill flag will use reaction confidence scores to prioritize the introduction of reactions/pathways with the most genetic evidence present in the contigs, i.e. not just adding random reactions until growth is achieved. Also, the log files of the carveme jobs will say how many metabolites/reactions are added from gapfilling, and from my experience it is usually none or in the low single digits.

Thanks for your interest in metaGEM and let me know if you have further questions!

from metagem.