Mycologia
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Blackwell, M.
Right arrow Articles by Jones, K. G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Blackwell, M.
Right arrow Articles by Jones, K. G.
Agricola
Right arrow Articles by Blackwell, M.
Right arrow Articles by Jones, K. G.
Mycologia, 95(6), 2003, pp. 987-992.
© 2003 by The Mycological Society of America

Extreme morphological divergence: phylogenetic position of a termite ectoparasite


Meredith Blackwell 1

     Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana 70803

Daniel A. Henk

     Biology Department, Duke University, Box 90338, Durham, North Carolina 27708

Kevin G. Jones

     Department of Biology, University of Virginia at Wise, Wise, Virginia 24293

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 

Species of Termitaria are lesion-forming ectoparasites occurring worldwide on a diverse group of termites. The reduced thallus consists of a basal cell layer from which haustorial cells penetrate the termite and a darkly pigmented sporodochium. One species, Termitaria snyderi, has been the subject of several morphological studies, but its phylogenetic position has remained enigmatic. Here we provide evidence of a close relationship between T. snyderi and the morphologically distinct ascomycetes, Kathistes analemmoides and K. calyculata, based on phylogenetic analysis of molecular characters derived from portions of the nuclear-encoded small-subunit ribosomal RNA gene (ssu rDNA) and supplemental evidence from the ß-tubulin gene. Trees were derived using parsimony and maximum-likelihood criteria. Bayesian analysis and parsimony bootstrap methods were used to assess support for the tree nodes.

Key words: asexual fungi, conidial fungi, insect-associated fungi, Isoptera, mitosporic fungi, ssu rDNA, Termitaria


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Termites are unique among insect groups in having a variety of minute ectoparasitic fungal specialists that attach to or penetrate their exoskeletons. Most of the species in 10 different genera are restricted to termites (Blackwell and Rossi 1986Citation). The life histories of the fungi are not well known, and hypotheses of phylogenetic relationships often have been based on scant morphological and life history data (Blackwell and Kimbrough 1976aCitation, bCitation, Blackwell et al 1980Citation, 1986Citation, Blackwell 1994Citation). In addition to the more elaborate thalli of several Laboulbeniomycetes found on termites, the thalli of some termite ectoparasites, such as species of Laboulbeniopsis and Coreomycetopsis, have thalli composed of 3–15 linearly superposed cells; single-celled Amphoromorpha species reported from termites probably are the capilliconidia of species of Basidiobolus (Blackwell and Malloch 1989Citation). None of these fungi has been cultured in agar medium.

A second group of termite fungi has more prominent, albeit reduced, thalli and includes species of Termitaria and Mattirolella. Termitaria snyderi conidia germinate and form thalli on the surface of the termite exoskeleton, and thalli can develop on any exposed body part, including the head, legs, antennae and abdomen. A thallus consists of a basal layer of cells, some or all of which may be haustorial mother cells. Each of the dark, thick-walled haustorial mother cells of the pseudoparenchymatous crust possesses a pale circular area at the base, actually an opening in the cell wall through which penetration pegs exit the cell to enter the host exoskeleton. From the penetration pegs the walled, nucleated haustoria penetrate through the exoskeleton to the basement membrane of the termite integument (Khan and Aldrich 1975Citation). Above the basal layer, the cells of the subhymenium give rise to a compact palisade of conidiogenous cells that are enclosed by a border of dark-walled cells forming a sporodochium (Khan and Kimbrough 1974aCitation). Khan and Aldrich (1973)Citation reported basipetal production of conidia from a fixed locus within a closed conidiogenous cell; conidium release occurred upon synchronous rupture of the conidiogenous cell apices.

