**Deliverable DWP5-6 (Task 5.9)**
Number	DWP5-6 (Task 5.9)
Due	M46 (Mar 2014)
Description	Database interface for risk-benefit assessment
Nature	P (Prototype)
Dissemination level	PP (Restricted to project partners and the Commission)
Responsible person	THL

Main message:

What did Plantlibra project produce as a database interface for risk-benefit assessment of plant food supplements?

Answer:

The present deliverable (DWP5-6) includes the results of Task 5.9 which is aimed at developing an internet based interface where users can access, combine, by plant and compound, and discuss information on risk and benefit assessment. A detailed description of the structure and functioning of the interface is provided in the present document. The Opasnet wiki-based platform (en.opasnet.org) was used to develop the internet based interface, named “Compound Intake Estimator”. It is a tool that allows the use to calculate the intakes of compounds, based on food or food supplement intake. In addition, the tool allows estimating the health risk of a compound at a defined level of intake. The risk estimated is obtained by combining the information on compound concentration in the plant food supplement with the specific health guidance value available for the compound. As health guidance values, the Acceptable Daily Intake or Tolerable Daily Intake are used when available. Otherwise, the Toxicological Threshold of Concern approach is applied. In case of compounds having a structural alert for genotoxicity, the Margin of Exposure concept is instead used. The interface is structured in such a way that the user can select the compound(s) of interest, the extraction technique, and specify the level of intake. In addition, the user can also add his own composition data. This feature is useful in case a user wants to compare the composition of new botanical preparation with data available from literature and evaluate if any of the compounds present might pose, or not, a health concern at a specific level of intake.

Background and objective of the work

The use of plants and plant derivatives to maintain health has been popular throughout Europe for many centuries. The consumption of teas, digestive drinks, juices, elixirs and extracts prepared from botanicals and used for health maintenance purposes has become part of European cultural heritage. Plant food supplements are a modern-day extension of this process.

PlantLIBRA project (PLANT food supplements: Levels of Intake, Benefit and Risk Assessment) aims to foster the safe use of food supplements containing plants or botanical preparations, by increasing science-based decision-making by regulators and food chain operators. To make informed decisions, competent authorities and food businesses need more quality-assured and accessible information and better tools (e.g. metadata banks). ^[1]

PlantLIBRA is structured to develop, validate and disseminate data and methodologies for risk and benefit assessment and implement sustainable international cooperation. International cooperation, on-spot and in-language capacity building are necessary to ensure the quality of the plants imported in the EU.

The Description of Work for PlantLIBRA project describes the objective of this deliverable in the following way:

Task 5.9: Integration of databases and guidance

Taking into account the database of WP6, THL will develop an Internet-based interface (DWP5-6) where users can access, combine, by plant and compound, and discuss information for risk-benefit assessment. In this context, THL will make selected case-studies available to the public through the Mediawiki platform, when possible. THL will develop guidance for conducting risk-benefit assessments of plants, extracts, preparations and individual products which will be made available online as Deliverable 5-6 (DWP5-6).

Compound intake estimator

Compound intake estimator calculates intakes of compounds based on food or food supplement intake, compound concentration in food, and guidance values for the compound. In the text, we talk about plant-based food supplements (PFSs), but the tool is generic and is usable for other products as well as long as the required concentration data is made available.

Question

How to estimate intakes of chemical compounds in plant-based food supplements and compare them to guidance values in an open collaborative web system?

Answer

The compound intake estimator calculates intakes and levels of health concern of plant-based food supplements based on product-specific concentration data and usage information provided by the user. The tool works online simply with a web browser. For user instructions and some data used, see below. The tools is openly available online at http://en.opasnet.org/w/Compound_intake_estimator.

Compound intakes and risk estimates

A concern indicator is calculated for each compound separately. If the intake is smaller than the threshold of concern, the indicator is < 1 and there is no concern. The threshold of concern is acceptable daily intake (ADI), tolerable daily intake (TDI), lower confidence limit of benchmark dose (BMDL10), or threshold of toxicological concern (TTC), depending on what data is available for each compound.

Compound concentrations

The following code allows to calculate the mean and median value of the concentration of each compound in the input data.

For seeing details of what the code actually does, click Show code. The code will appear, with documentation of each part on rows starting with ###. In brief, the code reads a piece of concentration and other data from Opasnet database, selects the rows defined by the user through the interface above, and shows the data on a summary table and a few graphs.

Rationale

For successful benefit-risk assessment of plant-based food supplements (PFSs), there is a need to estimate intakes of relevant compounds in a product. The health impacts of this exposure are then estimated based on available exposure-response functions. Optimally, all health impacts can be aggregated into e.g. disability-adjusted life years (DALYs) using disability weights or similar importance weights. However, the tool presented here does not go so far. Rather, it compares the intakes of each relevant compound to a guidance value, which is thought to separate safe intakes from non-safe intakes. Non-safe does not mean that the compound is causing health problems; it merely tells about compromised safety margins that are inherently embedded in the guidance values. Small exceedances are most likely harmless, but if the intake is clearly beyond its guidance value, careful consideration of potential risk is warranted.

This web tool, Compound intake estimator, was developed in Plantlibra project to promote the use of quantitative estimates of compound intakes and health risks. The tool consists of open data about guidance values (in the Opasnet database), an open source model (running on R software) to compute intakes and risk estimates, product-specific concentration data that is only partially available, and a user interface where the user can choose modelling options and add case-specific compound concentration data as input to the model. Detailed user guide, and the data and the model are available at http://en.opasnet.org/w/Compound_intake_estimator.

Users who are skilful with R can download the whole model to their own computers and use and develop it further as they please. This requires an R package called OpasnetUtils, which is available from the CRAN repository http://www.cran.r-project.org. It is also possible to develop a new page in Opasnet for a new product and its concentration data. Those who are interested in such work should contact the tool developer Jouni Tuomisto at THL, Finland.

Guidance values

There are several possible approaches for deriving guidance values. In this approach, we use acceptable daily intake (ADI), tolerable daily intake (TDI), lower confidence limit for benchmark dose 10 (BMDL10), and threshold for toxicological concern (TTC). If there are several guidance values available for a single compound, they are used in the preference order mentioned here.

Acceptable daily intake (ADI): An estimate of the amount of a chemical or toxin that a person can ingest daily over a human’s lifetime without appreciable health risk, divided by an average person’s lifespan. The concept of the ADI was developed by the WHO and FAO for chemicals such as additives to foods, pesticide residues and veterinary drugs in foods; ADIs are derived from lab toxicity data and from human experiences of such chemicals when available, and incorporate a safety factor; the ADI is an estimate of the amount of a substance in food that can be ingested over a lifetime by humans without significant risk to health.^[2]

Tolerable daily intake (TDI): Refers to the daily amount of a chemical that has been assessed safe for human being on long-term basis (usually whole lifetime). Originally acceptable daily intake (ADI) was introduced in 1961 to define the daily intake of a food additive which, during the entire lifetime, appears to be without appreciable risk. For contaminants and other foreign chemicals not used intentionally, the term TDI is often preferred.^[3]

