Genetic diversity patterns in farmed rainbow trout (Oncorhynchus mykiss) populations using genome-wide SNP and haplotype data.

Rainbow trout is one of the most popular aquaculture species worldwide, with a long history of domestication. However, limited information exists about the genetic diversity of farmed rainbow trout populations globally, with most available reports relying on low-throughput genotyping technologies. Notably, no information exists about the genetic diversity status of farmed rainbow trout in Sweden. Double-digest restriction-site-associated DNA sequencing was performed on more than 500 broodfish from two leading producers in Sweden and from the country's national breeding program. Following the detection of single nucleotide polymorphisms (SNPs), genetic diversity was studied by using either individual SNPs (n = 8680; one SNP retained per 300 bp sequence reads) or through SNP haplotypes (n = 20 558; all SNPs retained in 300 bp sequence reads). Similar amounts of genetic diversity were found amongst the three populations when individual SNPs were used. Furthermore, principal component analysis and discriminant analysis of principal components suggested two genetic clusters with the two industry populations grouped together. Genetic differentiation based on the FST fixation index was ~0.01 between the industry populations and ~0.05 when those were compared with the breeding program. Preliminary estimates of effective population size (Ne ) and inbreeding (based on runs of homozygosity; FROH ) were similar amongst the three populations (Ne ≈ 50-80; median FROH ≈ 0.11). Finally, the haplotype-based analysis suggested that animals from the breeding program had higher shared coancestry levels than those from the other two populations. Overall, our study provides novel insights into the genetic diversity and structure of Sweden's three main farmed rainbow trout populations, which could guide their future management.

Seasonal and age-related changes in sperm quality of farmed arctic charr (Salvelinus alpinus).

BACKGROUND: Substantial variation in male fertility is regularly observed in farmed Arctic charr. However, detailed investigations of its fluctuation during a reproductive season and across years are lacking. Furthermore, information about the effect of underlying genetic factors influencing sperm quality is scarce. The current study focused on seasonal and age-related factors that may affect sperm quality characteristics in males reared in natural and delayed photoperiods. Animals were sampled three times for two consecutive years, and sperm quality parameters were recorded using a computer-assisted sperm analysis (CASA) system. Thereafter, high-throughput sequencing technologies were applied, aiming to identify genomic regions related to the variation of sperm quality throughout the reproductive season. RESULTS: An across-season variation in the recorded sperm quality parameters was evident. Overall, 29% and 42% of males from the natural and delayed spawning groups had a highly variable total progressive motility. Males at four years of age showed significantly higher sperm motility and velocities during the early October and November recordings compared to the following year when the same animals were five years of age. On the other hand, the opposite was observed regarding sperm concentration during the last sampling. A genome-wide FST scan detected SNP differentiation among males with high and low variability in total progressive motility (PM) on eight chromosomes (FST > 0.17), Genome wide windows with the highest FST contained SNPs in proximity (within 250 kb up- and downstream distance) to 16 genes with sperm quality biological functions in mammalian species. CONCLUSION: Our findings provide a detailed view of seasonal, age-related, and genetic effects on sperm quality and can be used to guide decisions on broodstock selection and hatchery management.

Whole-genome re-sequencing provides key genomic insights in farmed Arctic charr (Salvelinus alpinus) populations of anadromous and landlocked origin from Scandinavia.

Arctic charr (Salvelinus alpinus) is a niche-market high-value species for Nordic aquaculture. Similar to other salmonids, both anadromous and landlocked populations are encountered. Whole-genome re-sequencing (22X coverage) was performed on two farmed populations of anadromous (Sigerfjord; n = 24) and landlocked (Arctic Superior; n = 24) origin from Norway and Sweden respectively. More than 5 million SNPs were used to study their genetic diversity and to scan for selection signatures. The two populations were clearly distinguished through principal component analysis, with the mean fixation index being ~0.12. Furthermore, the levels of genomic inbreeding estimated from runs of homozygosity were 6.23% and 8.66% for the Norwegian and the Swedish population respectively. Biological processes that could be linked to selection pressure associated primarily with the anadromous background and/or secondarily with domestication were suggested. Overall, our study provided insights regarding the genetic composition of two main strains of farmed Arctic charr from Scandinavia. At the same time, ample genomic resources were produced in the magnitude of millions of SNPs that could assist the transition of Nordic Arctic charr farming in the genomics era.