In addition to T. snyderi, four other species of Termitaria are known from termite hosts (Thaxter 1920Citation, Kimbrough and Lenz 1982Citation), and the genus is present in most regions where termites occur (Blackwell and Rossi 1986Citation, Hojo et al 2001Citation). Other species of termite parasites in Mattirolella (Kimbrough and Thorne 1982Citation) and Termitariopsis cavernosa on legionary ants (Blackwell et al 1980Citation) have a thallus morphology and ecology similar to those of T. snyderi (Blackwell and Rossi 1986Citation), and we suspect that they are all related. It is, however, the relationships of the morphologically distinguished group to other fungi that is the question needing attention.

Termitaria and Mattirolella were suggested as relatives of several fungi based on morphological features: the flat stromatic growth of the primary thallus has been likened to that of the Asterinales (Thaxter 1920Citation, Kimbrough and Thorne 1982Citation); the darkly pigmented sporodochia with tightly held conidiogenous cells have been compared to Bloxamia; conidium ontogeny in Termitaria and Mattirolella was thought to be similar to the processes observed in Bloxamia truncata, Thielaviopsis, Chalara and Chalaropsis (Khan and Kimbrough 1974aCitation, bCitation). However, morphological and developmental studies merely have served to highlight the unique morphology of species of Termitaria and Mattirolella when compared to other fungi.

The distinctive morphology caused Khan and Kimbrough (1974b)Citation to place Termitaria and Mattirolella in a new order of the Deuteromycotina, Termitariales. The classification of asexual fungi, however, has been de-emphasized since that time because it now is possible to use molecular characters to place asexual fungi among their closest sexual relatives. Here we report information on T. snyderi that provides strong support for a new hypothesis that this species is related to two species of arthropod-associated ascomycetes, Kathistes analemmoides and K. calyculata (Malloch and Blackwell 1990Citation).


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Collection and examination – Thalli of Termitaria snyderi on Reticulitermes flavipes (termites) were collected from loblolly pine (USA, LA: Baton Rouge, Louisiana State University Burden Research Plantation, 4560 Essen Lane, dead branch of Pinus taeda, 9 Jul 1997, Daniel A. Henk 14). Excess material was stored in alcohol as a voucher in LSUM. The collection of the termites was by tapping chunks of termite-infested wood, so that termites rained from the galleries into plastic bags. Termites were killed by freezing and examined for external fungi using a dissecting microscope (100x) (Blackwell 1980Citation).

Culture – Infected termite parts were removed aseptically and placed on 1/2; strength cornmeal agar (7 g Difco cornmeal agar, 12 g DifcoBacto agar, 1 L tap water) or Molisch's agar (15 g agar, 20 g sucrose, 10 g peptone, 0.25 g Mg2SO4, 0.25 g K2HPO4, 1 L tap water), which is used specifically for entomopathogens (Harry Evans pers comm, 1994). A small amount of chlortetracycline hydrochloride (Sigma, St. Louis, Missouri) was sprinkled on the agar surface to retard bacterial growth. The fungal lesions on the agar were observed daily for 2 wk to monitor the slowly growing T. snyderi hyphae and to determine whether contaminants were present.

DNA extraction – DNA extraction from fungal lesions on termites was carried out by adaptation of the protocols of Cenis (1992)Citation and Ferreira and Glass (1996)Citation. Three thalli were dissected from a single infected termite and placed in a 1.5 mL microfuge tube. The dry fungal material was irradiated at 700 W for 5 min in a domestic microwave oven and then dispersed in 100 µL extraction buffer (0.2 M Tris-HCl, pH 8.0; 0.25 M NaCl; 25 mM EDTA; 0.5% SDS) by using a vortex. Further homogenization of the material was performed with a disposable micropestle (Kontes, Vineland, New Jersey), and the extract was incubated at 65 C for 5 min. One-half volume of 7.5 M ammonium acetate was added to the homogenate, followed by incubation at -20 C for 10 min. Cellular debris and precipitated proteins were pelleted by centrifugation (15 min, 16 000x g), and the supernatant, transferred to a clean tube. Nucleic acids were precipitated by addition of an equal volume of propan-2-ol, and recovered by centrifugation, as described above. DNA was rinsed twice with 500 µL 70% ethanol, air dried, and dissolved in 100 µL sterile water. Aliquots of the DNA solution were used undiluted in PCR reactions. To extract DNA from cultured material, the sparse aerial hyphae were collected under a dissecting microscope by contact with a toothpick dipped in sterile water. Adherent hyphae were dispersed in 100 µL of extraction buffer in a 1.5 mL microfuge tube. Subsequent processing of the sample was performed as described above with the omission of the microwave treatment. DNA samples were diluted 1/10 for use in PCR reactions.

