A Megafauna’s Microfauna: Gastrointestinal Parasites of New Zealand’s Extinct Moa (Aves: Dinornithiformes)

Abstract

We perform the first multidisciplinary study of parasites from an extinct megafaunal clade using coprolites from the New Zealand moa (Aves: Dinornithiformes). Ancient DNA and microscopic analyses of 84 coprolites deposited by four moa species (South Island giant moa, Dinornis robustus; little bush moa, Anomalopteryx didiformis; heavy-footed moa, Pachyornis elephantopus; and upland moa, Megalapteryx didinus) reveal an array of gastrointestinal parasites including coccidians (Cryptosporidium and members of the suborder Eimeriorina), nematodes (Heterakoidea, Trichostrongylidae, Trichinellidae) and a trematode (Echinostomida). Parasite eggs were most prevalent and diverse in coprolites from lowland sites, where multiple sympatric moa species occurred and host density was therefore probably higher. Morphological and phylogenetic evidence supports a possible vicariant Gondwanan origin for some of the moa parasites. The discovery of apparently host-specific parasite taxa suggests paleoparasitological studies of megafauna coprolites may provide useful case-studies of coextinction.

About the Authors:

Jamie R. Wood

* E-mail: [email protected]

Affiliation: Landcare Research, Lincoln, Canterbury, New Zealand

Janet M. Wilmshurst

Affiliation: Landcare Research, Lincoln, Canterbury, New Zealand

Nicolas J. Rawlence

Current address: Department of Earth and Ocean Sciences, University of Waikato, Hamilton, New Zealand

Affiliation: Australian Centre for Ancient DNA, School of Earth and Environmental Sciences, University of Adelaide, Adelaide, South Australia, Australia

Karen I. Bonner

Affiliation: Landcare Research, Lincoln, Canterbury, New Zealand

Trevor H. Worthy

Current address: School of Earth and Environmental Sciences, University of Adelaide, Adelaide, South Australia, Australia

Affiliation: School of Biological, Earth and Environmental Sciences, University of New South Wales, Sydney, Australia

John M. Kinsella

Affiliation: HelmWest Laboratory, Missoula, Montana, United States of America

Alan Cooper

Affiliation: Australian Centre for Ancient DNA, School of Earth and Environmental Sciences, University of Adelaide, Adelaide, South Australia, Australia

Introduction

Coprolites (preserved faecal boli) are rich sources of paleoecological information. Historically, the main focus of coprolite studies has been dietary reconstruction [1]–[5]. However, paleoparasitological analyses have also been reported for a large number of coprolites, both archaeological and paleontological in origin [6], and from a broad range of vertebrate taxa including dinosaurs [7], rodents [8], ground sloth [9], lizards [10], carnivores [11], [12], birds [13] and humans [14]–[16]. Identification of parasites in coprolites can provide information on the ecology [17], population dynamics [18] and diseases of prehistoric animals, and the evolution of host-parasite relationships [19], [20].

During the past fifteen years, advances in ancient DNA (aDNA) techniques have provided increased potential for the application of molecular paleoparasitology to coprolites, yet implementation has so far been limited [20], [21]. The potential benefits of aDNA analysis, when used in conjunction with conventional microscopic techniques, include better taxonomic resolution (particularly if only eggs are present) and detection of very small parasites (e.g. Apicomplexa) or those with fragile thin-walled eggs (e.g. Strongylus) that may not preserve intact in coprolites [6]. Paleoparasitological analysis of extinct animal coprolites can not only inform us about host-parasite relationships, but may also be able to shed new light on biodiversity loss due to the process of coextinction where parasitic and mutualistic taxa disappear along with their host taxon. The relative importance of coextinction in total biodiversity loss is poorly understood due to a lack of empirical data [22], but analysis of coprolites may provide a way to obtain such data.

Recently, an ideal resource on which to perform a broad-scale study of gastrointestinal parasites from an extinct megafauna group has been uncovered in New Zealand. Here, accumulations of Holocene coprolites have been excavated from several cave and rock overhang sites on across the South Island [5], [23]–[24]. Identification of the coprolites using aDNA analysis has revealed they were deposited by moa (Aves: Dinornithiformes), a group of large avian herbivores that formerly occurred throughout New Zealand. Nine species of moa [25], [26], ranging from c. 30 to >200 kg [27], all underwent rapid extermination following initial settlement of New Zealand in the 13th Century AD [28]. Analyses of several of the coprolites has already provided answers to some questions regarding moa biology, including diet, habitat-use and niche partitioning [5], [23]–[24]. Here, we use both microscopic and aDNA analyses of parasites from moa (Aves: Dinornithiformes) coprolites to examine host-parasite relationships, and determine whether paleoparasitological analysis of coprolites may provide a suitable method for detecting coextinction events.

Materials and Methods

Coprolite Samples

