Acute phase cytokines, TAC1, and toll-like receptor4 mRNA expression and health associated with group size in veal calves

Chronic stressors are a major health and well-being issue in animals. Immune status of animals under chronic stress is compromised, thus reducing disease resistance and compromising well-being of the animal. The objective of this study was to determine the influence of group size of veal calves on immune status and leukocyte mRNA expression of acute phase cytokines, toll-like receptor 4 (TLR4) and tachykinin 1 (TAC1) over a five-month finishing period. Holstein bull calves (n=168), 44±3 days of age were assigned to one of three treatments; 2, 4, or 8 calves/pen (pen space allowance of 1.82m(2)/calf). Jugular blood samples were collected at the day of grouping and then monthly for 4 months. The differential leukocyte counts were determined and mRNA was extracted from the leukocytes. Reverse transcription-qPCR was used to measure the gene expression of interleukin-1 (IL-1β), IL-1 receptor antagonist (IL-1Ra), tumor necrosis factor (TNF-α), TLR4, and TAC1 in leukocytes. Health was evaluated before grouping and monthly for 4 months. On the 1st month after grouping, veal calves that were housed in groups of 8 have greater expression of IL-1β mRNA than calves housed in groups of 4 or 2 (treatment×month, P=0.04). Also at 1 month, groups of 8 had greater TAC1 expression (P<0.05) than calves housed in groups of 4 or 2. However, the expression of IL-1Ra, TNF-α, and TLR4 were not influenced by group size. In the first month of the trial, calves in groups of 8 coughed more (P<0.05) than calves in groups of 2 and coughed more than calves in groups of 4 and 2 during the 2nd month (treatment×month, P=0.03). Calves housed in groups of 8 tended to have greater neutrophil percentage (P=0.09), neutrophil to lymphocyte ratio (P=0.06), and had lower lymphocyte percentage (P=0.06) than those housed in groups of 4 or 2. In conclusion, the number of veal calves in a group, given the same space during the finishing period did not alter IL-1Ra, TNF-α, and TLR4 mRNA expression. However, housing of calves in groups of 8 was associated with greater expression of IL-1β and TAC1 mRNA in peripheral blood leukocytes, and coughing during the first 2 months after grouping. Therefore, housing of veal calves in larger groups may lead to greater susceptibility to respiratory disease and stress.

Effects of rubber flooring during the first 2 lactations on production, locomotion, hoof health, immune functions, and stress

Some housing systems on dairy farms can result in long-term chronic pain. The effects of acute pain on immunity have been explored, but chronic pain's influence on immune responses is still poorly understood. Therefore, the objective of this research was to determine chronic effects of flooring on immune responses and production in freestall housing for dairy cows. Thirty heifers were studied from before calving as first-calf heifers until d 180 of their second lactation. Treatments were rubber (Kraiburg; Agromatic Inc., Fond du Lac, WI) flooring or concrete with diamond grooves in a freestall barn, each in 2 quadrants of the barn. Heifers entered the treatments after calving, so the system was dynamic and each cow was considered an experimental unit. At the end of the first lactation, cows were housed in a bedded pack barn with pasture access until calving was imminent. At that time, they returned to their assigned treatment, but not necessarily into the same quadrant. Production, reproduction, cortisol, acute-phase proteins, and health data were recorded throughout lactation 1, locomotion was scored weekly, and hoof scoring and care was conducted on d 60 and 180 of lactations 1 and 2, and quantitative real-time-PCR of blood leukocytes was analyzed in mid lactation of lactation 1. Mature-equivalent milk fat, milk protein, and protein percentages during the first lactation were greater for cows on the rubber flooring. Hoof and leg therapy treatments per cow were fewer for rubber floor-housed cows. Locomotion scores were less for cows housed on rubber during the second lactation. White blood cell counts were less for cows housed on rubber, and caused by greater lymphocyte counts for cows housed on concrete. The possibility of chronic inflammation was substantiated by less IL-1β and more IL-1 receptor antagonists for cows housed on rubber at d 150 in the second lactation. Cortisol and acute-phase proteins did not differ between the treatments. Interferon-γ, IL-12, the modulator of tissue reconstruction (B-cell-transforming growth factor 1), and pain-modulating neurokinin (tachykinin 1) were not different at d 105. These data show indicators of chronic inflammation for cows housed on the concrete flooring compared with those housed on rubber. Implications for the use of rubber flooring in freestall barns are broader than just lameness and may affect many aspects of cow physiology and production.

