Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Lay Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (24)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Kotiaho, J. S.
Right arrow Articles by Tomkins, J. L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kotiaho, J. S.
Right arrow Articles by Tomkins, J. L.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Behavioral Ecology Vol. 12 No. 5: 553-557
© 2001 International Society for Behavioral Ecology

The discrimination of alternative male morphologies

Janne S. Kotiaho and Joseph L. Tomkins

Department of Zoology, University of Western Australia, Nedlands, WA 6907, Australia

Address correspondence to J.S. Kotiaho, who is now at the Department of Biological and Environmental Sciences, University of Jyväskylä, PO Box 35, FIN-40351, Jyväskylä, Finland. E-mail: jkotiaho{at}jyu.fi . J.L. Tomkins is now at the Department of Environmental and Evolutionary Biology, University of St. Andrews, St. Andrews, Fife, Scotland.

Received 10 March 2000; revised 4 August 2000; accepted 5 November 2000.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 REFERENCES
 
Male dimorphisms represent alternative selective regimes within a sex. As such, they can be used as a powerful tool in testing evolutionary theory. However, to realize this potential, we need to be able to accurately discriminate individuals into two separate morphs. In this article we discuss the existing methods and propose a new one. We test our method with data from three dimorphic species and compare these results to results with existing methods. We conclude that existing methods often mis-classify a large proportion of individuals, but applying our method notably reduces these errors.

Key words: male dimorphism, morph discrimination, sexual selection, alternative reproductive strategies.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 REFERENCES
 
Alternative reproductive behaviors are frequently associated with two or more distinct morphologies that occur within a sex (Alcock, 1996bGo; Emlen, 1997Go; Gross, 1985Go, 1991Go). In a few species, the morphologies associated with the different behavioral tactics are extreme, such as the order of magnitude difference in body mass of par and anadromous salmon (Gage et al., 1995Go; Gross, 1985Go). However, in most species dimorphic variability is limited to specific traits that vary discontinuously with body size. Discriminating between one morph and another in these species can be difficult, and this stimulated early researchers to formulate explicit methods of detecting dimorphisms. The first of these was the inspection of the modality of the character distribution; bimodality in character length (e.g., forceps lengths in the earwig Forficula auricularia) associated with unimodality in a linear measure of body size was considered indicative of a dimorphism (Bateson and Brindley, 1892Go; Huxley, 1932Go).

Because male dimorphisms represent alternative selective regimes within a sex, they can be used as powerful tests of evolutionary theory (Gage et al., 1995Go; Gross, 1996Go; Simmons et al., 1999Go; Tomkins and Simmons, 1996Go). However, in modern evolutionary biology, visual examination of the modality of character distributions are frequently inadequate and imprecise. This is because in many cases the character distributions are not strictly discontinuous, but instead all character sizes are expressed by at least some individuals (Figures 1,2,3; see also Eberhard and Gutiérrez, 1991Go; Emlen, 1996Go; Simmons et al., 1999Go; Tomkins, 1999Go; Tomkins and Simmons, 1996Go). Thus, to test evolutionary theories with dimorphic species, we need accurate statistical methods to detect dimorphisms and discriminate between alternative phenotypes.



View larger version (23K):
[in this window]
[in a new window]
 
Figure 1 The relationship between pronotum width and horn length from (A) field and (B) laboratory populations of O. taurus. (A) The vertical line illustrates the switch point (5.100) derived from the Eberhard and Gutiérrez (1991Go) original model (model 1), and the horizontal line illustrates the switch point (0.22) derived from our modification of the model (model 2). (B) The vertical line illustrates the switch point (5.135) derived from model 1, and the horizontal line illustrates the switch point (0.31) derived from model 2.

 


View larger version (24K):
[in this window]
[in a new window]
 
Figure 2 The relationship between pronotum width and horn length in O. binodis. The vertical line illustrates the switch point (5.565) derived from the Eberhard and Gutiérrez (1991Go) original model, and the horizontal line illustrates the switch point (0.48) derived from our modification of model 1 (model 2).

 