We used coprolites for which the depositing moa species had previously been determined by aDNA analysis [5, 24, unpublished data] (Table S1). The coprolites were all intact boli and in an excellent state of preservation (figured in 5, and supplementary material of 24). We examined the presence and abundance of parasites in the coprolites (n = 84) using microscopic techniques, and used molecular analyses on a subset of the total (n = 16). The coprolites represent four moa species and three geographic regions of New Zealand’s South Island (Fig. 1): the relatively high rainfall Dart River Valley (c. 500 m elevation), in the mountains of West Otago (heavy-footed moa, Pachyornis elephantopus, n = 8; South Island giant moa, Dinornis robustus, n = 17; little bush moa, Anomalopteryx didiformis, n = 3; upland moa, Megalapteryx didinus, n = 19); the Kawarau and Roxburgh River gorges (c. 200 m elevation) in the semi-arid region of Central Otago (P. elephantopus, n = 2); and the subalpine Euphrates Cave (c. 1000 m elevation), located on the Garibaldi Ridge, Northwest Nelson (M. didinus, n = 35). The Dart River and Central Otago coprolites are of late Holocene age (<3,000 BP) [5], [29], and the Euphrates Cave coprolites extend from the late to mid Holocene (oldest dated sample 6,368±31 radiocarbon years BP) until the approximate time of moa extinction [24]. Unprocessed remainders from the coprolites have been accessioned into the collections of Canterbury Museum for permanent storage (Table S1).

[Figure omitted. See PDF.]

Figure 1. Location and taxonomic representation of moa coprolite study sites.

(a), Moa coprolite from Dart River Valley (b), Moa coprolite sites on the South Island of New Zealand, showing moa taxa represented in coprolite assemblages.

https://doi.org/10.1371/journal.pone.0057315.g001

Microscopic Examination of Parasites

Subsamples (0.31–0.10 g) from each coprolite were boiled in KOH for ten minutes. A known number of exotic Lycopodium spores [30] (1–2 tablets, batch number 483216, mean of 18583 spores per tablet) were added to each sample. The resulting sediment was pipette mixed, and drops were mounted on microscope slides (2–3 per sample) in glycerol jelly. Slides were systematically scanned at 200×magnification, and helminth eggs and Lycopodium spores were counted, allowing quantification of egg abundance. Because of a change in laboratory operating procedures that occurred during this study, samples that were prepared earlier are quantified according to volume (mL), while those prepared later are quantified according to mass (g). The slides are held by Landcare Research, Lincoln, New Zealand.

aDNA Analysis

DNA extraction, polymerase chain reaction (PCR), cloning, and sequencing were carried out following the methods of Wood et al. [24] at the Australian Centre for Ancient DNA. Ancient DNA (aDNA) extraction and PCR setup was carried out in a geographically and physically isolated dedicated aDNA laboratory located 15 minutes from the University of Adelaide campus, where downstream procedures were performed in a modern DNA laboratory. Protocols to control for contamination [31] were strictly followed, including the use of Shrimp DNase to eliminate potential contaminants in PCR reagents.

We designed two sets of primers: (1) Nem18SF (5′-ATTCCGATAACGARCGAGAC-3′) and Nem18SR (5′-CCGCTKRTCCCTCTAAGAAGT-3′); (2) Nem18SlongF (5′-CAGGGCAAGTCTGGTGCCAGCAGC-3′) and Nem18SlongR (5′-GACTTTCGTTCTTGATTAATGAA-3′). Both sets of primers bind to regions that are conserved across a broad range of invertebrates (including the common parasitic helminth groups of apicomplexans, nematodes, and trematodes), and amplify variable regions of the 18S gene (c. 40–120 bp with Nem18S primers; c. 350–400 bp with Nem18Slong primers) (Figs. S1, S2, S3, S4).

Preliminary identification of clone sequences was performed using BLAST. If sequences were obtained using both primer sets, from the same coprolite, and returned identical taxonomic matches in BLAST, then these sequences were assumed to represent the same taxa and were concatenated for the phylogenetic analysis. Sequences from the coprolites were sorted into apicomplexans, nematodes, and trematodes based on nearest BLAST matches. These were aligned (using MUSCLE in Geneious) with 18S sequences (c.1700–1900 bp) from a range of representative taxa from each of these parasite groups (Table S2). The alignments were imported into BEAUti v.1.6.1 and the resulting xml. file was analysed using BEAST v.1.6.1. Our analyses incorporated a HKY model with estimated base frequencies and an age-independent transitions only sequence error model, a relaxed lognormal clock and a Yule tree prior (lognormal birth rate). Maximum credibility trees were produced from MCMC chain lengths of 25 (Apicomplexa, Trematoda) or 50 (Nematoda) million generations (parameters logged every 1000) and assessed for robustness using Tracer v.1.5. Tree output files were summarized using Tree Annotator (10% burnin). Sequences >50 bp in length were deposited in GenBank (Accession numbers KC405320– KC405484).

Results

Microscopic Examination of Parasites

Nematode eggs were observed on microscope slides, although no larvae were seen. The taxonomic usefulness of many helminth eggs is limited and identification can prove difficult [6], [32]. Here we describe the egg types present in the moa coprolites and suggest their likely taxonomic affinities. Prevalence and abundance of the different nematode egg types are shown in Table 1. Coprolites from low altitude sites (Dart River Valley, Kawarau Gorge and Roxburgh Gorge) had higher parasite egg diversity and prevalence overall (mean egg types per moa species = 2.25; eggs in 51% of coprolites) compared with those from the higher altitude Euphrates Cave (mean egg types per moa species = 1; eggs in 17.1% of coprolites). It was not possible to test whether these patterns were significant overall, due to two methods having been used to quantify egg abundance (eggs mL−1 and eggs g−1). However, for coprolites where egg abundance was calculated as eggs g−1, there was a significantly lower prevalence at Euphrates Cave (1000 m elevation) compared with Dart River (500 m elevation) (p = 0.037, t-test).

