Almudena Fernández

Epizootic hemorrhagic disease (EHD) is a non-contagious viral infection that can cause important economic losses in dairy farms. This study aimed to identify epidemiological and genetic factors influencing the susceptibility and severity of EHD in Holstein dairy cattle during the 2023 outbreak in Spain. Data from 2852 animals in 7 affected farms from 5 Spanish provinces were used. Symptoms were categorized in 5 categories: no symptoms, mild symptoms, severe symptoms and recovered, severe symptoms with aftereffects (sequelae) and severe symptoms followed by death. All animals were genotyped using the Illumina EuroG MD SNP array, and imputation to whole genome sequenced was carried out using the 1000 bull genomes data set as reference. Risk factors for EHD were explored using linear mixed effects models, as well as the loss of milk yield and culling risk probability due to severe EHD. Around 66% of animals showed EHDV antibodies, although only 25% of the sample in this study showed severe symptoms and a death rate of 2-3%. Results indicated that age was the main risk factor for severe EHD, with older cows showing higher susceptibility. Production losses were significant in cows with moderate to severe symptoms, especially if the outbreak occurred mid-lactation (up to -9 kg/d). Cows exhibiting severe symptoms demonstrated a markedly increased likelihood of being culled, with odds ratio of 10.86 (95% CI: 6.08-19.41) for cows with severe symptoms. The genetic component of EHD was evaluated using REML procedures, and a genome-wide association study was conducted to investigate genomic regions associated to the disease. A heritability of 0.08 was estimated, with some genes (AOC1, BST1, CD38, DPP6, DPYS, HOGA1, KMT2C, PIK2A, PIK3, PI4K2A) associated with the disease that were involved in immunity processes and the development of the pulmonary tree. This study underscores the genetic and epidemiological factors influencing EHD severity in dairy cattle, providing insights for managing future outbreaks and mitigating economic losses in affected regions.

Introduction: The low pregnancy rate by artificial insemination in sheep represents a fundamental challenge for breeding programs. In this species, oestrus synchronization is carried out by manipulating hormonal regimens through the insertion of progestogen intravaginal devices. This reproductive strategy may alter the vaginal microbiota affecting the artificial insemination outcome.Methods: In this study, we analyzed the vaginal microbiome of 94 vaginal swabs collected from 47 ewes with alternative treatments applied to the progesterone-releasing intravaginal devices (probiotic, maltodextrin, antibiotic and control), in two sample periods (before placing and after removing the devices). To our knowledge, this is the first study using nanopore-based metagenome sequencing for vaginal microbiome characterization in livestock.Results: Our results revealed a significant lower abundance of the genera Oenococcus (Firmicutes) and Neisseria (Proteobacteria) in pregnant compared to non-pregnant ewes. We also detected a significant lower abundance of Campylobacter in the group of samples treated with the probiotic.Discussion: Although the use of probiotics represents a promising practice to improve insemination results, the election of the suitable species and concentration requires further investigation. In addition, the use of progestogen in the synchronization devices seemed to increase the alpha-diversity and decrease the abundance of harmful microorganisms belonging to Gammaproteobacteria and Fusobacteriia classes, suggesting a beneficial effect of their use.

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

Background: A complete understanding of the genetic basis for sexual determination and differentiation is necessary in order to implement efficient breeding schemes at early stages of development. Atlantic salmon belongs to the family Salmonidae of fishes and represents a species of great commercial value. Although the species is assumed to be male heterogametic with XY sex determination, the precise genetic basis of sexual development remains unclear. The complexity is likely associated to the relatively recent salmonid specific whole genome duplication that may be responsible for certain genome instability. This instability together with the capacity of the sex-determining gene to move across the genome as reported by previous studies, may explain that sexual development genes are not circumscribed to the same chromosomes in all members of the species. In this study, we have used a 220 K SNP panel developed for Atlantic salmon to identify the chromosomes explaining the highest proportion of the genetic variance for sex as well as candidate regions and genes associated to sexual development in this species.Results: Results from regional heritability analysis showed that the chromosomes explaining the highest proportion of variance in these populations were Ssa02 (heritability = 0.42, SE = 0.12) and Ssa21 (heritability = 0.26, SE = 0.11). After pruning by linkage disequilibrium, genome-wide association analyses revealed 114 SNPs that were significantly associated with sex, being Ssa02 the chromosome containing a greatest number of regions. Close examination of the candidate regions evidenced important genes related to sex in other species of Class Actinopterygii, includingSDY, genes from familySOX, RSPO1, ESR1, U2AF2A, LMO7, GNRH-R, DNDandFIGLA.Conclusions: The combined results from regional heritability analysis and genome-wide association have provided new advances in the knowledge of the genetic regulation of sex determination in Atlantic salmon, supporting that Ssa02 is the candidate chromosome for sex in this species and suggesting an alternative population lineage in Spanish wild populations according to the results from Ssa21.

Selective breeding for improving host responses to infectious pathogens is a promising option for disease control. In fact, disease resilience, the ability of a host to survive or cope with infectious challenge, has become a highly desirable breeding goal. However, resilience is a complex trait composed of two different host defence mechanisms, namely resistance (the ability of a host to avoid becoming infected or diseased) and endurance (the ability of an infected host to survive the infection). While both could be targeted for genetic improvement, it is currently unknown how they contribute to survival, as reliable estimates of genetic parameters for both traits obtained simultaneously are scarce. A difficulty lies in obtaining endurance phenotypes for genetic analyses. In this study, we present the results from an innovative challenge test carried out in turbot whose design allowed disentangling the genetic basis of resistance and endurance to Philasterides dicentrarchi, a parasite causing scuticociliatosis that leads to substantial economic losses in the aquaculture industry. A noticeable characteristic of the parasite is that it causes visual signs that can be used for disentangling resistance and endurance. Our results showed the existence of genetic variation for both traits (heritability = 0.26 and 0.12 for resistance and endurance, respectively) and for the composite trait resilience (heritability = 0.15). The genetic correlation between resistance and resilience was very high (0.90) indicating that both are at a large extent the same trait, but no significant genetic correlation was found between resistance and endurance. A total of 18,125 SNPs obtained from 2b-RAD sequencing enabled genome-wide association analyses for detecting QTLs controlling the three traits. A candidate QTL region on linkage group 19 that explains 33% of the additive genetic variance was identified for resilience. The region contains relevant genes related to immune response and defence mechanisms. Although no significant associations were found for resistance, the pattern of association was the same as for resilience. For endurance, one significant association was found on linkage group 2. The accuracy of genomic breeding values was also explored for resilience, showing that it increased by 12% when compared with the accuracy of pedigree-based breeding values. To our knowledge, this is the first study in turbot disentangling the genetic basis of resistance and endurance to scuticociliatosis.

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

One of the major limitation for the application of QTL results in pig breeding and QTN identification has been the limited number of QTL effects validated in different animal material. The aim of the current work was to validate QTL regions through joint and specific genome wide association and haplotype analyses for growth, fatness and premier cut weights in three different genetic backgrounds, backcrosses based on Iberian pigs, which has a major role in the analysis due to its high productive relevance. The results revealed nine common QTL regions, three segregating in all three backcrosses on SSC1, 0–3 Mb, for body weight, on SSC2, 3–9 Mb, for loin bone-in weight, and on SSC7, 3 Mb, for shoulder weight, and six segregating in two of the three backcrosses, on SSC2, SSC4, SSC6 and SSC10 for backfat thickness, shoulder and ham weights. Besides, 18 QTL regions were specifically identified in one of the three backcrosses, five identified only in BC_LD, seven in BC_DU and six in BC_PI. Beyond identifying and validating QTL, candidate genes and gene variants within the most interesting regions have been explored using functional annotation, gene expression data and SNP identification from RNA-Seq data. The results allowed us to propose a promising list of candidate mutations, those identified in PDE10A, DHCR7, MFN2 and CCNY genes located within the common QTL regions and those identified near ssc-mir-103-1 considered PANK3 regulators to be further analysed.

Meat quality depends on tissue composition which is in turn influenced by different factors, such as diet, genotype, age, or sex. We evaluated the effects of breed, 24 h fasting, and dietary energy source (HO: oleic acid versus CH: carbohydrates) on the expression of candidate genes involved in adipogenesis, lipogenesis, and lipolysis in the adipose tissue from Iberian and Duroc growing pigs. The Iberian pigs showed greater feed intake, backfat thickness, and saturated fatty acids (SFA) content in the subcutaneous fat, whereas the Duroc pigs had greater ham weight and polyunsaturated fatty acids (PUFA) content. In both breeds, the diet induced changes in the fatty acid (FA) composition of subcutaneous fat samples. The HO group had higher monounsaturated fatty acids (MUFA) and oleic acid, and lower SFA than the CH group. Regarding gene expression, breed and feeding status (fasting versus postprandial) had significant effects on gene expression, with quantitative interactions between them, while diet showed negligible effects. In general, adipogenic and lipogenic genes were upregulated in the Iberian pigs and in postprandial samples. In contrast, the expression of lipolytic genes showed complex interaction effects. Our results agree with the phenotypic differences between the Iberian and Duroc breeds and with the inhibition of lipogenesis by fasting. Quantitative interactions between breed and feeding status effects were observed, which indicates a different response to fasting of the two breeds, with the obese Iberian breed showing a more stable expression of lipogenic genes. These results highlight the complexity of lipid metabolism regulation, especially in relation to lipolysis processes.