Single meal: Baobab, Honey, and Antelope¶

Here we will build a medium based on the anticipated metabolites present in the human gut for a meal matched to the Hadza people in Tanzania. Please be aware that this is unlikely to represent the full diversity of the diet of the Hadza.

The composition of the meal is the following:

200g Baobab
50g Baobab seeds (mixed with water)
100g of Honey
100g of Antelope meat

This sums up to 1051 kcal.

We will start by reading individual tables for the specific foods and scale the abundances.

In [86]:

import pandas as pd

meal = {
    "Baobab pulp": 200,
    "Baobab seed": 50,
    "Honey": 100,
    "Antelope": 100
}

foods = []
for food in meal:
    mets = pd.read_excel("../data/foods_diets.xlsx", food)
    mets["amount_g"] = mets.relative_abundance / mets.relative_abundance.sum() * meal[food]
    mets["concentration"] = mets["amount_g"] / mets["mw"] * 1000.0  # to yield mmol/meal
    mets["flux"] = mets["concentration"] / 8.0
    mets["food"] = food
    foods.append(mets)
foods = pd.concat(foods)

Now we combine the data.

In [87]:

foods.loc[foods.metabolite == "h2o", "flux"] += 1000.0 / 18.01 * 1000.0 / 8.0 # add 1L water
diet = foods.dropna(subset=["metabolite"]).groupby("metabolite").flux.sum().reset_index()
diet.head()

Out[87]:

	metabolite	flux
0	4abut	0.724400
1	Lcystin	0.002091
2	ab14lac_L	0.258269
3	acgal	1.051659
4	acon_C	1.122785

Adjust for intestinal adsorption¶

To achieve this we will load the Recon3 human model. AGORA and Recon IDs are very similar so we should be able to match them. We just have to adjust the Recon3 ones a bit. We start by identifying all available exchanges in Recon3 and adjusting the IDs.

In [88]:

from cobra.io import read_sbml_model
import pandas as pd

recon3 = read_sbml_model("../data/Recon3D.xml.gz")
exchanges = pd.Series([r.id for r in recon3.exchanges])
exchanges = exchanges.str.replace("__", "_").str.replace("_e$|EX_", "", regex=True)
exchanges.head()

Out[88]:

0     5adtststerone
1    5adtststerones
2             5fthf
3             5htrp
4             5mthf
dtype: object

In [89]:

diet["dilution"] = 1.0
diet.loc[diet.metabolite.isin(exchanges), "dilution"] = 0.2
diet["flux"] = diet["flux"] * diet["dilution"] 
diet[["metabolite", "dilution"]].drop_duplicates().dilution.value_counts()

Out[89]:

0.2    52
1.0    19
Name: dilution, dtype: int64

Adding host supplied components¶

Finally we add the host metabolites such as primary bile acids and mucins and a minuscule amount of oxygen.

In [90]:

diet.set_index("metabolite", inplace=True)
annotations = pd.read_csv("../data/agora_metabolites.csv")

# mucin
for met in annotations.loc[annotations.metabolite.str.contains("core"), "metabolite"]:
    diet.loc[met, "flux"] = 1

# primary BAs
for met in ["gchola", "tchola"]:
    diet.loc[met, "flux"] = 1

# anaerobic
diet.loc["o2", "flux"] = 0.001

diet.reset_index(inplace=True)
diet["reaction"] = "EX_" + diet.metabolite + "(e)"
diet

Out[90]:

	metabolite	flux	dilution	reaction
0	4abut	0.144880	0.2	EX_4abut(e)
1	Lcystin	0.002091	1.0	EX_Lcystin(e)
2	ab14lac_L	0.258269	1.0	EX_ab14lac_L(e)
3	acgal	0.210332	0.2	EX_acgal(e)
4	acon_C	1.122785	1.0	EX_acon_C(e)
...	...	...	...	...
82	gncore2_rl	1.000000	NaN	EX_gncore2_rl(e)
83	core7	1.000000	NaN	EX_core7(e)
84	gchola	1.000000	NaN	EX_gchola(e)
85	tchola	1.000000	NaN	EX_tchola(e)
86	o2	0.001000	NaN	EX_o2(e)

