Genome-wide association study of age at puberty and its (co)variances with fertility and stature in growing and lactating Holstein-Friesian dairy cattle.

Reproductive performance is a key determinant of cow longevity in a pasture-based, seasonal dairy system. Unfortunately, direct fertility phenotypes such as intercalving interval or pregnancy rate tend to have low heritabilities and occur relatively late in an animal's life. In contrast, age at puberty (AGEP) is a moderately heritable, early-in-life trait that may be estimated using an animal's age at first measured elevation in blood plasma progesterone (AGEP4) concentrations. Understanding the genetic architecture of AGEP4 in addition to genetic relationships between AGEP4 and fertility traits in lactating cows is important, as is its relationship with body size in the growing animal. Thus, the objectives of this research were 3-fold. First, to estimate the genetic and phenotypic (co)variances between AGEP4 and subsequent fertility during first and second lactations. Second, to quantify the associations between AGEP4 and height, length, and BW measured when animals were approximately 11 mo old (standard deviation = 0.5). Third, to identify genomic regions that are likely to be associated with variation in AGEP4. We measured AGEP4, height, length, and BW in approximately 5,000 Holstein-Friesian or Holstein-Friesian × Jersey crossbred yearling heifers across 54 pasture-based herds managed in seasonal calving farm systems. We also obtained calving rate (CR42, success or failure to calve within the first 42 d of the seasonal calving period), breeding rate (PB21, success or failure to be presented for breeding within the first 21 d of the seasonal breeding period) and pregnancy rate (PR42, success or failure to become pregnant within the first 42 d of the seasonal breeding period) phenotypes from their first and second lactations. The animals were genotyped using the Weatherby's Versa 50K SNP array (Illumina, San Diego, CA). The estimated heritabilities of AGEP4, height, length, and BW were 0.34 (90% credibility interval [CRI]: 0.30, 0.37), 0.28 (90% CRI: 0.25, 0.31), 0.21 (90% CRI: 0.18, 0.23), and 0.33 (90% CRI: 0.30, 0.36), respectively. In contrast, the heritabilities of CR42, PB21 and PR42 were all <0.05 in both first and second lactations. The genetic correlations between AGEP4 and these fertility traits were generally moderate, ranging from 0.11 to 0.60, whereas genetic correlations between AGEP4 and yearling body-conformation traits ranged from 0.02 to 0.28. Our GWAS highlighted a genomic window on chromosome 5 that was strongly associated with variation in AGEP4. We also identified 4 regions, located on chromosomes 14, 6, 1, and 11 (in order of decreasing importance), that exhibited suggestive associations with AGEP4. Our results show that AGEP4 is a reasonable predictor of estimated breeding values for fertility traits in lactating cows. Although the GWAS provided insights into genetic mechanisms underpinning AGEP4, further work is required to test genomic predictions of fertility that use this information.

Genome-wide association study of anogenital distance and its (co)variances with fertility in growing and lactating Holstein-Friesian dairy cattle.

Anogenital distance (AGD) is a moderately heritable trait that can be measured at a young age that may provide an opportunity to indirectly select for improved fertility in dairy cattle. In this study, we characterized AGD and its genetic and phenotypic relationships with a range of body stature and fertility traits. We measured AGD, shoulder height, body length, and body weight in a population of 5,010 Holstein-Friesian and Holstein-Friesian × Jersey crossbred heifers at approximately 11 mo of age (AGD1). These animals were born in 2018 across 54 seasonal calving, pasture-based dairy herds. A second measure of AGD was collected in a subset of herds (n = 17; 1,956 animals) when the animals averaged 29 mo of age (AGD2). Fertility measures included age at puberty (AGEP), then time of calving, breeding, and pregnancy during the first and second lactations. We constructed binary traits reflecting the animal's ability to calve during the first 42 d of their herd's seasonal calving period (CR42), be presented for breeding during the first 21 d of the seasonal breeding period (PB21) and become pregnant during the first 42 d of the seasonal breeding period (PR42). The posterior mean of sampled heritabilities for AGD1 was 0.23, with 90% of samples falling within a credibility interval (90% CRI) of 0.20 to 0.26, whereas the heritability of AGD2 was 0.29 (90% CRI 0.24 to 0.34). The relationship between AGD1 and AGD2 was highly positive, with a genetic correlation of 0.89 (90% CRI 0.82 to 0.94). Using a GWAS analysis of 2,460 genomic windows based on 50k genotype data, we detected a region on chromosome 20 that was highly associated with variation in AGD1, and a second region on chromosome 13 that was moderately associated with variation in AGD1. We did not detect any genomic regions associated with AGD2 which was measured in fewer animals. The genetic correlation between AGD1 and AGEP was 0.10 (90% CRI 0.00 to 0.19), whereas the genetic correlation between AGD2 and AGEP was 0.30 (90% CRI 0.15 to 0.44). The timing of calving, breeding, and pregnancy (CR42, PB21, and PR42) during first or second lactations exhibited moderate genetic relationships with AGD1 (0.19 to 0.52) and AGD2 (0.46 to 0.63). Genetic correlations between AGD and body stature traits were weak (≤0.16). We conclude that AGD is a moderately heritable trait, which may have value as an early-in-life genetic predictor for reproductive success during lactation.

Variance parameter estimation for age at puberty phenotypes under 2 levels of phenotype censorship.

Age at puberty (AGEP) is a moderately heritable trait in cattle that may be predictive of an animal's genetic merit for reproductive success later in life. In addition, under some mating strategies (for example, where mating begins before all animals have attained puberty) animals that attain puberty at a relatively young age will also likely conceive earlier than their herd mates, and thus begin their productive life earlier. Unfortunately, AGEP is challenging to measure because animals must be observed over a period of several months. Our objectives for this study were twofold. The first objective was to produce variance components for AGEP. The second objective was to investigate the implications of a simplified phenotyping strategy for AGEP, when the interval between repeated blood plasma progesterone measures was extended from weekly to monthly, increasing the extent of left, interval, and right censoring. We measured AGEP in a closely monitored population of around 500 Holstein-Friesian heifers, born in 2015 and managed under a seasonal, pasture-based dairy system. Animals were blood tested weekly from approximately 240 to 440 d of age and were deemed to have reached puberty when blood plasma progesterone elevation (>1 ng/mL) was detected in 2 of 3 consecutive blood tests (AGEP_Weekly). To simulate a simplified phenotyping strategy based on monthly herd visits (AGEP_Monthly), we selectively disregarded data from all but 3 blood test events, when animals were around 300, 330, and 360 d of age (standard deviation = 14.5 d). The posterior mean of estimated heritabilities for AGEP_Weekly was 0.54, with a 90% credibility interval (90% CRI) of 0.41 to 0.66, whereas it was 0.44 (90% CRI 0.32 to 0.57) for AGEP_Monthly. The correlation between EBVs for AGEP_Weekly and AGEP_Monthly was 0.87 (90% CRI, 0.84 to 0.89). We conclude that in this population, AGEP is a moderately heritable trait. Further, increasing phenotype censorship from weekly to monthly observations would not have altered the main conclusions of this analysis. Our results support the strategic use of censoring to reduce costs and animal ethics considerations associated with collection of puberty phenotypes.

Erratum to "Cumulative dairy cow genetic change from selection and crossbreeding over the last 2 decades in New Zealand closely aligns to model-based predictions published in 2000" (JDS Commun. 2:51-54).

[This corrects the article DOI: 10.3168/jdsc.2020-0043.].

Pregnancy status predicted using milk mid-infrared spectra from dairy cattle.

Accurate and timely pregnancy diagnosis is an important component of effective herd management in dairy cattle. Predicting pregnancy from Fourier-transform mid-infrared (FT-MIR) spectroscopy data is of particular interest because the data are often already available from routine milk testing. The purpose of this study was to evaluate how well pregnancy status could be predicted in a large data set of 1,161,436 FT-MIR milk spectra records from 863,982 mixed-breed pasture-based New Zealand dairy cattle managed within seasonal calving systems. Three strategies were assessed for defining the nonpregnant cows when partitioning the records according to pregnancy status in the training population. Two of these used records for cows with a subsequent calving only, whereas the third also included records for cows without a subsequent calving. For each partitioning strategy, partial least squares discriminant analysis models were developed, whereby spectra from all the cows in 80% of herds were used to train the models, and predictions on cows in the remaining herds were used for validation. A separate data set was also used as a secondary validation, whereby pregnancy diagnosis had been assigned according to the presence of pregnancy-associated glycoproteins (PAG) in the milk samples. We examined different ways of accounting for stage of lactation in the prediction models, either by including it as an effect in the prediction model, or by pre-adjusting spectra before fitting the model. For a subset of strategies, we also assessed prediction accuracies from deep learning approaches, utilizing either the raw spectra or images of spectra. Across all strategies, prediction accuracies were highest for models using the unadjusted spectra as model predictors. Strategies for cows with a subsequent calving performed well in herd-independent validation with sensitivities above 0.79, specificities above 0.91 and area under the receiver operating characteristic curve (AUC) values over 0.91. However, for these strategies, the specificity to predict nonpregnant cows in the external PAG data set was poor (0.002-0.04). The best performing models were those that included records for cows without a subsequent calving, and used unadjusted spectra and days in milk as predictors, with consistent results observed across the training, herd-independent validation and PAG data sets. For the partial least squares discriminant analysis model, sensitivity was 0.71, specificity was 0.54 and AUC values were 0.68 in the PAG data set; and for an image-based deep learning model, the sensitivity was 0.74, specificity was 0.52 and the AUC value was 0.69. Our results demonstrate that in pasture-based seasonal calving herds, confounding between pregnancy status and spectral changes associated with stage of lactation can inflate prediction accuracies. When the effect of this confounding was reduced, prediction accuracies were not sufficiently high enough to use as a sole indicator of pregnancy status.

Animal medical genetics: a historical perspective on more than 50 years of research into genetic disorders of animals at Massey University.

Over the last 50 years, there have been major advances in knowledge and technology regarding genetic diseases, and the subsequent ability to control them in a cost-effective manner. This review traces these advances through research into genetic diseases of animals at Massey University (Palmerston North, NZ), and briefly discusses the disorders investigated during that time, with additional detail for disorders of major importance such as bovine α-mannosidosis, ovine ceroid-lipofuscinosis, canine mucopolysaccharidosis IIIA and feline hyperchylomicronaemia. The overall research has made a significant contribution to veterinary medicine, has provided new biological knowledge and advanced our understanding of similar disorders in human patients, including testing various specific therapies prior to human clinical trials.

Cumulative dairy cow genetic change from selection and crossbreeding over the last 2 decades in New Zealand closely aligns to model-based predictions published in 2000.

A deterministic model was developed in 1998 to evaluate the concurrent effects of selection and crossbreeding on the rate of genetic gain and increases in productivity of New Zealand dairy cattle over the ensuing 25-yr period. The predictions of today's breed composition of the national dairy herd and genetic trends for body weight and lactation yields of milk, fat, and protein are compared with today's actual values. Selection was assumed to use an index that included live weight and lactation yields of milk, fat, and protein. Mating strategies involving the Holstein-Friesian (F), Jersey (J), and Ayrshire (A) breeds were evaluated. Effects of heterosis and age were included to calculate phenotypic live weight and yields of milk, fat, and protein per cow. At the time the model was developed, New Zealand had an across-breed evaluation system, but only straightbred bulls were used after progeny testing based on records of straightbred and crossbred daughters. The model predicted that if crossbred cows and bulls could be considered as bull parents, faster rates of genetic gain may result because of increased selection intensities in the cow to breed bull selection pathway. This scenario transpired, and the best bulls and cows for farm profit were used regardless of breed. Under that mating strategy for the 2018 birthyear, the model predicted the national breed composition would be 11% F, 34% J, 52% F×J, and 2% A; the actual breed composition was 36% F, 9% J, 53% F×J, and 1% A. The model-predicted annual genetic gains would be 16.7 L of milk, 1.2 kg of fat, 1.5 kg of protein, and -0.7 kg of body weight; the realized annual improvements were 13.6 L of milk, 1.31 kg of fat, 1.17 kg of protein, and -0.36 kg of body weight. Predicted long-term responses to selection can closely mirror realized improvements, confirming the value of modeling to inform animal breeding decision-making.

Accuracy of genomic prediction of shell quality in a White Leghorn line.

Several genomic methods were applied for predicting shell quality traits recorded at 4 different hen ages in a White Leghorn line. The accuracies of genomic prediction of single-step GBLUP and single-trait Bayes B were compared with predictions of breeding values based on pedigree-BLUP under single-trait or multitrait models. Breaking strength (BS) and dynamic stiffness (Kdyn) measurements were collected on 18,524 birds from 3 consecutive generations, of which 4,164 animals also had genotypes from an Affymetrix 50K panel containing 49,591 SNPs after quality control edits. All traits had low to moderate heritability, ranging from 0.17 for BS to 0.34 for Kdyn. The highest accuracies of prediction were obtained for the multitrait single-step model. The use of marker information resulted in higher prediction accuracies than pedigree-based models for almost all traits. A genome-wide association study based on a Bayes B model was conducted to detect regions explaining the largest proportion of genetic variance. Across all 8 shell quality traits analyzed, 7 regions each explaining over 2% of genetic variance and 54 regions each explaining over 1% of genetic variance were identified. The windows explaining a large proportion of genetic variance overlapped with several potential candidate genes with biological functions linked to shell formation. A multitrait repeatability model using a single-step method is recommended for genomic evaluation of shell quality in layer chickens.

The evolving role of Fourier-transform mid-infrared spectroscopy in genetic improvement of dairy cattle.

Over the last 100 years, significant advances have been made in the characterisation of milk composition for dairy cattle improvement programs. Technological progress has enabled a shift from labour intensive, on-farm collection and processing of samples that assess yield and fat levels in milk, to large-scale processing of samples through centralised laboratories, with the scope extended to include quantification of other traits. Fourier-transform mid-infrared (FT-MIR) spectroscopy has had a significant role in the transformation of milk composition phenotyping, with spectral-based predictions of major milk components already being widely used in milk payment and animal evaluation systems globally. Increasingly, there is interest in analysing the individual FT-MIR wavenumbers, and in utilising the FT-MIR data to predict other novel traits of importance to breeding programs. This includes traits related to the nutritional value of milk, the processability of milk into products such as cheese, and traits relevant to animal health and the environment. The ability to successfully incorporate these traits into breeding programs is dependent on the heritability of the FT-MIR predicted traits, and the genetic correlations between the FT-MIR predicted and actual trait values. Linking FT-MIR predicted traits to the underlying mutations responsible for their variation can be difficult because the phenotypic expression of these traits are a function of a diverse range of molecular and biological mechanisms that can obscure their genetic basis. The individual FT-MIR wavenumbers give insights into the chemical composition of milk and provide an additional layer of granularity that may assist with establishing causal links between the genome and observed phenotypes. Additionally, there are other molecular phenotypes such as those related to the metabolome, chromatin accessibility, and RNA editing that could improve our understanding of the underlying biological systems controlling traits of interest. Here we review topics of importance to phenotyping and genetic applications of FT-MIR spectra datasets, and discuss opportunities for consolidating FT-MIR datasets with other genomic and molecular data sources to improve future dairy cattle breeding programs.

Sire selection and genetic improvement of dairy cattle assuming pure market competition.

Breeding companies and farmers rely on selection indices to identify sires they expect to improve production system profitability. Such indices combine estimates of genetic merit for individual traits with corresponding economic values that are fixed based on circumstances at a particular time. Perfect market competition has also been proposed as an economic basis to value alternative sires. The objective of this study was to propose an economic model and develop it for the evaluation of dairy sires. The pure competition model (PUC) was used to evaluate the relationship between a profitability index based on the PUC approach versus the traditional selection index approach for 330 dairy sires comprising Holstein-Friesians, Jerseys, and Ayrshires. The correlation between these 2 selection indices was only 0.56, indicating that the conventional selection index did not correlate well with an index based on the PUC model. In particular, the higher ranking bulls were overvalued using the conventional selection index. Our study concluded that the use of fixed economic values is problematic for the delivery of consistent rankings in selection indices. In contrast, sire rankings based on PUC are more reliable because the sires are evaluated on the basis of efficiency gains rather than production while accounting for market prices and marginal values of dairy outputs over time.

Increasing the accuracy of genomic prediction in pure-bred Limousin beef cattle by including cross-bred Limousin data and accounting for an F94L variant in MSTN.

Explicitly fitting effects for major genes or QTL that account for a large percentage of variation in a whole genomic prediction model may increase prediction accuracy. This study compared approaches to account for a major effect of an F94L variant in the MSTN gene within the genomic prediction using bovine whole-genomic SNP markers. Among the beef cattle breeds, Limousin have been known to have an F94L variant that is not present in Angus. The reference population in this study consisted of 3060 beef cattle including pure-bred Limousin (PL), cross-bred Limousin with Angus (LF) and pure-bred Angus, genotyped using a BovineSNP50 BeadChip and directly for the MSTN-F94L variant. We compared prediction accuracies in PL animals using the three datasets from only the PL population, admixed PL and LF (AL) or multibreed analysis using all of the PL, LF and Angus (MB) population according to four-fold cross-validation after K-means clustering. The MSTN-F94L variant was the most strongly associated with five traits (birth weight, calving ease direct, milk, weaning weight and yield grade) among the 13 measured traits in PL and AL populations. Fitting the MSTN-F94L variant as a random effect, the genomic prediction accuracies for birth weight increased by 2.7% in PL, by 2.2% in AL and by 3.2% in MB. Prediction accuracies for five traits increased in the MB analysis. Fitting MSTN-F94L as a fixed effect in PL, AL and MB analyses resulted in increased prediction accuracy in PL for eight traits. Prediction accuracies can be improved by including a causal variant in genomic evaluation compared with simply using whole-genome SNP markers. Fitting the causal variant as a fixed effect along with markers fitted as random effects resulted in greater prediction accuracies for most traits. Causal variants should be genotyped along with SNP markers.

Strategies for noise reduction and standardization of milk mid-infrared spectra from dairy cattle.

The use of Fourier-transform mid-infrared (FTIR) spectroscopy is of interest to the dairy industry worldwide for predicting milk composition and other novel traits that are difficult or expensive to measure directly. Although there are many valuable applications for FTIR spectra, noise from differences in spectral responses between instruments is problematic because it reduces prediction accuracy if ignored. The purpose of this study was to develop strategies to reduce the impact of noise and to compare methods for standardizing FTIR spectra in order to reduce between-instrument variability in multiple-instrument networks. Noise levels in bands of the infrared spectrum caused by the water content of milk were characterized, and a method for identifying and removing outliers was developed. Two standardization methods were assessed and compared: piecewise direct standardization (PDS), which related spectra on a primary instrument to spectra on 5 other (secondary) instruments using identical milk-based reference samples (n = 918) analyzed across the 6 instruments; and retroactive percentile standardization (RPS), whereby percentiles of observed spectra from routine milk test samples (n = 2,044,094) were used to map and exploit primary- and secondary-instrument relationships. Different applications of each method were studied to determine the optimal way to implement each method across time. Industry-standard predictions of milk components from 2,044,094 spectra records were regressed against predictions from spectra before and after standardization using PDS or RPS. The PDS approach resulted in an overall decrease in root mean square error between industry-standard predictions and predictions from spectra from 0.190 to 0.071 g/100 mL for fat, from 0.129 to 0.055 g/100 mL for protein, and from 0.143 to 0.088 g/100 mL for lactose. Reductions in prediction error for RPS were similar but less consistent than those for PDS across time, but similar reductions were achieved when PDS coefficients were updated monthly and separate primary instruments were assigned for the North and South Islands of New Zealand. We demonstrated that the PDS approach is the most consistent method to reduce prediction errors across time. We also showed that the RPS approach is sensitive to shifts in milk composition but can be used to reduce prediction errors, provided that secondary-instrument spectra are standardized to a primary instrument with samples of broadly equivalent milk composition. Appropriate implementation of either of these approaches will improve the quality of predictions based on FTIR spectra for various downstream applications.

Variance component estimates for post-thaw sperm variables measured by computer assisted semen analyzer for inbred and non-inbred Hereford bulls.

Inbreeding has been shown to have negative effects on animal reproductive performance. It, however, has not been documented as to the extent that inbreeding impacts sperm cell characteristics especially after cells have undergone cryopreservation. In the present study, therefore, there was evaluation of sperm motion and size characteristics from two inbred lines and non-inbred Hereford bulls using Computer Assisted Sperm Analysis (CASA). Pre-frozen semen samples were obtained from Line 1, Prospector, and a non-inbred Hereford bulls. Average inbreeding of each group was 27%, 45%, and 0%, for Line 1, Prospector, and non-inbred groups, respectively. For each bull, two straws (0.5 ml) were evaluated, within each straw two slides with ten fields per slide analyzed with the CASA. Cells were analyzed by mixed model analysis using a model that contained line as fixed effect and individual bull, straw, slide, and fields are random effects. No significant differences between lines of bulls were found for motility or progressive motility. For the head elongation (ELON), and head area (AREA) significant differences existed between lines. Mean separation indicated that Line 1 group had a larger head area and the heads were not as elongated as the Prospector and non-inbred group cells. The lack of significant differences for motility and progressive motility across lines was surprising and indicates there was no inbreeding depression for post-thaw semen characteristics. One explanation for this could be that indirect selection pressure was placed on these characteristics as inbreeding was increased.