The effect of floor surface on dairy cow immune function and locomotion score

This study evaluated the effect of 2 dairy cow housing systems on cow locomotion, immune status, and expression of genes associated with lameness during the dry and periparturient periods. Cows were assigned to freestall housing with either rubber (RUB; n = 13) or concrete (CON; n = 14) at the feed-bunk and alley immediately after their first calving, and managed on this system during all subsequent lactations. At dry off, cows were moved to a straw bedded-pack dry cow pen, and remained there until about 2 d before subsequent calving. To investigate whether greater exposure to RUB or CON increased the differences between cows on each treatment, cows at the end of either their first (n = 16) or second (n = 11) lactations were included in the experiment. Locomotion scores and blood samples were obtained at -60 (beginning of the dry period), -30, 0 (after calving), +10 and +18 d relative to calving. Leukocyte counts were obtained by using an automated cell counter. Phagocytic activity, and cells positive for CD14 and CD18 expression were measured by flow cytometry using labeled microbeads and antibodies. Expression of tachikinin 1(TAC1), histamine receptor 1 (H1), and metalloproteinase (MMP)13 in blood leukocytes was estimated using quantitative real-time PCR. Treatment effects were determined using a repeated measures model. Provision of rubber flooring did not improve dairy cow locomotion during the subsequent study period. However, time relative to calving had an effect on locomotion score and speed, which were worst on d 0, probably because of the discomfort associated with calving. An interaction occurred between treatment and time for neutrophil and lymphocyte counts. The RUB cows had greater neutrophil and lesser lymphocyte numbers postpartum than CON. These cows also had more cells positive for CD14 postpartum compared with prepartum. Moreover, RUB cows showed upregulation of MMP13 and TAC1 compared with CON. These genes are associated with lameness and pain detection respectively. Greater neutrophil to lymphocyte ratios and CD14 expression are associated with physiological stress or with activated immunity. Rubber flooring is associated with an increase in activity and standing. This may have resulted in indications of physiological stress and upregulation of genes associated with lameness and pain for RUB cows. However, this study did not take into account the long-term effects of concrete or rubber flooring; for instance, occurrence of lameness or survivability within the herd.

Potential for estimation of body condition scores in dairy cattle from digital images

Body condition scoring, an indirect measure of the level of subcutaneous fat in dairy cattle, has been widely adopted for research and field assessment or for management purposes on farms. The feasibility of utilizing digital images to determine body condition score (BCS) was assessed for lactating dairy cows at the Scottish Agricultural College Crichton Royal Farm. Two measures of BCS were obtained by using the primary systems utilized in the United Kingdom (UK-BCS) and the United States (USBCS). Means were 2.12 (+/-0.35) and 2.89 (+/-0.40), modes were 2.25 and 2.75, and ranges were 1.0 to 3.5 and 1.5 to 4.5 for the UKBCS (n = 2,346) and USBCS (n = 2,571), respectively. Up to 23 anatomical points were manually identified on images captured automatically as cows passed through a weigh station. Points around the hooks were easier to identify on images than points around pins and the tailhead. All identifiable points were used to define and formulate measures describing the cow's contour. For both BCS systems, hook angle, posterior hook angle, and tailhead depression were significant predictors of BCS. When the full data set testing only the angles around the hooks was used, 100% of predicted BCS were within 0.50 points of actual USBCS and 92.79% were within 0.25 points; and 99.87% of predicted BCS were within 0.50 points of actual UKBCS and 89.95% were within 0.25 points. In a reduced data set considering only observations in which the tailhead depression angle was available, adding the tailhead depression to models did not improve model predictions. The relationships of the calculated angles with USBCS were stronger than those with UKBCS. This research demonstrates the potential for using digital images for assessing BCS. Future efforts should explore ways to automate this process by using a larger number of animals to predict scores accurately for cows across all levels of body condition.

Genotype by Environment Interaction for Production Traits While Accounting for Heteroscedasticity

