The forgotten variable? Does the euthanasia method and sample storage condition influence an organisms transcriptome - a gene expression analysis on multiple tissues in pigs.

BACKGROUND: Transcriptomic studies often require collection of fresh tissues post euthanasia. The chosen euthanasia method might have the potential to induce variations in gene expressions that are unlinked with the experimental design. The present study compared the suitability of 'nitrogen gas in foam' (ANOXIA) in comparison to a non-barbiturate anaesthetic, T-61® (T61), for euthanizing piglets used in transcriptome research. Further, the effect of common tissue storage conditions, RNAlater™ (RL) and snap freezing in liquid nitrogen (LN2), on gene expression profiles were also analysed. RESULTS: On comparison of the 3'mRNA-Seq data generated from pituitary, hypothalamus, liver and lung tissues, no significant differential expression in the protein coding genes were detected between the euthanasia methods. This implies that the nitrogen anoxia method could be a suitable alternative for euthanasia of piglets used in transcriptomic research. However, small nuclear RNAs (snRNAs) that constitute the eukaryotic spliceosomal machinery were found to be significantly higher (log2fold change ≥ 2.0, and adjusted p value ≤ 0.1) in pituitary samples collected using ANOXIA. Non-protein coding genes like snRNAs that play an important role in pre-mRNA splicing can subsequently modify gene expression. Storage in RL was found to be superior in preserving RNA compared to LN2 storage, as evidenced by the significantly higher RIN values in representative samples. However, storage in RL as opposed to LN2, also influenced differential gene expression in multiple tissues, perhaps as a result of its inability to inhibit biological activity during storage. Hence such external sources of variations should be carefully considered before arriving at research conclusions. CONCLUSIONS: Source of biological variations like euthanasia method and storage condition can confound research findings. Even if we are unable to prevent the effect of these external factors, it will be useful to identify the impact of these variables on the parameter under observation and thereby prevent misinterpretation of our results.

Differential gene expression in narrow- and broad-headed European glass eels (Anguilla anguilla) points to a transcriptomic link of head shape dimorphism with growth rate and chemotaxis.

One of the major challenges in evolutionary biology is to understand the mechanisms underlying morphological dimorphism and plasticity, including the genomic basis of traits and links to ecology. At the yellow eel stage of the European eel (Anguilla anguilla), two morphotypes are found: broad- and narrow-heads. This dimorphism has been linked to dietary differences, with broad-heads feeding on harder, larger prey than narrow-heads. However, recent research showed that both morphotypes could be distinguished at the glass eel stage, the nonfeeding predecessor of the yellow eel stage, implying that nondietary factors play a role in the development of this head shape dimorphism. Here, we used transcriptome profiling (RNAseq) to identify differentially expressed genes between broad- and narrow-headed glass eels. We found 260 significantly differentially expressed genes between the morphotypes, of which most were related to defence and immune responses. Interestingly, two genes involved in growth (soma and igf2) were significantly upregulated in narrow-heads, while nine genes involved in chemotaxis showed significant differential expression. Thus, we found support for the observation that head shape is associated with somatic growth, with fast-growing eels developing a narrower head. Additionally, observations in the wild have shown that slow-growers prefer freshwater, while fast-growers prefer brackish water. The differential expression of genes involved in chemotaxis seems to indicate that glass eel growth rate and habitat choice are linked. We hypothesize that two levels of segregation could take place in the European eel: first according to habitat choice and second according to feeding preference.

Genome-wide single-generation signatures of local selection in the panmictic European eel.

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus).

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene-linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high-resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms.

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Temporal genetic stability and high effective population size despite fisheries-induced life-history trait evolution in the North Sea sole.

Heavy fishing and other anthropogenic influences can have profound impact on a species' resilience to harvesting. Besides the decrease in the census and effective population size, strong declines in mature adults and recruiting individuals may lead to almost irreversible genetic changes in life-history traits. Here, we investigated the evolution of genetic diversity and effective population size in the heavily exploited sole (Solea solea), through the analysis of historical DNA from a collection of 1379 sole otoliths dating back from 1957. Despite documented shifts in life-history traits, neutral genetic diversity inferred from 11 microsatellite markers showed a remarkable stability over a period of 50 years of heavy fishing. Using simulations and corrections for fisheries induced demographic variation, both single-sample estimates and temporal estimates of effective population size (N(e) ) were always higher than 1000, suggesting that despite the severe census size decrease over a 50-year period of harvesting, genetic drift is probably not strong enough to significantly decrease the neutral diversity of this species in the North Sea. However, the inferred ratio of effective population size to the census size (N(e) /N(c) ) appears very small (10(-5) ), suggesting that overall only a low proportion of adults contribute to the next generation. The high N(e) level together with the low N(e) /N(c) ratio is probably caused by a combination of an equalized reproductive output of younger cohorts, a decrease in generation time and a large variance in reproductive success typical for marine species. Because strong evolutionary changes in age and size at first maturation have been observed for sole, changes in adaptive genetic variation should be further monitored to detect the evolutionary consequences of human-induced selection.

