Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir
© Codeço; licensee BioMed Central Ltd. 2001
Received: 13 July 2000
Accepted: 2 February 2001
Published: 2 February 2001
In the last decades, attention to cholera epidemiology increased, as cholera epidemics became a worldwide health problem. Detailed investigation of V. cholerae interactions with its host and with other organisms in the environment suggests that cholera dynamics is much more complex than previously thought. Here, I formulate a mathematical model of cholera epidemiology that incorporates an environmental reservoir of V. cholerae. The objective is to explore the role of the aquatic reservoir on the persistence of endemic cholera as well as to define minimum conditions for the development of epidemic and endemic cholera.
The reproduction rate of cholera in a community is defined by the product of social and environmental factors. The importance of the aquatic reservoir depends on the sanitary conditions of the community. Seasonal variations of contact rates force a cyclical pattern of cholera outbreaks, as observed in some cholera-endemic communities.
Further development on cholera modeling requires a better understanding of V. cholerae ecology and epidemiology. We need estimates of the prevalence of V. cholerae infection in endemic populations as well as a better description of the relationship between dose and virulence.
Since Snow's seminal work (1855), when cholera epidemics were first associated with a contaminated water supply, this disease has been the focus of intense study. Cholera is a waterborne disease characterized by severe diarrhea. The etiological agent is V. cholerae O1 (more recently also V. cholerae O139), which colonizes the small intestine and produces an enterotoxin responsible for a watery diarrhea. Without prompt treatment, a person with cholera may die of dehydration in a matter of hours after infection. Cholera outbreaks are generally associated to contaminated food and water supplies. Appropriate sanitation and safe water are the main weapons against this disease.
Until the 19th century, cholera was confined to the Indian sub-continent . There, cholera outbreaks are seasonal with one or two peaks per year . From this region, cholera has spread throughout the world seven times since 1817 . The last pandemics began in 1961 in Indonesia, spread through the Asian continent during the 60's, reached Africa in 1970  and Latin America in 1991. The occurrence of successive cholera outbreaks throughout Africa and Latin America during the 90's raised the concern that cholera had established itself in these regions as an endemic disease.
Once cholera arrives into a new region, either carried by an infected person or by contaminated water and food, we may expect one of three possible outcomes: no outbreak, an outbreak possibly followed by few waves; or a cholera outbreak followed by subsequent outbreaks that may assume a persistent seasonal pattern. I call these, respectively, cholera-free; epidemic and endemic patterns. What demographic, environmental and sociological factors drive these outcomes? Studies suggest that the number of susceptibles; exposure to untreated water and sewage; and the presence of an aquatic reservoir of V. cholerae are important factors   . How important they are and how they interact with each other and with other variables to drive cholera epidemiology, are open questions.
Recently, it was found that toxigenic V. cholerae can survive in some aquatic environments for months to years, in association with zooplankton and other aquatic organisms . Under stress, this pathogen assumes a viable but not culturable state, becoming undetectable to traditional bacteriological techniques. These findings suggest that the aquatic environment may be a reservoir of toxigenic V. cholerae in endemic regions. What is the role of the aquatic reservoir in promoting epidemic and endemic cholera? What regions are more likely to maintain endemic cholera? What are the best approaches to prevent and control cholera outbreaks? What are the best predictors of the fate of a community after the introduction of cholera?
Here, I start to address these questions by proposing a mathematical model that incorporates what are considered important determinants of cholera dynamics. Despite the great effort towards understanding cholera dynamics, the mechanisms driving V. cholerae dynamics in water is still poorly known. Thus, it is appropriate to start a modeling exercise with a simple model (with few variables) and then add complexity as needed. Following this methodology, I started with a basic Susceptible-Infective-Recovered model coupled to an aquatic population of V. cholerae. The objective is to explore the role of the aquatic reservoir on the persistence of endemic cholera as well as to define minimum conditions for the development of epidemic and endemic cholera. Mathematical models are powerful tools for the development of epidemiological theories . They can synthesize the current empirical knowledge about the disease into a coherent mechanistic framework. These models may help us to infer causal relationships and to suggest experimental designs to test alternative hypotheses.
This work is structured in three parts. First, the mathematical model for epidemic and endemic cholera is presented. Secondly, the proposed model is applied to three hypothetical communities to simulate cholera-free, epidemic and endemic situations. Thirdly, I consider possible causes of endemic oscillations. I finish this work discussing insights brought from model analyses into cholera dynamics.
The model proposed here is an extension of Capasso's model , used to describe the 1973's cholera epidemics in Italy. In Capasso's version, two equations describe the dynamics of infected people in the community and the dynamics of the aquatic population of pathogenic bacteria. In our formulation, the dynamics of the susceptible population is included since I wish to study long term dynamics . The mathematical model is:
Symbols used in the model
number of susceptibles
number of infected
concentration of toxigenic V. cholerae in water (cells/ml)
total human population
Human birth and death rates (day-1)
rate of exposure to contaminated water (day-1)
concentration of V. cholerae in water that yields 50% chance of
catching cholera (cells/ml)
rate at which people recover from cholera (day-1)
growth rate of V. cholerae in the aquatic environment (day-1)
loss rate of V. cholerae in the aquatic environment (day-1)
contribution of each infected person to the population of V. cholerae
in the aquatic environment (cell/ml day-1 person-1)
where K is the concentration of V. cholerae in water that yields 50% chance of catching cholera. I assume that the only route for infection is the ingestion of contaminated water from non-treated sources.
Equation 1b describes the dynamics of infected people in the community. This category includes not only cholera cases but also those with asymptomatic and mild infections. In reality, only 1 to 30% of V. cholerae infections actually develops into severe cholera cases . By combining all infection types into a single model compartment, I am assuming that they follow the same dynamics, i.e., that the case-to-infection ratio remains constant through the epidemics. This is reasonable if infection virulence is strongly determined by host factors (sensitivity to cholera toxin, blood type, etc) and bacterial factors (biotype, etc), that are not expected to change during the period of interest   . A potential problem with this assumption is that inoculum size may affect the case-to-infection ratio . It is not clear, however, if inoculum size is directly related to virulence or if it just increases the chance of intestine colonization. Here, I assumed that inoculum size affects the per capita infection rate (i.e., the probability of colonization), but not the severity of symptoms. Equation 1b states that the infected population increases as susceptibles become infected (first term in the equation) and decreases as they recover from the disease or die.
The third equation describes the dynamics of pathogenic V. cholerae in the aquatic reservoir, in this case, the set of untreated waters consumed by the population. Environmental V. cholerae is found in ponds, wells, rivers, estuaries and coast waters. These environments have quite distinct physical, chemical and biological characteristics and V. cholerae dynamics is probably controlled by different factors in each of them. Classical studies on the survival of culturable V. cholerae in aquatic environments suggest that V. cholerae cannot maintain a stable population in the environment. Observed rates of population extinction in microcosms are variable, ranging from 0.02 day-1 to >3 day-1 . Islam , however, showed that population decay does not necessarily imply death, but also the transition towards a non-culturable state. Using fluorescent antibody techniques, he found V. cholerae surviving for more than 15 months within the mucilagenous sheaths of a filamentous alga (extinction rate <0.002). Here, equation 1c states that bacterial density in the water results from the balance between local birth and death processes and the inflow of contaminated sewage. The parameter e defines the average contribution of each infected person to the aquatic population of V. cholerae.
At last, equations (2) specify the initial conditions (all individuals are initially susceptible).
Parameters for three hypothetical communities.
Expected cholera dynamics
1 × 106
1 × 106
1 × 106
First case: cholera-free population
Consider a community that does not experience cholera for generations. All individuals are susceptibles. There are neither infective or immune individuals nor toxigenic bacteria in the water.
The S c threshold increases proportionally to the barriers to cholera infection (K), and recovery (r) as well as the net mortality rate of V. cholerae in the water. It decreases, on the other hand, as contamination of water supplies as well as contact with these waters increase. In other words, the better water quality and sewage treatment are, the greater must be the pool of susceptibles in order to trigger a cholera outbreak.
Rearranging equation 5, we find the maximum degree of contamination each infected person may cause to the water reservoir without causing public health hazards:
Asymptomatic cases of V. cholerae infection yield 102-105 V. cholerae per gram of feces . In severe cases, yield increases to 106-109 cells/ml of rice-water stool . If we know the proportion of infections that are asymptomatic, mild and severe, we could obtain a rough estimation of e as
Consider a hypothetical community with 10,000 members sharing the same water source. Suppose that the water reservoir is directly contaminated with sewage. If a = 1, cholera outbreak in this community will develop only if individual contribution to water contamination exceeds ca. 7 cells/ml day-1. Now, consider that each infected person yields 1 × 104 cells person-1 day-1. This water will remain infective unless it is diluted to a value less than 7 cells/ml. Such dilution requires 4 × 104/7 = ca. 5700 liters of water per infected person per day.
A real pond, however, is not perfectly mixed as this model assumes. If water used for consumption is taken from a point close to the excretion site, the risk of acquiring the infection will be higher even if the pond is large. Moreover, V. cholerae is found associated to phytoplankton, macrophyte, zooplankton, crustacea and other aquatic organisms. On the surface of these organisms, density of bacteria may be 100 to 1000 times greater than in the aquatic medium. These organisms are not evenly distributed within the water body. Phytoplankton tends to concentrate on the water surface, zooplankton migrates along the water column on a daily basis. If water is taken from phyto- or zooplankton-rich patches, risk of catching cholera will increase, enhancing the probability of triggering an outbreak in the community.
This model is too simple to provide quantitative predictions on cholera dynamics. Nevertheless, its qualitative results verifies known alternative approaches to the prevention of cholera outbreaks:
1. Minimize water contamination as well as consumption of untreated waters: This is the classical approach to cholera control. Good sanitation reduces the parameter e and water treatment reduces parameter a. The smaller these parameters are, the larger must be the susceptible pool in order to a cholera outbreak to develop (equation 5).
2. Dilute cholera diarrhea with large amounts of water to make water uninfective.
Second case: Epidemic cholera
Third case: Endemic cholera
Consider our hypothetical community 3 (Table 2). Its parameters are identical to those set for community 2, except that n is greater. Since S 0 >S C, the introduction of infectives in community 3 starts a cholera outbreak. In this case, however, cholera does not vanish after the first peak. It actually returns in successive waves and eventually converges to a positive endemic equilibrium (Figure 5). Setting the derivatives of equation system 1 to zero and solving it algebraically, we obtain the endemic equilibrium:
The equilibrium number of infected individuals in the population is a fraction φ of the surplus population (H-S C), where φ is
Fraction φ increases as the susceptible turnover rate n increases and tends to a at large values of n. This fraction is a function of human parameters, only.
Although in theory any community with S 0 >S C and n > 0 is expected to become endemic, this is not likely to occur in practice. To maintain an endemic state, I* must be much greater than 1, otherwise infection will fade-out due to stochastic processes .
Stability analysis (see additional material: Appendix) indicate that the endemic equilibrium is stable to small perturbations if S C >S 0 and mb >nb. The route to the endemic equilibrium involves dumped oscillations (Figure 5). Cholera waves arise when the susceptible population grows above the threshold S C, triggering a new outbreak. Eventually, the number of infectives stabilizes into a fixed fraction of the population. Bacterial dynamics in the water is also oscillatory and follows the human excretion pattern.
In summary, our model predicts that cholera outbreaks should occur in communities whose susceptible pool is greater than the susceptible threshold (S C). The magnitude of the threshold depends on a combination of environmental, sociological and strain-specific factors. The environmental reservoir reduces the susceptible threshold. In the presence of a permanent reservoir of V. cholerae (nb tending to mb), the threshold tend to zero, and community of any size in contact with the reservoir should be subject to an outbreak.
Environmental triggers of cholera outbreaks
In endemic regions, cholera outbreaks are often associated to climatic events . In some parts of Africa, outbreaks occur during the dry season or right after heavy rainfalls . In the Americas, disasters caused by the El Niño preceded large cholera outbreaks . In Bangladesh, cholera season coincides with the post-monsoon period . Flooding and drought are likely to affect cholera dynamics in a complex way. Flooding of streets and cities washes contaminated feces and sewage into the rivers. It may also disrupt water distribution service and aggravate hygiene conditions . Drought, on the other hand, reduces the availability of potable water; aggravates hygiene conditions; increases the number of people sharing the same water supply; and may increase per capita water contamination. All these factors contribute to a greater rate of contact with unsafe waters (parameter a) and increased (floods) or reduced (floods or droughts) dilution of the per capita bacterial yield (parameter e).
Besides weather changes, seasonal variations of water temperature has been also associated to cholera outbreaks. Cholera outbreaks in Peru  and Bangladesh  are preceded by periods of warm waters, during which toxigenic V. cholerae is detected in the environment. Colwell  explains this pattern by an association between zooplankton seasonal growth (that generally peaks during the warm season) and V. cholerae dynamics.
Floods, droughts, temperature and biotic interactions have potentially quite different impacts on cholera seasonality. To untangle the potential roles of these factors, I conducted a series of simulations of equation system 1 where the following parameters are assumed to vary seasonally: 1) the contact rate with contaminated water (simulating seasonal variations in water quality); 2) the per capita excretion rate (simulating seasonal variations in per capita water contamination); 3) the growth rate of V. cholerae (simulating seasonal variations in autochtonous V. cholerae abundance in the water). In all simulations, a sine function with period equal to 365 days was used to simulate the seasonal oscillations.
Scenario 1. Seasonal oscillation of the contact rate (a)
Here, I substituted the parameter a in equation 1 by the sin function:
Parameter values used in the sensitivity analysis.
1,000; 20,000; 40,000; 60,000; 80,000; 100,000
0.0001; 0.002; 0.004; 0.006; 0.008; 0.01
0.1; 0.2; 0.4;.0.6; 0.8; 1
1,000; 20,000; 40,000; 60,000; 80,000; 100,000
0.02; 0.05; 0.09; 0.13; 0.16; 0.2
-0.02; -0.13; -0.25; -0.37; -0.48; -0.6
1; 20; 40; 60; 80; 100
Scenario 2. Seasonal oscillation of the per capita contamination rate (e)
In this scenario, the parameter e in equation system 1 was replaced by the periodic function:
A prevalence "plateau" often follows the annual cholera outbreak. The dynamics is very similar to that observed for the scenario 1. For some parameter combinations, however, this plateau rises to form a second outbreak of minor intensity. Conditions favoring the occurrence of this second peak include high K, high extinction rate and fast recover rate.
Scenario 3. Seasonal oscillation of V. cholerae extinction rate (nb-mb)
In this last scenario, the term netnb = (nb-mb) in equation system 1 was replaced by the periodic function:
The goal is to simulate, in a phenomenological way, the seasonal variation of V. cholerae abundance in the environment due to autochtonous growth. Seasonal growth of V. cholerae in the environment has been linked to factors as temperature, copepod abundance and chlorophyll. However, the mechanistic relationships between these variables are still not clearly defined and I opted for not including them explicitly in the model.
In the last decades, attention to cholera epidemiology increased, as cholera epidemics became a worldwide health problem. Detailed investigation of V. cholerae interactions with its host and with other organisms in the environment suggests that cholera dynamics is much more complex than it was previously thought. We know now that cholera toxin is codified by a lisogenic phage that can be transferred horizontally between bacteria . We also know that V. cholerae can survive in the water for a long time, assuming a non-culturable form, living in association with aquatic organisms. These findings increase the complexity of the cholera system, that passes from a two-factor system (human - pathogen) to a three or four factor system (virus-bacteria-aquatic organisms-humans?). Despite the great effort towards understanding cholera dynamics, we are still scratching the surface of this problem. Most studies associating cholera dynamics to environmental events are based on correlations. The causal mechanisms behind these correlations are still poorly understood. Unfortunately, many experimental studies performed in the past, on the tolerance and survival of V. cholerae in the environment (important for the establishment of causal relationships), may be misleading because they used bacteriological techniques for detecting V. cholerae. These findings may require further considerations on the light of the existence of unculturable forms of V. cholerae.
Experimental effort towards understanding cholera epidemiology has been accompanied by a theoretical effort. Authors have proposed new conceptual models of cholera epidemiology where the aquatic environment has an important role on the emergence of new strains of toxigenic V. cholerae  and on the seasonality of cholera outbreaks  .
The model introduced here is based on a previous cholera model , modified to incorporate an environmental reservoir. We still know very little about the dynamics of V. cholerae in the environment. Consequently, I opted for modeling the population dynamics of V. cholerae in the water in the simplest way possible, i.e., with density independent growth and death rates. This model is quite abstract but provides some insights about the role of an environmental reservoir on cholera epidemiology.
According to the model, a permanent environmental reservoir of toxigenic V. cholerae (mb = nb) should increase R0 towards infinity (equation 8). It means that any susceptible human population exposed to this water would be prone to cholera outbreaks. The magnitude of the problem, in terms of public health, would vary according to the probability of secondary transmission. In a community with good sanitation, cholera would show up as occasional primary cases without secondary transmission. In poor communities, on the other hand, endemism would result from the close contact between susceptibles and the source of contamination. An infinite R0 is an overestimation of this simplistic deterministic model. The probability of caching a V. cholerae infection is dose-dependent. It is likely a minimal reservoir should be required for the occurrence of cholera cases in the population.
Endemic cholera, however, can be maintained even in the absence of a permanent reservoir (mb >nb). This is a less obvious result. Endemism is maintained if bacteria net loss rate (mb-nb) is sufficiently low to maintain an aquatic population of V. cholerae until the susceptible pool crosses the threshold S C again. However, the best the sanitation condition (low a and e), the longer-lived the reservoir must be. On the other hand, the worse the sanitation (high a and e), the shorter can be the bacterial residence time in the water. In the limit, we may say that endemism in sanitized communities requires a permanent reservoir while endemism in poor communities requires just transient reservoirs (and a sufficiently high turnover of susceptibles).
This result suggests that cholera endemism in the US Golf Coast, for example, may be due to a permanent reservoir of V. cholerae. In some African regions, on the other hand, poor conditions may maintain endemism even in the absence of an environmental reservoir  .
This work is an effort towards the formulation of a cholera theory that incorporates a V. cholerae environmental reservoir. The model proposed is very simple and does not include many features of this complex system. Nonetheless, this study brings some new insights into cholera epidemiology. Most studies on cholera epidemiology concentrate on either social or environmental factors. This work, however, shows that the reproduction rate of cholera is a function of social and environmental factors. It is necessary to determine the relative weights of each one of these components in order to develop appropriate control strategies.
Further developments on cholera modeling require a better understanding of V. cholerae ecology and epidemiology. We need estimates of the prevalence of V. cholerae infection in endemic populations. We also need better estimates of the required infection dose as well as a better description of the relationship between dose and virulence.
I would like to thank Dr. Bradley Sack and Dr. Russek-Cohen for providing critiques of this manuscript. Support was provided by FAPERJ, FINEP and PRONEX.
- Epstein PR: Algal blooms in the spread and persistence of cholera. Biosystems. 1993, 31: 209-221. 10.1016/0303-2647(93)90050-M.View ArticlePubMedGoogle Scholar
- Islam MS, Miah MA, Hasan MK, Sack RB, Albert MJ: Detection of non-culturable Vibrio cholerae O1 associated with a cyanobacterium from an aquatic environment in Bangladesh. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1994, 88: 298-299.View ArticlePubMedGoogle Scholar
- Faruque SM, Albert MJ, Mekalanos JJ: Epidemiology, genetics and ecology of toxigenic Vibrio cholerae. Microbiology and Molecular Biology Reviews. 1998, 62: 1301-1314.PubMedPubMed CentralGoogle Scholar
- Morillon M, De Pina JJ, Husser JA, Baundet JM, Bertherat E, Martet G: Djibouti, histoire de deux épidémies de choléra: 1993-1994. Bull Soc Path Ex. 1998, 91: 407-411.Google Scholar
- Islam MS, Drasar B, Bradley SR: Probable role of blue-green algae in maintaining endemicity and seasonality of cholera in Bangladesh: a hypothesis. J Diarrhoeal Dis Res. 1994, 12: 245-256.PubMedGoogle Scholar
- Colwell RR, Huq A: Environmental reservoir of Vibrio cholerae, the causative agent of cholera. Annals of the New York Academy Society. 1994, 740: 44-54.View ArticleGoogle Scholar
- Anderson RM, May RM: Infectious diseases of humans. Oxford: Oxford University Press;. 1991Google Scholar
- Capasso V, Paveri-Fontana SL: A mathematical model for the 1973 cholera epidemic in the european mediterranean region. Rev Epidém et Santé Pub. 1979, 27: 121-132.Google Scholar
- Sengupta TK, Nandy RK, Mukhopadyay S, Hall RH, Sathyamoorthy V, Ghose AC: Characterization of a 20-k Da pilus protein expressed by a diarrheogenic strain of non-O1/non-O139 Vibrio cholerae. FEMS Microbiol. Letters. 1998, 160: 183-189. 10.1016/S0378-1097(98)00012-3.View ArticleGoogle Scholar
- Sack DA, Cadoz M: Cholera vaccines. In Plotkin AS, Orenstein WA. Vaccines. Philadelphia: WB Saunders Co;. 1999Google Scholar
- Sack FH, Pierce NF, Hennessey KN: Gastric acid in cholera and non-cholera diarrhea. Bulletin of the World Heath Organization. 1972, 47: 31-36.Google Scholar
- Clemens JD, Sack DA, Harris JR: ABO blood groups and cholera: new observations on specificity of risk and modification of vaccine efficacy. The Journal of Infectious Diseases. 1989, 159: 770-773.View ArticlePubMedGoogle Scholar
- Glass RI, Holmgren J, Haley CE: Predisposition for cholera of individuals with O blood group. Possible evolutionary significance. Am J Epidemiol. 1985, 121: 791-796.PubMedGoogle Scholar
- Hornick RB, Music SI, Wenzel RP, Cash RA, Libonati JP, Snyder MJ, et al: The Broad Street pump revisited: response of volunteers to ingested cholera vibrios. Bulletin of the New York Academy of Medicine. 1971, 47: 1181-1191.PubMedPubMed CentralGoogle Scholar
- Feachem R, Bradley DJ, Garelick H, Mara DD: Vibrio cholerae and cholera. In Sanitation and disease. Health aspects of excreta and wastewater menagement. John Wiley & Sons;. 1983, 297-325.Google Scholar
- Islam MS, Drasar B, Bradley DJ: Survival of toxigenic Vibrio cholerae 01 on a duckweed, Lemna minor. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1990, 84: 422-424.View ArticlePubMedGoogle Scholar
- Lobitz B, Beck L, Huq B, Wood B, Fuchs G, Faruque SM, et al: Climate and infectious diseases: use of remote sensing for detection of Vibrio cholerae by indirect measurement. Proceedings of the National Academic Society USA. 2000, 97: 1438-1443. 10.1073/pnas.97.4.1438.View ArticleGoogle Scholar
- Lawoyin TO, Ogunbodede NA, Olumide EAA, Onadeko MO: Outbreak of cholera in Ibadan, Nigeria. European Journal of Epidemiology. 1999, 15: 367-370. 10.1023/A:1007547117763.View ArticlePubMedGoogle Scholar
- WHO: Cholera in 1998. Weekly Epidemiological Recort. 1999, 74: 257-263.Google Scholar
- Franco AA, Fix AD, Prada A, Paredes E, Palomino JC, Wright AC, Johnson JÁ, McCarter R, Guerra H, Morris JG: Cholera in Lima, Peru, correlates with prior isolation of Vibrio cholerae from the environment. Am J Epidemiol. 1997, 148: 1067-1075.View ArticleGoogle Scholar
- Colwell R: Global climate and infectious disease: the cholera paradigm. Science. 1996, 274: 2025-10.1126/science.274.5295.2025.View ArticlePubMedGoogle Scholar
- Levin BR, Tauxe R: Cholera: nice bacteria and bad viruses. Current Biology. 1996, 6: 1389-1391.View ArticlePubMedGoogle Scholar
- Taylor RK: Virus on virus infects bacterium. Nature. 1999, 399: 312-313. 10.1038/20565.View ArticlePubMedGoogle Scholar
- Utsalo SJ, Eko FO, Umoh F, Asindi AA: Faecal excretion of Vibrio cholerae during convalescence of cholera patients in Calabar, Nigeria. European Journal of Epidemiology. 1999, 15: 379-381. 10.1023/A:1007541317104.View ArticlePubMedGoogle Scholar
- Walker I: Amazonian streams and small rivers. In Limnology in Brazil. Edited by Tundisi JG, Bicudo CEM, Tundisi TM. Rio de Janeiro: ABC/SBL;. 1995, 167-194.Google Scholar
- Nisbet RM, Gurney WSC: Modelling fluctuating populations. New York: John Wiley & Sons;. 1982Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/content/backmatter/1471-2334-1-1-b1.pdf
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.