Grazing (G) provides an alternative management system for dairy production. Heteroscedasticity (HV) of the data may bias estimates of genetic correlations of yield traits between environments, an indicator of genotype-by-environment interaction (GxE). The objective of this study was to investigate the effect of HV on estimates of heritabilities and genetic correlations for mature-equivalent milk, protein, and fat yield, and lactation-average somatic cell scores of daughters, and to determine if HV affects the ability of sire's predicted transmitting ability (PTA) to predict daughter production in G and confinement (C) herds. Data consisted of 72,489 records from 35,674 cows in 366 G herds from 11 states, and 117,629 records from 50,963 cows in 373 C herds from the same 11 states plus 1 geographically contiguous state. Herds were divided into variance quartiles (Q(V)1-Q(V)4) based on milk yield. A transformation was used to reduce HV by standardizing the within-herd standard deviation to the average across-herd standard deviation of a base year for each parity, and was similar to the method used in current USDA-DHIA genetic evaluations. Regression of daughter yield on sire PTA showed that PTA overestimated production of all traits in Q(V)1-Q(V)3 and of milk in Q(V)4 of G herds. For C herds, yields of milk in Q(V)1 and Q(V)2, and of protein and fat in Q(V)1 were overestimated, and protein was underestimated in Q(V)4. Reducing HV had little effect on G herds, but for C herds, regression did not differ from unity for milk and protein in Q(V)1 and Q(V)2. For milk, protein, and fat in G, heritabilities were approximately 0.17, 0.17, and 0.19, respectively. The heritabilities for milk, protein, and fat in C herds were approximately 0.16, 0.17, and 0.21, respectively. Genetic correlations between C and G did not suggest a GxE in 3 upper quartiles, but a possible GxE (correlation = 0.21, estimated standard error = 0.22) for the lowest quartile. Reducing HV did not affect estimates of heritabilities or genetic correlations. Results indicated that modest evidence for existence of GxE did not arise solely from HV.

Prepartum Milking of Heifers Influences Future Production and Health

Transition heifers face multiple stressors during the periparturient period, including first exposure to milking, that may adversely impact dry matter intake (DMI), reduce milk production, compromise immune function, and increase susceptibility to disease. It was hypothesized that reducing the combined stressors experienced at calving would improve the periparturient performance, health, and well-being of heifers. The objective of this study was to determine the effect of initiating the milking procedure 3 wk before expected calving on production, DMI, body weight, energy balance, udder health, calving traits, and health status, as indicated by plasma acute phase protein concentrations. Twenty-two primigravid heifers, blocked by expected calving date, were assigned randomly either to a prepartum milking (PM) group or control group. The PM heifers were milked twice daily beginning at 21 d before expected calving, and control heifers were not milked until after calving. All heifers had access to the same precalving and post-calving diets. Results indicated that PM heifers produced more milk during the first 2 wk after calving and had greater DMI as a percentage of body weight during the first month after calving than did control heifers, although energy balance was more negative for PM heifers. The PM heifers had reduced somatic cell counts through the first month after calving and lower average somatic cell scores during lactation despite having more quarters with mastitis infection at calving. The PM heifers had less udder edema at the third milking postcalving, and had reduced concentrations of haptoglobin in blood sooner than did control heifers. These results indicate that prepartum milking is an alternative management practice that has beneficial effects on the production, health, and well-being of first-lactation cows.

Short Communication: Behavioral and Physiological Indicators of Sensitivity or Chronic Pain Following Tail Docking

Docking the tails of dairy cattle causes mild to moderate behavior changes and physiological indicators of acute pain, but no studies have investigated the possibility that tail docking may lead to chronic pain. In human amputees, an incidence of increased limb surface temperature is associated with phantom limb pain, a central nervous system representation that survives peripheral loss. The objectives of this study were to assess indicators of sensitivity or chronic pain in heifers by using behavioral indicators and thermography. We tested 14 Holstein heifers, 7 docked and 7 intact, from a previous neonatal tail-docking experiment. All 14 animals were videotaped during a test sequence of alternating cold (-9 degrees C), hot (54 degrees C), and neutral packs applied to the underside of the tail. Packs were placed approximately 30.5 cm from the tail head on all animals. A thermal image of the tail was taken using infrared imagery prior to and after temperature sensitivity testing. Docked heifers tended to have greater changes in surface temperatures following the test sequence than did nondocked heifers. In docked heifers, temperatures on the underside of the tail were higher than those at the tip of the tail, both prior to and following the test sequence. Docked heifers also showed substantially higher stomping activity following application of the cold pack. Shifting increased in intact heifers after application of the hot pack, but shifting of the docked heifers did not change. Greater changes were observed in the tail surface temperatures of the docked heifers following temperature manipulation, similar to human amputees who are experiencing phantom limb pain, indicating that similar mechanisms are present in the stump of the docked tail. The behaviors of docked heifers indicated changes in their sensitivity to heat and cold.