DNA amplification, polymerase chain reaction and sequencing – Sequence data were obtained from the 5' end of the ssu rRNA gene and from a portion of the gene encoding ß-tubulin. PCR amplifications were performed with 100 µL PCR Supermix (Gibco BRL, Bethesda, MD) containing 10 µM each of the appropriate primers (see below). For ssu rDNA, primers NS17 (Gargas and Taylor 1992Citation) and NS4 (White et al 1990Citation) were used in a cycling protocol consisting of 94 C for 5 min, followed by 35 cycles of 95 C for 1 min; 54 C for 1 min; and 72 C for 1.5 min. Tubulin amplifications were performed using primers ßt2a and ßt2b (Glass and Donaldson 1995Citation, Donaldson et al 1995Citation) in a cycling procedure consisting of 94 C for 5 min, followed by 35 cycles of 94 C for 1 min; 65 C for 1 min; 72 C for 1 min with a 5 s increase in extension time per cycle. Aliquots of the reactions were subjected to electrophoresis through 1% agarose gels, and PCR products were recovered with a Prep-A-Gene DNA Purification Kit (Bio-Rad Laboratories, Hercules, California). DNA sequencing reactions were performed with primers NS17, NS2, NS3 and NS4 (18S rDNA), and ßt2a and ßt2b for tubulin products, with the Taq dideoxy Terminator Cycle Sequencing Kit (Perkin-Elmer GeneAmp PCR System 2400, PE Applied Biosystems, Foster City, California). Reactions were analyzed on an ABI Prism Model 377 automated sequencer (PE Applied Biosystems, Foster City, California).

Phylogenetic analysis – Sequences of 43 taxa accessioned in GenBank are: Alternaria alternata (U05194), Ambrosiozyma platypodis (L36984), Aspergillus fumigatus (M60300), Balansia sclerotica (U32399), Blastomyces dermititidis (M55624), Capronia mansonii (X79318), Candida albicans (M60302), Candida tropicalis (M55527), Ceratocystis fimbriata (U32418), Ceratocystis virescens (U32419), Chaetomium globosum (U20379), Claviceps paspali (U32401), Cryphonectria parasitica (L42441), Daldinia concentrica (U32402), Diaporthe phaseolorum (L36985), Haematonectria haematococca (U32413), Halosphaeriopsis mediosetigera (U32420), Hypomyces polyporinus (U32410), Hypoxylon atroroseum (U32411), Kathistes analemmoides (AF313767), Kathistes calyculata (AF313768), Leptosphaeria maculans (U04238), Leucostoma persoonii (M83259), Microascus trigonosporus (L36987), Morchella esculenta (L36998), Neurospora crassa (X04971), Ophiostoma ulmi (M83261), Ophiostoma piliferum (U20377), Protomyces inouyei (D11377), Petriella setifera (U32421), Peziza badia (L37539), Pleospora herbarum (U05201), Pyxidiophora sp. 1 (IMI-1989) (AF313769), Pyxidiophora sp. 2 (Baton Rouge, Chalara anamorph) (AY212811*), Saccharomyces cerevisiae (J01353), Sphaerostilbella aureonitens (U32416), Subbaromyces splendens (U63552), Talaromyces flavus (M83262), Taphrina deformans (U00971), Termitaria snyderi (culture) (AY280724* and AY280725, AY278207*), Termitaria snyderi (lesion) (AY212812*, AY278206*), Xylaria curta (U32417) and Xylaria hypoxylon (U20378). Sequences acquired in this study that were aligned with the multialignment program Clustal W (Thompson et al 1994Citation) are marked with an asterisk. The alignment (TreeBASE matrix SN1427–4151) was optimized visually, and ambiguous regions were excluded from the analyses. Ambrosiozyma platypodis, Candida albicans, Candida tropicalis, Protomyces inouyei, Saccharomyces cereviseae and Taphrina deformans were designated outgroup taxa based on the results of many previous studies showing that these groups lie outside the clade of ascocarpic ascomycetes.

