Inbreeding depression load for litter size in Entrepelado and Retinto Iberian pig varieties1

Individual-specific hidden inbreeding depression load (IDL) can be accounted for in livestock populations by appropriate best linear unbiased prediction approaches. This genetic effect has a recessive pattern and reveals when inherited in terms of identity-by-descent. Nevertheless, IDL inherits as a pure additive genetic background and can be selected using standard breeding values. The main target of this research was to evaluate IDL for litter size in 2 Iberian pig varieties (Entrepelado and Retinto) from a commercial breeding-stock. Analyses were performed on the total number of piglets born (both alive and dead) and used data from 3,200 (8.02 ± 0.04 piglets/litter) Entrepelado and 4,744 Retinto litters (8.40 ± 0.03 piglets/litter). Almost 50% of Entrepelado sows were inbred (1.7% to 25.0%), whereas this percentage reduced to 37.4% in the Retinto variety (0.2% to 25.0%). The analytical model was solved by Bayesian inference and accounted for 2 systematic effects (sow age and breed/variety of the artificial insemination boar), 2 permanent environmental effects (herd-year-season and sow), and 2 genetic effects (IDL and infinitesimal additive). In terms of posterior means (PM), additive genetic and IDL variances were similar in the Entrepelado variety (PM, 0.68 vs. 0.76 piglets2, respectively) and their 95% credibility intervals (95CI) overlapped, although without including zero (0.38 to 0.94 vs. 0.15 to 1.31 piglets2, respectively). The same pattern revealed in the Retinto variety, with IDL variance (PM, 0.41 piglets2; 95CI, 0.07 to 0.88 piglets2) slightly larger than the additive genetic variance (PM, 0.37 piglets2; 95CI, 0.16 to 0.59 piglets2). The relevance of IDL was also checked by a Bayes factor and the deviance information criterion, the model including this effect being clearly favored in both cases. Although the analysis assumed null genetic covariance between IDL and infinitesimal additive effects, a moderate negative correlation (-0.31) was suggested when plotting the PM of breeding values in the Entrepelado variety; a negative genetic trend for IDL was also revealed in this Iberian pig variety (-0.25 piglets for 100% inbred offspring of individuals born in 2014), whereas no trend was detected in Retinto breeding-stock. Those were the first estimates of IDL in a commercial livestock population, they giving evidence of a relevant genetic background with potential consequences on the reproductive performance of Iberian sows.

Low genome-wide homozygosity in 11 Spanish ovine breeds

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (F_ROH ) are generally small (F_ROH  < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

Background

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

Results

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110–114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

Conclusions

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5557-9) contains supplementary material, which is available to authorized users.

Implementation of Bayesian methods to identify SNP and haplotype regions with transmission ratio distortion across the whole genome: TRDscan v.1.0

Realized deviations from the expected Mendelian inheritance of alleles from heterozygous parents have been previously reported in a broad range of organisms (i.e., transmission ratio distortion; TRD). Various biological mechanisms affecting gametes, embryos, fetuses, or even postnatal offspring can produce patterns of TRD. However, knowledge about its prevalence and potential causes in livestock species is still scarce. Specific Bayesian models have been recently developed for the analyses of TRD for biallelic loci, which accommodated a wide range of population structures, enabling TRD investigation in livestock populations. The parameterization of these models is flexible and allows the study of overall (parent-unspecific) TRD and sire- and dam-specific TRD. This research aimed at deriving Bayesian models for fitting TRD on the basis of haplotypes, testing the models for both haplotype- and SNP-based methods in simulated data and actual Holstein genotypes, and developing a specific software for TRD analyses. Results obtained on simulated data sets showed that the statistical power of the analysis increased with sample size of trios (n), proportion of heterozygous parents, and the magnitude of the TRD. On the other hand, the statistical power to detect TRD decreased with the number of alleles at each loci. Bayesian analyses showed a strong Pearson correlation coefficient (≥0.97) between simulated and estimated TRD that reached the significance level of Bayes factor ≥10 for both single-marker and haplotype analyses when n ≥ 25. Moreover, the accuracy in terms of the mean absolute error decreased with the increase of the sample size and increased with the number of alleles at each loci. Using real data (55,732 genotypes of Holstein trios), SNP- and haplotype-based distortions were detected with overall TRD, sire-TRD, or dam-TRD, showing different magnitudes of TRD and statistical relevance. Additionally, the haplotype-based method showed more ability to capture TRD compared with individual SNP. To discard possible random TRD in real data, an approximate empirical null distribution of TRD was developed. The program TRDscan v.1.0 was written in Fortran 2008 language and provides a powerful statistical tool to scan for TRD regions across the whole genome. This developed program is freely available at http://www.casellas.info/files/TRDscan.zip.

Screening for epistatic selection signatures: A simulation study

Detecting combinations of alleles that diverged between subpopulations via selection signature statistics can contribute to decipher the phenomenon of epistasis. This research focused on the simulation of genomic data from subpopulations under divergent epistatic selection (ES). We used D’_IS² and F_ST statistics in pairs of loci to scan the whole-genome. The results showed the ability to identify loci under additive-by-additive ES (ES_aa) by reporting large statistical departures between subpopulations with a high level of divergence, while it did not show the same advantage in the other types of ES. Despite this, limitations such as the difficulty to distinguish between the quasi-complete fixation of one locus by ES_aa from other events were observed. However, D’_IS² can detect loci under ES_aa by defining a minimum boundary for the minor allele frequency on a multiple subpopulation analysis where ES only takes place in one subset. Even so, the major limitation was distinguishing between ES and single-locus selection (SS); therefore, we can conclude that divergent locus can be also a result of ES. The test conditions with D-statistics of both Ohta (1982a, 1982b) and Black and Krafsur (1985) did not provide evidence to differentiate ES in our simulation framework of isolated subpopulations.

Freemartinism in replacement ewe-lambs of the Ripollesa sheep breed

The freemartinism syndrome affects almost all female calves born as co-twins to male calves, whereas little is known about this phenomenon in female sheep. Within this context, 1,185 ewe-lambs from the Ripollesa sheep breed were genotyped for the presence of oY1 polymorphism (a non-autosomal region of the Y chromosome). Neither ewe-lambs from single births (856) nor ewe-lambs from all-female multiple births (170) were revealed as freemartins, whereas five of 159 ewe-lambs from multiple births with male co-twins were freemartins (3.15 ± 1.38%). All freemartin ewe-lambs were confirmed by physical examination of external genitalia. The results confirm a low incidence of freemartinism from heterosexual twin pregnancies in Ripollesa sheep.

12 Effect of additional prostaglandin F2α during the Ovsynch protocol applied in different postpartum intervals in lactating dairy cows: Preliminary results

Timed AI (TAI) programs have been developed to increase the service risk. However, when Ovsynch is used as a synchronization protocol, lack of regression of the corpus luteum after prostaglandin F2α (PGF) administration has been demonstrated, leading to reductions in fertility. In this way, a second PGF administration might increase the frequency of achieving complete luteolysis. This study aimed to evaluate the effect of an additional treatment with PGF on fertility when an Ovsynch synchronization protocol was applied at different postpartum intervals. Multiparous high-production dairy cows (n=471) of a commercial farm were stratified by postpartum intervals [&lt;120 (n=136), 120-150 (n=131), 150-180 (n=92), 180-210 (n=52), and &gt;210 days (n=60)] and randomly allocated to receive either the Ovsynch (n=108) or Ovsynch with second PGF protocol (n=243). Cows inseminated 12h after observed heat detection (AI) served as a control (n=120). The Ovsynch protocol consisted of an injection of gonadotropin-releasing hormone (100 µg; Fertagyl®, Intervet, Millsboro, MD, USA) on Day 0, an injection of PGF (0.5 mg; Estrumate®, Schering-Plough Animal Health, Montréal, QC, Canada) on Day 7, another injection of gonadotropin-releasing hormone (100µg) on Day 9, and timed insemination on Day 10. The Ovsynch with second PGF protocol consisted of the same hormone injection schedule as the Ovsynch protocol, but a second injection of PGF (0.5 mg; Estrumate®) was administered 24h after the first PGF. Data were analysed using the chi-square test. The level of significance was taken as P&lt;0.05. When data were pooled at 120 to 150 days of postpartum interval, pregnancy rate was higher, but not significantly different (P&gt;0.05), for cows that received the Ovsynch (37.5%, n=32) or Ovsynch with second PGF protocol (33.8%, n=77) compared to nonsynchronized cows (22.8%, n=22). After 180 days of postpartum interval, Ovsynch with second PGF protocol triggered pregnancy rates similar to those of nonsynchronized cows (33.3%, n=21 and 28.6%, n=14), whereas the Ovsynch protocol resulted in lower percentages (11.8%, n=17; P&gt;0.05). Similar pregnancy rates (P&gt;0.05) were observed when TAI protocols were administered before 120 days (26.5 and 30.2% for Ovsynch and Ovsynch with second PGF, respectively) or between 150 and 180 postpartum days (31.8 and 33.9% for Ovsynch and Ovsynch with second PGF, respectively) compared with AI cows (25.7 and 28.8% for &lt;120 days and between 150 and 180 postpartum days, respectively). These preliminary results showed a tendency of a better pregnancy rate when both TAI synchronization protocols were applied between 120 and 150 postpartum days when compared with AI. Also, the application of a second dose of PGF seems recommendable when the Ovsynch protocol is used to synchronize cows of &gt;180 postpartum days. Ongoing research will allow an increase in the number of the samples in each postpartum interval in order to confirm these results.