Genotype x environment interaction for grazing versus confinement. I. Production traits

The objective of this study was to investigate the possible existence of a genotype x environment interaction (GxE) for production traits of US Holsteins in grazing versus confinement herds. Grazing herds were defined as those that utilized grazing for at least 6 mo and were enrolled in dairy herd improvement (DHI). Control herds were confinement DHI herds of comparable size in similar regions. The performance of daughters in grazing herds and control herds was examined using linear regression of mature equivalent milk, fat, and protein yield on the November 2000 USDA-DHI predicted transmitting abilities (PTA) of their sires for those traits. Heritabilities and genetic correlations were estimated using restricted maximum likelihood in a bivariate animal model that considered the same trait in different environments as different traits. Product-moment and rank correlations were calculated between sires' estimated breeding values, estimated separately in both environments. For grazing herds, the coefficient of regression of milk, fat and protein on PTA were 0.78, 0.76, and 0.78, respectively. Corresponding coefficients in the control herds were 0.99, 0.96, and 0.98. Estimates of heritability for the traits ranged from 0.2 to 0.25, and differences between grazing and control environments were small. Estimates of the genetic correlations for the traits in both environments were 0.89, 0.88, and 0.91 for milk, fat, and protein, respectively. Within-quartile analyses revealed a lower correlation for milk and protein between the upper and lower grazing quartiles, while the same quartiles for the control herds did not differ from unity. Rank correlation coefficients between sire estimated breeding values from the 2 environments were 0.59, 0.63, and 0.66 for milk, fat, and protein, respectively. The mean rank change for the top 100 sires between the two environments was 27. The regression coefficients indicate that expected daughter differences may be overstated by current sire PTA in grazing herds. Genetic correlations less than unity suggests that there is, at least, some reranking among sires in both environments, while the rank correlations indicate the possibility of sire reranking when evaluations were performed within management system. However, differences are not so large as to justify separate genetic evaluations for each system.

Genotype × Environment Interaction for Grazing vs. Confinement. II. Health and Reproduction Traits

Continual selection for increased milk yield for more than 40 yr, combined with the antagonistic relationship between increasing yield, somatic cell count, and fertility, have resulted in sires that may not be optimal for producing daughters for grazing systems where seasonal calving is very important. The objective of this study was to investigate the possible existence of a genotype x environment interaction (G x E) in grazing vs. confinement herds within the United States for lactation average somatic cell score (LSCS), days open (DO), days to first service (DFS), and number of services per conception (SPC). Grazing herds were defined as those that utilized grazing for at least 6 mo each year and were enrolled in Dairy Herd Improvement (DHI). Control herds were confinement DHI herds of similar size in the same states. For LSCS, the performance of daughters in grazing and control herds was examined using linear regression of LSCS on the November 2000 USDA-DHIA sire predicted transmitting abilities (PTA) for SCS. Genetic parameters for all traits were estimated using REML in a bivariate animal model that treated the same trait in different environments as different traits. Rank correlations were calculated between sires' estimated breeding values for LSCS, calculated separately for sires in both environments. The coefficient of regression of daughter LSCS on sire PTA was less in grazing herds than in control herds. The coefficient of regression for control herds was closer to expectation. Estimates of heritability were approximately 0.12 for LSCS, and less than 0.05 for the reproduction traits. Heritabilities for DO, DFS, and SPC were slightly higher for control herds. Estimates of genetic correlation for each reproductive trait between the 2 environments were high and not significantly different from unity. Generally, these traits appear to be under similar genetic control, but a lower coefficient of regression of LSCS on sire PTA for SCS in grazing herds suggests differences in daughter performance in grazing herds may be overstated based on current PTA for SCS.

Genotype x environment interactions in conventional versus pasture-based dairies in Canada

The objective of this study was to evaluate effects of genotype x environmental interaction (GxE) on dairy traits among herds in Canada that practice intensive rotational grazing versus conventional methods that rely on stored feeds. Based on responses to questionnaires, 22 herds were selected for the grazing group and the conventional group consisted of 34 large free-stall dairies. Production data consisted of 6,749 lactations of 2,817 cows in the grazing herds and 29,371 lactations of 12,774 cows in the control herds. Subsets of data for reproduction and conformation also were created. Per cow production was greater in the conventionally managed herds (9,947 kg of milk) than in the grazing herds (9,400 kg). Phenotypic data were regressed on sire predicted transmitting abilities (PTA) and genetic correlations for all traits were estimated between environments. A scaling effect was observed across environments for yield traits, as phenotypic differences among cows were greater than predicted by sire PTA in conventional herds and consistent with sire PTA in grazing herds. Heritabilities for yield also were greater in conventional herds. Other effects of GxE were minor, with genetic correlations of near 0.90 or greater for all traits. Among yield traits, genetic correlation was lowest for fat (0.88 with SE = 0.04). These results demonstrate that effects of GxE are small between the two types of management systems analyzed and that graziers can accurately select sires based on national rankings.

Standardization of lactation means of somatic cell scores for calculation of genetic evaluations

Sample day records of means of somatic cell scores were analyzed to develop adjustments to standardize records for length of lactation. Estimates for effect of DIM were from a model that included random lactation and cow effects for lactations to date of 753,929 Holsteins and 21,842 Jerseys calving in Pennsylvania or Wisconsin from 1987 to 1991. Lactation, cow, and residual variances were estimated using REML. Lactation and cow variances relative to a phenotypic variance of 1.00 were .57 and .31, respectively, for Holsteins and .52 and .35, respectively, for Jerseys. Estimates of effect of DIM were used to compute additive adjustments. Final lactation mean of somatic cell score at < or = 305 DIM for 1,857,532 Holsteins and 113,998 Jerseys from all participating states were standardized for lactation length and analyzed to determine the national effects of calving age and the regional effects of calving month. Multiplicative adjustments were developed for calving age and additive adjustments for calving month. Sample day records of lactation means of somatic cell scores were used to estimate weights based on number of somatic cell sample days to account for the lower accuracy of short records for genetic evaluation.

Age and seasonal effects on Holstein yield for four regions of the United States over time

Multiplicative factors to adjust Holstein yield for age and season of calving were calculated and analyzed from 20 data files for region and time. Regions were 1) California, 2) Minnesota and Wisconsin, 3) New York and Pennsylvania, and 4) six southern states (Texas, Louisiana, Mississippi, Alabama, Georgia, and Florida). Time periods were calving years of 1964 to 1968, 1971 to 1975, 1976 to 1980, 1981 to 1985, and 1986 to 1990. The same statistical model was used that had been used to derive current USDA factors from 1964 to 1968 data. The number of lactation records for each data file ranged from 81,394 to 2,238,201 and increased with time. Effect of calving season on milk yield generally decreased with time; the largest decrease was for the southern region. Effect of calving age also generally decreased for recent data, and the largest decreases were for California. Updated factors to adjust yield for calving age and season are indicated for some regions.

Genetic evaluation of somatic cell scores for United States dairy cattle

Increases in milk yield from genetic selection may be accompanied by correlated increases in genetic susceptibility to clinical mastitis and somatic cells. Unlike clinical mastitis, somatic cell scores can be easily determined and recorded and are related to milk loss from subclinical mastitis. Selection against high somatic cell scores should decrease incidence of clinical mastitis and provide direct economic benefits through higher milk quality premiums. Genetic evaluation for lactation means of linear somatic cell scores has been implemented by USDA and parallels that for yield traits. Because additive genetics accounts for only about 10% of differences in somatic cell scores among cows, more information is needed for the same degree of confidence in genetic estimates as for yield. Only 80% of DHIA cows currently have somatic cell records. Thus, reliabilities of somatic cell evaluations are smaller than those for yield traits. Most progress in selection for lower somatic cell scores will come through sires of cows considered as bull dams. Somatic cell evaluations may best be reported through an economic index with a small amount of emphasis on somatic cell score relative to yield traits. Greater emphasis on somatic cell scores would decrease genetic gain in yield traits, which are economically more important.

Selection on somatic cell score to improve resistance to mastitis in the United States