Conservation of the introgressed European water frog complex using molecular tools.

In Belgium, the Pelophylax esculentus complex has recently been subjected to multiple introductions of non-native water frogs, increasing the occurrence of hybridisation events. In the present study, we tested the reliability of morphometric and recently developed microsatellite tools to identify introgression and to determine the origin of exotic Belgian water frogs. By analysing 150 individuals of each taxon of the P. esculentus complex and an additional 60 specimens of the introduced P. cf. bedriagae, we show that neither of the currently available tools appears to have sufficient power to reliably distinguish all Belgian water frog species. We therefore aimed at increasing the discriminatory power of a microsatellite identification tool by developing a new marker panel with additional microsatellite loci. By adding only two new microsatellite loci (RlCA5 and RlCA1b20), all taxa of the P. esculentus complex could be distinguished from each other with high confidence. Three more loci (Res3, Res5 and Res17) provided a powerful discrimination of the exotic species.

Morphological and genetic seasonal dynamics of European eel Anguilla anguilla recruitment in southern France.

The fine scale morphological and genetic dynamics of successive waves of Anguilla anguilla glass eel recruitment was studied over a 2 year period at a southern European Mediterranean location (Camargue, France) with continuous recruitment. Using morphometric [total length (L(T)), mass (M), condition (K) and pigmentation stage] as well as genetic (allozyme) markers, the aim was to test for the existence of temporally separated spawning groups and explore the relation between genetic variability and morphological heterogeneity of recruits. The results showed that L(T), M and K varied over time, being highest from the end of summer to winter (peaking in December) and lowest in spring (lowest in April). The pigmentation stages within monthly samples were highly diverse with a heterogeneous seasonal pattern. Allozyme data showed high genetic variability values within samples, but low genetic differentiation among samples (F(ST) = 0.003, P < 0.05). Pairwise comparisons between samples indicated a positive correlation between genetic differentiation and difference in recruitment time (days), with a marked increase in genetic differentiation around 250 days between monthly recruitment samples. Furthermore, genetic diversity increased with the number of pigmentation stages per sample and was negatively correlated with the North Atlantic Oscillation (NAO) index during the putative year of trans-oceanic migration. No correlation, however, was found between the level of multilocus heterozygosity (MLH) and growth variables. A situation of genetic patchiness with fluctuating parental contribution can thus best explain the patterns observed, although the existence of two separate spawning periods cannot be excluded. More discriminatory and sensitive genetic markers, such as (neutral and adaptive) microsatellites, could probably provide additional insights into the most probable hypothesis explaining the population structure and recruitment heterogeneity of A. anguilla.

Response of mixed methanotrophic consortia to different methane to oxygen ratios.

Methane (CH4) and oxygen (air) concentrations affect the CH4 oxidation capacity (MOC) and mixed methanotrophic community structures in compost (fresh) and landfill (age old) top cover soils. A change in the mixed methanotrophic community structure in response has implications for landfill CH4 bio-filter remediation and possible bio-product outcomes (i.e., fatty acid methyl esters (FAME) content and profiles and polyhydroxybutyrate (PHB) contents). Therefore the study aimed to evaluate the effect of variable CH4 to oxygen ratios (10-50% CH4 in air) on mixed methanotrophic community structures enriched from landfill top cover (LB) and compost soils (CB) and to quantify flow on impacts on MOC, total FAME contents and profiles, and PHB accumulation. A stable consortium developed achieving average MOCs of 3.0±0.12, 4.1±0.26, 6.9±0.7, 7.6±1.3 and 9.2±1.2mgCH4g-1DWbiomassh-1 in LB and 2.9±0.04, 5.05±0.32, 6.7±0.31, 7.9±0.61 and 8.6±0.48mgCH4g-1DWbiomassh-1 in CB for a 20day cultivation period at 10, 20, 30, 40 and 50% CH4, respectively. CB at 10% CH4 had a maximal FAME content of 40.5±0.8mgFAMEg-1DWbiomass, while maximal PHB contents (25mgg-1DWbiomass) was observed at 40% CH4 in LB. Despite variable CH4/O2 ratios, the mixed methanotrophic community structures in both LB and CB were relatively stable, dominated by Methylosarcina sp., and Chryseobacterium, suggesting that a resilient consortium had formed which can now be tested in bio-filter operations for CH4 mitigations in landfills.