Equally weighted maximum-parsimony (MP) analysis was carried out with PAUP 4.0b10 (Swofford 2001Citation) to perform a heuristic search starting from 100 random-addition sequence replicates and swapping with the TBR algorithm. Support using MP was assessed with 1000 bootstrap replicates (Felsenstein 1985Citation). Equally parsimonious trees were compared with maximum-likelihood scores using the Shimodaira-Hasegawa test. Analyses were done with a variety of samplings among diverse ascocarpic ascomycetes, but additional information on the taxa of interest was not obtained. Bayesian Markov Chain Monte Carlo (B-MCMC) analysis was carried out with the program MrBAYES 2.01 (Huelsenbeck and Ronquist 2001Citation). The general time reversible model with gamma shape estimate of rate heterogeneity and estimated base frequencies and invariant sites was implemented in the B-MCMC analysis. The analysis consisted of 200 000 generations of four chains sampling every 100 generations. The first 24 000 generations were discarded as burn in time, leaving 1760 trees as the sample representing the posterior probability distribution.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
DNA obtained from T. snyderi lesions and a putative T. snyderi culture from a different termite in the same collection (DAH 14) were identical for all readable bases except one for about 1000 bp of the ssu rRNA gene and all of the 448 bp of the ß-tubulin gene regions obtained by sequencing. Within the 448 bp segment amplified with ßt2a and ßt2b exons 3 (27 bp, 9 a.a.), 4 (42 bp, 14 a.a.), 5 (54 bp, 18 a.a.), and part of 6 (165 bp, 55 a.a.) and the introns 3 (55 bp), 4 (52 bp), and 5 (53 bp) were found. Comparison of partial sequences of the two genes, one including three highly divergent ß-tubulin intron regions, was the criterion used to establish the identity of the culture as that of T. snyderi. Information on the culture will be discussed elsewhere, but it was essential to this study as a check on the identity of the DNA extracted from the termite lesion.

Parsimony analysis resulted in two most-parsimonious trees that were incongruent only in the placement of Morchella esculenta and Peziza badia. The most-likely tree (-ln L = 6129.43615) (Fig. 1) showed a sister relationship between M. esculenta and P. badia, while the other tree (-ln L = 6129.58604; P = 0.445) showed these taxa as a grade leading to the remaining Euascomycetes. B-MCMC analysis and parsimony bootstrap provided similar levels of support for nodes throughout the tree. Both of the trees placed T. snyderi and the two species of Kathistes within a single clade excluded from the perithecial ascomycetes with 100% Bayesian analysis and bootstrap support.



View larger version (43K):
[in this window]
[in a new window]
 
FIG. 1. One of two most-parsimonious trees. Upper or single numbers indicate support above 50% in 1000 bootstrap replicates with parsimony. Lower numbers represent probability of nodes with support greater than 90% in Bayesian analysis. This tree was considered the better tree by comparison with maximum likelihood

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The close relationship of K. analemmoides and K. calyculata to T. snyderi was not predicted on the basis of morphological traits or life-cycle studies, but rather was discovered by the inclusion of DNA sequences from a variety of ascomycetes targeted in part for comparison because of their associations with arthropods.

