Technical advance | Open | Open Peer Review | Published:
Worldwide food recall patterns over an eleven month period: A country perspective
BMC Public Healthvolume 8, Article number: 308 (2008)
Following the World Health Organization Forum in November 2007, the Beijing Declaration recognized the importance of food safety along with the rights of all individuals to a safe and adequate diet. The aim of this study is to retrospectively analyze the patterns in food alert and recall by countries to identify the principal hazard generators and gatekeepers of food safety in the eleven months leading up to the Declaration.
The food recall data set was collected by the Laboratory of the Government Chemist (LGC, UK) over the period from January to November 2007. Statistics were computed with the focus reporting patterns by the 117 countries. The complexity of the recorded interrelations was depicted as a network constructed from structural properties contained in the data. The analysed network properties included degrees, weighted degrees, modularity and k-core decomposition. Network analyses of the reports, based on 'country making report' (detector) and 'country reported on' (transgressor), revealed that the network is organized around a dominant core.
Ten countries were reported for sixty per cent of all faulty products marketed, with the top 5 countries having received between 100 to 281 reports. Further analysis of the dominant core revealed that out of the top five transgressors three made no reports (in the order China > Turkey > Iran). The top ten detectors account for three quarters of reports with three > 300 (Italy: 406, Germany: 340, United Kingdom: 322).
Of the 117 countries studied, the vast majority of food reports are made by 10 countries, with EU countries predominating. The majority of the faulty foodstuffs originate in ten countries with four major producers making no reports. This pattern is very distant from that proposed by the Beijing Declaration which urges all countries to take responsibility for the provision of safe and adequate diets for their nationals.
The worldwide importance of a safe diet is reflected in the in depth discussions and reports leading up to and following the Beijing Declaration on Food Safety [1–11]. Following the World Health Organization Forum in November 2007, the Declaration was adopted by over 50 countries and identified the right of the individual to a safe and adequate diet along with guidance to countries for its provision. The Declaration stipulates the essential public health function of food safety controls and, moreover, the application of equal food safety measures within and between countries. It urges all countries to "establish food safety authorities ...within a comprehensive production-to-consumption legislative framework".
A great deal of discussion has occurred regarding the development of a new concerted model for food safety controls in Europe. Following decades of piecemeal regulation, the European Food Safety Authority (EFSA) was established in 2002 to provide a platform for scientific advice and the commitment required to ensure customer protection . The final format veered from the US Food and Drug Administration approach, which has a focus on the three strands of risk-assessment, -management and -communication. In contrast, the EU model has separated risk assessment from management to ensure that "the control must be at the heart of the Commission's risk management process".
In 1979, the European Commission set up the Rapid Alert System for Food and Feed (RASFF), a key instrument for customer protection . Although, it has the principal aim of alerting countries to immediate hazards, the RASFF provides a useful database for studying historical trends in food safety issues, along with the potential to predict future risks . Whatever the nature of the potential hazard, be it a chemical, microbiological, parasite, packaging or contaminant of human origin, the Beijing Declaration fully recognizes the need for improved and continued annual reporting systems. A key aspect of the Declaration is that each country should be actively engaged in the process. Numerous reports exist on the breakdown of food types, hazard types and countries of origin and reporting body [12, 13]. However, despite the Beijing Declaration, there is a paucity of reported studies of analyses of patterns from a nation perspective – both the reporting countries and countries of origin for hazardous foods.
The aim of this study is to retrospectively analyze the patterns in risk reporting within and between countries to identify the principal hazard generators and gatekeepers of food safety. The reporting period has been chosen as the eleven months prior to the Beijing Declaration with a view to a future comparative study being possible to monitor enhanced or reduced adherence to the Declaration over time.
The food recall data set was collected by the LGC (UK) over the period from January to November 2007, inclusive. The data set contained detailed information which included: Date; Agency; Alert Reference; Type of Alert; Country Notified by; Company; Reasons for recall; Category; Coded category; Metal categories; Country of Origin; Type of Control; and Status. Descriptive statistics and statistical analyses were computed by using Excel, SPSS 15.0 and the R statistical framework . Network properties such as degrees, weighted degrees, modularity and k-core decomposition were performed and graph visualisations were drawn by using the igraph library for network science . Country codes with full country names are listed in Additional file 1: (Weighted and unweighted in-degree (transgressor) and out-degree (detector) for all countries). A brief explanation and calculation of modularity is provided in Additional file 2: (Definition of modularity).
Results and discussion
With the rapid improvement in computational power, network analysis has become a useful tool to analyze complex information for underlying structure or patterns, otherwise undetectable with descriptive analyses . Information on the network size, connections, and structural properties such as the number of layers or cluster formations is capable of revealing important collective information about the system . Mathematical manipulations utilised in this paper can be equally applied to data containing hundreds (such as the present data), thousands or millions of data points (nodes) which can be visualised by graphical methods. Owing to the recently discovered similarities in structural properties among different systems, ranging from social groups through scientific collaboration, traffic, and neural networks to metabolic networks and food webs, describing these systems by their network properties has turned to be highly instrumental . The deconvolution approach to complex interrelationships using network analyses has been applied to diverse systems including i) drug-therapy interactions, ii) evolution-ecology investigations, and iii) elucidating mechanisms of disease aetiology [19–22].
The complexity of the recorded interrelations for the food recall data is best to be pictured as a network, which can be constructed from structural properties (relationships) contained in the data. Figure 1a depicts all 117 countries involved in food alerts either as detector or transgressor (or both), connected by 574 links. Two countries are linked if one of them reported on the other at least once during the examined time period where the direction of the arrow reflects the roles of the countries: arrows originate from detectors and are oriented towards transgressors. Mutual arrows are present between two countries if they both reported on each other. The weight assigned to an arrow reflects the number of reports with the width of the arrow proportional to the logarithm of its weight. Thus, the role of a country in the network is reflected by the number, direction and width of its adjacent arrows.
The average number of reports between countries (i.e. the average edge weight) was 4.28 ± 9.11 (maximum weight = 126) and only 29 (5.1%) links had weights above the mean plus one standard deviation (Figure 1b). These links will be called strong links from now on. It is notable that the distribution of the weights was positively skewed with 456 links below the mean. This suggests that the full graph is predominantly connected by relatively weak links. The weight distribution is further discussed below in the 'breakdown of reporting pattern' section. Interestingly, the largest distance in the network (the so-called diameter of the network) is only 5, even though there are a lot more countries in the dataset.
Countries connected with strong links (Fig 1b) were: TUR, ITA, GER, FRA, USA, MLT, NED, NOR, THA, NGR, IND, GBR, CHN, GRE, POL, ESP, IRJ, DEN, CZE, SVK. Of these 20 countries, the UK, ESP and GER made alerts/reports on products from their own countries (as indicated by the 'loops' in Fig 1b). For clarity, the UK is labeled as GBR. To unveil the inner structure of the network, attempts were made to decompose it to several smaller clusters (densely connected subgroups that are connected by a small number of edges). A common approach to decompose a network is to find a partition of the nodes that maximizes the measure called modularity . (For a more detailed description of modularity, refer to Additional file 2: Definition of modularity). Modularity maximization was performed by the heuristic introduced by Latapy and Pons . As a rule of thumb, Newman stated that a modularity value < 0.3 indicates the absence of a strong clustered structure . The best modularity we could achieve was only 0.214 (which is already likely to be enhanced by the presence of loop edges, i.e. countries reporting on themselves), which suggests that the network is organized around a dominant 'core' surrounded by nodes that connect to the core only via weak links. Assuming that important countries are well connected and globally centred (as supported by our cluster analysis), further analyses were conducted to identify these countries based on their centrality, using k-core decomposition .
The k-core analysis is a step by step process where the least connected nodes are gradually removed from the graph. The whole network is called the 0-core. Removing isolated nodes (having no links) result in the so-called 1-core. The process continues by removing the nodes having only a single link (thus their degree is 1). This may have created other nodes that are left with only a single link, so the process is repeated until there are no more nodes with degree 1. The remaining network is called the 2-core, since the degree of each node is at least 2. Removing all nodes having degree 2 yields the 3-core and so on. The process ends when we removed all nodes. The coreness or shell index of a node is the largest k for which the node is the member of the k-core but not the member of the k+1-core. In summary, this iterative process results in a series of subgraphs that gradually reveal the globally central region of the original network and stops when no more layers can be peeled off.
The k-core analysis revealed 12 layers (i.e. no countries could be eliminated so the remaining ones would have 13 connections). The innermost subgraph (the 12-core) is depicted in Figure 2.
This 12-core is the largest subgraph where each node has at least 12 interconnections. The average edge weight within the subgraph is 7.35 ± 13.98 (max = 126.00). This subgraph consists of 21 countries with 190 connections, with 16 connections having a weight larger than mean + SD. These are: USA → USA (126), GER → TUR (73), ITA → ESP (68), ITA → CHN (58), GBR → GBR (51), GRE → CHN (45), GER → CHN (36), NED → USA (31), DEN → GER (30), ESP → ESP (28), GER → GER (27), ESP → CHN (27), ITA → TUR (27), DEN → FRA (26), GBR → CHN (23), FRA → TUR (22).
In practical terms, countries in this 12-core subgraphs are the ones that have played an important role in food safety during the time period when the data were collected (January – November 2007, inclusive). However, this 'importance' may come from two very distinct sources: detector or transgressor which can be accounted for in the k-core decomposition. When considering this additional information, the 117 countries form two distinct 5-layer graphs. The transgressor 'layer 5' contains all countries on which at least 5 reports were made among each other whereas 'layer 0' shows countries on which no reports were made during the time period of this study. The detector 'layer 5' contains all countries who issued at least 5 reports on other countries in the same layer whereas 'layer 0' is formed by countries that made no reports (but intuitively, they are the ones that had been reported on, otherwise they would not have been included). Countries by layers are summarized in Table 1.
The 5th layer for both detectors and transgressors are depicted separately in Fig 3, with a notable size difference. Only 11 European countries appear to play an important role in policing industrial food safety, mostly from the European Union (EU). Among the transgressors, only 36% were European and 30% from the EU.
The number of transgressors (of the 5th layer) appeared to be independent of the population (hence remotely to the market size). Population statistics were obtained from the CIA's 2008 World Factbook . The rank correlation between the number of reports received and population was small and non-significant (Kendall tau = .168, p = .2025).
Interestingly, there was a noticeable overlap between the transgressor and detector graphs for the UK, France, Italy, Germany and Belgium. These countries seem to guard 'food safety' and report all problems/hazards, regardless of the origin. Not surprisingly, three of these countries (UK, France and Germany) not only reported faulty food products from other countries but also made a considerable number of reports on their own products (Fig 1).
Owing to the noticeable difference between the set of transgressors and detectors, this aspect was further analyzed focusing exclusively on countries that are most involved in testing and those with a defective food product. All data are provided in Additional file 1: (Weighted and unweighted in-degree (transgressor) and out-degree (detector) for all countries).
Breakdown of reporting patterns
An overview of reporting patterns demonstrates that the reporting behaviour is heavily skewed towards a small number of countries, both in terms of reporting and those being extensively reported (Table 1). Figure 4 shows that the number of reports logged or received by the 117 counties decreases rapidly. The highest number of reports made during the eleven-month period (406) exceeded the highest number of reports received (281) considerably.
Reports were made against products produced by the 115 countries, with the top 5 countries having received between 100 to 281 reports in the 11 month period. The next twenty countries received between 25 to 100 reports against their products in this period. The numbers further decline steadily from 23 to 1 report being made for products from the next 88 countries, with no reports against products from Israel and Luxembourg. Ten countries that are most frequent transgressors provided sixty per cent of all faulty products marketed.
In contrast, the pattern for the 30 detectors is that the top three make over 300 reports each (Italy 406, Germany 340, UK 322) as detailed in Additional file 1: (Weighted and unweighted in-degree (transgressor) and out-degree (detector) for all countries). The next five countries make between 99 and 142 reports with the final 20 countries making between 71 and 1 reports. In all, the top ten detectors account for three quarters of reports.
Twelve countries constitute the most active level with the reporting activities highlighted in Figure 5a for numbers of transgressors. Relationships are given for both detectors (in red) and transgressors (in blue). Considerable variations in reports occur when both the quantities and countries are examined (Fig 5). Large variations exist for the top twelve countries between those reporting and those reported. A lesser degree of variation was observed whether reports relate to countries or products (i.e. Fig 5a versus 5c).
Reporting rates varied enormously between countries with no clear relationship with population size. The twelve-country segment selected with optimal rates of reporting differs from those for reported activities which may reflect differences in market size. EU countries predominate with each of Italy, Germany and the UK reporting products from over 50 countries (Fig 5a). This pattern remains when the total number of reports is considered (Fig 5c). However, the pattern of detectors changes when total reports are considered as the USA enters at number 4 signifying it makes a lot of reports from against fewer countries. In fact the USA made 126 reports against its own products during the period which may reflect the market size.
As expected for this segment, countries with larger populations often predominate, reflecting market size. However, this generalization frequently breaks down when the data are analyzed for individual countries as detectors or as transgressors. The order of reports against a country's products is given in Figure 5b in terms of numbers of transgressors. Of the core twelve countries, four are present owing only to reports against their products (China, Turkey, India and Brazil) and they have no reports against products arising from other countries. In contrast, a relatively large number of products from other countries have been reported by Italy and Germany, with a significantly reduced number of countries reporting them. This trend is followed by four countries, albeit to a lesser degree (Spain, Netherlands, France and Belgium).
In terms of total numbers of reports against a country, China leads, closely followed by the USA and Turkey. The pattern continues with countries from the EU in the order: Spain, Ireland, Germany, France and the UK.
On consideration of the total number of reports, the pattern exhibited for numbers of transgressors is reflected in the data for total number of products reported (Figure 5b and 5d). The same four countries have no reports of products from other countries with China and Turkey having a large number of reports against their products (> 200). Italy and Germany remain at the top in terms of foreign product reporting but with few reports against their products by other countries. Similarly, there are a large number of reports against products from the USA, which has in turn reported many products.
Limitations and extensions
The analyses were limited to approximately one year period. By expanding the timeframe of the analysis, trends and effects of various legislative and regulatory changes, such as the Beijing Declaration, can be investigated. As the reasons for food alerts/recalls vary greatly from illegal import (1.1%) through unauthorized food additives or ingredient (11.2%) to mycotoxin contamination (23.7%), future research would benefit from focused analyses. The methodology used for this paper can also be applied to various groups of countries based on the reasons for food alerts/recalls. Particular interests are reports owing to i) processing, ii) microbiological or iii) chemical contaminations.
The vast majority of food reports are made by 10 countries out of the 117 studied, with EU countries predominating. The majority of the faulty foodstuffs originate in ten countries with four major food producers making no reports at all. A substantial worldwide improvement in food testing and reporting is required in response to the Beijing Declaration.
World Health Organization. [http://www.who.int/foodsafety/fs_management/meetings/forum07/en/index.html]
Kok EJ, Keijer J, Meter GA, Kuiper HA: Comparative safety assessment of plant-derived foods. Regul Toxicol Pharmacol. 2008, 50 (1): 98-113.
Hugas M, Tsigarida E, Robinson T, Calistri P: Risk assessment of biological hazards in the European Union. Int J Food Microbiol. 2007, 120 (1–2): 131-135.
Sergent T, Ribonnet L, Kolosova A, Garsou S, Schaut A, De Saeger S, Van Peteghem C, Larondelle Y, Pussemier L, Schneider Y-J: Molecular and cellular effects of food contaminants and secondary plant components and their plausible interactions at the intestinal level. Food Chem Toxicol. 2008, 46 (3): 813-841.
Birchard K: Europe tackles consumers fears over food safety. Lancet. 2001, 357 (9264): 1274-
Millstone E, Lang T, Naska A, Eames M, Barling D, van Zwanenberg P, Trichopoulou A: 'European Policy on Food Safety': Comments and suggestions on the White Paper on Food Safety. Trends Food Sci Tech. 2000, 11 (12): 458-466.
Bergeaud-Blackler F, Paola Ferretti M: More politics, stronger consumers? A new division of responsibility for food in the European Union. Appetite. 2006, 47 (2): 134-142.
Coppens P, Fernandes da Silva M, Pettman S: European regulations on nutraceuticals, dietary supplements and functional foods: A framework based on safety. Toxicology. 2006, 221 (1): 59-74.
Paulsen P, Luf W, Smulders FJM: Different legislations on toxicants in foodstuffs. Food Toxicants Analysis. Edited by: Pico Y. 2007, Amsterdam: Elsevier, 11-31.
Halkier B, Holm L: hifting responsibilities for food safety in Europe: An introduction. Appetite. 2006, 47 (2): 127-33.
Marvin HJP, Kleter GA, Prandini A, Dekkers S, Bolton DJ: Early identification systems for emerging foodborne hazards. Food Chem Toxicol. 2007
Rapid Alert System for Food and Feed (RASFF). [http://ec.europa.eu/food/food/rapidalert/index_en.htm]
Kleter GA, Prandini A, Filippi L, Marvin HJP: Identification of potentially emerging food safety issues by analysis of reports published by the European Community's Rapid Alert System for Food and Feed (RASFF) during a four-year period. Food Chem Toxicol. 2007
R Development Core Team: R: A language and environment for statistical computing. 2006, R Foundation for Statistical Computing, Vienna, Austria, ISBN 3-900051-07-0, [http://www.r-project.org]
Csárdi G, Nepusz T: The igraph software package for complex network research. InterJournal Complex Systems. 2006, 1695-
Newman MEJ: The structure and function of complex networks. SIAM Review. 2003, 45: 167-256.
Palumbo MC, Farina L, Colosimo A, Tun K, Dhar PK, Giuliani A: Networks everywhere? Some general implications of an emergent metaphor. Current Bioinformatics. 2006, 1: 219-234.
Watts DJ, Strogatz SH: Collective dynamics of "small world" networks. Nature. 1998, 393: 440-442.
Ma'ayan A, Jenkins SL, Goldfarb J, Iyengar R: Network analysis of FDA approved drugs and their targets. Mt Sinai J Med. 2007, 74 (1): 27-32.
Nacher JC, Schwartz JM: A global view of drug-therapy interactions. BMC Pharmacol. 2008, 8: 5-
Prolx SR, Promislow DEL, Phillips PC: Network thinking in ecology and evolution. Trends Ecol Evol. 2005, 20 (6): 345-353.
Kleemann R, Verschuren L, van Erk MJ, Nikolsky Y, Cnubben NHP, Verheij ER, Smilde AK, Hendriks HFJ, Zadelaar S, Smith GJ, Kaznacheev V, Nikolskaya T, Melnikov A, Hurt-Camejo E, Greef van der J, van Ommen B, Kooistra T: Atherosclerosis and liver inflammation induced by increased dietary cholesterol intake: a combined transcriptomics and metabolomics analysis. Genome Biol. 2007, 8: R200-
Newman MEJ: Fast algorithm for detecting community structure in networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2004, 69 (Pt 6 - 2): 066133-
Latapy M, Pons P: Computing communities in large networks using random walks. J Graph Algorithms Appl. 2006, 10 (2): 191-218.
Seidman SB: Network structure and minimum degree. Soc Networks. 1983, 5: 269-287.
Central Intelligence Agency:2008 The World Factbook. ISSN 1553-8133. [https://www.cia.gov/library/publications/the-world-factbook/index.html]
Newman MEJ, Girvan M: Finding and evaluating community structure in networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2004, 69 (2 Pt 2): 026113-
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/8/308/prepub
The food recall data set was collected by LGC (UK) over the period from January to November 2007, supported under contract with the UK Department for Innovation, Universities & Skills as part of the Government Chemist programme (see http://www.governmentchemist.org.uk).
The authors declare that they have no competing interests.
TN conducted the network analyses, prepared the network visualisations and assisted in writing the manuscript. AP prepared the data for network analyses, conducted statistical analyses and assisted in writing the manuscript. DN suggested the study, sought the data, assisted with the data interpretation and contributed to writing the manuscript. All authors have read and approved the final version of the manuscript.