Matthew L. Spangler

Intraspecific Polymorphism, Interspecific Divergence, and the Origins of Function-Altering Mutations in Deer Mouse Hemoglobin

Major challenges for illuminating the genetic basis of phenotypic evolution are to identify causative mutations, to quantify their functional effects, to trace their origins as new or preexisting variants, and to assess the manner in which segregating variation is transduced into species differences. Here, we report an experimental analysis of genetic variation in hemoglobin (Hb) function within and among species of Peromyscus mice that are native to different elevations. A multilocus survey of sequence variation in the duplicated HBA and HBB genes in Peromyscus maniculatus revealed that function-altering amino acid variants are widely shared among geographically disparate populations from different elevations, and numerous amino acid polymorphisms are also shared with closely related species. Variation in Hb-O2 affinity within and among populations of P. maniculatus is attributable to numerous amino acid mutations that have individually small effects. One especially surprising feature of the Hb polymorphism in P. maniculatus is that an appreciable fraction of functional standing variation in the two transcriptionally active HBA paralogs is attributable to recurrent gene conversion from a tandemly linked HBA pseudogene. Moreover, transpecific polymorphism in the duplicated HBA genes is not solely attributable to incomplete lineage sorting or introgressive hybridization; instead, it is mainly attributable to recurrent interparalog gene conversion that has occurred independently in different species. Partly as a result of concerted evolution between tandemly duplicated globin genes, the same amino acid changes that contribute to variation in Hb function within P. maniculatus also contribute to divergence in Hb function among different species of Peromyscus. In the case of function-altering Hb mutations in Peromyscus, there is no qualitative or quantitative distinction between segregating variants within species and fixed differences between species.

Comparison of molecular breeding values based on within- and across-breed training in beef cattle.

BackgroundAlthough the efficacy of genomic predictors based on within-breed training looks promising, it is necessary to develop and evaluate across-breed predictors for the technology to be fully applied in the beef industry. The efficacies of genomic predictors trained in one breed and utilized to predict genetic merit in differing breeds based on simulation studies have been reported, as have the efficacies of predictors trained using data from multiple breeds to predict the genetic merit of purebreds. However, comparable studies using beef cattle field data have not been reported.MethodsMolecular breeding values for weaning and yearling weight were derived and evaluated using a database containing BovineSNP50 genotypes for 7294 animals from 13 breeds in the training set and 2277 animals from seven breeds (Angus, Red Angus, Hereford, Charolais, Gelbvieh, Limousin, and Simmental) in the evaluation set. Six single-breed and four across-breed genomic predictors were trained using pooled data from purebred animals. Molecular breeding values were evaluated using field data, including genotypes for 2227 animals and phenotypic records of animals born in 2008 or later. Accuracies of molecular breeding values were estimated based on the genetic correlation between the molecular breeding value and trait phenotype.ResultsWith one exception, the estimated genetic correlations of within-breed molecular breeding values with trait phenotype were greater than 0.28 when evaluated in the breed used for training. Most estimated genetic correlations for the across-breed trained molecular breeding values were moderate (> 0.30). When molecular breeding values were evaluated in breeds that were not in the training set, estimated genetic correlations clustered around zero.ConclusionsEven for closely related breeds, within- or across-breed trained molecular breeding values have limited prediction accuracy for breeds that were not in the training set. For breeds in the training set, across- and within-breed trained molecular breeding values had similar accuracies. The benefit of adding data from other breeds to a within-breed training population is the ability to produce molecular breeding values that are more robust across breeds and these can be utilized until enough training data has been accumulated to allow for a within-breed training set.

Estimation of breed-specific heterosis effects for birth, weaning, and yearling weight in cattle.