The use of telomeric length as authenticity marker in fish and seafood - a new perspective in the detection of adulteration.

In this review we propose the use of telomeric length (TL) as an authenticity marker that could provide an alternative method for differentiating fish and seafood samples or detecting fraud. Considering the ever-growing number of incidents of economically motivated fish and seafood adulteration using even more sophisticated methods to overcome current authenticity markers, the need to identify novel authenticity markers becomes essential. The TL of fish and seafood depends on individual characteristics (e.g., sex, age) and the environmental stimuli (e.g., temperature, water quality) to which these are exposed. Hence, both wild marine and freshwater populations occupying different geographical origin habitats might differ substantially because of the environmental cues affecting them. Moreover, the implementation of various rearing practices in aquaculture, such as different levels of fish and seafood density and increased ambient noise combined with site-specific environmental cues could affect TL, providing regulatory authorities with valuable information by distinguishing wild from reared populations and organic from conventional ones. In the present review the effects of both the environmental conditions and individual characteristics on the telomeric stability of fish and seafood telomeres are discussed, suggesting TL as a potential prospect authenticity marker that could be used to prevent fish and seafood adulteration.

A genomic-based vision on the genetic diversity and key performance traits in selectively bred Arctic charr (Salvelinus alpinus).

Routine implementation of genomic information for guiding selection decisions is not yet common in the majority of aquaculture species. Reduced representation sequencing approaches offer a cost-effective solution for obtaining genome-wide information in species with a limited availability of genomic resources. In the current study, we implemented double-digest restriction site-associated DNA sequencing (ddRAD-seq) on an Arctic charr strain with the longest known history of selection (approximately 40 years) aiming to improve selection decisions. In total, 1730 animals reared at four different farms in Sweden and spanning from year classes 2013-2017 were genotyped using ddRAD-seq. Approximately 5000 single nucleotide polymorphisms (SNPs) were identified, genetic diversity-related metrics were estimated, and genome-wide association studies (GWAS) for body length at different time points and age of sexual maturation were conducted. Low genetic differentiation amongst animals from the different farms was observed based on both the results from pairwise F st values and principal component analysis (PCA). The existence of associations was investigated between the mean genome-wide heterozygosity of each full-sib family (year class 2017) and the corresponding inbreeding coefficient or survival to the eyed stage. A moderate correlation (-0.33) was estimated between the mean observed heterozygosity of each full-sib family and the corresponding inbreeding coefficient, while no linear association was obtained with the survival to the eyed stage. GWAS did not detect loci with major effect for any of the studied traits. However, genomic regions explaining more than 1% of the additive genetic variance for either studied traits were suggested across 14 different chromosomes. Overall, key insights valuable for future selection decisions of Arctic charr have been obtained, suggesting ddRAD as an attractive genotyping platform for obtaining genome-wide information in a cost-effective manner.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes.

BACKGROUND: The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations. METHODS: Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time. RESULTS: Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe. CONCLUSIONS: Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

Growth performance of five different strains of Nile tilapia (Oreochromis niloticus) introduced to Tanzania reared in fresh and brackish waters.

