Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed-invected.

Polymorphic Batesian mimics exhibit multiple protective morphs that each mimic a different noxious model. Here, we study the genomic transitions leading to the evolution of different mimetic wing patterns in the polymorphic Mocker Swallowtail Papilio dardanus. We generated a draft genome (231 Mb over 30 chromosomes) and re-sequenced individuals of three morphs. Genome-wide single nucleotide polymorphism (SNP) analysis revealed elevated linkage disequilibrium and divergence between morphs in the regulatory region of engrailed, a developmental gene previously implicated in the mimicry switch. The diverged region exhibits a discrete chromosomal inversion (of 40 kb) relative to the ancestral orientation that is associated with the cenea morph, but not with the bottom-recessive hippocoonides morph or with non-mimetic allopatric populations. The functional role of this inversion in the expression of the novel phenotype is currently unknown, but by preventing recombination, it allows the stable inheritance of divergent alleles enabling geographic spread and local coexistence of multiple adaptive morphs.

The contribution of mitochondrial metagenomics to large-scale data mining and phylogenetic analysis of Coleoptera.

A phylogenetic tree at the species level is still far off for highly diverse insect orders, including the Coleoptera, but the taxonomic breadth of public sequence databases is growing. In addition, new types of data may contribute to increasing taxon coverage, such as metagenomic shotgun sequencing for assembly of mitogenomes from bulk specimen samples. The current study explores the application of these techniques for large-scale efforts to build the tree of Coleoptera. We used shotgun data from 17 different ecological and taxonomic datasets (5 unpublished) to assemble a total of 1942 mitogenome contigs of >3000 bp. These sequences were combined into a single dataset together with all mitochondrial data available at GenBank, in addition to nuclear markers widely used in molecular phylogenetics. The resulting matrix of nearly 16,000 species with two or more loci produced trees (RAxML) showing overall congruence with the Linnaean taxonomy at hierarchical levels from suborders to genera. We tested the role of full-length mitogenomes in stabilizing the tree from GenBank data, as mitogenomes might link terminals with non-overlapping gene representation. However, the mitogenome data were only partly useful in this respect, presumably because of the purely automated approach to assembly and gene delimitation, but improvements in future may be possible by using multiple assemblers and manual curation. In conclusion, the combination of data mining and metagenomic sequencing of bulk samples provided the largest phylogenetic tree of Coleoptera to date, which represents a summary of existing phylogenetic knowledge and a defensible tree of great utility, in particular for studies at the intra-familial level, despite some shortcomings for resolving basal nodes.

The phylogeny of Galerucinae (Coleoptera: Chrysomelidae) and the performance of mitochondrial genomes in phylogenetic inference compared to nuclear rRNA genes.

With efficient sequencing techniques, full mitochondrial genomes are rapidly replacing other widely used markers, such as the nuclear rRNA genes, for phylogenetic analysis but their power to resolve deep levels of the tree remains controversial. We studied phylogenetic relationships of leaf beetles (Chrysomelidae) in the tribes Galerucini and Alticini (root worms and flea beetles) based on full mitochondrial genomes (103 newly sequenced), and compared their performance to the widely sequenced nuclear rRNA genes (full 18S, partial 28S). Our results show that: (i) the mitogenome is phylogenetically informative from subtribe to family level, and the per-nucleotide contribution to nodal support is higher than that of rRNA genes, (ii) the Galerucini and Alticini are reciprocally monophyletic sister groups, if the classification is adjusted to accommodate several 'problematic genera' that do not fit the dichotomy of lineages based on the presence (Alticini) or absence (Galerucini) of the jumping apparatus, and (iii) the phylogenetic results suggest a new classification system of Galerucini with eight subtribes: Oidina, Galerucina, Hylaspina, Metacyclina, Luperina, Aulacophorina, Diabroticina and Monoleptina.

Soup to Tree: The Phylogeny of Beetles Inferred by Mitochondrial Metagenomics of a Bornean Rainforest Sample.

