Eric Normandeau

An Improved Genotyping by Sequencing (GBS) Approach Offering Increased Versatility and Efficiency of SNP Discovery and Genotyping

Highly parallel SNP genotyping platforms have been developed for some important crop species, but these platforms typically carry a high cost per sample for first-time or small-scale users. In contrast, recently developed genotyping by sequencing (GBS) approaches offer a highly cost effective alternative for simultaneous SNP discovery and genotyping. In the present investigation, we have explored the use of GBS in soybean. In addition to developing a novel analysis pipeline to call SNPs and indels from the resulting sequence reads, we have devised a modified library preparation protocol to alter the degree of complexity reduction. We used a set of eight diverse soybean genotypes to conduct a pilot scale test of the protocol and pipeline. Using ApeKI for GBS library preparation and sequencing on an Illumina GAIIx machine, we obtained 5.5 M reads and these were processed using our pipeline. A total of 10,120 high quality SNPs were obtained and the distribution of these SNPs mirrored closely the distribution of gene-rich regions in the soybean genome. A total of 39.5% of the SNPs were present in genic regions and 52.5% of these were located in the coding sequence. Validation of over 400 genotypes at a set of randomly selected SNPs using Sanger sequencing showed a 98% success rate. We then explored the use of selective primers to achieve a greater complexity reduction during GBS library preparation. The number of SNP calls could be increased by almost 40% and their depth of coverage was more than doubled, thus opening the door to an increase in the throughput and a significant decrease in the per sample cost. The approach to obtain high quality SNPs developed here will be helpful for marker assisted genomics as well as assessment of available genetic resources for effective utilisation in a wide number of species.

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher FST values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.

The genetic consequences of spatially varying selection in the panmictic American eel (Anguilla rostrata).

Our understanding of the genetic basis of local adaptation has recently benefited from the increased power to identify functional variants associated with environmental variables at the genome scale. However, it often remains challenging to determine whether locally adaptive alleles are actively maintained at intermediate frequencies by spatially varying selection. Here, we evaluate the extent to which this particular type of balancing selection explains the retention of adaptive genetic variation in the extreme situation of perfect panmixia, using the American eel (Anguilla rostrata) as a model. We first conducted a genome scan between two samples from opposite ends of a latitudinal environmental gradient using 454 sequencing of individually tagged cDNA libraries. Candidate SNPs were then genotyped in 992 individuals from 16 sampling sites at different life stages of the same cohort (including larvae from the Sargasso Sea, glass eels, and 1-year-old individuals) as well as in glass eels of the following cohort. Evidence for spatially varying selection was found at 13 loci showing correlations between allele frequencies and environmental variables across the entire species range. Simulations under a multiple-niche Levene’s model using estimated relative fitness values among genotypes rarely predicted a stable polymorphic equilibrium at these loci. Our results suggest that some genetic-by-environment interactions detected in our study arise during the progress toward fixation of a globally advantageous allele with spatially variable effects on fitness.

Mapping phenotypic, expression and transmission ratio distortion QTL using RAD markers in the Lake Whitefish (Coregonus clupeaformis)

The evolution of reproductive isolation in an ecological context may involve multiple facets of species divergence on which divergent selection may operate. These include variation in quantitative phenotypic traits, regulation of gene expression, and differential transmission of particular allelic combinations. Thus, an integrative approach to the speciation process involves identifying the genetic basis of these traits, in order to understand how they are affected by divergent selection in nature and how they ultimately contribute to reproductive isolation. In the Lake Whitefish (Coregonus clupeaformis), dwarf and normal species pairs sympatrically occur in several North American postglacial lakes. The limnetic dwarf whitefish distinguishes from its normal benthic relative by numerous life history, behavioural, morphological and gene expression traits, in relation with the exploitation of distinct ecological niches. Here, we have applied the RAD‐Sequencing method to a hybrid backcross family to reconstruct a high‐density genetic linkage map and perform QTL mapping in the Lake Whitefish. The 3061 cM map encompassed 3438 segregating RAD markers distributed over 40 linkage groups, for an average resolution of 0.89 cM. We mapped phenotypic and expression QTL underlying ecologically important traits as well as transmission ratio distortion QTL, and identified genomic regions harbouring clusters of such QTL. A narrow genomic region strongly associated with sex determination was also evidenced. Positional and functional information revealed in this study will be useful in ongoing population genomic studies to illuminate our understanding of the genomic architecture of reproductive isolation between whitefish species pairs.

Effects of chronic metal exposure on wild fish populations revealed by high-throughput cDNA sequencing