View larger version (25K):
[in this window]
[in a new window]
 
Figure 3 The relationship between pronotum width and forceps length in F. auricularia. The vertical line illustrates the switch point (1.960) derived from the Eberhard and Gutiérrez (1991Go) original model, and the horizontal line illustrates the switch point (4.50) derived from our modification of model 1 (model 2).

 

Eberhard and Gutiérrez (1991Go) pioneered a method of statistically detecting and testing for dimorphisms. Before this method only a handful of studies had attempted any statistical determination between alternative morphologies (Cook, 1987Go; Eberhard, 1987Go; Goldsmith, 1985Go; Ollason, 1972Go). Eberhard and Gutiérrez's (1991Go) model (which will be referred to as model 1) is as follows:

(1)
in which Y and X are a linear measure of character size and a linear measure of body size, respectively; XD is the proposed switch point (the body size at which one phenotype gives way to another); D = 0 if X < XD, D = 1 if X >= XD; {alpha} is constant; ß is the regression coefficient, and {epsilon} is the error (Eberhard and Gutiérrez, 1991Go).

This model provides a statistical test for the existence of dimorphic variation in a character associated with body size. First, model 1 tests whether there is a body size switch point, XD, at which the distribution of the character size Y becomes discontinuous. This switch point can be found by first iterating the XD that gives the best fit (highest R2) for model 1 and then testing if this XD fitted in model 1 gives a regression coefficient, ß3, that is significantly different from zero (Eberhard and Gutiérrez, 1991Go).

Second, model 1 tests if there is a body size switch point, XD, at which the linear slope between body size, X, and character length, Y, changes. If the ß3 of the previous test was significant, then the change in the slope at the body size XD may be examined by testing the regression coefficient ß2 against zero. A significant deviation from zero indicates that, in addition to the discontinuity of the character length, there is also a change in the slope. If the ß3 of the previous test was not significant, then the term ß3D in model 1 may be left out, and the change in the slope at the body size XD may be examined with the reduced model by testing the regression coefficient ß2 against zero. The above methods enable the statistical determination of whether a dimorphism exists, and they also establish the exact body size at which individuals are most likely to switch from one morph to the other.

Male dimorphisms (in insects particularly) tend to be conditionally expressed (Eberhard, 1982Go), or what Gross (1996Go) has described as "status dependent." The condition for the expression of alternative behaviors and morphologies is generally body size, which is linked to the status of the individual. For example, in the dung beetle Onthophagus acuminatus the unconscious strategy of beetles can be summarized as two tactics: if small, sneak copulations and do not grow horns, but if large, grow horns and guard females (Emlen, 1997Go). Thus, body size is usually the variable underlying the expression of many dimorphic traits. Model 1 therefore provides the important information about where the body size switch point lies, at which one morph changes to another. The accuracy of this methodology is, however, dependent on a tight relationship between body size and the dimorphic trait (Eberhard and Gutiérrez, 1991Go; Tomkins, 1999Go; Tomkins and Simmons, 1996Go). When there is variance in the switch points of individuals within a population, either environmental or genetic in origin (Tomkins, 1999Go), a large overlap in body size will exist, and many individuals will be misclassified by model 1 (Eberhard and Gutiérrez, 1991Go; Tomkins, 1999Go; Tomkins and Simmons, 1996Go).

For example, in Figure 1, the vertical lines in the graphs indicate the best switch point based on model 1. To the left of the line, individuals are classified as minor morphs, and to the right of the line as major morphs. It is evident that some of the individuals even with the longest horns are erroneously classified as minors (top left of the lines). In addition, some individuals with very small horns are classified as majors (bottom right of the lines). Similar misclassification of minors as majors is evident in Figures 2 and 3. Unfortunately, the overlaps between morphs are very common (Alcock, 1996aGo; Eberhard and Gutiérrez, 1991Go; Emlen, 1996Go; Hunt and Simmons, 1998aGo; Moczek and Emlen, 1999Go; Simmons et al., 1999Go; Tomkins and Simmons, 1999Go) and are the reason that further statistical discrimination is needed.