Combining multi-OMICs information to identify key-regulator genes for pleiotropic effect on fertility and production traits in beef cattle

The identification of biological processes related to the regulation of complex traits is a difficult task. Commonly, complex traits are regulated through a multitude of genes contributing each to a small part of the total genetic variance. Additionally, some loci can simultaneously regulate several complex traits, a phenomenon defined as pleiotropy. The lack of understanding on the biological processes responsible for the regulation of these traits results in the decrease of selection efficiency and the selection of undesirable hitchhiking effects. The identification of pleiotropic key-regulator genes can assist in developing important tools for investigating biological processes underlying complex traits. A multi-breed and multi-OMICs approach was applied to study the pleiotropic effects of key-regulator genes using three independent beef cattle populations evaluated for fertility traits. A pleiotropic map for 32 traits related to growth, feed efficiency, carcass and meat quality, and reproduction was used to identify genes shared among the different populations and breeds in pleiotropic regions. Furthermore, data-mining analyses were performed using the Cattle QTL database (CattleQTLdb) to identify the QTL category annotated in the regions around the genes shared among breeds. This approach allowed the identification of a main gene network (composed of 38 genes) shared among breeds. This gene network was significantly associated with thyroid activity, among other biological processes, and displayed a high regulatory potential. In addition, it was possible to identify genes with pleiotropic effects related to crucial biological processes that regulate economically relevant traits associated with fertility, production and health, such as MYC, PPARG, GSK3B, TG and IYD genes. These genes will be further investigated to better understand the biological processes involved in the expression of complex traits and assist in the identification of functional variants associated with undesirable phenotypes, such as decreased fertility, poor feed efficiency and negative energetic balance.

Genetic evaluation of age at first calving for Guzerá beef cattle using linear, threshold, and survival Bayesian models

Age at first calving (AFC) is characterized as a censored trait due to missing values provided by recording mistakes and nonoccurrence or delay in calving communication. In this context, we aimed to compare several statistical methods for genetic evaluation of AFC in Guzerá beef cattle under a Bayesian approach. Seven different methods were used for this purpose. The traditional linear mixed model (LM), which considers only uncensored records; the LM with simulated records (SM), which is based on data augmentation framework; the penalty method, in which a constant of 21 d was added to censored records; the bivariate threshold-linear method considering (TLcens) or not (TLmiss) censored information; and the piecewise Weibull proportional hazards model considering (PWPHcens) or not (PWPH) censored records. Heritability estimates ranged from 0.19 (TLcens) to 0.28 (SM) in nonsurvival approaches; and 0.40 and 0.46 to PWPH and PWPHcens methods, respectively. In general, breeding values correlations between different methods and the percentage of selected bulls in common indicated reranking, with these correlation ranging from -0.28 (between SM and PWPH) to 0.99 (between TLmiss and LM). The traditional LM, which considers only uncensored records, should be preferred due to its robustness and simplicity. Based on cross-validation analyses, we conclude that the TLmiss could be also a suitable alternative for breeding value prediction, and censored methods did not improve the analysis.

Skin mites in mice (Mus musculus): high prevalence of Myobia sp. (Acari, Arachnida) in Robertsonian mice