Lower confidence limit for benchmark dose 10 (BMDL10): For the dose-response function of a health impact, it is possible to estimate the dose that causes on average 10 % response compared with the highest achievable response. This is called the benchmark dose (BMD10). BMDL10 is the lower confidence limit of BMD10, typically the 2.5th percentile of the probability distribution of the BMD10 estimate.^[4]

Threshold for toxicological concern (TTC): According to the TTC concept, a "safe" level of exposure can be identified for many chemicals based on their chemical structure and the known toxicity of chemicals that share similar structural characteristics. The TTC approach is exclusively designed as a substitute for substance-specific information in situations where there is limited or no information on the toxicity of the compound and information on exposure indicates that human exposure is very low.^[5]

If a compound in a PFS does not have any information on these guidance values, it is ignored in this analysis. The same happens to compounds that exist in the product but do not have any concentration information.

Intake and guidance value calculations

The total intake of essential oil is calculated like this:

intake_oil = amount_suppl * fraction_oil + amount_oil

where

intake is the total intake of essential oil (or another product part with the compounds of interest)
amount_suppl is the intake of food supplement
fraction_oil is the fraction of essential oil in the food supplement
amount_oil is the intake of the essential oil consumed directly, i.e. in addition to what is in the food supplement.

The intakes of individual compounds are calculated like this:

intake_comp (mg /kg /d) = intake_oil (g /d) * conc_comp (mg /g) / 60 (kg)

where

intake_comp are the intakes of each compound
intake_oil is the total intake of essential oil (or another product part with the compounds of interest)
conc_comp are the concentrations of each compound in the essential oil.
60 kg is the assumed body weight of the person

The concern indicator is calculated like this:

concern indicator_comp,i = intake_comp (mg /kg /d) / guidance value_comp,i (mg /kg /d) * safety factor_i

where

concern indicator_comp,i are the preferred indicators for each compound,
guidance value_comp are the guidance values for each compound,
safety factor_i are the safety factors used for each indicator (10000 for BMDL10, 1 for others, which already include a safety factor)

The concern indicator is calculated for each compound separately. If the intake is smaller than the guidance value, the indicator is < 1 and there is no concern. Only one guidance value indicator i is used for each compound, the preference order (when there are several guidance values for a compound) is explained in the previous section.

For BMDL10, the typical approach is to divide BMDL10 by the intake and interpret values > 10000 as safe. However, it would be confusing for the user to treat one guidance value differently than others. Therefore, this tool applies the same equation to all guidance values and uses 10000 as the default safety factor for BMDL10. For all other guidance values, the default safety factor is 1, because they already have embedded safety factors for interspecies and intraspecies variation.

User interface

You can fill in the following information:

Name of plant food supplement (PFS) product. If this is given, it is shown on the graphs and outputs as title.
Fraction of essential oil in the supplement. How much (as a percentage) does the total product contain the essential oil?
Amount of supplement per day. How much is the supplement (the total product) consumed per day (given as grams of product per day)?
The compound concentration data is taken from a table in the Opasnet database. However, if you know that the product contains something that is not listed in the table, you may add the information here.
- Name of additional compound. The name is used to combine the concentration data with toxicological data. So, the name should be found in one of these two tables on this page: Cramer classes of compounds, or Guidance values. If the name is not found, the new piece of information is ignored.
- Concentration of the additional compound in the supplement. This is given as mg of compound per g of PFS product.
- Technical note: You should enter the name of a single compound between quotation marks, e.g. 'Estragol'. If you add several compounds, you should use vector formatting for R, e.g. c('Estragole', 'Fenchone'). The concentration value can be without quotation marks unless it is a distribution or a vector. Acceptable entries for concentrations are e.g.: 2.1; '21.4-24.5'; and c('2.1', '21.4-24.5').
Show intermediate results? If Yes is selected, you will get two tables of a) concentration and b) intake estimates with all compounds in the product. Otherwise, only a default output is shown (see below).
Which extraction techniques to consider? Different extraction methods may result in different concentration estimates even from the same product. Some of the toxicological information is given having a particular extraction method in mind, and it is not always possible to use toxicological information with any concentration information. Therefore, if you know that some extraction methods should not be used to estimate toxicology, you can unselect them from the list. The default is to use all data, but you can select Yes, let me choose to the extraction techniques, and then select from the list that appears.
The default outcome shows a table and a graph where the level of concern is shown for each compound. However, because there are dozens of compounds in any product, the graph shows only those that have the level more than one, i.e. those compounds that need further scrutiny.

Concentrations

The concentration code allows to calculate the mean and median value of the concentration of each compound that appear in table. In addition, it is possible to plot the concentration (expressed as percentage of essential oil) of each compound reported in the data table. Each single value is reported in the chart.

Data

This is an example data to show what kind of information is used by the tool. It is also used to calculate the example results. Some of the actual data for PFSs is proprietary and is currently only available to the partners of the Plantlibra project.

Example data(%)
Obs	Extraction.technique	Compound	Result	Reference
1	Superheated water extraction, 100 °C	2-Carene	0.08; 0.1; 0.12	Jayawardena & Smith, 2010
2	Superheated water extraction, 100 °C	Z-Cinnamaldehyde	1.63; 2.1; 2.57	Jayawardena & Smith, 2010
3	Superheated water extraction, 100 °C	Cinnamaldehyde	81.14; 83.7; 86.26	Jayawardena & Smith, 2010
4	Superheated water extraction, 100 °C	Eugenol	0.45; 0.8; 1.15	Jayawardena & Smith, 2010
5	Hydrodistillation	1,8-Cineole	0.2	Chericoni et al, 2005
6	Superheated water extraction, 100 °C	Cinnamyl acetate	5.47; 7.2; 8.93	Jayawardena & Smith, 2010

Calculations

This is the code for the actual compound intake estimator model. You need to run the code below only if you update the guidance value tables on this page.

+ Show code - Hide code


library(OpasnetUtils)

################## PRE-DEFINED INPUTS

# compound.conc is the concentration of a compound in the critical constituent (or food, if fr.const == 1) 
# fr.const is the fraction of the critical constituent in the whole food
# extraction is a criteria to select specific extraction techniques
# addcomp is an additional compound not in the product. Data about its concentration in food is given by the user.

# compound.intake is the intakes of the compounds via food or food supplement, each given as mg /kg /d.