Growth rate correlates to individual heterozygosity in the European eel, Anguilla anguilla L.

Heterozygosity-fitness correlations (HFCs) have been reported in populations of many species. We provide evidence for a positive correlation between genetic variability and growth rate at 12 allozyme loci in a catadromous marine fish species, the European eel (Anguilla anguilla L.). More heterozygous individuals show a significantly higher length and weight increase and an above average condition index in comparison with more homozygous individuals. To a lesser extent, six microsatellite loci show a similar pattern, with positive but not significant correlations between heterozygosity and growth rate. The HFCs observed could be explained by an effect of either direct allozyme over-dominance or associative overdominance. Selection affecting some of the allozyme loci would explain the greater strength of the HFCs found at allozymes in comparison with microsatellites and the lack of correlation between MLH at allozymes and MLH at microsatellites. Associative overdominance (where allozyme loci are merely acting as neutral markers of closely linked fitness loci) might provide an explanation for the HFCs if we consider that allozyme loci have a higher chance than microsatellites to be in linkage disequilibrium with fitness loci.

The catadromous European eel Anguilla anguilla (L.) as a model for freshwater evolutionary ecotoxicology: relationship between heavy metal bioaccumulation, condition and genetic variability.

Understanding the effects of pollutants on the genome is of crucial importance to preserve the evolutionary potential of endangered natural populations. The highly vagile European eel (Anguilla anguilla L.) has suffered a dramatic decline in recruitment since two decades, urging for a better understanding of the genetic impact of pollution. Its catadromous life history constitutes a model to assess local selection of pollutants on condition and genetic variability, as juveniles recruit in European rivers without appreciable pollution load or interfering genetic background. Because of its high fat content and local benthic feeding behaviour, the feeding stage is considered extremely prone to the bioaccumulation of pollutants. We studied the relationship between heavy metal bioaccumulation, fitness (condition) and genetic variability in the European eel. The muscle tissues of 78 sub-adult eels, originating from three Belgian river basins (Scheldt, Meuse and Yser), were examined for nine heavy metal pollutants (Hg, Cd, Pb, Cu, Zn, Ni, Cr, As and Se), while in total 123 individuals were genotyped at 12 allozyme and 8 microsatellite loci. A significant negative correlation between heavy metal pollution load and condition was observed, suggesting an impact of pollution on the health of sub-adult eels. In general, we observed a reduced genetic variability in strongly polluted eels, as well as a negative correlation between level of bioaccumulation and allozymatic multi-locus heterozygosity (MLH). Microsatellite genetic variability did not show any pollution related differences, suggesting a differential response at metabolic enzymes and possibly direct overdominance of heterozygous individuals.

Gene transcription reflects poor health status of resident European eel chronically exposed to environmental pollutants.

Understanding the effects of chronic exposure to pollutants on the genome and transcriptome of diadromous fish populations is crucial for their resilience under combined anthropogenic and environmental selective pressures. The catadromous European eel (Anguilla anguilla L.) has suffered a dramatic decline in recruitment for three decades, necessitating a thorough assessment of the transcriptional effects of environmental pollutants on resident and migrating eels in natural systems. We investigated the relationship between muscular bioaccumulation levels of metals (Hg, Cd, Pb, Cu, Zn, Ni, Cr, As and Se), PCBs and organochlorine pesticides (DDTs), the health status (condition factor and lipid reserves) and the associated transcriptional response in liver and gill tissues for genes involved in metal detoxification (metallothionein, MT) and oxidative metabolism (cytochrome P4501A, CYP1A) of xenobiotic compounds. In total 84 resident eels originating from three Belgian river basins (Scheldt, Meuse and Yzer) were analyzed along with five unpolluted aquaculture samples as control group. There was a large spatial variation in individual contaminant intensity and profile, while tissue pollution levels were strongly and negatively associated with condition indices, suggesting an important impact of pollution on the health of sub-adult resident eels. Gene transcription patterns revealed a complex response mechanism to a cocktail of pollutants, with a high variation at low pollution levels, but strongly down-regulated hepatic and gill gene transcription in highly polluted eels. Resident eels clearly experience a high pollution burden and seem to show a dysfunctional gene transcription regulation of detoxification genes at higher pollutant levels, correlated with low energy reserves and condition. To fully understand the evolutionary implications of pollutants on eel reproductive fitness, analyses of mature migrating eels and the characterization of their transcriptome-wide gene transcription response would be appropriate to unveil the complex responses associated with multiple interacting stressors and the long-term consequences at the entire species level. In the meanwhile, jointly monitoring environmental and tissue pollution levels at a European scale should be initiated, while preserving high quality habitats to increase the recovery chance of European eel in the future.

