Parentage assignment with genomic markers: a major advance for understanding and exploiting genetic variation of quantitative traits in farmed aquatic animals

APIS: an updated parentage assignment software managing triploids induced from diploid parents

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields, and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximizes the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96-192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1,232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low false discovery rate (<0.01), and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

PMSeeker: A Scheme Based on the Greedy Algorithm and the Exhaustive Algorithm to Screen Low-Redundancy Marker Sets for Large-Scale Parentage Assignment with Full Parental Genotyping

Parentage assignment is a genetic test that utilizes genetic characteristics, such as molecular markers, to identify the parental relationships within populations, which, in commercial fish farming, are almost always large and where full information on potential parents is known. To accurately find the true parents, the genotypes of all loci in the parentage marker set (PMS) are required for each individual being tested. With the same accuracy, a PMS containing a smaller number of markers will undoubtedly save experimental costs. Thus, this study established a scheme to screen low-redundancy PMSs using the exhaustive algorithm and greedy algorithm. When screening PMSs, the greedy algorithm selects markers based on the parental dispersity index (PDI), a uniquely defined metric that outperforms the probability of exclusion (PE). With the conjunctive use of the two algorithms, non-redundant PMSs were found for more than 99.7% of solvable cases in three groups of random sample experiments in this study. Then, a low-redundancy PMS can be composed using two or more of these non-redundant PMSs. This scheme effectively reduces the number of markers in PMSs, thus conserving human and experimental resources and laying the groundwork for the widespread implementation of parentage assignment technology in economic species breeding.

Potential of low-density genotype imputation for cost-efficient genomic selection for resistance to Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Flavobacterium columnare is the pathogen agent of columnaris disease, a major emerging disease that affects rainbow trout aquaculture. Selective breeding using genomic selection has potential to achieve cumulative improvement of the host resistance. However, genomic selection is expensive partly because of the cost of genotyping large numbers of animals using high-density single nucleotide polymorphism (SNP) arrays. The objective of this study was to assess the efficiency of genomic selection for resistance to F. columnare using in silico low-density (LD) panels combined with imputation. After a natural outbreak of columnaris disease, 2874 challenged fish and 469 fish from the parental generation (n = 81 parents) were genotyped with 27,907 SNPs. The efficiency of genomic prediction using LD panels was assessed for 10 panels of different densities, which were created in silico using two sampling methods, random and equally spaced. All LD panels were also imputed to the full 28K HD panel using the parental generation as the reference population, and genomic predictions were re-evaluated. The potential of prioritizing SNPs that are associated with resistance to F. columnare was also tested for the six lower-density panels.

RESULTS

The accuracies of both imputation and genomic predictions were similar with random and equally-spaced sampling of SNPs. Using LD panels of at least 3000 SNPs or lower-density panels (as low as 300 SNPs) combined with imputation resulted in accuracies that were comparable to those of the 28K HD panel and were 11% higher than the pedigree-based predictions.

CONCLUSIONS

Compared to using the commercial HD panel, LD panels combined with imputation may provide a more affordable approach to genomic prediction of breeding values, which supports a more widespread adoption of genomic selection in aquaculture breeding programmes.

Characterization of the genetic pool of the Canadienne dairy cattle breed

BACKGROUND

Canadienne cattle are the oldest breed of dairy cattle in North America. The Canadienne breed originates from cattle that were brought to America by the mid-seventeenth century French settlers. The herd book was established in 1886 and the current breed characteristics include dark coat color, small size compared to the modern Holstein breed, and overall rusticity shaped by the harsh environmental conditions that were prevalent during the settlement of North America. The Canadienne breed is an invaluable genetic resource due to its high resilience, longevity and fertility. However, it is heavily threatened with a current herd limited to an estimated 1200 registered animals, of which less than 300 are fullblood. To date, no effort has been made to document the genetic pool of this heritage breed in order to preserve it.

RESULTS

In this project, we used genomic data, which allow a precise description of the genetic makeup of a population, to provide valuable information on the genetic diversity of this heritage breed and suggest management options for its long-term viability. Using a panel that includes 640,000 single nucleotide polymorphisms (SNPs), we genotyped 190 animals grouped into six purity ranges. Unsupervised clustering analyses revealed three genetically distinct groups among those with the higher levels of purity. The observed heterozygosity was higher than expected even in the 100% purebreds. Comparison with Holstein genotypes showed significantly shorter runs of homozygosity for the Canadienne breed, which was unexpected due to the high inbreeding value calculated from pedigree data.

CONCLUSIONS

Overall, our data indicate that the fullblood gene pool of the Canadienne breed is more diversified than expected and that bloodline management could promote breed sustainability. In its current state, the Canadienne is not a dead-end breed but remains highly vulnerable due to its small population size.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus)

