DNA barcoding was used to investigate laboratory and commercial stocks of Eisenia species from four provinces of South Africa. The COI gene was partially amplified and sequenced in selected earthworms from eight local populations (focal groups) and two European laboratory stocks (non-focal groups). Only nine COI haplotypes were identified from the 224 sequences generated. One of these haplotypes was found to belong to the megascolecid Perionyx excavatus. The remaining eight haplotypes belonged to the genus Eisenia, although only a single E. fetida haplotype, represented by six specimens, was found in one of the European populations. The other seven haplotypes, all occurring in South Africa, were E. andrei. One of the commercial stocks from South Africa and a laboratory culture from Europe were mixtures of E. andrei-P. excavatus and E. andrei-E. fetida, respectively. Previous allozyme studies have helped to suggest that some of the populations included in this study may be suffering from inbreeding depression, which could result in adverse consequences for both the vermiculture industry and ecotoxicological research in South Africa.
Eisenia fetida (Savigny, 1826) and Eisenia andrei Bouché, 1972 have become cosmopolitan earthworm species because of their worldwide use in ecotoxicological testing and vermicomposting. Originating from Palaearctic Europe, they have been successfully introduced to other ecozones mainly because of their wide temperature tolerance and robustness (Hendrix et al. 2008). Both E. fetida and E. andrei are the earthworm species recommended by the Organisation for Economic Co-operation and Development (OECD 1984, 2004) and the International Organization for Standardization (ISO 2008, 2012) for the testing of chemicals.
Historically, Savigny only described E. fetida (.E. foetida), which was later suspected of harbouring a cryptic sister species. Bouché (1972) divided E. fetida into two subspecies, E. foetida foetida (current E. fetida) and E. foetida unicolour (current E. andrei). Using allozyme polymorphism, Jaenicke (1982) and Øien and Stenersen (1984) indicated that these subspecies are different species. Their findings were supported by Domínguez et al. (2005) and Pérez-Losada et al. (2005), who concluded that E. fetida and E. andrei are different biological and phylogenetic species, as judged by their reproductive isolation and DNA divergence.
In South Africa and world-wide, E. andrei and E. fetida are used in the vermiculture industry and scientific research. Two South African research laboratories in the field of terrestrial ecotoxicology, at Stellenbosch University and North-West University, respectively, have used E. fetida and E. andrei for decades, and the output of their research has been published in the local and international scientific literature (Reinecke & Viljoen 1991; Reinecke & Reinecke 1997; Prinsloo et al. 1999; Reinecke et al. 2001; Reinecke et al. 2002; Maboeta & van Rensburg 2003a, b; Maboeta et al. 2008; Owojori et al. 2009; Voua Otomo & Reinecke 2010).
Despite the interest in both species, there has been no molecular study of populations introduced locally or of cultures of E. andrei and E. fetida used in South Africa. Molecular work on selected laboratory and field populations have focused on the toxicological effects of particular toxicants on DNA integrity and allozyme polymorphism in these earthworms (Reinecke & Reinecke 2004; Voua Otomo et al. 2011). Voua Otomo et al. (2009) conducted a DNA barcoding study on an Eisenia sp. laboratory stock housed in the Zoology Department of Stellenbosch University as a means of researching its taxonomic identity.
The need for molecular studies on these earthworms is critical for several reasons. Being economically and scientifically important, basic information such as species identity and the genetic differentiation between Eisenia spp. stocks should be relevant to the breeders, potential buyers and researchers alike. The ecotoxicological literature reveals that countless researchers worldwide rely upon informally identified commercial earthworm stocks for laboratory bioassays (e.g., Beyer 1996; Fitzpatrick et al. 1996; Saint-Denis et al. 1998; Krauss et al. 2000; Gevao et al. 2001; Miyazaki et al. 2002; Gambi et al. 2007; Lin et al. 2010).
Moreover, earthworm cultures kept isolated for many generations may, with time, suffer from inbreeding depression characterized by low heterozygosity (Voua Otomo et al. 2011). This may undermine sustainable earthworm breeding and quality research.
The aim of this study was to conduct a DNA barcode investigation of earthworm stocks from selected vermiculture establishments and research laboratories in South Africa in order to confirm their taxonomic status, and assess their levels of genetic richness and differentiation.
MATERIAL AND METHODS
In the present study, the term “population” is used in an inclusive manner and thus may refer to a free-living “wild” population or to a captive breeding stock. A total of eight focal and two non-focal populations were the subject of this study. Focal populations included two vermiculture stocks from Johannesburg (Gauteng, South Africa), two vermiculture stocks and a laboratory culture from Potchefstroom (North West, South Africa), a free-living population and a laboratory culture from Stellenbosch (Western Cape, South Africa) and a vermiculture stock from Port Elizabeth (Eastern Cape, South Africa). Two non-focal laboratory cultures were acquired from Brno (Czech Republic) and Southampton (UK). Table 1 provides the presumed identities (as given by the owners) of the respective earthworm groups, their geographical locality, their function/use and, when applicable, an excerpt from the list of recent publications based upon research work carried out on the populations concerned.
Because of the economic importance of Eisenia spp. and considering that potential earthworm buyers are mostly unable to distinguish between different earthworm species, we decided not to sort the randomly picked local specimens according to phenotypic features, thus allowing us to identify possible mixed cultures.
Localities, presumed identity, use and publications record for the earthworm groups included in the present study. Abbreviations: n — the number of specimens used for COI genotyping from the respective groups, ANs — the Genbank accession numbers for the respective sequences.
Total genomic DNA was extracted from 224 earthworms using the NucleoSpin® Tissue kit (Macherey-Nagel). Samples of five to ten milligrams of the tail section of the selected specimens were treated according to the manufacturer's instructions. The universal primers LCO1490 and HCO2198 (Folmer et al. 1994) were used to amplify 683 bp of the cytochrome oxidase I (COI) gene.
PCR reactions consisted of 0.3 µl (∼30 ng) DNA template, 12.5 µl PCR Master Mix (Fermentas), 11 µl nuclease-free water (Fermentas) and 10 pmol (∼1 µl) of each of the primers. PCR cycling comprised an initial denaturation step at 94 °C for 5 min followed by 35 cycles at 94 °C for 30 s, 50 °C for 30 s and 72 °C for 45 s. A final extension step at 72 °C for 5 min completed the reactions. Successful amplification was verified by electrophoretic means using agarose gels (0.75 g SeaKem® LE Agarose, Lonza, in 50 ml TAE buffer, 1.5% (w/v) stained with 5 µl ethidium bromide). Sequencing reactions were performed using the ABI v3.1 BigDye® kit. Purified sequences were run on an ABI 3500XL Genetic Analyser.
All the barcodes generated in the present study were deposited in GenBank (Table 1). They were tentatively identified using the BOLD (Barcode of Life Data Systems) Identification System and compared to published COI sequences of E. andrei, E. fetida and Allolobophoridella eiseni deposited in GenBank by Pérez-Losada et al. (2005).
All the sequences were aligned, edited and analysed in MEGA v5 (Tamura et al. 2011) using the Kimura-2-parameter (K2P) method (Kimura 1980). A neighbour-joining tree was subsequently constructed. Bootstrap support was obtained from 1000 iterations. Since COI diversity is highly dependent on effective population size and because of the uneven sample sizes of the groups included in this study, we used the Contrib software of Petit et al. (1998) to assess haplotypic richness and diversity contribution after rarefaction. The software package NETWORK 18.104.22.168 (Fluxus Technology Ltd) was used to compute a haplotype network of the distinct Eisenia spp. COI sequences occurring in South Africa, using the Median-joining method.
Nine distinct sequences of the COI gene were identified amongst the 224 worms included in this study. The haplotype distributions across the populations revealed that H1 (haplotype 1) was the most widespread and H2 the most frequent, representing more than 70 % of all the COI sequences (Table 2). Five haplotypes were unique to their population of origin, viz. H4 (JNB; Johannesburg), H6 (SUN; Stellenbosch University), H7 (NWU; North-West University), H8 (PE; Port Elizabeth) and H9 (ENG; Southampton).
The analysis of all the haplotypes together with previously published COI sequences of E. andrei and E. fetida revealed that the nine distinct sequences of COI identified in the ten groups could represent four different earthworm species. Haplotypes HI to H6 grouped with previously identified sequences of E. andrei (K2P ≤8.28%) (Fig. 1). H7 grouped with BOLD sequences identified as E. andrei. However, K2P distances revealed that sequence divergence between H7 and the other E. andrei haplotypes was as high as 31.10%. The identity of H7 is therefore uncertain, especially considering the fact that it grouped with unpublished (i.e. potentially unverified), alleged E. andrei sequences from BOLD (EWSJC613-10, EWSJC614-10) (Fig. 1). H8 grouped with GenBank sequences of the megascolecid Perionyx excavatus (K2P ≤1.2%). The BOLD system also identified H8 as P. excavatus. H9 grouped with previously identified sequences of E. fetida (K2P ≤11.7%). The earthworm cultures from Port Elizabeth and Southampton were mixtures of E. andrei - P. excavatus and E. andrei - E. fetida, respectively.
Genetic richness and differentiation of local populations
Eight of the nine COI haplotypes (H1-H8) occurred in the selected South African earthworm stocks. H8, as established above, does not belong to the genus Eisenia. Consequently, only seven Eisenia COI haplotypes were found to occur in local populations. All of these, with the exception of H7, grouped with conclusively identified specimens of E. andrei. Table 3 provides the Kimura 2-parameter distance matrix between these haplotypes. Prior to rarefaction analyses, PE (H = 0.426), JOZ (H = 0.378) and NWU (H = 0.377) had, in order, the three highest haplotype diversities (Table 4). After rarefaction to a common sample size of 10, this order changed to JOZ (r(10) = 2), PE (r(10) = 1.658) and JNB (r(10) = 1.5). MID also contributed more to the total genetic diversity amongst populations (HT= 0.4498), as indicated by the only positive CT (CT = 0.322), which was mostly due to the strong divergence (CD = 0.38) of MID from the other populations (Table 4). Differentiation indices DHT and DGST >0.75 for MID revealed that this population was indeed the most divergent of the local populations included in the present study. Negative CD values for the other populations reflected a lack of significant differentiation between them. This was confirmed by conventional population pairwise FSTs that showed non-significant differentiation amongst these populations.
Figure 2 represents a network of the E. andrei haplotypes found in local South African populations. The dubious haplotype H7 was excluded from this analysis. The least number of mutations found was between H1 and H2 (a single mutation) and the highest number of mutations was between H2 and H4 (31 mutations).
Haplotype distribution and frequency across all the populations investigated. H2 was the most frequent haplotype, representing more than 70 % of all the COI sequences.
DNA analysis reveals that the sequences generated from South African-based Eisenia populations grouped unequivocally with known sequences of E. andrei. Earthworm breeders and researchers have assumed that these local groups represent cultures and populations of E. fetida. Reinecke and Viljoen (1991) stated that local Eisenia populations could be a mixture of E. andrei and E. fetida. To date, no locally occurring E. fetida specimen has been formally identified using DNA markers. The occurrence of mixed local populations of E. andrei and E. fetida cannot be excluded as it is acknowledged that both species commonly occur in mixed colonies and that E. andrei could outcompete E. fetida during periods of food abundance (Elvira et al. 1996). Domínguez et al. (2005) noted that E. andrei is the predominant species in commercial vermiculture establishments, while E. fetida is mostly found in free-living populations. Considering that seven out the eight local earthworm groups investigated were bred in captivity, perhaps the inclusion of more field populations would have helped to detect the presence of E. fetida.
Kimura 2-parameter distance matrix (%) between the E. andrei COI haplotypes (H1-H7) found in the studied South African earthworm groups. Distances between H7 and the rest of the haplotypes vary between 23.28 and 31.10%; the identity of H7 remains uncertain.
The vermiculture stock from Port Elizabeth was a mixture of E. andrei and P. excavatus; the oriental compost worm known to be able to reproduce parthenogenetically and to thrive in similar living conditions as E. andrei and E. fetida (Hallatt et al. 1990). These results suggest that the untrained buyer seeking to purchase E. fetida in South Africa has a greater likelihood of acquiring E. andrei; and occasionally together with individuals of another species such as P. excavatus.
The unique COI sequence (H7) identified as an E. andrei sequence through the BOLD system was extremely divergent from the other E. andrei sequences. Using the K2P method, the accepted threshold for species delimitation on the basis of DNA barcode data is 15% K2P (Chang & James 2011). The divergence between H7 and the other E. andrei haplotypes was consistently more than 23 % K2P. An increasing number of cryptic oligochaete species have been reported in the literature since the recent advent of earthworm molecular studies (King et al. 2008; Pérez-Losada et al. 2009; Blakemore et al. 2010; James et al. 2010; Novo et al. 2010). H7 could represent an as yet undescribed species. However, additional molecular and morphological investigations would be required to shed further light on the matter.
COI haplotype numbers were limited to two or three distinct sequences within each of the local groups. This translated into a haplotype diversity (H) lower than 0.45 in all the populations. When compared to other such molecular studies in which COI polymorphism in earthworms has been investigated, the present haplotype diversity is proportionally very low. King et al. (2008) sequenced the COI gene in selected lineages of the European earthworm Allolobophora chlorotica and found H values as high as 0.95. Similarly, Novo et al. (2009) obtained H values as high as 0.92 in populations of the hormogastrid earthworm Hormogaster elisae from the central Iberian Peninsula. Equally high haplotypic richness has been reported in several other species of earthworms such as Dendrobaena octaedra (Cameron et al. 2008; Knott & Haimi 2010), Amynthas wulinensis (Chang et al. 2007), Aporrectodea rosea, Octolasion lacteum, and Lumbricus rubellus (Klarica et al. 2012).
Moreover, laboratory and vermicomposting cultures are susceptible to the founder effect (Mayr 1942) as they are usually started with a limited number of individuals. This may explain the comparatively poor haplotype diversity observed in South African E. andrei stocks. For Eisenia spp., the phenomenon could be compounded by the fact that known habitats of these species (compost heaps, manure, rich soils, etc.) are naturally fragmented. Despite their status as standard laboratory test species, molecular studies of free-living E. andrei and E. fetida are rare. The population genetics of these species has yet to be thoroughly investigated in Europe, where they originated.
Measure of genetic diversity and divergence for each South African population of Eisenia andrei based on COI sequence data after rarefaction to a common sample size of ten. Abbreviations: n — number of specimens included per population; Nb Hap. — number of haplotypes; H (SE) — haplotype diversity with standard error in brackets; π (SE) — nucleotide diversity with standard error in brackets; r (10) — allelic richness after rarefaction to a common size of ten specimens per sample; DHS, DHT, DGST — divergence indices from the other populations; CT, CS, CD — contribution indices to total diversity; CrT, CrS, CrD — contribution indices to total allelic richness (see Petit et al. (1998) for more details).
Being a species introduced to South Africa, E. andrei also suffered the effects of another significant factor upon being brought into the country; the propagule pressure, which stipulates that species introduced in large and consistent quantities are more likely to persist in their new environment compared to those introduced in limited numbers and involving relatively few release events (Lockwood et al. 2005). This particular factor may also help to explain the local predominance of E. andrei over E. fetida by assuming that larger and more consistent introduction events may have occurred for E. andrei.
Of all the local groups investigated, MID was the only significantly divergent population. The haplotype distributions across the populations (Table 2) show that MID was the only population not harbouring H2, the haplotype which represented 75 % or more of the COI sequences within the local populations. This perhaps indicates that H2 is rare in free-living populations of E. andrei or that this particular haplotype is selected against under relatively harsh environmental conditions.
Finally, Voua Otomo et al. (2011) established, using allozyme polymorphism, that the mean observed heterozygosity per locus (Ho) in two of the earthworm groups investigated in this study (SUN and MID - previously thought to be E. fetida) was zero. It is suspected that inbreeding could be occurring in these populations.
This may have significant implications for both the research sector and the vermiculturing industry. The SUN and MID groups have for instance been used in ecotoxicological research (Table 1). If the genetic diversity of laboratory populations is drastically reduced, the reliability of results from laboratory testing could be compromised. The lack of genetic variation has been associated with decreased fitness, often affecting traits such as growth, reproduction and survival (Charlesworth & Charlesworth 1987; Reed & Frankham 2003). Velando et al. (2006) researched the deleterious effects of inbreeding on the reproduction of E. andrei and reported that inbreeding causes a “strong reduction of cocoon production”
The use of DNA barcoding has helped to show that E. fetida may be rarer in South Africa than previously assumed. E. andrei is the main species used in both the vermiculture industry and laboratory research. Most of these captive stocks are genetically homogenous and may in some instances suffer from inbreeding depression.
This work is based upon research supported financially by the National Research Foundation of South Africa. The authors wish to acknowledge A.J. and S.A. Reinecke, P. Theron, H. Bouwman, K. Reid and G. Heron for providing some of the materials used in this study. The authors also wish to acknowledge Dr S. James and an anonymous reviewer for comments on the manuscript.
- W.N. Beyer 1996. Accumulation of chlorinated benzenes in earthworms. Bulletin of Environmental Contamination and Toxicology 57: 729–736. Google Scholar
- R. Blakemore , E. Kupriyanova & M. Grygier 2010. Neotypification of Drawida hattamimizu Hatai, 1930 (Annelida, Oligochaeta, Megadrili, Moniligastridae) as a model linking mtDNA (COI) sequences to an earthworm type, with a response to the ‘can of worms’ theory of cryptic species. ZooKeys 41: 1–29. Google Scholar
- M.B. Bouché 1972. Lombriciens de France: écologie et systématique. Paris: Institut National des Recherche Agronomique. Google Scholar
- E.K. Cameron , E.M. Bayne & D.W. Coltman 2008. Genetic structure of invasive earthworms Dendrobaena octaedra in the boreal forest of Alberta: insights into introduction mechanisms. Molecular Ecology 17: 1189–1197. Google Scholar
- C.-H. Chang & S. James 2011. A critique of earthworm molecular phylogenetics. Pedobiologia 54(suppl.): S3–S9. Google Scholar
- C.-H. Chang , Y.-H. Lin , I.-H. Chen , S.-C. Chuang & J.-H. Chen 2007. Taxonomic re-evaluation of the Taiwanese montane earthworm Amynthas wulinensis Tsai, Shen & Tsai, 2001 (Oligochaeta: Megascolecidae): polytypic species or species complex? Organisms Diversity and Evolution 7: 231–240. Google Scholar
- D. Charlesworth & B. Charlesworth 1987. Inbreeding depression and its evolutionary consequences. Annual Review of Ecology and Systematics 18: 237–268. Google Scholar
- J. Domínguez , A. Velando & A. Ferreiro 2005. Are Eisenia fetida (Savigny, 1826) and Eisenia andrei Boulé, 1972 (Oligochaeta, Lumbricidae) different biological species? Pedobiologia 49: 81–87. Google Scholar
- C. Elvira , J. Domínguez & M.J.I. Briones 1996. Growth and reproduction of Eisenia andrei and Eisenia fetida (Oligochaeta, Lumbricidae) in different organic residues. Pedobiologia 40: 377–384. Google Scholar
- L.C. Fitzpatrick , J.F. Muratti-Ortiz , B.J. Venables & A.J. Goven 1996. Comparative toxicity in earthworms Eisenia fetida and Lumbricus terrestris exposed to cadmium nitrate using artificial soil and filter paper protocols. Bulletin of Environmental Contamination and Toxicology 57: 63–68. Google Scholar
- O. Folmer , M. Black , W. Hoeh , R. Lutz & R. Vrijenhoek 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology 3: 294–299. Google Scholar
- N. Gambi , A. Pasteris & E. Fabbri 2007. Acetylcholinesterase activity in the earthworm Eisenia andrei at different conditions of carbaryl exposure. Comparative Biochemistry and Physiology Part C: Toxicology and Pharmacology 145: 678–685. Google Scholar
- B. Gevao , C. Mordaunt , K.T. Semple , T.G. Piearce & K.C. Jones 2001. Bioavailability of nonextractable (bound) pesticide residues to earthworms. Environmental Science and Technology 35: 501–507. Google Scholar
- L. Hallatt , A.J. Reinecke & S.A. Viljoen 1990. Life cycle of the oriental compost worm Perionyx excavatus (Oligochaeta). South African Journal of Zoology 25: 41–45. Google Scholar
- PF. Hendrix , M.A. Callaham , J.M. Drake , C.-Y Huang , S.W. James , B.A. Snyder & W. Zhang 2008. Pandora's box contained bait: the global problem of introduced earthworms. Annual Review of Ecology, Evolution, and Systematics 39: 593–613. Google Scholar
- ISO. 2008. ISO 17 512–1. Soil quality - Avoidance test for determining the quality of soils and effects of chemicals on behaviour - Part 1: Test with earthworms (Eisenia fetida and Eisenia andrei). Geneva, Switzerland: International Organization for Standardization. Google Scholar
- ISO. 2012. ISO 11268–2. Soil quality - Effects of pollutants on earthworms - Part 2. Determination of effects on reproduction of Eisenia fetida/Eisenia andrei. Geneva, Switzerland: International Organization for Standardization. Google Scholar
- J. Jaenicke 1982. “Eisenia foetida” is two biological species. Megadrilogica 4: 6–8. Google Scholar
- S.W. James , D. Porco , T. Decaëns , B. Richard , R. Rougerie & C. Erséus 2010. DNA barcoding reveals cryptic diversity in Lumbricus terrestris L.. 1758 Clitellata): resurrection of L. herculeus (Savigny, 1826). PLoS ONE 5: e15629 [1–8], (doi:10.1371/journal.pone.0015629) Google Scholar
- M. Kimura 1980. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution 16: 111–120. Google Scholar
- R.A. King , A.L. Tibble & W.O.C. Symondson 2008. Opening a can of worms: unprecedented sympatric cryptic diversity within British lumbricid earthworms. Molecular Ecology 17: 4684–4698. Google Scholar
- J. Klarica , A. Kloss-Brandstätter , M. Traugott & A. Juen 2012. Comparing four mitochondrial genes in earthworms - implications for identification, phylogenetics, and discovery of cryptic species. Soil Biology and Biochemistry 45: 23–30. Google Scholar
- K.E. Knott & J. Haimi 2010. High mitochondrial DNA sequence diversity in the parthenogenetic earthworm Dendrobaena octaedra. Heredity 105: 341–347. Google Scholar
- M. Krauss , W. Wilcke & W. Zech 2000. Availability of polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) to earthworms in urban soils. Environmental Science and Technology 34: 4335–4340. Google Scholar
- D. Lin , Q. Zhou , X. Xie & Y. Liu 2010. Potential biochemical and genetic toxicity of triclosan as an emerging pollutant on earthworms (Eisenia fetida). Chemosphere 81: 1328–1333. Google Scholar
- J.L. Lockwood , P. Cassey & T. Blackburn 2005. The role of propagule pressure in explaining species invasions. Trends in Ecology and Evolution 20: 223–228. Google Scholar
- M.S. Maboeta , S.A. Reinecke & A.J. Reinecke 2004. The relationship between lysosomal biomarker and organismal responses in an acute toxicity test with Eisenia fetida (Oligochaeta) exposed to the fungicide copper oxychloride. Environmental Research 96: 95–101. Google Scholar
- M.S. Maboeta & L. van Rensburg 2003a. Bioconversion of sewage sludge and industrially produced woodchips. Water, Air, and Soil Pollution 150: 219–233. Google Scholar
- M.S. Maboeta & L. van Rensburg 2003b. Vermicomposting of industrially produced woodchips and sewage sludge utilizing Eisenia fetida. Ecotoxicology and Environmental Safety 56: 265–270. Google Scholar
- M.S. Maboeta , L. van Rensburg & P.J.J. van Rensburg 2008. Earthworm (Eisenia fetida) bioassay to assess the possible effects of platinum tailings disposal facilities on the environment along a gradient. Applied Ecology and Environmental Research 6: 13–19. Google Scholar
- R.A. Maleri , A.J. Reinecke & S.A. Reinecke 2007. A comparison of nickel toxicity to pre-exposed earthworms (Eisenia fetida, Oligochaeta) in two different test substrates. Soil Biology and Biochemistry.59: 2849–2853. Google Scholar
- R.A. Maleri , A.J. Reinecke & S.A. Reinecke 2008. Metal uptake of two ecophysiologically different earthworms (Eisenia fetida and Aporrectodea caliginosa) exposed to ultramafic soils. Applied Soil Ecology 38: 42–50. Google Scholar
- E. Mayr 1942. Systematics and the origin of species. New York: Columbia University Press. Google Scholar
- A. Miyazaki , T. Amano , H. Saito & Y. Nakano 2002. Acute toxicity of chlorophenols to earthworms using a simple paper contact method and comparison with toxicities to fresh water organisms. Chemosphere 47: 65–69. Google Scholar
- M. Novo , A. Almodóvar & D.J. Díaz-Cosín 2009. High genetic divergence of hormogastrid earthworms (Annelida, Oligochaeta) in the central Iberian Peninsula: evolutionary and demographic implications. Zoologica scripta 38: 537–552. Google Scholar
- M. Novo , A. Almodóvar , R. Fernández , D. Trigo & D.J. Díaz-Cosín 2010. Cryptic speciation of hormogastrid earthworms revealed by mitochondrial and nuclear data. Molecular Phylogenetics and Evolution 56: 507–512. Google Scholar
- OECD. 1984. Guideline for the testing of chemicals No. 207. Earthworm, acute toxicity tests. Paris, France: Organisation for Economic Co-operation and Development. Google Scholar
- OECD. 2004. Guideline for the testing of chemicals No. 222. Earthworm reproduction test (Eisenia fetida/ Eisenia andrei). Paris, France: Organisation for Economic Co-operation and Development. Google Scholar
- N. Øien & J. Stenersen 1984. Esterases of earthworms—III. Electrophoresis reveals that Eisenia fetida (Savigny) is two species. Comparative Biochemistry and Physiology Part C: Comparative Pharmacology 78: 277–282. Google Scholar
- O.J. Owojori , A.J. Reinecke & A.B. Rozanov 2010. Influence of clay content on bioavailability of copper in the earthworm Eisenia fetida. Ecotoxicology and Environmental Safety 73: 407–414. Google Scholar
- O.J. Owojori , A. J. Reinecke , P. Voua-Otomo & S.A. Reinecke 2009. Comparative study of the effects of salinity on life-cycle parameters of four soil-dwelling species (Folsomia candida, Enchytraeus doerjesi, Eisenia fetida and Aporrectodea caliginosa). Pedobiologia 52: 351–360. Google Scholar
- M. Pérez-Losada , J. Eiroa , S. Mato & J. Domínguez 2005. Phylogenetic species delimitation of the earthworms Eisenia fetida (Savigny, 1826) and Eisenia andrei Bouché, 1972 (Oligochaeta, Lumbricidae) based on mitochondrial and nuclear DNA sequences. Pedobiologia 49: 317–324. Google Scholar
- M. Pérez-Losada , M. Ricoy , J.C. Marshall & J. Domínguez 2009. Phylogenetic assessment of the earthworm Aporrectodea caliginosa species complex (Oligochaeta: Lumbricidae) based on mitochondrial and nuclear DNA sequences. Molecular Phylogenetics and Evolution 52: 293–302. Google Scholar
- R.J. Petit , A. El Mousadik & O. Pons 1998. Identifying populations for conservation on the basis of genetic markers. Conservation Biology 12: 844–855. Google Scholar
- M.W. Prinsloo , S.A. Reinecke , W.J. Przybylowicz , J. Mesjasz-Przybylowicz & A.J. Reinecke 1999. Micro-PIXE studies of Cd distribution in the nephridia of the earthworm Eisenia fetida (Oligochaeta). Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 158: 317–322. Google Scholar
- D.H. Reed & R. Frankham 2003. Correlation between fitness and genetic diversity. Conservation Biology 17: 230–237. Google Scholar
- A.J. Reinecke , M.S. Maboeta , L.A. Vermeulen & S.A. Reinecke 2002. Assessment of lead nitrate and mancozeb toxicity in earthworms using the avoidance response. Bulletin of Environmental Contamination and Toxicology 68: 779–786. Google Scholar
- A.J. Reinecke & S.A. Reinecke 2003. The influence of exposure history to lead on the lysosomal response in Eisenia fetida (Oligochaeta). Ecotoxicology and Environmental Safety 55: 30–37. Google Scholar
- A.J. Reinecke , S.A. Reinecke & M.S. Maboeta 2001. Cocoon production and viability as endpoints in toxicity testing of heavy metals with three earthworm species. Pedobiologia 45: 61–68. Google Scholar
- A.J. Reinecke & S.A. Viljoen 1991. A comparison of the biology of Eisenia fetida and Eisenia andrei (Oligochaeta). Biology and Fertility of Soils 11: 295–300. Google Scholar
- S.A. Reinecke & A.J. Reinecke 1997. The influence of lead and manganese on spermatozoa of Eisenia fetida (Oligochaeta). Soil Biology and Biochemistry 29: 737–742. Google Scholar
- S.A. Reinecke & A.J. Reinecke 2004. The comet assay as biomarker of heavy metal genotoxicity in earthworms. Archives of Environmental Contamination and Toxicology46: 208–215. Google Scholar
- M. Saint-Denis , F. Labrot , J.F. Narbonne & D. Ribera 1998. Glutathione, glutathione-related enzymes, and catalase activities in the earthworm Eisenia fetida andrei. Archives of Environmental Contamination and Toxicology 35: 602–614. Google Scholar
- K. Tamura , D. Peterson , N. Peterson , G. Stecher , M. Nei & S. Kumar 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution 28: 2731–2739. Google Scholar
- A. Velando , J. Domínguez & A. Ferreiro 2006. Inbreeding and outbreeding reduces cocoon production in the earthworm Eisenia andrei. European Journal of Soil Biology 42 (suppl. 1): S354–S357. Google Scholar
- P. Voua Otomo , B. Jansen van Vuuren & S.A. Reinecke 2009. Usefulness of DNA barcoding in ecotoxicological investigations: resolving taxonomic uncertainties using Eisenia Malm 1877 as an example. Bulletin of Environmental Contamination and Toxicology 82: 261–264. Google Scholar
- P. Voua Otomo , O.J. Owojori , S.A. Reinecke , S. Daniels & A.J. Reinecke 2011. Using estimates of metal bioavailability in the soil and genetic variation of allozymes to investigate heavy metal tolerance in the earthworm Eisenia fetida (Oligochaeta). Ecotoxicology and Environmental Safety 74: 2070–2074. Google Scholar
- P. Voua Otomo & S.A. Reinecke 2010. Increased cytotoxic and genotoxic tolerance of Eisenia fetida (Oligochaeta) to cadmium after long-term exposure. Ecotoxicology 19: 362–368. Google Scholar