Thus, for establishing the body size switch point for a dimorphic trait, which is essential for understanding the ontogeny and evolution of dimorphisms, model 1 is the appropriate model. However, model 1 is not as effective in identifying individuals as either belonging to one morph or another. We suggest a modified methodology that specifically aims for the discrimination of morphs based on the dimorphic character itself, rather than on body size. By modifying model 1 to substitute Y with X and X with Y, we find the switch point in the dimorphic character size that defines the dimorphism, rather than the switch point in body size that defines the dimorphism. In this model (which will be referred to as model 2),

(2)
The terms are as defined in the original model 1 with the exception that the suggested switch point is YD and D = 0 if Y < YD, and D = 1 if Y >= YD. Repeating the calculations discussed with model 1 allows one to statistically test the existence of dimorphism and to find the best switch point that determines the dimorphism similarly as with the model 1. The only difference between the results is that now the switch point is found directly for the dimorphic character itself rather than for a correlate of it—namely, body size.

It may seem inappropriate to use linear regression models and regress X on Y (i.e., independent body size on dependent character size). However, because in the majority of morphometric studies, Y is not strictly dependent on X, and X is not measured without error (the two basic assumptions of linear regression models; Sokal and Rohlf, 1981Go; Zar, 1996Go), there is statistical justification for regressing X on Y as well as Y on X.

We illustrate the utility of our modification of the model 1 in the discrimination between morphs by comparing the two models with data from three species: two species of dung beetles with a horn dimorphism (Onthophagus taurus [two samples] and Onthophagus binodis) and a species of earwig with a forceps dimorphism (Forficula auricularia). The first sample of O. taurus (Figure 1A) consists of individuals from a field population in Margaret River in southwestern Western Australia. These data were first presented in Simmons et al. (1999Go). The second sample of O. taurus (Figure 1B) consists of second generation, laboratory-reared individuals also originating from Margaret River. The sample of O. binodis consists of first-generation laboratory-reared individuals originating from Walpole in southwestern Western Australia. The F. auricularia were collected from West Wideopen Island in the Farnes group on the northeastern coast of England, and the data were presented first in Tomkins and Simmons (1996Go).

The best fitting switch point for field-collected O. taurus, using model 1, occurred at the pronotum width of 5.100 mm, and the dimorphism at this point was best characterized as being discontinuous (ß3 significantly different from zero) and having a change in the slope (ß2 significantly different from zero; Tables 1 and 2). Using model 2, the best fitting switch point occurred at the horn length of 0.201 mm, and the dimorphism at this point was best characterized as being discontinuous (ß3 significantly different from zero) and having a change in the slope (ß2 significantly different from zero) (Tables 1 and 2). Using model 1, the best fitting switch point for the laboratory-reared O. taurus occurred at the pronotum width of 5.135 mm, and the dimorphism at this point was best characterized as being discontinuous 3 significantly different from zero) and having a change in the slope (ß2 significantly different from zero; Tables 1 and 2). Using model 2, the best fitting switch point occurred at the horn length of 0.310 mm, and the dimorphism at this point was best characterized as being discontinuous 3 significantly different from zero) and having a change in the slope (ß2 significantly different from zero; Tables 1 and 2).


View this table:
[in this window]
[in a new window]
 
Table 1 Statistical tests for dimorphism in O. taurus from a field and laboratory population: multiple regressions
 

View this table:
[in this window]
[in a new window]
 
Table 2 Statistical tests for dimorphism in O. taurus from field and laboratory populations: ANOVAs
 

The best fitting switch point for O. binodis using model 1 occurred at the pronotum width of 5.565 mm, and the dimorphism at this point was best characterized as being continuous (ß3 not significantly different from zero) but having a change in the slope (ß2 significantly different from zero) (Tables 3 and 4). Using model 2, the best fitting switch point occurred at the horn length of 0.48 mm, and the dimorphism at this point was best characterized as being continuous (ß3 not significantly different from zero) but having a change in the slope (ß2 significantly different from zero; Tables 3 and 4).


