Martijn J.T.N. Timmermans

Genetic variation in toxicant-stressed populations: an evaluation of the "genetic erosion" hypothesis.

The question whether environmental pollution affects genetic diversity in natural populations remains unanswered to date despite the fact that genetic variation is one of the three pillars of biodiversity recognized in the Rio convention of 1993. The loss of genetic diversity in populations subjected to anthropogenic stress can be designated as “genetic erosion” and may be considered as a factor of concern in risk assessment of toxic chemicals. Theoretically there are four different ways in which toxicants can affect genetic variation: (i) by increasing mutation rates, (ii) by directional selection on tolerant genotypes, (iii) by causing bottleneck events, and (iv) by altering migration. This paper reviews studies that have documented genetic change in animal populations exposed to environmental pollution. In these studies, genetic variation is measured in a variety of ways: heritability of quantitative characters, heterozygosity of allozyme loci, haplotype diversity in mitochondrial DNA, and variability in RAPD fingerprints. Studies on cadmium tolerance of Collembola living in metal-contaminated soil suggest that strong directional selection pressure may decrease genetic variability of traits immediately linked to tolerance. Allozyme studies in fish have documented a similar decrease of genetic variation in populations living in strongly acidified waters. A correlation between RAPD-PCR-based genetic similarity and site contamination has been documented in crayfish. Overall, there is significant support for the genetic erosion hypothesis, but the issue cannot be considered settled. In most studies insufficient attention is given to factors such as population size, bottlenecks and mutation, which may influence genetic variability in addition to the toxicant selection regime. At the moment, there does not seem to be a sound scientific basis for incorporating genetic diversity measurements into risk assessment, despite the variety of easily applicable molecular techniques available. It is often not known what kind of variation is measured by these techniques (neutral or selectable) and how the markers are inherited. Given the importance of the issue, as stressed by the Rio Convention, a concentrated research effort is necessary to better define the question and find a general approach to evaluate its importance in ecological risk assessment.

Why barcode? High-throughput multiplex sequencing of mitochondrial genomes for molecular systematics

Mitochondrial genome sequences are important markers for phylogenetics but taxon sampling remains sporadic because of the great effort and cost required to acquire full-length sequences. Here, we demonstrate a simple, cost-effective way to sequence the full complement of protein coding mitochondrial genes from pooled samples using the 454/Roche platform. Multiplexing was achieved without the need for expensive indexing tags (‘barcodes’). The method was trialled with a set of long-range polymerase chain reaction (PCR) fragments from 30 species of Coleoptera (beetles) sequenced in a 1/16th sector of a sequencing plate. Long contigs were produced from the pooled sequences with sequencing depths ranging from ∼10 to 100× per contig. Species identity of individual contigs was established via three ‘bait’ sequences matching disparate parts of the mitochondrial genome obtained by conventional PCR and Sanger sequencing. This proved that assembly of contigs from the sequencing pool was correct. Our study produced sequences for 21 nearly complete and seven partial sets of protein coding mitochondrial genes. Combined with existing sequences for 25 taxa, an improved estimate of basal relationships in Coleoptera was obtained. The procedure could be employed routinely for mitochondrial genome sequencing at the species level, to provide improved species ‘barcodes’ that currently use the cox1 gene only.

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.

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.

Genetic structure in Orchesella cincta (Collembola): strong subdivision of European populations inferred from mtDNA and AFLP markers.

Population genetic structure is determined both by current processes and historical events. Current processes include gene flow, which is largely influenced by the migration capacity of a species. Historical events are, for example, glaciation periods, which have had a major impact on the distribution of many species. Species with a low capacity or tendency to move about or disperse often exhibit clear spatial genetic structures, whereas mobile species mostly show less spatial genetic differentiation. In this paper we report on the genetic structure of a small, wingless arthropod species (Orchesella cincta: Collembola) in Europe. For this purpose we used mtDNA COII sequences and AFLP markers. We show that large genetic differences exist between populations of O. cincta, as expected from O. cinctas winglessness and sedentary lifestyle. Despite the fact that most variability was observed within populations (59%), a highly significant amount of AFLP variation (25%) was observed between populations from northwestern Europe, central Europe and Italy. This suggests that gene flow among regions is extremely low, which is additionally supported by the lack of shared mtDNA alleles between regions. Based on the genetic variation and sequence differences observed we conclude that the subdivision occurred long before the last glaciation periods. Although the populations still interbreed in the lab, we assume that in the long term the genetic isolation of these regions may lead to speciation processes.

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

BackgroundGenomic 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.ResultsEight 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.ConclusionWe 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, Ca2+‐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.

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.

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

BackgroundEnvironmental 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.DescriptionCollembase 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 ~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 (~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).ConclusionCollembase 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.

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.