[Figure omitted. See PDF.]

Table 1. Prevalence and abundance of nematode egg types in moa coprolites.

https://doi.org/10.1371/journal.pone.0057315.t001

Egg type 1 (Fig. 2a).

Description: Ovoid to slightly barrel-shaped (approximately parallel sides). Egg wall relatively thick (c. 5–6 µm), transparent, equal thickness around entire egg, with smooth outer surface. Inner mass often golden colour and shrunken away from outer wall. Dimensions: 50–67.5×31.5–45 µm (mean 61×39 µm) (n = 16). Likely affinity: Eggs are similar to those of Heterakoidea spp. (Nematoda). Hosts: A. didiformis, D. robustus, P. elephantopus, M. didinus. Present in coprolites from the Dart River Valley and Kawarau Gorge.

[Figure omitted. See PDF.]

Figure 2. Helminth eggs from moa coprolites.

(a), egg type 1, cf. Heterakoidea; (b–c), egg type 2, undetermined Nematoda; (d–e), egg type 3, Trichinellidae cf. Capillaria.

https://doi.org/10.1371/journal.pone.0057315.g002

Egg type 2 (Fig. 2b, c).

Description: Ovoid, light brown. Wall c. 5 µm thick. Outer layer with short, evenly spaced spinules protruding up to 3 µm from the surface. The spinulose layer is often absent around the poles. Dimensions: 70–72.5×40–47.5 µm (n = 4). Likely affinity: undetermined Nematoda. Appears to be similar to egg type 1 but with a spinulose layer on the external surface. Hosts: D. robustus, P. elephantopus, M. didinus. Present in coprolites from the Dart River Valley and Kawarau Gorge.

Egg type 3 (Fig. 2d, e).

Description: Elongate ovoid, orange-brown, surface grades from reticulate to longitudinally striate. Obvious polar pores. Dimensions: 52–60×30–35 µm (n = 2). Likely affinity: Trichinellidae (c.f. Capillaria) (Nematoda), due to the presence of polar pores. Hosts: P. elephantopus and M. didinus. Present in coprolites from the Dart River Valley and Euphrates Cave.

aDNA Analysis

In total, 167/233 (71.7%) of the clone sequences using the Nem18S primers and 42/61 (68.9%) of the clone sequences using the Nem18Slong primers were identified as being from parasites. Six distinct clusters of parasite clone sequences, likely to represent identical or very closely related taxa, were identified from sequence alignments and BLAST matches (Figs. S5, S6). Three of the groups had nearest BLAST matches within Apicomplexa, two within Nematoda and one within Trematoda. Non-target sequences included moa, plant, fungi and soil micro-organisms. Fish sequences obtained in initial PCRs were eliminated by using shrimp DNAase, and therefore likely reflected contaminants in the PCR reagents.

Maximum credibility trees.

Single representative sequences from each of the six groups were used in this analysis, and are listed in Table S3.

Apicomplexa

Sequence 1, obtained from a single Dart River M. didinus coprolite (Table 2), was well-supported (posterior value 1.0) within the clade including Cryptosporidium species, and sister (100% bootstrap support) to C. ‘struthionis’, an undescribed strain sequenced from ostrich (Struthio camelus) (Figs. 3, S7). Sequences 2 (M. didinus from Dart River and Euphrates Cave) and 3 (D. robustus from Dart River) were both placed within the suborder Eimeriorina with posterior values of 1.0 and 0.43 respectively (sequence 2 as sister to Calyptosporidae and 3 as sister to all non-Cryptosporidiidae Eimeriorina) (Fig. S7). It should be noted that several families within Eimeriorina (Elleipsisomatidae, Selenococcidiidae, and Spirocystidae) were not represented in the analysis due to unavailability of 18S sequences on Genbank, and this may have affected the exact placement of clones within this group.

[Figure omitted. See PDF.]

Figure 3. Phylogenetic position of 18S sequence groups obtained from moa coprolites.

(a), apicomplexa; (b), nematoda; (c), trematoda. The complete maximum credibility trees are provided as Figs. S7, S8, S9.

https://doi.org/10.1371/journal.pone.0057315.g003

[Figure omitted. See PDF.]

Table 2. Prevalence of six helminth taxa identified by aDNA analysis of moa coprolites.

https://doi.org/10.1371/journal.pone.0057315.t002

Nematoda

Both Nematoda sequences were placed within well-supported clades (Fig. S8). Sequence 4 is nested within the superfamily Heterakoidea (posterior value 1.0) (Figs. 3, S8), and was the most widespread sequence obtained, occurring in coprolites from all four moa species and both Dart River Valley and Euphrates Cave (Table 2). Sequence 5 is nested within the family Trichostrongylidae (posterior value 0.98) (Figs. 3, S8) and was recorded from just a single M. didinus coprolite from Euphrates Cave (Table 2).

Trematoda

Sequence 6 (M. didinus from Dart River and Euphrates Cave; Table 2) is well-supported as being nested within the Trematode order Echinostomida (posterior value 1.0), and sister to Notocotylidae (Figs. 3, S9).

Discussion

Parasite Assemblage