87 rows × 4 columns

And we will merge this table with some annotations to make it more accessible.

In [91]:

skeleton = pd.merge(diet, annotations, on="metabolite")

skeleton["global_id"] = skeleton.reaction
skeleton["reaction"] = "EX_" + skeleton.metabolite + "_m"
skeleton.head()

Out[91]:

	metabolite	flux	dilution	reaction	name	hmdb	kegg.compound	pubchem.compound	inchi	chebi	global_id
0	4abut	0.144880	0.2	EX_4abut_m	4-Aminobutanoate	HMDB00112	C00334	119.0	InChI=1S/C4H9NO2/c5-3-1-2-4(6)7/h1-3,5H2,(H,6,7)	NaN	EX_4abut(e)
1	Lcystin	0.002091	1.0	EX_Lcystin_m	L-cystine	HMDB00192	C00491	NaN	InChI=1S/C6H12N2O4S2/c7-3(5(9)10)1-13-14-2-4(8...	NaN	EX_Lcystin(e)
2	ab14lac_L	0.258269	1.0	EX_ab14lac_L_m	L-Arabinono-1,4-lactone	NaN	NaN	NaN	NaN	NaN	EX_ab14lac_L(e)
3	acgal	0.210332	0.2	EX_acgal_m	N-acetyl-D-galactosamine	HMDB00212	C01132	35717.0	NaN	NaN	EX_acgal(e)
4	acon_C	1.122785	1.0	EX_acon_C_m	cis-Aconitate	NaN	NaN	NaN	NaN	NaN	EX_acon_C(e)

Complete the medium¶

Great we now have a pretty good skeleton. One issue that this will never be fully complete. There will always be some components missing that are essential for microbial growth. Fortunately, we provide a algorithm in MICOM to complete a medium with the smallest set of additional components to provide growth to all intestinal taxa.

In [92]:

from micom.workflows.db_media import complete_db_medium

manifest, imports = complete_db_medium(
    "../data/agora103_strain.qza", 
    medium=skeleton, 
    growth=0.1, 
    threads=12, 
    max_added_import=10,  # do not add more than 10 mmol/h of flux per component
    strict=["EX_o2(e)"],  # force anaerobic environment
    weights="mass"        # minimize added molecular weight
)

Completing 818 strain-level models on a medium with 87 components (1 strict).

Output()

Obtained growth for 815 models adding additional flux of 0.82/2940.27 on average.

In [93]:

manifest.can_grow.value_counts()

Out[93]:

True     815
False      3
Name: can_grow, dtype: int64

In [94]:

filled = imports.max()
added = filled.sum() - skeleton.loc[skeleton.reaction.isin(filled.index), "flux"].sum()

print(f"Added flux is {added.sum():.2f}/{filled.sum():.2f} mmol/h.")

Added flux is 153.99/3143.57 mmol/h.

Let's see what was added in large amounts.

In [95]:

added_fluxes = filled.copy()
shared = added_fluxes.index[added_fluxes.index.isin(skeleton.reaction)]
added_fluxes[shared] -= skeleton.flux[shared]
added_fluxes.sort_values(ascending=False)[0:20]

Out[95]:

EX_h_m         10.000000
EX_h2_m        10.000000
EX_glyleu_m    10.000000
EX_no_m        10.000000
EX_ser_L_m     10.000000
EX_glyc_m       8.883167
EX_succ_m       8.277121
EX_co2_m        7.558247
EX_asn_L_m      6.665649
EX_for_m        5.926777
EX_urea_m       4.965212
EX_mal_L_m      4.867167
EX_acald_m      4.838348
EX_cytd_m       4.458517
EX_ph2s_m       4.322594
EX_acac_m       4.121076
EX_fum_m        3.904283
EX_nh4_m        3.869570
EX_no3_m        3.851188
EX_arg_L_m      3.519534
dtype: float64

Looks okay. So we will now assemble the final medium. For this we add the new components to each sample and rebuild the annotations for a nicely formatted medium.

In [96]:

added_df = filled[filled > 1e-8].reset_index() 
added_df.iloc[:, 0] = added_df.iloc[:, 0].str.replace("EX_|_m$", "", regex=True)
added_df.columns = ["metabolite", "flux"]

completed = pd.merge(added_df, annotations, on="metabolite", how="left")
completed["reaction"] = "EX_" + completed.metabolite + "_m"
completed["global_id"] = "EX_" + completed.metabolite + "(e)"
completed

Out[96]:

	metabolite	flux	name	hmdb	kegg.compound	pubchem.compound	inchi	chebi	reaction	global_id
0	ala_L	0.427871	L-alanine	HMDB00161	C00041	5950.0	InChI=1S/C3H7NO2/c1-2(4)3(5)6/h2H,4H2,1H3,(H,5...	NaN	EX_ala_L_m	EX_ala_L(e)
1	arg_L	3.769624	L-argininium(1+)	HMDB00517	C00062	6322.0	InChI=1S/C6H14N4O2/c7-4(5(11)12)2-1-3-10-6(8)9...	NaN	EX_arg_L_m	EX_arg_L(e)
2	asn_L	6.665649	L-asparagine	HMDB00168	C00152	6267.0	InChI=1S/C4H8N2O3/c5-2(4(8)9)1-3(6)7/h2H,1,5H2...	NaN	EX_asn_L_m	EX_asn_L(e)
3	asp_L	3.163472	L-aspartate(1-)	HMDB00191	C00049	5960.0	InChI=1S/C4H7NO4/c5-2(4(8)9)1-3(6)7/h2H,1,5H2,...	NaN	EX_asp_L_m	EX_asp_L(e)
4	ca2	2.710538	calcium(2+)	HMDB00464	C00076	NaN	InChI=1S/Ca/q+2	NaN	EX_ca2_m	EX_ca2(e)
...	...	...	...	...	...	...	...	...	...	...
173	MGlcn54	0.009248	mucin-type O-glycan No 54	NaN	NaN	NaN	NaN	NaN	EX_MGlcn54_m	EX_MGlcn54(e)
174	n2	0.028071	Nitrogen	NaN	NaN	NaN	NaN	NaN	EX_n2_m	EX_n2(e)
175	cmp	0.008410	CMP	HMDB00095	C00055	6131.0	NaN	NaN	EX_cmp_m	EX_cmp(e)
176	datp	0.002391	dATP	HMDB01532	C00131	15993.0	NaN	NaN	EX_datp_m	EX_datp(e)
177	so3	1.072348	Sulfite	HMDB00240	C00094	1100.0	InChI=1S/H2O3S/c1-4(2)3/h(H2,1,2,3)/p-2	NaN	EX_so3_m	EX_so3(e)

178 rows × 10 columns

Validate the medium¶

And we will now validate whether the medium works.

In [97]:

from micom.workflows.db_media import check_db_medium

check = check_db_medium("../data/agora103_strain.qza", medium=completed, threads=12)

Checking 818 strain-level models on a medium with 178 components.

Output()

In [98]:

check.growth_rate.describe()

Out[98]:

count    818.000000
mean       0.190551
std        0.105812
min        0.036577
25%        0.100000
50%        0.156764
75%        0.256409
max        0.646666
Name: growth_rate, dtype: float64

In [99]:

import qiime2 as q2

arti = q2.Artifact.import_data("MicomMedium[Global]", completed)
arti.save("../media/baobab_honey_antelope.qza")

Out[99]:

'../media/baobab_honey_antelope.qza'

In [ ]: