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.
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ᵒN, we examined the patterns of genome‐wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using FST‐based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (FST = 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.
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 (LN
), on gene expression profiles were also analysed.
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 LN
storage, as evidenced by the significantly higher RIN values in representative samples. However, storage in RL as opposed to LN
, 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.
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.
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.
Chronic obstructive pulmonary disease (COPD) is a highly prevalent respiratory disease characterized by airflow limitation and chronic inflammation. MiR-155 is described as an ancient regulator of ...the immune system. Our objective was to establish a role for miR-155 in cigarette smoke (CS)-induced inflammation and COPD. We demonstrate increased miR-155 expression by RT-qPCR in lung tissue of smokers without airflow limitation and patients with COPD compared to never smokers and in lung tissue and alveolar macrophages of CS-exposed mice compared to air-exposed mice. In addition, we exposed wild type and miR-155 deficient mice to CS and show an attenuated inflammatory profile in the latter. Alveolar macrophages were sorted by FACS from the different experimental groups and their gene expression profile was analyzed by RNA sequencing. This analysis revealed increased expression of miR-155 targets and an attenuation of the CS-induced increase in inflammation-related genes in miR-155 deficient mice. Moreover, intranasal instillation of a specific miR-155 inhibitor attenuated the CS-induced pulmonary inflammation in mice. Finally, elastase-induced emphysema and lung functional changes were significantly attenuated in miR-155 deficient mice. In conclusion, we highlight a role for miR-155 in CS-induced inflammation and the pathogenesis of COPD, implicating miR-155 as a new therapeutic target in COPD.
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.
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.
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.
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 (Ne) 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 (Ne/Nc) appears very small (10−5), suggesting that overall only a low proportion of adults contribute to the next generation. The high Ne level together with the low Ne/Nc 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.
► We investigated the relation between pollutants bioaccumulation levels and the health status in European eel. ► We measured the associated transcriptional response in liver and gill tissues for ...metallothionein and cytochrome P450. ► We found a large spatial variation in contaminant profile, correlated with low energy reserves and condition. ► Gene transcription patterns reveal a complex response mechanism to a pollutant cocktail. ► Yellow eels show a dysfunctional gene transcription regulation of detoxification genes at higher pollutant levels.
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.