In spite of the growth of molecular ecology, systematics and next-generation sequencing, the discovery and analysis of diversity is not currently integrated with building the tree-of-life. Tropical arthropod ecologists are well placed to accelerate this process if all specimens obtained through mass-trapping, many of which will be new species, could be incorporated routinely into phylogeny reconstruction. Here we test a shotgun sequencing approach, whereby mitochondrial genomes are assembled from complex ecological mixtures through mitochondrial metagenomics, and demonstrate how the approach overcomes many of the taxonomic impediments to the study of biodiversity. DNA from approximately 500 beetle specimens, originating from a single rainforest canopy fogging sample from Borneo, was pooled and shotgun sequenced, followed by de novo assembly of complete and partial mitogenomes for 175 species. The phylogenetic tree obtained from this local sample was highly similar to that from existing mitogenomes selected for global coverage of major lineages of Coleoptera. When all sequences were combined only minor topological changes were induced against this reference set, indicating an increasingly stable estimate of coleopteran phylogeny, while the ecological sample expanded the tip-level representation of several lineages. Robust trees generated from ecological samples now enable an evolutionary framework for ecology. Meanwhile, the inclusion of uncharacterized samples in the tree-of-life rapidly expands taxon and biogeographic representation of lineages without morphological identification. Mitogenomes from shotgun sequencing of unsorted environmental samples and their associated metadata, placed robustly into the phylogenetic tree, constitute novel DNA "superbarcodes" for testing hypotheses regarding global patterns of diversity.

Bulk de novo mitogenome assembly from pooled total DNA elucidates the phylogeny of weevils (Coleoptera: Curculionoidea).

Complete mitochondrial genomes have been shown to be reliable markers for phylogeny reconstruction among diverse animal groups. However, the relative difficulty and high cost associated with obtaining de novo full mitogenomes have frequently led to conspicuously low taxon sampling in ensuing studies. Here, we report the successful use of an economical and accessible method for assembling complete or near-complete mitogenomes through shot-gun next-generation sequencing of a single library made from pooled total DNA extracts of numerous target species. To avoid the use of separate indexed libraries for each specimen, and an associated increase in cost, we incorporate standard polymerase chain reaction-based "bait" sequences to identify the assembled mitogenomes. The method was applied to study the higher level phylogenetic relationships in the weevils (Coleoptera: Curculionoidea), producing 92 newly assembled mitogenomes obtained in a single Illumina MiSeq run. The analysis supported a separate origin of wood-boring behavior by the subfamilies Scolytinae, Platypodinae, and Cossoninae. This finding contradicts morphological hypotheses proposing a close relationship between the first two of these but is congruent with previous molecular studies, reinforcing the utility of mitogenomes in phylogeny reconstruction. Our methodology provides a technically simple procedure for generating densely sampled trees from whole mitogenomes and is widely applicable to groups of animals for which bait sequences are the only required prior genome knowledge.

Phylogenetic community ecology of soil biodiversity using mitochondrial metagenomics.

High-throughput DNA methods hold great promise for the study of taxonomically intractable mesofauna of the soil. Here, we assess species diversity and community structure in a phylogenetic framework, by sequencing total DNA from bulk specimen samples and assembly of mitochondrial genomes. The combination of mitochondrial metagenomics and DNA barcode sequencing of 1494 specimens in 69 soil samples from three geographic regions in southern Iberia revealed >300 species of soil Coleoptera (beetles) from a broad spectrum of phylogenetic lineages. A set of 214 mitochondrial sequences longer than 3000 bp was generated and used to estimate a well-supported phylogenetic tree of the order Coleoptera. Shorter sequences, including cox1 barcodes, were placed on this mitogenomic tree. Raw Illumina reads were mapped against all available sequences to test for species present in local samples. This approach simultaneously established the species richness, phylogenetic composition and community turnover at species and phylogenetic levels. We find a strong signature of vertical structuring in soil fauna that shows high local community differentiation between deep soil and superficial horizons at phylogenetic levels. Within the two vertical layers, turnover among regions was primarily at the tip (species) level and was stronger in the deep soil than leaf litter communities, pointing to layer-mediated drivers determining species diversification, spatial structure and evolutionary assembly of soil communities. This integrated phylogenetic framework opens the application of phylogenetic community ecology to the mesofauna of the soil, among the most diverse and least well-understood ecosystems, and will propel both theoretical and applied soil science.

The evolutionary genetics of highly divergent alleles of the mimicry locus in Papilio dardanus.