View this table:
[in this window]
[in a new window]
 
Table 3 Statistical tests for dimorphism in O. binodis: multiple regressions
 

View this table:
[in this window]
[in a new window]
 
Table 4 Statistical tests for dimorphism in O. binodis: ANOVAs
 

The best fitting switch point for F. auricularia, based on model 1, occurred at the pronotum width of 1.959 mm, and the dimorphism at this point was best characterized as being discontinuous 3 significantly different from zero) and having a change in the slope (ß2 significantly different from zero; Tables 5 and 6). Using model 2, the best fitting switch point occurred at the forceps length of 4.50 mm, and the dimorphism at this point was best characterized as being discontinuous (ß3 significantly different from zero) and having a change in the slope (ß2 significantly different from zero; Tables 5 and 6).


View this table:
[in this window]
[in a new window]
 
Table 5 Statistical tests for dimorphism in F. auricularia: multiple regressions
 

View this table:
[in this window]
[in a new window]
 
Table 6 Statistical tests for dimorphism in F. auricularia: ANOVAs
 

In Figures 1,2,3, we have plotted the horn lengths of O. taurus and O. binodis and the forceps length of F. auricularia on their respective pronotum widths. As pointed out earlier, it is immediately clear from the figures that the switch points calculated by using model 1 misclassify a substantial proportion of the males. However, by using our model 2, misclassification is notably reduced.

As morphological dimorphisms are frequently associated with alternative reproductive behaviors (Alcock et al., 1977Go; Cook, 1990Go; Emlen, 1997Go; Goldsmith, 1987Go; Hunt and Simmons, 1998bGo, 2000Go; Hunt et al., 1999Go), only a combination of behavioral and morphological studies can ultimately confirm the correct position of a switch point. In dung beetles, the dimorphism in horn lengths has been shown to be associated with different behaviors. In O. acuminatus, major morphs with large horns defend the tunnels in which a female constructs a brood mass, while minor morphs try to sneak copulations (Emlen, 1997Go). In O. taurus and O. binodis, major and minor morphs use different tactics in helping the female in brood mass construction; major males help females, whereas minor males do not help (Cook, 1990Go; Hunt and Simmons, 1998bGo). However, only one study has examined whether the change in male morphology from one morph to another coincides with the change from one behavior to another. Hunt and Simmons (2000Go) showed that in O. taurus the start of horn development coincided with a behavioral change from a nonhelping minor to a helping major morph; furthermore, the change in helping was not continuous but showed an abrupt discontinuity. This abrupt change in behavior accompanied with much less abrupt change in horn lengths (Figure 1) suggests that the biologically correct switch point from minor to major morph in O. taurus occurs immediately at the start of horn development and not across a body size threshold per se.

From Hunt and Simmons's study (2000Go), it is obvious that model 1 sometimes misclassifies many males of O. taurus (Figure 1). However, our model 2 performs better and provides us with a horn length switch point that coincides well with the behavioral switch described by Hunt and Simmons (2000Go).

Other methods of describing dimorphisms in Onthophagus have been used. For example, Emlen (1996Go) modified logistic regressions and nonlinear regressions to determine a switch point for both the pronotum width and the horn length (inflection point and curve height, respectively, in Emlen, 1996Go). However, in these models the switch point is located at the inflection point of the sigmoidal distribution of horn lengths plotted on pronotum width (Emlen, 1996Go). Although this model is useful for calculating relative horn size, it may be unsuitable for finding the correct switch point. It is likely that the positions of behavioral and morphological switch points in O. taurus are also applicable to other Onthophagus species with approximately similar sigmoidal relationships between horn length and pronotum width. Hence, the determination of switch points based on inflection of the sigmoid (Emlen, 1996Go) is likely to overstimate the position of the real switch point and misclassify a proportion of small majors as minors.