Heterosis, assumed proportional to expected breed heterozygosity, was calculated for 6834 individuals with birth, weaning and yearling weight records from Cycle VII and advanced generations of the U.S. Meat Animal Research Center (USMARC) Germplasm Evaluation (GPE) project. Breeds represented in these data included: Angus, Hereford, Red Angus, Charolais, Gelbvieh, Simmental, Limousin and Composite MARC III. Heterosis was further estimated by proportions of British × British (B × B), British × Continental (B × C) and Continental × Continental (C × C) crosses and by breed-specific combinations. Model 1 fitted fixed covariates for heterosis within biological types while Model 2 fitted random breed-specific combinations nested within the fixed biological type covariates. Direct heritability estimates (SE) for birth, weaning ,and yearling weight for Model 1 were 0.42 (0.04), 0.22 (0.03), and 0.39 (0.05), respectively. The direct heritability estimates (SE) of birth, weaning, and yearling weight for Model 2 were the same as Model 1, except yearling weight heritability was 0.38 (0.05). The B × B, B × C, and C × C heterosis estimates for birth weight were 0.47 (0.37), 0.75 (0.32), and 0.73 (0.54) kg, respectively. The B × B, B × C, and C × C heterosis estimates for weaning weight were 6.43 (1.80), 8.65 (1.54), and 5.86 (2.57) kg, respectively. Yearling weight estimates for B × B, B × C, and C × C heterosis were 17.59(3.06), 13.88 (2.63), and 9.12 (4.34) kg, respectively. Differences did exist among estimates of breed-specific heterosis for weaning and yearling weight, although the variance component associated with breed-specific heterosis was not significant. These results illustrate that there are differences in breed-specific heterosis and exploiting these differences can lead to varying levels of heterosis among mating plans.

Effects of shade and feeding zilpaterol hydrochloride to finishing steers on performance, carcass quality, heat stress, mobility, and body temperature

Steers ( = 480; 22% with black hides and 78% with red hides) were used to study the effects of shade and feeding zilpaterol hydrochloride (ZH) on performance, carcass quality, heat stress, mobility, and body temperature (BT). A randomized block design with a 2 × 2 factorial treatment arrangement was used with 4 replicates per treatment. Factors included housing type (open or shaded pens) and the feeding of ZH (0 or 8.33 mg/kg DM) the last 21 d on feed with a 3-d withdrawal. Cattle were blocked by BW into a heavy or light block and randomly assigned to pen within each block. Rumen boluses to record BT were inserted before ZH feeding. Respiration rate and panting scores were recorded daily during the ZH feeding period. Mobility scores were collected at various time points from before ZH feeding through harvest. Interactions between ZH and housing type were not significant ( > 0.26) for animal performance, carcass characteristics, and respiration or panting score. No differences ( > 0.44) were observed for DMI, ADG, or G:F on a live basis due to ZH; however, cattle fed in open pens tended ( = 0.08) to have a greater ADG than cattle in shaded pens. Cattle fed ZH had 14 kg heavier carcasses with larger LM area ( < 0.01) than control cattle. Respiration rates for cattle fed ZH were greater ( = 0.05) with no differences ( = 0.88) due to housing. Time affected ( < 0.01) mobility scores, with observations on the morning of harvest at the abattoir being the worst for all groups of cattle. An interaction ( < 0.01) was observed between ZH and housing type for BT. Cattle fed ZH, in both shaded and open pens, had lower ( < 0.05) average, maximum, and area under the curve BT than control cattle fed in the same housing type. However, the observed reduction in BT due to ZH was greater for cattle fed ZH in open pens than for cattle fed ZH in shaded pens. From these results, we conclude that ZH improved HCW with little impact on heat stress or mobility, suggesting that animal welfare was not affected by feeding ZH for 21 d at the end of the feeding period.

Genome-wide association study for feed efficiency and growth traits in U.S. beef cattle

BackgroundSingle nucleotide polymorphism (SNP) arrays for domestic cattle have catalyzed the identification of genetic markers associated with complex traits for inclusion in modern breeding and selection programs. Using actual and imputed Illumina 778K genotypes for 3887 U.S. beef cattle from 3 populations (Angus, Hereford, SimAngus), we performed genome-wide association analyses for feed efficiency and growth traits including average daily gain (ADG), dry matter intake (DMI), mid-test metabolic weight (MMWT), and residual feed intake (RFI), with marker-based heritability estimates produced for all traits and populations.ResultsModerate and/or large-effect QTL were detected for all traits in all populations, as jointly defined by the estimated proportion of variance explained (PVE) by marker effects (PVE ≥ 1.0%) and a nominal P-value threshold (P ≤ 5e-05). Lead SNPs with PVE ≥ 2.0% were considered putative evidence of large-effect QTL (n = 52), whereas those with PVE ≥ 1.0% but < 2.0% were considered putative evidence for moderate-effect QTL (n = 35). Identical or proximal lead SNPs associated with ADG, DMI, MMWT, and RFI collectively supported the potential for either pleiotropic QTL, or independent but proximal causal mutations for multiple traits within and between the analyzed populations. Marker-based heritability estimates for all investigated traits ranged from 0.18 to 0.60 using 778K genotypes, or from 0.17 to 0.57 using 50K genotypes (reduced from Illumina 778K HD to Illumina Bovine SNP50). An investigation to determine if QTL detected by 778K analysis could also be detected using 50K genotypes produced variable results, suggesting that 50K analyses were generally insufficient for QTL detection in these populations, and that relevant breeding or selection programs should be based on higher density analyses (imputed or directly ascertained).ConclusionsFourteen moderate to large-effect QTL regions which ranged from being physically proximal (lead SNPs ≤ 3Mb) to fully overlapping for RFI, DMI, ADG, and MMWT were detected within and between populations, and included evidence for pleiotropy, proximal but independent causal mutations, and multi-breed QTL. Bovine positional candidate genes for these traits were functionally conserved across vertebrate species.

Genome wide association study of cholesterol and poly- and monounsaturated fatty acids, protein, and mineral content of beef from crossbred cattle

The objectives were to determine the variation explained by the BovineSNP50v2 BeadChip for cholesterol (CH), polyunsaturated fatty acids (PUFA), monounsaturated fatty acids (MUFA), protein, and minerals in beef cattle, and to identify chromosomal regions that harbor major allelic variants underlying the variation of these traits. Crossbred steers and heifers (n=236) segregating at the inactive myostatin allele on BTA2 were harvested and steaks were sampled from the M. semitendinosus and the M. longissimus thoracis et lumborum for nutrient analysis. A Bayes C algorithm was employed in genome-wide association analysis. The resulting posterior heritability (SD) estimates ranged from 0.43 (0.10) to 0.71 (0.08) for lipid traits and 0.05 (0.08) to 0.75 (0.06) for mineral traits. Across cuts, correlations between genomic estimated breeding values (GEBV) were similar for CH, MUFA and PUFA. The top 0.5% 1-Mb windows for all traits explained up to 9.93% of the SNP variance. Slight differences did exist between cuts and between different measurement scales of fatty acids.

Breed effects and genetic parameter estimates for calving difficulty and birth weight in a multibreed population

Birth weight (BWT) and calving difficulty (CD) were recorded on 4,579 first-parity females from the Germplasm Evaluation Program at the U.S. Meat Animal Research Center (USMARC). Both traits were analyzed using a bivariate animal model with direct and maternal effects. Calving difficulty was transformed from the USMARC scores to corresponding -scores from the standard normal distribution based on the incidence rate of the USMARC scores. Breed fraction covariates were included to estimate breed differences. Heritability estimates (SE) for BWT direct, CD direct, BWT maternal, and CD maternal were 0.34 (0.10), 0.29 (0.10), 0.15 (0.08), and 0.13 (0.08), respectively. Calving difficulty direct breed effects deviated from Angus ranged from -0.13 to 0.77 and maternal breed effects deviated from Angus ranged from -0.27 to 0.36. Hereford-, Angus-, Gelbvieh-, and Brangus-sired calves would be the least likely to require assistance at birth, whereas Chiangus-, Charolais-, and Limousin-sired calves would be the most likely to require assistance at birth. Maternal breed effects for CD were least for Simmental and Charolais and greatest for Red Angus and Chiangus. Results showed that the diverse biological types of cattle have different effects on both BWT and CD. Furthermore, results provide a mechanism whereby beef cattle producers can compare EBV for CD direct and maternal arising from disjoined and breed-specific genetic evaluations.

The effect of myostatin genotype on body temperature during extreme temperature events.

