Closed nicola-debernardini closed 3 years ago
Hi Nicola, this looks odd, indeed. I ran the whole workflow but can't reproduce the issue. The exact workflow I used:
model="LT608329.1"
gapseq find -p all -m Archaea -t Archaea $model.fna.gz
gapseq find-transport $model.fna.gz
gapseq draft -r $model-all-Reactions.tbl \
-t $model-Transporter.tbl \
-p $model-all-Pathways.tbl \
-c $model.fna.gz \
-u 200 -l 100
# gapfill with minimal media
gapseq fill -m $model-draft.RDS \
-n MM_anaerobic_CO2_H2.csv \
-c $model-rxnWeights.RDS \
-g $model-rxnXgenes.RDS \
-b 100 -e highH2
# gapfill with minimal media plus ALLmed
gapseq fill -m $model-draft.RDS \
-n MM_anaerobic_CO2_H2_plus_ALLmed.csv \
-c $model-rxnWeights.RDS \
-g $model-rxnXgenes.RDS \
-b 100 -e highH2
With the minimal medium I get: Final growth rate: 0.1019939
With the extended medium: Final growth rate: 0.09577336
I am wondering if this could be due to different LP-solvers. In my case I used cplex version 1210. What LP-solver are you using?
Best Silvio
Hi Silvio,
Many thanks for your answer. I realized I used the version 1.1 6ae3a6a to run only the last step ("fill") and version 1.1 d6522aa to run "find", "find-transport" and "draft". You probably solved the problem in the latter versions of one of these three commands.
I tried to run your procedure using version 1.1 d6522aa in all the steps and the results I obtained are identical to yours. The solver I used is cplex version 12.8 but the different version of the solver didn't cause any change.
Also, thanks for sharing you procedure. I never used the "-m" option during the "find" step but I saw it runs much faster when I used it.
Thanks again for the answer!
Cheers, Nicola
Hi, I'm performing the gap filling of Methanothermobacter wolfeii - LT608329.1 using the minimal medium "MM_anaerobic_CO2_H2" as shown in the issue 54 but I obtain different results when try to add compounds to that medium.
The medium used is: MM_anaerobic_CO2_H2.csv
The medium used has this composition: MM_anaerobic_CO2_H2_plus_ALLmed.csv
Do you have any clue why does this happen? Why does the same reference genome gapfilled on medium containing the same compounds and with the same uptake bounds lead to difference growth rates when I add compounds to one of them (in this case not growing at all)?
The models were gapfilled with version (v.1.1 6ae3a6a) - models.zip.
I appreciate any help you can provide.
Best regards, Nicola