Kathistes is a relatively recently described and little known ascomycete genus that contains three species (Malloch and Blackwell 1990Citation). The species are known only from dung and are characterized by hyaline perithecial ascomata with light brown necks, basally clustered evanescent asci, several-septate hyaline ascospores and unique multicellular structures of undetermined function (sporidiomata) associated with the base of mature perithecia. The sporidiomata contain minute conidium-like cells that presumably escape through a "nozzle" formed by a single elongated wall cell, but there is no other information on the function or fate of these cells. Two species, K. analemmoides and K. calyculata, have been studied in culture, and they are characterized by ascospore germination by budding to produce yeast-like growth that may darken in color to tan; only scant mycelial growth occurs in culture with age. Species of Kathistes are associated with arthropods for dispersal in dung habitats, but they are not known to be arthropod parasites or associates of termites.

By contrast T. snyderi has a darkly pigmented sporodochium enclosing tightly packed conidiogenous cells as its only known state and apparently is restricted to termite hosts. Thus, based on our current knowledge, species of Kathistes and Termitaria do not have comparable states in their life cycles. Because the morphological traits could not be compared among the taxa to increase confidence in the grouping of these fungi, B-MCMC was used as a model-based method by which we specifically addressed the potential spurious results of long-branch attraction that are possible under some circumstances in parsimony analysis (Felsenstein 1978Citation). B-MCMC, in fact, did support the hypothesis.

Our strategy to find close relatives of T. snyderi was dependent on the inclusion of a diverse group of insect-associated ascomycetes in our phylogenetic database, and, although we have discovered that K. analemmoides and K. calyculata are its closest relatives, we have not been able to resolve the clade among other ascomycetes, a recurring problem in attempts to discover the relationships of some major ascomycete clades. For example, such a problem occurred in the placement of the Laboulbeniomycetes (represented in Fig. 1 by Pyxidiophora), a group that includes another previously problematic termite parasite, Laboulbeniopsis (Henk et al 2003). Lack of resolution has been particularly frustrating in our attempts to track evolutionary transformations in morphology and life histories in insect-associated lineages with highly divergent morphologies. Our inability to determine the near relatives of the Termitaria-Kathistes clade severely reduces the information we can infer from the tree on morphological traits and life history strategies.

At this time we are unable to suggest additional ascomycetes for sampling but hope for the accidental discovery of other members of the Termitaria-Kathistes clade from among the remaining, largely unsampled, diversity of ascocarpic ascomycetes. However, the common problem associated with reconstructing rapid ancient radiations, probably represented here by the deeper, poorly supported branches in both maximum-parsimony and maximum-likelihood analyses, might require that additional genes be applied to achieve a solution (Lanyon 1993Citation).


    ACKNOWLEDGMENTS
 
This research was supported in part by the National Science Foundation (Phylogeny of Laboulbeniales, DEB-9615520 with REU supplements and DEB-0090301) and a Howard Hughes Medical Institute grant through the Undergraduates Biological Sciences Education Program to Louisiana State University. Drs. Sung-Oui Suh and Alex Weir provided valuable support throughout the study.


    FOOTNOTES
 
1 Corresponding author. E-mail: mblackwell{at}lsu.edu Back

Accepted for publication May 8, 2003.


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Blackwell M., 1980 New records of termite-associated fungi from Georgia. J Invert Path 35:101-104

———. 1994 Minute mycological mysteries: the influence of arthropods on the lives of fungi. Mycologia 86:1-17

———, Kimbrough JW., 1976a A developmental study of the termite-associated fungus Coreomycetopsis oedipus. Mycologia 68:551-558

———, ———. 1976b Ultrastructure of the termite-associated fungus Laboulbeniopsis termitarius. Mycologia 68:541-550

———, Malloch D., 1989 Similarity of Amphoromorpha and secondary capilliconidia of Basidiobolus. Mycologia 81:735-741

———, Rossi W., 1986 Biogeography of fungal ectoparasites of termites. Mycotaxon 25:581-601

———, Samson R, Kimbrough JW., 1980 Termitariopsis cavernosa, gen. et sp. nov., a sporodochial fungus ectoparasitic on ants. Mycotaxon 12:97-104