BACKGROUND: The phylogenetic history of genes underlying phenotypic diversity can offer insight into the evolutionary origin of adaptive traits. This is especially true where single genes have large phenotypic effects, for example in determining polymorphic mimicry in butterflies. Here, we characterise the evolutionary history of two candidate genes for the mimicry switch in the polymorphic Batesian mimic Papilio dardanus coding for the transcription factors engrailed and invected. RESULTS: We show that phased haplotypes associated with the dominant morphs f. poultoni and f. planemoides are phylogenetically highly divergent, in particular at non-synonymous sites. Some non-synonymous changes are shared between the divergent alleles suggesting either convergence or a shared ancestry. Gene trees for invected do not show this pattern. Despite their great divergence, all engrailed alleles of P. dardanus were monophyletic with respect to alleles of closely related species. Phylogenetic analyses therefore reveal no evidence for introgression from other species. A McDonald-Kreitman test conducted on a population sample from South Africa confirms a significant excess of intraspecific non-synonymous diversity in P. dardanus engrailed, suggesting long-term balanced polymorphism at this locus. CONCLUSIONS: The divergence between engrailed haplotypes suggests an evolutionary history distorted by selection with multiple changes reflecting recurrent selective sweeps. The high level of intraspecific polymorphism observed is characteristic of balancing selection on this locus, as expected if the gene engrailed is under phenotypic selection for the maintenance of multiple mimetic morphs. Non-synonymous changes in key functional portions of a major transcription factor are likely to be deleterious but if maintained in a dominant allele at low frequency, heterozygosity would reduce the associated genetic load.

A functional isopenicillin N synthase in an animal genome.

Horizontal transfer of genes is widespread among prokaryotes, but is less common between microorganisms and animals. Here, we present evidence for the presence of a gene encoding functional isopenicillin N synthase, an enzyme in the ß-lactam antibiotics biosynthesis pathway, in the genome of the soil-living collembolan species, Folsomia candida (FcIPNS). At present, this gene is only known from bacteria and fungi, as is the capacity to produce ß-lactam antibiotics. The FcIPNS gene was located on two genomic contigs, was physically linked to a predicted insect ATP-binding cassette transporter gene, and contained three introns each flanked by eukaryotic splicing recognition sites (GT/AG). Homology searches revealed no similarity between these introns and the FcIPNS regions of bacteria or fungi. All amino acids conserved across bacteria and fungi were also conserved in F. candida. Recombinant FcIPNS was able to convert its substrate amino δ-(l-α-aminoadipyl)-l-cysteinyl-d-valine into isopenicillin N, providing strong evidence that FcIPNS is functional. Phylogenetic analysis clustered FcIPNS outside the bacterial IPNS clade, and also outside the fungal IPNS clade, suggesting an ancient gene transfer followed by divergence in the F. candida genome. In conclusion, the data suggest that the soil-living collembolan F. candida has assimilated the capacity for antibacterial activity by horizontal gene transfer, which may be an important adaptive trait in the microbe-dominated soil ecosystem.

Comparative genomics of the mimicry switch in Papilio dardanus.

The African Mocker Swallowtail, Papilio dardanus, is a textbook example in evolutionary genetics. Classical breeding experiments have shown that wing pattern variation in this polymorphic Batesian mimic is determined by the polyallelic H locus that controls a set of distinct mimetic phenotypes. Using bacterial artificial chromosome (BAC) sequencing, recombination analyses and comparative genomics, we show that H co-segregates with an interval of less than 500 kb that is collinear with two other Lepidoptera genomes and contains 24 genes, including the transcription factor genes engrailed (en) and invected (inv). H is located in a region of conserved gene order, which argues against any role for genomic translocations in the evolution of a hypothesized multi-gene mimicry locus. Natural populations of P. dardanus show significant associations of specific morphs with single nucleotide polymorphisms (SNPs), centred on en. In addition, SNP variation in the H region reveals evidence of non-neutral molecular evolution in the en gene alone. We find evidence for a duplication potentially driving physical constraints on recombination in the lamborni morph. Absence of perfect linkage disequilibrium between different genes in the other morphs suggests that H is limited to nucleotide positions in the regulatory and coding regions of en. Our results therefore support the hypothesis that a single gene underlies wing pattern variation in P. dardanus.

Towards a mitogenomic phylogeny of Lepidoptera.