The parasite groups identified from the moa coprolites have all been recorded previously from New Zealand birds [33] and are typical of parasites recorded from extant ratites [34]–[35]. Although there was some overlap in the parasite taxa identified by both microscopic and DNA analysis (Heterakoidea), each analysis method detected taxa that the other did not, highlighting the usefulness of a multidisciplinary approach in paleoparasitological studies. Whereas microscopic analysis was used on more samples and may have detected some of the less prevalent parasite taxa, DNA analysis detected coccidians, which due to their small size may not have preserved as well as larger, thick-walled eggs.

Within each of the six parasite DNA sequence groups (Figs. S5, S6) there is some genetic variation, which may partly be due to ancient DNA damage [36], but may also represent different parasite haplotypes or closely-related species. An example of the latter is evident in group 6 (Echinostomida), where sequences with an A in position 99 were recovered from both Dart River Valley and Euphrates Cave coprolites, but sequences with a G in that position were present only in coprolites from Euphrates Cave (Fig. S5).

Our results for nematode egg counts (Table 1) indicate that moa at lower altitudes had higher parasite prevalence and diversity. This does not appear to be due to a higher number of moa species represented in the lowland coprolite assemblages. For example, all three nematode egg types were detected in M. didinus coprolites from Dart River Valley (overall prevalence 0.31), whereas just one type was present in the M. didinus coprolites from Euphrates Cave (overall prevalence 0.17). This effect is likely due to host density, with has been shown to positively correlate with parasite abundance [37]. In prehuman New Zealand, moa density was probably relatively high in lowland sites, where multiple sympatric species often coexisted (i.e. remains of four moa species from the Dart River Valley, one from Euphrates Cave). However, the lower parasite prevalence observed at Euphrates Cave may also be partly due to the extended temporal range of the coprolites from the site. Whereas coprolites from the low altitude sites are all late Holocene (a period of relative climatic stability), coprolites from Euphrates Cave extend back into the early-mid Holocene period, when climatic conditions in that region were relatively warm and dry [38]. A larger sample of radiocarbon dated coprolites from Euphrates Cave would be required to test any potential temporal affect on parasite prevalence.

Gondwanan Vicariance

The hypothesis that New Zealand’s indigenous fauna and flora represent recent dispersal events following complete submergence of the current landmass during the Oligocene has been a topic of contentious debate during the last decade [39], [40]. Moa have long been regarded as having vicariant Gondwanan origins [41], although recent molecular studies using mitochondrial DNA have presented a range of widely varying taxonomic relationships and evolutionary histories for palaeognaths, some inferring repeated dispersal events and loss of flight (reviewed by Allentoft and Rawlence [42]). However, a recent re-evaluation of morphological relationships suggests that a vicariant Gondwanan origin for moa cannot be ruled out [43]. This is in line with new fossil evidence providing support for the presence of multiple ancient vicariant lineages in New Zealand, including sphenodontine reptiles [44], acanthisittid wrens [45] and freshwater limpets (Latia) [46]. Our results provide some evidence that moa parasites may also have Gondwanan affinities.

Within the Heterakoidea, Inglis & Harris [47] erected the family Kiwinematidae to include the genera Hatterianema and Kiwinema described from New Zealand’s indigenous tuatara (Sphenodon punctatus) and kiwi (Apteryx sp.) respectively. Subsequently, a newly erected African genus Mammalakis has also been placed within Kiwinematidae [48]. Kiwinematidae are characterized by primitive features from which features in other Heterakoidea families could have hypothetically been derived, and thus may represent a Gondwanan relict [47]. Although we cannot definitely attribute the Heterakoidea from moa coprolites to Kiwinematidae because this family is not represented on Genbank, the DNA sequence is sufficiently divergent from both Ascaridia and Heterakis to suggest that it belongs to neither of these widespread genera (pairwise % identities for Nem18Slong fragment: 97.4% Ascaridia galli : moa taxon; 95.9% Heterakis gallinarum : moa taxon; 98.4% Ascaridia galli : Heterakis gallinarum). Further DNA work on extant Kiwinematidae may help resolve the placement of the moa coprolite taxon within Heterakoidea.

Another potential Gondwanan link may lie in the Cryptosporidium sequenced from a moa coprolite. Of seven Cryptosporidium spp. included in the phylogenetic analysis, the coprolite sequence formed a well-supported ‘ratite’ clade with Cryptosporidium ‘struthionis’, basal to the two Cryptosporidium clades reported by Xiao et al. [49] from mammals, snakes and lizards (Figs. 3, S7). A diverse array of avian Cryptosporidium genotypes have been recognized using the SSU rRNA locus [50], and future phylogenetic analyses of this parasite group may provide an interesting complimentary data set with which to understand evolutionary relationships between bird groups.

Coextinction and Habitat Fragmentation

The process of coextinction (loss of parasitic and mutualistic taxa) makes a significant contribution to biodiversity loss during extinction events [51], and could potentially account for the majority of species losses [52]. Therefore, an understanding of the process could assist with making more accurate estimates of the total numbers of species at risk of extinction [51]–[52]. A lack of empirical data has so far precluded an accurate assessment of the importance of coextinction in overall biodiversity loss [52]. Paleoparasitological analysis of coprolites from extinct animals such as moa may provide important quantitative case studies of coextinction events, at least of gastrointestinal parasites.

