Single nucleotide polymorphism discovery from expressed sequence tags in the waterflea Daphnia magna
- Luisa Orsini1Email author,
- Mieke Jansen1Email author,
- Erika L Souche†2, 3,
- Sarah Geldof1, 2 and
- Luc De Meester1
© Orsini et al; licensee BioMed Central Ltd. 2011
Received: 21 June 2010
Accepted: 13 June 2011
Published: 13 June 2011
Daphnia (Crustacea: Cladocera) plays a central role in standing aquatic ecosystems, has a well known ecology and is widely used in population studies and environmental risk assessments. Daphnia magna is, especially in Europe, intensively used to study stress responses of natural populations to pollutants, climate change, and antagonistic interactions with predators and parasites, which have all been demonstrated to induce micro-evolutionary and adaptive responses. Although its ecology and evolutionary biology is intensively studied, little is known on the functional genomics underpinning of phenotypic responses to environmental stressors. The aim of the present study was to find genes expressed in presence of environmental stressors, and target such genes for single nucleotide polymorphic (SNP) marker development.
We developed three expressed sequence tag (EST) libraries using clonal lineages of D. magna exposed to ecological stressors, namely fish predation, parasite infection and pesticide exposure. We used these newly developed ESTs and other Daphnia ESTs retrieved from NCBI GeneBank to mine for SNP markers targeting synonymous as well as non synonymous genetic variation. We validate the developed SNPs in six natural populations of D. magna distributed at regional scale.
A large proportion (47%) of the produced ESTs are Daphnia lineage specific genes, which are potentially involved in responses to environmental stress rather than to general cellular functions and metabolic activities, or reflect the arthropod's aquatic lifestyle. The characterization of genes expressed under stress and the validation of their SNPs for population genetic study is important for identifying ecologically responsive genes in D. magna.
Single Nucleotide Polymorphisms (SNPs), defined as single-base changes or short insertion or deletion mutations (indels), are the most abundant class of genetic variation found in eukaryotic genomes. SNPs are widespread, and present in both coding and non coding regions [1–3]. Until few years ago, the use of SNP markers was limited to model organisms with sequenced genomes, mostly because of the costs associated with SNP discovery. Methods for indirect SNP discovery detect heteroduplexes on the basis of mismatch-induced altered DNA characteristics [4–8], whereas typical direct SNP discovery strategies [9, 10] involve sequencing of locus-specific amplification (LSA) products from multiple individuals. One of the most common strategies adopted for SNP development in non-model organisms is the use of Expressed Sequence Tags (ESTs) as a resource for SNP marker detection [11–14]. This is a consequence of the increasing availability of EST libraries for non-model organisms. As a result, SNPs are becoming increasingly important for research on non-model organisms. SNPs offer the potential for genome wide scans of selectively neutral as well as adaptive variation [15, 16], with simple mutation models and powerful analytical methods , and with application to noninvasive analysis and historical DNA .
Most genetic surveys of natural populations focus on neutral loci. Whereas this provides valuable insights into the historical demography and evolution of populations (see ), it does not allow to understand the dynamics of genes that affect fitness along environmental gradients. Daphnia (Crustacea: Cladocera) offers a unique opportunity to study neutral and selective variation in natural populations with a known ecological background. Daphnia play a pivotal role in the ecology of standing waters, are widely used in population studies and environmental risk assessments, and are supported by a large community of ecologists, evolutionary biologists and ecotoxicologists [19–21]. This is in part due to their ease of culture, convenient size, short generation time and cyclic parthenogenetic reproduction, which make them very suitable for laboratory and field experiments, experimental evolution, and quantitative genetic analyses in multiple environments. Thanks to the sustained efforts of the Daphnia Genomics Consortium (http://daphnia.cgb.indiana.edu, DGC), Daphnia is also regarded as a model invertebrate in ecological genomics [21, 22]. Daphnia magna is, especially in Europe, intensively used to study stress responses to pollutants, climate change, and antagonistic interactions with predators and parasites [23–25]. Daphnia has also been subject to intensive population genetic study, with a strong focus on the impact of its peculiar reproduction mode, cyclical parthenogenesis, population genetic structure and among-population genetic differentiation [26–29]. Yet, although some knowledge has recently been acquired in the fields of functional responses to parasite infection  and proteomics , precious little is known about the complex interaction between neutral genetic variation, reflecting population genetic structure and demography, and the functional genomics underpinning phenotypic responses to environmental stressors. Among other reasons, the lack of suitable markers for functional traits has been one of the main limitations.
We developed three EST libraries using clonal lineages of D. magna exposed to standardized selection pressures, namely fish predation, exposure to parasites (Pasteuria ramosa), and exposure to pesticides (carbaryl). These environmental stressors are known to induce pronounced micro-evolutionary responses in D. magna[23–25]. Our EST sequences as well as EST sequences published in NCBI at the time of the analysis were mined for SNP markers targeting synonymous and non-synonymous polymorphisms. An in silico discovery tool purposely designed to mine EST sequences (Souche et al, in prep) was used for SNP discovery. The newly developed SNPs were validated by genotyping individuals from six natural populations of D. magna distributed at regional scale. The identification and characterization of genes differentially expressed in stress conditions and the validation of SNP mutations that could be linked to specific environmental stressors opens new interesting perspectives in the study of functional polymorphism in natural populations of D. magna.
Development of EST libraries
Exposure of clones to environmental stressors
The clones used to develop the EST libraries were two genetic isolates (M10 and Mu11) exposed independently to three environmental stressors that are known to induce micro-evolutionary responses in D. magna: fish predation, parasite infection and pesticide exposure. Clone M10 was hatched from Oud-Heverlee pond in Belgium (see  for information on habitat), whereas clone Mu11 was isolated from a pond in Germany (see ). Animals were grown for two generations under standard light/food regime to eliminate residual maternal effects after hatching. Clone Mu11 was exposed to fish kairomones enriched medium, while clone M10 was exposed independently to parasites and pesticides. All exposed animals were grown in a climate chamber (20+/-1°C) in 100 ml ADaM medium  and daily fed with 1 × 105 cells of Scenedesmus obliquus/ml. Exposure to the environmental stressors was performed on juveniles not older than 24 h released from the second clutch. Juveniles were exposed in a concentration of 10 animals/50 ml, daily fed with 1 × 108 cells of S. obliquus/ml, and kept in daily refreshed medium.
The exposure to fish predation was mimicked by culturing D. magna in fish kairomone enriched ADaM medium . The medium was obtained by filtering (0.450 μm mesh) and diluting water (5 times) from a 20L aquarium where three gold orfes (Leuciscus idus) were kept for 24 h.
The common endoparasite of D. magna, Pasteuria ramosa, was used to induce parasite infection. P. ramosa infection happens via horizontal transmission of spores released from dead infected hosts. Monoclonal cultures of the clone M10 were exposed to a spore solution consisting of squashed and filtered (60 μm nylon filter) infected animals. A final concentration of 1 × 106 mature spores/ml was used to infect the Daphnia host. To avoid that adaptation to specific strains of the parasite spores would affect the response to infection, the host and the spores were collected from two distinct ponds in Belgium [hosts were hatched from Oud-Heverlee pond, whereas P. ramosa spores were collected from the sediment of OM2 pond (for details on the ponds see ). A detailed description of the infection method is reported elsewhere .
We used carbaryl (1-napthyl methylcarbarmate), a commonly used pesticides in agriculture, as a model pesticide. Different concentrations of carbaryl (CAS 63-25-2, purity 99.8%, Sigma-Aldrich, Germany) were used: 5.6 μg/L, 8.0 μg/L and 11.4 μg/L. These concentrations were determined by earlier work on pesticide exposure to be sublethal but at the same time affect Daphnia fitness [35, 36].
After exposure, animals were transferred to liquid nitrogen for storage. Juveniles exposed to fish predation and parasite infection were stored after 48 h and 96 h. Carbaryl exposed animals were stored after 48 h, 96 h and 144 h exposure time. Depending on the age and the survival of the exposed Daphnia, between 18 and 53 individuals were pooled for RNA extraction.
Construction of cDNA libraries and sequencing of ESTs
Total RNA was extracted from the pool of exposed and non-exposed (control) samples using the Trizol® extraction method (Invitrogen, Life technologies, Belgium), following the manufacturer's instructions. After extraction and DNAse treatment (Fermentas, Germany; ), RNA purity was checked with a spectrophotometer (NanoDrop Technologies, USA).
Pooled RNA samples (three time points -48, 96 and 144 h- for the carbaryl exposures and two -48 and 96 h- for parasite and fish kairomone exposures) from individual stressors were used to create three separate libraries. Each library represents the genes differentially expressed in a pool of exposed individuals versus a control (non exposed individuals). The construction of the subtractive libraries was done according to Soetaert et al. . In brief, we used a combination of the SMART™ PCR cDNA synthesis kit (Clontech, USA) and the PCR-Select™ cDNA subtraction kit (Clontech, USA) to obtain libraries by means of Suppression Subtractive Hybridization . Differently expressed genes in the three subtraction libraries were cloned using the pGEM®-T Easy Vector System II (Promega, USA) following the manufacturer's instructions. The gene fragments were PCR amplified with the cloning vectors, M13F (5' CGA CTG TGT AAA ACG ACG GCC AG 3') and M13R (5'CAG GAA ACA GCT ATG ACC ATG ATT ACG CC 3'). The PCR products were purified with Exo-SAP (Fermentas, USA) at a concentration of 10 u/5 μl of PCR product prior sequencing. PCR products were single-strand sequenced using Big Dye terminator chemistry on a CEQ™ 8000 automated sequencer (Beckman Coulter). The PCR fragments were sequenced with the vector primers SP6 (5' ATT TAG GTG ACA CTA TAG 3') and T7 (5' TAA TAC GAC TCA CTA TAG GG 3). Sequence assembly and editing were done using CodonCode Aligner http://www.codoncode.com.
EST assembly, annotation and processing for SNP development
cDNA fragments were trimmed from residual vector sequences using CodonCode Aligner. Low quality and short sequences fragments (fragments <50 bp) were removed from the total set of sequences. Base calling and attribution of quality values were performed using Phred , using default parameters. The software cross_match (http://www.phrap.org/phredphrapconsed.html#block_phrap) was used to mask oligonucleotides, primers and adapters used during the library construction on the remaining ESTs. Masked sequences, polyA tails and low complexity sequences were trimmed with SeqClean using default parameters and a minimum accepted length of 50bp (http://compbio.dfci.harvard.edu/tgi/software/). The trimmed fragments and EST sequences obtained from NCBI GenBank at the time of the analysis  (Additional file 1) were assembled in clusters based on their similarities using TGICL  and default parameters, except for the overlap percent identity cut-off that was set to 93 instead of 80. A quality value of 20 was arbitrarily assigned to all EST nucleotides downloaded from NCBI GenBank in order to process sequences with unknown quality in the SNP pipeline (see below). The group of ESTs assembled in one consensus sequence and the remaining unclustered ESTs (referred hereafter as contigs and singlets) were used for sequence similarity searches against the NCBI and the http://wfleabase.org/ genome databases using Blast . These contigs and singlets were subsequently annotated using Open Reading Frame (ORF) prediction by OrfPredictor , Gene Ontology assignments by Blast2GO  and Prosite functional protein domains by ScanProsite . Orthologs and paralogs were identified using the OrthoMCL database that includes 1,270,853 sequences from 138 genomes .
SNP markers development
The clusters obtained by assembling EST sequences from GenBank and our ESTs were mined for SNPs using a pipeline that integrates six freely available SNP discovery tools: SNPserver , PolyBayes , Quality SNP , PolyFreq , MiraEST , and ssahaSNP . Stringent quality criteria were used for SNP mining. We excluded clusters containing singlet sequences for obvious reasons; we considered only one contig per cluster when multiple contigs were present to avoid redundancy; a mismatch was regarded as SNP candidate if it appeared at least twice in the alignment. Consequently, only the contigs containing four overlapping ESTs were searched for SNPs. The ACE file produced by TGICL was used as input of the pipeline. We also designed SNPs probes and PCR primers for five nuclear genes, including one nuclear receptor gene [Ultraspiracle (Usp)] and four metabolic enzyme genes [Phosphoglucose isomerase (Pgi), Mannose-phosphate isomerase (Mpi), l-lactate dehydrogenase (Ldh), Enolase (Enol)]. These genes were previously used in a study involving D. magna and D. pulex. The sequences were kindly provided by Christoph R. Haag.
SNP validation and characterization
SNP loci were validated by genotyping natural populations of D. magna. We sampled the surface sediment layer of six populations (185 dormant eggs) of D. magna, scattered across Belgium (three from the coast and three from the inland, Additional file 2). By hatching dormant eggs from the superficial sediment layers, we ensured to capture the genetic diversity of the populations without it being affected by clonal erosion typically occurring in the active populations of cyclical parthenogenetic species [55, 56]. Our approach aims at sampling genetic diversity as it occurs at hatching, at the beginning of the growing season. Genotyping of SNPs was performed using the MassARRAY platform from Sequenom (San Diego, USA) at the VIB Vesalius Research Center, Leuven, Belgium. Genomic positions of the genetic variants were selected and 70 bp-up and downstream sequences were used for primer design with the Sequenom MassARRAY Assay Design 3.1 software with default parameters. Respective forward PCR, reverse PCR and extension primers for the Sequenom genotyping assays can be found in supplementary material (Additional file 3). Multiplexes levels were between 17 and 40. The genotyping was performed according to the iPLEX protocol from Sequenom (available at http://www.sequenom.com/). Quality control criteria were adopted (water as negative control and inter-plate duplicates testing).
SNP markers were characterized by counting the number of synonymous and non-synonymous mutations, and by identifying the codon position responsible for the non-synonymous changes.
Population genetics analysis
A test for outlier SNPs was performed to identify the strictly neutral SNPs to use in population genetic statistics in order to validate the newly developed markers. The outlier loci were identified with two methods, Lositan  (a selection detection workbench constructed around the program fdist ) and BAYESCAN . The simulation conditions were as follows. For Lositan , we used 10,000 simulation replicates (infinite allele model, IAM) for a sample size of 50 individuals, with forced mean FST calculated on the real data set. Multiple runs (generally 3) were run to avoid spurious results. For BAYESCAN, 10 pilot runs of 5,000 iterations and an additional burn-in of 50,000 iterations were followed by 100,000 iterations (sample size of 5,000 and thinning interval of 20) to identify loci under selection from locus specific Bayes factors. A Bayes factor of 3, corresponding to a posterior probability of 0.5 (substantial selection) was considered as the minimal threshold for a locus to be considered under selection. Multiple runs (generally 3) were run to avoid spurious results.
The strictly neutral loci were used to estimate the among-population differentiation with the Microsatellite Analyser (MSA) software v 3.12 . They were also used in an analysis of population genetic structure and inference of population demography. The spatial genetic structure was described using the Bayesian analysis implemented in the BAPS software [61–63]. We performed a non-spatial genetic mixture analysis using individuals and populations as basic units to be clustered. Different starting values of K (2 to 12) were used to verify the robustness of the results. A Mantel test (10,000 permutations) was performed using IBD online software  to test for Isolation by Distance.
The three EST libraries yielded 1,698 gene fragments. Of these fragments, 480 were from the parasite library, 480 from the fish library and 738 from the pesticide library. After a quality trimming, 1,070 fragments remained (301 from the fish library, 366 from the parasite library and 403 were from the pesticide library; GenBank accession numbers: HO045245-HO046616). Of the total gene fragments, 685 showed significant similarity (Blast e < 10-5) with genes within the NCBI non-redundant protein database (64% of the total sequenced fragments, Additional file 4). Of these, 4 matched nuclear eukaryotic genes, 16 matched mitochondrial genes, 82 matched ribosomal genes, and 196 showed homology to known proteins in other organisms (Additional file 4). An additional 60 ESTs aligned to annotated D. pulex and D. magna genes, and 327 ESTs aligned to predicted or conserved hypothetical proteins. Of the 685 gene fragments, 287 have a gene prediction based on homology with known genes in other organism. The annotated gene fragments showed significant similarity with a number of species, including a large proportion of insects (65%) (Additional file 5).
Annotation of assembled ESTs
The sequences considered for SNP mining, including GeneBank and our sequences for assembling into a non-redundant unigene set had an average length of 454 bp and a standard deviation of 137 bp. Over 80% of all ESTs were longer than 350 bp. From the total number of ESTs, 1,674 clusters of sequences were produced. A total of 10,737 ESTs, representing 75.3% of the 14,253 processed ESTs, were assembled in 1,812 contigs, with a redundancy of 75.3%. Only 2,446 ESTs (17.2% of processed ESTs) remained singletons. Of the total number of obtained contigs, 878 (48.5%) contained 2 ESTs and 143 (7.9%) contained more than 10 ESTs. Out of the 1,812 contigs, 55 (3.0%) contained only ESTs produced in this study, 1,483 (81.9%) contained only ESTs downloaded from GenBank, and 274 (15.1%) contained a mix of both ESTs. 1,183 (65%) of these contigs showed significant sequence similarity (Blast e < 10-5) with genes within the NCBI non-redundant protein database. An additional 421 contigs were aligned to annotated D. pulex proteins. Therefore, 88.5% of the chosen contigs are given putative functions based on sequence similarity to annotated proteins from genome of other species; 26% are unique to the lineage leading to Daphnia. The 574 contigs containing more than 4 overlapping sequences were mined for SNPs: 414 (72%) showed significant sequence similarity with genes within NCBI, 122 additional contigs were aligned to annotated D. pulex proteins. Therefore, 93% of the chosen contigs are given putative functional annotations; 23% are unique to the lineage leading to Daphnia.
In vitro development of SNP markers
The contigs containing more than 4 overlapping sequences were mined for SNPs. The discovery tools used in the pipeline identified different numbers of SNPs (AutoSNP: 778; PolyBayes: 757, QualitySNP: 641, PolyFreq: 51; MiraEST: 333) depending on the specific quality parameters. The number of SNP candidates detected per polymorphic contig ranged between 27 and 1, with 96 contigs showing only one candidate SNP, 163 contigs showing between 2 and 9 candidate SNPs and 17 contigs showing more than 10 candidate SNPs. A total of 986 candidate SNPs were discovered. Excluding the candidate SNPs that were contiguous, close to an indel or to the 5' and 3' ends of the ESTs, and the ones that presented high polymorphism in the flanking regions, we had 159 candidate SNPs left, including more than one candidate per EST in some contigs. PCR primers and oligonucleotide probes were designed only for SNPs indentified by a minimum of 3 discovery tools, in order to avoid false positives. Moreover, we targeted only one SNP per contig to use in subsequent genotyping optimization assay. The use of a single candidate SNP per EST was chosen to avoid complications in the subsequent population genetic analyses caused by linkage among SNPs within the same sequence. Excluding the unsuitable candidates based on all the criteria outlined above, we tested 138 SNPs with Sequenom.
SNP validation and characterization
Population genetic structure using neutral SNPs
Pairwise FST values among six natural populations of D. magna
The 12 loci that departed from neutral expectation showing reduced polymorphism were present within ESTs. These EST showed similarity mostly to transmembrane proteins (Additional file 9). Four of the 12 loci were linked to biotic stress in either Daphnia magna or D. pulex. The genomic regions containing the outlier loci will be object of further analysis to indentify functional polymorphism linked to specific environmental stressors.
EST linked SNPs
SNPs are rapidly becoming the marker of choice for many applications in population ecology, evolution and conservation genetics, because of the potential for high genotyping efficiency, data quality, genome-wide coverage and analytical simplicity (e.g. in modeling mutational dynamics, ). However, the use of SNPs in population genetics and genomics is proceeding slower for non-model organisms than for model organisms. One reason is the need of genomics resources for the focal species in order to be able to develop a large panel of SNPs. The second reason is the ascertainment bias due to a limited number of individuals from which the SNP panel is generally developed. In this study we capitalize on rapidly developing genomics resources for Daphnia and produce a first set of SNP markers for population genomic studies. The limitations associated with ascertainment bias are alleviated by using individuals originating from different geographic locations and by using different EST libraries (ESTs produced in this study and ESTs retrieved from NCBI GenBank) from which the SNPs are developed.
In comparison to SNPs developed from genomic sequences, EST-derived SNPs have several advantages [65–67]. Since ESTs are transcribed sequences, EST-derived SNPs are associated with actual genes allowing the use of gene-associated SNPs for mapping and comparative genome studies ( and references therein). EST-derived SNPs are also a rich source of candidate polymorphisms underlying important traits leading to the identification of quantitative trait nucleotides (QTN) (e.g. ) linked to ecologically relevant genes. With this study we have achieved a first step toward the identification of candidate point mutations that can be linked to phenotypic responses to environmental stress in future experimental studies.
Validation of SNPs in population genetic analysis
Our screen for outlier loci identified a large portion (12%) of loci showing reduced polymorphism as compared to neutral expectations and therefore potentially under selection. The proportion of loci showing reduced genetic variation was higher than generally observed in literature (2.6 to 5.5% in humans ; 1.4-3.2% in lake whitefish ecotypes ; 2.6 to 3.3.% in Norway spruce ; 5.5% in white spruce ; 9.5% in salmon ; 1.3-3.6% in common frog ; 5% in sticklebacks ). Most of the mentioned studies are based on anonymous markers or markers falling outside genic regions, whereas our markers are designed within expressed genomic regions. This is likely to have increased the percentage of loci potentially under direct effect of selection in our analysis. In addition, the pre-exposure of animals used to generate the EST libraries to stressors known to induce microevolutionary and adaptive responses in D. magna may have increased the efficiency in targeting loci under selection. Further studies on the genes where the outlier SNPs are located will be needed to identify candidate genes underlying adaptive responses to environmental stressors.
The analysis of the population structure based on the neutral loci in six natural populations of D. magna revealed patterns that are consistent with discrete populations distributed at regional scale. This is consistent with earlier work on population genetics of Daphnia populations. However, the level of genetic differentiation among our populations is lower than previously observed in D. magna[77, 78]. This can be explained by the use of different genetic markers. The use of allozymes [77, 78] vs. SNPs (present study) may have contributed to the somewhat larger values of genetic differentiation observed in previous studies compared to the current study. Even though we observe an overall lower genetic differentiation than previously observed, our results conform to the emerging pattern from literature, with significant genetic differentiation among populations even at relatively short geographic distances (e.g. [27, 29, 80, 81]). The use of SNPs in ecology and population genetics is still uncommon [2, 3], therefore it is difficult to draw comparisons between the range of genetic differentiation observed in our study and literature studies. However, recent simulation studies on the use of SNPs in ecology  illustrate that ~30 SNPs should be sufficient to detect moderate (FST = 0.01) levels of differentiation, while studies aimed at detecting demographic independence (e.g. FST < 0.005) may require 80 or more SNPs. The number of SNPs used in this study is dense for a species with a relatively small genome size (~200Mb). Additionally, the FST values measured in our populations range between 0.06 and 0.20. This conforms to high among-populations genetic differentiation previously observed using other genetic markers [77, 78].
We do not observe isolation-by-distance, which is in agreement with other studies at similar regional geographic scales for D. magna (e.g. ). The absence of isolation-by-distance reflects high genetic differentiation among neighboring populations rather than the absence of genetic differentiation at larger geographic distances. Earlier reports on isolation-by-distance in Daphnia refer to an increase in genetic differentiation at near-continental scales (phylogeographic patterns; see e.g.). Furthermore, the absence of isolation by distance coupled with a strong spatial structure of the populations conform to patterns encountered in other metapopulation systems and can be explained with high genetic drift in small populations .
Building upon the rapid development of genomic tools for D. magna we produced and validated a first set of SNP markers from EST sequences to use in future population genomic studies. By exposing animals that were used to produce the EST libraries to standardized selection pressures, we identified candidate EST potentially underlying responses to environmental stress. The SNPs developed in the present study represent an important first step toward the identification of candidate genes underpinning stress responses in natural populations.
We thank John Colbourne for helping with the genes annotation. We thank Anja Coors for providing the fish and parasite exposed samples and Tine Vandenbrouck for the supporting information on EST library development. Verification of carbaryl concentrations were conducted by Prof. W.P. de Voogt (University of Amsterdam). We thank Christoph R. Haag for providing the sequences of five nuclear genes that were used for SNP mining, and Alexander Triantafyllidis and Joost Raeymaekers for comments on an earlier version of the manuscript. MJ has a PhD fellowship of the IWT (Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen); ES was funded by the EU network of Excellence Marine Genomics Europe and the EU project Aquafirst (STREP-2004-513692); LO is a postdoctoral researcher for the FWO (Research Fund, Flanders). This project was financially supported by the GOA/08/06 project of the KULeuven Research Fund, and FWO project G.0229.09. The project is part of ESF EEFG EUROCORES project STRESSFLEA. Our work benefits from and contributes to the Daphnia Genomics Consortium.
- Avise JC: Phylogeography. The history and formation of species. 2000, Harvard University PressGoogle Scholar
- Morin PA, Luikart G, Wayne RK, group Sw: SNPs in ecology, evolution and conservation. TREE. 2004, 19: 208-216.Google Scholar
- Morin PA, Martien KK, Taylor BL: Assessing statistical power of SNPs for population structure and conservation studies. Mol Ecol Resources. 2009, 9: 66-73. 10.1111/j.1755-0998.2008.02392.x.View ArticleGoogle Scholar
- Fisher S, Lerman LS: Lenght-independent separation of DNA restriction fragments in two-dimensional gel electrophoresis. Cell. 1979, 16: 191-200. 10.1016/0092-8674(79)90200-9.View ArticleGoogle Scholar
- Goldrick MM: RNase cleavage-based methods for mutation/SNP detection, past and present. Hum Mutat. 2001, 18 (3): 190-204. 10.1002/humu.1175.View ArticlePubMedGoogle Scholar
- Orita M, Iwahana H, Kanazawa H, Hayashi K, Sekiya T: Detection of polymorphisms of human DNA by gel electrophoresis as single-strand conformation polymorphisms. PNAS. 1989, 86: 2766-2770. 10.1073/pnas.86.8.2766.View ArticlePubMedPubMed CentralGoogle Scholar
- Vijg J, van Orsouw NJ: Two-dimensional gene scanning: exploring human genetic variability. Electrophoresis. 1999, 20: 1239-1249. 10.1002/(SICI)1522-2683(19990101)20:6<1239::AID-ELPS1239>3.0.CO;2-T.View ArticlePubMedGoogle Scholar
- Xiao W, Oefner PJ: Denaturing high-performance liquid chromatography: A review. Hum Mutat. 2001, 17: 439-474. 10.1002/humu.1130.View ArticlePubMedGoogle Scholar
- Suh Y, Vijg J: SNP discovery in associating genetic variation with human disease phenotypes. Mutat Res. 2005, 573: 41-53.View ArticlePubMedGoogle Scholar
- Twyman RM: SNP discovery and typing technologies for pharmacogenomics. Curr Top Med Chem. 2004, 4: 1423-1431.View ArticlePubMedGoogle Scholar
- Coates BS, Sumerford DV, Hellmich RL, Lewis LC: Mining an Ostrinia nubilalis midgut expressed sequence tag (EST) library for candidate genes and single nucleotide polymorphisms (SNPs). Insect Mol Biol. 2008, 17: 607-620. 10.1111/j.1365-2583.2008.00833.x.View ArticlePubMedGoogle Scholar
- Colesa ND, Colemana CE, Christensena SA, Jellena EN, Stevensa MR, Bonifaciob A, Rojas-Beltranb JA, Fairbanksa DJ, Maughana PJ: Development and use of an expressed sequenced tag library in quinoa (Chenopodium quinoa Willd.) for the discovery of single nucleotide polymorphisms. Plant Science. 2005, 168: 439-447. 10.1016/j.plantsci.2004.09.007.View ArticleGoogle Scholar
- Colson I, Du Pasquier L, Ebert D: Intragenic tandem repeats in Daphnia magna: structure, function and distribution. BMC Res Notes. 2009, 2: 206-10.1186/1756-0500-2-206.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim DW, Jung TS, Nam SH, Kwon HR, Kim A, Chae SH, Choi SH, Kim DW, Kim RN, Park HS: GarlicESTdb: an online database and mining tool for garlic EST sequences. BMC Plant Biology. 2009, 9: 61-66. 10.1186/1471-2229-9-61.View ArticlePubMedPubMed CentralGoogle Scholar
- Luikart G, phillip R, Tallmon D, Jordan S, Taberlet P: The power and promise of population genomics: from genotyping to genome typing. Nat Rev Genet. 2003, 4: 981-994.View ArticlePubMedGoogle Scholar
- Wayne RK, Morin PA: Conservation genetics in the new molecular age. Front Ecol Environ. 2004, 2: 89-97. 10.1890/1540-9295(2004)002[0089:CGITNM]2.0.CO;2.View ArticleGoogle Scholar
- Ryman N, Palm S, Andre C, Carvalho GR, Dahlgren TG, Jorde PE, Laikre L, Larsson LC, Palme A, Ruzzante DE: Power for detecting genetic divergence: differences between statistical methods and marker loci. Mol Ecol. 2006, 15: 2031-2045. 10.1111/j.1365-294X.2006.02839.x.View ArticlePubMedGoogle Scholar
- Morin PA, McCarthy M: Highly accurate SNP genotyping from historical and low quality samples. Mol Ecol Resources. 2007, 7: 937-946.View ArticleGoogle Scholar
- Colbourne JK, Singan VRa, Gilbert D: wFleaBase: the Daphnia genome database. BMC Bioinformatics. 2005, 6: 45-10.1186/1471-2105-6-45.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaw JR, Colbourne JK, Davey JC, Glaholt SP, Hampton TH, Chen CY, Folt CL, Hamilton JW: Gene response profiles for Daphnia pulex exposed to the environmental stressor cadmium reveal a novel crustacean metallothionein. BMC genomics. 2007, 8: 477-10.1186/1471-2164-8-477.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaw JR, Pfrender M, Eads BD, Klaper R, Callaghan A, Colson I, Jansen B, Gilbert D, Colbourne JK: Daphnia as an emerging model for toxicological genomics. Advances in experimental biology on toxicogenomics. Edited by: Hogstrand C, Kille P. 2008, Elsevier Press, 165-219.Google Scholar
- Colbourne JK, Pfrender ME, Gilbert D, Thomas WK, Tucker A, Oakley TH, Tokishita S, Aerts A, Arnold GJ, Basu MK, et al: The ecoresponsive genome of Daphnia pulex. Science. 2011, 331: 555-561. 10.1126/science.1197761.View ArticlePubMedPubMed CentralGoogle Scholar
- Coors A, Vanoverbeke J, De Bie T, De Meester L: Land use, genetic diversity and toxicant tolerance in natural populations of Daphnia magna. Aquat Toxicol. 2009, 95: 71-79. 10.1016/j.aquatox.2009.08.004.View ArticlePubMedGoogle Scholar
- Cousyn C, De Meester L, Colbourne JK, Brendonck L, Verchuren D, Volckaert F: Rapid, local adaptation of zooplankton behavior to changes in predation pressure in the absence of neutral genetic changes. PNAS. 2001, 98: 6256-6260. 10.1073/pnas.111606798.View ArticlePubMedPubMed CentralGoogle Scholar
- Decaestecker E, Gaba S, Raeymaekers J, Stoks R, Van Kerckhoven L, Ebert D, De Meester L: Host-parasite Red Queen dynamics archived in pond sediment. Nature. 2007, 450: 870-874. 10.1038/nature06291.View ArticlePubMedGoogle Scholar
- Carvalho GR: Genetics of aquatic clonal organisms. Genetic and evolution of aquatic organisms. Edited by: Beaumont AR. 1994, London: Chapman and Hall, 291-323.Google Scholar
- De Meester L, Vanoverbeke J, De Gelas K, Ortells R, Spaak P: Genetic structure of cyclic parthenogenetic zooplankton populations- a conceptual framework. Arch Hydrobiol. 2006, 167: 217-244. 10.1127/0003-9136/2006/0167-0217.View ArticleGoogle Scholar
- Haag CR, D. E: Genotypic selection in Daphnia populations consisting of inbred sibships. J Evol Biol. 2007, 20: 881-891. 10.1111/j.1420-9101.2007.01313.x.View ArticlePubMedGoogle Scholar
- Hebert PDN: Genetics of Daphnia. The Biology of Daphnia. Edited by: Peters R, de Bernardi R. 1987, Occas. Pub. Ital. Inst. Hydrobiol, 45: 439-460.Google Scholar
- Routtu J, Jansen B, Colson I, De Meester L, Ebert D: The first-generation Daphnia magna linkage map. BMC genomics. 11: 508-
- Pauwels K, Stoks R, De Meester L: Enhanced anti-predator defence in the presence of food stress in the water flea Daphnia magna. Funct Ecol. 2010, 322-329. 24
- Ben-Ami F, Regoes RR, Ebert D: A quantitative test of the relationship between parasite dose and infection probability across different host-parasite combinations. Proc Biol Sci. 2008, 275: 853-859. 10.1098/rspb.2007.1544.View ArticlePubMedPubMed CentralGoogle Scholar
- Klüttgen BUD, Engels MTRH: Combined effects of 3,4-dichloroaniune and food concentration on life-table data of two related cladocerans, Daphnia magna and Ceriodaphnia quadrangula. Chemosphere. 1994, 32: 2015-2028.View ArticleGoogle Scholar
- Jansen M, Stoks R, Decaestecker E, Coors A, Van De Meutter V, De Meester L: Local exposure shapes spatial patterns in infectivity and community structure of Daphnia parasites. J Anim Ecol. 2010, 79: 1023-1033. 10.1111/j.1365-2656.2010.01718.x.View ArticlePubMedGoogle Scholar
- Coors A, De Meester L: Synergistic, antagonistic and additive effects of multiple stressors: predation threat, parasitism and pesticide exposure in Daphnia magna. J Appl Ecol. 2008, 45: 1820-1828. 10.1111/j.1365-2664.2008.01566.x.View ArticleGoogle Scholar
- Coors A, Decaestecker E, Jansen M, De Meester L: Pesticide exposure strongly enhances parasite virulence in an invertebrate host model. Oikos. 2008, 117: 1840-1846. 10.1111/j.1600-0706.2008.17028.x.View ArticleGoogle Scholar
- Sambrook J, Russell D: Molecular Cloning: A Laboratory Manual. 2000, Cold Spring Harbor Laboratory PressGoogle Scholar
- Soetaert A, Moens LN, Van der Ven K, Van Leemput K, Naudts B, Blust R, De Coen WM: Molecular impact of propiconazole on Daphnia magna using a reproduction-related cDNA array. Comp Biochem Physiol C: Toxicol Pharmacol. 2006, 142: 66-76. 10.1016/j.cbpc.2005.10.009.Google Scholar
- Moens LN, Van der Ven K, Caturla M, De Coen WM: Toxicogenomics applications in environmental toxicology. Vitro Methods in Aquatic Toxicology. 2003, Chichester, UK Springer-Praxis Publishing, 185-220.Google Scholar
- Ewing B, Hillier L, Wendl M, Green P: Basecalling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.View ArticlePubMedGoogle Scholar
- Watanabe H, Tatarazako N, Oda S, Nishide H, Uchiyama I, Morita M, Iguchi T: Analysis of expressed sequence tags of the water flea Daphnia magna. Genome. 2005, 48 (4): 606-609. 10.1139/g05-038.View ArticlePubMedGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al: TIGR Gene Indices clustering tools (TGILC): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.View ArticlePubMedPubMed CentralGoogle Scholar
- Min XJ, Butler G, Storms R, Tsang A: OrfPredictor: predicting protein-coding regions in EST-derived sequences. Nucleic Acids Res. 2005, 33: 677-680. 10.1093/nar/gki394.View ArticleGoogle Scholar
- Conesa A, Gotz S: Blast2GO: A Comprehensive Suite for Functional Analysis in Plant Genomics. Int J Plant Genomics. 2008, 2008: 619832-View ArticlePubMedGoogle Scholar
- De Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N: ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006, 34: 362-365. 10.1093/nar/gkl124.View ArticleGoogle Scholar
- Chen F, Mackey AJ, Stoeckert CJ, Roos DS: OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res. 2006, 34: D363-368. 10.1093/nar/gkj123.View ArticlePubMedGoogle Scholar
- Savage D, Batley J, Erwin T, Logan E, Love CG, Lim GA, Mongin E, Barker G, Spangenberg GC, Edwards D: SNPServer: a real-time SNP discovery tool. Nucleic Acids Res. 2005, 33: W493-495. 10.1093/nar/gki462.View ArticlePubMedPubMed CentralGoogle Scholar
- Marth GT, Korf I, Yandell MD, Yeh RT, Gu Z, Zakeri H, Stitziel NO, Hillier L, Kwok PY, Gish WR: A general approach to single-nucleotide polymorphism discovery. Nat Genet. 1999, 23: 452-456. 10.1038/70570.View ArticlePubMedGoogle Scholar
- Tang J, Vosman B, Voorrips RE, van der Linden CG, Leunissen JA: QualitySNP: a pipeline for detecting single nucleotide polymorphisms and insertions/deletions in EST data from diploid and polyploid species. BMC Bioinformatics. 2006, 7: 438-10.1186/1471-2105-7-438.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Huang X: A method for finding single-nucleotide polymorphisms with allele frequencies in sequences of deep coverage. BMC Bioinformatics. 2005, 6: 220-10.1186/1471-2105-6-220.View ArticlePubMedPubMed CentralGoogle Scholar
- Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WE, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004, 14: 1147-1159. 10.1101/gr.1917404.View ArticlePubMedPubMed CentralGoogle Scholar
- Ning Z, Cox AJ, Mullikin JC: SSAHA: a fast search method for large DNA databases. Genome Res. 2001, 11: 1725-1729. 10.1101/gr.194201.View ArticlePubMedPubMed CentralGoogle Scholar
- Haag CR, McTaggart SJ, Didier A, Little TJ, Charlesworth D: Nucleotide polymorphism and within-gene recombination in Daphnia magna and D. pulex, two cyclical parthenogens. Genetics. 2009, 182: 313-323. 10.1534/genetics.109.101147.View ArticlePubMedPubMed CentralGoogle Scholar
- De Meester L, Gomez A, Okamura B, Schwenk K: The Monopolization Hypothesis and the dispersal-gene flow paradox in aquatic organisms. Acta Oecol. 2002, 23: 121-135. 10.1016/S1146-609X(02)01145-1.View ArticleGoogle Scholar
- Vanoverbeke J, De Meester L: Clonal erosion and genetic drift in cyclical parthenogens - the interplay between neutral and selective processes. J Evol Biol. 2010, 23: 997-1012. 10.1111/j.1420-9101.2010.01970.x.View ArticlePubMedGoogle Scholar
- Antao T, Lopes A, Lopes RJ, Beja-Pereira A, Luikart G: LOSITAN: a workbench to detect molecular adaptation based on a Fst-outlier method. BMC Bioinformatics. 2008, 9: 323-10.1186/1471-2105-9-323.View ArticlePubMedPubMed CentralGoogle Scholar
- Beaumont MA, Nichols RA: Evaluating loci for use in the genetic analysis of population structure. P Roy Soc B. 1996, 363: 1619-1626.View ArticleGoogle Scholar
- Foll M, Gaggiotti O: A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics. 2008, 180: 977-993. 10.1534/genetics.108.092221.View ArticlePubMedPubMed CentralGoogle Scholar
- Dieringer D, Schlötterer C: Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes. 2003, 3: 167-169. 10.1046/j.1471-8286.2003.00351.x.View ArticleGoogle Scholar
- Corander J, Marttinen P: Bayesian identification of admixture events using multi-locus molecular markers. Mol Ecol. 2006, 15: 2833-2843. 10.1111/j.1365-294X.2006.02994.x.View ArticlePubMedGoogle Scholar
- Corander J, Sirén J, Arjas E: Bayesian spatial modeling of genetic population structure. Computational Stat. 2008, 23: 111-129. 10.1007/s00180-007-0072-x.View ArticleGoogle Scholar
- Guillot G, Estoup A, Mortier F, Cosson JF: A spatial statistical model for landscape genetics. Genetics. 2005, 170: 1261-1280. 10.1534/genetics.104.033803.View ArticlePubMedPubMed CentralGoogle Scholar
- Bohonak AJ: IBD (Isolation By Distance): A program for analyses of isolation by distance. J Hered. 2002, 93: 153-154. 10.1093/jhered/93.2.153.View ArticlePubMedGoogle Scholar
- Hayes B, Laerdahl JK, Lien S, Moen T, Berg P, Hindar K, Davidson WS, Koop BF, Adzhubei A, Hoyheim B: An extensive resource of single nucleotide polymorphism markers associated with Atlantic salmon (Salmo salar) expressed sequences. Aquaculture. 2007, 265: 82-90. 10.1016/j.aquaculture.2007.01.037.View ArticleGoogle Scholar
- Akey JM, Zhang K, Xiong M, Jin L: The effect of single nucleotide polymorphism identification strategies on estimates of linkage disequilibrium. Mol Biol Evol. 2003, 20: 232-242. 10.1093/molbev/msg032.View ArticlePubMedGoogle Scholar
- Picoult-Newberg L, Ideker TE, Pohl MG, Taylor SL, Donaldson MA, Nickerson DA, Boyce-Jacino M: Mining SNPs from EST databases. Genome Res. 1999, 9: 167-174.PubMedPubMed CentralGoogle Scholar
- Wang S, Sha Z, Sonstegard TS, Liu H, Xu P, Somridhivej B, Peatman E, Kucuktas H, Liu Z: Quality assessment parameters for EST-derived SNPs from catfish. BMC genomics. 2008, 9: 450-10.1186/1471-2164-9-450.View ArticlePubMedPubMed CentralGoogle Scholar
- Jalving R, van't Slot R, van Oost BA: Chicken single nucleotide polymorphism identification and selection for genetic mapping. Poult Sci. 2004, 83: 1925-1931.View ArticlePubMedGoogle Scholar
- Akey JM, Zhang G, Zhang K, Jin L, Shriver MD: Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 2002, 12: 1805-1814. 10.1101/gr.631202.View ArticlePubMedPubMed CentralGoogle Scholar
- Campbell D, Bernatchez L: Generic scan using AFLP markers as a means to assess the role of directional selection in the divergence of sympatric whitefish ecotypes. Mol Biol Evol. 2004, 21: 945-956. 10.1093/molbev/msh101.View ArticlePubMedGoogle Scholar
- Achere V, Favre JM, Besnard G, Jeandroz S: Genomic organization of molecular differentiation in Norway spruce (Picea abies). Mol Ecol. 2005, 14: 3191-3201. 10.1111/j.1365-294X.2005.02646.x.View ArticlePubMedGoogle Scholar
- Namroud MC, Beaulieu J, Juge N, Laroche J, Bousquet J: Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce. Mol Ecol. 2008, 17: 3599-3613. 10.1111/j.1365-294X.2008.03840.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Vasemagi A, Nilsson J, Primmer CR: Expressed sequence tag-linked microsatellites as a source of gene-associated polymorphisms for detecting signatures of divergent selection in atlantic salmon (Salmo salar L.). Mol Biol Evol. 2005, 22: 1067-1076. 10.1093/molbev/msi093.View ArticlePubMedGoogle Scholar
- Bonin A, Taberlet P, Miaud C, Pompanon F: Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria). Mol Biol Evol. 2006, 23: 773-783. 10.1093/molbev/msj087.View ArticlePubMedGoogle Scholar
- Makinen HS, Cano JM, Merila J: Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations. Mol Ecol. 2008, 17: 3565-3582. 10.1111/j.1365-294X.2008.03714.x.View ArticlePubMedGoogle Scholar
- Haag CR, Riek M, Hottinger JW, Pajunen VI, Ebert D: Founder events as determinants of within-island and among-island genetic structure of Daphnia metapopulations. Heredity. 2006, 96: 150-158. 10.1038/sj.hdy.6800774.View ArticlePubMedGoogle Scholar
- Vanoverbeke J, De Meester L: Among-population genetic differentiation in the cyclical parthenogen Daphnia magna (Crustacea, Anomopoda) and its relation to geographic distance and conal diversity. Hydrobiologia. 1997, 360: 135-142. 10.1023/A:1003160903708.View ArticleGoogle Scholar
- Hedrick PW: Perspective: highly variable loci and their interpretation in evolution and conservation. Evolution. 1999, 53: 313-318. 10.2307/2640768.View ArticleGoogle Scholar
- Amend A, Garbelotto M, Fang Z, Keeley S: Isolation by landscape in populations of a prized edible mushroom Tricholoma matsutake. Conserv Genet. 2010, 11: 795-802. 10.1007/s10592-009-9894-0.View ArticleGoogle Scholar
- Orsini L, Corander J, Alasentie A, Hanski I: Genetic spatial structure in a butterfly metapopulation correlates better with past than present demographic structure. Mol Ecol. 2008, 17: 2629-2642. 10.1111/j.1365-294X.2008.03782.x.View ArticlePubMedGoogle Scholar
- De Gelas K, De Meester L: Phylogeography of Daphnia magna in Europe. Mol Ecol. 2005, 14: 753-764. 10.1111/j.1365-294X.2004.02434.x.View ArticlePubMedGoogle Scholar