compound.intake <- Ovariable("compound.intake", 
	dependencies = data.frame(Names = c(
		"food.amount", # amount of food supplement consumed
		"compound.conc", # concentration of compounds in essential oil
		"oil.amount", # amount of essential oil consumed in addition to the supplement product
		"fr.const", # fraction of essential oil in the supplement product
		"extraction", # Which extraction methods the user wants to consider
		"addcomp" # concentrations of compounds in the essential oil (missing from default data)
		# conc is the temporary ovariable for compound.conc used in calculations
	)), 
	formula = function(...) {

		conc <- compound.conc * 10 * extraction 
		# remove extractions that are not considered.
		# compound.conc from w/w% to mg/g

#		oapply would be an alternative to remove extra indices, but there is a bizarre error.
#		compound.conc <- oapply(compound.conc, cols = "Extraction.technique", fun = mean)

		test <- colnames(conc@output) != "Extraction.technique"
		conc@marginal	<- conc@marginal[test]
		conc@output 	<- conc@output[test] # unit: mg /g
		# removes index Extraction.technique, which is no longer needed.

		conc@output <- orbind(conc, addcomp * 1)
		# addcomp must not contain columns that are not in conc. It must contain Compound.
		# It is multiplied by 1 to make the result column name "Result".

		out <- (food.amount * fr.const / 100 + oil.amount) * conc / 60 
		# calculates the oil intake from food supplement and then adds direct oil intake.
		# scale to 60 kg person
		# fr.const from % to fraction

		test2 <- unique(out@output$Compound[result(out) > 0]) # remove all compounds where intake = 0
		out@output <- out@output[out@output$Compound %in% test2 , ]

		return(out)
	}
)

# guidance.values is the guidance value (ADI, BMDL10, TTC or similar) for compounds. unit: mg /kg /d
# cramer is cramer classes of compounds (temporary)
# cramer.values are the thresholds of toxicological concern for each Cramer class (temporary)

cramer <- tidy(opbase.data("Op_en6193", subset = "cramer_classes_of_compounds"))
cramer$Result <- NULL
cramer$Guidance <- "TTC" # Threshold for toxicological concern

cramer.values <- tidy(opbase.data("Op_en6193", subset = "cramer_values"))
cramer.values <- merge(cramer, cramer.values)
cramer.values[["Cramer class"]] <- NULL
cramer.values$Result <- cramer.values$Result / 1000 # from ug /kg /d to mg /kg /d

guidance.values <- tidy(opbase.data("Op_en6193", subset = "guidance_values"))
guidance.values <- orbind(cramer.values, guidance.values)
guidance.values <- Ovariable("guidance.values", data = guidance.values)

safety.factor <- Ovariable("safety.factor", ddata = "Op_en6193.safety_factors") #  [[Compound intake estimator]]

# food.risk calculates "hazard quotients" or risks relative to guidance values. Value < 1 are of minimal concern

food.risk <- Ovariable(name = "food.risk",
	dependencies = data.frame(Name = c("compound.intake", "guidance.values", "safety.factor")),
	formula = function(...) {

		# guidance.risk is intake compared with guidance values.

		out <- compound.intake / (guidance.values / safety.factor)

		removepoorguidance <- function(dat, goodg) { # Remove inferior guidance values

			# Compounds that have good guidance values (goodg).
			hasgoodg <- unique(dat[dat$Guidance == goodg , "Compound"])

			# The row contains data about a chemical with a good guidance value
			# but the data is about an inferior guidance value.

			out <- dat[!(dat$Compound %in% hasgoodg & dat$Guidance != goodg) , ]

			return(out)
		}

		out <- removepoorguidance(out@output, "ADI")
		out <- removepoorguidance(out, "TDI")
		out <- removepoorguidance(out, "BMDL10")

		test2 <- unique(out$Compound[out$Result > 0]) # remove all compounds where risk = 0
		out <- out[out$Compound %in% test2 , ]

		return(out)
	}
)

extraction <- Ovariable("extraction", ddata = "Op_en6193.extraction_techniques")
extraction@data$extractionResult <- 1

### Define formatting for the concentration graphs that will be shown.

concplot <- function(condition) {
	dat <- compound.conc@output[condition, ]
	if(nrow(dat) != 0) {
		out <- ggplot(dat, 
			aes(
				x = Compound, 
				y = compound.concResult, 
				colour = Extraction.technique
			)) + 
			geom_point(size=5) +
			theme_grey(base_size = 24) +
			theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
			labs(y = "Concentration (%)")
		return(out)
	} else {
		 warning("No data to show")
	}
}

objects.store(compound.intake, guidance.values, food.risk, extraction, concplot, safety.factor)

