Stephane Boyer

Spider: an R package for the analysis of species identity and evolution, with particular reference to DNA barcoding.

Spider: SPecies IDentity and Evolution in R is a new R package implementing a number of useful analyses for DNA barcoding studies and associated research into species delimitation and speciation. Included are functions essential for generating important summary statistics from DNA barcode data, assessing specimen identification efficacy, and for testing and optimizing divergence threshold limits. In terms of investigating evolutionary and taxonomic questions, techniques for assessing diagnostic nucleotides and probability of reciprocal monophyly are also provided. Additionally, a sliding window function offers opportunities to analyse information across a gene, essential for marker design in degraded DNA studies. Spider capitalizes on R’s extensible ethos and offers an integrated platform ideal for the analysis of both nucleotide and morphological data. The program can be obtained from the comprehensive R archive network (CRAN, http://cran.r‐project.org) and from the R‐Forge package development site (http://spider.r‐forge.r‐project.org/).

Sliding window analyses for optimal selection of mini-barcodes, and application to 454-pyrosequencing for specimen identification from degraded DNA.

DNA barcoding remains a challenge when applied to diet analyses, ancient DNA studies, environmental DNA samples and, more generally, in any cases where DNA samples have not been adequately preserved. Because the size of the commonly used barcoding marker (COI) is over 600 base pairs (bp), amplification fails when the DNA molecule is degraded into smaller fragments. However, relevant information for specimen identification may not be evenly distributed along the barcoding region, and a shorter target can be sufficient for identification purposes. This study proposes a new, widely applicable, method to compare the performance of all potential ‘mini-barcodes’ for a given molecular marker and to objectively select the shortest and most informative one. Our method is based on a sliding window analysis implemented in the new R package SPIDER (Species IDentity and Evolution in R). This method is applicable to any taxon and any molecular marker. Here, it was tested on earthworm DNA that had been degraded through digestion by carnivorous landsnails. A 100 bp region of 16 S rDNA was selected as the shortest informative fragment (mini-barcode) required for accurate specimen identification. Corresponding primers were designed and used to amplify degraded earthworm (prey) DNA from 46 landsnail (predator) faeces using 454-pyrosequencing. This led to the detection of 18 earthworm species in the diet of the snail. We encourage molecular ecologists to use this method to objectively select the most informative region of the gene they aim to amplify from degraded DNA. The method and tools provided here, can be particularly useful (1) when dealing with degraded DNA for which only small fragments can be amplified, (2) for cases where no consensus has yet been reached on the appropriate barcode gene, or (3) to allow direct analysis of short reads derived from massively parallel sequencing without the need for bioinformatic consolidation.

Using Next-Generation Sequencing to Analyse the Diet of a Highly Endangered Land Snail (Powelliphanta augusta) Feeding on Endemic Earthworms

Predation is often difficult to observe or quantify for species that are rare, very small, aquatic or nocturnal. The assessment of such species’ diet can be conducted using molecular methods that target prey DNA remaining in predators’ guts and faeces. These techniques do not require high taxonomic expertise, are applicable to soft-bodied prey and allow for identification at the species level. However, for generalist predators, the presence of mixed prey DNA in guts and faeces can be a major impediment as it requires development of specific primers for each potential prey species for standard (Sanger) sequencing. Therefore, next generation sequencing methods have recently been applied to such situations. In this study, we used 454-pyrosequencing to analyse the diet of Powelliphantaaugusta , a carnivorous landsnail endemic to New Zealand and critically endangered after most of its natural habitat has been lost to opencast mining. This species was suspected to feed mainly on earthworms. Although earthworm tissue was not detectable in snail faeces, earthworm DNA was still present in sufficient quantity to conduct molecular analyses. Based on faecal samples collected from 46 landsnails, our analysis provided a complete map of the earthworm-based diet of P . augusta . Predated species appear to be earthworms that live in the leaf litter or earthworms that come to the soil surface at night to feed on the leaf litter. This indicates that P . augusta may not be selective and probably predates any earthworm encountered in the leaf litter. These findings are crucial for selecting future translocation areas for this highly endangered species. The molecular diet analysis protocol used here is particularly appropriate to study the diet of generalist predators that feed on liquid or soft-bodied prey. Because it is non-harmful and non-disturbing for the studied animals, it is also applicable to any species of conservation interest.

Nationwide inventory of mosquito biodiversity (Diptera: Culicidae) in Belgium, Europe

To advance our restricted knowledge on mosquito biodiversity and distribution in Belgium, a national inventory started in 2007 (MODIRISK) based on a random selection of 936 collection points in three main environmental types: urban, rural and natural areas. Additionally, 64 sites were selected because of the risk of importing a vector or pathogen in these sites. Each site was sampled once between May and October 2007 and once in 2008 using Mosquito Magnet Liberty Plus traps. Diversity in pre-defined habitat types was calculated using three indices. The association between species and environmental types was assessed using a correspondence analysis. Twenty-three mosquito species belonging to traditionally recognized genera were found, including 21 indigenous and two exotic species. Highest species diversity (Simpson 0.765) and species richness (20 species) was observed in natural areas, although urban sites scored also well (Simpson 0.476, 16 species). Four clusters could be distinguished based on the correspondence analysis. The first one is related to human modified landscapes (such as urban, rural and industrial sites). A second is composed of species not associated with a specific habitat type, including the now widely distributed Anopheles plumbeus. A third group includes species commonly found in restored natural or bird migration areas, and a fourth cluster is composed of forest species. Outcomes of this study demonstrate the effectiveness of the designed sampling scheme and support the choice of the trap type. Obtained results of this first country-wide inventory of the Culicidae in Belgium may serve as a basis for risk assessment of emerging mosquito-borne diseases.

Noninvasive molecular methods to identify live scarab larvae: an example of sympatric pest and nonpest species in New Zealand

Despite the negative impact that many scarab larvae have on agro‐ecosystems, very little attention has been paid to their taxonomy. Their often extremely similar morphological characteristics have probably contributed to this impediment, which has also meant that they are very difficult to identify in the field. Molecular methods can overcome this challenge and are particularly useful for the identification of larvae to enable management of pest species occurring sympatrically with nonpest species. However, the invasive collection of DNA samples for such molecular methods is not compatible with subsequent behavioural, developmental or fitness studies. Two noninvasive DNA sampling and DNA analysis methods suitable for the identification of larvae from closely related scarab species were developed here. Using the frass and larval exuviae as sources of DNA, field‐collected larvae of Costelytra zealandica (White) and Costelytra brunneum (Broun) (Scarabaeidae: Melolonthinae) were identified by multiplex PCR based on the difference in size of the resulting PCR products. This study also showed that small quantities of frass can be used reliably even 7 days after excretion. This stability of the DNA is of major importance in ecological studies where timeframes rarely allow daily monitoring. The approach developed here is readily transferable to the study of any holometabolous insect species for which morphological identification of larval stages is difficult.

Un‐nesting DNA Russian dolls – the potential for constructing food webs using residual DNA in empty aphid mummies

Constructing food‐web assemblages comprising parasitoid wasps involves large field collections of hosts followed by labour‐intensive rearing of the insects to evaluate the rates of parasitism along with morphological or molecular identification of the parasitoid species. This article presents research towards a new molecular method for the practical and accurate construction of aphid‐based food webs. We hypothesize that parasitoid and hyperparasitoid DNA left inside aphid mummies after emergence of these third and fourth trophic‐level guilds can be simultaneously detected using universal polymerase chain reaction (PCR) primers for nonspecific DNA amplification in combination with single‐stranded conformation polymorphism (SSCP) analysis. Such a protocol theoretically allows food‐web construction to be performed with no a priori knowledge of the species present. Moreover, the use of empty mummies circumvents rearing and minimizes labour and time in the field and laboratory. To test our hypothesis, we conducted DNA analyses on laboratory‐produced parasitized aphids (mummies) from Myzus persicae and Brevicoryne brassicae (two important aphid pest species) after exposure to the parasitoid Diaeretiella rapae and the hyperparasitoid Asaphes vulgaris. DNA is amplifiable in empty aphid mummies for as long as 3 weeks after parasitoid emergence. However, the simultaneous identification of several species in a single mummy sample was rare, which hinders the accurate inference of trophic links. DNA quality and relative quantity, together with preferential amplification, are potential explanations of current results. Technical refinements are needed to ensure full reliability and detection of complex trophic links. The use of PCR‐SSCP for food‐web construction is novel, and its potential to include an important number of different species is yet to be fully explored.

Diet Composition and Variability of Wild Octopus vulgaris and Alloteuthis media (Cephalopoda) Paralarvae: A Metagenomic Approach

The high mortality of cephalopod early stages is the main bottleneck to grow them from paralarvae to adults in culture conditions, probably because the inadequacy of the diet that results in malnutrition. Since visual analysis of digestive tract contents of paralarvae provides little evidence of diet composition, the use of molecular tools, particularly next generation sequencing (NGS) platforms, offers an alternative to understand prey preferences and nutrient requirements of wild paralarvae. In this work, we aimed to determine the diet of paralarvae of the loliginid squid Alloteuthis media and to enhance the knowledge of the diet of recently hatched Octopus vulgaris paralarvae collected in different areas and seasons in an upwelling area (NW Spain). DNA from the dissected digestive glands of 32 A. media and 64 O. vulgaris paralarvae was amplified with universal primers for the mitochondrial gene COI, and specific primers targeting the mitochondrial gene 16S gene of arthropods and the mitochondrial gene 16S of Chordata. Following high-throughput DNA sequencing with the MiSeq run (Illumina), up to 4,124,464 reads were obtained and 234,090 reads of prey were successfully identified in 96.87 and 81.25% of octopus and squid paralarvae, respectively. Overall, we identified 122 Molecular Taxonomic Units (MOTUs) belonging to several taxa of decapods, copepods, euphausiids, amphipods, echinoderms, molluscs, and hydroids. Redundancy analysis (RDA) showed seasonal and spatial variability in the diet of O. vulgaris and spatial variability in A. media diet. General Additive Models (GAM) of the most frequently detected prey families of O. vulgaris revealed seasonal variability of the presence of copepods (family Paracalanidae) and ophiuroids (family Euryalidae), spatial variability in presence of crabs (family Pilumnidae) and preference in small individual octopus paralarvae for cladocerans (family Sididae) and ophiuroids. No statistically significant variation in the occurrences of the most frequently identified families was revealed in A. media. Overall, these results provide new clues about dietary preferences of wild cephalopod paralarvae, thus opening up new scenarios for research on trophic ecology and digestive physiology under controlled conditions.

Using molecular tools to identify New Zealand endemic earthworms in a mine restoration project

Abstract A restoration ecology project was commenced on the West Coast region of New Zealand to re-establish the local fauna of endemic Powelliphanta spp. carnivorous landsnails at an opencast coal mine site after mining activities. The aim of the current research is to provide recommendations for the use of earthworms to improve the restoration of ecological communities, especially the landsnails. To provide such recommendations, different aspects of the ecology and bio-systematics of the New Zealand endemic earthworm fauna have been studied using molecular techniques. About 1,500 earthworm individuals have been collected across 17 sampling sites in the Stockton mine area. In New Zealand, 173 endemic earthworm species are known. Only minor revisions to the earthworms’ taxonomic status have been made since 1959. Species identification was performed by morphological analysis (following Lee’s taxonomic key) and molecular analysis (using the mitochondrial 16S gene). The latter analyses conducted on a selection of 83 individuals revealed the existence of at least 17 different taxa, most of which are probably undescribed species. Some of these earthworm species are predated by an endangered carnivorous landsnail, Powelliphanta augusta Walker, Trewick & Barker. Because the conservation of P. augusta may rely greatly on the understanding of their diet, earthworm DNA was sought after in the snails’ feces, using molecular analyses. Molecular analyses have been helpful in establishing an inventory of the species present in the study site, facilitating new species taxonomic descriptions and elucidating the predatorprey relationship.

New records of springtails in New Zealand pasture: how well are our pastoral invertebrates known?

Abstract Twenty-four collembolan species are recorded from improved pastures and clovers in New Zealand, of which 17 can be named to species or probable species, the others only to genus. Of the 17 named species, nine have been recorded before from New Zealand but the other eight are new records for the country. All named species are considered as introduced to New Zealand, probably originally from Europe and are unlikely to colonise native habitats. As all named species reported as new records can be abundant at times, this indicates poor knowledge of a major part of New Zealands agricultural fauna. Collembola are a group of important microarthropod detritivores that make a significant contribution to ecosystem services. The absence so far of quantification of the contribution this and other soil groups make to ecological resilience and function is a serious problem.

Invasion success of a scarab beetle within its native range: host range expansion versus host-shift

Only recently has it been formally acknowledged that native species can occasionally reach the status of ‘pest’ or ‘invasive species’ within their own native range. The study of such species has potential to help unravel fundamental aspects of biological invasions. A good model for such a study is the New Zealand native scarab beetle, Costelytra zealandica (White), which even in the presence of its natural enemies has become invasive in exotic pastures throughout the country. Because C. zealandica still occurs widely within its native habitat, we hypothesised that this species has only undergone a host range expansion (ability to use equally both an ancestral and new host) onto exotic hosts rather than a host shift (loss of fitness on the ancestral host in comparison to the new host). Moreover, this host range expansion could be one of the main drivers of its invasion success. In this study, we investigated the fitness response of populations of C. zealandica from native and exotic flora, to several feeding treatments comprising its main exotic host plant as well as one of its ancestral hosts. Our results suggest that our initial hypothesis was incorrect and that C. zealandica populations occurring in exotic pastures have experienced a host-shift rather than simply a host-range expansion. This finding suggests that an exotic plant introduction can facilitate the evolution of a distinct native host-race, a phenomenon often used as evidence for speciation in phytophagous insects and which may have been instrumental to the invasion success of C. zealandica.