The backbone phylogeny of Lepidoptera remains unresolved, despite strenuous recent morphological and molecular efforts. Molecular studies have focused on nuclear protein coding genes, sometimes adding a single mitochondrial gene. Recent advances in sequencing technology have, however, made acquisition of entire mitochondrial genomes both practical and economically viable. Prior phylogenetic studies utilised just eight of 43 currently recognised lepidopteran superfamilies. Here, we add 23 full and six partial mitochondrial genomes (comprising 22 superfamilies of which 16 are newly represented) to those publically available for a total of 24 superfamilies and ask whether such a sample can resolve deeper lepidopteran phylogeny. Using recoded datasets we obtain topologies that are highly congruent with prior nuclear and/or morphological studies. Our study shows support for an expanded Obtectomera including Gelechioidea, Thyridoidea, plume moths (Alucitoidea and Pterophoroidea; possibly along with Epermenioidea), Papilionoidea, Pyraloidea, Mimallonoidea and Macroheterocera. Regarding other controversially positioned higher taxa, Doidae is supported within the new concept of Drepanoidea and Mimallonidae sister to (or part of) Macroheterocera, while among Nymphalidae butterflies, Danainae and not Libytheinae are sister to the remainder of the family. At the deepest level, we suggest that a tRNA rearrangement occurred at a node between Adeloidea and Ditrysia+Palaephatidae+Tischeriidae.

Mitogenome sequences stabilize the phylogenetics of weevils (Curculionoidea) and establish the monophyly of larval ectophagy.

The weevils and their relatives (superfamily Curculionoidea) constitute a huge radiation of beetles, but basal relationships in this group remain controversial, in particular within the largest family, Curculionidae. We used next-generation sequencing to generate mitochondrial genome data comprising 12 protein coding genes for 27 taxa. Together with two published ingroup sequences, (sub)family relationships were assessed with Bayesian and ML searches under various models and partitioning schemes. Forward and reverse strands of the mitochondrial genome differed in nucleotide skew. Consequently synonymous codon usage differed substantially on either strand for each amino acid, whereby codons ending in AT and GC were favored on the forward and reverse strands, respectively. Data partitioning by forward/reverse strand and codon position greatly improved likelihood scores and nodal support, whereas the tree topology was largely stable. The analysis generally supports the basal position of several 'orthocerus' lineages with straight antennae and male genitalia of an ancestral type; a paraphyletic mixed group of Heteromorphi exhibiting mixed ancestral and derived antennal and genitalic characters; and the derived 'gonathocerous' lineages with kinked antennae corresponding to the strongly supported Curculionidae. The latter did not include the wood boring Platypodinae that was recovered as sister to Dryophthoridae, while the Scolytinae and Cossoninae formed two independent lineages of wood borers within Curculionidae. A basal split in Curculionidae placed the Entiminae and Hyperinae as sister to all other subfamilies with high support, which provides a new ecological concept for structuring the Curculionidae according to the ectophagous larval life style in the former versus endophagous larvae in all others. This basal split is also supported by gene order rearrangements in a tRNA cluster. Recent studies supporting the monophyly of wood boring weevils may be attributable to long-branch attraction, as molecular rates in their mitochondrial genomes were found to be higher than in other lineages, but this did not confound tree searches under combined analysis of mitochondrial protein coding genes.

Phylogenetically informative rearrangements in mitochondrial genomes of Coleoptera, and monophyly of aquatic elateriform beetles (Dryopoidea).

Mitochondrial gene order in Coleoptera has been thought to be conservative but a survey of 60 complete or nearly complete genomes revealed a total of seven different gene rearrangements (deletions, gene order reversals), mainly affecting tRNA genes. All of these were found to be limited to a single taxon or a subclade of Coleoptera. The phylogenetic distribution of a translocation of tRNA(Pro) in three species of elateriform beetles was investigated further by sequencing three nearly complete mitochondrial genomes (Dascillidae, Byrrhidae, Limnichidae) and ten additional individuals for a ∼1370 bp diagnostic fragment spanning the relevant region. Phylogenetic analysis consistently recovered the monophyly of families previously grouped in the contentious superfamily Dryopoidea, a group of approximately 10 beetle families with mainly aquatic lifestyles. The Byrrhidae (moss beetles) were not part of this lineage, although they may be its sister group, to recover the widely accepted Byrrhoidea. The tRNA(Pro) translocation was present in all members of Dryopoidea, but not in any other Elateriformia, providing independent support for this lineage and for a single origin of aquatic habits.

High throughput nano-liter RT-qPCR to classify soil contamination using a soil arthropod.