A resource of genome-wide single-nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel.

Reduced representation genome sequencing such as restriction-site-associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single-nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome-wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome-wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19,703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long-term effective population size was estimated to range between 132,000 and 1,320,000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82,425 loci and 376,918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.

Detecting genome-wide gene transcription profiles associated with high pollution burden in the critically endangered European eel.

The European eel illustrates an example of a critically endangered fish species strongly affected by human stressors throughout its life cycle, in which pollution is considered to be one of the factors responsible for the decline of the stock. The objective of our study was to better understand the transcriptional response of European eels chronically exposed to pollutants in their natural environment. A total of 42 pre-migrating (silver) female eels from lowly, highly and extremely polluted environments in Belgium and, for comparative purposes, a lowly polluted habitat in Italy were measured for polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs) and brominated flame retardants (BFRs). Multipollutant level of bioaccumulation was linked to their genome-wide gene transcription using an eel-specific array of 14,913 annotated cDNAs. Shared responses to pollutant exposure were observed when comparing the highly polluted site in Belgium with the relatively clean sites in Belgium and Italy. First, an altered pattern of transcription of genes was associated with detoxification, with a novel European eel CYP3A gene and gluthatione S-transferase transcriptionally up-regulated. Second, an altered pattern of transcription of genes associated with the oxidative phosphorylation pathway, with the following genes involved in the generation of ATP being transcriptionally down-regulated in individuals from the highly polluted site: NADH dehydrogenase, succinate dehydrogenase, ubiquinol-cytochrome c reductase, cytochrome c oxidase and ATP synthase. Although we did not measure metabolism directly, seeing that the transcription level of many genes encoding enzymes involved in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated in the highly polluted site suggests that pollutants may have a significant effect on energy metabolism in these fish.

Permanent genetic resources added to Molecular Ecology Resources Database 1 August 2011-30 September 2011.

This article documents the addition of 299 microsatellite marker loci and nine pairs of single-nucleotide polymorphism (SNP) EPIC primers to the Molecular Ecology Resources (MER) Database. Loci were developed for the following species: Alosa pseudoharengus, Alosa aestivalis, Aphis spiraecola, Argopecten purpuratus, Coreoleuciscus splendidus, Garra gotyla, Hippodamia convergens, Linnaea borealis, Menippe mercenaria, Menippe adina, Parus major, Pinus densiflora, Portunus trituberculatus, Procontarinia mangiferae, Rhynchophorus ferrugineus, Schizothorax richardsonii, Scophthalmus rhombus, Tetraponera aethiops, Thaumetopoea pityocampa, Tuta absoluta and Ugni molinae. These loci were cross-tested on the following species: Barilius bendelisis, Chiromantes haematocheir, Eriocheir sinensis, Eucalyptus camaldulensis, Eucalyptus cladocalix, Eucalyptus globulus, Garra litaninsis vishwanath, Garra para lissorhynchus, Guindilla trinervis, Hemigrapsus sanguineus, Luma chequen. Guayaba, Myrceugenia colchagüensis, Myrceugenia correifolia, Myrceugenia exsucca, Parasesarma plicatum, Parus major, Portunus pelagicus, Psidium guayaba, Schizothorax richardsonii, Scophthalmus maximus, Tetraponera latifrons, Thaumetopoea bonjeani, Thaumetopoea ispartensis, Thaumetopoea libanotica, Thaumetopoea pinivora, Thaumetopoea pityocampa ena clade, Thaumetopoea solitaria, Thaumetopoea wilkinsoni and Tor putitora. This article also documents the addition of nine EPIC primer pairs for Euphaea decorata, Euphaea formosa, Euphaea ornata and Euphaea yayeyamana.

Application of SNPs for population genetics of nonmodel organisms: new opportunities and challenges.

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 December 2009-31 January 2010.

