JWD Subscribe to eTOC alerts
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Journal of Wildlife Diseases, 41(4), 2005, pp. 701-706
© Wildlife Disease Association  2005
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
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 Feldman, S. H.
Right arrow Articles by Green, D. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Feldman, S. H.
Right arrow Articles by Green, D. E.

PHYLOGENETIC CLASSIFICATION OF THE FROG PATHOGEN AMPHIBIOTHECUM (DERMOSPORIDIUM) PENNERI BASED ON SMALL RIBOSOMAL SUBUNIT SEQUENCING

Sanford H. Feldman1,3, Jeffrey H. Wimsatt1 and D. Earl Green2

1 Center for Comparative Medicine, University of Virginia, Charlottesville, Virginia 22908, USA
2 National Wildlife Health Center, United States Geological Survey, Madison, Wisconsin 53711, USA

3 Corresponding author (email: shf2b{at}virginia.edu)


   ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
ABSTRACT:   We determined 1,600 base pairs of DNA sequence in the 18S small ribosomal subunit from two geographically distinct isolates of Dermosporidium penneri. Maximum likelihood and parsimony analysis of these sequences place D. penneri in the order Dermocystida of the class Mesomycetozoea. The 18S rRNA sequences from these two isolates only differ within a single region of 16 contiguous nucleotides. Based on the distant phylogenetic relationship of these organisms to Amphibiocystidium ranae and similarity to Sphaerothecum destruens we propose the organism be renamed Amphibiothecum penneri.
  Key words:  Amphibiothecum, dermocystida, Dermosporidium, frog diseases, mesomycetozoea.


   INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Species of the genus Dermosporidium are infrequently studied anuran parasites that were previously characterized based on morphologic observations using light microscopy. Dermosporidium hylae was the first species identified and was described from superficial cysts within the dorsal skin of the tree frog Hyla rubra in Brazil (Carini, 1940). Several decades later, other organisms of similar morphology, clinical appearance, and specificity for amphibians were described including Dermosporidium granulosum from Rana temporaria (Broz and Privora, 1952), Dermosporidium multigranulare from Rana esculenta (Broz and Kulda, 1954), and Dermosporidium penneri from the toad Bufo americanus (Jay and Pohley, 1981). However, the assignment of the genus Dermosporidium was recently challenged based on the lack of supportive data for the morphologic criteria used in the initial descriptions. A single genus "Amphibiocystidium" was proposed to encompass microorganisms previously defined in the genera Dermocycoides, Dermosporidium, and amphibian Dermocystidium (Pascolini et al., 2003). Moreover, morphologic and phylogenetic analysis of the small ribosomal subunit of Amphibiocystidium ranae revealed the taxon of this amphibian pathogen to be closely related to the Dermocystida infecting fish (Pereira et al., 2005). The sole use of morphologic criteria has led to some confusion regarding the classification of these microorganism(s). Taxonomic genomic characterization of closely related organisms including those in the genus Dermosporidium should help to clarify these key relationships. Based on morphologic similarities of microorganisms classified as members of the genus Dermosporidium to characteristics of members of the genus Dermocystidium, the former was assumed to belong to the class Mesomycetozoea (Mendoza et al., 2002). In our current study, we present molecular evidence of the relationship of D. penneri to members of the Order Dermocystida based on phylogenetic analysis of the DNA sequence in the gene encoding the 18S small subunit ribosomal. Based on our findings and the recommendations of Pereira et al. (2005) for changing the nomenclature of amphibian disease organisms in the class Dermocystida, we propose changing the name of Dermosporidium penneri to Amphibiothecum penneri.


   MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