The question of whether parasite coextinctions occurred in New Zealand is not new; in 1994, Bush and Kennedy [53] pondered whether “when the moas of New Zealand went extinct over a century ago, did they take with them parasite metapopulations or are those parasites found today in other ratites (e.g. the kiwi, emu, cassowary, rhea and ostrich)?”. Although a complete 18S DNA survey of gastrointestinal parasites in extant New Zealand birds would be required to prove whether the taxa present in moa coprolites are now extinct, some of the taxa were identified only from coprolites of Megalapteryx, providing some evidence for host-specificity and likely coextinction (Fig. 4). These taxa include species in the order Echinostomida and suborder Eimeriorina, which were both present in all Megalapteryx coprolites from Dart River and in coprolites from Euphrates Cave, but were not present in coprolites of the three other moa species (Table 2). The Echinostomida sequence was in a well-supported clade with members of the Notocotylidae (Catatropis and Notocotylus), a family of trematodes that inhabit the digestive tracts (commonly caeca) of mammals and birds [54]. Both Catatropis and Notocotylus have been recorded from wild birds in New Zealand, in particular avian herbivores in close association with water (ducks, geese) and wading birds [33]. Aquatic snails, such as the native Potamopyrgus antipodarum are the intermediate hosts [55]. The presence of such a parasite in Megalapteryx would not be unexpected, as evidence from coprolites show that this moa species occasionally fed around margins of alpine tarns and lakes, and grazed aquatic plants [56]. Eimeriorina sp. 2, Trichostrongylidae and Cryptosporidium were also only identified from Megalapteryx, but from single coprolites, so further samples would need to be analysed to test the host-specificity of these taxa.

[Figure omitted. See PDF.]

Figure 4. Network diagrams of identified moa - parasite interactions.

(a) interactions between moa species and parasite egg types, based on morphological identification of parasite eggs from 84 coprolites; (b) interactions between moa species and DNA sequence groups, based on DNA identifications of parasites from a subset of 16 coprolites. Parasite boxes are scaled relative to overall prevalence in analysed coprolite assemblages and moa bones are proportional to the number of coprolites analysed from each species/locality. Localities (in parentheses) are: O, Central Otago; D, Dart River Valley; E, Euphrates Cave.

https://doi.org/10.1371/journal.pone.0057315.g004

The apparent differences in parasite diversity between moa species (Tables 1, 2) could be due to several factors. The evolutionary histories of hosts and parasites are often closely mirrored [57], and the basal position of Megalapteryx within moa [25] may explain the apparent host specificity of several parasites identified from coprolites of this species. However Dinornis is the second most basal moa genus, yet we found similar parasite diversity to Pachyornis and Anomalopteryx, which represent more recent splits within the moa phylogeny [25]. Another possibility is that the parasite diversity may relate to the species’ ecology. For example, Megalapteryx may have favoured feeding near water sources such as small alpine tarns [56] and therefore been more susceptible to waterborne parasites (e.g., Cryptosporidium) and parasites with aquatic intermediate hosts.

Analysis of parasites in coprolites from other New Zealand extant bird species offers the potential to examine how parasite communities deal with severe habitat fragmentation. For example, the extensive pre-European kakapo (Strigops habroptilus) coprolite record across New Zealand [58], [59] could provide a means to contrast past kakapo parasite diversity with that in modern birds, which now have severely contracted population size (<150 individual birds) and distribution (few offshore islands).

Conclusions

We have confirmed the presence of apicomplexan, nematode and trematode gastrointestinal parasites in the coprolites of New Zealand’s extinct moa. Several of these parasites appear to have been host-specific, and therefore are likely to have become extinct with the moa. This study has shown the potential for multidisciplinary paleoparasitological studies of coprolites to contribute to our understanding of evolutionary histories of both parasites and hosts, and to provide case studies of coextinction. The relatively young ages of moa coprolites, and the recent discovery of many specimens from multiple sites, offer an ideal sample for such a study. Questions relating to parasite-host evolutionary histories, and coextinction, could also be answered by detailed paleoparasitological analyses of many Late Quaternary coprolites known from around the world, including ground sloth [9] and mammoth [60].

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Alignment of Nematoda, Trematoda, and Apicomplexa 18S sequences used for designing the Nem18SF and Nem18SR primers.

https://doi.org/10.1371/journal.pone.0057315.s001

(DOC)

Figure S2.

Alignment of Nematoda 18S sequences used for designing the Nem18SlongF and Nem18SlongR primers.

https://doi.org/10.1371/journal.pone.0057315.s002

(DOC)

FigureS3.

Alignment of Trematoda 18S sequences for designing the Nem18SlongF and Nem18SlongR primers.

https://doi.org/10.1371/journal.pone.0057315.s003

(DOC)

Figure S4.

Alignment of Apicomplexa 18S sequences for designing the Nem18SlongF and Nem18SlongR primers.

https://doi.org/10.1371/journal.pone.0057315.s004

(DOC)

Figure S5.

Alignment of clone sequences obtained from moa coprolites using Nem18SF and Nem18SR primers.

https://doi.org/10.1371/journal.pone.0057315.s005

(DOC)

Figure S6.

Alignment of clone sequences obtained from moa coprolites using Nem18SlongF and Nem18SlongR primers.