Five introduced strains of Nile tilapia (Oreochromis niloticus) were tested for growth performance both in fresh- and brackish-water (2 salinity units) environments for 56 days. The BIG NIN, GIFT, Chitralada, "Ruvu Farm" and Silver YY strains with initial mean average weight (± standard error) of 96.4 ± 6.90 g, 104.1 ± 7.19 g, 137.2 ± 7.21 g, 53.2 ± 6.98 g and 95.3 ± 7.11 g, respectively were used. Individuals were tagged and pooled in hapas (12 m × 8.5 m × 2 m each), aligned into different ponds (20 m × 20 m each). Stocking density of 5 fish/m2 and 350 g/kg crude protein diet were used. Overall, the average weight gain for GIFT strain was 7.5%, 32%, 45% and 86.5% higher than BIG NIN, Chitralada, "Ruvu Farm" and Silver YY strains, respectively, across both environments. All strains performed significantly better (p < 0.05) when reared in brackish-water than their respective counterparts in freshwater, except for the BIG NIN strain. The morphometric correlations for all strains in both environments ranged from moderate (0.50) to strong positive (0.92). The GIFT strain demonstrated superior growth and genotype by environment interaction was weak and not important to be prioritized in breeding programs.

Genotyping Strategies Using ddRAD Sequencing in Farmed Arctic Charr (Salvelinus alpinus).

Incorporation of genomic technologies into fish breeding programs is a modern reality, promising substantial advances regarding the accuracy of selection, monitoring the genetic diversity and pedigree record verification. Single nucleotide polymorphism (SNP) arrays are the most commonly used genomic tool, but the investments required make them unsustainable for emerging species, such as Arctic charr (Salvelinus alpinus), where production volume is low. The requirement to genotype a large number of animals for breeding practices necessitates cost effective genotyping approaches. In the current study, we used double digest restriction site-associated DNA (ddRAD) sequencing of either high or low coverage to genotype Arctic charr from the Swedish national breeding program and performed analytical procedures to assess their utility in a range of tasks. SNPs were identified and used for deciphering the genetic structure of the studied population, estimating genomic relationships and implementing an association study for growth-related traits. Missing information and underestimation of heterozygosity in the low coverage set were limiting factors in genetic diversity and genomic relationship analyses, where high coverage performed notably better. On the other hand, the high coverage dataset proved to be valuable when it comes to identifying loci that are associated with phenotypic traits of interest. In general, both genotyping strategies offer sustainable alternatives to hybridization-based genotyping platforms and show potential for applications in aquaculture selective breeding.

Assessing the potential of improving growth and survival to the eyed stage in selectively bred Arctic charr (Salvelinus alpinus).

The Arctic charr breeding programme has been a main driving force for developing the aquaculture industry in Sweden. Selection has been performed for almost 40 years using animals from a closed breeding nucleus. The aim of the current study was to evaluate the potential of further improving growth-related traits taking into account the existence of genotype-by-environment interaction. Furthermore, we investigated the magnitude of the genetic component associated with survival to the eyed stage and potential associations with inbreeding coefficients. A preliminary heritability estimate of 0.23 (SE 0.20) was obtained for survival to the eyed stage using records spanning from 2000 to 2017 (n = 230). Moreover, moderate-to-high heritability estimates (0.27-0.49) were obtained for growth-related traits (body weight and length), using animals from the latest generation of selection (year class 2017). Those animals (n = 2,776), originating from 55 full-sib families, were split into two groups and reared in separate land-based facilities of commercial fish farms in Sweden. The growth-related traits were recorded twice in both sites when animals were of >1 and >2 years of age. Existence of sexual growth dimorphism was indicated with the males having on average 6%-8% higher total length and 22%-34% higher body weight. Furthermore, high genetic correlations regarding growth traits were obtained amongst animals reared at the two different sites (0.82-0.95). In addition, we assessed the accuracy of best linear unbiased prediction (BLUP)-derived estimated breeding values (EBVs) when phenotypes from each rearing site were subsequently masked and used as a validation set. A mean prediction accuracy of 0.60 (length) and 0.64 (weight) were derived for both rearing sites. Overall, our results suggest that further growth improvements should be possible in the subsequent generations of selection. Finally, even though indications for the existence of an underlying genetic component(s) involved in survival to the eyed stage were obtained additional data will be required for elucidating its magnitude.