We developed oligonucleotide primers unique to the class Mesomycetozoea by sequence alignment of the 18S rRNA subunit genes of Ichthyophonus hoferi (GenBank accession U25637; Ragan et al., 1996), Ichthyophonus irregularis (GenBank accession AF232303; Rand et al., 2000), Dermocystidium salmonis (GenBank accession U21337; Jostensen et al., 2002), Dermocystidium sp. (Gen-Bank accession U21336; Ragan et al., 1996), Rhinosporidium seeberi (GenBank accession AF118851; Herr et al., 1999), Pseudoperkinsus tapetis (GenBank accession AF192386; Figueras et al., 2000), and Bufo valliceps (GenBank accession M59386; Hedges et al., 1990). The primers were designed to amplify as much of the 18S rRNA gene as possible, recognizing members of the order Mesomycetozoea without concurrently amplifying 18S rRNA gene sequence found in toads. The primer sequences are 5'-TAAGCCATGCATGTCAAGTATAA (forward) and 5'-ACTAGGAATTCCTCGTBA AGATS (reverse).

Samples of A. penneri microdissected cysts from infected skin of Bufo americanus used in this study had been previously collected at New London Submarine Base, Connecticut (designated 18617-002) and Petersburg National Battlefield, Virginia (designated 4825-004) and had been stored frozen in ethanol. Histologically, A. penneri elicited a granulatous inflammatory cell response in which spores were extracellular, 10 µm in diameter, and contained approximately twenty 1–4 µm inclusions and a central nucleus (data not shown). A more comprehensive description of the morpohology of A. penneri isolates and the pathology elicited by its presence has been published elsewhere (Green and Converse, 2005). DNA was isolated using DNeasy tissue kit (Qiagen Corporation, Valencia, California, USA) according to the manufacturer’s recommendations for animal tissues, and sampling emphasized microdissected cysts within infected tissues. Five microliters of the isolated DNA was subjected to amplification in 50 µl polymerase chain reactions (PCRs) using Ex-Taq polymerase reagents (Takara Mirus Bio Corp., Madison, Wisconsin, USA) according to manufacturer’s recommendations and 0.5 µM of each primer. A 9600 gradient Robocycler (Strategene, La Jolla, California, USA) was set to the following PCR parameters: 95 C hot start for 90 sec, 35 cycles of 95 C (45 sec)/55 C (60 sec)/72 C (140 sec), and a final extension of 72 C for 5 min. Five microliters of each reaction was analyzed by agarose gel electrophoresis. The resulting 1.6 kB amplicons were ligated into TOPO pCR-4 (Invitrogen Corp., Carlsbad, California, USA) and used to transform chemically competent Top 10 E. coli (Invitrogen) that were plated onto LB agar containing 100 µg/ml ampicillin. Ten colonies from each sample were selected and grown separately overnight in LB-ampicillin broth. The plasmids were isolated using a Plasmid minikit (Qiagen), eluted with water, and analyzed by EcoRI restiction digest (New England Biolabs, Waverly, Massachusetts, USA). Six plasmids containing the amplicon insert from each isolate were submitted for DNA sequencing to Davis Sequencing (Davis, California, USA) in the forward (from the T3 site) and reverse (from the T7 site) directions. The entire procedure of amplification, cloning, and sequence analysis was repeated once with the same sample DNAs and the sequence results combined with those determined initially for unambiguous characterization of the amplified DNA sequence.

