Genetic Parameters and Genotype × Diet Interaction for Body Weight Performance and Fat in Gilthead Seabream.

There has been thorough research on the genotype by diet interaction and the extent of its impact on the genetic evaluation, using a partly replaced marine animal protein on the gilthead seabream. To do that, 8356 individuals were gathered from two batches and followed different diets: a high-plant-protein diet containing 85% plant proteins and a standard commercial one containing 30% marine animal protein. During the experiment, body weight, growth and fat content were recorded. High heritability estimates were detected for the body weight performance and fat content. A small effect of genotype by diet interaction was detected in all phenotypes (presented as the genetic correlations from 0.95 to 0.97) but a medium-high ranking correlation between the breeding values for each trait was estimated (0.72-0.70). A higher expected response to selection for the body weight performance was detected using the standard commercial rather than the plant-based diet. Based on the findings, the establishment of a plant-based diet breeding strategy can be achieved provided the reduction of the cost of aquafeed is attained, though a lower genetic gain is expected.

Genetic diversity and substructuring of the Hungarian merino sheep breed using microsatellite markers.

The Hungarian Merino sheep breed (Ovis aries) is the most significant animal resource of the Hungarian sheep sector which, unfortunately, has gone through a huge reduction in number during the last decades and became endangered in 2014. A modern molecular genetic survey is now becoming more than necessary in order to characterize the within-breed genetic diversity and structure. For that reason, six Hungarian Merino flocks were genotyped in 16 microsatellite markers. In total, 144 different alleles were found and the mean values of observed and expected heterozygosity were 0.714 and 0.705, respectively, suggesting a noticeable genetic variability of the breed. The genetic differentiation of the Hungarian flocks was generally low, as reflected by the estimated total FST value (0.036), the extended pattern of admixture in Structure analysis, as well as, by the noticeable level of genetic clustering in UPGMA and FCA analyses. However, two out of the six studied flocks tended to be genetically more distant. The outcome of our study could be a starting point for a planned breeding strategy of the Hungarian Merino breed, by keeping the within-flock genetic variability in priority, as well as, by preserving the potential genetic uniqueness with close monitoring of the inbreeding.

Genomic analysis of the rare British Lop pig and identification of distinctive genomic markers.

Concentration of production on a few commercial pig breeds has led to the marginalization of many native, numerically small breeds, increasing their risk of endangerment. In the UK, one such rare breed is the British Lop, a lop-eared breed, of similar origin to the Welsh breed. The objective of the present study was to address the genomic status of the British Lop and its relationship with other breeds and identify a small set of genomic markers that uniquely characterize and distinguish British Lop animals. Results have shown that the British Lop is a relatively distinct population with reduced genomic diversity and effective size consistent with its status as a rare breed. Furthermore, we demonstrated the genetic closeness of the British Lop to phenotypically similar breeds such as Landrace and Welsh as well Large White, Middle White and Pietrain. Finally, a set of 75 Single Nucleotide Polymorphisms distributed across multiple chromosomes were identified and validated as markers that can consistently distinguish British Lops from other closely related breeds. Results may inform breeding and management strategies aiming to enhance diversity as well as the development of a breed purity test.

Genomic Selection and Genome-Wide Association Analysis for Stress Response, Disease Resistance and Body Weight in European Seabass.

The majority of the genetic studies in aquaculture breeding programs focus on commercial traits such as body weight, morphology, and resistance against diseases. However, studying stress response in European seabass may contribute to the understanding of the genetic component of stress and its future use to select broodstock whose offspring may potentially be less affected by handling. A total of 865 European seabass offspring were used to measure body weight and stress response. Moreover, a disease challenge experiment with Vibrio anguillarum was conducted in a subset (332) of the above fish to study disease resistance. Fish were genotyped with a 57k SNP array, and a Genome-Wide Association study (GWAS) was performed. Five SNPs were found to be statistically significant, three of which affect stress indicators and body weight (in a subgroup of the population), and a putative SNP affects growth performance, while no SNP associated with resistance to Vibrio was found. A moderate to high genomic heritability regarding stress indicators and body weight was estimated using the Restricted Maximum Likelihood (REML) process. Finally, the accuracy, along with the correlation between Estimated Breeding Values (EBVs) and Genomic Estimated Breeding Values (GEBVs), were calculated for all the traits.

Chromosome genome assembly for the meagre, Argyrosomus regius, reveals species adaptations and sciaenid sex-related locus evolution.

The meagre, Argyrosomus regius, has recently become a species of increasing economic interest for the Mediterranean aquaculture and there is ongoing work to boost production efficiency through selective breeding. Access to the complete genomic sequence will provide an essential resource for studying quantitative trait-associated loci and exploring the genetic diversity of different wild populations and aquaculture stocks in more detail. Here, we present the first complete genome for A. regius, produced through a combination of long and short read technologies and an efficient in-house developed pipeline for assembly and polishing. Scaffolding using previous linkage map data allowed us to reconstruct a chromosome level assembly with high completeness, complemented with gene annotation and repeat masking. The 696 Mb long assembly has an N50 = 27.87 Mb and an L50 = 12, with 92.85% of its length placed in 24 chromosomes. We use this new resource to study the evolution of the meagre genome and other Sciaenids, via a comparative analysis of 25 high-quality teleost genomes. Combining a rigorous investigation of gene duplications with base-wise conservation analysis, we identify candidate loci related to immune, fat metabolism and growth adaptations in the meagre. Following phylogenomic reconstruction, we show highly conserved synteny within Sciaenidae. In contrast, we report rapidly evolving syntenic rearrangements and gene copy changes in the sex-related dmrt1 neighbourhood in meagre and other members of the family. These novel genomic datasets and findings will add important new tools for aquaculture studies and greatly facilitate husbandry and breeding work in the species.

A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep.

BACKGROUND: The management of farm animal genetic resources and the adaptation of animals to climate change will probably have major effects on the long-term sustainability of the livestock sector. Genomic data harbour useful relevant information that needs to be harnessed for effectively managing genetic resources. In this paper, we report the genome characterization of the highly productive Mediterranean Chios dairy sheep and focus on genetic diversity measures related with local adaptation and selection and the genetic architecture of animal resilience to weather fluctuations as a novel adaptative trait linked to climate change. RESULTS: We detected runs of homozygosity (ROH) and heterozygosity (ROHet) that revealed multiple highly homozygous and heterozygous hotspots across the Chios sheep genome. A particularly highly homozygous region was identified on chromosome 13 as a candidate of directional genetic selection associated with milk traits, which includes annotated genes that were previously shown to be linked to local adaptation to harsh environmental conditions. Favourable heterozygosity related with a potentially protective role against livestock diseases and enhanced overall fitness was revealed in heterozygous-rich regions on sheep chromosomes 3, 10, 13 and 19. Furthermore, genomic analyses were conducted on sheep resilience phenotypes that display changes in milk production in response to weather variation. Sheep resilience to heat stress was a significantly heritable trait (h2 = 0.26) and genetically antagonistic to milk production. Genome-wide association and regional heritability mapping analyses revealed novel genomic markers and regions on chromosome 5 that were significantly associated with sheep resilience to climate change. Subsequently, an annotation analysis detected a set of genes on chromosome 5 that were associated with olfactory receptor complexes that could participate in heat stress mitigation through changes in respiration rate and respiratory evaporation. Other genes were grouped in previously reported biological processes relevant to livestock heat dissipation, including stress and immune response. CONCLUSIONS: Our results may contribute to the optimal management of sheep genetic resources and inform modern selective breeding programmes that aim at mitigating future environmental challenges towards sustainable farming, while better balancing animal adaptation and productivity. Our results are directly relevant to the studied breed and the respective environmental conditions; however, the methodology may be extended to other livestock species of interest.

Understanding the seasonality of performance resilience to climate volatility in Mediterranean dairy sheep.

As future climate challenges become increasingly evident, enhancing performance resilience of farm animals may contribute to mitigation against adverse weather and seasonal variation, and underpin livestock farming sustainability. In the present study, we develop novel seasonal resilience phenotypes reflecting milk production changes to fluctuating weather. We evaluate the impact of calendar season (autumn, winter and spring) on animal performance resilience by analysing 420,534 milk records of 36,908 milking ewes of the Chios breed together with relevant meteorological data from eastern Mediterranean. We reveal substantial seasonal effects on resilience and significant heritable trait variation (h2 = 0.03-0.17). Resilience to cold weather (10 °C) of animals that start producing milk in spring was under different genetic control compared to autumn and winter as exemplified by negative genetic correlations (- 0.09 to - 0.27). Animal resilience to hot weather (25 °C) was partially under the same genetic control with genetic correlations between seasons ranging from 0.43 to 0.86. We report both favourable and antagonistic associations between animal resilience and lifetime milk production, depending on calendar season and the desirable direction of genetic selection. Concluding, we emphasise on seasonal adaptation of animals to climate and the need to incorporate the novel seasonal traits in future selective breeding programmes.

QTL for Stress and Disease Resistance in European Sea Bass, Dicentrarhus labrax L.

There is a growing interest in selective breeding in European sea bass (Dicentrarchus labrax), especially regarding family selection based on growth performance. In particular, quantitative trait loci (QTL) identification in sea bass enhances the application of marker-assisted breeding for the genetic improvement of the production traits. The aims of the study were to identify potential QTL affecting stress and immunological indicators, body weight, and mortality after vibriosis injection in sea bass as well as to estimate heritability and genetic/phenotypic correlations for the aforementioned traits. To this end, stress test was performed on 960 offspring and a sub-group of them (420) was selected to explore the mortality after vibrio injection. Selective genotyping was performed in 620 offspring for 35 microsatellite markers and distributed into 6 linkage groups. The length of the genetic linkage map was 283.6 cM and the mean distance between the markers was 8.1 cM. QTL affecting body weight in three different growth periods detected on linkage groups LG1, LG4, LG6, and LG14. A QTL associated with weight in early growth stages (290-306 days post-hatching) was also identified on LG3. QTL analysis confirmed the existence of QTL affecting cortisol levels, on LG3 and LG14. Moreover, new QTL affecting only cortisol and glucose levels were detected on LG1 and LG23. No QTL affecting hormonal or biochemical marks was found on LG4 and LG6. Heritability of cortisol, lysozyme levels, and mortality were high (0.36, 0.55, and 0.38, respectively).

A practical approach to detect ancestral haplotypes in livestock populations.

BACKGROUND: The effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle. RESULTS: The identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content. CONCLUSION: In agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations.

Exploring a Nonmodel Teleost Genome Through RAD Sequencing-Linkage Mapping in Common Pandora, Pagellus erythrinus and Comparative Genomic Analysis.

Common pandora (Pagellus erythrinus) is a benthopelagic marine fish belonging to the teleost family Sparidae, and a newly recruited species in Mediterranean aquaculture. The paucity of genetic information relating to sparids, despite their growing economic value for aquaculture, provides the impetus for exploring the genomics of this fish group. Genomic tool development, such as genetic linkage maps provision, lays the groundwork for linking genotype to phenotype, allowing fine-mapping of loci responsible for beneficial traits. In this study, we applied ddRAD methodology to identify polymorphic markers in a full-sib family of common pandora. Employing the Illumina MiSeq platform, we sampled and sequenced a size-selected genomic fraction of 99 individuals, which led to the identification of 920 polymorphic loci. Downstream mapping analysis resulted in the construction of 24 robust linkage groups, corresponding to the karyotype of the species. The common pandora linkage map showed varying degrees of conserved synteny with four other teleost genomes, namely the European seabass (Dicentrarchus labrax), Nile tilapia (Oreochromis niloticus), stickleback (Gasterosteus aculeatus), and medaka (Oryzias latipes), suggesting a conserved genomic evolution in Sparidae. Our work exploits the possibilities of genotyping by sequencing to gain novel insights into genome structure and evolution. Such information will boost the study of cultured species and will set the foundation for a deeper understanding of the complex evolutionary history of teleosts.

Second generation genetic linkage map for the gilthead sea bream Sparus aurata L.

An updated second linkage map was constructed for the gilthead sea bream, Sparus aurata L., a fish species of great economic importance for the Mediterranean aquaculture industry. In contrast to the first linkage map which mainly consisted of genomic microsatellites (SSRs), the new linkage map is highly enriched with SSRs found in Expressed Sequence Tags (EST-SSRs), which greatly facilitates comparative mapping with other teleosts. The new map consists of 321 genetic markers in 27 linkage groups (LGs): 232 genomic microsatellites, 85 EST-SSRs and 4 SNPs; of those, 13 markers were linked to LGs but were not ordered. Eleven markers (5 SSRs, 5 EST-SSRs and 1 SNP) are not assigned to any LG. The total length of the sex-averaged map is 1769.7cM, 42% longer than the previously published one, and the number of markers in each LG ranges from 2 to 30. The inter-marker distance varies from 0 to 75.6cM, with an average of 5.75cM. The male and female maps have a length of 1349.2 and 2172.1cM, respectively, and the average distance between markers is 4.38 and 7.05cM, respectively. Comparative mapping with the three-spined stickleback (Gasterosteus acuulatus) chromosomes and scaffolds showed conserved synteny with 132 S. aurata markers (42.9% of those mapped) having a hit on the stickleback genome.

Quantitative trait loci involved in sex determination and body growth in the gilthead sea bream (Sparus aurata L.) through targeted genome scan.

Among vertebrates, teleost fish exhibit a considerably wide range of sex determination patterns that may be influenced by extrinsic parameters. However even for model fish species like the zebrafish Danio rerio the precise mechanisms involved in primary sex determination have not been studied extensively. The zebrafish, a gonochoristic species, is lacking discernible sex chromosomes and the sex of juvenile fish is difficult to determine. Sequential protandrous hermaphrodite species provide distinct determination of the gender and allow studying the sex determination process by looking at the mechanism of sex reversal. This is the first attempt to understand the genetic basis of phenotypic variation for sex determination and body weight in a sequential protandrous hermaphrodite species, the gilthead sea bream (Sparus aurata). This work demonstrates a fast and efficient strategy for Quantitative Trait Loci (QTL) detection in the gilthead sea bream, a non-model but target hermaphrodite fish species. Therefore a comparative mapping approach was performed to query syntenies against two other Perciformes, the European sea bass (Dicentrarchus labrax), a gonochoristic species and the Asian sea bass (Lates calcarifer) a protandrous hermaphrodite. In this manner two significant QTLs, one QTL affecting both body weight and sex and one QTL affecting sex, were detected on the same linkage group. The co-segregation of the two QTLs provides a genomic base to the observed genetic correlation between these two traits in sea bream as well as in other teleosts. The identification of QTLs linked to sex reversal and growth, will contribute significantly to a better understanding of the complex nature of sex determination in S. aurata where most individuals reverse to the female sex at the age of two years through development and maturation of the ovarian portion of the gonad and regression of the testicular area. [Genomic sequences reported in this manuscript have been submitted to GenBank under accession numbers HQ021443-HQ021749.].