Assessing the genetic diversity of farmed and wild Rufiji tilapia (Oreochromis urolepis urolepis) populations using ddRAD sequencing.

Rufiji tilapia (Oreochromis urolepis urolepis) is an endemic cichlid in Tanzania. In addition to its importance for biodiversity conservation, Rufiji tilapia is also attractive for farming due to its high growth rate, salinity tolerance, and the production of all-male hybrids when crossed with Nile tilapia (Oreochromis niloticus). The aim of the current study was to assess the genetic diversity and population structure of both wild and farmed Rufiji tilapia populations in order to inform conservation and aquaculture practices. Double-digest restriction-site-associated DNA (ddRAD) libraries were constructed from 195 animals originating from eight wild (Nyamisati, Utete, Mansi, Mindu, Wami, Ruaha, Kibasira, and Kilola) and two farmed (Bwawani and Chemchem) populations. The identified single nucleotide polymorphisms (SNPs; n = 2,182) were used to investigate the genetic variation within and among the studied populations. Genetic distance estimates (F st) were low among populations from neighboring locations, with the exception of Utete and Chemchem populations (F st = 0.34). Isolation-by-distance (IBD) analysis among the wild populations did not detect any significant correlation signal (r = .05; p-value = .4) between the genetic distance and the sampling (Euclidean distance) locations. Population structure and putative ancestry were further investigated using both Bayesian (Structure) and multivariate approaches (discriminant analysis of principal components). Both analysis indicated the existence of three distinct genetic clusters. Two cross-validation scenarios were conducted in order to test the efficiency of the SNP dataset for discriminating between farmed and wild animals or predicting the population of origin. Approximately 95% of the test dataset was correctly classified in the first scenario, while in the case of predicting for the population of origin 68% of the test dataset was correctly classified. Overall, our results provide novel insights regarding the population structure of Rufiji tilapia and a new database of informative SNP markers for both conservation management and aquaculture activities.

Application of Low Coverage Genotyping by Sequencing in Selectively Bred Arctic Charr (Salvelinus alpinus).

Arctic charr (Salvelinus alpinus) is a species of high economic value for the aquaculture industry, and of high ecological value due to its Holarctic distribution in both marine and freshwater environments. Novel genome sequencing approaches enable the study of population and quantitative genetic parameters even on species with limited or no prior genomic resources. Low coverage genotyping by sequencing (GBS) was applied in a selected strain of Arctic charr in Sweden originating from a landlocked freshwater population. For the needs of the current study, animals from year classes 2013 (171 animals, parental population) and 2017 (759 animals; 13 full sib families) were used as a template for identifying genome wide single nucleotide polymorphisms (SNPs). GBS libraries were constructed using the PstI and MspI restriction enzymes. Approximately 14.5K SNPs passed quality control and were used for estimating a genomic relationship matrix. Thereafter a wide range of analyses were conducted in order to gain insights regarding genetic diversity and investigate the efficiency of the genomic information for parentage assignment and breeding value estimation. Heterozygosity estimates for both year classes suggested a slight excess of heterozygotes. Furthermore, FST estimates among the families of year class 2017 ranged between 0.009 - 0.066. Principal components analysis (PCA) and discriminant analysis of principal components (DAPC) were applied aiming to identify the existence of genetic clusters among the studied population. Results obtained were in accordance with pedigree records allowing the identification of individual families. Additionally, DNA parentage verification was performed, with results in accordance with the pedigree records with the exception of a putative dam where full sib genotypes suggested a potential recording error. Breeding value estimation for juvenile growth through the usage of the estimated genomic relationship matrix clearly outperformed the pedigree equivalent in terms of prediction accuracy (0.51 opposed to 0.31). Overall, low coverage GBS has proven to be a cost-effective genotyping platform that is expected to boost the selection efficiency of the Arctic charr breeding program.

Species-Specific Marker Discovery in Tilapia.

Tilapias (family Cichlidae) are of importance in aquaculture and fisheries. Hybridisation and introgression are common within tilapia genera but are difficult to analyse due to limited numbers of species-specific genetic markers. We tested the potential of double digested restriction-site associated DNA (ddRAD) sequencing for discovering single nucleotide polymorphism (SNP) markers to distinguish between 10 tilapia species. Analysis of ddRAD data revealed 1,371 shared SNPs in the de novo-based analysis and 1,204 SNPs in the reference-based analysis. Phylogenetic trees based on these two analyses were very similar. A total of 57 species-specific SNP markers were found among the samples analysed of the 10 tilapia species. Another set of 62 species-specific SNP markers was identified from a subset of four species which have often been involved in hybridisation in aquaculture: 13 for Oreochromis niloticus, 23 for O. aureus, 12 for O. mossambicus and 14 for O. u. hornorum. A panel of 24 SNPs was selected to distinguish among these four species and validated using 91 individuals. Larger numbers of SNP markers were found that could distinguish between the pairs of species within this subset. This technique offers potential for the investigation of hybridisation and introgression among tilapia species in aquaculture and in wild populations.

Optimizing Genomic Prediction of Host Resistance to Koi Herpesvirus Disease in Carp.

Genomic selection (GS) is increasingly applied in breeding programs of major aquaculture species, enabling improved prediction accuracy and genetic gain compared to pedigree-based approaches. Koi Herpesvirus disease (KHVD) is notifiable by the World Organization for Animal Health and the European Union, causing major economic losses to carp production. GS has potential to breed carp with improved resistance to KHVD, thereby contributing to disease control. In the current study, Restriction-site Associated DNA sequencing (RAD-seq) was applied on a population of 1,425 common carp juveniles which had been challenged with Koi herpes virus, followed by sampling of survivors and mortalities. GS was tested on a wide range of scenarios by varying both SNP densities and the genetic relationships between training and validation sets. The accuracy of correctly identifying KHVD resistant animals using GS was between 8 and 18% higher than pedigree best linear unbiased predictor (pBLUP) depending on the tested scenario. Furthermore, minor decreases in prediction accuracy were observed with decreased SNP density. However, the genetic relationship between the training and validation sets was a key factor in the efficacy of genomic prediction of KHVD resistance in carp, with substantially lower prediction accuracy when the relationships between the training and validation sets did not contain close relatives.

Population Structure and Genetic Diversity of Nile Tilapia (Oreochromis niloticus) Strains Cultured in Tanzania.

Understanding population structure and genetic diversity within and between local Nile tilapia lines cultured in Tanzania is important for sustainable aquaculture production. This study investigated the genetic structure and diversity among seven Nile tilapia populations in Tanzania (Karanga, Igunga, Ruhila, Fisheries Education and Training Agency, Tanzania Fisheries Research Institute, Kunduchi, and Lake Victoria). Double-digest restriction site-associated DNA (ddRAD) libraries were prepared from 140 individual fish (20 per population) and sequenced using an Illumina HiSeq 4000 resulting in the identification of 2,180 informative single nucleotide polymorphisms (SNPs). Pairwise Fst values revealed strong genetic differentiation between the closely related populations; FETA, Lake Victoria, and Igunga and those from TAFIRI and Karanga with values ranging between 0.45 and 0.55. Population structure was further evaluated using Bayesian model-based clustering (STRUCTURE) and discriminant analysis of principal components (DAPC). Admixture was detected among Karanga, Kunduchi, and Ruhila populations. A cross-validation approach (25% of individual fish from each population was considered of unknown origin) was conducted in order to test the efficiency of the SNP markers to correctly assign individual fish to the population of origin. The cross-validation procedure was repeated 10 times resulting in 77% of the tested individual fish being allocated to the correct population. Overall our results provide a new database of informative SNP markers for both conservation management and aquaculture activities of Nile tilapia strains in Tanzania.

Applications of genotyping by sequencing in aquaculture breeding and genetics.

Selective breeding is increasingly recognized as a key component of sustainable production of aquaculture species. The uptake of genomic technology in aquaculture breeding has traditionally lagged behind terrestrial farmed animals. However, the rapid development and application of sequencing technologies has allowed aquaculture to narrow the gap, leading to substantial genomic resources for all major aquaculture species. While high-density single-nucleotide polymorphism (SNP) arrays for some species have been developed recently, direct genotyping by sequencing (GBS) techniques have underpinned many of the advances in aquaculture genetics and breeding to date. In particular, restriction-site associated DNA sequencing (RAD-Seq) and subsequent variations have been extensively applied to generate population-level SNP genotype data. These GBS techniques are not dependent on prior genomic information such as a reference genome assembly for the species of interest. As such, they have been widely utilized by researchers and companies focussing on nonmodel aquaculture species with relatively small research communities. Applications of RAD-Seq techniques have included generation of genetic linkage maps, performing genome-wide association studies, improvements of reference genome assemblies and, more recently, genomic selection for traits of interest to aquaculture like growth, sex determination or disease resistance. In this review, we briefly discuss the history of GBS, the nuances of the various GBS techniques, bioinformatics approaches and application of these techniques to various aquaculture species.

Mapping and Sequencing of a Significant Quantitative Trait Locus Affecting Resistance to Koi Herpesvirus in Common Carp.

Cyprinids are the most highly produced group of fishes globally, with common carp being one of the most valuable species of the group. Koi herpesvirus (KHV) infections can result in high levels of mortality, causing major economic losses, and is listed as a notifiable disease by the World Organization for Animal Health. Selective breeding for host resistance has the potential to reduce morbidity and losses due to KHV. Therefore, improving knowledge about host resistance and methods of incorporating genomic data into breeding for resistance may contribute to a decrease in economic losses in carp farming. In the current study, a population of 1,425 carp juveniles, originating from a factorial cross between 40 sires and 20 dams was challenged with KHV. Mortalities and survivors were recorded and sampled for genotyping by sequencing using Restriction Site-Associated DNA sequencing (RADseq). Genome-wide association analyses were performed to investigate the genetic architecture of resistance to KHV. A genome-wide significant QTL affecting resistance to KHV was identified on linkage group 44, explaining approximately 7% of the additive genetic variance. Pooled whole genome resequencing of a subset of resistant (n = 60) and susceptible animals (n = 60) was performed to characterize QTL regions, including identification of putative candidate genes and functional annotation of associated polymorphisms. The TRIM25 gene was identified as a promising positional and functional candidate within the QTL region of LG 44, and a putative premature stop mutation in this gene was discovered.

Genome-wide association and genomic prediction of resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax) using RAD sequencing.

BACKGROUND: European sea bass (Dicentrarchus labrax) is one of the most important species for European aquaculture. Viral nervous necrosis (VNN), commonly caused by the redspotted grouper nervous necrosis virus (RGNNV), can result in high levels of morbidity and mortality, mainly during the larval and juvenile stages of cultured sea bass. In the absence of efficient therapeutic treatments, selective breeding for host resistance offers a promising strategy to control this disease. Our study aimed at investigating genetic resistance to VNN and genomic-based approaches to improve disease resistance by selective breeding. A population of 1538 sea bass juveniles from a factorial cross between 48 sires and 17 dams was challenged with RGNNV with mortalities and survivors being recorded and sampled for genotyping by the RAD sequencing approach. RESULTS: We used genome-wide genotype data from 9195 single nucleotide polymorphisms (SNPs) for downstream analysis. Estimates of heritability of survival on the underlying scale for the pedigree and genomic relationship matrices were 0.27 (HPD interval 95%: 0.14-0.40) and 0.43 (0.29-0.57), respectively. Classical genome-wide association analysis detected genome-wide significant quantitative trait loci (QTL) for resistance to VNN on chromosomes (unassigned scaffolds in the case of 'chromosome' 25) 3, 20 and 25 (P < 1e06). Weighted genomic best linear unbiased predictor provided additional support for the QTL on chromosome 3 and suggested that it explained 4% of the additive genetic variation. Genomic prediction approaches were tested to investigate the potential of using genome-wide SNP data to estimate breeding values for resistance to VNN and showed that genomic prediction resulted in a 13% increase in successful classification of resistant and susceptible animals compared to pedigree-based methods, with Bayes A and Bayes B giving the highest predictive ability. CONCLUSIONS: Genome-wide significant QTL were identified but each with relatively small effects on the trait. Tests of genomic prediction suggested that incorporating genome-wide SNP data is likely to result in higher accuracy of estimated breeding values for resistance to VNN. RAD sequencing is an effective method for generating such genome-wide SNPs, and our findings highlight the potential of genomic selection to breed farmed European sea bass with improved resistance to VNN.