Alignment of all forward sequences for each geographic isolate, and separately for all reverse sequences, was performed using Omiga 1.1 (Oxford Molecular Ltd., Oxford, UK). Separate consensus sequences were determined for the Connecticut and the Virginia isolates, and the two consensus sequences aligned for comparison. Phylogenetic analysis of the consensus sequences used the same GenBank sequences used for primer development plus Amphibiocystidium ranae (GenBank accession AY550245 and AY692319; Pereira et al., 2005), Dermocystidium percae (GenBank accession AF533941; Pascolini et al., 2003), Amoebidium parasiticum (GenBank accession Y19155; Ustinova et al., 2000), Anurofeca richardsi (GenBank accession AF070445; Baker et al., 1999), Ichthyophonida gen. sp. LKM51 (GenBank accession AJ130859; Van Hannen et al., 1999), Psorospermium haeckeli (GenBank accession U33180; Ragan et al., 1996), Sphaeroforma arctica (GenBank accession Y16260; Jostensen et al., 2002), Sphaerothecum destruens (GenBank accession AY267346; Arkush et al., 2003), Capsaspora owczarzaki (GenBank accession AF436886; Hartel et al., 2002). Microciona prolifera (GenBank accession L10825; Wainright et al., 1993), and Mnemiopsis leidyi (GenBank accession L10826; Wainright et al., 1993) were used to root the phylograms in the phylogenetic analyses. All GenBank sequences were converted to GCG sequence and truncated at the 5' and 3' ends to coincide with the sequences determined for the A. penneri amplicons. All 20 sequences were aligned and analyzed by parsimony analysis and again by maximum likelihood analysis using PHYLIP inference software (Felsenstein, 2004).

The significance of branch points (nodes) for the consensus phylogram generated was determined by bootstrap analysis of 1,000 datasets, each dataset searched for the most parsimonious tree. Bootstrap maximum likelihood analysis of the same 1,000 datasets was performed (with global rearrangement, the rough option turned off, and with the jumble set and random number seed again set to three). A single phylogram is presented using TREEVIEW (Page, 1996) that represents a consensus of trees from maximum likelihood analysis showing the number of trees in agreement at each node (Fig. 1Go). Mnemiopsis leidyi and M. prolifera were used together as an outgroup to root the consensus tree. Phylograms were generated using maximum likelihood analysis with and without imposing a molecular clock (i.e., the tips of all branches were made equidistant from the root). The log likelihood of each maximum likelihood method was analyzed by Chi Square goodness-of-fit testing against the assumption that the small ribosomal subunit gene can function as a molecular clock.


Figure 1
View larger version (29K):
[in this window]
[in a new window]
 
FIGURE 1. This phylogram from maximum likelihood analysis incorporates a molecular clock assumption. The scale shows the branch length equivalent to the occurrence of substitution of 10% of the nucleotide sequence. The relatedness of two microorganisms in this phylogram is the sum of the length of all horizontal branches separating the two species (neglect vertical distances between species). The number to the left at each node is the consensus among 1,000 trees by bootstrap maximum likelihood analysis. This maximum likelihood model allows changes to occur multiple times at any single nucleotide site, thus a branch twice as long as another may not necessarily contain twice as many nucleotides differences. The phylogram shows definitively that A. penneri is a member of the order Dermocystida.

 

   RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Amplified 18S rRNA sequences of 18617-002 (Connecticut) and 4825-004 (Virginia) were 1,607 and 1,599 bases in length, respectively, and have been deposited in GenBank (accession numbers AY772000 and AY772001). Alignment of sequences from the two A. penneri isolates shows complete sequence agreement with the exception of a single contiguous DNA region. This regional difference spans 16 nucleotides wherein eight nucleotides were deleted in the 18S rRNA gene of the Virginia isolate and two nucleotide substitutions occurred among the remaining eight nucleotides (Fig. 2Go). The clones from both isolates showed the same variability (A, C, or G) at nucleotide position 1,590 that may have occurred as an artifact resulting from degeneracy in the design of the reverse primer at that position.


Figure 2
View larger version (6K):
[in this window]
[in a new window]
 
