Methods for identifying and interpreting sex-linked SNP markers and carrying out sex assignment: application to thornback ray (Raja clavata).

Sex-determining modes remain unknown in numerous species, notably in fishes, in which a variety of modalities have been reported. Additionally, noninvasive individual sexing is problematic for species without external sex attributes or for early life stages, requiring cytogenetic or molecular analyses when sex chromosomes or sex-linked markers have been characterized. Genomics now provide a means to achieve this. Here, we review common sex-determination systems and corresponding statistical methods for identifying sex-linked genetic markers and their use for sex assignment, focusing on single nucleotide polymorphism (SNP) markers derived from reduced representation sequencing methods. We demonstrate the dependence of expected sex assignment error on the number of sex-linked SNPs and minor allele frequency. The application of three methods was made here: (a) identification of heterozygote excess in one sex, (b) FST outlier analysis between the two sexes and (c) neuronal net modelling. These methods were applied to a large SNP data set (4604 SNPs) for 1680 thornback rays (Raja clavata). Using method (a), nineteen putative sex-linked SNPs were identified. Comparison with the reference genome of a related species (Amblyraja radiata) indicated that all 19 SNPs are probably located on the same chromosome. These results suggest that thornback ray has a XX/XY sex-determination system. Method (b) identified eight SNPs probably located on different chromosomes. Method (a) led to the lowest sex assignment error among the three methods (4.2% error for females and 3.7% for males).

Estimating effective population size using RADseq: Effects of SNP selection and sample size.

Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.

Concentrations of mercury and other trace elements in two offshore skates: sandy ray Leucoraja circularis and shagreen ray L. fullonica.

Trace metal concentrations in muscle and liver tissues from two offshore species of skate were examined. Concentrations of mercury in muscle of Leucoraja circularis (n=20; 23-110.5cm total length, 157-490m water depth) and L. fullonica (n=24; 28.5-100cm total length, 130-426m water depth) were 0.02-1.8 and 0.04-0.61mgkg-1, respectively. Concentrations of both As and Hg increased with total length. Only the largest specimen had a concentration of Hg in muscle >1.0mgkg-1. Data were limited for specimens>90cm long, and further studies on contaminants in larger-bodied skates could usefully be undertaken.

Interannual variability of fisheries economic returns and energy ratios is mostly explained by gear type.

According to portfolio theory applied to fisheries management, economic returns are stabilised by harvesting in a portfolio stocks of species whose returns are negatively correlated and for which the portfolio economic return variance is smaller than the sum of stock specific return variances. Also, variability is expected to decrease with portfolio width. Using a range of indicators, these predictions were tested for the French fishing fleets in the Bay of Biscay (Northeast Atlantic) during the period 2001-2009. For this, vessels were grouped into eight fishing fleets based on the gears used and exploited species were grouped into five functional groups. The portfolio width of fleets ranged from 1-3 functional groups, or 4-19 species. Economic fleet returns (sale revenues minus fishing costs) varied strongly between years; the interannual variability was independent of portfolio width (species or functional groups). Energy ratio expressed by the ratio between fuel energy used for fishing and energy contained in landings varied from 0.3 for purse seines to 9.7 for trawlers using bottom trawls alone or in combination with pelagic trawls independent of portfolio width. Interannual variability in total sale revenues was larger than the sum of species specific sales revenue variability, except for fleets using hooks and pelagic trawlers; it increased with the number of species exploited. In conclusion, the interannual variability of economic returns or energy ratios of French fisheries in the Bay of Biscay did not decrease with the number of species or functional groups exploited, though it varied between fleets.

Does presence of a mid-ocean ridge enhance biomass and biodiversity?

In contrast to generally sparse biological communities in open-ocean settings, seamounts and ridges are perceived as areas of elevated productivity and biodiversity capable of supporting commercial fisheries. We investigated the origin of this apparent biological enhancement over a segment of the North Mid-Atlantic Ridge (MAR) using sonar, corers, trawls, traps, and a remotely operated vehicle to survey habitat, biomass, and biodiversity. Satellite remote sensing provided information on flow patterns, thermal fronts, and primary production, while sediment traps measured export flux during 2007-2010. The MAR, 3,704,404 km(2) in area, accounts for 44.7% lower bathyal habitat (800-3500 m depth) in the North Atlantic and is dominated by fine soft sediment substrate (95% of area) on a series of flat terraces with intervening slopes either side of the ridge axis contributing to habitat heterogeneity. The MAR fauna comprises mainly species known from continental margins with no evidence of greater biodiversity. Primary production and export flux over the MAR were not enhanced compared with a nearby reference station over the Porcupine Abyssal Plain. Biomasses of benthic macrofauna and megafauna were similar to global averages at the same depths totalling an estimated 258.9 kt C over the entire lower bathyal north MAR. A hypothetical flat plain at 3500 m depth in place of the MAR would contain 85.6 kt C, implying an increase of 173.3 kt C attributable to the presence of the Ridge. This is approximately equal to 167 kt C of estimated pelagic biomass displaced by the volume of the MAR. There is no enhancement of biological productivity over the MAR; oceanic bathypelagic species are replaced by benthic fauna otherwise unable to survive in the mid ocean. We propose that globally sea floor elevation has no effect on deep sea biomass; pelagic plus benthic biomass is constant within a given surface productivity regime.