Given the inherent variability of aquatic systems, predicting the in situ effects of contaminants on such ecosystems still represents a major challenge for ecotoxicology. In this context, transcriptomic tools can help identify and investigate the mechanisms of toxicity beyond the traditional morphometric, physiological and population-level endpoints. In this study, we used the 454 sequencing technology to examine the in situ effects of chronic metal (Cd, Cu) exposure on the yellow perch (Perca flavescens) transcriptome. Total hepatic mRNA from fish sampled along a polymetallic gradient was extracted, reverse transcribed, labeled with unique barcode sequences and sequenced. This approach allowed us to identify correlations between the transcription level of single genes and the hepatic concentrations of individual metals; 71% of the correlations established were negative. Chronic metal exposure was thus associated with a decrease in the transcription levels of numerous genes involved in protein biosynthesis, in the immune system, and in lipid and energy metabolism. Our results suggest that this marked decrease could result from an impairment of bile acid metabolism by Cd and energy restriction but also from the recruitment of several genes involved in epigenetic modifications of histones and DNA that lead to gene silencing.

Population-specific gene expression responses to hybridization between farm and wild Atlantic salmon

Because of intrinsic differences in their genetic architectures, wild populations invaded by domesticated individuals could experience population‐specific consequences following introgression by genetic material of domesticated origin. Expression levels of 16 000 transcripts were quantified by microarrays in liver tissue from farm, wild, and farm‐wild backcross (i.e. F1 farm‐wild hybrid × wild; total n = 50) Atlantic salmon (Salmo salar) raised under common environmental conditions. The wild populations and farm strain originated from three North American rivers in eastern Canada (Stewiacke, Tusket, and Saint John rivers, respectively). Analysis of variance revealed 177 transcripts with different expression levels among the five strains compared. Five times more of these transcripts were differentiated between farmed parents and Tusket backcrosses (n = 53) than between Stewiacke backcrosses and their farmed parents (n = 11). Altered biological processes in backcrosses also differed between populations both in number and in the type of processes impacted (metabolism vs immunity). Over‐dominant gene expression regulation in backcrosses varied considerably between populations (23% in Stewiacke vs 44% in Tusket). Hence, the consequences of introgression of farm genetic material on gene expression depended on population‐specific genetic architectures. These results support the need to evaluate impacts of farm‐wild genetic interactions at the population scale.

RNA-seq Reveals Transcriptomic Shock Involving Transposable Elements Reactivation in Hybrids of Young Lake Whitefish Species

Identifying the molecular basis of reproductive isolation among diverging lineages represents an essential step toward understanding speciation in natural populations. Postzygotic barriers can lead to hybrid breakdown, a syndrome that has been documented in several systems, potentially involving the reactivation of transposable elements. In northeastern North America, two lake whitefish lineages have repeatedly colonized postglacial lakes ~12,000 years ago, and a dwarf limnetic species has evolved multiple times from the normal benthic species. Reproductive isolation is incomplete between them; viable hybrids can be generated in the laboratory but significant mortality occurs and is associated with a malformed phenotype in backcross embryos, thus revealing a hybrid breakdown syndrome. By means of RNA-seq analyses, the objective of this study was to determine which genes were misregulated in hybrids and rigorously test the hypothesis of transposable element reactivation. We compared the transcriptomic landscape in pure embryos, F1-hybrids, and healthy and malformed backcrosses at the late embryonic stage. Extensive expression differences consistent with previously documented adaptive divergence between pure normal and dwarf embryos were identified for the first time. Pronounced transcriptome-wide deregulation in malformed backcrosses was observed, with over 15% of transcripts differentially expressed in all comparisons, compared with 1.5% between pure parental forms. Convincing evidence of transposable elements and noncoding transcripts reactivation in malformed backcrosses is presented. We propose that hybrid breakdown likely results from extensive genomic incompatibilities, plausibly encompassing transposable elements. Combined with previous studies, these results reveal synergy among many reproductive barriers, thus maintaining divergence between these two young whitefish species.

Differences in transcription levels among wild, domesticated, and hybrid Atlantic salmon (Salmo salar) from two environments

Escaped domesticated individuals can introduce disadvantageous traits into wild populations due to both adaptive differences between population ancestors and human‐induced changes during domestication. In contrast to their domesticated counterparts, some endangered wild Atlantic salmon populations encounter during their marine stage large amounts of suspended sediments, which may act as a selective agent. We used microarrays to elucidate quantitative transcriptional differences between a domesticated salmon strain, a wild population and their first‐generation hybrids during their marine life stage, to describe transcriptional responses to natural suspended sediments, and to test for adaptive genetic variation in plasticity relating to a history of natural exposure or nonexposure to suspended sediments. We identified 67 genes differing in transcription level among salmon groups. Among these genes, processes related to energy metabolism and ion homoeostasis were over‐represented, while genes contributing to immunity and actin‐/myosin‐related processes were also involved in strain differentiation. Domestic–wild hybrids exhibited intermediate transcription patterns relative to their parents for two‐thirds of all genes that differed between their parents; however, genes deviating from additivity tended to have similar levels to those expressed by the wild parent. Sediments induced increases in transcription levels of eight genes, some of which are known to contribute to external or intracellular damage mitigation. Although genetic variation in plasticity did not differ significantly between groups after correcting for multiple comparisons, two genes (metallothionein and glutathione reductase) tended to be more plastic in response to suspended sediments in wild and hybrid salmon, and merit further examination as candidate genes under natural selection.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs.

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.