A. González-Rodríguez

A Bayesian Model for the Analysis of Transgenerational Epigenetic Variation

Epigenetics has become one of the major areas of biological research. However, the degree of phenotypic variability that is explained by epigenetic processes still remains unclear. From a quantitative genetics perspective, the estimation of variance components is achieved by means of the information provided by the resemblance between relatives. In a previous study, this resemblance was described as a function of the epigenetic variance component and a reset coefficient that indicates the rate of dissipation of epigenetic marks across generations. Given these assumptions, we propose a Bayesian mixed model methodology that allows the estimation of epigenetic variance from a genealogical and phenotypic database. The methodology is based on the development of a T matrix of epigenetic relationships that depends on the reset coefficient. In addition, we present a simple procedure for the calculation of the inverse of this matrix (T−1) and a Gibbs sampler algorithm that obtains posterior estimates of all the unknowns in the model. The new procedure was used with two simulated data sets and with a beef cattle database. In the simulated populations, the results of the analysis provided marginal posterior distributions that included the population parameters in the regions of highest posterior density. In the case of the beef cattle dataset, the posterior estimate of transgenerational epigenetic variability was very low and a model comparison test indicated that a model that did not included it was the most plausible.

Genetic diversity and divergence among Spanish beef cattle breeds assessed by a bovine high-density SNP chip

The availability of SNP chips for massive genotyping has proven to be useful to genetically characterize populations of domestic cattle and to assess their degree of divergence. In this study, the Illumina BovineHD BeadChip genotyping array was used to describe the genetic variability and divergence among 7 important autochthonous Spanish beef cattle breeds. The within-breed genetic diversity, measured as the marker expected heterozygosity, was around 0.30, similar to other European cattle breeds. The analysis of molecular variance revealed that 94.22% of the total variance was explained by differences within individuals whereas only 4.46% was the result of differences among populations. The degree of genetic differentiation was small to moderate as the pairwise fixation index of genetic differentiation among breeds (F) estimates ranged from 0.026 to 0.068 and the Neis D genetic distances ranged from 0.009 to 0.016. A neighbor joining (N-J) phylogenetic tree showed 2 main groups of breeds: Pirenaica, Bruna dels Pirineus, and Rubia Gallega on the one hand and Avileña-Negra Ibérica, Morucha, and Retinta on the other. In turn, Asturiana de los Valles occupied an independent and intermediate position. A principal component analysis (PCA) applied to a distance matrix based on marker identity by state, in which the first 2 axes explained up to 17.3% of the variance, showed a grouping of animals that was similar to the one observed in the N-J tree. Finally, a cluster analysis for ancestries allowed assigning all the individuals to the breed they belong to, although it revealed some degree of admixture among breeds. Our results indicate large within-breed diversity and a low degree of divergence among the autochthonous Spanish beef cattle breeds studied. Both N-J and PCA groupings fit quite well to the ancestral trunks from which the Spanish beef cattle breeds were supposed to derive.

On the performance of tests for the detection of signatures of selection: a case study with the Spanish autochthonous beef cattle populations

Background Procedures for the detection of signatures of selection can be classified according to the source of information they use to reject the null hypothesis of absence of selection. Three main groups of tests can be identified that are based on: (1) the analysis of the site frequency spectrum, (2) the study of the extension of the linkage disequilibrium across the length of the haplotypes that surround the polymorphism, and (3) the differentiation among populations. The aim of this study was to compare the performance of a subset of these procedures by using a dataset on seven Spanish autochthonous beef cattle populations.ResultsAnalysis of the correlations between the logarithms of the statistics that were obtained by 11 tests for detecting signatures of selection at each single nucleotide polymorphism confirmed that they can be clustered into the three main groups mentioned above. A factor analysis summarized the results of the 11 tests into three canonical axes that were each associated with one of the three groups. Moreover, the signatures of selection identified with the first and second groups of tests were shared across populations, whereas those with the third group were more breed-specific. Nevertheless, an enrichment analysis identified the metabolic pathways that were associated with each group; they coincided with canonical axes and were related to immune response, muscle development, protein biosynthesis, skin and pigmentation, glucose metabolism, fat metabolism, embryogenesis and morphology, heart and uterine metabolism, regulation of the hypothalamic–pituitary–thyroid axis, hormonal, cellular cycle, cell signaling and extracellular receptors.ConclusionsWe show that the results of the procedures used to identify signals of selection differed substantially between the three groups of tests. However, they can be classified using a factor analysis. Moreover, each canonical factor that coincided with a group of tests identified different signals of selection, which could be attributed to processes of selection that occurred at different evolutionary times. Nevertheless, the metabolic pathways that were associated with each group of tests were similar, which suggests that the selection events that occurred during the evolutionary history of the populations probably affected the same group of traits.