This article documents the addition of 220 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Allanblackia floribunda, Amblyraja radiata, Bactrocera cucurbitae, Brachycaudus helichrysi, Calopogonium mucunoides, Dissodactylus primitivus, Elodea canadensis, Ephydatia fluviatilis, Galapaganus howdenae howdenae, Hoplostethus atlanticus, Ischnura elegans, Larimichthys polyactis, Opheodrys vernalis, Pelteobagrus fulvidraco, Phragmidium violaceum, Pistacia vera, and Thunnus thynnus. These loci were cross-tested on the following species: Allanblackia gabonensis, Allanblackia stanerana, Neoceratitis cyanescens, Dacus ciliatus, Dacus demmerezi, Bactrocera zonata, Ceratitis capitata, Ceratitis rosa, Ceratits catoirii, Dacus punctatifrons, Ephydatia mülleri, Spongilla lacustris, Geodia cydonium, Axinella sp., Ischnura graellsii, Ischnura ramburii, Ischnura pumilio, Pistacia integerrima and Pistacia terebinthus.

Isolation and characterization of expressed sequence tag-linked microsatellite loci for the European eel (Anguilla anguilla).

A European eel (Anguilla anguilla) expressed sequence tag database consisting of 795 contigs and 4008 singletons was screened for microsatellites sequences. Primers were designed to amplify 96 repeats, of which 86 gave good quality amplification products. Twenty-eight microsatellites were selected for further microsatellite genotyping. Only two loci were found to be monomorphic; out of the 26 polymorphic loci, number of alleles per locus ranged from two to 14, while the observed and expected heterozygosities ranged from 0.05 to 0.93, and from 0.05 to 0.95, respectively. All 28 primer sets tested revealed positive amplification in American eel (Anguilla rostrata).

Influence of DNA isolation from historical otoliths on nuclear-mitochondrial marker amplification and age determination in an overexploited fish, the common sole (Solea solea L.).

Historical otolith collections are crucial in assessing the evolutionary consequences of natural and anthropogenic changes on the demography and connectivity of commercially important fish species. Hence, it is important to define optimal protocols for purifying DNA from such valuable information sources while avoiding any damage to the physical structure of the otolith. Before being able to conclude on the harmlessness of a method, it is important to validate protocols on different kinds of otoliths by testing purification methodologies under standardized conditions. Here we compare the effect of two DNA extraction methods on the success in identifying the age in an overexploited marine fish, the common sole (Solea solea L.). To ensure optimal future population genetic and demographic analyses, we assessed DNA quantity and tested the DNA quality by investigating the amplification success of a mitochondrial and nuclear marker. Our results show that the choice of the DNA extraction method had a significant effect on the success of using these otoliths in age and growth analyses. Standard commercial and published protocols resulted in a severe damaging of the otolith structure, hampering accurate preparation and analyses of the morphological structures of the otoliths. Shortening the lysis time and lowering the EDTA (ethylene diamine tetraacetic acid) and SDS (sodium dodecylsulphate) concentration turned out to be beneficial for the stability of otolith structure, while maintaining an overall high DNA quality measured through polymerase chain reaction amplification success. We therefore recommend that care should be taken when choosing the extraction method for a molecular study on archived samples, in order to enable the maximal use of information embedded in historical material.

High genetic diversity in cryptic populations of the migratory sutchi catfish Pangasianodon hypophthalmus in the Mekong River.

The detection and conservation of spawning units is of crucial importance in highly migratory species. The sutchi catfish Pangasianodon hypophthalmus (Pangasiidae; Teleostei) is a common large-sized tropical fish, which migrates annually to several upstream spawning sites on the Lower Mekong River and feeds on the huge floodplain of the Lower Mekong and Tonle Sap for the other half of the year. We hypothesised that because of the relative size of the feeding and spawning habitat, genetic variability would be high and homogeneous in foraging populations, but that spawning stocks would be distinct in space and time. To test these predictions, 567 individuals from 10 geographic locations separated by up to 1230 km along the Lower Mekong River were genotyped at seven microsatellite loci. The level of genetic diversity was much higher than other freshwater fish and reached values comparable to marine species (mean H(e)=0.757). All samples collected at the potential spawning sites deviated from Hardy-Weinberg expectations, suggesting admixture. Individual-based clustering methods revealed genetic heterogeneity and enabled the detection of three genetically distinct sympatric populations. There was no evidence of recent reduction in effective population size in any population. Contrasting with the vast extent of the feeding grounds, the shortage of spawning grounds seems to have moved sutchi catfish towards diachronous spawning. Hence the sustainable exploitation of this natural resource hinges on the conservation of the limited spawning grounds and open migration routes between the spawning and feeding grounds.