Myobia sp. and Demodex sp. are two skin mites that infest mice, particularly immunodeficient or transgenic lab mice. In the present study, wild house mice from five localities from the Barcelona Roberstonian system were analysed in order to detect skin mites and compare their prevalence between standard (2n = 40) and Robertsonian mice (2n > 40). We found and identified skin mites through real-time qPCR by comparing sequences from the mitochondrial 16S rRNA and the nuclear 18S rRNA genes since no sequences are available so far using the mitochondrial gene. Fourteen positive samples were identified as Myobia musculi except for a deletion of 296 bp out to 465 bp sequenced, and one sample was identified as Demodex canis. Sampling one body site, the mite prevalence in standard and Robertsonian mice was 0 and 26%, respectively. The malfunction of the immune system elicits an overgrowth of skin mites and consequently leads to diseases such as canine demodicosis in dogs or rosacea in humans. In immunosuppressed mice, the probability of developing demodicosis is higher than in healthy mice. Since six murine toll-like receptors (TLRs) are located in four chromosomes affected by Robertsonian fusions, we cannot dismiss that differences in mite prevalence could be the consequence of the interruption of TLR function. Although ecological and/or morphological factors cannot be disregarded to explain differences in mite prevalence, the detection of translocation breakpoints in TLR genes or the analysis of TLR gene expression are needed to elucidate how Robertsonian fusions affect the immune system in mice.

On individual-specific prediction of hidden inbreeding depression load

Inbreeding depression is caused by increased homozygosity in the genome and merges two genetic mechanisms, a higher impact from recessive mutations and the waste of overdominance contributions. It is of major concern for the conservation of endangered populations of plants and animals, as major abnormalities are more frequent in inbred families than in outcrosses. Nevertheless, we lack appropriate analytical methods to estimate the hidden inbreeding depression load (IDL) in the genome of each individual. Here, a new mixed linear model approach has been developed to account for the inbreeding depression-related background of each individual in the pedigree. Within this context, inbred descendants contributed relevant information to predict the IDL contained in the genome of a given ancestor; moreover, known relationships spread these predictions to the remaining individuals in the pedigree, even if not contributing inbred offspring. Results obtained from the analysis of weaning weight in the MARET rabbit population demonstrated that the genetic background of inbreeding depression distributed heterogeneously across individuals and inherited generation by generation. Moreover, this approach was clearly preferred in terms of model fit and complexity when compared with classical approaches to inbreeding depression. This methodology must be viewed as a new tool for a better understanding of inbreeding in domestic and wild populations.

Technical note: PaGELL v.1.5: A flexible parametric program for the Bayesian analysis of longevity data within the context of animal breeding

This technical note presents the program PaGELL v.1.5 (Parametric Genetic Evaluation of Lifespan in Livestock), a flexible software program to analyze (right-censored) longevity data in livestock populations, with a special emphasis on the genetic evaluation of the breeding stock. This software relies on a parametric generalization of the proportional hazard model; more specifically, the baseline hazard function follows a Weibull process and flexibility is gained by including an additional time-dependent effect with the number of change points defined by the user. The program can accommodate 3 different sources of variation (i.e., systematic, permanent environmental, and additive genetic effects) and both fixed and time-dependent patterns (only for systematic and permanent environmental effects). Analyses are performed within a Bayesian context by sampling from the joint posterior distribution of the model, and model fit can be easily determined by the calculation of the deviance information criterion. Although this software has already been used on field data sets, its performance has been double-checked on simulated data set, and results are presented in this technical note. PaGELL v.1.5 was written in Fortran 95 language and, after compiling with the GNU Fortran Compiler v.4.7 and later, it has been tested in Windows, Linux, and MacOS operating systems (both 32- and 64-bit platforms). This program is available at http://www.casellas.info/files/pageII.zip.

Linkage disequilibrium, persistence of phase, and effective population size in Spanish local beef cattle breeds assessed through a high-density single nucleotide polymorphism chip

Linkage disequilibrium (LD) and persistence of phase are fundamental approaches for exploring the genetic basis of economically important traits in cattle, including the identification of QTL for genomic selection and the estimation of effective population size () to determine the size of the training populations. In this study, we have used the Illumina BovineHD chip in 168 trios of 7 Spanish beef cattle breeds to obtain an overview of the magnitude of LD and the persistence of LD phase through the physical distance between markers. Also, we estimated the time of divergence based on the persistence of the LD phase and calculated past from LD estimates using different alternatives to define the recombination rate. Estimates of average (as a measure of LD) for adjacent markers were close to 0.52 in the 7 breeds and decreased with the distance between markers, although in long distances, some LD still remained (0.07 and 0.05 for markers 200 kb and 1 Mb apart, respectively). A panel with a lower boundary of 38,000 SNP would be necessary to launch a successful within-breed genomic selection program. Persistence of phase, measured as the pairwise correlations between estimates of in 2 breeds at short distances (10 kb), was in the 0.89 to 0.94 range and decreased from 0.33 to 0.52 to a range of 0.01 to 0.08 when marker distance increased from 200 kb to 1 Mb, respectively. The magnitude of the persistence of phase between the Spanish beef breeds was similar to those found in dairy breeds. For across-breed genomic selection, the size of the SNP panels must be in the range of 50,000 to 83,000 SNP. Estimates of past showed values ranging from 26 to 31 for 1 generation ago in all breeds. The divergence among breeds occurred between 129 and 207 generations ago. The results of this study are relevant for the future implementation of within- and across-breed genomic selection programs in the Spanish beef cattle populations. Our results suggest that a reduced subset of the SNP panel would be enough to achieve an adequate precision of the genomic predictions.