Genomic differentiation between Asturiana de los Valles, Avileña-Negra Ibérica, Bruna dels Pirineus, Morucha, Pirenaica, Retinta and Rubia Gallega cattle breeds

The Spanish local beef cattle breeds have most likely common origin followed by a process of differentiation. This particular historical evolution has most probably left detectable signatures in the genome. The objective of this study was to identify genomic regions associated with differentiation processes in seven Spanish autochthonous populations (Asturiana de los Valles (AV), Avileña-Negra Ibérica (ANI), Bruna dels Pirineus (BP), Morucha (Mo), Pirenaica (Pi), Retinta (Re) and Rubia Gallega (RG)). The BovineHD 777K BeadChip was used on 342 individuals (AV, n=50; ANI, n=48; BP, n=50; Mo, n=50; Pi, n=48; Re, n=48; RG, n=48) chosen to be as unrelated as possible. We calculated the fixation index (F ST ) and performed a Bayesian analysis named SelEstim. The output of both procedures was very similar, although the Bayesian analysis provided a richer inference and allowed us to calculate significance thresholds by generating a pseudo-observed data set from the estimated posterior distributions. We identified a very large number of genomic regions, but when a very restrictive significance threshold was applied these regions were reduced to only 10. Among them, four regions can be highlighted because they comprised a large number of single nucleotide polymorphisms and showed extremely high signals (Kullback-Leiber divergence (KLD)>6). They are located in BTA 2 (5 575 950 to 10 152 228 base pairs (bp)), BTA 5 (17 596 734 to 18 850 702 bp), BTA 6 (37 853 912 to 39 441 548 bp) and BTA 18 (13 345 515 to 15 243 838 bp) and harbor, among others, the MSTN (Myostatin), KIT-LG (KIT Ligand), LAP3 (leucine aminopeptidase 3), NAPCG (non-SMC condensing I complex, subunit G), LCORL (ligand dependent nuclear receptor corepressor-like) and MC1R (Melanocortin 1 receptor) genes. Knowledge on these genomic regions allows to identify potential targets of recent selection and helps to define potential candidate genes associated with traits of interest, such as coat color, muscle development, fertility, growth, carcass and immunological response.

Performance of genomic selection under a single-step approach in autochthonous Spanish beef cattle populations

This study evaluated different strategies for implementing a single-step genomic selection programme in two autochthonous Spanish beef cattle populations (Pirenaica-Pi and Rubia Gallega-RG). The strategies were compared in terms of accuracy attained under different scenarios by simulating genomic data over the known genealogy. Several genotyping approaches were tested, as well as, other factors like marker density, effective population size, mutation rate and heritability of the trait. The results obtained showed gains in accuracy with respect to pedigree BLUP evaluation in all cases. The greatest benefit was obtained when the candidates to selection had their genotypes included in the evaluation. Moreover, genotyping the individuals with the most accurate predictions maximized the gains but other suboptimal strategies also yielded satisfactory results. Furthermore, the gains in accuracy increased with the marker density reaching a plateau at around 50,000 markers. Likewise, the effective population size and the mutation rate have also shown an effect, both increasing the accuracy with decreasing values of these population parameters. Finally, the results obtained for the RG population showed greater gains compared to the Pi population, probably attributed to the wider implantation of artificial insemination.

Non-linear recursive models for growth traits in the Pirenaica beef cattle breed.

One of the main goals of selection schemes in beef cattle populations is to increase carcass weight at slaughter. Live weights at different growth stages are frequently used as selection criteria under the hypothesis that they usually have a high and positive genetic correlation with weight at slaughter. However, the presence of compensatory growth may bias the prediction ability of early weights for selection purposes. Recursive models may represent an interesting alternative for understanding the genetic and phenotypic relationship between weight traits during growth. For the purposes of this study, the analysis was performed for three different set of data from the Pirenaica beef cattle breed: weight at 120 days (W120) and at 210 days (W210); W120 and carcass weight at slaughter at 365 days (CW365); W210 and CW365. The number of records for each analysis was 8592, 4648 and 3234, respectively. A pedigree composed of 56323 individuals was also included. The statistical model comprised sex, year-season of birth, herd and slaughterhouse, plus a non-linear recursive dependency between traits. The dependency was modeled as a polynomial up to the 4th degree and models were compared using a Logarithm of Conditional Predictive Ordinates. The results of model comparison suggest that the best models were the 3rd degree polynomial for W120-W210 and W120-CW365 and the 2nd degree polynomial for W210-CW365. The posterior mean estimates for heritabilities ranged between 0.29 and 0.44 and the posterior mean estimates of the genetic correlations were null or very low, indicating that the relationship between traits is fully captured by the recursive dependency. The results imply that the predictive ability of the performance of future growth is low if it is only based on records of early weights. The usefulness of slaughterhouse records in beef cattle breeding evaluation is confirmed.