BACKGROUND

Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species.

RESULTS

From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip.

CONCLUSION

This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Animal board invited review: Widespread adoption of genetic technologies is key to sustainable expansion of global aquaculture

•

The extent of application of genetic technologies to aquaculture production varies widely by species and geography.

•

Achieving a more universal application of seed derived from scientifically based breeding programmes is an important goal in order to meet increasing global demands for seafood production.

•

This article reviews the status of genetic technologies across the world’s top 10 highly produced species.

•

Opportunities and barriers to achieving broad-scale uptake of genetic technologies in global aquaculture are discussed.

•

A future outlook for potential disruptive genetic technologies and how they might affect global aquaculture production is given.

Aquaculture production comprises a diverse range of species, geographies, and farming systems. The application of genetics and breeding technologies towards improved production is highly variable, ranging from the use of wild-sourced seed through to advanced family breeding programmes augmented by genomic techniques. This technical variation exists across some of the most highly produced species globally, with several of the top ten global species by volume generally lacking well-managed breeding programmes. Given the well-documented incremental and cumulative benefits of genetic improvement on production, this is a major missed opportunity. This short review focusses on (i) the status of application of selective breeding in the world’s most produced aquaculture species, (ii) the range of genetic technologies available and the opportunities they present, and (iii) a future outlook towards realising the potential contribution of genetic technologies to aquaculture sustainability and global food security.

What Can Genetics Do for the Control of Infectious Diseases in Aquaculture?

Simple Summary

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. This state-of-the-art review describes the application of genetics and genomics as novel tools to control infectious disease in aquaculture.

Abstract

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. An innovative approach to mitigate their impact on production is breeding for disease resistance: selection for domestication, family-based selection, marker-assisted selection, and more recently, genomic selection. Advances in genetics and genomics approaches to the control of infectious diseases are key to increasing aquaculture efficiency, profitability, and sustainability and to reducing its environmental footprint. Interaction and co-evolution between a host and pathogen can, however, turn breeding to boost infectious disease resistance into a potential driver of pathogenic change. Parallel molecular characterization of the pathogen and its virulence and antimicrobial resistance genes is therefore essential to understand pathogen evolution over time in response to host immunity, and to apply appropriate mitigation strategies.

Genomic Selection in Aquaculture Species

To date, genomic prediction has been conducted in about 20 aquaculture species, with a preference for intra-family genomic selection (GS). For every trait under GS, the increase in accuracy obtained by genomic estimated breeding values instead of classical pedigree-based estimation of breeding values is very important in aquaculture species ranging from 15% to 89% for growth traits, and from 0% to 567% for disease resistance. Although the implementation of GS in aquaculture is of little additional investment in breeding programs already implementing sib testing on pedigree, the deployment of GS remains sparse, but could be boosted by adaptation of cost-effective imputation from low-density panels. Moreover, GS could help to anticipate the effect of climate change by improving sustainability-related traits such as production yield (e.g., carcass or fillet yields), feed efficiency or disease resistance, and by improving resistance to environmental variation (tolerance to temperature or salinity variation). This chapter synthesized the literature in applications of GS in finfish, crustaceans and molluscs aquaculture in the present and future breeding programs.

Linking Pedigree Information to the Gene Expression Phenotype to Understand Differential Family Survival Mechanisms in Highly Fecund Fish: A Case Study in the Larviculture of Pacific Bluefin Tuna

Mass spawning in fish culture often brings about a marked variance in family size, which can cause a reduction in effective population sizes in seed production for stock enhancement. This study reports an example of combined pedigree information and gene expression phenotypes to understand differential family survival mechanisms in early stages of Pacific bluefin tuna, Thunnus orientalis, in a mass culture tank. Initially, parentage was determined using the partial mitochondrial DNA control region sequence and 11 microsatellite loci at 1, 10, 15, and 40 days post-hatch (DPH). A dramatic proportional change in the families was observed at around 15 DPH; therefore, transcriptome analysis was conducted for the 15 DPH larvae using a previously developed oligonucleotide microarray. This analysis successfully addressed the family-specific gene expression phenotypes with 5739 differentially expressed genes and highlighted the importance of expression levels of gastric-function-related genes at the developmental stage for subsequent survival. This strategy demonstrated herein can be broadly applicable to species of interest in aquaculture to comprehend the molecular mechanism of parental effects on offspring survival, which will contribute to the optimization of breeding technologies.

Optimization of Genomic Selection to Improve Disease Resistance in Two Marine Fishes, the European Sea Bass (Dicentrarchus labrax) and the Gilthead Sea Bream (Sparus aurata)

