Estimating breeding values with molecular relatedness and reconstructed pedigrees in natural mating populations of common sole, Solea solea

Genetic evaluations of dairy goats with few pedigree data: different approaches to use molecular information

One of the limitations of implementing animal breeding programs in small-scale or extensive production systems is the lack of production records and genealogical records. In this context, molecular markers could help to gain information for the breeding program. This study addresses the inclusion of molecular data into traditional genetic evaluation models as a random effect by molecular pedigree reconstruction and as a fixed effect by Bayesian clustering. The methods were tested for lactation curve traits in 14 dairy goat herds with incomplete phenotypic data and pedigree information. The results showed an increment of 37.3% of the relationships regarding the originals with MOLCOAN and clustering into five genetic groups. Data leads to estimating additive variance, error variance, and heritability with four different models, including pedigree and molecular information. Deviance Information Criterion (DIC) values demonstrate a greater fitting of the models that include molecular information either as fixed (genetic clusters) or as random (molecular matrix) effects. The molecular information of simple markers can complement genetic improvement strategies in populations with little information.

Impact of genotypic errors with equal and unequal family contribution on accuracy of genomic prediction in aquaculture using simulation

Genotypic errors, conflict between recorded genotype and the true genotype, can lead to false or biased population genetic parameters. Here, the effect of genotypic errors on accuracy of genomic predictions and genomic relationship matrix are investigated using a simulation study based on population and genomic structure comparable to black tiger prawn, Penaeus monodon. Fifty full-sib families across five generations with phenotypic and genotypic information on 53 K SNPs were simulated. Ten replicates of different scenarios with three heritability estimates, equal and unequal family contributions were generated. Within each scenario, four SNP densities and three genotypic error rates in each SNP density were implemented. Results showed that family contribution did not have a substantial impact on accuracy of predictions across different datasets. In the absence of genotypic errors, 3 K SNP density was found to be efficient in estimating the accuracy, whilst increasing the SNP density from 3 to 20 K resulted in a marginal increase in accuracy of genomic predictions using the current population and genomic parameters. In addition, results showed that the presence of even 10% errors in a 10 and 20 K SNP panel might not have a severe impact on accuracy of predictions. However, below 10 K marker density, even a 5% error can result in lower accuracy of predictions.

Genetic Estimates for Growth and Shape-Related Traits in the Flatfish Senegalese Sole

Simple Summary

To increase competitiveness, the aquaculture flatfish industry demands animals with optimal growth rates and a high shape quality. Genetic breeding is an essential tool to achieve these goals but it requires the estimation of the genetic components of these traits under industrial conditions. The current study provides phenotypic data and genetic parameters of eight traits related to growth and shape quality. The high heritabilities and correlations obtained support that genetic breeding programs can be successfully implemented in Senegalese sole to optimize production.

Abstract

Shape quality is very important in flatfish aquaculture due to the impact on commercialization. The Senegalese sole (Solea senegalensis) is a valuable flatfish with a highly elliptic body that slightly changes with age and size, and it is prone to accumulating malformations during the production cycle. The present study aims to investigate the genetic parameters of two growth traits (weight and standard length) and six shape quality predictors (ellipticity, three body heights (body height at the pectoral fin base [BHP], body maximum height [BMH] and caudal peduncle height [CPH]) and two ratios (BMH/BHP and BMH/CPH)). These traits were measured before the on-growing stage (age ~400 days (d)) and at harvest (~800 d). Phenotypic data, heritabilities and genetic and phenotypic correlations between the traits are presented and discussed. High or very high heritabilities (0.433–0.774) were found for growth traits, body heights and ellipticity and they were higher at 400 than 800 d. In contrast, the ratios of BMH/BHP and BMH/CPH were less heritable (0.144–0.306). Positive and very high (>0.95) correlations between growth traits and the three heights were found and decreased with age. In contrast, ellipticity had negative and medium-high genetic correlations with growth traits and heights, indicating fish selected for bigger size would also become rounder. The ratio of BMH/CPH showed low genetic correlations with all traits and provided complementary information to ellipticity for a better fitting to the expected lanceolate body morphology of sole. The genetic correlations for all traits at both ages were very high, indicating that selection before entering the growth-out stage in recirculation aquaculture systems is recommended to accelerate genetic gains.

