Abstract
Mutations are normally classified according to their proximal effect on an organism’s fitness, whether beneficial, deleterious, or neutral. While this is a very useful first-pass categorization of mutations, the realization that mutations are not always haphazard, but in fact may be part of a regulated design, means that creationists should be looking for a deeper classification of mutations based on whether or not they conform to their organism’s design. Design-consistent mutations are those which occur within the pattern expected by the genome’s architecture, and design-inconsistent mutations are those which occur outside of the genome’s architecture. Features such as metabolic consistency, mutational mechanism, mutation rate, reversibility, and preservation of genome semantics can be used by biologists to assess whether or not a mutation is design-consistent or design-inconsistent.
Keywords: mutation, design, molecular biology, genome, genome architecture, simple sequence repeats, recombination, transposable elements, semantics, information theory
Introduction
Modern biochemistry has shown that the cell is a much more fascinating piece of machinery than ever would have been expected a century earlier. It has moved in our understanding from being a blob of protoplasm to an intricate wonder of nanotechnology.
Likewise, our understanding of the genome and its intricacies has increased by leaps and bounds over the last few decades. While it was previously thought that only protein-coding genes would be relevant, it is now known that the regulation of protein-coding genes is just as important, if not more so. While it was previously thought to be intellectually irresponsible to consider a biological function for transposable elements, we are now understanding their ubiquity and importance in shaping the genome (Sternberg 2002).
Our understanding of the process of mutation is undergoing a similar revolution. Historically, creationists and evolutionists have been in agreement that mutations which occur in organisms are haphazard—that is, there is no designed purpose for them. However, this understanding is beginning to change. It is becoming increasingly apparent that the cell itself can induce mutational processes in the right genes to produce beneficial changes at appropriate times. The most well-studied of these systems is the somatic hypermutation (SMH) system in the vertebrate immune system. In order to increase the binding potential of immunoglobulins to antigens, the cells direct mutations to a specific region of a specific gene in order to produce immunoglobulins that have a higher affinity to the antigens (Papavasiliou and Schatz 2002).
This mutational process skips the region of the gene which attaches to the B-cell entirely, and focuses only on the region of the gene which binds to the antigen. It is not deterministic—that is, the specific changes which are made seem to be stochastic—but the changes are focused to the right gene in the right situation, bypassing well over 99.99% of the genome and focusing on the correct few hundred base pairs which would matter.
More and more examples of focused mutations have been explored. Some transposable elements are triggered in direct response to specific cell stressors. For instance, Hall (1999) showed that E. coli can use insertion sequences to activate the gene required to metabolize beta-glucoside sugars.
Simple sequence repeats (SSRs) have shown to mutate primarily in copy-number, acting as a genomic tuning knob or state switch. King, Trifonov, and Kashi (2006) enumerate several, including an AC repeat in a promoter which causes variation in body weight in Angus beef cattle, and an AC repeat in tilapia fish with similar effects. Bayless and Moxon (2006) report that a 4-nucleotide repeat (CAAC) can cause the lic2A gene of H. influenzae to switch between three states—off, low expression, and high expression—as it alters the reading frame on which the ATG start codon is found.
Semi-palindromic DNA often points to potentially beneficial areas where mutations might take place. For instance, one study of the genetic adaptation of E. coli to low glucose concentrations found identical modifications of the mgl operator sequence (mglO) in multiple populations (Notley-McRobb and Ferenci 1999). These were later found to be in loops that were near stem-loop structures in DNA, formed by semi-palindromic sequences (Wright 2004).
Caporale has synthesized this new research into what she terms an “implicit genome.” That is, a genome has an implicit range of mutations which are likely to occur, and these mutations are part of evolutionary “strategies” for organisms to survive changing environments (Caporale 1999; Caporale 2006).
While Caporale and others attribute the creation of implicit genomes to indirect selection, this idea also seems to play well into a creationary understanding of the way that genomes should work. In fact, understanding cells to have internal mechanisms for large-scale genetic adaptation has been steadily growing in creation thinking over the last decade (Anderson and Purdom 2008; Ashcraft 2003; Bartlett 2006, 2008; Borger 2009; Lightner 2008; Wood and Cavanaugh 2001).
The Need for a Second-Order Classification of Mutations
In the current literature of both creationists and evolutionists, mutations are often classified according to their effect on an organism’s survival within a specific environment (Baumgardner et al 2008; Eyre-Walker and Keightley 2007). A mutation is considered “beneficial” if it helps the organism survive, “deleterious” if it hinders the organism, and “neutral” or “nearly-neutral” if there is little observed effect. This convention of using beneficial/neutral/ deleterious for categorization is useful because it can often be measured directly and quantitatively.
Unfortunately, many people take these categorizations to signify more than can be inferred from them. For example, those who view all mutations as haphazard often view beneficial mutations as evidence of evolution. Likewise, those who view mutations as being possibly internally generated by a guided process within an organism may view the fact that a mutation is beneficial as being prima facie evidence that a mutation was internally generated, and likewise view the fact that a mutation is deleterious to be prima facie evidence that a mutation was haphazard.
Upon closer examination, however, knowing whether a mutation is “beneficial,” “deleterious,” or “neutral” does not by itself tell us whether or not a given mutation occurred according to design or not.
For instance, in order to survive large-scale environmental changes, populations of organisms may keep a supply of organisms with alternative biochemical configurations through a mutational process (Bartlett 2008). Compared to the overall population of organisms, these mutations would actually be phenotypically deleterious, even though they are part of the overall biological design for hedging against possible environmental changes.
Likewise, a phenotypically beneficial mutation is not necessarily part of an overall design. Behe has termed this sort of event “trench warfare”—the mutation may give a phenotypic advantage within a competitive environment, but at the cost of debilitating some important function of the organism. An easy example of this would be sickle-cell anemia—while it may be beneficial in some circumstances because it prevents malaria, the way it debilitates the person who has it overall leads to the conclusion that this was not a designed feature (Behe 2007). Many other such locally beneficial mutations which have an overall deleterious effect on the complexity of a cell’s biochemistry have been documented (Anderson 2005).
Therefore, while the scale of beneficial/neutral/ deleterious works well for a first-order classification of mutations, creation biologists should be looking deeper into a second-order classification based on its consistency with the design of the organism. This classification separates mutations into one of two possibilities—“design-consistent” mutations, and “design-inconsistent” mutations. A “design-consistent” mutation is one which appears to have occurred within the genome’s implicit range, and a “design-inconsistent” mutation is one which appears to be haphazard (that is, philosophically random as described by Bartlett (2008).
Guidelines for Determining Second-Order Classification
The problem with a second-order classification system is that since we do not have total knowledge of the original plan, it makes it difficult to determine whether a mutation was consistent with that design or not. However, creationists can apply theological concepts to achieve a basic understanding of the plan, which can then illuminate our investigations, even in absence of full prior knowledge of the full plan. The notion of a Genesis “created kind” (called baramins in creation biology) is a key theological notion which will aid our investigation. Because God created the animals according to their kinds (Genesis 1:11; 1:21; 1:24; 1:25), it can be presumed that the biological plan is a sort of dynamic stasis, where basic patterns are preserved, yet variance is allowed to aid in both survival of the baramin and the fulfillment of the baramin’s role. The extent of the intended dynamic and static elements of the baramin are not known a priori. The criteria presented here should help creation biologists take the basic theological concepts provided by scripture, and combine these with the data of biology to achieve a fuller understanding of the pattern of life that God created.
The following are several parameters which can be considered which will help make the second-order classification of mutations based on the assumed dynamic stasis of the original kinds.
Metabolic consistency
One of the major ways of determining whether or not a mutation is design-consistent is to look at the internal operation of the organism, examining the consistency of function of organisms with and without the mutation. Design-consistent mutations should maintain internal consistency whether or not they are beneficial within the current environment. It is hard to separate problems which occur from internal metabolic problems as opposed to those caused from the outside. However, a decent test of metabolic consistency in single-celled organisms would be whether or not a mutation caused metabolic problems for an organism which was growing in a nutrient-rich environment free of competition, predation, and toxic compounds. Similar stress-free environments could be constructed for testing the metabolic consistency of mutations in multicellular organisms.
Sickle-cell anemia, for instance, while it is beneficial as far as preventing malaria, causes large-scale functional problems for the organism even in the best environments. Therefore, the mutation causing sickle-cell anemia is metabolically inconsistent, even though it may provide benefit in certain circumstances.
Anderson (2005) pointed out that most drug-resistant mutations of bacteria involve a fitness cost in most normal environments. However, for some of the mutations cited, the fitness cost was not severe, and therefore, by this criteria, would be considered metabolically consistent. In Mycobacterium tuberculosis, for example, some mutations which confer rifampin resistance also do well in normal cultures. In the three mutations isolated by Billington, McHugh, and Gillespie (1999), one mutation type had no relative fitness decreases, and another one had only moderate relative fitness decreases, thus indicating that they are metabolically consistent. One of the mutations had a drastic reduction in relative fitness, indicating that this mutation was probably not metabolically consistent.
Mutational mechanism
A mutation which is in response to a specific stress or group of stresses, or is timed to occur with a particular stage of life for an organism, or for which there is an enzyme whose core function is to produce such a mutation, is likely to be design-consistent. Williams uses mutational mechanism as the primary differentiator between design-consistent mutations, which he terms as “recombination,” and design-inconsistent mutations, which he terms as simply “mutations” (Williams 2005). Other creationists have followed this approach as well (Ashcraft 2003). This criteria by itself is not sufficient, especially considering that the physical cause of many mutations is unknown. In addition, some design-consistent mutational events may be the result of DNA sequence alone, and may occur without specialized enzyme assistance. In addition, some mutational mechanisms may become mistargeted due to either a problem in the mutational mechanism or the sites which are targeted.
An obvious example of a mutational mechanism would be the V(D)J recombination system used in the production of immunoglobulins. This system recombines three types of gene fragments (variable, diversity, and joining) into a nearly-limitless array of immunoglobulins. This system uses a combination of a pair of enzymes (RAG1 and RAG2) and a recombination signal sequence (RSS), to cut and splice gene fragments at the appropriate locations (Market and Papavasilious 2003). However, this targeting mechanism can also lead to cancers when cryptic RSS sequences (sequences similar to the RSS used by V(D)J reactions, but in other parts of the genome) become accessible to the recombination enzymes (Schlissel, Kaffer, and Curry 2006).
As mentioned earlier in the paper, Hall (1999) showed that E. coli can use insertion sequences to activate the gene required to metabolize beta-glucoside sugars. Because this mutation only occurs under the conditions where the mutation is needed, it can be inferred that this is the result of a cellular mechanism.
Mutation rate
A mutation which occurs at a significantly higher rate than the average mutation rate for the organism is likely to be design-consistent. This points to the mutation being part of the phenotype of the organism, rather than the mutation occurring arbitrarily. Even though we are living in a post-Fall world, we assume that most of our biological systems function properly according to their design on a daily basis. Thus, a high mutation rate, especially across an entire population, is suggestive that a mutation is design-consistent.
As with other considerations, this one must be used with care. Problematic mutations can cause a site to become mutationally active when it should not be. One test suggested by Bartlett (2008) in differentiating between a design-consistent and a design-inconsistent hotspot would be to compare the average fitness effects of mutations at that hotspot with the average fitness effects of induced arbitrary mutations. If the mutations in the hotspot tend to be metabolically consistent, and their relative fitness is greater than the relative fitness of organisms with arbitrary mutations, then this is evidence of a mutational hotspot being design-consistent.
In the human genome, this appears to hold as a general pattern. Chuang and Li (2004) have noted that mutational hotspots tend to occur in genomic regions involved in extracellular communication, while mutational coldspots tend to occur in cellular housekeeping functions. Thus, the mutations seem to be targeted at areas that would provide adaptation to new environments, and silenced at areas that would more likely cause metabolic inconsistencies.
Reversibility
If a mutation is easily reversible (that is, the frequency of reversion is significantly greater than what would be expected from the overall mutation rate), this is good evidence that the mutation is design-consistent. If one of the purposes of mutations is to provide a hedge against environmental changes, then it would be reasonable to think that if the hedge is successful, the organism needs to be able to make a future hedge of going back to the original configuration. Therefore, reversibility is a key indicator of design-consistent mutations.
SSRs are quite interesting subjects because, in general, they are both highly mutable as well as being easily reversible. Historically, they have been viewed as evolutionary byproducts, or junk DNA. However, current research is continually finding new ways in which SSRs allow for the genome to adapt to changing circumstances (Kashi and King 2006).
Preservation of genomic semantics
Every day we are learning more and more about the semantics of the genome. The genome’s semantics can be considered its internal logic—how different sequences serve to regulate and format the genome’s function and adaptation (Sternberg 2002). In computer programming, a program’s semantics is the set of rules, conventions, and axioms which are assumed to hold true as the program progresses.
Many mutations that lead to sickness are those which cause certain regions of the genome to operate in a semantically different way than before. For instance, some mutations have transformed a non-SSR sequence into an SSR. This causes the cell to modify the SSR’s copy number where it would not have before (King, Trifonov, and Kashi 2006). One example of this is in a heritable form of colorectal cancer. In this cancer, a T to A mutation creates a novel SSR. In later generations, the expansion of this SSR leads to colorectal cancer (Laken et al 1997). Other types of signals may also be created improperly due to mutation, which could cause other systems to act on the wrong site (see for instance the previous discussion about V(D)J recombination).
Therefore, in many cases, the altering of the semantics of the genome often points to a design-inconsistent mutation. However, this criteria should be used with care. DNA might contain a large quantity of meta-information. Meta-information is “information about information.” If, for instance, the cell has sufficient meta-information about the roles of different DNA sequences (as opposed to simply the sequences themselves), there is no reason to think that it could not add or remove semantic elements as needed, using the meta-information as a guide. Dipterans, for example, have a class of transposable elements known as mini-me elements (microsatellite initiating mobile elements), which are retrotransposons that contain primers for SSRs (Wilder and Hollocher 2001). The functions of these are not well-characterized, but their abundance within Diptera (comprising, for instance, 1.2% of the Drosophila melanogaster genome), indicates that they are probably a part of the organism’s design.
Since scientific knowledge of genome semantics is still in its infancy, other factors such as metabolic consistency and precision should be considered in determining whether or not a violation of genomic semantics has taken place, or if the mutation is simply a part of a higher-level or undiscovered semantic within the cell.
Other Considerations
These heuristic guidelines are certainly not complete, nor do any of them stand on their own. Nonetheless, when used in combination, they can open up a new way for creationists to look at mutational processes within cells. Hopefully as this categorization is applied, additional considerations will be supplied to this list.
Also note that the weighting of the various criteria are dependent on our understanding of God’s general plan for organisms. If the goal is a dynamic stasis, then reversibility should be a heavily weighted factor. If the goal is for life to proceed according to a specific direction (for instance, see Gene 20091), then mutation rate should be weighted higher.
Conclusions
While the first-level classifications of mutations as beneficial, neutral, or deleterious is certainly useful, the perspective of creation biology can add additional depth by looking at whether or not a mutation is consistent with the cell’s overall architecture, or if the mutation goes against that architecture. Only distinguishing between mutations on the basis of immediate fitness gains or losses can blind biologists to the architecture in which those mutations are made. By differentiating between design-consistent and design-inconsistent mutations, creation biologists can better understand the overall architecture of each baramin, and ultimately understand both God’s purposes for these organisms as well as how humans can best cooperate with their designs in ecological and biotechnical endeavors.
Appendix—Experimental Methodology
A basic outline of how an experimenter could test each of these criteria is given below. Because experimental mutational studies are generally easier on single-celled organisms which can be grown in a laboratory setting, these methodologies are tailored to that environment. It is the hope of this author that further work can expand upon these methodologies and demonstrate the most effective ways of experimentally determining whether mutations are design-consistent or design-inconsistent.
- Metabolic Consistency—to determine metabolic consistency, plate both the wild-type and the mutant-type on several different media, and examine the relative fitness of the mutant-type on each. Compare the best fitness scores of both the wild-type and mutant-type. If the relative fitness of the mutant-type on its best media is significantly lower than the fitness of the wild-type on its best media, then the mutation is probably metabolically inconsistent. Further experimental work comparing relative fitness values of mutations in the context of the other criteria is required for determining what the relative fitness value ranges should be for a mutation to be considered metabolically consistent.
- Mutational Mechanism—this criteria is easy to rule in (by knowing a mutational mechanism capable of producing the mutation) but hard to rule out. If no known mutational mechanism causes the mutation, but a specific environmental inducer causes a specific or semi-specific mutational response, then it is reasonable to infer a mutational mechanism even without knowing what it is.
- Mutation Rate—the rate of mutation should be compared to other sites within this organism’s genome. If the production of this mutation is significantly higher than the average site for the organism, then a high mutation rate is established.
- Reversibility—the reversion rate for a given mutation should be compared to the mutation rate for the mutant’s original production. It is best if the specific gene sequences can be obtained for the wild-type, the mutant, and the reversion, in order to determine if the reversion is a true genetic reversion, or if a different mutation supplied the missing functionality. If the original mutation was a spontaneous mutation by the fluctuation test (Luria and Delbrück 1943) or the Lederberg test (Lederberg and Lederberg 1952), then it is reasonable to think that the reversions should be spontaneous, too. If the original mutation is adaptive, however, then the reversion mutation might require a specific environmental signal to occur, and it might take some effort for the experimenter to determine what this signal is.
- Preservation of Genomic Semantics—the interplay of the mutant sequence and the surrounding genetic context should be examined. Known genetic sequence motifs should be identified, as well as the mutation’s impact on those motifs. If the mutation corrupts an existing motif, then it is likely that the mutation is violating genomic semantics. However, if a set of mutations are consistently altering motifs in narrow ways, then it is possible that the mutations are following an as-yet-unknown semantic rule within the cell.
References
Anderson, K. L. 2005. Is bacterial resistance to antibiotics an appropriate example of evolutionary change? Creation Research Society Quarterly 41:318–326. Retrieved from, http://www.creationresearch.org/crsq/articles/41/41_4/bact_resist.htm.
Anderson, K. L. and G. Purdom. 2008. A creationist perspective of beneficial mutations in bacteria. In Proceedings of the Sixth International Conference on Creationism, ed. A. A. Snelling, pp. 73–86. Pittsburgh, Pennsylvania: Creation Science Fellowship. Retrieved from, http://www.answersingenesis.org/articles/aid/v4/n1/beneficial-mutations-in-bacteria.
Ashcraft, C. A. 2003. Genetic variability by design. TJ 18, no. 2:98–104. Retrieved from, http://creation.com/images/pdfs/tj/j18_2/j18_2_98-104a.pdf.
Bartlett, J. L. 2006. Metaprogramming and genomics. Occasional Papers of the BSG 8:17–18. Retrieved from, http://www.creationbiology.org/content.aspx?page_id=22&club_id=201240&module_id=36813.
Bartlett, J. L. 2008. Statistical and philosophical notions of randomness in creation biology. Creation Research Society Quarterly 45, no. 2:91–99.
Baumgardner, J., J. Sanford, W. Brewer, P. Gibson, and W. ReMine. 2008. Mendel’s Accountant: A new population genetics simulation tool for studying mutation and natural selection. In Proceedings of the Sixth International Conference on Creationism, ed. A. A. Snelling, pp. 87–98. Pittsburgh, Pennsylvania: Creation Science Fellowship. Retrieved from, http://www.icr.org/i/pdf/technical/Mendels-Accountant.pdf.
Bayless, C. D. and E. R. Moxon. 2006. Repeats and variation in pathogen selection. In The implicit genome, ed. L. H. Caporale, pp. 54–76. New York: Oxford University Press.
Behe, M. J. 2007. The edge of evolution. New York: Free Press.
Billington, O. J., T. D. McHugh, and S. H. Gillespie. 1999. Physiological cost of rifampin resistance induced in vitro in Mycobacterium tuberculosis. Antimicrobial Agents and Chemotherapy 43:1866–1869. Retrieved from, http://aac.asm.org/cgi/content/abstract/43/8/1866.
Borger, P. 2009. The design of life: Part 3—an introduction to variation-inducing genetic elements. Journal of Creation 23, no. 1:99–106. Retrieved from, http://www.evoinfo.org/Publications/A/Borger3.pdf.
Caporale, L. H. 1999. Chance favors the prepared genome. Annals of the New York Academy of Sciences 870:1–21.
Caporale, L. H. 2006. An overview of the implicit genome. In The implicit genome, ed. L. H. Caporale, pp. 3–25. New York: Oxford University Press.
Chuang, J. H. and H. Li. 2004. Functional bias and spatial organization of genes in mutational hot and cold regions in the human genome. PLoS Biology 2, no. 2:e29. doi:10.1371/journal.pbio.0020029. Retrieved from, http://www.plosbiology.org/article/info:doi%2F10.1371%2Fjournal.pbio.0020029.
Eyre-Walker, A. and P. D. Keightley. 2007. The distribution of fitness effects of new mutations. Nature Reviews Genetics 8:610–618. Retrieved from, http://www.lifesci.sussex.ac.uk/home/Adam_Eyre-Walker/Website/Publications_files/EWNRG07.pdf.
Gene, M. 2009. The increasing hydrophobicity effect. Retrieved from, http://designmatrix.wordpress.com/2009/06/27/theincreasing-hydrophobicity-effect/.
Hall, B. G. 1999. Transposable elements as activators of cryptic genes in E. coli. Genetica 107:181–187.
Kashi, Y. and D. G. King. 2006. Simple sequence repeats as advantageous mutators in evolution. TRENDS in Genetics 22:253–259. Available online at doi:10.1016/j.tig.2006.03.005.
King, D. G., E. N. Trifonov, and Y. Kashi. 2006. Tuning knobs in the genome: Evolution of simple sequence repeats by indirect selection. In The implicit genome, ed. L. H. Caporale, pp. 77–90. New York: Oxford University Press.
Laken, S. J., G. M. Petersen, S. B. Gruber, C. Oddoux, H. Ostrer, F. M. Giardiello, S. R. Hamilton, H. Hampel, A. Markowitz, D. Klimstra, S. Jhanwar, S. Winawer, K. Offit, M. C. Luce, K. W. Kinzler, and B. Vogeistein. 1997. Familial colorectal cancer in Ashkenazim due to a hypermutable tract in APC. Nature Genetics 17:79–83.
Lederberg, J. and E. M. Lederberg. 1952. Replica plating and indirect selection of bacterial mutants. Journal of Bacteriology 63:399–406. Retrieved from, http://profiles.nlm.nih.gov/BB/A/B/F/J/_/bbabfj.pdf.
Luria, S. E. and M. Delbrück. 1943. Mutations of bacteria from virus sensitivity to virus resistance. Genetics 28:491–511.
Lightner, J. K. 2008. Karyotype variability within the cattle baramin. Answers Research Journal 1:77–88. Retrieved from, http://www.answersingenesis.org/articles/arj/v1/n1/karyotype-variability-cattle.
Market, E. and F. N. Papavasilious. 2003. V(D)J recombination and the evolution of the adaptive immune system. PLoS Biology 1, no. 1:24–27. Retrieved from, http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.0000016.
Notley-McRobb, L. and T. Ferenci. 1999. Adaptive mglregulatory mutations and genetic diversity evolving in glucose-limited Escherichia coli populations. Environmental Microbiology 1:33–43.
Papavasilious, F. N. and D. G. Schatz. 2002. Somatic hypermutation of immunoglobulin genes: Merging mechanisms for genetic diversity. Cell 109(2): S35–S44. Retrieved from, http://dx.doi.org/10.1016/S0092-8674(02)00706-7.
Schlissel, M. S., C. R. Kaffer, and J. D. Curry. 2006. Leukemia and lymphoma: A cost of doing business for adaptive immunity. Genes and Development 20:1539–1544. Retrieved from, http://genesdev.cshlp.org/content/20/12/1539.full.
Sternberg, R. V. 2002. On the roles of repetitive DNA elements in the context of a unified genomic-epigenetic system. Annals of the New York Academy of Sciences 981:154–188.
Wilder, J. and H. Hollocher. 2001. Mobile elements and the genesis of microsatellites in Dipterans. Molecular Biology and Evolution 18:384–392. Retrieved from, http://mbe.oxfordjournals.org/cgi/content/full/18/3/384.
Williams, A. 2005. Inheritance of biological information Part III: Control of information transfer and change. TJ 19(3):21–28. Retrieved from, http://creation.com/images/pdfs/tj/j19_3/j19_3_21-28.pdf.
Wood, T. C. and D. P. Cavanaugh. 2001. A baraminological analysis of subtribe Flaveriinae (Asteracae: Helenieae) and the Origin of Biological Complexity. Origins 52:7–27. Retrieved from, http://www.grisda.org/origins/52007.htm.
Wright, B. E. 2004. Stress-directed adaptive mutations and evolution. Molecular Microbiology 52, no. 3:643–650. Retrieved from, http://www3.interscience.wiley.com/cgibin/fulltext/118789850/HTMLSTART.