Disease outbreaks are a major threat to the aquaculture industry, and can be controlled by selective breeding. With the development of high-throughput genotyping technologies, genomic selection may become accessible even in minor species. Training population size and marker density are among the main drivers of the prediction accuracy, which both have a high impact on the cost of genomic selection. In this study, we assessed the impact of training population size as well as marker density on the prediction accuracy of disease resistance traits in European sea bass (Dicentrarchus labrax) and gilthead sea bream (Sparus aurata). We performed a challenge to nervous necrosis virus (NNV) in two sea bass cohorts, a challenge to Vibrio harveyi in one sea bass cohort and a challenge to Photobacterium damselae subsp. piscicida in one sea bream cohort. Challenged individuals were genotyped on 57K-60K SNP chips. Markers were sampled to design virtual SNP chips of 1K, 3K, 6K, and 10K markers. Similarly, challenged individuals were randomly sampled to vary training population size from 50 to 800 individuals. The accuracy of genomic-based (GBLUP model) and pedigree-based estimated breeding values (EBV) (PBLUP model) was computed for each training population size using Monte-Carlo cross-validation. Genomic-based breeding values were also computed using the virtual chips to study the effect of marker density. For resistance to Viral Nervous Necrosis (VNN), as one major QTL was detected, the opportunity of marker-assisted selection was investigated by adding a QTL effect in both genomic and pedigree prediction models. As training population size increased, accuracy increased to reach values in range of 0.51-0.65 for full density chips. The accuracy could still increase with more individuals in the training population as the accuracy plateau was not reached. When using only the 6K density chip, accuracy reached at least 90% of that obtained with the full density chip. Adding the QTL effect increased the accuracy of the PBLUP model to values higher than the GBLUP model without the QTL effect. This work sets a framework for the practical implementation of genomic selection to improve the resistance to major diseases in European sea bass and gilthead sea bream.

Social network analysis of the genealogy of strawberry: retracing the wild roots of heirloom and modern cultivars

The widely recounted story of the origin of cultivated strawberry (Fragaria × ananassa) oversimplifies the complex interspecific hybrid ancestry of the highly admixed populations from which heirloom and modern cultivars have emerged. To develop deeper insights into the three-century-long domestication history of strawberry, we reconstructed the genealogy as deeply as possible—pedigree records were assembled for 8,851 individuals, including 2,656 cultivars developed since 1775. The parents of individuals with unverified or missing pedigree records were accurately identified by applying an exclusion analysis to array-genotyped single-nucleotide polymorphisms. We identified 187 wild octoploid and 1,171 F. × ananassa founders in the genealogy, from the earliest hybrids to modern cultivars. The pedigree networks for cultivated strawberry are exceedingly complex labyrinths of ancestral interconnections formed by diverse hybrid ancestry, directional selection, migration, admixture, bottlenecks, overlapping generations, and recurrent hybridization with common ancestors that have unequally contributed allelic diversity to heirloom and modern cultivars. Fifteen to 333 ancestors were predicted to have transmitted 90% of the alleles found in country-, region-, and continent-specific populations. Using parent–offspring edges in the global pedigree network, we found that selection cycle lengths over the past 200 years of breeding have been extraordinarily long (16.0-16.9 years/generation), but decreased to a present-day range of 6.0-10.0 years/generation. Our analyses uncovered conspicuous differences in the ancestry and structure of North American and European populations, and shed light on forces that have shaped phenotypic diversity in F. × ananassa.

Maximum likelihood parentage assignment using quantitative genotypes

The cost of parentage assignment precludes its application in many selective breeding programmes and molecular ecology studies, and/or limits the circumstances or number of individuals to which it is applied. Pooling samples from more than one individual, and using appropriate genetic markers and algorithms to determine parental contributions to pools, is one means of reducing the cost of parentage assignment. This paper describes and validates a novel maximum likelihood (ML) parentage-assignment method, that can be used to accurately assign parentage to pooled samples of multiple individuals-previously published ML methods are applicable to samples of single individuals only-using low-density single nucleotide polymorphism (SNP) 'quantitative' (also referred to as 'continuous') genotype data. It is demonstrated with simulated data that, when applied to pools, this 'quantitative maximum likelihood' method assigns parentage with greater accuracy than established maximum likelihood parentage-assignment approaches, which rely on accurate discrete genotype calls; exclusion methods; and estimating parental contributions to pools by solving the weighted least squares problem. Quantitative maximum likelihood can be applied to pools generated using either a 'pooling-for-individual-parentage-assignment' approach, whereby each individual in a pool is tagged or traceable and from a known and mutually exclusive set of possible parents; or a 'pooling-by-phenotype' approach, whereby individuals of the same, or similar, phenotype/s are pooled. Although computationally intensive when applied to large pools, quantitative maximum likelihood has the potential to substantially reduce the cost of parentage assignment, even if applied to pools comprised of few individuals.

Multi-year pair-bonding in Murray cod (Maccullochella peelii)

Mating strategies in fishes are known to include polygyny, polyandry and monogamy and provide valuable insights regarding powerful evolutionary forces such as sexual selection. Monogamy is a complex of mating systems that has been relatively neglected. Previous work on mating strategies in fishes has often been based on observation and focused on marine species rather than freshwater fishes. SNPs are increasingly being used as a molecular ecology tool in non-model organisms, and methods of probabilistic genetic analysis of such datasets are becoming available for use in the absence of parental genotypes. This approach can be used to infer mating strategies. The long-term pair bonding seen in mammals, reptiles and birds has not been recorded in freshwater fishes-in every other respect an extremely diverse group. This study shows that multi-year pair bonding occurs in an Australian Percichthyid fish that exhibits paternal care of eggs and larvae. Using SNPs, full sibling pairs of larvae were found over multiple years in a three-year study. Stable isotope signatures of the larvae support the genetic inference that full sibling pairs shared a common mother, the ultimate source of that isotopic signature during oogenesis. Spatial and temporal clustering also suggests that the full sibling larvae are unlikely to be false positive identifications of the probabilistic identification of siblings. For the first time, we show multi-year pair bonding in a wild freshwater fish. This will have important conservation and management implications for the species. This approach could provide insights into many behavioural, ecological and evolutionary questions, particularly if this is not a unique case. Our findings are likely to initiate interest in seeking more examples of monogamy and alternative mating strategies in freshwater fishes, particularly if others improve methods of analysis of SNP data for identification of siblings in the absence of parental genotypes.

Genetic architecture and genomic selection of female reproduction traits in rainbow trout

BACKGROUND

Rainbow trout is a significant fish farming species under temperate climates. Female reproduction traits play an important role in the economy of breeding companies with the sale of fertilized eggs. The objectives of this study are threefold: to estimate the genetic parameters of female reproduction traits, to determine the genetic architecture of these traits by the identification of quantitative trait loci (QTL), and to assess the expected efficiency of a pedigree-based selection (BLUP) or genomic selection for these traits.

RESULTS

A pedigreed population of 1343 trout were genotyped for 57,000 SNP markers and phenotyped for seven traits at 2 years of age: spawning date, female body weight before and after spawning, the spawn weight and the egg number of the spawn, the egg average weight and average diameter. Genetic parameters were estimated in multi-trait linear animal models. Heritability estimates were moderate, varying from 0.27 to 0.44. The female body weight was not genetically correlated to any of the reproduction traits. Spawn weight showed strong and favourable genetic correlation with the number of eggs in the spawn and individual egg size traits, but the egg number was uncorrelated to the egg size traits. The genome-wide association studies showed that all traits were very polygenic since less than 10% of the genetic variance was explained by the cumulative effects of the QTLs: for any trait, only 2 to 4 QTLs were detected that explained in-between 1 and 3% of the genetic variance. Genomic selection based on a reference population of only one thousand individuals related to candidates would improve the efficiency of BLUP selection from 16 to 37% depending on traits.

CONCLUSIONS

Our genetic parameter estimates made unlikely the hypothesis that selection for growth could induce any indirect improvement for female reproduction traits. It is thus important to consider direct selection for spawn weight for improving egg production traits in rainbow trout breeding programs. Due to the low proportion of genetic variance explained by the few QTLs detected for each reproduction traits, marker assisted selection cannot be effective. However genomic selection would allow significant gains of accuracy compared to pedigree-based selection.

Development and validation of a RAD-Seq target-capture based genotyping assay for routine application in advanced black tiger shrimp (Penaeus monodon) breeding programs

Background

The development of genome-wide genotyping resources has provided terrestrial livestock and crop industries with the unique ability to accurately assess genomic relationships between individuals, uncover the genetic architecture of commercial traits, as well as identify superior individuals for selection based on their specific genetic profile. Utilising recent advancements in de-novo genome-wide genotyping technologies, it is now possible to provide aquaculture industries with these same important genotyping resources, even in the absence of existing genome assemblies. Here, we present the development of a genome-wide SNP assay for the Black Tiger shrimp (Penaeus monodon) through utilisation of a reduced-representation whole-genome genotyping approach (DArTseq).

Results

Based on a single reduced-representation library, 31,262 polymorphic SNPs were identified across 650 individuals obtained from Australian wild stocks and commercial aquaculture populations. After filtering to remove SNPs with low read depth, low MAF, low call rate, deviation from HWE, and non-Mendelian inheritance, 7542 high-quality SNPs were retained. From these, 4236 high-quality genome-wide loci were selected for baits-probe development and 4194 SNPs were included within a finalized target-capture genotype-by-sequence assay (DArTcap). This assay was designed for routine and cost effective commercial application in large scale breeding programs, and demonstrates higher confidence in genotype calls through increased call rate (from 80.2 ± 14.7 to 93.0% ± 3.5%), increased read depth (from 20.4 ± 15.6 to 80.0 ± 88.7), as well as a 3-fold reduction in cost over traditional genotype-by-sequencing approaches.