BACKGROUND: To incorporate genomics data into environmental assessments a mechanistic perspective of interactions between chemicals and induced biological processes needs to be developed. Since chemical compounds with structural similarity often induce comparable biological responses in exposed animals, gene expression signatures can serve as a starting point for the assessment of chemicals and their toxicity, but only when relevant and stable gene panels are available. To design such a panel, we isolated differentially expressed gene fragments from the soil arthropod Folsomia candida, a species often used for ecotoxicological testing. Animals were exposed to two chemically distinct compounds, being a metal (cadmium) and a polycyclic aromatic hydrocarbon (phenanthrene). We investigated the affected molecular responses resulting from either treatment and developed and validated 44 qPCR assays for their responses using a high throughput nano-liter RT-qPCR platform for the analysis of the samples. RESULTS: Suppressive subtractive hybridization (SSH) was used to retrieve stress-related gene fragments. SSH libraries revealed pathways involved in mitochondrial dysfunction and protein degradation for cadmium and biotransformation for phenanthrene to be overrepresented. Amongst a small cluster of SSH-derived cadmium responsive markers were an inflammatory response protein and an endo-glucanase. Conversely, cytochrome P450 family 6 or 9 was specifically induced by phenanthrene. Differential expressions of these candidate biomarkers were also highly significant in the independently generated test sample set. Toxicity levels in different training samples were not reflected by any of the markers' intensity of expressions. Though, a model based on partial least squares differential analysis (PLS-DA) (with RMSEPs between 9 and 22% and R(2)s between 0.82 and 0.97) using gene expressions of 25 important qPCR assays correctly predicted the nature of exposures of test samples. CONCLUSIONS: For the application of molecular bio-indication in environmental assessments, multivariate analyses obviously have an added value over univariate methods. Our results suggest that compound discrimination can be achieved by PLS-DA, based on a hard classification of the within-class rankings of samples from a test set. This study clearly shows that the use of high throughput RT-qPCR could be a valuable tool in ecotoxicology combining high throughput with analytical sensitivity.

Transcriptional plasticity of a soil arthropod across different ecological conditions.

Ecological functional genomics, dealing with the responses of organisms to their natural environment is confronted with a complex pattern of variation and a large number of confounding environmental factors. For gene expression studies to provide meaningful information on conditions deviating from normal, a baseline or normal operating range (NOR) response needs to be established which indicates how an organism's transcriptome reacts to naturally varying ecological factors. Here we determine the transcriptional plasticity of a soil arthropod, Folsomia candida, exposed to various natural environments, as part of a first attempt in establishing such a NOR. Animals were exposed to 26 different field soils after which gene expression levels were measured. The main factor found to regulate gene expression was soil-type (sand or clay). Cell homeostasis and DNA replication were affected in collembolans exposed to sandy soil, indicating general stress. Multivariate analysis identified soil fertility as the main factor influencing gene expression. Regarding land-use, only forest soils showed an expression pattern deviating from the others. No significant effect of land-use, agricultural practice or soil type on fitness was observed, but arsenic concentration was negatively correlated with reproductive output. In conclusion, transcriptional responses remained within a limited range across the different land-uses but were significantly affected by soil-type. This may be caused by the contrasting soil physicochemical properties to which F. candida strongly responds. The broad range of conditions over which this soil-living detritivore is able to survive and reproduce, indicates a strategy of high plasticity, which comes with extensive gene expression regulation.

Reference genes for QRT-PCR tested under various stress conditions in Folsomia candida and Orchesella cincta (Insecta, Collembola).

BACKGROUND: Genomic studies measuring transcriptional responses to changing environments and stress currently make their way into the field of evolutionary ecology and ecotoxicology. To investigate a small to medium number of genes or to confirm large scale microarray studies, Quantitative Reverse Transcriptase PCR (QRT-PCR) can achieve high accuracy of quantification when key standards, such as normalization, are carefully set. In this study, we validated potential reference genes for their use as endogenous controls under different chemical and physical stresses in two species of soil-living Collembola, Folsomia candida and Orchesella cincta. Treatments for F. candida were cadmium exposure, phenanthrene exposure, desiccation, heat shock and pH stress, and for O. cincta cadmium, desiccation, heat shock and starvation. RESULTS: Eight potential reference genes for F. candida and seven for O. cincta were ranked by their stability per stress factor using the programs geNorm and Normfinder. For F. candida the succinate dehydrogenase (SDHA) and eukaryotic transcription initiation factor 1A (ETIF) genes were found the most stable over the different treatments, while for O. cincta, the beta actin (ACTb) and tyrosine 3-monooxygenase (YWHAZ) genes were the most stable. CONCLUSION: We present a panel of reference genes for two emerging ecological genomic model species tested under a variety of treatments. Within each species, different treatments resulted in differences in the top stable reference genes. Moreover, the two species differed in suitable reference genes even when exposed to similar stresses. This might be attributed to dissimilarity of physiology. It is vital to rigorously test a panel of reference genes for each species and treatment, in advance of relative quantification of QRT-PCR gene expression measurements.