cat("The following objects were successfully stored:
compound.intake, guidance.values, food.risk, extraction, concplot, safety.factor\n")

Using different data on different layers: [7].

Extraction techniques

Extraction techniques(-)
Obs	Extraction.technique	Notes
1	Not known
2	Hydrodistillation
3	Solvent extraction
4	Steam distillation
5	Supercritical fluid extraction, 120 bar, 40 °C
6	Supercritical fluid extraction, 90 bar, 40 °C
7	Supercritical fluid extraction, 90 bar, 50 °C
8	Superheated water extraction, 100 °C
9	Superheated water extraction, 150 °C
10	Superheated water extraction, 200 °C
11	Solid phase
12	70% ethanol (1:5)
13	Acetone and water mixture (35 to 67:1)
14	Commercial product
15	Methanol (40%) and water
16	Not reported
17	Ultrasound assisted extraction
18	70% methanol, Ultrasound assisted extraction
19	ddH2O and Na2HPO4,partitioning with ethyl acetate, organic phases dried and redissolved in methanol
20	Maceration
21	Soxhlet extraction

Data about guidance values

Generic data about guidance values

Safety factors(-)
Obs	Guidance	Safety factor
1	ADI	1
2	BMDL10	10000
3	TDI	1
4	TTC	1

Cramer values(µg /kg-BW /d)
Obs	Cramer class	Threshold of toxicological concern	Description
1	1	30
2	2	1.5
3	3	1.5
4	4	0.0025	This value is the TTC for compounds having a structural alert for genotoxicity

Values 30, 1.5 and 0.025 µg /kg /d come from EFSA (2012)^[6]

Compound-specific data about guidance values

^[7]

Guidance values(mg /kg-BW /d)
Obs	Compound	Guidance	Result
1	1,8-Cineole	ADI	2.8
2	4-Terpineol	ADI	1.2
3	Benzylbenzoate	ADI	5
4	Cinnamaldehyde	ADI	0.7
5	Coumarin	TDI	0.1
6	Estragole	BMDL10	3.3-6.5
7	Eugenol	ADI	2.5
8	Fenchone	ADI	10.64
9	Limonene	ADI	1
10	Linalool	ADI	0.5
11	Myrcene	ADI	2.5
12	Trans-anethole	ADI	2
13	Safrole	BMDL10	1.9-5.1
14	α-Terpineol	ADI	1.2

Cramer classes of compounds

Cramer classes of compounds(-)
Obs	Compound	Cramer class	Notes
1	1,8-Cineole	3	A new line was created for 1.8-Cineole based on "ß-Phellandrene + 1,8-Cineole: 1+3"
2	2-Carene	1
3	2-Methoxycinnamaldehyde	1
4	3-Phenylpropyl acetate	1
5	4-Methoxyphenylacetone	1
6	4-Terpineol	3
7	5-Caffeoylquinic Acid	3
8	Acacetin	3
9	Allo-ocimene	1
10	Apigenin	3
11	Apiole	4
12	Aromadendrene	1
13	Aromatic compound		Class could not be assigned
14	Benzaldehyde	1
15	Benzenepropanal	1
16	Benzenepropanol acetate	1
17	Benzenepropanol	1
18	Benzylbenzoate	1
19	Bilobalide	3
20	Borneol	1
21	Cadinene	1
22	Calamenene	2
23	Camphene	1
24	Camphor	3
25	Carvacrol	1
26	Carveol isomere	1	Carveol (2-methyl-5-prop-1-en-2-ylcyclohex-2-en-1-ol) is 1. Isomere is not specified so the class was assigned to the general compound
27	Caryophyllene alcohol	3
28	Caryophyllene oxide	3
29	Caryophyllene	1
30	Chavicol	1
31	Cinnamaldehyde	1
32	Cinnamic acid	1
33	Cinnamyl acetate	1
34	Cinnamyl alcohol	1
35	Cis-anethole	2
36	Cis-limonene oxide	2
37	Cis-P-menth-2-en-1-ol	3
38	Cis-sabinene hydrate	2
39	Cis-ß-Ocimene	1
40	Cis-α-ocimene	1
41	Cis-β-ocimene	1
42	Coumarin	3
43	Cryptone	2
44	Cymene isomere	1	cymene is 1. Isomere is not specified so the class was assigned to the general compound
45	E-Anethole	1
46	e-Cadinene + Pioncarvone		Class could not be assigned
47	E-Cinnamyl acetate	1
48	E-o-Methoxycinnamaldehye	1
49	Epoxy-6.7-Humulene + Mw=220		Class could not be assigned
50	Estragole	1	Test compound, not real value
51	Eugenol acetate	3
52	Eugenol	1
53	Fenchyl acetate	1
54	Ginkgolide A	3
55	Ginkgolide B	3
56	Ginkgolide C	3
57	Ginkgolide J	3
58	Ginkgolide A, B, C	3
59	Ginkgolide A, B, C, J	3
60	Humulene epoxide II	3
61	Isoledene	1
62	Isorhamnetin	3
63	Kaempferol	3
64	Limonene	1
65	Linalool	3
66	Luteolin	3
67	Methoxy eugenol	3
68	Methyl cinnamaldehyde	1
69	Methyl eugenol	1
70	Monoterpenes		Class could not be assigned
71	Myricetin	3
72	N-Hexadecanoic acid	1
73	N-Octadecanoic acid	1
74	o-Cymene	1
75	p-Anisaldehyde	1
76	Para-Methoxy cinnamic aldehyde	1
77	p-Cymene	1
78	P-Cymene-8-ol	3
79	Phenyl cyclohexene	2
80	Phenylmethyle isovalerate	1
81	Quercitin	3
82	Rutin	1	Test compound, not real value
83	Sabinene	1
84	Safrole	3
85	Sesquiterpene epoxide		Class could not be assigned. Sesquierpene is class 3, but don't now if something changes with epoxide
86	Sesquiterpenol Mw=220		Class could not be assigned
87	Sesquiterpenol Mw=222		Class could not be assigned
88	Spathulenol	3
89	ß-Bisabolene	1
90	ß-Cubebene	1
91	ß-Elemene	1
92	ß-Myrcene	1
93	ß-Phellandrene	1	ß-Phellandrene + 1,8-Cineole: 1+3
94	ß-Pinene	1
95	Styrene	1
96	Terpinen-4-ol	3
97	Terpinolene isomere	1	terpinolene is 1 , isomere not known. Isomere is not specified so the class was assigned to the general compound
98	Terpinolene	1
99	Tetradecanal	1
100	Trans-2-methoxycinnamaldehyde	1
101	Trans-limonene oxide	2
102	Trans-pinene hydrate	2
103	Trans-Pinocarveol	1
104	Trans-Piperitol	1
105	Trans-P-menth-2-en-1-ol	3
106	Trans-ß-Ocimene	1
107	Z-Cinnamaldehyde	1
108	α-Cadinol	3
109	α-Calacorene	2
110	α-Copaene	1
111	α-Cubebene	1
112	α-Famesene		Class could not be assigned
113	α-Gurjunene	1
114	α-Humulene	1
115	α-Muurolene	1
116	α-Muurolol	3
117	α-P-Dimethylstyrene	1
118	α-Phellandrene	1
119	α-Pinene	1
120	α-Selinene	1
121	α-Terpinene	1
122	α-Terpineol	3
123	α-Thujene	1
124	β-Anise aldehyde	1
125	β-Phellandrene	1
126	β-Pinene	1
127	γ-Cadinene	1
128	γ-terpinene	1
129	γ-Terpinene	1
130	Δ3-Carene	1
131	δ-Cadinene	1

Cramer class data from User:Anto.guzzon.

References

↑ PlantLIBRA website
↑ Segen's Medical Dictionary. 2012 Farlex, Inc. [1] accessed 15th April, 2014.
↑ Wikipedia: Tolerable daily intake accessed 15th April 2014
↑ U.S.EPA: Benchmark Dose (BMD) Methodology [2] accessed 15th April, 2014.
↑ Opinion on the Use of the Threshold of Toxicological Concern (TTC) Approach for Human Safety Assessment of Chemical Substances with focus on Cosmetics and Consumer Products. European Union (2012) SCCP/1171/08. [3] accessed 15th April 2014.
↑ EFSA (2012). Scientific Opinion on Exploring options for providing advice about possible human health risks based on the concept of Threshold of Toxicological Concern (TTC). EFSA Journal 2012;10(7):2750
↑ van den Berg SJPL, Restani P, Boersma MG, Delmulle L, Rietjens IMCM (2011) Levels of Genotoxic and Carcinogenic Ingredients in Plant Food Supplements and Associated Risk Assessment. Food and Nutrition Sciences, 2(9): 989-1010.

Compound concentration data for PlantLIBRA cases

PlantLIBRA WP5 (HYLO) collected data about compound concentrations in three case study plants. The methods for the data collection are described elsewhere.

Cinnamon (dried bark)
Bitter fennel (fruit)
Gingko (leaf)

The data were uploaded to Opasnet database. The database is a versatile data storage for all kinds of data tables related to benefit-risk assessment, health, food, and other topics. It has two parts: an open part that can be accessed using a web browser or through a machine-readable JSON interface; and a password-protected part that is technically similar but is not visible to the Internet.

Data on composition were collected running a systematic search on Pubmed, using information available from previous assessments of authoritative bodies (e.g. EFSA, EMA, etc.) and the ePlantLIBRA database, which was developed in this project. ePlantLIBRA is based on eBASIS from the EuroFIR Network. ^[1] eBASIS was produced by the EuroFIR network and is now being further developed in a co-operation between the PlantLIBRA Project and EuroFIR AISBL. EuroFIR is funded by the European Commission under the Food Quality and Safety Priority of the FP7 for Research and Technological Development (CT).

eBASIS is a unique database that contains critically evaluated published data on the content and biological effects of bioactive constituents in plant based foods and an up-to-date list of plant and plant part names in 15 EU languages. Bioactive Compounds are defined as inherent non-nutrient constituents of food plants and edible mushrooms with anticipated health promoting, beneficial and/or toxic effects when ingested. ^[2]

In an early stage there was a plan to store the compound data only in the ePlantLIBRA database and access that dynamically through a machine-readable interface from Opasnet. However, this plan was abandoned for several reasons. The main reason was that ePlantLIBRA database is actually only a part of the larger eBASIS database, which contains also proprietary, non-public information. It would have been a large work to build an additional security functionality to allow access to some but not all parts of the database. In addition, the tool we are describing here needs a moderate amount of data, which can easily be copied and uploaded to another system (Opasnet database) for this purpose. And thirdly, when new compound concentration data becomes available (either from ePlantLIBRA or elsewhere), it is straightforward to upload also the new data to Opasnet database as a part of Compound intake estimator. Thus, the use of Compound intake estimator is not dependent on ePlantLIBRA database but can benefit from its future development.

Case study results

The following figures demonstrate some output graphs from the model. The graphs give a user-friendly way to get an overview of the compounds in a product. It is especially interesting to compare concentrations and levels of concern of different compounds. As can be seen from the pair of graphs of cinnamon dried bark essential oil, most of the product is cinnamaldehyde, which is only marginally above a guidance value and thus not likely a concern. In contrast, safrole seems to be somewhat of concern or at least a compound to pay attention to, as it goes beyond a guidance value by a factor of 10 - 80.

Bitter fennel is a different story. The clearly largest concern is estragole, which goes beyond the BMDL10 by a factor of 50 - 600. Interestingly, this compound is present in low concentrations when trans-anethole is present at high concentrations and vice versa. Thus, in the first case estragole might be causing a clear risk and fennel products would therefore need a a benefit-risk assessment. There is also another interesting compound, namely apiole, which goes above the level of concern by a factor of 10 - 200, even when its concentration is so low (0.1 %) that it is not even shown on the graph. With this compound, the first question in mind is whether there is actually evidence supporting its much higher putative potency compared with the other compounds.

Discussion

Concentration and intake data are not easily understood by consumers. There are dozens or hundreds of compounds in a single PFS product. It is therefore very difficult to get any idea which compounds are important and whether a consumer should be interested in further information about these compounds. There is a need for tools visualising these data to promote informed discussion and consumer decisions. The model presented here offers a solution to this need.

This model will make it easy to compare two critical things: exposure to compounds relative to health-based guidance values, and concentrations in the product. By comparing these two graphs, it was possible to identify from each example product two critical compounds that would warrant further scrutiny. The reasons for need of further work were different with each compound, as described in the previous section.

The compound intake estimator model can be used online simply with a web browser. It can easily be copied to other projects or uses in Opasnet. It can also be used on one's own computer, but that requires some skills of R software. It is a small task to add new data sets to the model when they become available. It is open source and thus it can be taken and developed further. It is a part of a large web workspace Opasnet that has several independent users and continuous, long-term maintenance. This flexibility and durability is important when considering the lifespan of the model. Even after the end of the PlantLIBRA project, the model will be maintained and kept available for users.

EFSA data collection and PFS

The critical thing with the compound intake estimator model is the availability of relevant compound concentration data about PFSs. EFSA has recently studied availability of such data. The EFSA contract CFT/EFSA/DCM/2011/03 provided EFSA with an updated food composition database covering approximately 1750 foods in combination with additional FoodEx2 facet descriptors included in the EFSA FoodEx2 classification system. A model for data transfer, compatible with the EuroFIR technical annex and CEN Food Data Standard, and the EFSA data structure was developed and tested. Fourteen national food database compiler organisations supplied initial food lists mapped to the EFSA food list. Eight of them contained data about food supplements and their composition. Datasets compatible with EFSA's data structure were produced based on the models. ^[3]

Therefore, there exists a large amount of potentially useful PFS compound concentration data in EFSA. However, EFSA has not made this data available for open use, and this will not probably happen before the end of PlantLIBRA project. In any case, it is important to acknowledge and remember this possibility and potential for improvement of the model described here.

**Supplement data provided by compilers.^[3]**
Country	Number of supplements listed
Finland	92
France	449
Italy	126
Netherlands	90²
Serbia	8
Slovenia	95
Sweden	122
UK	103

²Does not include supplements for children.

References

↑ ePlantLIBRA website, Eurofir, accessed 15th April, 2014.
↑ eBASIS (Bioactive Substances in Food Information System), Eurofir, accessed 15th April, 2014.
↑ ^3.0 ^3.1 M.A. Roe, S. Bell, M. Oseredczuk, T. Christensen, S. Westenbrink, H. Pakkala, K. Presser and P.M. Finglas, 2013. Updated food composition database for nutrient intake. EFSA supporting publication 2013:EN-355, 21 pp. Available online: [4] [5] [6]

Appendix 1: Guidance values and other inputs

Move the tables about guidance value to this location.

Appendix 2: Computer code for the tool

Compound intakes and risks

Show details

<rcode label='Intakes and risks' embed=1 graphics=1
variables="
	name:name|description:Name of plant food supplement|type:text|default:Cinnamon essential oil|
category:Intake information|
	name:food.amount|description:Amount of supplement per day (g /d)|type:text|default:1|
	name:fr.const|description:Fraction of essential oil in the supplement (%)|type:text|default:10|
	name:oil.amount|description:Amount of essential oil consumed directly per day (g /d)|type:text|default:0|
	name:intermediates|description:Show intermediate results?|
		type:selection|options:FALSE;No;TRUE;Yes|default:FALSE|
	name:extract|description:Do you want to limit to some extraction techniques?|type:selection|
		options:FALSE;No, analyze all techniques;TRUE;Yes, let me choose|default:FALSE|
	name:add|description:Do you want to add own compound data?|
		type:selection|options:FALSE;No;TRUE;Yes|default:FALSE|
	name:addcompound|description:Name(s) of additional compound(s)|default:' '|
category:Additional compound information|
category_conditions:add;TRUE|
	name:addconc|description:Concentration(s) of the compound(s) (mg of compound /g product)|default:0|
	name:extracted|description:Which extraction techniques to consider?|type:checkbox|options:
		1;Not known;
		2;Hydrodistillation;
		3;Solvent extraction;
		4;Steam distillation;
		5;Supercritical fluid extraction, 120 bar, 40 °C;
		6;Supercritical fluid extraction, 90 bar, 40 °C;
		7;Supercritical fluid extraction, 90 bar, 50 °C;
		8;Superheated water extraction, 100 °C;
		9;Superheated water extraction, 150 °C;
		10;Superheated water extraction, 200 °C;
		11;Solid phase|
		default:1;2;3;4;5;6;7;8;9;10;11|
	category:Extraction technique|
	category_conditions:extract;TRUE
">

library(OpasnetUtils)
library(ggplot2)

################## SETUP MODEL

openv.setN(100) # Sets the number of model iterations to 100.

# Fetches the ovariables needed in the model from [[Compound intake estimator]].
# Ovariables downloaded: guidance.values, compound.intake, food.risk, extraction, concplot.

objects.latest("Op_en6193", code_name="initiate") 

# compound.conc is the concentration of compounds in the critical constituent.
# Data comes from [[Compound intake estimator]].

compound.conc <- Ovariable(name = "compound.conc", ddata = "Op_en6193.example_data")

################## USER INPUTS

# name is the name of product. It is used in the output graphs.

# food.amount is the amount that the consumer consumes the defined product (g /day)

food.amount <- Ovariable("food.amount", data = data.frame(Result = food.amount))

# fr.const is the fraction of a critical constituent in the food or food supplement (given as percentage)

fr.const <- Ovariable("fr.const", data = data.frame(Result = fr.const))

# oil.amount is the amount that the consumer consumes essential oil in addition to the product (g /day)

oil.amount <- Ovariable("oil.amount", data = data.frame(Result = oil.amount))

# extraction is a list of extraction techniques to consider. (by default, all extraction techniques)

if(exists('extracted')) {
	extraction@data <- extraction@data[extracted , ] # Order of extraction techniques must not change!
} else {
	extraction@data <- data.frame(extractionResult = 1) # Remove column Extraction.techniques
}

# addcomp are user-defined additional compounds for a supplement, if they are not in the table.
## addcompound is the name and addconc is the concentration of the additional compound.
# addcomp and addcompound can be vectors.

addcompound <- ifelse(! exists('addcompound') || addcompound == ' ', NA, addcompound)
addconc <- ifelse(! exists('addconc') || addconc == ' ', 0, addconc)

addcomp <- Ovariable("addcomp", data = data.frame(Compound = addcompound, Result = addconc))


############## ACTUAL MODEL

food.risk <- EvalOutput(food.risk)

################ PRINTING MODEL OUTPUTS

if(intermediates) {

	cat("Concentrations of compounds in ", name, " (mg /kg DW)\n")

	oprint(summary(compound.conc, marginals = c("Compound")))

	cat("Intakes of different compounds (mg /kg BW /d)\n")

	oprint(summary(compound.intake, marginals = c("Compound")), digits = 3)

}

cat("Risks relative to guidance values (ADI, TDI, BMDL10, or TTC). Values < 1 are of minimal concern.\n")

oprint(summary(food.risk, marginals = c("Compound")))

test <- unique(food.risk@output[food.risk@output$food.riskResult >= 1 , "Compound"])

temp <- food.risk@output[food.risk@output$Compound %in% test , ]

if(nrow(temp) == 0) {
	cat("There are no compounds that exceed a guidance value.\n")
	stop("There are no compounds that exceed a guidance value.\n")
}

cat("Risks relative to guidance values (ADI, TDI, BMDL10, or TTC). Values < 1 are of minimal concern and are not shown in the 
graph. Note! BMDL10 is used to calculate margin of exposure (MOE), but instead of MOE, (10000 / MOE) is shown on this graph 
for consistency.\n")

ggplot(temp, aes(x = Compound, y = food.riskResult, fill = Guidance)) + 
	geom_boxplot() + 
	theme_grey(base_size = 24) +
	labs(title = paste("Exposure to compounds in", name), y = "Ratio (Exposure /Guidance)") +
	theme(axis.text.x = element_text(angle = 90, hjust = 1))

</rcode>

Compound concentrations

Show details

<rcode label='Plot concentrations' graphics=1 embed=1
variables="
	name:limit|description:Plot values above this limit (%)|default:10|
	name:allcompounds|description:Do you want to plot all compounds?|type:selection|options:TRUE;Yes;FALSE;No, choose from the list below|default:TRUE|
	name:allextractions|description:Do you want to plot all extraction techniques?|type:selection|options:TRUE;Yes;FALSE;No, choose from the list below|default:TRUE|
	name:compound|description:If not, what compounds do you want to plot?|type:checkbox|options:
		'1,8-Cineole';1,8-Cineole;
		'2-Carene';2-Carene;
		'2-Methoxycinnamaldehyde';2-Methoxycinnamaldehyde;
		'2-Methoxycinnamaldehye';2-Methoxycinnamaldehye;
		'2-Methoxycinnamaldehye cinnamalaldehyde';2-Methoxycinnamaldehye cinnamalaldehyde;
		'3-Phenylpropyl acetate';3-Phenylpropyl acetate;
		'4-Terpineol';4-Terpineol;
		'a-Copaene';a-Copaene;
		'a-Muurolene';a-Muurolene;
		'a-Terpineol';a-Terpineol;
		'Aromadendrene';Aromadendrene;
		'Aromatic compound';Aromatic compound;
		'Benzaldehyde';Benzaldehyde;
		'Benzaldeyhde';Benzaldeyhde;
		'Benzene propanol';Benzene propanol;
		'Benzenepropanal';Benzenepropanal;
		'Benzenepropanol acetate';Benzenepropanol acetate;
		'Benzylbenzoate';Benzylbenzoate;
		'Borneol';Borneol;
		'Cadinene';Cadinene;
		'Calamenene';Calamenene;
		'Camphene';Camphene;
		'Camphor';Camphor;
		'Carvacrol';Carvacrol;
		'Carveol isomere';Carveol isomere;
		'Caryophyllene';Caryophyllene;
		'Caryophyllene alcohol';Caryophyllene alcohol;
		'Caryophyllene oxide';Caryophyllene oxide;
		'Caryophyllene oxyde';Caryophyllene oxyde;
		'Chavicol';Chavicol;
		'Cinamyl acetate';Cinamyl acetate;
		'Cinnamaldehyde';Cinnamaldehyde;
		'Cinnamic acid';Cinnamic acid;
		'Cinnamyl acetate';Cinnamyl acetate;
		'Cinnamyl alchol';Cinnamyl alchol;
		'Cinnamyl alcohol';Cinnamyl alcohol;
		'Cinnamyl alcool';Cinnamyl alcool;
		'Cis-P-menth-2-en-1-ol';Cis-P-menth-2-en-1-ol;
		'Cis-ß-Ocimene';Cis-ß-Ocimene;
		'Coumarin';Coumarin;
		'Cryptone';Cryptone;
		'Cymene isomere';Cymene isomere;
		'E-Anethole';E-Anethole;
		'e-Cadinene + Pioncarvone';e-Cadinene + Pioncarvone;
		'E-Cinnamyl acetate';E-Cinnamyl acetate;
		'E-o-Methoxycinnamaldehye';E-o-Methoxycinnamaldehye;
		'Epoxy-6.7-Humulene + Mw=220';Epoxy-6.7-Humulene + Mw=220;
		'Eugenol';Eugenol;
		'Eugenol acetate';Eugenol acetate;
		'Humulene epoxide II';Humulene epoxide II;
		'Isoledene';Isoledene;
		'Limonene';Limonene;
		'Linalool';Linalool;
		'Linaool';Linaool;
		'Methoxy eugenol';Methoxy eugenol;
		'Methyl cinnamaldehyde';Methyl cinnamaldehyde;
		'Methyl eugenol';Methyl eugenol;
		'Monoterpenes';Monoterpenes;
		'N-Hexadecanoic acid';N-Hexadecanoic acid;
		'N-Octadecanoic acid';N-Octadecanoic acid;
		'P-Cymene';P-Cymene;
		'P-Cymene-8-ol';P-Cymene-8-ol;
		'Para-Methoxy cinnamic aldehyde';Para-Methoxy cinnamic aldehyde;
		'Phenyl cyclohexene';Phenyl cyclohexene;
		'Phenylmethyle isovalerate';Phenylmethyle isovalerate;
		'Sabiene';Sabiene;
		'Sabinene';Sabinene;
		'Safrole';Safrole;
		'Sesquiterpene epoxide';Sesquiterpene epoxide;
		'Sesquiterpenol Mw=220';Sesquiterpenol Mw=220;
		'Sesquiterpenol Mw=222';Sesquiterpenol Mw=222;
		'Spathulenol';Spathulenol;
		'ß-Bisabolene';ß-Bisabolene;
		'ß-Cubebene';ß-Cubebene;
		'ß-Elemene';ß-Elemene;
		'ß-Myrcene';ß-Myrcene;
		'ß-Phellandrene';ß-Phellandrene;
		'ß-Phellandrene + 1,8-Cineole';ß-Phellandrene + 1,8-Cineole;
		'ß-Pinene';ß-Pinene;
		'Styrene';Styrene;
		'Terpinen-4-ol';Terpinen-4-ol;
		'Terpiolene';Terpiolene;
		'Terpiolene isomere';Terpiolene isomere;
		'Tetradecanal';Tetradecanal;
		'Trans-2-methoxycinnamaldehyde';Trans-2-methoxycinnamaldehyde;
		'Trans-P-menth-2-en-1-ol';Trans-P-menth-2-en-1-ol;
		'Trans-Pinocarveol';Trans-Pinocarveol;
		'Trans-Piperitol';Trans-Piperitol;
		'Trans-ß-Ocimene';Trans-ß-Ocimene;
		'Z-Cinnamaldehyde';Z-Cinnamaldehyde;
		'α-Cadinol';α-Cadinol;
		'α-Calacorene';α-Calacorene;
		'α-Copaene';α-Copaene;
		'α-Cubebene';α-Cubebene;
		'α-Famesene';α-Famesene;
		'α-Gurjunene';α-Gurjunene;
		'α-Humulene';α-Humulene;
		'α-Muurolene';α-Muurolene; 
		'α-Muurolol';α-Muurolol;
		'α-P-Dimethylstyrene';α-P-Dimethylstyrene;
		'α-Phellandrene';α-Phellandrene;
		'α-Pinene';α-Pinene;
		'α-Selinene';α-Selinene;
		'α-Terpinene';α-Terpinene;
		'α-Terpineol';α-Terpineol;
		'α-Thujene';α-Thujene;
		'γ-Cadinene';γ-Cadinene;
		'γ-Terpinene';γ-Terpinene;
		'δ-Cadinene';δ-Cadinene;
		'Δ3-Carene';Δ3-Carene|
	category:Chemical compounds|
	category_conditions:allcompounds;FALSE|
	name:extraction|type:checkbox|options:
		'Hydrodistillation';Hydrodistillation;
		'Not known';Not known;
		'Solvent extraction';Solvent extraction;
		'Steam distillation';Steam distillation;
		'Supercritical fluid extraction, 120 bar, 40 °C';Supercritical fluid extraction, 120 bar, 40 °C;
		'Supercritical fluid extraction, 90 bar, 40 °C';Supercritical fluid extraction, 90 bar, 40 °C;
		'Supercritical fluid extraction, 90 bar, 50 °C';Supercritical fluid extraction, 90 bar, 50 °C;
		'Superheated water extraction, 100 °C';Superheated water extraction, 100 °C;
		'Superheated water extraction, 150 °C';Superheated water extraction, 150 °C;
		'Superheated water extraction, 200 °C';Superheated water extraction, 200 °C|
	category:Extraction technique|
	category_conditions:allextractions;FALSE
">
### Load some R packages needed to implement the code.

library(OpasnetUtils)
library(ggplot2)

################ SETUP THE MODEL

openv.setN(100) # Set model iterations

# Fetches the ovariables needed in the model from [[Compound intake estimator]].
# Ovariables downloaded: guidance.values, compound.intake, food.risk, extraction, concplot

objects.latest("Op_en6193", code_name="initiate") 

### Download the data from the database. You can also see the raw data on this page below,
### or access the database via this link: http://heande.opasnet.org/wiki/Special:Opasnet_Base?id=Heande3857

# compound.conc is the concentration of compounds in the critical constituent.
# Data comes from [[Compound intake estimator]].

compound.conc <- Ovariable(name = "compound.conc", ddata = "Op_en6193.example_data")

################ USER INPUTS

### Select only those rows that fulfill the criteria of compound and extraction method set by the user (if any).

dat <- compound.conc@data
if(!allcompounds) {dat <- dat[dat$Compound %in% compound , ]}
if(!allextractions) {dat <- dat[dat$Extraction.technique %in% extraction , ]}
compound.conc@data <- dat

### Create an information object (called ovariable) out of the data for easy handling.

compound.conc <- EvalOutput(compound.conc)

### Print a summary table of the output.

oprint(summary(compound.conc), digits = 3)

### Use the graph formatting for several graphs. The graphs only differ by the range of concentration that is shown on Y axis.

concplot(compound.conc@output$compound.concResult >= limit)
concplot(compound.conc@output$compound.concResult <= 10 & compound.conc@output$compound.concResult > 1)
concplot(compound.conc@output$compound.concResult <= 1 )

</rcode>

Model initiation

This is the code for the actual compound intake estimator model. This part of the model is stored online in Opasnet, and it is automatically fetched when the other parts of the model are run.

Show details

<rcode name="initiate" label="Initiate model" embed=1>

library(OpasnetUtils)

################## PRE-DEFINED INPUTS

# compound.conc is the concentration of a compound in the critical constituent (or food, if fr.const == 1) 
# fr.const is the fraction of the critical constituent in the whole food
# extraction is a criteria to select specific extraction techniques
# addcomp is an additional compound not in the product. Data about its concentration in food is given by the user.

# compound.intake is the intakes of the compounds via food or food supplement, each given as mg /kg /d.

compound.intake <- Ovariable("compound.intake", 
	dependencies = data.frame(Names = c(
		"food.amount", # amount of food supplement consumed
		"compound.conc", # concentration of compounds in essential oil
		"oil.amount", # amount of essential oil consumed in addition to the supplement product
		"fr.const", # fraction of essential oil in the supplement product
		"extraction", # Which extraction methods the user wants to consider
		"addcomp" # concentrations of compounds in the essential oil (missing from default data)
		# conc is the temporary ovariable for compound.conc used in calculations
	)), 
	formula = function(...) {

		conc <- compound.conc * 10 * extraction 
		# remove extractions that are not considered.
		# compound.conc from w/w% to mg/g

#		oapply would be an alternative to remove extra indices, but there is a bizarre error.
#		compound.conc <- oapply(compound.conc, cols = "Extraction.technique", fun = mean)

		test <- colnames(conc@output) != "Extraction.technique"
		conc@marginal	<- conc@marginal[test]
		conc@output 	<- conc@output[test] # unit: mg /g
		# removes index Extraction.technique, which is no longer needed.

		conc@output <- orbind(conc, addcomp * 1)
		# addcomp must not contain columns that are not in conc. It must contain Compound.
		# It is multiplied by 1 to make the result column name "Result".

		out <- (food.amount * fr.const / 100 + oil.amount) * conc / 60 
		# calculates the oil intake from food supplement and then adds direct oil intake.
		# scale to 60 kg person
		# fr.const from % to fraction

		test2 <- unique(out@output$Compound[result(out) > 0]) # remove all compounds where intake = 0
		out@output <- out@output[out@output$Compound %in% test2 , ]

		return(out)
	}
)

# guidance.values is the guidance value (ADI, BMDL10, TTC or similar) for compounds. unit: mg /kg /d
# cramer is cramer classes of compounds (temporary)
# cramer.values are the thresholds of toxicological concern for each Cramer class (temporary)

cramer <- tidy(opbase.data("Op_en6193", subset = "cramer_classes_of_compounds"))
cramer$Result <- NULL
cramer$Guidance <- "TTC" # Threshold for toxicological concern

cramer.values <- tidy(opbase.data("Op_en6193", subset = "cramer_values"))
cramer.values <- merge(cramer, cramer.values)
cramer.values[["Cramer class"]] <- NULL
cramer.values$Result <- cramer.values$Result / 1000 # from ug /kg /d to mg /kg /d

guidance.values <- tidy(opbase.data("Op_en6193", subset = "guidance_values"))
guidance.values <- orbind(cramer.values, guidance.values)
guidance.values <- Ovariable("guidance.values", data = guidance.values)

safety.factor <- Ovariable("safety.factor", ddata = "Op_en6193.safety_factors") #  [[Compound intake estimator]]

# food.risk calculates "hazard quotients" or risks relative to guidance values. Value < 1 are of minimal concern

food.risk <- Ovariable(name = "food.risk",
	dependencies = data.frame(Name = c("compound.intake", "guidance.values", "safety.factor")),
	formula = function(...) {

		# guidance.risk is intake compared with guidance values.

		out <- compound.intake / (guidance.values / safety.factor)

		removepoorguidance <- function(dat, goodg) { # Remove inferior guidance values

			# Compounds that have good guidance values (goodg).
			hasgoodg <- unique(dat[dat$Guidance == goodg , "Compound"])

			# The row contains data about a chemical with a good guidance value
			# but the data is about an inferior guidance value.

			out <- dat[!(dat$Compound %in% hasgoodg & dat$Guidance != goodg) , ]

			return(out)
		}

		out <- removepoorguidance(out@output, "ADI")
		out <- removepoorguidance(out, "TDI")
		out <- removepoorguidance(out, "BMDL10")

		test2 <- unique(out$Compound[out$Result > 0]) # remove all compounds where risk = 0
		out <- out[out$Compound %in% test2 , ]

		return(out)
	}
)

extraction <- Ovariable("extraction", ddata = "Op_en6193.extraction_techniques")
extraction@data$extractionResult <- 1

### Define formatting for the concentration graphs that will be shown.

concplot <- function(condition) {
	dat <- compound.conc@output[condition, ]
	if(nrow(dat) != 0) {
		out <- ggplot(dat, 
			aes(
				x = Compound, 
				y = compound.concResult, 
				colour = Extraction.technique
			)) + 
			geom_point(size=5) +
			theme_grey(base_size = 24) +
			theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
			labs(y = "Concentration (%)")
		return(out)
	} else {
		 warning("No data to show")
	}
}

objects.store(compound.intake, guidance.values, food.risk, extraction, concplot, safety.factor)

cat("The following objects were successfully stored:
compound.intake, guidance.values, food.risk, extraction, concplot, safety.factor\n")

</rcode>

[1] PlantLIBRA website

[2] Segen's Medical Dictionary. 2012 Farlex, Inc. [1] accessed 15th April, 2014.

[3] Wikipedia: Tolerable daily intake accessed 15th April 2014

[4] U.S.EPA: Benchmark Dose (BMD) Methodology [2] accessed 15th April, 2014.

[5] Opinion on the Use of the Threshold of Toxicological Concern (TTC) Approach for Human Safety Assessment of Chemical Substances with focus on Cosmetics and Consumer Products. European Union (2012) SCCP/1171/08. [3] accessed 15th April 2014.

[EFSA_2012-6] EFSA (2012). Scientific Opinion on Exploring options for providing advice about possible human health risks based on the concept of Threshold of Toxicological Concern (TTC). EFSA Journal 2012;10(7):2750

[van_den_Berg_2011-7] van den Berg SJPL, Restani P, Boersma MG, Delmulle L, Rietjens IMCM (2011) Levels of Genotoxic and Carcinogenic Ingredients in Plant Food Supplements and Associated Risk Assessment. Food and Nutrition Sciences, 2(9): 989-1010.

[8] PlantLIBRA website, Eurofir, accessed 15th April, 2014.

[9] BASIS (Bioactive Substances in Food Information System), Eurofir, accessed 15th April, 2014.

[roe-10] 3.0 ^3.1 M.A. Roe, S. Bell, M. Oseredczuk, T. Christensen, S. Westenbrink, H. Pakkala, K. Presser and P.M. Finglas, 2013. Updated food composition database for nutrient intake. EFSA supporting publication 2013:EN-355, 21 pp. Available online: [4] [5] [6]

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[1]

[2]

[3]

Plantlibra deliverable DWP5-6

Contents

Background and objective of the work

Compound intake estimator

Question

Answer

Compound intakes and risk estimates

Compound concentrations

Rationale

Guidance values

Intake and guidance value calculations

User interface

Concentrations

Data

Calculations

Extraction techniques

Data about guidance values

See also

References

Compound concentration data for PlantLIBRA cases

Case study results

Discussion

EFSA data collection and PFS

See also

References

Appendix 1: Guidance values and other inputs

Appendix 2: Computer code for the tool

Compound intakes and risks

Compound concentrations

Model initiation

Navigation menu