———, Bridges JR, Moser JC, Perry TJ., 1986 Hyperphoretic dispersal of a Pyxidiophora anamorph. Science 232:993-995[Abstract/Free Full Text]

Cenis JL., 1992 Rapid extraction of fungal DNA for PCR amplification. Nucl Acids Res 20:2380.[Free Full Text]

Donaldson GC, Ball LA, Axelrood PE, Glass NL., 1995 Primer sets developed to amplify conserved genes from filamentous ascomycetes are useful in differentiating Fusarium species associated with conifers. Appl Env Microbiol 61:1331-1340[Abstract]

Felsenstein J., 1978 Cases in which parsimony or compatibility methods will be positively misleading. Syst Zool 27:401-410

———. 1985 Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783-791

———. 1995 PHYLIP (Phylogeny Inference Package) Version 3.572c. Seattle: University of Washington

Ferreira AVB, Glass NL., 1996 PCR from fungal spores after microwave treatment. Fungal Genetics Newsletter 43:http://www.fgsc.net/fgn43/ferreir.html

Gargas A, Taylor JW., 1992 Polymerase chain reaction (PCR) primers for amplifying and sequencing 18S rDNA from lichenized fungi. Mycologia 84:589-592

Glass NL, Donaldson GC., 1995 Development of primer sets designed for use with the PCR to amplify conserved genes from filamentous ascomycetes. Appl Env Microbiol 61:1323-1330[Abstract]

Henk DA, Weir A, Blackwell M., 2003 Laboulbeniopsis termitarius, an ectoparasite of termites newly recognized as a member of the Laboulbeniomycetes. Mycologia 95:561-564[Abstract/Free Full Text]

Hojo M, Miura T, Maekawa K, Iwata R, Yamane A., 2001 Termitaria species (Termitariales, Deuteromycetes) found on Japanese termites (Isoptera). Sociobiology 38:327-342

Huelsenbeck JP, Ronquist F., 2001 MrBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17:754-755[Abstract/Free Full Text]

Khan SR, Aldrich HC., 1973 Conidiogenesis in Termitaria snyderi (Fungi Imperfecti). Can J Bot 51:2307-2314

———, ———. 1975 The haustoria of Termitaria snyderi. J Invert Path 25:247-260

———, Kimbrough JW., 1974a Morphology and development of Termitaria snyderi. Mycologia 66:446-462

———, ———. 1974b Taxonomic position of Termitaria and Mattirolella (entomogenous Deuteromycetes). Amer J Bot 61:395-399

Kimbrough JW, Lenz M., 1982 New species of Termitaria (Termitariales, Deuteromycetes) on Australian termites (Isoptera). Bot Gaz (Crawfordsville) 143:262-272

———, Thorne BL., 1982 Structure and development of Mattirolella crustosa (Termitariales, Deuteromycetes) on Panamanian termites. Mycologia 74:201-209

Lanyon SM., 1993 Phylogenetic frameworks—towards a firmer foundation for the comparative approach. Biol J Linn Soc 49:45-61

Malloch D, Blackwell M., 1990 Kathistes, a new genus of pleomorphic ascomycetes. Can J Bot 68:1712-1721

Swofford DL., 2001 PAUP* 4.0. Phylogenetic analysis using parsimony (*and other methods). Version 4d65. Sunderland, Massachusetts: Sinauer Associates Inc

Thompson JD, Higgins DG, Gibson TG., 1994 CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673-4680[Abstract/Free Full Text]

Thaxter R., 1920 Second note on certain peculiar fungus-parasites of living insects. Bot Gaz (Crawfordsville) 69:1-27

White TJ, Bruns T, Lee S, Taylor J., 1990 Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JH, White TJ, eds. PCR protocols. A guide to methods and applications. San Diego, California: Academic Press. p 315–322





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Blackwell, M.
Right arrow Articles by Jones, K. G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Blackwell, M.
Right arrow Articles by Jones, K. G.
Agricola
Right arrow Articles by Blackwell, M.
Right arrow Articles by Jones, K. G.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS