Genetic analysis of the effects of heat stress before and after lambing on pre‐weaning live weight in Spanish Merino lambs

Background

Heat stress (HS) is a major environmental effect on sheep production. Hereby, we estimated the genetic (co)variance component of HS on the pre‐weaning performance of 19,022 Merino lambs by analysing the climatological index of temperature and relative humidity (recorded 30 days before lambing and after lambing) using transversal and longitudinal mixed linear models.

Methods and Results

The global impact of HS during the last 30 days of pregnancy was −17% for birthweight and ranged between −4% and −8% for live weight at 15, 30 days of age (W30), and average daily gain from birth at 30 days. The results from both statistical approaches showed very similar heritabilities (h²), ranging from 0.192 to 0.237 for the direct genetic (D) effects and from 0.072 to 0.082 for the maternal genetic (M) effects, but the antagonism between (D) and (M) was higher when a longitudinal model was used. A significant genotype‐environmental effect was also found regardless of whether the climatological covariables were considered in the model. In addition, we employed D and M breeding values for W30 as an example to create a new subjacent index by first using a principal component analysis and employing the leading eigenvalues as a weighted factor that provides the information needed to identify those genotypes that maximise the response for both genetic effects over a wide range of climate–environment levels.

Conclusions

Our study revealed that the HS indexes of the mother during the gestation period have a significant effect on the growth of the lambs during the early stages of life.

BREEDING AND GENETICS SYMPOSIUM: Breeding for resilience to heat stress effects in dairy ruminants. A comprehensive review

Selection for heat tolerant (HT) animals in dairy production has been so far linked to estimation of declines in production using milk recording and meteorological information on the day of control using reaction norms. Results from these models show that there is a reasonable amount of genetic variability in the individual response to high heat loads, which makes feasible selection of HT animals at low costs. However, the antagonistic relationship between level of production and response to heat stress (HS) implies that selection for HT animals under this approach must be done with caution so that productivity is not damaged. Decomposition of the genetic variability in principal components (PC) can provide selection criteria independent of milk production level although biological interpretation of PC is difficult. Moreover, given that response to heat stress for each animal is estimated with very sparse information collected under different physiological and management circumstances, biased (normally underestimation) and lack of accuracy may be expected. Alternative phenotypic characterization of HT can come from the use of physiological traits, which have also shown moderate heritability. However, costs of a large scale implementation based on physiological characteristics has precluded its use. Another alternative is the use of biomarkers that define heat tolerance. A review of biomarkers of HS from more recent studies is provided. Of particular interest are milk biomarkers, which together with infrared spectra prediction equations can provide useful tools for genetic selection. In the 'omics' era, genomics, transcriptomics, proteomics and metabolomics have been already used to detect genes affecting HT. A review of findings in these areas is also provided. Except for the slick hair gene, there are no other genes for which variants have been clearly associated with HT. However, integration of omics information could help in pointing at knots of the HS control network and, in the end, to a panel of markers to be used in the selection of HT animals. Overall, HT is a complex phenomenon that requires integration of fine phenotypes and omics information to provide accurate tools for selection without damaging productivity. Technological developments to make on-farm implementation feasible and with greater insight into the key biomarkers and genes involved in HT are needed.

Short-term glutamate administration positively affects the number of antral follicles and the ovulation rate in cyclic adult goats

The acute effects of short-term glutamate administration on the number of antral follicles and ovulation rate were examined in adult goats. Neither live weight (44.5±1.3 kg) nor body condition (3.3±0.8 units) differed between the control (untreated) and glutamate-treated (0.175 mg/kg) animals (p>0.05). However, the number of antral follicles (3.4 vs. 2.1, p=0.05) and ovulation rate (2.5 vs. 1.5, p=0.05) was higher in the glutamate-administered group than in the controls.

Glutamate supply positively affects serum cholesterol concentrations without increases in total protein and urea around the onset of puberty in goats

Different neurotransmitter and neuromodulatory systems regulate synthesis and secretion of GnRH. Whereas the endocrine and neural systems are activated in response to the metabolic status and the circulating levels of specific blood metabolites, glutamate receptors have been reported at hepatic level. This study evaluated the possible effect of glutamate supplementation upon changes in serum concentrations across time for total protein (TP), urea (UR) and cholesterol (CL) around the onset of puberty in goats. Prepuberal female goats (n=18) were randomly assigned to: (1) excitatory amino acids group, GLUT, n=10; 16.52±1.04kg live weight (LW), 3.4±0.12 body condition score (BCS) receiving an i.v. infusion of 7mgkg(-1) LW of l-glutamate, and (2) Control group, CONT, n=8; 16.1±1.04kg LW, 3.1±0.12 BCS. General averages for LW (23.2±0.72kg), BCS (3.37±0.10 units), serum TP (65.28±2.46mgdL(-1)), UR (23.42±0.95mgdL(-1)), CL (77.89±1.10mgdL(-1)) as well as the serum levels for TP and UR across time did not differ (P>0.05) between treatments. However, while GLUT positively affected (P<0.05) both the onset (207±9 vs. 225±12 d) and the percentage (70 vs. 25%) of females showing puberty, a treatment×time interaction effect (P<0.05) was observed in the GLUT group, with increases in serum cholesterol, coincident with the onset of puberty. Therefore, in peripuberal glutamate supplemented goats, serum cholesterol profile could act as a metabolic modulator for the establishment of puberty, denoting also a potential role of glutamate as modulator of lipid metabolism.

Identification of two paralogous caprine CD36 genes that display highly divergent mRNA expression profiles

The CD36 molecule plays a pivotal role in a variety of immunological and cellular processes, including pathogen recognition, inflammation and apoptosis. Herein, we demonstrate that this gene is duplicated in goats, with two copies (CD36 and CD36-like) that display highly divergent mRNA expression profiles. In this way, CD36 mRNA is mostly expressed in the adipose tissue and heart whilst CD36-like mRNA shows a high expression in the liver. We have also found evidence of the presence of two paralogous CD36 and CD36-like genes in the bovine genome, suggesting that CD36 duplication took place before goat-cattle radiation, i.e. at least 20 MYR ago. Finally, we have characterized the polymorphism of the coding regions of the goat CD36 and CD36-like genes. In doing so, we have identified one synonymous polymorphism at the CD36-like gene (c.390A>C) that displays a significant association (P=0.04) with milk somatic cell count, a parameter often used to diagnose mastitis in domestic ruminants.

Genetic variation of adaptation to heat stress in two Spanish dairy goat breeds

Monthly test-day records of milk yield and composition in Murciano-Granadina (MG) and Payoya (PYA) dairy goats were combined with weather data from meteorological stations, to analyse the effects of heat stress on dairy traits, measured with an index of temperature and relative humidity (THI). A 'repeatability animal model' and a 'reaction norm animal model' were used to estimate genetic (co)variance components. Estimates obtained with both models were very similar. The h(2) of daily yields in MG did not vary throughout the THI scale, but the pattern of variation of content traits showed negative trends for increasing THI values. In PYA, a slight positive tendency throughout the THI scale was observed for the same traits. The genetic correlations between subsequent points in the THI scale were lower than 0.80 when they were computed between low and high THI points. The same reaction norm was observed for all traits. Using the 'reaction norm animal model', it was possible to identify those animals that show the same performance throughout the THI trajectory (robust) and those with varying performances (tolerant and non-tolerant to heat stress). Results in this study also show that heat tolerance decreases, while the genetic level for milk traits increases. Losses because of heat stress were equivalent to 1.9 and 3.1% in the yearly yield of fat and protein for MG and PYA, respectively.

Pleiotropic effects of the goat prolactin receptor genotype on milk fatty acid composition

In the lactating mammary gland, prolactin (PRL) stimulates the synthesis of lactose as well as fatty acid uptake, lipogenesis, and triacylglycerol synthesis. Associations between bovine PRL receptor (PRLR) genotype and fat yield have been reported, which illustrates the role of PRL in conveying lipids toward the udder as well as in stimulating their local synthesis during lactation. Conversely, and to the best of our knowledge, the effects of PRLR genotype on milk fatty acid content have not been studied so far in any mammalian species. In this study, we sequenced most of the coding region of the caprine PRLR gene in several individuals from the Malagueña and Murciano-Granadina breeds. This approach allowed us to identify 2 long and short mRNA isoforms, produced by alternative splicing, and 4 single-nucleotide polymorphisms (SNPs), namely, c.177T>C, c.1131G>A, c.1201G>A and c.1355C>T. Two of these SNPs are nonsynonymous and involve G401R (c.1201G>A) and T452I (c.1355C>T) amino acid substitutions in the cytoplasmic domain of PRLR, which plays a fundamental role in signal transduction. Performance of an association analysis with milk composition traits in a Murciano-Granadina goat population revealed highly suggestive effects on palmitoleic acid content, whereas suggestive effects were detected for other fatty acids, such as palmitic and linoleic. These results are consistent with the pleiotropic effects of PRL on mammary gland lipid metabolism and milk composition.

Alloprocoptic selection: A mode of natural selection promoting polymorphism

The fitness effect of genetic variation at three loci coding for enzymes is studied in Drosophila melanogaster. The fertility of a female is determined by the female genotype as well as by the genotype of the male with which she mates. Significant interactions exist between female and male genotype, so that the fertility of a given mating combination cannot be predicted from the average fertility of the two genotypes involved. Multiple stable equilibria are possible when such interactions exist. At two loci, the fertility is greater than expected when the two mating individuals are homozygous for different alleles and smaller than expected when they are homozygous for the same allele. This mode of selection in which association of opposites increases their fitness is herein named alloprocoptic selection. It will contribute to maintaining genetic polymorphism in nature.

Short communication: genetic variability in the predicted microRNA target sites of caprine casein genes

The main goal of the current work was to identify single nucleotide polymorphisms (SNP) that might create or disrupt microRNA (miRNA) target sites in the caprine casein genes. The 3' untranslated regions of the goat alpha(S1)-, alpha(S2)-, beta-, and kappa-casein genes (CSN1S1, CSN1S2, CSN2, and CSN3, respectively) were resequenced in 25 individuals of the Murciano-Granadina, Cashmere, Canarian, Saanen, and Sahelian breeds. Five SNP were identified through this strategy: c.175C>T at CSN1S1; c.109T>C, c.139G>C, and c.160T>C at CSN1S2; and c.216C>T at CSN2. Analysis with the Patrocles Finder tool predicted that all of these SNP are located within regions complementary to the seed of diverse miRNA sequences. These in silico results suggest that polymorphism at miRNA target sites might have some effect on casein expression. We explored this issue by genotyping the c.175C>T SNP (CSN1S1) in 85 Murciano-Granadina goats with records for milk CSN1S1 concentrations. This substitution destroys a putative target site for miR-101, a miRNA known to be expressed in the bovine mammary gland. Although TT goats had higher levels (6.25 g/L) of CSN1S1 than their CT (6.05 g/L) and CC (6.04 g/L) counterparts, these differences were not significant. Experimental confirmation of the miRNA target sites predicted in the current work and performance of additional association analyses in other goat populations will be an essential step to find out if polymorphic miRNA target sites constitute an important source of variation in casein expression.

Random regression model of growth during the first three months of age in Spanish Merino sheep1,2

A total of 88,727 individual BW records of Spanish Merino lambs, obtained from 30,214 animals between 2 and 92 d of age, were analyzed using a random regression model (RRM). These animals were progeny of 546 rams and 15,586 ewes raised in 30 flocks, between 1992 and 2002, with a total of 45,941 animals in the pedigree. The contemporary groups (animals of the same flock, year, and season, with 452 levels), the lambing number (11 levels), the combination sex of lambs with type of litter (4 levels), and a fixed regression coefficient of age on BW were included as fixed effects. A total of 7 RRM were compared, and the best fit was obtained for a model of order 3 for the direct and maternal genetic effects and for the individual permanent environmental effect. For the maternal permanent environmental effect the best model had an order 2. The residual variance was assumed to be heterogeneous with 10 age classes; the covariance between both genetic effects was included. According to the results of the selected RRM, the heritability for both genetic effects (h(a)2 and h(m)2) increased with age, with estimates of 0.123 to 0.186 for h(a)2 and of 0.059 to 0.108 for h(m)2. The correlations between direct and genetic maternal effects were -0.619 to -0.387 during the first 45 d of age and decreased as age increased, until reaching values from -0.366 to -0.275 between 45 to 75 d of age. Important changes in ranking of the animals were found based on the breeding value estimation with the current method and with the random regression procedure. The use of RRM to analyze the genetic trajectory of growth in this population of Merino sheep is highly recommended.

Short communication: Identification of two polymorphisms in the goat lipoprotein lipase gene and their association with milk production traits