Estimation of genetic parameters for growth and carcass traits in turbot (Scophthalmus maximus)

Information on phenotypic and genetic (co)variance for production traits in turbot is required to improve breeding programs. So far, information on morphometric growth traits is sparse and completely lacking on quality carcass traits like fillet weight or fillet yield for turbot. As part of a long-term study we explored the phenotypic and genetic (co)variance of 16 biometrical and carcass traits of three different European turbot strains. Fish were reared under commercial grow-out conditions, including size grading. We used molecular relatedness (MR) methods based on genotyping with 96 microsatellite markers and animal models. We included an adapted condition factor for Pleuronectiformes (FCIPLN) and average daily weight gain (ADG) between the ages of 300 and 500 d post-hatch (dph) for their potential correlation with body weight at harvest. Heritability estimates for all traits were low to medium (0.04–0.29) when strains were jointly analyzed. Separate analysis of strains yielded higher heritability estimates (0.12–0.43). Genetic correlations between weight-related traits were highly positive (0.70–0.99), while runs with yield and ratio traits often resulted in unreliable estimates of genetic correlation due to high standard errors. Body weight (h2=0.19), fillet yield (h2=0.15), and dressing percentage (h2=0.17) are particularly promising selection traits for turbot breeding.

'apparent': a simple and flexible R package for accurate SNP-based parentage analysis in the absence of guiding information

Background

The accurate determination of parent-progeny relationships within both in situ natural populations and ex situ genetic resource collections can greatly enhance plant breeding/domestication efforts and support plant genetic resource conservation strategies. Although a range of parentage analysis tools are available, none are designed to infer such relationships using genome-wide single nucleotide polymorphism (SNP) data in the complete absence of guiding information, such as generational groups, partial pedigrees, or genders. The R package (‘apparent’) developed and presented here addresses this gap.

Results

‘apparent’ adopts a novel strategy of parentage analysis based on a test of genetic identity between a theoretically expected progeny (EP_ij), whose genotypic state can be inferred at all homozygous loci for a pair of putative parents (i and j), and all potential offspring (PO_k), represented by the k individuals of a given germplasm collection. Using the Gower Dissimilarity metric (GD), genetic identity between EP_ij and PO_k is taken as evidence that individuals i and j are the true parents of offspring k. Significance of a given triad (parental pair_ij + offspring_k) is evaluated relative to the distribution of all GD_ij|k values for the population. With no guiding information provided, ‘apparent’ correctly identified the parental pairs of 15 lines of known pedigree within a test population of 77 accessions of Actinidia arguta, a performance unmatched by five other commonly used parentage analysis tools. In the case of an inconclusive triad analysis due to the absence of one parent from the test population, ‘apparent’ can perform a subsequent dyad analysis to identify a likely single parent for a given offspring. Average dyad analysis accuracy was 73.3% in the complete absence of pedigree information but increased to 100% when minimal generational information (adults vs. progeny) was provided.

Conclusions

The ‘apparent’ R package is a fast and accurate parentage analysis tool that uses genome-wide SNP data to identify parent-progeny relationships within populations for which no a priori knowledge of family structure exists.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2662-3) contains supplementary material, which is available to authorized users.

A genetic linkage map of sole (Solea solea): a tool for evolutionary and comparative analyses of exploited (flat)fishes

Linkage maps based on markers derived from genes are essential evolutionary tools for commercial marine fish to help identify genomic regions associated with complex traits and subject to selective forces at play during exploitation or selective breeding. Additionally, they allow the use of genomic information from other related species for which more detailed information is available. Sole (solea solea L.) is a commercially important flatfish species in the North Sea, subject to overexploitation and showing evidence of fisheries-induced evolutionary changes in growth- and maturation-related traits. Sole would definitely benefit from a linkage map to better understand how evolution has shaped its genome structure. This study presents a linkage map of sole based on 423 single nucleotide polymorphisms derived from expressed sequence tags and 8 neutral microsatellite markers. The total map length is 1233.8 cM and consists of 38 linkage groups with a size varying between 0 to 92.1 cM. Being derived from expressed sequence tags allowed us to align the map with the genome of four model fish species, namely medaka (Oryzias latipes), Nile tilapia (Oreochromis niloticus), three-spined stickleback (Gasterosteus aculeatus) and green spotted pufferfish (Tetraodon nigroviridis). This comparison revealed multiple conserved syntenic regions with all four species, and suggested that the linkage groups represent 21 putative sole chromosomes. The map was also compared to the linkage map of turbot (Scophthalmus maximus), another commercially important flatfish species and closely related to sole. For all putative sole chromosomes (except one) a turbot homolog was detected, confirming the even higher degree of synteny between these two flatfish species.

Fatty acid composition of eggs and its relationships to egg and larval viability from domesticated common sole (Solea solea) breeders

The study of lipids and fatty acids (FAs) has been used in the assessment of egg quality because their composition can influence the fertilization rate, hatching, survival and growth of marine fish larvae. For these reasons, the lipid content (TL) and fatty acid composition of common sole (Solea solea) eggs were measured and correlated to egg and larval viability parameters throughout an entire reproductive season. Seventeen batches of fertile eggs obtained from natural spawning of captive breeders were characterized for the TL, FA profile, hatching rate (HR) and survival rate of larvae (SR) at 0-6 days post-hatching (dph). The egg FA composition reflected the composition of the feed supplied to the broodstock during summer and autumn (before and during vitellogenesis) rather than that supplied during the spawning season. In general, the egg FA profile showed minimal differences among the early-, mid- and late-spawning periods (possibly due to the change of the diet and/or water temperature) indicating that it is possible to obtain a similar egg quality in terms of egg FA profile over 2 months of spawning. Saturated FAs and monounsaturated FAs (MUFA) were positively correlated with HR, while TL, 22 : 6n-3 (DHA), 20 : 4n-6 (ARA), polyunsaturated FAs of the (n-3) series (n-3 PUFA) and polyunsaturated FAs of the (n-6) series were negatively correlated (p ≤ 0.05). MUFA, 20 : 5n-3 (EPA), n-6/n-3 were positively correlated with SR, while DHA, n-3 PUFA, DHA/EPA were negatively correlated (p ≤ 0.05). In conclusion, the feed supplied before and during vitellogenesis has a major role in determining the egg FA profile in common sole. The relationships found between TL and FAs with egg and larval viability parameters differ from many other farmed marine fish species, which may suggest the need for a specific broodstock feed for this species.

Pedigree-free estimates of heritability in the wild: promising prospects for selfing populations

Estimating the genetic variance available for traits informs us about a population's ability to evolve in response to novel selective challenges. In selfing species, theory predicts a loss of genetic diversity that could lead to an evolutionary dead-end, but empirical support remains scarce. Genetic variability in a trait is estimated by correlating the phenotypic resemblance with the proportion of the genome that two relatives share identical by descent ('realized relatedness'). The latter is traditionally predicted from pedigrees (Φ A : expected value) but can also be estimated using molecular markers (average number of alleles shared). Nevertheless, evolutionary biologists, unlike animal breeders, remain cautious about using marker-based relatedness coefficients to study complex phenotypic traits in populations. In this paper, we review published results comparing five different pedigree-free methods and use simulations to test individual-based models (hereafter called animal models) using marker-based relatedness coefficients, with a special focus on the influence of mating systems. Our literature review confirms that Ritland's regression method is unreliable, but suggests that animal models with marker-based estimates of relatedness and genomic selection are promising and that more testing is required. Our simulations show that using molecular markers instead of pedigrees in animal models seriously worsens the estimation of heritability in outcrossing populations, unless a very large number of loci is available. In selfing populations the results are less biased. More generally, populations with high identity disequilibrium (consanguineous or bottlenecked populations) could be propitious for using marker-based animal models, but are also more likely to deviate from the standard assumptions of quantitative genetics models (non-additive variance).

Populus trichocarpa cell wall chemistry and ultrastructure trait variation, genetic control and genetic correlations

The increasing ecological and economical importance of Populus species and hybrids has stimulated research into the investigation of the natural variation of the species and the estimation of the extent of genetic control over its wood quality traits for traditional forestry activities as well as the emerging bioenergy sector. A realized kinship matrix based on informative, high-density, biallelic single nucleotide polymorphism (SNP) genetic markers was constructed to estimate trait variance components, heritabilities, and genetic and phenotypic correlations. Seventeen traits related to wood chemistry and ultrastructure were examined in 334 9-yr-old Populus trichocarpa grown in a common-garden plot representing populations spanning the latitudinal range 44° to 58.6°. In these individuals, 9342 SNPs that conformed to Hardy-Weinberg expectations were employed to assess the genomic pair-wise kinship to estimate narrow-sense heritabilities and genetic correlations among traits. The range-wide phenotypic variation in all traits was substantial and several trait heritabilities were > 0.6. In total, 61 significant genetic and phenotypic correlations and a network of highly interrelated traits were identified. The high trait variation, the evidence for moderate to high heritabilities and the identification of advantageous trait combinations of industrially important characteristics should aid in providing the foundation for the enhancement of poplar tree breeding strategies for modern industrial use.

DNA barcodes of Asian Houbara Bustard (Chlamydotis undulata macqueenii)

Populations of Houbara Bustards have dramatically declined in recent years. Captive breeding and reintroduction programs have had limited success in reviving population numbers and thus new technological solutions involving molecular methods are essential for the long term survival of this species. In this study, we sequenced the 694 bp segment of COI gene of the four specimens of Asian Houbara Bustard (Chlamydotis undulata macqueenii). We also compared these sequences with earlier published barcodes of 11 individuals comprising different families of the orders Gruiformes, Ciconiiformes, Podicipediformes and Crocodylia (out group). The pair-wise sequence comparison showed a total of 254 variable sites across all the 15 sequences from different taxa. Three of the four specimens of Houbara Bustard had an identical sequence of COI gene and one individual showed a single nucleotide difference (G > A transition at position 83). Within the bustard family (Otididae), comparison among the three species (Asian Houbara Bustard, Great Bustard (Otis tarda) and the Little Bustard (Tetrax tetrax)), representing three different genera, showed 116 variable sites. For another family (Rallidae), the intra-family variable sites among the individuals of four different genera were found to be 146. The COI genetic distances among the 15 individuals varied from 0.000 to 0.431. Phylogenetic analysis using 619 bp nucleotide segment of COI clearly discriminated all the species representing different genera, families and orders. All the four specimens of Houbara Bustard formed a single clade and are clearly separated from other two individuals of the same family (Otis tarda and Tetrax tetrax). The nucleotide sequence of partial segment of COI gene effectively discriminated the closely related species. This is the first study reporting the barcodes of Houbara Bustard and would be helpful in future molecular studies, particularly for the conservation of this threatened bird in Saudi Arabia.

Social and genetic interactions drive fitness variation in a free-living dolphin population

The evolutionary forces that drive fitness variation in species are of considerable interest. Despite this, the relative importance and interactions of genetic and social factors involved in the evolution of fitness traits in wild mammalian populations are largely unknown. To date, a few studies have demonstrated that fitness might be influenced by either social factors or genes in natural populations, but none have explored how the combined effect of social and genetic parameters might interact to influence fitness. Drawing from a long-term study of wild bottlenose dolphins in the eastern gulf of Shark Bay, Western Australia, we present a unique approach to understanding these interactions. Our study shows that female calving success depends on both genetic inheritance and social bonds. Moreover, we demonstrate that interactions between social and genetic factors also influence female fitness. Therefore, our study represents a major methodological advance, and provides critical insights into the interplay of genetic and social parameters of fitness.