https://doi.org/10.1371/journal.pone.0057315.s006

(DOC)

Figure S7.

Maximum-credibility tree for 18S sequences of representative Apicomplexa (from Genbank), and moa coprolite sequences 1–3. The tree is rooted with Gymnodium (Dinoflagellata).

https://doi.org/10.1371/journal.pone.0057315.s007

(JPG)

Figure S8.

Maximum-credibility tree for 18S sequences of representative Nematoda (from Genbank), and moa coprolite sequences 4–5. The tree is rooted with Gordius (Nematomorpha). All sequences are correctly resolved within clades representing 5 major orders of parasitic nematodes, except *Litomosoides (Spiruda) and Rondonia (Rhabditida).

https://doi.org/10.1371/journal.pone.0057315.s008

(JPG)

Figure S9.

Maximum-credibility tree for 18S sequences of representative Trematoda (from Genbank), and moa coprolite sequence 6. The tree is rooted with Notocaryoplana (Turbullaria).

https://doi.org/10.1371/journal.pone.0057315.s009

(JPG)

Table S1.

Moa coprolite specimens that were used in this study. Specimen numbers relate to the Australian Centre for Ancient DNA sample database (A). Where voucher specimens exist, the museum registration numbers are also given (CM, Canterbury Museum, New Zealand; OM, Otago Museum, New Zealand).

https://doi.org/10.1371/journal.pone.0057315.s010

(DOC)

Table S2.

18S sequences from GenBank that were used in the phylogenetic analyses.

https://doi.org/10.1371/journal.pone.0057315.s011

(DOC)

Table S3.

Six distinct parasite sequences obtained from moa coprolites that were used in the phylogenetic analyses. Note that the analyses of sequences 1, 2 and 4 used the concatenated sequences from both primer sets (see methods section for explanation).

https://doi.org/10.1371/journal.pone.0057315.s012

(DOC)

Acknowledgments

We thank J. Austin and G. Rattray for laboratory assistance, S. Richardson, A. Camens and C. Wood for fieldwork assistance, and S. Wagstaff and D. Tompkins for helpful discussion. The New Zealand Department of Conservation permitted the collection of coprolites from Dart River and Euphrates Cave (Global Concession Permit Number: CA-31615-OTH).

Author Contributions

Conceived and designed the experiments: JRW JMW THW AC. Performed the experiments: JRW KIB NJR. Analyzed the data: JRW JMK. Wrote the paper: JRW JMW THW AC.

Citation: Wood JR, Wilmshurst JM, Rawlence NJ, Bonner KI, Worthy TH, Kinsella JM, et al. (2013) A Megafauna’s Microfauna: Gastrointestinal Parasites of New Zealand’s Extinct Moa (Aves: Dinornithiformes). PLoS ONE 8(2): e57315. https://doi.org/10.1371/journal.pone.0057315

References

1. James HF, Burney DA (1997) The diet and ecology of Hawaii’s extinct flightless waterfowl: evidence from coprolites. Biol J Linn Soc 62: 279–297.

2. Poinar HN, Hofreiter M, Spaulding WG, Martin PS, Stankiewicz BA, et al. (1998) Molecular coproscopy: Dung and diet of the extinct ground sloth Nothrotheriops shastensis. Science 218: 402–406.

3. Akeret O, Haas JN, Leuzinger U, Jacomet S (1999) Plant macrofossils and pollen in goat/sheep faeces from the Neolithic lake-shore settlement Arbon Bleiche 3, Switzerland. The Holocene 9: 175–182.

4. Kropf M, Mead JI, Anderson RS (2007) Dung, diet, and the paleoenvironment of the extinct shrub-ox (Euceratherium collinum) on the Colorado Plateau, USA. Quat Res 67: 143–151.

5. Wood JR, Rawlence NJ, Rogers GM, Austin JJ, Worthy TH, et al. (2008) Coprolite deposits reveal the diet and ecology of the extinct New Zealand megaherbivore moa (Aves: Dinornithiformes). Quat Sci Rev 27: 2593–2602.

6. Bouchet F, Guidon N, Dittmar K, Harter S, Ferreira LF, et al. (2003) Parasite remains in archaeological sites. Mem Inst Oswald Cruz 98 (supplement 1)47–52.

7. Poinar G, Boucot AJ (2006) Evidence of intestinal parasites of dinosaurs. Parasitology 133: 245–249.

8. Sardella NH, Fugassa MH (2009) Paleoparasitological analysis of rodent coprolites in Holocenic samples from Patagonia, Argentina. J Parasitol 95: 646–651.

9. Schmidt GD, Duszynski DW (1992) Parasites of the extinct Shasta ground sloth, Nothrotheriops shastensis, in Rampart Cave, Arizona. J Parasitol 78: 811–816.

10. Araujo AJG, Confalonieri UEC, Ferreira LF (1982) Oxyurid (Nematoda) egg from coprolites from Brazil. J Parasitol 68: 511–512.

11. Ferreira LF, Araujo A, Duarte AN (1993) Nematode larvae in fossilized animal coprolites from Lower and Middle Pleistocene sites, Central Italy. J Parasitol 79: 440–442.

12. Fugassa MH, Denegri GM, Sardella NH, Araujo A, Guichon RA, et al. (2006) Paleoparasitological records in a canid coprolite from Patagonia, Argentina. J Parasitol 92: 1110–1113.