FIGURE 2. The small ribosomal subunit sequences for samples 4825-004 and 18617-002 when aligned demonstrate complete consensus, with the exception one region of 16 nucleotides shown here between positions 146 and 162. Sample 4825-004 (Virginia) lacks eight nucleotides and contains two nucleotide substitutions within this region.

 
The topology of the consensus phylograms generated by bootstrap analyses using parsimony and maximum likelihood methods was the same (data not shown). Using the maximum likelihood method to compare phylograms within the Dermocystida clade, all nodes are well supported (greater than 50% consensus among tree topologies) except at the following nodes: the node distinguishing Rhinosporidium seeberi and Dermocystidium sp.; distinguishing A. penneri from S. destruens, the node distinguishing D. salmonis from A. ranae, and the node distinguishing C. owczarzaki from P. haeckeli. This analysis showed that A. penneri is a member of the order Dermocystida within the class Mesomycetozoea and further suggests that the two A. penneri isolates analyzed are distinct from each other. The phylogram generated is identical with others published previously for the class Mesomycetozoea (Wainright et al., 1993; Ragan et al., 1996; Herr et al., 1999; Hertel et al., 2002; Mendoza et al., 2002; Pekkarinen et al., 2003) and similar to the phylogram generated by Pereira et al. (2005). The rosette agent, which is included in many of these phylograms was recently renamed Sphaerothecum destruens (Arkus et al., 2003).

Phylograms generated using maximum likelihood methods, with and without a molecular clock assumption, were of the same (unrooted) topology (data not shown). Imposing the molecular clock (natural log likelihood=–8769) when compared with no clock assumption (natural log likelihood=–8679) during analysis of the 18S small ribosomal subunit sequence was not statistically supported ({chi}2=180, d.f.=18, P<0.005). Thus, maximum likelihood tree analysis with molecular clock assumption does not provide a relative estimate of the divergence of species within the class Mesomycetozoea.


   DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The number of parasitic and saprophytic microorganisms assigned to the monophylitic class Mesomycetozoea is continuing to increase (Ragan et al., 1996; Mendoza et al., 2002; Pereira et al., 2005). The 18S rRNA sequence analysis presented in our study supports the categorization of A. penneri as an organism closely related to the Dermocystidia parasites of fish, as originally deduced using morphologic characters (Jay and Pohley, 1981). Interestingly, A. penneri forms a clade that includes S. destruens (since the node separating these is not well supported) presenting a taxonomic quandry given the recent suggestion to use the epithet Dermotheca for organisms with phylogenetic affinity to S. destruens (Pereira et al., 2005). The morphologic character of A. penneri cysts is round to oval 4–6 mm in diameter located in divided chambers in the ventral subcutis of the toad B. americanus (Jay and Pohley, 1981), unlike the giant elongate-cylindrical intracellular cysts found in liver cells of salmonids characteristic of S. destruens (Arkush et al., 2003). The spores of A. penneri are 10–12 µm in diameter, whereas the spores of S. destruens are 2–6 µm. Thus, the phylogenetic similarity of A. penneri and S. destruens in the present study is not reflected in host species, or cyst location and morphology. Therefore, we propose assigning a new genus name, Amphibiothecum, which reflects both the host species and the phylogenetic similarity to S. destruens. The topology of our phylogram is similar to that reported by Pereira et al. (2005), but the node distinguishing A. ranae from D. salmonis is not well supported in our study whereas the node distinguishing the two A. ranae isolates is very well supported. Although our study, as well as studies by others (Pereira et al., 2005), demonstrates that the 18S sequence of D. percae is (phylogenetically) more rudimentary when compared with the same region of gene sequence of other members in the class Dermocystida, we have chosen to not address altering the genus designation in this report.

Amphibian population declines and species extinctions have been occurring worldwide over the past 25 yr (Green and Sherman, 2001; Green et al., 2002). In the United States amphibian population declines associated with emerging infectious diseases such as Batrachochytrium dendrobatidis (chytrid), iridovirus, and a dermocystidium-like infection (Green et al., 2002; 2003) have adversely affected natural populations. Although infection with A. penneri has not been shown to contribute directly to amphibian die-offs (Green et al., 2002), its importance may increase as ecologic stressors change over time. The 18S rRNA sequence determined for Amphibiothecum penneri isolates in the present study can be used to develop PCR protocols that are specific for this micro-organism; this will allow for rapid, sensitive, and specific identification of this organism from amphibian specimens. These sequences also provide much needed information to better understand the taxonomy of the Mesomycetozoea.


   ACKNOWLEDGMENTS
 