Conclusion

Importantly, this assay equips the P. monodon industry with the ability to simultaneously assign parentage of communally reared animals, undertake genomic relationship analysis, manage mate pairings between cryptic family lines, as well as undertake advance studies of genome and trait architecture. Critically this assay can be cost effectively applied as P. monodon breeding programs transition to undertaking genomic selection.

Harnessing genomics to fast-track genetic improvement in aquaculture

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

Potential for Genetic Improvement of Resistance to Perkinsus olseni in the Manila Clam, Ruditapes philippinarum, Using DNA Parentage Assignment and Mass Spawning

The Manila clam Ruditapes philippinarum, a major cultured shellfish species, is threatened by infection with the microparasite Perkinsus olseni, whose prevalence increases with high water temperatures. Under the current trend of climate change, the already severe effects of this parasitic infection might rapidly increase the frequency of mass mortality events. Treating infectious diseases in bivalves is notoriously problematic, therefore selective breeding for resistance represents a key strategy for mitigating the negative impact of pathogens. A crucial step in initiating selective breeding is the estimation of genetic parameters for traits of interest, which relies on the ability to record parentage and accurate phenotypes in a large number of individuals. Here, to estimate the heritability of resistance against P. olseni, a field experiment mirroring conditions in industrial clam production was set up, a genomic tool was developed for parentage assignment, and parasite load was determined through quantitative PCR. A mixed-family cohort of potentially 1,479 clam families was produced in a hatchery by mass spawning of 53 dams and 57 sires. The progenies were seeded in a commercial clam production area in the Venice lagoon, Italy, where high prevalence of P. olseni had previously been reported. Growth and parasite load were monitored every month and, after 1 year, more than 1,000 individuals were collected for DNA samples and phenotype recording. A pooled sequencing approach was carried out using DNA samples from the hatchery broodstock and from a Venice lagoon clam population, providing candidate markers used to develop a 245-SNP panel. Parentage assignment for 246 F1 individuals showed sire and dam representation were high (75 and 85%, respectively), indicating a very limited risk of inbreeding. Moderate heritability (0.23 ± 0.11-0.35 ± 0.13) was estimated for growth traits (shell length, shell weight, total weight), while parasite load showed high heritability, estimated at 0.51 ± 0.20. No significant genetic correlations were found between growth-associated traits and parasite load. Overall, the preliminary results provided by this study show high potential for selecting clams resistant to parasite load. Breeding for resistance may help limit the negative effects of climate change on clam production, as the prevalence of the parasite is predicted to increase under a future scenario of higher temperatures. Finally, the limited genetic correlation between resistance and growth suggests that breeding programs could incorporate dual selection without negative interactions.

First Evidence of Realized Selection Response on Fillet Yield in Rainbow Trout Oncorhynchus mykiss, Using Sib Selection or Based on Correlated Ultrasound Measurements

Fillet yield, the proportion of edible fillet relative to body weight, is a major trait to improve in fish sold processed, as it has a direct impact on profitability and can simultaneously decrease the environmental impact of producing a given amount of fillet. However, it is difficult to improve by selective breeding, because it cannot be measured on live breeding candidates, its phenotypic variation is low, and, as a ratio, it is not normally distributed and a same change in fillet yield can be the result of different changes in fillet weight and body weight. Residual headless gutted carcass weight (rHGCW) is heritable and highly genetically correlated to Fillet% in rainbow trout, and can be predicted by the ratio of abdominal wall thickness to depth of the peritoneal cavity (E8/E23), measured on live fish by ultrasound tomography. We selected broodstock based on rHGCW, measured on sibs of the selection candidates, on ultrasound measurements (E8/E23) measured on the selection candidates, or a combination of both. Seven broodstock groups were selected: fish with 15% highest (rHGCW+) or lowest (rHGCW−) EBV for rHGCW, with 15% highest (E8/E23+) or lowest (E8/E23−) EBV for E8/E23, with both rHGCW+ and E8/E23+ (Both+) or rHGCW− and E8/E23− (Both−), or with close to zero EBVs for both traits (Mid). Seven corresponding groups of offspring were produced and reared communally. At harvest size (1.5 kg mean weight), 1,561 trout were slaughtered, measured for the traits of interest, and pedigreed with DNA fingerprinting. Offspring from groups Both+, rHGCW+ and E8/E23+ had a higher EBV for rHGCW than the control group, while down-selected groups had a lower EBV. Looking at the phenotypic mean for Fillet% (correlated response), up-selected fish had more fillet than down-selected fish. The highest difference was between Both+ (69.36%) and Both− (68.20%), a 1.16% units difference in fillet percentage. The change in Fillet% was explained by an opposite change in Viscera%, while Head% remained stable. Selection using sib information on rHGCW was on average more efficient than selection using the candidates’ own E8/E23 phenotypes, and downward selection (decreasing Fillet%) was more efficient than upward selection.