13. Fugassa MH, Sardella NH, Denegri GM (2007) Paleoparasitological analysis of a raptor pellet from southern Patagonia. J Parasitol 93: 421–422.

14. Leuzinger ML, Bouchet F (2003) Dioctophymidae eggs in coprolite from Neolithic site of Arbon-Bleiche 3 (Switzerland). J Parasitol 89: 1073–1076.

15. Sianto L, Reinhard JJ, Chame M, Chaves S, Mendonca S, et al. (2005) The finding of Echinostoma (Trematoda: Digenea) and Hookworm eggs in coprolites collected from a Brazilian mummified body dated 600–1,200 years before present. J Parasitol 91: 972–975.

16. Fugassa MH, Beltrame MO, Sardella NH, Civalero MT, Aschero C (2010) Paleoparasitological results from coprolites dated at the Pleistocene-Holocene transition as source of paleoecological evidence in Patagonia. J Arch Sci 37: 880–884.

17. Watve MG, Sukumar R (1995) Parasite abundance and diversity in mammals: Correlates with host ecology. Proc Natl Acad Sci USA 92: 8945–8949.

18. Holmstad PR, Hudson PJ, Skorping A (2005) The influence of a parasite community on the dynamics of a host population: a longitudinal study on willow ptarmigan and their parasites. Oikos 111: 377–391.

19. Araujo A, Jansen AM, Bouchet F, Reinhard K, Ferreira LF (2003) Parasitism, the diversity of life, and paleoparasitology. Mem Inst Oswaldo Cruz 98 (supplement 1)5–11.

20. Iniguez AM, Reinhard KJ, Araujo A, Ferreira LF, Vicente ACP (2003) Enterobius vermicularis: ancient DNA from north and south American human coprolites. Mem Inst Oswaldo Cruz 98 (supplement 1)67–69.

21. Loreille O, Roumat E, Verneau O, Bouchet F, Hänni C (2001) Ancient DNA from Ascaris: extraction amplification and sequences from eggs collected in coprolites. Int J Parasitol 31: 1101–1106.

22. Dunn RR, Harris NC, Colwell RK, Koh LP, Sodhi NS (2009) The sixth mass coextinction: are most endangered species parasites and mutualists? Proc R Soc London Ser B 276: 3037–3045.

23. Wood JR, Wilmshurst JM, Worthy TH, Cooper A (2012) First coprolite evidence for the diet of Anomalopteryx didiformis, an extinct forest ratite from New Zealand. N Z J Ecol 36: 164–170.

24. Wood JR, Wilmshurst JM, Wagstaff SJ, Worthy TH, Rawlence NJ, et al. (2012) High-reolution coproecology: using coprolites to reconstruct the habits and habitats of New Zealand’s extinct upland moa (Megalapteryx didinus). PLOS ONE 7: e40025.

25. Bunce M, Worthy TH, Phillips MJ, Holdaway RN, Willerslev W, et al. (2009) The evolutionary history of the extinct ratite moa and New Zealand Neogene paleogeography. Proc Natl Acad Sci USA 106: 20646–20651.

26. Gill BJ, Bell BD, Chambers GK, Medway DG, Palma RL, et al.. (2010) Checklist of the birds of New Zealand, Norfolk and Macquarie Islands, and the Ross Dependency, Antarctica. Wellington: Ornithological Society of New Zealand and Te Papa Press. 464 p.

27. Dickison MJ (2007) Allometry of giant flightless birds. Unpublished PhD thesis, Duke University, USA. 125p.

28. Wilmshurst JM, Anderson AJ, Higham TFG, Worthy TH (2008) Dating the late prehistoric dispersal of Polynesians to New Zealand using the commensal Pacific rat. Proc Natl Acad Sci USA 105: 7676–7680.

29. Wood JR, Wilmshurst JM, Rawlence NJ (2011) Radiocarbon-dated faunal remains correlate very large rock avalanche deposit with prehistoric Alpine fault rupture. N Z J Geol Geophys 54: 431–434.

30. Moore PD, Collinson M, Webb JA (1991) Pollen Analysis. Oxford: Blackwell Scientific. 216 p.

31. Cooper A, Poinar HN (2000) Ancient DNA: Do it right or not at all. Science 289: 1139.

32. Reinhard KJ (1998) Parasitology. In: Cockburn A, Cockburn E, Reyman TA, editors. Mummies, disease and ancient cultures. Cambridge: Cambridge University Press. 377–380.

33. McKenna P (2010) An updated checklist of helminth and protozoan parasites of birds of New Zealand. WebmedCentral Parasitology 1, WMC00705.

34. Craig TM, Diamond PL (1996) Parasites of ratites. In: Tully TN, Shane SM, editors. Ratite management, medicine and surgery. Malabar, Flordia: Krieger Publishing Company. 115–126.

35. Ponce Gordo F, Herrera S, Castro AT, Garcia Duran B, Martinez Diaz RA (2002) Parasites from farmed ostriches (Struthio camelus) and rheas (Rhea americana) in Europe. Vet Parasitol 107: 137–160.

36. Mitchell D, Willerslev E, Hansen A (2005) Damage and repair of ancient DNA. Mutat Res 571: 265–276.