This research was funded through a cooperative agreement with the United States Geological Survey, National Wildlife Health Center.


   LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
ARKUSH, K. D., L. MENDOZA, M. A. ADKISON, AND R. P. HEDRICK. 2003. Observations on the life stages of Sphaerothecum destruens n.g., n. sp., a mesomycetozoean fish pathogen formally referred to as the Rosette agent. Journal of Eukaryotic Microbiology 50: 430–438.[Medline]

BAKER, G. C., T. J. BEEBEE, AND M. A. RAGAN. 1999. Prototheca richardsi, a pathogen of anuran larvae, is related to a clade of protistan parasites near the animal-fungal divergence. Microbiology 145: 1777–1784.[Abstract]

BROZ, O., AND M. PRIVORA. 1952. Two skin parasites of Rana temporaria: Dermocystidium ranae Guyenot & Naville and Dermosporidium granulosum n. sp. Parasitology 42: 65–69.[Medline]

———, AND J. KULDA. 1954. Dermosporidium multigranulare n. sp. parasit z kuze Rana esculenta. Acta Societatis Zoologicae Bohemicae 18: 91–97.

CARINI, A. 1940. Sobre um parasito semelhante ao "Rhinosporidium," encontrado em quistos da pele de uma "Hyla." Arquivos de Institute Biologie Sao Paulo. 11:93–98.

FELSENSTEIN, J. 2004. PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by the author. Department of Genome Sciences, University of Washington, Seattle.

FIGUERAS, A., G. LORENZO, M. C. ORDAS, M. GOUY, AND B. NOVOA. 2000. Sequence of the small subunit ribosomal DNA gene of Pseudoperkinsus tapetis isolated from carpet shell clam in Galicia (NW of Spain). Marine Biotechnology 2: 419–428.

GREEN, D. E., AND K. A. CONVERSE. 2005. Diseases of frogs and toads. In Wildlife diseases: Landscape epidemiology, spatial distribution, and utilization of remote sensing technology, S. K. Majumdar, J. Huffman, and F. J. Brenner, and A. I. Panah (eds.). Pennsylvania Academy of Science, Philadelphia, pp. 89–117.

———, AND C. K. SHERMAN. 2001. Diagnostic histological findings in the Yosemite toads (Bufo canorus) from a die-off in the 1970s. Journal of Herpetology 35: 92–103.

———, K. A. CONVERSE, AND A. K. SCHRADER. 2002. Epizootiology of sixty-four amphibian morbidity and mortality events in the USA, 1996–2001. Annals New York Acadamy of Sciences 969: 323–339.

———, S. H. FELDMAN, AND J. H. WIMSATT. 2003. Emergence of a Perkinsus-like agent in anuran liver during die-offs of local populations: PCR detection and phylogenetic characterization. American Association of Zoo Veterinarians Conference Proceedings. Milwaukee, Wisconsin.

HEDGES, S. B., K. D. MOBERG, AND L. R. MAXSON. 1990. Tetrapod phylogeny inferred from 18S and 28S ribosomal RNA sequences and a review of the evidence for amniote relationships. Molecular Biology of Evolution 7: 607–633.

HERR, R. A., L. AJELLO, J. W. TAYLOR, S. N. ARSECULERATNE, AND L. MENDOZA. 1999. Phylogenetic analysis of Rhinosporidium seeberi’s 18S small subunit ribosomal DNA groups this pathogen among members of the protoctistan Mesomycetozoa clade. Journal of Clinical Microbiology 37: 2750–2754.[Abstract/Free Full Text]

HERTEL, L. A., C. J. BAYNE, AND E. S. LOKER. 2002. The symbiont Capasapora owczarzaki, nov. gen. Nov. sp., isolated from three strains of the pulmonate snail Biomphalaria glabrata is related to members of the Mesomycetozoea. International Journal for Parasitology 32: 1183–1191.[Medline]