Accuracy of Genomic Evaluations of Juvenile Growth Rate in Common Carp (Cyprinus carpio) Using Genotyping by Sequencing.

Cyprinids are the most important group of farmed fish globally in terms of production volume, with common carp (Cyprinus carpio) being one of the most valuable species of the group. The use of modern selective breeding methods in carp is at a formative stage, implying a large scope for genetic improvement of key production traits. In the current study, a population of 1,425 carp juveniles, originating from a partial factorial cross between 40 sires and 20 dams, was used for investigating the potential of genomic selection (GS) for juvenile growth, an exemplar polygenic production trait. RAD sequencing was used to identify and genotype SNP markers for subsequent parentage assignment, construction of a medium density genetic map (12,311 SNPs), genome-wide association study (GWAS), and testing of GS. A moderate heritability was estimated for body length of carp at 120 days (as a proxy of juvenile growth) of 0.33 (s.e. 0.05). No genome-wide significant QTL was identified using a single marker GWAS approach. Genomic prediction of breeding values outperformed pedigree-based prediction, resulting in 18% improvement in prediction accuracy. The impact of reduced SNP densities on prediction accuracy was tested by varying minor allele frequency (MAF) thresholds, with no drop in prediction accuracy until the MAF threshold is set <0.3 (2,744 SNPs). These results point to the potential for GS to improve economically important traits in common carp breeding programs.

High-resolution mapping of the recombination landscape of the phytopathogen Fusarium graminearum suggests two-speed genome evolution.

Recombination is a major evolutionary force, increasing genetic diversity and permitting efficient coevolution of fungal pathogen(s) with their host(s). The ascomycete Fusarium graminearum is a devastating pathogen of cereal crops, and can contaminate food and feed with harmful mycotoxins. Previous studies have suggested a high adaptive potential of this pathogen, illustrated by an increase in pathogenicity and resistance to fungicides. In this study, we provide the first detailed picture of the crossover events occurring during meiosis and discuss the role of recombination in pathogen evolution. An experimental recombinant population (n = 88) was created and genotyped using 1306 polymorphic markers obtained from restriction site-associated DNA sequencing (RAD-seq) and aligned to the reference genome. The construction of a high-density linkage map, anchoring 99% of the total length of the reference genome, allowed the identification of 1451 putative crossovers, positioned at a median resolution of 24 kb. The majority of crossovers (87.2%) occurred in a relatively small portion of the genome (30%). All chromosomes demonstrated recombination-active sections, which had a near 15-fold higher crossover rate than non-active recombinant sections. The recombination rate showed a strong positive correlation with nucleotide diversity, and recombination-active regions were enriched for genes with a putative role in host-pathogen interaction, as well as putative diversifying genes. Our results confirm the preliminary analysis observed in other F. graminearum strains and suggest a conserved 'two-speed' recombination landscape. The consequences with regard to the evolutionary potential of this major fungal pathogen are also discussed.

Gene-centromere mapping in meiotic gynogenetic European seabass.

BACKGROUND: Fully isogenic lines in fish can be developed using "mitotic" gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families. RESULTS: From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly. CONCLUSION: Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.