Disparate population and holobiont structure of pocilloporid corals across the Red Sea gradient demonstrate species-specific evolutionary trajectories.

Global habitat degradation heightens the need to better understand patterns of genetic connectivity and diversity of marine biota across geographical ranges to guide conservation efforts. Corals across the Red Sea are subject to pronounced environmental differences, but studies so far suggest that animal populations are largely connected, excepting evidence for a genetic break between the northern-central and southern regions. Here, we investigated population structure and holobiont assemblage of two common pocilloporid corals, Pocillopora verrucosa and Stylophora pistillata, across the Red Sea. We found little evidence for population differentiation in P. verrucosa, except for the southernmost site. Conversely, S. pistillata exhibited a complex population structure with evidence for within-reef and regional genetic differentiation, in line with differences in their reproductive mode (P. verrucosa is a broadcast spawner and S. pistillata is a brooder). Analysis for genomic loci under positive selection identified 85 sites (18 of which were in coding sequences) that distinguished the southern P. verrucosa population from the remainder of the Red Sea population. By comparison, we found 128 loci (24 of which were residing in coding sequences) in S. pistillata with evidence for local adaptation at various sites. Functional annotation of the underlying proteins revealed putative roles in the response to stress, lipid metabolism, transport, cytoskeletal rearrangement, and ciliary function (among others). Microbial assemblages of both coral species showed pervasive association with microalgal symbionts from the genus Symbiodinium (former clade A) and bacteria from the genus Endozoicomonas that exhibited significant differences according to host genotype and environment. The disparity of population genetic and holobiont assemblage patterns even between closely related species (family Pocilloporidae) highlights the need for multispecies investigations to better understand the role of the environment in shaping evolutionary trajectories. It further emphasizes the importance of networks of reef reserves to achieve conservation of genetic variants critical to the future survival of coral ecosystems.

Low levels of sibship encourage use of larvae in western Atlantic bluefin tuna abundance estimation by close-kin mark-recapture.

Globally, tunas are among the most valuable fish stocks, but are also inherently difficult to monitor and assess. Samples of larvae of Western Atlantic bluefin tuna Thunnus thynnus (Linnaeus, 1758) from standardized annual surveys in the northern Gulf of Mexico provide a potential source of "offspring" for close-kin mark-recapture (CKMR) estimates of abundance. However, the spatial patchiness and highly skewed numbers of larvae per tow suggest sampled larvae may come from a small number of parents, compromising the precision of CKMR. We used high throughput genomic profiling to study sibship within and among larval tows from the 2016 standardized Gulf-wide survey compared to targeted sampling carried out in 2017. Full- and half-siblings were found within both years, with 12% of 156 samples in 2016 and 56% of 317 samples in 2017 having at least one sibling. There were also two pairs of cross cohort half-siblings. Targeted sampling increased the number of larvae collected per sampling event but resulted in a higher proportion of siblings. The combined effective sample size across both years was about 75% of the nominal size, indicating that Gulf of Mexico larval collections could be a suitable source of juveniles for CKMR in Western Atlantic bluefin tuna.

Rethinking the influence of hydroelectric development on gene flow in a long-lived fish, the Lake Sturgeon Acipenser fulvescens.

Many hydroelectric dams have been in place for 50 - >100 years, which for most fish species means that enough generations have passed for fragmentation induced divergence to have accumulated. However, for long-lived species such as Lake Sturgeon, Acipenser fulvescens, it should be possible to discriminate between historical population structuring and contemporary gene flow and improve the broader understanding of anthropogenic influence. On the Winnipeg River, Manitoba, two hypotheses were tested: 1) Measureable quantities of former reservoir dwelling Lake Sturgeon now reside downstream of the Slave Falls Generating Station, and 2) genetically differentiated populations of Lake Sturgeon occur upstream and downstream, a result of historical structuring. Genetic methods based on ten microsatellite markers were employed, and simulations were conducted to provide context. With regards to contemporary upstream to downstream contributions, the inclusion of length-at-age data proved informative. Both pairwise relatedness and Bayesian clustering analysis substantiated that fast-growing outliers, apparently entrained after residing in the upstream reservoir for several years, accounted for ~15% of the Lake Sturgeon 525-750 mm fork length captured downstream. With regards to historical structuring, upstream and downstream populations were found to be differentiated (FST = 0.011, and 0.013-0.014 when fast-growing outliers were excluded), and heterozygosity metrics were higher for downstream versus upstream juveniles. Historical asymmetric (downstream) gene flow in the vicinity of the generating station was the most logical explanation for the observed genetic structuring. In this section of the Winnipeg River, construction of a major dam does not appear to have fragmented a previously panmictic Lake Sturgeon population, but alterations to habitat may be influencing upstream to downstream contributions in unexpected ways.

Linking genomics and population genetics with R.

Population genetics and genomics have developed and been treated as independent fields of study despite having common roots. The continuous progress of sequencing technologies is contributing to (re-)connect these two disciplines. We review the challenges faced by data analysts and software developers when handling very big genetic data sets collected on many individuals. We then expose how r, as a computing language and development environment, proposes some solutions to meet these challenges. We focus on some specific issues that are often encountered in practice: handling and analysing single-nucleotide polymorphism data, handling and reading variant call format files, analysing haplotypes and linkage disequilibrium and performing multivariate analyses. We illustrate these implementations with some analyses of three recently published data sets that contain between 60 000 and 1 000 000 loci. We conclude with some perspectives on future developments of r software for population genomics.

Salmonid Chromosome Evolution as Revealed by a Novel Method for Comparing RADseq Linkage Maps.

Whole genome duplication (WGD) can provide material for evolutionary innovation. Family Salmonidae is ideal for studying the effects of WGD as the ancestral salmonid underwent WGD relatively recently, ∼65 Ma, then rediploidized and diversified. Extensive synteny between homologous chromosome arms occurs in extant salmonids, but each species has both conserved and unique chromosome arm fusions and fissions. Assembly of large, outbred eukaryotic genomes can be difficult, but structural rearrangements within such taxa can be investigated using linkage maps. RAD sequencing provides unprecedented ability to generate high-density linkage maps for nonmodel species, but can result in low numbers of homologous markers between species due to phylogenetic distance or differences in library preparation. Here, we generate a high-density linkage map (3,826 markers) for the Salvelinus genera (Brook Charr S. fontinalis), and then identify corresponding chromosome arms among the other available salmonid high-density linkage maps, including six species of Oncorhynchus, and one species for each of Salmo, Coregonus, and the nonduplicated sister group for the salmonids, Northern Pike Esox lucius for identifying post-duplicated homeologs. To facilitate this process, we developed MapComp to identify identical and proximate (i.e. nearby) markers between linkage maps using a reference genome of a related species as an intermediate, increasing the number of comparable markers between linkage maps by 5-fold. This enabled a characterization of the most likely history of retained chromosomal rearrangements post-WGD, and several conserved chromosomal inversions. Analyses of RADseq-based linkage maps from other taxa will also benefit from MapComp, available at: https://github.com/enormandeau/mapcomp/

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.

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.