Bayesian analysis of parent-specific transmission ratio distortion in seven Spanish beef cattle breeds

Transmission ratio distortion (TRD) is the departure from the expected Mendelian ratio in offspring, a poorly investigated biological phenomenon in livestock species. Given the current availability of specific parametric methods for the analysis of segregation data, this study focused on the screening of TRD in 602 402 single nucleotide polymorphisms covering all autosomal chromosomes in seven Spanish beef cattle breeds. On average, 0.13% (n = 786) and 0.01% (n = 29) of genetic markers evidenced sire- or dam-specific TRD respectively. There were no single nucleotide polymorphisms accounting for both sire- and dam-specific TRD at the same time, and only one marker (rs43147474) accounted for (sire-specific) TRD in all seven breeds. It must be noted that rs43147474 is located in the fourth intronic region of the GTP-binding protein 10 gene, and this locus has been previously linked to the maintenance of mitochondria and nucleolar architectures. Alternatively, other candidate genes surround this hot-spot for sire-specific TRD in the cattle genome, and they are related to embryonic and postnatal lethality as well as prostate cancer, among others. This research characterized the distribution of TRD in the bovine genome, highlighting heterogeneous results when comparing across breeds.

Genetic and environmental relationships of detailed milk fatty acids profile determined by gas chromatography in Brown Swiss cows

The aim of this study was to characterize the profile of 47 fatty acids, including conjugated linoleic acid (CLA), 13 fatty acid groups, and 5 Δ(9)-desaturation indices in milk samples from Brown Swiss cows. The genetic variation was assessed and the statistical relevance of the genetic background for each trait was evaluated using the Bayes factor test. The additive genetic, herd-date, and residual relationships were also estimated among all single fatty acids and groups of fatty acids. Individual milk samples were collected from 1,158 Italian Brown Swiss cows and a detailed analysis of fat percentages and milk fatty acid compositions was performed by gas chromatography. Bayesian animal models were used for (co)variance components estimation. Exploitable genetic variation was observed for most of the de novo synthesized fatty acids and saturated fatty acids, except for C4:0 and C6:0, whereas long-chain fatty acids and unsaturated fatty acids (including CLA) were mainly influenced by herd-date effects. Herd-date effect explained large portions of the total phenotypic variance for C18:2 cis-9,cis-12 (0.668), C18:3 cis-9,cis-12,cis-15 (0.631), and the biohydrogenation and elongation products of these fatty acids. The desaturation ratios showed higher heritability estimates than the individual fatty acids, except for CLA desaturation index (0.098). Among the medium-chain fatty acids, C12:0 had greater heritability than C14:0 (0.243 vs. 0.097, respectively). Both C14:0 and C16:0 showed negative additive genetic correlations with the main monounsaturated and polyunsaturated fatty acids of milk fat, suggesting that their synthesis in the mammary gland may be influenced by the presence of unsaturated fatty acids. No correlation was observed between C4:0 and the other short-chain fatty acids (except for C6:0), confirming the independence of C4:0 from de novo mammary fatty acid synthesis. Among the genetic correlations dealing with potentially beneficial fatty acids, C18:0 was positively correlated with vaccenic and rumenic acids and negatively with linoleic acid. Finally, fatty acids C6:0 through C14:0 showed relevant correlations due to unknown environmental effects, suggesting the potential existence of genetic variances in micro-environmental sensitivity. This study allowed us to acquire new knowledge about the genetic and the environmental relationships among fatty acids. Likewise, the existence of genetic variation for most of de novo synthetized fatty acids and saturated fatty acids was also observed. Overall, these results provide useful information to combine feeding with genetic selection strategies for obtaining a desirable milk fatty acids profile, depending on the origin of fatty acids in milk.