Adaptive differences in gene expression associated with heavy metal tolerance in the soil arthropod Orchesella cincta.

Field-selected tolerance to heavy metals has been reported for Orchesella cincta (Arthropoda: Collembola) populations occurring at metal-contaminated mining sites. This tolerance correlated with heritable increase in metal excretion efficiency, less pronounced cadmium (Cd)-induced growth reduction and overexpression of the metallothionein gene. We applied transcriptomics to determine differential gene expression caused by this abiotic stress in reference and Cd-tolerant populations. Many cDNAs responded to Cd exposure in the reference population. Significantly fewer clones were Cd responsive in tolerant animals. Analysis of variance revealed transcripts that interact between Cd exposure and population. Hierarchical cluster analysis of these clones identified two major groups. The first one contained cDNAs that were up-regulated by Cd in the reference culture but non-responsive or down-regulated in tolerant animals. This cluster was also characterized by elevated constitutive expression in the tolerant population. Gene ontology analysis revealed that these cDNAs were involved in structural integrity of the cuticle, anti-microbial defence, calcium channel-blocking, sulphur assimilation and chromatin remodelling. The second group consisted of cDNAs down-regulated in reference animals but not responding or slightly up-regulated in tolerant animals. Their functions involved carbohydrate metabolic processes, Ca(2+)-dependent stress signalling, redox state, proteolysis and digestion. The reference population showed a strong signature of stress-induced genome-wide perturbation of gene expression, whereas the tolerant animals maintained normal gene expression upon Cd exposure. We confirmed the micro-evolutionary processes occurring in soil arthropod populations and suggest a major contribution of gene regulation to the evolution of a stress-adapted phenotype.

Collembase: a repository for springtail genomics and soil quality assessment.

BACKGROUND: Environmental quality assessment is traditionally based on responses of reproduction and survival of indicator organisms. For soil assessment the springtail Folsomia candida (Collembola) is an accepted standard test organism. We argue that environmental quality assessment using gene expression profiles of indicator organisms exposed to test substrates is more sensitive, more toxicant specific and significantly faster than current risk assessment methods. To apply this species as a genomic model for soil quality testing we conducted an EST sequencing project and developed an online database. DESCRIPTION: Collembase is a web-accessible database comprising springtail (F. candida) genomic data. Presently, the database contains information on 8686 ESTs that are assembled into 5952 unique gene objects. Of those gene objects approximately 40% showed homology to other protein sequences available in GenBank (blastx analysis; non-redundant (nr) database; expect-value < 10-5). Software was applied to infer protein sequences. The putative peptides, which had an average length of 115 amino-acids (ranging between 23 and 440) were annotated with Gene Ontology (GO) terms. In total 1025 peptides (approximately 17% of the gene objects) were assigned at least one GO term (expect-value < 10-25). Within Collembase searches can be conducted based on BLAST and GO annotation, cluster name or using a BLAST server. The system furthermore enables easy sequence retrieval for functional genomic and Quantitative-PCR experiments. Sequences are submitted to GenBank (Accession numbers: EV473060 - EV481745). CONCLUSION: Collembase http://www.collembase.org is a resource of sequence data on the springtail F. candida. The information within the database will be linked to a custom made microarray, based on the Agilent platform, which can be applied for soil quality testing. In addition, Collembase supplies information that is valuable for related scientific disciplines such as molecular ecology, ecogenomics, molecular evolution and phylogenetics.

Uncovering Trophic Interactions in Arthropod Predators through DNA Shotgun-Sequencing of Gut Contents.