Increased genetic susceptibility to mastitis has accompanied the rapid genetic increase in milk yield, and genetic selection for mastitis resistance should be considered. Somatic cell score is recommended as an indicator trait to achieve genetic improvement for mastitis resistance. Heritability of somatic cell score is around 10%, and genetic correlation between somatic cell score and clinical mastitis is around .6 to .8. Selection for lower somatic cell score is consistent with the goal of maximizing genetic improvement for total economic merit and should be included in breeding programs. National genetic evaluations for somatic cells scores will use the same animal models and methods as are currently used for milk yield traits. Reliabilities of PTA for somatic cell scores will be smaller than for yield traits because of lower heritability and availability of records from fewer cows. Several forms are proposed for reporting genetic evaluations of somatic cells to producers, and advantages and disadvantages are discussed. Using somatic cell scores for breeding decisions would marginally decrease genetic gain for milk yield and increase total economic merit. Optimal selection indexes would slow the rate of increase in mastitis, rather than decrease its incidence.

Genetic variation in lactation means of somatic cell scores for six breeds of dairy cattle

Lactation means of somatic cell scores from sample days were used to estimate the components of variation for additive genetic, permanent environmental, and herd-sire interaction effects. Data included records of 4314 Ayrshire, 7845 Brown Swiss, 18,115 Guernsey, 1,135,752 Holstein, 67,862 Jersey, and 787 Milking Shorthorn cows from across the US. Records were preadjusted for length of lactation. Fixed effects of herd-year, calving age, and calving month were included in animal models for estimation of variance components. Additive genetic estimates from REML relative to a phenotypic variance of 1.00 were .07 for Ayrshires, .07 for Brown Swiss, .11 for Guernseys, .09 for Holsteins, .09 for Jerseys, and .08 for Milking Shorthorns; permanent environmental estimates were .25, .26, .22, .21, .20, and .35; and herdsire interaction estimates were .04, .02, .00, .02, .02, and .01. Effects of calving age were similar for all regions of the US but differed for Jerseys and Holsteins. Effects of calving month were similar for all breeds. Cows calving during summer had the highest lactation means for somatic cell score from sample days. Impact of calving month was greatest in the Southeast.

The importance of maternal lineage on milk yield traits of dairy cattle

Maternal lineage effects on milk yield traits, considered indicative of cytoplasmic inheritance, were evaluated with animal models. Cattle were from a selection experiment begun in 1968. Maternal pedigrees were traced to the first female member in the Holstein-Friesian Herdbook; purchased cows entering the herd, considered foundation females, were assigned to maternal lineage groups. All models accounted for year-season of calving, parity, and selection lines. Maternal lineage effects were included in a repeated records model with cow effects and preadjustment for sire and maternal grandsire transmitting abilities. Maternal lineage accounted for 5.2, 4.1, and 10.5% of phenotypic variation of preadjusted records of milk yield, fat yield, and fat percentage, respectively. Maternal lineage was evaluated as a fixed effect in an animal model including random animal and permanent environmental effects. Maternal lineage significantly affected fat percentage but not milk yield. Maternal genetic (nuclear) effects and their covariance with additive animal effects did not significantly account for additional variation nor did they influence maternal lineage estimates. Maternal lineage affected calculated net energy of milk but was not important for SNF yield or concentration. Maternal lineage influenced fat percentage, energy concentration, and, to a lesser extent, fat yield in milk of dairy cattle.

Molecular analysis of cytoplasmic genetic variation in Holstein cows

Mitochondrial DNA from Holstein maternal lineages implicated to express cytoplasmic genetic effects on lactation traits was subcloned and screened for molecular polymorphisms. Sixteen of 35 lineages sampled differed from the most common mitochondrial DNA form by at least one restriction endonuclease cleavage site in the 4.3 kilobase segment examined. Variation existed in the region that regulates DNA replication and transcription as well as in transfer and ribosomal ribonucleic acid coding regions of the DNA. The index of nucleotide diversity calculated from polymorphism frequencies indicated that the minimum extent of variation between two random lineages was 1.16 x 10(-4) nucleotide differences per base pair in the segment examined. Presence of a HpaII marker near nucleotide 360 was associated with lower (P less than .001) milk fat percentages. Molecular markers indicated that pedigrees may not be sufficient to separate true cytoplasmic lineages for quantitative genetic analyses. These findings provide a molecular genetic basis for further study of cytoplasmic effects on phenotypic variation in Holstein cattle.