Eberhard and Gutiérrez's (1991Go) model (model 1) is based on the notion of conditional expression determined by body size. This model is not redundant because it determines the lateral position of the point across which most individuals are likely to change phenotypes, providing a vital component of understanding the evolution and ontogeny of dimorphisms. Our modification (model 2) is not based on body size but rather directly on the dimorphic character of interest and thus is not affected by overlaps in body size between morphs. In all of the four data sets that we analyzed, our modification (model 2) performed better in classifying the morphs than the original model. Our modification provides additional information about the dimorphism and has its primary advantage in discriminating accurately between the alternative phenotypes.


    ACKNOWLEDGEMENTS
 
We thank John Hunt and Leigh W. Simmons for comments on the manuscript. J.S.K. was funded by the Academy of Finland and J.L.T. was funded by a University Postdoctoral Research Fellowship from the University of Western Australia.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 REFERENCES
 
Alcock J, 1996a. Provisional rejection of three alternative hypotheses on the maintenance of a size dichotomy in Dawson's burrowing bee, Amegilla dawsoni (Apidae, Apinae, Anthophorini). Behav Ecol Sociobiol 39: 181-188.[Web of Science]

Alcock J, 1996b. The relation between male body size, fighting, and mating success in Dawson's burrowing bee, Amegilla dawsoni (Apidae, Apinae, Anthophorini). J Zool Soc Lond 239: 663-675.

Alcock J, Jones E, Buchmann SL, 1977. Male mating strategies in the bee Centris pallida (Anthrophoridae: Hymenoptera). Am Nat 111: 145-155.[Web of Science]

Bateson W, Brindley HH, 1892. On some cases of variation in secondary sexual characters, statistically examined. Proc R Soc Lond 1892: 585-594.

Cook DF, 1987. Sexual selection in dung beetles I. A multivariate study of the morphological variation in two species of Onthophagus (Scarabaeidae: Onthophagini). Aust J Zool 35: 123-132.[Web of Science]

Cook DF, 1990. Differences in courtship, mating and postcopulatory behavior between male morphs of the dung beetle Onthophagus binodis Thunberg (Coleoptera: Scarabaeidae). Anim Behav 40: 428-436.[Web of Science]

Eberhard WG, 1982. Beetle horn dimorphism: making the best of a bad lot. Am Nat 119: 420-426.[Web of Science]

Eberhard WG, 1987. Use of horns in fights by the dimorphic males of Ageopsis nigricollis (Coleoptera, Scarabeidae, Dynastinae). J Kans Entomol Soc 60: 504-509.

Eberhard WG, Gutiérrez EE, 1991. Male dimorphisms in beetles and earwigs and the question of developmental constraints. Evolution 45: 18-28.

Emlen DJ, 1996. Artificial selection on horn body-length size allometry in the horned beetle Onthophagus acuminatus (Coleoptera: Scarabidae). Evolution 50: 1219-1230.[Web of Science]

Emlen DJ, 1997. Alternative reproductive tactics and male dimorphism in the horned beetle Onthophagus acuminatus (Coleoptera: Scarabaeidae). Behav Ecol Sociobiol 41: 335-341.[Web of Science]

Gage MJG, Stockley P, Parker GA, 1995. Effects of alternative male mating strategies on characteristics of sperm production in the Atlantic salmon (Salmo salar): Theoretical and empirical investigations. Phil Trans R Soc Lond B 350: 391-399.

Goldsmith SK, 1985. Male dimorphism in Dendrobias mandibolaris Audinet-Serville (Coleoptera: Cerambycidae). J Kans Entomol Soc 58: 534-538.

Goldsmith SK, 1987. The mating system and alternative reproductive behaviours of Dendrobias mandibularis (Coleoptera: Cerambycidae). Behav Ecol Sociobiol 20: 111-115.

Gross MR, 1985. Disruptive selection for alternative life histories in salmon. Nature 313: 47-48.

Gross MR, 1991. Evolution of alternative reproductive strategies: frequency-dependent sexual selection in male bluegill sunfish. Phil Trans R Soc Lond B 332: 59-66.