Characterizing trophic networks is fundamental to many questions in ecology, but this typically requires painstaking efforts, especially to identify the diet of small generalist predators. Several attempts have been devoted to develop suitable molecular tools to determine predatory trophic interactions through gut content analysis, and the challenge has been to achieve simultaneously high taxonomic breadth and resolution. General and practical methods are still needed, preferably independent of PCR amplification of barcodes, to recover a broader range of interactions. Here we applied shotgun-sequencing of the DNA from arthropod predator gut contents, extracted from four common coccinellid and dermapteran predators co-occurring in an agroecosystem in Brazil. By matching unassembled reads against six DNA reference databases obtained from public databases and newly assembled mitogenomes, and filtering for high overlap length and identity, we identified prey and other foreign DNA in the predator guts. Good taxonomic breadth and resolution was achieved (93% of prey identified to species or genus), but with low recovery of matching reads. Two to nine trophic interactions were found for these predators, some of which were only inferred by the presence of parasitoids and components of the microbiome known to be associated with aphid prey. Intraguild predation was also found, including among closely related ladybird species. Uncertainty arises from the lack of comprehensive reference databases and reliance on low numbers of matching reads accentuating the risk of false positives. We discuss caveats and some future prospects that could improve the use of direct DNA shotgun-sequencing to characterize arthropod trophic networks.

Metagenome Skimming of Insect Specimen Pools: Potential for Comparative Genomics.

Metagenomic analyses are challenging in metazoans, but high-copy number and repeat regions can be assembled from low-coverage sequencing by "genome skimming," which is applied here as a new way of characterizing metagenomes obtained in an ecological or taxonomic context. Illumina shotgun sequencing on two pools of Coleoptera (beetles) of approximately 200 species each were assembled into tens of thousands of scaffolds. Repeated low-coverage sequencing recovered similar scaffold sets consistently, although approximately 70% of scaffolds could not be identified against existing genome databases. Identifiable scaffolds included mitochondrial DNA, conserved sequences with hits to expressed sequence tag and protein databases, and known repeat elements of high and low complexity, including numerous copies of rRNA and histone genes. Assemblies of histones captured a diversity of gene order and primary sequence in Coleoptera. Scaffolds with similarity to multiple sites in available coleopteran genome sequences for Dendroctonus and Tribolium revealed high specificity of scaffolds to either of these genomes, in particular for high-copy number repeats. Numerous "clusters" of scaffolds mapped to the same genomic site revealed intra- and/or intergenomic variation within a metagenome pool. In addition to effect of taxonomic composition of the metagenomes, the number of mapped scaffolds also revealed structural differences between the two reference genomes, although the significance of this striking finding remains unclear. Finally, apparently exogenous sequences were recovered, including potential food plants, fungal pathogens, and bacterial symbionts. The "metagenome skimming" approach is useful for capturing the genomic diversity of poorly studied, species-rich lineages and opens new prospects in environmental genomics.

Family-Level Sampling of Mitochondrial Genomes in Coleoptera: Compositional Heterogeneity and Phylogenetics.

Mitochondrial genomes are readily sequenced with recent technology and thus evolutionary lineages can be densely sampled. This permits better phylogenetic estimates and assessment of potential biases resulting from heterogeneity in nucleotide composition and rate of change. We gathered 245 mitochondrial sequences for the Coleoptera representing all 4 suborders, 15 superfamilies of Polyphaga, and altogether 97 families, including 159 newly sequenced full or partial mitogenomes. Compositional heterogeneity greatly affected 3rd codon positions, and to a lesser extent the 1st and 2nd positions, even after RY coding. Heterogeneity also affected the encoded protein sequence, in particular in the nad2, nad4, nad5, and nad6 genes. Credible tree topologies were obtained with the nhPhyML ("nonhomogeneous") algorithm implementing a model for branch-specific equilibrium frequencies. Likelihood searches using RAxML were improved by data partitioning by gene and codon position. Finally, the PhyloBayes software, which allows different substitution processes for amino acid replacement at various sites, produced a tree that best matched known higher level taxa and defined basal relationships in Coleoptera. After rooting with Neuropterida outgroups, suborder relationships were resolved as (Polyphaga (Myxophaga (Archostemata + Adephaga))). The infraorder relationships in Polyphaga were (Scirtiformia (Elateriformia ((Staphyliniformia + Scarabaeiformia) (Bostrichiformia (Cucujiformia))))). Polyphagan superfamilies were recovered as monophyla except Staphylinoidea (paraphyletic for Scarabaeiformia) and Cucujoidea, which can no longer be considered a valid taxon. The study shows that, although compositional heterogeneity is not universal, it cannot be eliminated for some mitochondrial genes, but dense taxon sampling and the use of appropriate Bayesian analyses can still produce robust phylogenetic trees.