JAY, J. M., AND W. J. POHLEY. 1981. Dermosporidium penneri sp. n. from the skin of the American toad, Bufo americanus (Amphibia: Bufonidae). Journal of Parasitology 67: 108–110.

JOSTENSEN, J.-P., S. SPERSTAD, S. JOHANSEN, AND B. LANDFALD. 2002. Molecular-phylogenetic, structural and biochemical features of a cold-adapted, marine ichthyosporean near the animal-fungal divergence, described from in vitro cultures. European Journal of Protistology. 38: 93–104.

MENDOZA, L., J. W. TAYLOR, AND L. AJELLO. 2002. The class Mesomycetozoea: a heterogeneous group of microorganisms at the animal-fungal boundary. Annual Reviews of Microbiology. 56: 315–344.

PAGE, R. D. M. 1996. TREEVIEW: An application to display phylogenetic trees on personal computers. Computer Applications in the Biological Sciences 12: 357–358.

PASCOLINI, R., P. DASZAK, A. A. CUNNINGHAM, S. TEI, D. VAGNETTI, S. BUCCI, A. FAGOTTI, AND I. DI ROSA. 2003. Parasitism by Dermocystidium ranae in a population of Rana esculenta complex in Central Italy and description of Amphibiocystidium n. gen. Diseases of Aquatic Organisms 56: 65–74.[Medline]

PEKKARINEN, M., J. LOM, C. A. MURPHY, M. A. RAGA, AND I. L. DYKOVA. 2003. Phylogenetic position and ultrastructure of two Dermocystidium species (Ichthyosporea) from the common perch (Perca fluviatilis). Acta Protozoologica 42: 287–307.

PEREIRA, C. N., I. DI ROSA, A. FAGOTTI, F. SIMONCELLI, R. PASCOLINI, AND L. MENDOZA. 2005. The pathogen of frogs Amphibiocystidium ranae is a member of the order dermocystida in the class Mesomycetozoea. Journal of Clinical Microbiology 34: 192–198.

RAGAN, M. A., C. L. GOGGIN, R. H. CAWTHORN, L. CERENIUS, W. V. C. JAMIESON, S. M. PLOURDE, T. G. RAND, K. SODERHALL, AND R. R. GUTELL. 1996. A novel clade of protistan parasites near the animal-fungal divergence. Proceedings of the National Academy of Sciences United States of America 93: 11907–11912.

RAND, T. G., K. WHITE, J. J. CANNONE, R. R. GUTELL, C. A. MURPHY, AND M. A. RAGAN. 2000. Ichthyophonus irregularis sp. nov. from the yellowtail flounder Limanda ferruginea from the Nova Scotia shelf. Diseases of Aquatic Organisms 41: 31–36.[Medline]

USTINOVA, I., L. KRIENITZ, AND V. A. R. HUSS. 2000. Hyaloraphidium curvatum is not a green alga but is a lower fungus; Amoebidium parasiticum is not a fungus, but a member of the DRIPs. Protist 151: 253–262.[Medline]

VAN HANNEN, E. J., W. MOOIJ, M. P. VAN AGTERVELD, H. J. GONS, AND H. J. LAANBROEK. 1999. Detritus-dependent development of the microbial community in an experimental system: qualitative analysis by denaturing gradient gel electrophoresis. Applied and Environmental. Microbiology 65: 2478–2484.

WAINRIGHT, P. O., G. HINKLE, M. L. SOGIN, AND S. K. STICKEL. 1993. Monophyletic origins of the metazoa: An evolutionary link with fungi. Science 260: 340–342.[Abstract/Free Full Text]

Received for publication 8 November 2004.




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
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 Feldman, S. H.
Right arrow Articles by Green, D. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Feldman, S. H.
Right arrow Articles by Green, D. E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS