Landscape genomics of the American lobster (Homarus americanus).

In marine species experiencing intense fishing pressures, knowledge of genetic structure and local adaptation represent a critical information to assist sustainable management. In this study, we performed a landscape genomics analysis in the American lobster to investigate the issues pertaining to the consequences of making use of putative adaptive loci to reliably infer population structure and thus more rigorously delineating biological management units in marine exploited species. Toward this end, we genotyped 14,893 single nucleotide polymorphism (SNPs) in 4190 lobsters sampled across 96 sampling sites distributed along 1000 km in the northwest Atlantic in both Canada and the USA. As typical for most marine species, we observed a weak, albeit highly significant genetic structure. We also found that adaptive genetic variation allows detecting fine-scale population structure not resolved by neutral genetic variation alone. Using the recent genome assembly of the American lobster, we were able to map and annotate several SNPs located in functional genes potentially implicated in adaptive processes such as thermal stress response, salinity tolerance and growth metabolism pathways. Taken together, our study indicates that weak population structure in high gene flow systems can be resolved at various spatial scales, and that putatively adaptive genetic variation can substantially enhance the delineation of biological management units of marine exploited species.

Sex matters in massive parallel sequencing: Evidence for biases in genetic parameter estimation and investigation of sex determination systems.

Using massively parallel sequencing data from two species with different life history traits, American lobster (Homarus americanus) and Arctic Char (Salvelinus alpinus), we highlight how an unbalanced sex ratio in the samples and a few sex-linked markers may lead to false interpretations of population structure and thus to potentially erroneous management recommendations. Here, multivariate analyses revealed two genetic clusters separating samples by sex instead of by expected spatial variation: inshore and offshore locations in lobster, or east and west locations in Arctic Char. To further investigate this, we created several subsamples artificially varying the sex ratio in the inshore/offshore and east/west groups and then demonstrated that significant genetic differentiation could be observed despite panmixia in lobster, and that FST values were overestimated in Arctic Char. This pattern was due to 12 and 94 sex-linked markers driving differentiation for lobster and Arctic Char, respectively. Removing sex-linked markers led to nonsignificant genetic structure in lobster and a more accurate estimation of FST in Arctic Char. The locations of these markers and putative identities of genes containing or nearby the markers were determined using available transcriptomic and genomic data, and this provided new information related to sex determination in both species. Given that only 9.6% of all marine/diadromous population genomic studies to date have reported sex information, we urge researchers to collect and consider individual sex information. Sex information is therefore relevant for avoiding unexpected biases due to sex-linked markers as well as for improving our knowledge of sex determination systems in nonmodel species.

Seascape genomics provides evidence for thermal adaptation and current-mediated population structure in American lobster (Homarus americanus).

Investigating how environmental features shape the genetic structure of populations is crucial for understanding how they are potentially adapted to their habitats, as well as for sound management. In this study, we assessed the relative importance of spatial distribution, ocean currents and sea surface temperature (SST) on patterns of putatively neutral and adaptive genetic variation among American lobster from 19 locations using population differentiation (PD) approaches combined with environmental association (EA) analyses. First, PD approaches (using bayescan, arlequin and outflank) found 28 outlier SNPs putatively under divergent selection and 9770 neutral SNPs in common. Redundancy analysis revealed that spatial distribution, ocean current-mediated larval connectivity and SST explained 31.7% of the neutral genetic differentiation, with ocean currents driving the majority of this relationship (21.0%). After removing the influence of spatial distribution, no SST were significant for putatively neutral genetic variation whereas minimum annual SST still had a significant impact and explained 8.1% of the putatively adaptive genetic variation. Second, EA analyses (using Pearson correlation tests, bayescenv and lfmm) jointly identified seven SNPs as candidates for thermal adaptation. Covariation at these SNPs was assessed with a spatial multivariate analysis that highlighted a significant temperature association, after accounting for the influence of spatial distribution. Among the 505 candidate SNPs detected by at least one of the three approaches, we discovered three polymorphisms located in genes previously shown to play a role in thermal adaptation. Our results have implications for the management of the American lobster and provide a foundation on which to predict how this species will cope with climate change.

RAD genotyping reveals fine-scale genetic structuring and provides powerful population assignment in a widely distributed marine species, the American lobster (Homarus americanus).

Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species' range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species' natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (FCT  = 0.0011; P-value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean FST  = 0.00185; CI: 0.0007-0.0021, P-value < 0.0002), providing strong evidence for weak, albeit fine-scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest FST . Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid 'high-grading bias', 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine-scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.

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.

Geographic variation of multiple paternity in the American lobster, Homarus americanus.

We studied the frequency of multiple paternity for American lobster (Homarus americanus) at three Canadian sites differing in exploitation rate and mean adult size. The probability of detecting multiple paternity using four microsatellite loci and 100 eggs per female was in excess of 99% under various scenarios of paternal contribution. Overall, 13% of the 108 examined females carried a clutch sired by two or three males. Multiple paternity was observed at the two most exploited sites (11% at Magdalen Islands and 28% at Grand Manan Island), whereas single paternity only was observed at the least exploited site (Anticosti Island). Within populations females with a clutch sired by more than one male tended to be smaller than females with a clutch sired by a single male. Based on these and other findings, we postulate a link between female promiscuity and sperm limitation in the American lobster.