Considerations for Maintaining Family Diversity in Commercially Mass-Spawned Penaeid Shrimp: A Case Study on Penaeus monodon

Skewed family distributions are common in aquaculture species that are highly fecund, communally (mass) spawned, and/or communally reared. The magnitude of skews pose challenges for maintaining family-specific genetic diversity, as increased resources are required to detect individuals from underrepresented families, or reliably determine relative survival as a measure of family performance. There is limited understanding of family skews or changes in family proportion of communally reared shrimp under commercial rearing conditions and particularly how this may affect genotyping strategies to recover family performance data in breeding programs. In this study, three separate batches of shrimp, Penaeus monodon, were communally spawned and reared, and then sampled as larvae when ponds were stocked at 30 days of culture (DOC) and as juveniles from commercial ponds during harvest at 150 DOC. A total of 199 broodstock contributed to the 5,734 progeny that were genotyped with a custom multiplex single nucleotide polymorphism (SNP) panel, and family assignments were cross-referenced using two parentage assignment methods, CERVUS and COLONY. A total of 121 families were detected, with some families contributing up to 11% of progeny at 30 DOC and up to 18% of progeny at harvest. Significant changes were detected for 20% of families from 30 to 150 DOC, with up to a 9% change in relative contribution. Family skew data was applied in several models to determine the optimal sample size to detect families, along with the ability to detect changes in relative family contribution over time. Results showed that an order of magnitude increase in sampling was required to capture the lowest represented 25% of families, as well as significantly improve the accuracy to determine changes in family proportion from 30 to 150 DOC. Practical measures may be implemented at the hatchery to reduce family skews; a cost-effective measure may be to address the initial magnitude differences in viable progeny produced among families, by pooling equal quantities of hatched larvae from each family. This study demonstrates the relationships between skews in families under commercial conditions, the ability to accurately detect families, and the balance of sampling effort and genotyping cost in highly fecund species such as shrimp.

Hatchery-reared enhancement program for silver carp (Hypophthalmichthys molitrix) in the middle Yangtze River: monitoring the effectiveness based on parentage analysis

Introduction

A hatchery-reared silver carp (Hypophthalmichthys molitrix) program has been intensively carried out since 2010 to enhance the rapidly declining fisheries production in the middle Yangtze River. However, only a little information regarding the effectiveness of the enhancement program has been reported. In this context, this study investigates on an enhancement program through monitoring the efficacy based on parentage analysis.

Methods

A total of 1,529 hatchery-reared fish and 869 larvae were sampled from the middle Yangtze River in 2016 and 2017 and were genotyped by thirteen microsatellite loci. Based on the results of parentage analysis the larvae were divided into three populations: (1) larvae population with both parents being hatchery-reared fish (=R), (2) larvae population with only a male or a female parent being hatchery-reared fish (=H), and (3) larvae population with no hatchery-reared fish parent (=W). The following analyses were also carried out: (1) assessing the contribution of hatchery-reared offspring to larval resources, and (2) evaluating the genetic effect of stock enhancement on the wild population.

Results

In total, 10.37% and 11.56% of larvae were identified as the offspring produced by hatchery-reared fish released in 2016 and 2017, respectively. In 2017, some of the larvae were assigned unambiguously to hatchery-reared fish released in 2016. In terms of the number of offspring produced, the hatchery-reared fish have shown significant variations. No significant differences were found among all the larvae populations concerning genetic parameters for diversity. High levels of genetic diversity of all larvae populations were obtained. Low F _STvalues obtained from pairwise F _ST analysis, as well as the analysis of molecular variance (AMOVA), revealed high genetic structural similarity among all the larvae populations. The genetic composition of the W larvae population in 2017 was different from that of all other larvae populations (all larvae populations in 2016, and R and H larvae populations in 2017), as demonstrated from the results of STRUCTURE and PCA analyses.

Conclusion

It was demonstrated that hatchery-reared fish are successful in producing the offspring in the natural environment during multiple years, which might assist in increasing the abundance of larvae. The hatchery-reared fish had variations in terms of the success rates on reproduction. Also, the hatchery-reared enhancement program had no significant effect on the genetic diversity or the genetic structure of wild populations. However, the genetic component of the W larvae population in 2017 was changed as compared to 2016, which was not due to the hatchery-reared enhancement program for silver carp. This could be due to flooding, but the specific causes need further studies. Our results clearly show the necessity to continuously inspect the genetic impact of the enhancement program so that historical information can be utilized for further research.