Gross MR, 1996. Alternative reproductive tactics: diversity within sexes. Trends Ecol Evol 11: 92-98.[Web of Science]

Hunt J, Kotiaho JS, Tomkins JL, 1999. Dung pad residence time covaries with male horn morphology in the dung beetle Onthophagus taurus (Coleoptera: Scarabidae). Ecol Entomol 24: 174-180.

Hunt J, Simmons LW, 1998a. Patterns of fluctuating asymmetry in beetle horns: no evidence for reliable signaling. Behav Ecol 9: 465-470.[Abstract/Free Full Text]

Hunt J, Simmons LW, 1998b. Patterns of parental provisioning covary with male morphology in a horned beetle (Onthophagus taurus) (Coleoptera: Scarabidae). Behav Ecol Sociobiol 42: 447-451.[Web of Science]

Hunt J, Simmons LW, 2000. Maternal and paternal effects on offspring phenotype in the dung beetle Onthophagus taurus. Evolution 54: 936-941.[Web of Science][Medline]

Huxley JS, 1932. Problems of relative growth. London: Methuen.

Moczek AP, Emlen DJ, 1999. Proximate determination of male horn dimorphisim in the beetle Onthophagus taurus (Coleoptera: Scarabidae). J Evol Biol 12: 27-37.

Ollason JG, 1972. A statistical description of structural variation in the cerci of the common earwig (Forficula auricularia). J Zool Lond 167: 153-160.

Simmons LW, Tomkins JL, Hunt J, 1999. Sperm competition games played by dimorphic male beetles. Proc R Soc Lond B 266: 145-150.[Abstract/Free Full Text]

Sokal RR, Rohlf FJ, 1981. Biometry. New York: Freeman.

Tomkins JL, 1999. Environmental and genetic determinants of the male forceps length dimorphism in the European earwig Forficula auricularia L. Behav Ecol Sociobiol 47: 1-8.

Tomkins JL, Simmons LW, 1996. Dimorphisms and fluctuating asymmetry in the forceps of male earwigs. J Evol Biol 9: 753-770.[Web of Science]

Tomkins JL, Simmons LW, 1999. Heritability of size but not symmetry in a sexually selected trait chosen by female earwigs. Heredity 82: 151-157.[Web of Science]

Zar JH, 1996. Biostatistical analysis. Englewood Cliffs, New Jersey: Prentice-Hall.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
ScienceHome page
J. M. Rowland and D. J. Emlen
Two Thresholds, Three Male Forms Result in Facultative Male Trimorphism in Beetles
Science, February 6, 2009; 323(5915): 773 - 776.
[Abstract] [Full Text] [PDF]


Home page
Behav EcolHome page
S. C. Cotter, M. Beveridge, and L. W. Simmons
Male morph predicts investment in larval immune function in the dung beetle, Onthophagus taurus
Behav. Ecol., March 1, 2008; 19(2): 331 - 337.
[Abstract] [Full Text] [PDF]


Home page
Proc R Soc BHome page
J. L Tomkins, J. S Kotiaho, and N. R LeBas
Phenotypic plasticity in the developmental integration of morphological trade-offs and secondary sexual trait compensation
Proc R Soc B, March 7, 2005; 272(1562): 543 - 551.
[Abstract] [Full Text] [PDF]


Home page
Behav EcolHome page
C. D. Kelly
Allometry and sexual selection of male weaponry in Wellington tree weta, Hemideina crassidens
Behav. Ecol., January 1, 2005; 16(1): 145 - 152.
[Abstract] [Full Text] [PDF]


Home page
Behav EcolHome page
J. S. Kotiaho
Sexual selection and condition dependence of courtship display in three species of horned dung beetles
Behav. Ecol., November 1, 2002; 13(6): 791 - 799.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Lay Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (24)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Kotiaho, J. S.
Right arrow Articles by Tomkins, J. L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kotiaho, J. S.
Right arrow Articles by Tomkins, J. L.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?