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.