Lipoprotein lipase (LPL) is a glycoprotein that plays a central role in plasma triglyceride metabolism by hydrolyzing triglyceride-rich chylomicrons and very low density lipoproteins. The activity of milk LPL has been shown to differ among several goat breeds, suggesting the existence of a genetic polymorphism influencing the functional properties of this enzyme. We have characterized the complete coding sequence of the goat LPL gene in 18 individuals belonging to 3 breeds. The coding region of the goat LPL cDNA was 1,437 bp long and encoded a protein of 478 amino acids. Moreover, we have identified 2 single nucleotide polymorphisms (SNP) including a G50C missense mutation, which involved a Ser-->Thr amino acid replacement at position 17 of the signal peptide, and a C2094T substitution in the 3' untranslated region. A univariate mixed model was used to evaluate the association between LPL genotypes and milk production and composition in 130 Murciano-Granadina goats. The G50C SNP was suggestively associated with milk fat content and tended to affect the milk dry weight basis. The C2094T SNP was not associated with any of the measured traits.

Goat acetyl-coenzyme A carboxylase alpha: molecular characterization, polymorphism, and association with milk traits

Acetyl-coenzyme A carboxylase alpha (ACACA), the major regulatory enzyme of fatty acid biosynthesis, catalyzes the conversion of acetyl-CoA to malonyl-CoA. We have sequenced 5.5 kb of the goat ACACA cDNA in 18 individuals. The translated cDNA sequence encodes 1,832 AA and shares a high AA identity (99%) with the ovine and bovine ACACA orthologous sequences. One silent single nucleotide polymorphism was identified at exon 45 (C5493T). The genotyping of this polymorphism in 4 goat breeds, Murciano-Granadina, Teramana, Majorera, and Malagueña, showed that it was segregating in all 4 breeds at variable frequencies. In addition, a mixed animal model was used to evaluate the association of C5493T genotypes with milk traits. Results from this study show that the C5493T mutation was suggestively associated with fat yield, lactose content, and somatic cell count.

Potential benefit from using the alphas1-casein genotype information in a selection scheme for dairy goats

A stochastic approach is proposed to predict responses to selection when using alpha(s1)-casein genotype information in a selection scheme of a Spanish breed of dairy goats. Two independent selection objectives were considered: protein yield (PY), where the major additive gene CSN1S1, which codes for alpha(s1)-casein, has a small effect, and protein content (P%), where this gene has a large effect on performances. Significant differences in response between using and ignoring information on the major gene were observed only when the major gene has a large effect. The main result was in the case of P%, the total genetic gain obtained in the early generations of selection was maintained in the long-term. Taking account of genotype information either in the evaluation model or in the selection criteria leads to a faster fixation of the favourable allele and a reduction of the total genetic variance over generations. The inbreeding rates varied across generations, the highest rates observed in later generations of selection and when the major gene has a large effect and its genotype was included in the genetic evaluation procedure. It is concluded that inclusion of the casein genotype as an additional selection criteria will improve gains for protein traits, in particular P%. Recommendations are also given in order to optimize the use of this molecular information in dairy goat selection programs.

Estimation of genetic parameters for milk yield and fat and protein contents of milk from murciano-granadina goats

A total of 10,289 lactations, corresponding to 4967 does distributed in 35 herds, was used to estimate heritabilities, repeatabilities, and correlations among milk yield (kilograms), fat contents, and protein contents (grams/100 g of milk) with single-trait and multiple-trait animal models using repeated records. Single-trait estimates were similar to those obtained with multiple-trait analysis. Heritabilities were 0.18, 0.16, and 0.25 with single-trait analysis and 0.17, 0.14, and 0.22 with multiple-trait analysis for milk yield, fat content, and protein content, respectively. Repeatabilities were 0.39, 0.36, and 0.47 with single-trait analysis and 0.36, 0.33, and 0.41 with multiple-trait analysis for milk yield, fat content, and protein content, respectively. Genetic correlations were -0.89 between milk yield and fat content, -0.65 between milk yield and protein content, and 0.93 between fat and protein contents. Phenotypic correlations were -0.48 between milk yield and fat content, -0.47 between milk yield and protein content, and 0.54 between fat and protein contents. Heritabilities and repeatabilities were lower, and correlations were higher, than those found in the literature. Some combination of both EBV of milk yield and protein content would be a good criterion for selection to improve milk yield while maintaining fat and protein contents.

Boundary-layer model for the population dynamics of single species

We develop a new discrete-time model, called the boundary-layer model, to describe the dynamics of single species that have a capacity for fast growth at very low population densities. The model explicitly separates the dynamics of the population at very low densities (within the "boundary layer") and at high densities. The boundary-layer model provides a better fit than other models such as the logistic or the theta model to data from experimental populations of Drosophila willistoni and D. pseudoobscura.