Tracing individuals and populations of the tambaqui, Colossoma macropomum (Cuvier, 1818), from Brazilian hatcheries using microsatellite markers

BACKGROUND

In recent years, tracing of alimentary produce of animal origin has become increasingly important, for economic, food safety and ecological reasons. The tambaqui, Colossoma macropomum, is the native fish most farmed in Brazil. The reliable identification of the origin of tambaquis (wild or farmed) offered for sale to the general public has become necessary to satisfy regulatory norms and uphold consumer confidence. Molecular methods based on the analysis of DNA sequences have often been used to evaluate the potential for tracing farmed fish, given their reliability and precision.

RESULTS

Full likelihood and Bayesian approaches proved to be the most efficient for the identification, respectively, of individuals and populations for most of the fish sampled from seven hatcheries and one wild stock. The exclusion method and genetic distances were the least effective approaches for the identification of individuals and populations. The Bayesian method identified correctly more than 99% of the fry from most stocks, except those of the Santarém hatchery and River Amazon wild stock, which presented the best results for individual identification.

CONCLUSIONS

The identification of populations was effective for most hatcheries, although the identification of individuals from most stocks was hampered by the reduced genetic variability. © 2018 Society of Chemical Industry.

The European sea bass: a key marine fish model in the wild and in aquaculture

The European sea bass (Dicentrarchus labrax L.) is a marine fish of key economic and cultural importance in Europe. It is now more an aquaculture than a fisheries species (>96% of the production in 2016), although modern rearing techniques date back only from the late 1980s. It also has high interest for evolutionary studies, as it is composed of two semispecies (Atlantic and Mediterranean lineages) that have come into secondary contact following the last glaciation. Based on quantitative genetics studies of most traits of interest over the past 10–15 years, selective breeding programs are now applied to this species, which is at the beginning of its domestication process. The availability of a good quality reference genome has accelerated the development of new genomic resources, including SNP arrays that will enable genomic selection to improve genetic gain. There is a need to improve feed efficiency, both for economic and environmental reasons, but this will require novel phenotyping approaches. Further developments will likely focus on the understanding of genotype‐by‐environment interactions, which will be important both for efficient breeding of farmed stocks and for improving knowledge of the evolution of natural populations. At the interface between both, the domestication process must be better understood to improve production and also to fully evaluate the possible impact of aquaculture escapees on wild populations. The latter is an important question for all large‐scale aquaculture productions.

Detailed insights into pan-European population structure and inbreeding in wild and hatchery Pacific oysters (Crassostrea gigas) revealed by genome-wide SNP data

Cultivated bivalves are important not only because of their economic value, but also due to their impacts on natural ecosystems. The Pacific oyster (Crassostrea gigas) is the world's most heavily cultivated shellfish species and has been introduced to all continents except Antarctica for aquaculture. We therefore used a medium-density single nucleotide polymorphism (SNP) array to investigate the genetic structure of this species in Europe, where it was introduced during the 1960s and has since become a prolific invader of coastal ecosystems across the continent. We analyzed 21,499 polymorphic SNPs in 232 individuals from 23 localities spanning a latitudinal cline from Portugal to Norway and including the source populations of Japan and Canada. We confirmed the results of previous studies by finding clear support for a southern and a northern group, with the former being indistinguishable from the source populations indicating the absence of a pronounced founder effect. We furthermore conducted a large-scale comparison of oysters sampled from the wild and from hatcheries to reveal substantial genetic differences including significantly higher levels of inbreeding in some but not all of the sampled hatchery cohorts. These findings were confirmed by a smaller but representative SNP dataset generated using restriction site-associated DNA sequencing. We therefore conclude that genomic approaches can generate increasingly detailed insights into the genetics of wild and hatchery produced Pacific oysters.

Genomic Selection in Aquaculture: Application, Limitations and Opportunities With Special Reference to Marine Shrimp and Pearl Oysters