Familial episodic ataxia in lambs is potentially associated with a mutation in the fibroblast growth factor 14 (FGF14) gene.

Familial episodic ataxia of lambs is a congenital transient autosomal dominant disorder of newborn lambs, with varying expressivity. Affected lambs show episodes of an asymmetric ataxic gait, base-wide extensor hypertonia of the thoracic limbs and flexor hypertonia of the pelvic limbs. The aim of the study was to determine the genetic variant causing familial episodic ataxia in lambs. Using whole genome sequencing of two half-sib affected lambs, their sire, and their two normal dams, a heterozygous C>T transition at OAR10:77593415 (Oar_v3.1) in exon 1 of the fibroblast growth factor 14 (FGF14) gene (c.46C>T) was identified. The c.46C>T transition resulted in a premature stop codon at position 16 of the 247 amino acid FGF14 protein (p.Q16*). PCR and Sanger sequencing was used to genotype an additional 20 clinically affected animals, demonstrating all lambs carried the c.46C>T variant but 1 clinically more severely affected inbred lamb was homozygous (TT). A further 11 unrelated normal ewes were positionally sequenced, none of which had the variant, while in 18 lambs of unknown status born over 2 years of breeding trials six lambs were found to have the c.46C>T variant, likely clinically unidentified heterozygotes due to the variable expressivity, while 12 did not. In conclusion, familial episodic ataxia of lambs is potentially associated with a c.46C>T variant in the FGF14 gene. Further research is required into the mechanism behind the apparent recovery of lambs.

A comparison of identity-by-descent and identity-by-state matrices that are used for genetic evaluation and estimation of variance components.

The genetic covariance matrix conditional on pedigree is proportional to the pedigree-based additive relationship matrix (PARM), which is twice the matrix of identity-by-descent (IBD) probabilities. In genomic prediction, IBD probabilities in the PARM, which are expected genetic similarities between relatives that are derived from the pedigree, are substituted by realized similarities that are derived from genotypes to obtain a genomic additive relationship matrix (GARM). Different definitions of similarity lead to different GARMs, and two commonly used GARMS are the matrix G, which is based on an allele substitution effect model, and the matrix T, which is based on an allele effect model. We show that although the two matrices T and G are not proportional, they give identical predictions of differences between breeding values. When genomic information is used for variance component estimation, the GARM Gx is computed from genotype covariates that have been standardized to have unit variance. That approach is equivalent to fitting a random regression model using the same standardized covariates. We show that under Hardy-Weinberg and linkage equilibria (LE) that the genetic variance is kσγ2, where σγ2 is the variance of a randomly sampled element from the vector of k substitution effects. However, if linkage disequilibrium (LD) has been generated through selection, covariances between genotypes at different loci will be negative, and therefore, the additive genetic variance will be lower than kσγ2. When the GARM Gx is assumed to be proportional to the genetic covariance matrix, the parameter being estimated is kσγ2. We have demonstrated by simulation that kσγ2 overestimates the additive genetic variance when LD is generated by selection. We argue that unlike the PARM, GARMs are not proportional to a genetic covariance matrix conditional on the observed causal genotypes. The objective here is to recognize the difference between these covariance matrices and its implications.

Genomic study and Medical Subject Headings enrichment analysis of early pregnancy rate and antral follicle numbers in Nelore heifers.

Zebu animals () are known to take longer to reach puberty compared with taurine animals (), limiting the supply of animals for harvest or breeding and impacting profitability. Genomic information can be a helpful tool to better understand complex traits and improve genetic gains. In this study, we performed a genomewide association study (GWAS) to identify genetic variants associated with reproductive traits in Nelore beef cattle. Heifer pregnancy (HP) was recorded for 1,267 genotyped animals distributed in 12 contemporary groups (CG) with an average pregnancy rate of 0.35 (±0.01). Disregarding one of these CG, the number of antral follicles (NF) was also collected for 937 of these animals, with an average of 11.53 (±4.43). The animals were organized in CG: 12 and 11 for HP and NF, respectively. Genes in linkage disequilibrium (LD) with the associated variants can be considered in a functional enrichment analysis to identify biological mechanisms involved in fertility. Medical Subject Headings (MeSH) were detected using the MESHR package, allowing the extraction of broad meanings from the gene lists provided by the GWAS. The estimated heritability for HP was 0.28 ± 0.07 and for NF was 0.49 ± 0.09, with the genomic correlation being -0.21 ± 0.29. The average LD between adjacent markers was 0.23 ± 0.01, and GWAS identified genomic windows that accounted for >1% of total genetic variance on chromosomes 5, 14, and 18 for HP and on chromosomes 2, 8, 11, 14, 15, 16, and 22 for NF. The MeSH enrichment analyses revealed significant ( < 0.05) terms associated with HP-"Munc18 Proteins," "Fucose," and "Hemoglobins"-and with NF-"Cathepsin B," "Receptors, Neuropeptide," and "Palmitic Acid." This is the first study in Nelore cattle introducing the concept of MeSH analysis. The genomic analyses contributed to a better understanding of the genetic control of the reproductive traits HP and NF and provide new selection strategies to improve beef production.

Epistatic interactions associated with fatty acid concentrations of beef from angus sired beef cattle.

BACKGROUND: Consumers are becoming increasingly conscientious about the nutritional value of their food. Consumption of some fatty acids has been associated with human health traits such as blood pressure and cardiovascular disease. Therefore, it is important to investigate genetic variation in content of fatty acids present in meat. Previously publications reported regions of the cattle genome that are additively associated with variation in fatty acid content. This study evaluated epistatic interactions, which could account for additional genetic variation in fatty acid content. RESULTS: Epistatic interactions for 44 fatty acid traits in a population of Angus beef cattle were evaluated with EpiSNPmpi. False discovery rate (FDR) was controlled at 5 % and was limited to well-represented genotypic combinations. Epistatic interactions were detected for 37 triacylglyceride (TAG), 36 phospholipid (PL) fatty acid traits, and three weight traits. A total of 6,181, 7,168, and 0 significant epistatic interactions (FDR < 0.05, 50-animals per genotype combination) were associated with Triacylglyceride fatty acids, Phospholipid fatty acids, and weight traits respectively and most were additive-by-additive interactions. A large number of interactions occurred in potential regions of regulatory control along the chromosomes where genes related to fatty acid metabolism reside. CONCLUSIONS: Many fatty acids were associated with epistatic interactions. Despite a large number of significant interactions, there are a limited number of genomic locations that harbored these interactions. While larger population sizes are needed to accurately validate and quantify these epistatic interactions, the current findings point towards additional genetic variance that can be accounted for within these fatty acid traits.

Impact of fitting dominance and additive effects on accuracy of genomic prediction of breeding values in layers.

Most genomic prediction studies fit only additive effects in models to estimate genomic breeding values (GEBV). However, if dominance genetic effects are an important source of variation for complex traits, accounting for them may improve the accuracy of GEBV. We investigated the effect of fitting dominance and additive effects on the accuracy of GEBV for eight egg production and quality traits in a purebred line of brown layers using pedigree or genomic information (42K single-nucleotide polymorphism (SNP) panel). Phenotypes were corrected for the effect of hatch date. Additive and dominance genetic variances were estimated using genomic-based [genomic best linear unbiased prediction (GBLUP)-REML and BayesC] and pedigree-based (PBLUP-REML) methods. Breeding values were predicted using a model that included both additive and dominance effects and a model that included only additive effects. The reference population consisted of approximately 1800 animals hatched between 2004 and 2009, while approximately 300 young animals hatched in 2010 were used for validation. Accuracy of prediction was computed as the correlation between phenotypes and estimated breeding values of the validation animals divided by the square root of the estimate of heritability in the whole population. The proportion of dominance variance to total phenotypic variance ranged from 0.03 to 0.22 with PBLUP-REML across traits, from 0 to 0.03 with GBLUP-REML and from 0.01 to 0.05 with BayesC. Accuracies of GEBV ranged from 0.28 to 0.60 across traits. Inclusion of dominance effects did not improve the accuracy of GEBV, and differences in their accuracies between genomic-based methods were small (0.01-0.05), with GBLUP-REML yielding higher prediction accuracies than BayesC for egg production, egg colour and yolk weight, while BayesC yielded higher accuracies than GBLUP-REML for the other traits. In conclusion, fitting dominance effects did not impact accuracy of genomic prediction of breeding values in this population.

Genetic parameters and genetic correlations among triacylglycerol and phospholipid fractions in Angus cattle.

The objective of this study was to estimate genetic parameters for intramuscular fatty acids from triacylglycerol (TAG) and phospholipid (PL) fractions in beef LM tissue. Longissimus muscle samples were obtained from 1,833 Angus cattle to determine the intramuscular fatty acid composition for 31 lipids and lipid classes from TAG and PL fractions and were classified by structure into saturated (SFA), monounsaturated (MUFA), polyunsaturated (PUFA), omega-3 (n-3), and omega-6 (n-6) fatty acids. An atherogenic index (AI) was also determined as a measure of the unsaturated fatty acid to SFA ratio. Restricted maximum likelihood methods combined with pedigree data were used to estimate variance components with the WOMBAT software package. Heritability estimates ranged from 0.00 to 0.63 for the major classes of fatty acids. Heritability estimates differed between the TAG and PL fractions, with higher estimates for TAG up to 0.64 and lower estimates for PL that ranged from 0.00 to 0.14. Phenotypic and genetic correlations among individual fatty acids were determined for the TAG fraction as well as among carcass traits, including rib eye area, numerical marbling score, yield grade, ether fat, and Warner-Bratzler shear force value. Strong negative or positive genetic correlations were observed among individual fatty acids in the TAG fraction, which ranged from -0.99 to 0.97 ( < 0.05). Moderate correlations between carcass traits and fatty acids from the TAG fraction ranged from -0.43 to 0.32 ( < 0.05). These results indicate that fatty acids prominent in beef tissues show significant genetic variation as well as genetic relationships with carcass traits.

Genetic parameters for sensory traits in longissimus muscle and their associations with tenderness, marbling score, and intramuscular fat in Angus cattle.

The objective of this study was to estimate heritabilities for sensory traits and genetic correlations among sensory traits and with marbling score (MS), Warner-Bratzler shear force (WBSF), and intramuscular fat content (IMFC). Samples of LM from 2,285 Angus cattle were obtained and fabricated into steaks for laboratory analysis and 1,720 steaks were analyzed by a trained sensory panel. Restricted maximum likelihood procedures were used to obtain estimates of variance and covariance components under a multitrait animal model. Estimates of heritability for MS, IMFC, WBSF, tenderness, juiciness, and connective tissue traits were 0.67, 0.38, 0.19, 0.18, 0.06, and 0.25, respectively. The genetic correlations of MS with tenderness, juiciness, and connective tissue were estimated to be 0.57 ± 0.14, 1.00 ± 0.17, and 0.49 ± 0.13, all positive and strong. Estimated genetic correlations of IMFC with tenderness, juiciness, and connective tissue were 0.56 ± 0.16, 1.00 ± 0.21, and 0.50 ± 0.15, respectively. The genetic correlations of WBSF with tenderness, juiciness, and connective tissue were all favorable and estimated to be -0.99 ± 0.08, -0.33 ± 0.30 and -0.99 ± 0.07, respectively. Strong and positive genetic correlations were estimated between tenderness and juiciness (0.54 ± 0.28) and between connective tissue and juiciness (0.58 ± 0.26). In general, genetic correlations were large and favorable, which indicated that strong relationships exist and similar gene and gene networks may control MS, IMFC, and juiciness or WBSF, panel tenderness, and connective tissue. The results from this study confirm that MS currently used in selection breeding programs has positive genetic correlations with tenderness and juiciness and, therefore, is an effective indicator trait for the improvement of tenderness and juiciness in beef. This study also indicated that a more objective measure, particularly WBSF, a trait not easy to improve through phenotypic selection, is an excellent candidate trait for genomic selection aimed at improving eating satisfaction.