Extreme heat and cold events can create deleterious physiological changes in cattle as they attempt to cope. The genetic background of animals can influence their response to these events. The objective of the current study was to determine the impact of myostatin genotype (MG) on body temperature during periods of heat and cold stress. Two groups of crossbred steers and heifers of unknown pedigree and breed fraction with varying percentages of Angus, Simmental, and Piedmontese were placed in a feedlot over 2 summers and 2 winters. Before arrival, animals were genotyped for the Piedmontese-derived myostatin mutation (C313Y) to determine their MG as either homozygous normal (0 copy; n = 84), heterozygous (1 copy; n = 96), or homozygous for inactive myostatin (2 copy; n = 59). Hourly tympanic and vaginal temperature measurements were collected for steers and heifers, respectively, for 5 d during times of anticipated heat and cold stress. Mean (±SD) ambient temperature for summer and winter stress events were 24.4 (±4.64) and -1.80 (±11.71), respectively. A trigonometric function (sine + cosine) with periods of 12 and 24 h was used to describe the diurnal cyclical pattern. Hourly body temperature was analyzed within a season, and fixed effects included MG, group, trigonometric functions nested within group, and interaction of MG with trigonometric functions nested within group; random effects were animal and residual (Model [I]). A combined analysis of season and group was also investigated with the inclusion of season as a main effect and the nesting of effects within both group and season (Model [C]). In both models, the residual was fitted using an autoregressive covariance structure. A 3-way interaction of MG, season, and trigonometric function periodicities of 24 h (P < 0.001) and 12 h (P < 0.02) for Model [C] indicate that a genotype × environment interaction exists for MG. For MG during summer stress events the additive estimate was 0.10°C (P < 0.01) and dominance estimate was -0.12°C (P < 0.001). During winter stress events the additive estimate was 0.10°C (P < 0.001) and dominance estimate was 0.054°C (P > 0.05). The current study illustrated that a genotype × environment interaction exists for MG and 1-copy animals were more robust to environmental extremes in comparison with 0- or 2-copy animals.

A simulation study for the analysis of uncertain binary responses: Application to first insemination success in beef cattle

A simulation was carried out to investigate the methods of analyzing uncertain binary responses for success or failure at first insemination. A linear mixed model that included, herd, year, and month of mating as fixed effects; and unrelated service sire, sire and residual as random effects was used to generate binary data. Binary responses were assigned using the difference between days to calving and average gestation length. Females deviating from average gestation length lead to uncertain binary responses. Thus, the methods investigated were the following: (1) a threshold model fitted to certain (no uncertainty) binary data (M1); (2) a threshold model fitted to uncertain binary data ignoring uncertainty (M2); and (3) analysis of uncertain binary data, accounting for uncertainty from day 16 to 26 (M3) or from day 14 to 28 (M4) after introduction of the bull, using a threshold model with fuzzy logic classification. There was virtually no difference between point estimates obtained from M1, M3, and M4 with true values. When uncertain binary data were analyzed ignoring uncertainty (M2), sire variance and heritability were underestimated by 22 and 24%, respectively. Thus, for noisy binary data, a threshold model contemplating uncertainty is needed to avoid bias when estimating genetic parameters.

Genetic parameters for docility, weaning weight, yearling weight, and intramuscular fat percentage in Hereford cattle.

Cattle behavior, including measures of docility, is important to beef cattle producers not only from a human safety perspective but also due to potential correlations to economically relevant traits. Field data from the American Hereford Association was used to estimate genetic parameters for chute score (CS; = 25,037), weaning weight (WW; = 24,908), yearling weight (YW; = 23,978), and intramuscular fat percentage (IMF; = 12,566). Single-trait and bivariate animal models were used to estimate heritabilities and genetic correlations. All models included fixed effects of sex and contemporary group, defined as herd-year-season, and direct genetic and residual components were included as random effects. For CS and WW, additional random effects of maternal genetic and maternal permanent environment were also fitted. For CS, WW, YW, and IMF, heritability estimates were 0.27 ± 0.02, 0.35 ± 0.03, 0.36 ± 0.02, and 0.27 ± 0.02, respectively. Genetic correlations between CS and WW, CS and YW, CS and IMF, WW and YW, WW and IMF, and YW and IMF were -0.12 ± 0.06, -0.10 ± 0.05, -0.08 ± 0.06, 0.47 ± 0.05, -0.19 ± 0.09, and -0.41 ± 0.05, respectively. Heritability estimates for all traits suggest that they would respond favorably to selection and that selection for increased WW or YW could decrease marbling. Genetic correlations between CS and WW, YW, and IMF were all favorable but weak, suggesting that selection for improved docility will not have negative consequences on growth or carcass quality. Furthermore, maternal additive and maternal permanent environmental variances for CS were near 0, suggesting that their inclusion in National Cattle Evaluations is not warranted.