Within aquaculture industries, selection based on genomic information (genomic selection) has the profound potential to change genetic improvement programs and production systems. Genomic selection exploits the use of realized genomic relationships among individuals and information from genome-wide markers in close linkage disequilibrium with genes of biological and economic importance. We discuss the technical advances, practical requirements, and commercial applications that have made genomic selection feasible in a range of aquaculture industries, with a particular focus on molluscs (pearl oysters, Pinctada maxima) and marine shrimp (Litopenaeus vannamei and Penaeus monodon). The use of low-cost genome sequencing has enabled cost-effective genotyping on a large scale and is of particular value for species without a reference genome or access to commercial genotyping arrays. We highlight the pitfalls and offer the solutions to the genotyping by sequencing approach and the building of appropriate genetic resources to undertake genomic selection from first-hand experience. We describe the potential to capture large-scale commercial phenotypes based on image analysis and artificial intelligence through machine learning, as inputs for calculation of genomic breeding values. The application of genomic selection over traditional aquatic breeding programs offers significant advantages through being able to accurately predict complex polygenic traits including disease resistance; increasing rates of genetic gain; minimizing inbreeding; and negating potential limiting effects of genotype by environment interactions. Further practical advantages of genomic selection through the use of large-scale communal mating and rearing systems are highlighted, as well as presenting rate-limiting steps that impact on attaining maximum benefits from adopting genomic selection. Genomic selection is now at the tipping point where commercial applications can be readily adopted and offer significant short- and long-term solutions to sustainable and profitable aquaculture industries.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

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.

Genomic Tools and Selective Breeding in Molluscs

The production of most farmed molluscs, including mussels, oysters, scallops, abalone, and clams, is heavily dependent on natural seed from the plankton. Closing the lifecycle of species in hatcheries can secure independence from wild stocks and enables long-term genetic improvement of broodstock through selective breeding. Genomic techniques have the potential to revolutionize hatchery-based selective breeding by improving our understanding of the characteristics of mollusc genetics that can pose a challenge for intensive aquaculture and by providing a new suite of tools for genetic improvement. Here we review characteristics of the life history and genetics of molluscs including high fecundity, self-fertilization, high genetic diversity, genetic load, high incidence of deleterious mutations and segregation distortion, and critically assess their impact on the design and effectiveness of selective breeding strategies. A survey of the results of current breeding programs in the literature show that selective breeding with inbreeding control is likely the best strategy for genetic improvement of most molluscs, and on average growth rate can be improved by 10% per generation and disease resistance by 15% per generation across the major farmed species by implementing individual or family-based selection. Rapid advances in sequencing technology have resulted in a wealth of genomic resources for key species with the potential to greatly improve hatchery-based selective breeding of molluscs. In this review, we catalog the range of genomic resources currently available for molluscs of aquaculture interest and discuss the bottlenecks, including lack of high-quality reference genomes and the relatively high cost of genotyping, as well as opportunities for applying genomics-based selection.

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.

Genomic relationships reveal significant dominance effects for growth in hybrid Eucalyptus

Non-additive genetic effects can be effectively exploited in control-pollinated families with the availability of genome-wide markers. We used 41,304 SNP markers and compared pedigree vs. marker-based genetic models by analysing height, diameter, basic density and pulp yield for Eucalyptus urophylla × E.grandis control-pollinated families represented by 949 informative individuals. We evaluated models accounting for additive, dominance, and first-order epistatic interactions (additive by additive, dominance by dominance, and additive by dominance). We showed that the models can capture a large proportion of the genetic variance from dominance and epistasis for growth traits as those components are typically not independent. We also showed that we could partition genetic variances more precisely when using relationship matrices derived from markers compared to using only pedigree information. In addition, phenotypic prediction accuracies were only slightly increased by including dominance effects for growth traits since estimates of non-additive variances yielded rather high standard errors. This novel result improves our current understanding of the architecture of quantitative traits and recommends accounting for dominance variance when developing genomic selection strategies in hybrid Eucalyptus.

Heritability of Boldness and Hypoxia Avoidance in European Seabass, Dicentrarchus labrax

To understand the genetic basis of coping style in European seabass, fish from a full factorial mating (10 females x 50 males) were reared in common garden and individually tagged. Individuals coping style was characterized through behavior tests at four different ages, categorizing fish into proactive or reactive: a hypoxia avoidance test (at 255 days post hatching, dph) and 3 risk-taking tests (at 276, 286 and 304 dph). We observed significant heritability of the coping style, higher for the average of risk-taking scores (h² = 0.45 ± 0.14) than for the hypoxia avoidance test (h² = 0.19 ± 0.10). The genetic correlations between the three risk-taking scores were very high (r_A = 0.96–0.99) showing that although their repeatability was moderately high (r_P = 0.64–0.72), successive risk-taking tests evaluated the same genetic variation. A mild genetic correlation between the results of the hypoxia avoidance test and the average of risk-taking scores (0.45 ± 0.27) suggested that hypoxia avoidance and risk-taking tests do not address exactly the same behavioral and physiological responses. Genetic correlations between weight and risk taking traits showed negative values whatever the test used in our population i.e. reactive individual weights were larger. The results of this quantitative genetic analysis suggest a potential for the development of selection programs based on coping styles that could increase seabass welfare without altering growth performances. Overall, it also contributes to a better understanding of the origin and the significance of individual behavioral differences.