37. Arneberg P, Skorping A, Grenfell B, Read AF (1998) Host densities as determinants of abundance in parasite communities. Proc R Soc London Ser B 265: 1283–1289.

38. Williams PW, King DNT, Zhao JX, Collerson KD (2005) Late Pleistocene to Holocene composite speleothem 18O and 13C chronologies from South Island, New Zealand – did a global Younger Dryas really exist? Earth Planet Sci. Lett. 230: 301–317.

39. Gibbs G (2006) Ghosts of Gondwana: The history of life in New Zealand. Nelson: Potton. 232p.

40. Waters JM, Craw D (2006) Goodbye Gondwana? New Zealand biogeography, geology, and the problem of cicrularity. Syst Biol 55: 351–356.

41. Fleming CA (1975) The geological history of New Zealand and its biota. In: Kuschel G, editor. Biogeography and ecology in New Zealand. The Hague: Junk. 1–86.

42. Allentoft ME, Rawlence NJ (2011) Moa’s Ark or volant ghosts of Gondwana? Insights from nineteen years of ancient DNA research on the extinct moa (Aves: Dinornithiformes) of New Zealand. Ann Anat 194: 36–51.

43. Worthy TH, Scofield RP (2012) 21st Century advances in knowledge of the biology of moa (Aves: Dinornithiformes), diagnoses revised, and a new morphological analysis. N Z J Zool 39: 87–153.

44. Jones ME, Tennyson AJD, Worthy JP, Evans SE, Worthy TH (2009) A sphenodontine (Rhynchocephalia) from the Miocene of New Zealand and palaeobiogeography of the tuatara (Sphenodon). Proc R Soc London Ser B 276: 1385–1390.

45. Worthy TH, Hand SJ, Nguyen J, Tennyson AJD, Worthy J, et al. (2010) Biogeographical and phylogenetic implications of an Early Miocene wren (Aves: Passeriformes: Acanthisittidae) from New Zealand. J Vert Pal 30: 479–498.

46. Marshall BA (2011) A new species of Latia Gray, 1850 (Gastropoda: Pulmonata: Hygrophila: Chilinoidea: Latiidae) from Miocene Palaeo-lake Manuherikia, southern New Zealand, and biogeographic implications. Mollusc Res 31: 47–52.

47. Inglis WG, Harris EA (1990) Kiwinematidae n. fam. (Nematoda) for Kiwinema n. g. and Hatterianema Chabaud & Dollfus, 1966: heterakoids of native New Zealand vertebrates. Syst Parasitol 15: 75–79.

48. Inglis WG (1991) Mammalakis n. g. and Mammalakinae n. subfam. (Nematoda: Heterakoidea: Kiwinematidae): parasites of mole rat (Rodentia: Bathyergidae and Splacidae). Syst Parasitol 20: 89–95.

49. Xiao L, Escalante L, Yang C, Sulaiman I, Escalante AE, et al. (1999) Phylogenetic analysis of Cryptosporidium parasites based on the small-subunit rRNA gene locus. Appl Environ Microbiol 65: 1578–1583.

50. Ryan U (2010) Cryptosporidium in birds, fish and amphibians. Exp Parasitol 124: 113–120.

51. Koh LP, Dunn RR, Sodhi NS, Colwell RK, Proctor HC, et al. (2004) Species coextinctions and the biodiversity crisis. Science 305: 1632–1634.

52. Dunn RR (2009) Coextinction: anecdotes, models, and speculation. In: Turvey S, editor. Holocene extinctions. Oxford: Oxford University Press. 167–180.

53. Bush AO, Kennedy CR (1994) Host fragmentation and helminth parasites: Hedging your bets against extinction. Int J Parasitol 24: 1333–1343.

54. Barton DP, Blair D (2005) Family Notocotylidae Lühe, 1909. In: Jones A, Bray RA, Gibson DI, editors. Keys to the Trematoda, Vol 2. London: CAB International and The Natural History Museum.

55. Morley NJ (2008) The role of the invasive snail Potamopyrgus antipodarum in the transmission of trematode parasites in Europe and its implications for ecotoxological studies. Aquat Sci 70: 107–114.

56. Horrocks M, D’Costa D, Wallace R, Gardner R, Kondo R (2004) Plant remains in coprolites: diet of a sub-alpine moa (Dinornithiformes) from southern New Zealand. Emu 104: 149–156.

57. Clayton DH, Al-Tamimi S, Johnson KP (2003). The ecological basis of coevolutionary history. In: Page RDM, editor. Tangled trees. Phylogeny, cospeciation, and coevolution. Chicago: University of Chicago Press. 310–341.

58. Horrocks M, Salter J, Braggins J, Nichol S, Moorhouse R, et al. (2008) Plant microfossil analysis of coprolites of the critically endangered kakapo (Strigops habroptilus) parrot from New Zealand. Rev Palaeobot Palyno 149: 229–245.

59. Wood JR, Wilmshurst JM, Worthy TH, Holzapfel AS, Cooper A (2012) A lost link between a flightless parrot and a parasitic plant and the potential role of coprolites in conservation paleobiology. Cons Biol 26: 1091–1099.

60. Davis OK, Agenbroad L, Martin PS, Mead JI (1984) The Pleistocene dung blanket of Bechan Cave, Utah. Spec Publ Carnegie Mus Nat Hist 8: 267–282.