Gene expression patterns during intramuscular fat development in cattle

Deposition of intramuscular fat, or "marbling," in beef cattle contributes significantly to meat quality variables, including juiciness, flavor, and tenderness. The accumulation of intramuscular fat is largely influenced by the genetic background of cattle, as well as their age and nutrition. To identify genes that can be used as early biomarkers for the prediction of marbling capacity, we studied the muscle transcriptome of 2 cattle crossbreeds with contrasting intramuscular fat content. The transcriptomes of marbling LM tissue of heifers from Wagyu x Hereford (WxH; n = 6) and Piedmontese x Hereford (PxH; n = 7) crosses were profiled by using a combination of complementary DNA microarray and quantitative reverse transcription-PCR. Five biopsies of LM were taken from each animal at approximately 3, 7, 12, 20, and 25 mo from birth. Tissue was also collected from the LM of each animal at slaughter (approximately 30 mo). Microarray experiments, conducted on the first 3 biopsies of 2 animals from each crossbreed, identified 97 differentially expressed genes. The gene expression results indicated that the LM transcriptome of animals with high marbling potential (WxH) could be reliably distinguished from less marbled animals (PxH) when the animals were as young as 7 mo of age. At this early age, one cannot reliably determine meaningful differences in intramuscular fat deposition. We observed greater expression of a set of adipogenesis- and lipogenesis-related genes in the LM of young WxH animals compared with their PxH contemporaries. In contrast, genes highly expressed in PxH animals were associated with mitochondrial oxidative activity. Further quantitative reverse transcription-PCR experiments revealed that the messenger RNA of 6 of the lipogenesis-related genes also peaked at the age of 20 to 25 mo in WxH animals. The messenger RNA expression of ADIPOQ, SCD, and THRSP was highly correlated with intramuscular fat content of an individual in WxH animals. Our study provides clear evidence of early molecular changes associated with marbling and also identifies specific time frames when intramuscular fat development in cattle muscle can be detected by using gene expression. This information could be used by animal scientists to design optimal nutrition for high marbling potential. In addition, the genes found to be highly expressed during development of marbling could be used to develop genetic markers or biomarkers to assist with beef production strategies.

A validated whole-genome association study of efficient food conversion in cattle

The genetic factors that contribute to efficient food conversion are largely unknown. Several physiological systems are likely to be important, including basal metabolic rate, the generation of ATP, the regulation of growth and development, and the homeostatic control of body mass. Using whole-genome association, we found that DNA variants in or near proteins contributing to the background use of energy of the cell were 10 times as common as those affecting appetite and body-mass homeostasis. In addition, there was a genic contribution from the extracellular matrix and tissue structure, suggesting a trade-off between efficiency and tissue construction. Nevertheless, the largest group consisted of those involved in gene regulation or control of the phenotype. We found that the distribution of micro-RNA motifs was significantly different for the genetic variants associated with residual feed intake than for the genetic variants in total, although the distribution of promoter sequence motifs was not different. This suggests that certain subsets of micro-RNA are more important for the regulation of this trait. Successful validation depended on the sign of the allelic association in different populations rather than on the strength of the initial association or its size of effect.

Accuracy of prediction of genomic breeding values for residual feed intake and carcass and meat quality traits in Bos taurus, Bos indicus, and composite beef cattle

The aim of this study was to assess the accuracy of genomic predictions for 19 traits including feed efficiency, growth, and carcass and meat quality traits in beef cattle. The 10,181 cattle in our study had real or imputed genotypes for 729,068 SNP although not all cattle were measured for all traits. Animals included Bos taurus, Brahman, composite, and crossbred animals. Genomic EBV (GEBV) were calculated using 2 methods of genomic prediction [BayesR and genomic BLUP (GBLUP)] either using a common training dataset for all breeds or using a training dataset comprising only animals of the same breed. Accuracies of GEBV were assessed using 5-fold cross-validation. The accuracy of genomic prediction varied by trait and by method. Traits with a large number of recorded and genotyped animals and with high heritability gave the greatest accuracy of GEBV. Using GBLUP, the average accuracy was 0.27 across traits and breeds, but the accuracies between breeds and between traits varied widely. When the training population was restricted to animals from the same breed as the validation population, GBLUP accuracies declined by an average of 0.04. The greatest decline in accuracy was found for the 4 composite breeds. The BayesR accuracies were greater by an average of 0.03 than GBLUP accuracies, particularly for traits with known genes of moderate to large effect mutations segregating. The accuracies of 0.43 to 0.48 for IGF-I traits were among the greatest in the study. Although accuracies are low compared with those observed in dairy cattle, genomic selection would still be beneficial for traits that are hard to improve by conventional selection, such as tenderness and residual feed intake. BayesR identified many of the same quantitative trait loci as a genomewide association study but appeared to map them more precisely. All traits appear to be highly polygenic with thousands of SNP independently associated with each trait.

Genome-wide association studies of female reproduction in tropically adapted beef cattle

The genetics of reproduction is poorly understood because the heritabilities of traits currently recorded are low. To elucidate the genetics underlying reproduction in beef cattle, we performed a genome-wide association study using the bovine SNP50 chip in 2 tropically adapted beef cattle breeds, Brahman and Tropical Composite. Here we present the results for 3 female reproduction traits: 1) age at puberty, defined as age in days at first observed corpus luteum (CL) after frequent ovarian ultrasound scans (AGECL); 2) the postpartum anestrous interval, measured as the number of days from calving to first ovulation postpartum (first rebreeding interval, PPAI); and 3) the occurrence of the first postpartum ovulation before weaning in the first rebreeding period (PW), defined from PPAI. In addition, correlated traits such as BW, height, serum IGF1 concentration, condition score, and fatness were also examined. In the Brahman and Tropical Composite cattle, 169 [false positive rate (FPR) = 0.262] and 84 (FPR = 0.581) SNP, respectively, were significant (P < 0.001) for AGECL. In Brahman, 41% of these significant markers mapped to a single chromosomal region on BTA14. In Tropical Composites, 16% of these significant markers were located on BTA5. For PPAI, 66 (FPR = 0.67) and 113 (FPR = 0.432) SNP were significant (P < 0.001) in Brahman and Tropical Composite, respectively, whereas for PW, 68 (FPR = 0.64) and 113 (FPR = 0.432) SNP were significant (P < 0.01). In Tropical Composites, the largest concentration of PPAI markers were located on BTA5 [19% (PPAI) and 23% (PW)], and BTA16 [17% (PPAI) and 18% (PW)]. In Brahman cattle, the largest concentration of markers for postpartum anestrus was located on BTA3 (14% for PPAI and PW) and BTA14 (17% PPAI). Very few of the significant markers for female reproduction traits for the Brahman and Tropical Composite breeds were located in the same chromosomal regions. However, fatness and BW traits as well as serum IGF1 concentration were found to be associated with similar genome regions within and between breeds. Clusters of SNP associated with multiple traits were located on BTA14 in Brahman and BTA5 in Tropical Composites.

Genetic and phenotypic characterisation of animal, carcass, and meat quality traits from temperate and tropically adapted beef breeds. 4. Correlations among animal, carcass, and meat quality traits

Beef cattle data from temperate (TEMP, n = 3947) and tropically (TROP, n = 4137) adapted breeds were analysed to compute estimates of genetic and phenotypic correlations between animal, abattoir carcass, and meat quality measures. Live animal traits included: liveweight (S2LWT), scanned subcutaneous rump fat depth (S2P8), scanned eye muscle area (S2EMA), flight time (S1FT), and finishing average daily gain (FADG). Carcass traits included: hot carcass weight (CWT), retail beef yield percentage (RBY), intramuscular fat percentage (IMF), subcutaneous rump fat depth (P8), eye muscle length by width (ELW), and meat colour score (MEATC). Meat quality measures taken on 2 muscles [M. longissimus thoracis et lumborum (LTL) and M. semitendinosus (ST)] included: shear force of LTL (LTL_SF) and ST (ST_SF); compression of the ST (ST_C); cooking loss % of the LTL (LTL_CL%) and ST (ST_CL%); Minolta LTL L* (LTL_L*), a* (LTL_a*), ST a* (ST_a*); and consumer-assessed LTL tenderness score (LTL_TEND). Genetic and phenotypic correlations between animal measures and related carcass traits were moderate to very high for TEMP and TROP. Genetic correlations between S2LWT and CWT were 0.89 and 0.82, between S2P8 and P8 0.80 and 0.88, and between S2EMA and ELW 0.62 and 0.68, for TEMP and TROP, respectively. Genetic correlations between animal measures and other carcass traits varied; moderate genetic correlations were estimated between S2P8 and RBY (–0.57, –0.19 for TEMP, TROP) and S2P8 and IMF (0.39, 0.23 for TEMP, TROP). Genetic correlations between animal and meat quality measures were moderate to low. For TEMP, moderate genetic correlations were estimated between S2P8 and LTL_TEND (0.38), FADG and ST_a* (–0.49), and FADG and LTL_TEND (0.45); and for TROP, S1FT and LTL_SF (–0.54), and S2EMA and LTL_L* (–0.46). Phenotypic correlations between animal and meat quality were generally low and close to zero. Several moderate to high genetic correlations existed between carcass and meat quality traits. In general, fatness measures were genetically correlated with tenderness (e.g. IMF and LTL_TEND 0.61, 0.31 for TEMP, TROP). CWT was genetically correlated with meat colour (CWT and LTL_L* 0.66, 0.60 for TEMP, TROP) and objective tenderness measures (CWT and ST_C –0.52, –0.22 for TEMP, TROP). Once again phenotypic correlations between carcass and meat quality were low, indicating that few phenotypic predictors of meat quality traits were identified. Several of the genetic correlations show that both animal and abattoir carcass traits may be of use as indirect measures for carcass and meat quality traits in multiple trait genetic evaluation systems.

Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle

A putative functional mutation (rs109231213) near PLAG1 (BTA14) associated with stature was studied in beef cattle. Data from 8199 Bos taurus, Bos indicus and Tropical Composite cattle were used to test the associations between rs109231213 and various phenotypes. Further, 23 496 SNPs located on BTA14 were tested for association with these phenotypes, both independently and fitted together with rs109231213. The C allele of rs109231213 significantly increased hip height, weight, net food intake, age at puberty in males and females and decreased IGF-I concentration in blood and fat depth. When rs109231213 was fitted as a fixed effect in the model, there was an overall reduction in associations between other SNPs and these traits but some SNPs remained associated (P < 10(-4) ). Frequency of the mutant C allele of rs109231213 differed among B. indicus (0.52), B. taurus (0.96) and Tropical Composite (0.68). Most chromosomes carrying the C allele had the same surrounding 10 SNP haplotype, probably because the C allele was introgressed into Brahman from B. taurus cattle. A region of reduced heterozygosity surrounds the C allele; this is small in B. taurus but 20 Mb long in Brahmans, indicating recent and strong selection for the mutant allele. Thus, the C allele appears to mark a mutation that has been selected almost to fixation in the B. taurus breeds studied here and introduced into Brahman cattle during grading up and selected to a frequency of 0.52 despite its negative effects on fertility.

Finding genes for economically important traits: Brahman cattle puberty

Age at puberty is an important component of reproductive performance in beef cattle production systems. Brahman cattle are typically late-pubertal relative to Bos taurus cattle and so it is of economic relevance to select for early age at puberty. To assist selection and elucidate the genes underlying puberty, we performed a genome-wide association study (GWAS) using the BovineSNP50 chip (~54 000 polymorphisms) in Brahman bulls (n = 1105) and heifers (n = 843) and where the heifers were previously analysed in a different study. In a new attempt to generate unbiased estimates of single-nucleotide polymorphism (SNP) effects and proportion of variance explained by each SNP, the available data were halved on the basis of year and month of birth into a calibration and validation set. The traits that defined age at puberty were, in heifers, the age at which the first corpus luteum was detected (AGECL, h2 = 0.56 ± 0.11) and in bulls, the age at a scrotal circumference of 26 cm (AGE26, h2 = 0.78 ± 0.10). At puberty, heifers were on average older (751 ± 142 days) than bulls (555 ± 101 days), but AGECL and AGE26 were genetically correlated (r = 0.20 ± 0.10). There were 134 SNPs associated with AGECL and 146 SNPs associated with AGE26 (P < 0.0001). From these SNPs, 32 (~22%) were associated (P < 0.0001) with both traits. These top 32 SNPs were all located on Chromosome BTA 14, between 21.95 Mb and 28.4 Mb. These results suggest that the genes located in that region of BTA 14 play a role in pubertal development in Brahman cattle. There are many annotated genes underlying this region of BTA 14 and these are the subject of current research. Further, we identified a region on Chromosome X where markers were associated (P < 1.00E–8) with AGE26, but not with AGECL. Information about specific genes and markers add value to our understanding of puberty and potentially contribute to genomic selection. Therefore, identifying these genes contributing to genetic variation in AGECL and AGE26 can assist with the selection for early onset of puberty.

Genome-wide association study for inhibin, luteinizing hormone, insulin-like growth factor 1, testicular size and semen traits in bovine species

The fertility of young bulls impacts on reproduction rates, farm profit and the rate of genetic progress in beef herds. Cattle researchers and industry therefore routinely collect data on the reproductive performance of bulls. Genome-wide association studies were carried out to identify genomic regions and genes associated with reproductive traits measured during the pubertal development of Tropical Composite bulls, from 4 to 24 months of age. Data from 1 085 bulls were collected for seven traits: blood hormone levels of inhibin at 4 months (IN), luteinizing hormone following a gonadotropin releasing hormone challenge at 4 months (LH), insulin-like growth factor 1 at 6 months (IGF1), scrotal circumference at 12 months (SC), sperm motility at 18 months (MOT), percentage of normal spermatozoa at 24 months (PNS) and age at a scrotal circumference of 26 cm (AGE26, or pubertal age). Data from 729 068 single-nucleotide polymorphisms were used in the association analysis. Significant polymorphism associations were discovered for IN, IGF1, SC, AGE26 and PNS. Based on these associations, INHBE, INHBC and HELB are proposed as candidate genes for IN regulation. Polymorphisms associated with IGF1 mapped to the PLAG1 gene region, validating a reported quantitative trait locus on chromosome 14 for IGF1. The X chromosome contained most of the significant associations found for SC, AGE26 and PNS. These findings will contribute to the identification of diagnostic genetic markers and informed genomic selection strategies to assist breeding of cattle with improved fertility. Furthermore, this work provides evidence contributing to gene function annotation in the context of male fertility.

A single nucleotide polymorphism-derived regulatory gene network underlying puberty in 2 tropical breeds of beef cattle

Harsh tropical environments impose serious challenges on poorly adapted species. In beef cattle, tropical adaptation in the form of temperature and disease resistance, coupled with acclimatization to seasonal and limited forage, comes at a cost to production efficiency. Prominent among these costs is delayed onset of puberty, a challenging phenotype to manipulate through traditional breeding mechanisms. Recently, system biology approaches, including gene networks, have been applied to the genetic dissection of complex phenotypes. We aimed at developing and studying gene networks underlying cattle puberty. Our starting material comprises the association results of ~50,000 SNP on 22 traits, including age at puberty, and 2 cattle breed populations: Brahman (n = 843) and Tropical Composite (n = 866). We defined age at puberty as the age at first corpus luteum (AGECL). By capturing the genes harboring mutations minimally associated (P < 0.05) to AGECL or to a set of traits related with AGECL, we derived a gene network for each breed separately and a third network for the combined data set. At the intersection of the 3 networks, we identified candidate genes and pathways that were common to both breeds. Resulting from these analyses, we identified an enrichment of genes involved in axon guidance, cell adhesion, ErbB signaling, and glutamate activity, pathways that are known to affect pulsatile release of GnRH, which is necessary for the onset of puberty. Furthermore, we employed network connectivity and centrality parameters along with a regulatory impact factor metric to identify the key transcription factors (TF) responsible for the molecular regulation of puberty. As a novel finding, we report 5 TF (HIVEP3, TOX, EYA1, NCOA2, and ZFHX4) located in the network intersecting both breeds and interacting with other TF, forming a regulatory network that harmonizes with the recent literature of puberty. Finally, we support our network predictions with evidence derived from gene expression in hypothalamic tissue of adult cows.

Genetic analyses of live-animal ultrasound and abattoir carcass traits in Australian Angus and Hereford cattle

In order to estimate genetic parameters, abattoir carcass data on 1,713 Angus and 1,007 Hereford steers and heifers were combined with yearling live-animal ultrasound measurements on 8,196 Angus and 3,405 Hereford individuals from seedstock herds. Abattoir measures included carcass weight (CWT), percentage of retail beefyield (RBY), near-infrared measured intramuscular fat percentage (CIMF), preslaughter scanned eye muscle area (CEMA), and subcutaneous fat depth at the 12th rib (CRIB) and at the P8 site (CP8). Ultrasound scans on yearling animals included 12th-rib fat depth (SRIB), rump fat depth at the P8 site (SP8), eye muscle area (SEMA), and percentage of intramuscular fat (SIMF). Records on CWT were adjusted to 650-d slaughter age, and the remaining abattoir traits were adjusted to 300-kg CWT. Scan data were analyzed treating records on males and females as different traits. Multivariate analyses were performed on a variety of trait combinations using animal model and REML algorithm. Heritability (h2) estimates for CWT, RBY, CIMF, CP8, CRIB, and CEMA were .31, .68, .43, .44, .28, and .26, respectively, for Angus and .54, .36, .36, .08, .27, .38, respectively, for Hereford. Pooled across sexes, h2 estimates for SIMF, SP8, SRIB, and SEMA were .33, .55, .51, and .42, respectively, for Angus and .20, .31, .18, and .38, respectively, for Hereford. Genetic correlations (r(g)) between the same pair of carcass traits measured at yearling through scanning and directly at the abattoir were moderate to strongly positive, suggesting that selection using yearling ultrasound measurements of seedstock cattle should result in predictable genetic improvement for abattoir carcass characteristics. Estimates of r(g) between the scanned fat measurements and RBY were negative, ranging from -.85 for Angus heifers to -.05 for Hereford heifers. Also, the estimates of r(g) between SEMA and the fat records measured at the abattoir were negative and ranged from -.94 in Hereford heifers to -.02 in Angus heifers.

Gene expression profiling of muscle tissue in Brahman steers during nutritional restriction

Expression profiling using microarrays allows for the detailed characterization of the gene networks that regulate an animal's response to environmental stresses. During nutritional restriction, processes such as protein turnover, connective tissue remodeling, and muscle atrophy take place in the skeletal muscle of the animal. These processes and their regulation are of interest in the context of managing livestock for optimal production efficiency and product quality. Here we expand on recent research applying complementary DNA (cDNA) microarray technology to the study of the effect of nutritional restriction on bovine skeletal muscle. Using a custom cDNA microarray of 9,274 probes from cattle muscle and s.c. fat libraries, we examined the differential gene expression profile of the LM from 10 Brahman steers under three different dietary treatments. The statistical approach was based on mixed-model ANOVA and model-based clustering of the BLUP solutions for the gene x diet interaction effect. From the results, we defined a transcript profile of 156 differentially expressed array elements between the weight loss and weight gain diet substrates. After sequence and annotation analyses, the 57 upregulated elements represented 29 unique genes, and the 99 downregulated elements represented 28 unique genes. Most of these co-regulated genes cluster into groups with distinct biological function related to protein turnover and cytoskeletal metabolism and contribute to our mechanistic understanding of the processes associated with remodeling of muscle tissue in response to nutritional stress.

Prospecting major genes in dairy buffaloes

BACKGROUND

Asian buffaloes (Bubalus bubalis) have an important socio-economic role. The majority of the population is situated in developing countries. Due to the scarce resources in these countries, very few species-specific biotechnology tools exist and a lot of cattle-derived technologies are applied to buffaloes. However, the application of cattle genomic tools to buffaloes is not straightforward and, as results suggested, despite genome sequences similarity the genetic polymorphisms are different.

RESULTS

The first SNP chip genotyping platform designed specifically for buffaloes has recently become available. Herein, a genome-wide association study (GWAS) and gene network analysis carried out in buffaloes is presented. Target phenotypes were six milk production and four reproductive traits. GWAS identified SNP with significant associations and suggested candidate genes that were specific to each trait and also genes with pleiotropic effect, associated to multiple traits.

CONCLUSIONS

Network predictions of interactions between these candidate genes may guide further molecular analyses in search of disruptive mutations, help select genes for functional experiments and evidence metabolism differences in comparison to cattle. The cattle SNP chip does not offer an optimal coverage of buffalo genome, thereafter the development of new buffalo-specific genetic technologies is warranted. An annotated reference genome would greatly facilitate genetic research, with potential impact to buffalo-based dairy production.

A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle

High intramuscular fat (IMF) awards price premiums to beef producers and is associated with meat quality and flavor. Studying gene interactions and pathways that affect IMF might unveil causative physiological mechanisms and inform genomic selection, leading to increased accuracy of predictions of breeding value. To study gene interactions and pathways, a gene network was derived from genetic markers associated with direct measures of IMF, other fat phenotypes, feedlot performance, and a number of meat quality traits relating to body conformation, development, and metabolism that might be plausibly expected to interact with IMF biology. Marker associations were inferred from genomewide association studies (GWAS) based on high density genotypes and 29 traits measured on 10,181 beef cattle animals from 3 breed types. For the network inference, SNP pairs were assessed according to the strength of the correlation between their additive association effects across the 29 traits. The co-association inferred network was formed by 2,434 genes connected by 28,283 edges. Topological network parameters suggested a highly cohesive network, in which the genes are strongly functionally interconnected. Pathway and network analyses pointed towards a trio of transcription factors (TF) as key regulators of carcass IMF: PPARGC1A, HNF4G, and FOXP3. Importantly, none of these genes would have been deemed as significantly associated with IMF from the GWAS. Instead, a total of 313 network genes show significant co-association with the 3 TF. These genes belong to a wide variety of biological functions, canonical pathways, and genetic networks linked to IMF-related phenotypes. In summary, our GWAS and network predictions are supported by the current literature and suggest a cooperative role for the 3 TF and other interacting genes including CAPN6, STC2, MAP2K4, EYA1, COPS5, XKR4, NR2E1, TOX, ATF1, ASPH, TGS1, and TTPA as modulators of carcass and meat quality traits in beef cattle.

The evolution of tropical adaptation: comparing taurine and zebu cattle

Beef cattle breeds consist of three major genetic subdivisions. The taurine group is adapted to temperate environments, and the zebu and Sanga groups are both adapted to tropical environments. With the advent of genotyping and sequencing technologies in agriculture, genome-wide exploration of the genetic basis for the differences in tropical adaptation has only just become possible. In this study, approximately 9000 single nucleotide polymorphism markers were genotyped on 317 animals of a selection of taurine, zebu, and composite breeds to characterize any systematic differences between these groups. We identified 91 intra-breed-class markers; 78 were polymorphic only within the zebu animals, while 13 were polymorphic only in the taurine animals. There were no fixed differences (fixed for alternate alleles between the two breed types) between zebu and taurine animals. We found 14 regions with significantly different allele frequencies between zebu and taurine animals indicative of variable selection pressure or genetic drift. We also found 12 independent regions of differential extended haplotype homozygosity (EHH), indicative of recent selection or rapid fixation of the alternate allele within a short period of time in one of the two breed classes. A preliminary functional genomics analysis of these regions pointed towards signatures of tropical attributes including keratins, heat-shock proteins and heat resistance genes. We anticipate this investigation to be a stepping-stone for future studies to identify genomic regions specific to the two cattle groups, and to subsequently assist in the discrimination between temperate and tropically adapted cattle.

Genetic and phenotypic characterisation of animal, carcass, and meat quality traits from temperate and tropically adapted beef breeds. 3. Meat quality traits

Meat quality measures, including objective measures of tenderness (shear force and compression), were taken on 2 muscles [M. longissimus thoracis et lumborum (LTL) and M. semitendinosus (ST)] from 7566 carcasses from temperate (TEMP) and tropically adapted (TROP) beef cattle breeds. Animals were finished to 1 of 3 market carcass weight end-points (220, 280, or 340 kg) either on pasture or in a feedlot, and in 2 different geographic regions for TROP. Both the phenotypic and genetic expression of the traits were estimated at each market weight and for each finishing regime. Heritabilities and correlations between the traits were estimated for TEMP and TROP separately. Smaller additive variances and heritabilities were observed for temperate breeds compared with tropically adapted breeds for most of the traits studied. For TROP, the heritability of traits measured on the ST muscle [compression (ST_C), shear force (ST_SF), and L* Minolta lightness value (ST_L*)] was 0.27, 0.42, and 0.16, respectively, and for traits measured on the LTL muscle [compression (LTL_C), shear force (LTL_SF), L* Minolta lightness value (LTL_L*), a* Minolta redness value (LTL_a*), cooking loss% (LTL_CL%), and consumer assessed tenderness score (LTL_TEND)] 0.19, 0.30, 0.18, 0.13, 0.20, and 0.31, respectively. For TEMP, the heritability of traits measured on the ST muscle [ST_C, ST_SF, ST_L*, a* Minolta redness value (ST_a*), cooking loss % (ST_CL%)] was 0.12, 0.11, 0.17, 0.13, and 0.15, respectively, and of traits measured on the LTL muscle (LTL_C, LTL_SF, LTL_L, and LTL_TEND) were 0.08, 0.09, 0.17 and 0.18 respectively. Genetic correlations were moderate to high for tenderness measures (shear force and compression) between muscles for the same tenderness measure (e.g. LTL_SF and ST_SF was 0.46 for TROP) and within a muscle for the different measures (e.g. ST_C and ST_SF was 0.83 for TROP). Phenotypic and genetic correlations between LTL_L* and all objective measures of tenderness were negative (e.g. LTL_SF and LTL_L* for TROP was –0.40). The genetic relationship between LTL_SF and LTL_TEND was –0.79 and –0.49 for TROP and TEMP, respectively. Finishing system affected the phenotypic expression of all traits. Pasture-finished, compared with feedlot-finished, animals had higher shear force and compression measures, darker meat colour, and lower sensory tenderness scores for both TEMP and TROP. For TROP, heifers had higher shear force and compression measures, lower sensory tenderness scores, and darker meat colour (lower L* values) than steers. Genetic correlations between markets were generally high and close to unity with the exception of the ST_L*, LTL_L*, ST_C, and ST_SF for TEMP. Geographic region had little effect on the phenotypic and genetic expression of meat quality traits for TROP. Genetic correlations between finishing regimes for all traits were positive and close to unity, with the exception of ST_C and LTL_SF for TEMP, and LTL_L* and LTL_CL% for TROP. Genetic improvement of meat quality traits is a possibility for tropically adapted breeds given the moderate heritabilities, adequate phenotypic variance, generally favourable genetic correlations between traits, and little evidence of genotype by environment interactions.

Gene expression profiling of bovine skeletal muscle in response to and during recovery from chronic and severe undernutrition1

Gene expression profiles of LM from beef cattle that underwent significant postweaning undernutrition were studied using complementary DNA (cDNA) microarrays. After 114 d of undernutrition, the RNA from LM showed 2- to 6-fold less expression of many genes from the classes of muscle structural proteins, muscle metabolic enzymes, and extracellular matrix compared with animals on a rapid growth diet. The expression levels of these genes had mostly returned to pretreatment levels after 84 d of realimentation. The gene expression changes associated with undernutrition and BW loss showed an emphasis on downregulation of gene expression specific to fast-twitch fibers, typical of starving mammals, with a preferential atrophy of glycolytic fast-twitch fibers. We also identified a small group of genes that showed 2- to 5-fold elevated expression in LM after 114 d of undernutrition. Putative roles for these genes in atrophying skeletal muscle are regulation of myogenic differentiation (CSRP3), maintenance of mesenchymal stem cells (CYR61), modulation of membrane function (TM4SF2), prevention of oxidative damage (SESN1), and regulation of muscle protein degradation (SQSTM1). A significant increase in stearoyl-CoA desaturase (SCD) gene expression was observed in atrophying muscle, suggesting either that increased fatty acid synthesis is part of the tissue response to caloric restriction, or that SCD plays another role in energy metabolism in the mixed cellular environment of bovine skeletal muscle.

Gene network analyses of first service conception in Brangus heifers: use of genome and trait associations, hypothalamic-transcriptome information, and transcription factors

Measures of heifer fertility are economically relevant traits for beef production systems and knowledge of candidate genes could be incorporated into future genomic selection strategies. Ten traits related to growth and fertility were measured in 890 Brangus heifers (3/8 Brahman × 5/8 Angus, from 67 sires). These traits were: BW and hip height adjusted to 205 and 365 d of age, postweaning ADG, yearling assessment of carcass traits (i.e., back fat thickness, intramuscular fat, and LM area), as well as heifer pregnancy and first service conception (FSC). These fertility traits were collected from controlled breeding seasons initiated with estrous synchronization and AI targeting heifers to calve by 24 mo of age. The BovineSNP50 BeadChip was used to ascertain 53,692 SNP genotypes for ∼802 heifers. Associations of genotypes and phenotypes were performed and SNP effects were estimated for each trait. Minimally associated SNP (P < 0.05) and their effects across the 10 traits formed the basis for an association weight matrix and its derived gene network related to FSC (57.3% success and heritability = 0.06 ± 0.05). These analyses yielded 1,555 important SNP, which inferred genes linked by 113,873 correlations within a network. Specifically, 1,386 SNP were nodes and the 5,132 strongest correlations (|r| ≥ 0.90) were edges. The network was filtered with genes queried from a transcriptome resource created from deep sequencing of RNA (i.e., RNA-Seq) from the hypothalamus of a prepubertal and a postpubertal Brangus heifer. The remaining hypothalamic-influenced network contained 978 genes connected by 2,560 edges or predicted gene interactions. This hypothalamic gene network was enriched with genes involved in axon guidance, which is a pathway known to influence pulsatile release of LHRH. There were 5 transcription factors with 21 or more connections: ZMAT3, STAT6, RFX4, PLAGL1, and NR6A1 for FSC. The SNP that identified these genes were intragenic and were on chromosomes 1, 5, 9, and 11. Chromosome 5 harbored both STAT6 and RFX4. The large number of interactions and genes observed with network analyses of multiple sources of genomic data (i.e., GWAS and RNA-Seq) support the concept of FSC being a polygenic trait.

A mixture model-based cluster analysis of DNA microarray gene expression data on Brahman and Brahman composite steers fed high-, medium-, and low-quality diets

The objective of this study is to explore aspects of the statistical analysis of gene expression response at the muscle tissue level to varying levels of energy and protein in the diet. Eleven Brahman and Brahman composite steers (weighing 302 +/- 9.8 kg, on average) were allocated randomly into high- (HIGH), medium- (MED), and low- (LOW) quality forage diets for 27 d. After this period, a biopsy of the longissimus dorsi muscle was taken from each animal and total RNA was extracted to generate the labeled target for microarray experimentation. These targets were hybridized to a complementary DNA (cDNA) microarray of 9,274 probes from cattle muscle and subcutaneous fat cDNA libraries. After edits, 151,904 expression intensity levels of 4,747 genes were analyzed. Emphasis was given to the choice of power transformation of the intensity channel readings and to the consistency of readings within each diet quality group. The statistical approach to isolate differentially expressed genes was based on model-based clustering via a mixture of normal distributions estimated through maximal likelihood. The base-2 logarithm was found to be the optimal power transformation to normalize gene intensity levels. A two-sample t-statistic was defined as a measure of possible differential expression. For each of the three diet contrasts, HIGH vs. LOW, HIGH vs. MED, and MED vs. LOW, three clusters were found, two of which contained more than 94% genes with almost no altered gene expression levels, whereas the third cluster contained the remaining genes with a differential expression. Results from the HIGH vs. LOW contrast identified 27 genes with a greater than 95% posterior probability of belonging to the cluster of differentially expressed genes.

Gene expression profiling of Hereford Shorthorn cattle following challenge with Boophilus microplus tick larvae

The ability of cattle to resist tick infestations is partly genetically determined. In order to better define the nature of Bos taurus resistance to the cattle tick Boophilus microplus, skin gene expression was studied using a cattle skin derived cDNA microarray. Expression profiles were determined in skin biopsies sampled from three highly tick resistant animals (HR) and two animals with lower tick resistance (LR) at time 0, immediately before challenge, and again 24 h after challenge. The analysis of the resulting expression data addressed two biological questions: first, for any animal exposed to ticks, which genes are differentially expressed in the 24 h following challenge; and second, which genes are differentially expressed between animals of high and low resistance at 24 h after challenge? In total, 214 genes were found to be differentially expressed in response to larval challenge across all the animals. Seventy-two genes were upregulated and 76 were downregulated at 24 h after challenge. Genes with significantly altered gene expression levels following tick infestation were predominantly keratin genes or mitochondrial genes, as well as odorant binding protein (OBP) and Bos taurus major allergen BDA20. In addition, we identified 66 genes with differential expression between HR and LR animals at 24 h. Of these, genes representing the extracellular matrix and immunoglobulin gene expression pathways were overrepresented. Three differentially expressed genes, OBP, Bos taurus major allergen BDA20 and dendritic cell protein HFL-B5 were further analysed by quantitative reverse transcription PCR (qRT-PCR). The qRT-PCR assay results closely mirrored the expression profiles found in the microarray experiment.

Gene expression phenotypes for lipid metabolism and intramuscular fat in skeletal muscle of cattle

Gene expression phenotypes were evaluated for intramuscular fat (IMF) in bovine skeletal muscle as an alternative to traditional estimates of IMF%. Gene expression data from a time course of LM development in high- and low-marbling Bos taurus cattle crosses were compared to identify genes involved in intramuscular adipocyte lipid metabolism with developmentally similar gene expression profiles. Three sets of genes were identified: triacylglyceride (TAG) synthesis and storage, fatty acid (FA) synthesis, and PPARγ-related genes. In an independent analysis in the LM of 48 Bos indicus cattle, TAG and FA gene sets were enriched in the top 100 genes of which expression was most correlated with IMF% (P = 1.2 × 10(-24) and 3.5 × 10(-9), respectively). In general, genes encoding enzymes involved in the synthesis of FA and TAG in the intramuscular adipocytes were present in the top 100 genes. In B. indicus, effects of a steroid hormone growth promotant (HGP), 2 experimental sites [New South Wales (NSW) and Western Australia (WA)], and 3 tenderness genotypes on the expression levels of genes in the TAG gene set and the correlation of gene expression with IMF% were investigated. Although correlation between expression of 12 individual TAG genes and IMF% was observed in HGP-treated animals in both experimental sites (mean r = 0.43), correlation was not observed for untreated animals at the NSW site (mean r = -0.07, P < 3 × 10(-6)). However, TAG genes showed an average 1.6-fold (P < 0.0004) reduction in expression in the LM of HGP-treated cattle relative to untreated cattle, an effect consistent across both experimental sites. Cattle possessing the favored tenderness calpain 1 and 3 and calpastatin alleles exhibited a greater (P = 0.008) reduction in expression in NSW (1.8-fold reduction, P = 0.0002) compared with WA (1.2-fold reduction, P = 0.03). Tenderness genotype had no impact (P > 0.05) on the correlation of TAG genes with IMF%. In general, the interactions among genotype, treatment and location, and TAG gene set gene expression were consistent with the interactions among the same factors and IMF% detected using 255 animals, of which the 48 in this study were a subset. Thus, the TAG gene set constitutes a gene expression phenotype able to predict effects of different genotypes and treatments on IMF% using much smaller groups than current approaches, even in animals with very low IMF%.

Joint analysis of multiple cDNA microarray studies via multivariate mixed models applied to genetic improvement of beef cattle

In functional genomic laboratories, it is common to use the same microarray slide across studies, each investigating a unique biological question, and each analyzed separately due to computational limitations and/or because there is no hybridization of samples from different studies on one slide. However, the question of analyzing data from multiple studies is a major current issue in microarray data analysis because there are gains to be made in the accuracy of estimated effects by exploiting a covariance structure between gene expression data across studies. We propose an approach for combining multiple studies using multivariate mixed models, with the assumption of a nonzero correlation among genes across experiments, while imposing a null residual covariance. We applied this method to jointly analyze three experiments in genetics of cattle with a total of 54 arrays, each with 19,200 spots and 7,638 elements. The resulting seven-variate model contains 752,476 equations and 56 covariances. To identify differentially expressed genes, we applied model-based clustering to a linear combination of the random gene x variety interaction effect. We enhanced the biological interpretation of the results by applying an iterative algorithm to identify the gene ontology classes that significantly changed in each experiment. We found 118 elements with coordinate expression that clustered into distinct biological functions such as adipogenesis and protein turnover. These results contribute to our understanding of the mechanistic processes involved in adipogenesis and nutrient partitioning.

Physiology and Endocrinology Symposium: How single nucleotide polymorphism chips will advance our knowledge of factors controlling puberty and aid in selecting replacement beef females

The promise of genomic selection is accurate prediction of the genetic potential of animals from their genotypes. Simple DNA tests might replace low-accuracy predictions for expensive or lowly heritable measures of puberty and fertility based on performance and pedigree. Knowing with some certainty which DNA variants (e.g., SNP) affect puberty and fertility is the best way to fulfill the promise. Several SNP from the BovineSNP50 assay have tentatively been associated with reproductive traits including age at puberty, antral follicle count, and pregnancy observed on different sets of heifers. However, sample sizes are too small and SNP density is too sparse to definitively determine genomic regions harboring causal variants affecting reproductive success. Additionally, associations between individual SNP and similar phenotypes are inconsistent across data sets, and genomic predictions do not appear to be globally applicable to cattle of different breeds. Discrepancies may be a result of different QTL segregating in the sampled populations, differences in linkage disequilibrium (LD) patterns such that the same SNP are not correlated with the same QTL, and spurious correlations with phenotype. Several approaches can be used independently or in combination to improve detection of genomic factors affecting heifer puberty and fertility. Larger samples and denser SNP will increase power to detect real associations with SNP having more consistent LD with underlying QTL. Meta-analysis combining results from different studies can also be used to effectively increase sample size. High-density genotyping with heifers pooled by pregnancy status or early and late puberty can be a cost-effective means to sample large numbers. Networks of genes, implicated by associations with multiple traits correlated with puberty and fertility, could provide insight into the complex nature of these traits, especially if corroborated by functional annotation, established gene interaction pathways, and transcript expression. Example analyses are provided to demonstrate how integrating information about gene function and regulation with statistical associations from whole-genome SNP genotyping assays might enhance knowledge of genomic mechanisms affecting puberty and fertility, enabling reliable DNA tests to guide heifer selection decisions.

Genetic and phenotypic characterisation of animal, carcass, and meat quality traits from temperate and tropically adapted beef breeds. 2. Abattoir carcass traits

A total of 11 abattoir carcass measures were recorded on 7854 carcasses in temperate (TEMP) and tropically adapted (TROP) beef breeds. Breeds for TEMP included Angus, Hereford, Murray Grey, and Shorthorn; Brahman, Belmont Red, and Santa Gertrudis accounted for TROP breeds. Measurements included carcass weight (CWT), retail beef yield percentage (RBY), intramuscular fat percentage (IMF), subcutaneous fat depth at the P8 site (P8) and at 12/13th rib (RIB), eye muscle length by width (ELW), deep butt temperature (DBTEMP), fat colour score (FATC), meat colour score (MEATC), marbling score (MARB), and carcass muscle score (MUSC). Animals were finished to 3 different market weight endpoints, either on pasture or in a feedlot, and in different geographic regions for the TROP breeds. Both the phenotypic and genetic expressions of the traits were estimated at each level of market weight endpoint and finishing regime. Heritabilities (h2), and genetic (rg) and phenotypic (rp) correlations between traits were estimated for TEMP and TROP separately. The design effects of market weight endpoint and finishing regimes were the most important sources of variation for continuously measured traits. Main effects for the scored traits were finishing regime for FATC and MEATC and market weight endpoint for MARB and MUSC. Feedlot finished cattle had the whitest FATC and the lightest MEATC. For TEMP, estimates of h2 for CWT, RBY, IMF, P8, RIB, ELW, DBTEMP, FATC, MEATC, MARB, and MUSC were 0.39, 0.57, 0.38, 0.36, 0.27, 0.30, 0.10, 0.05, 0.11, 0.17, and 0.14, respectively. In comparison, h2 for the same order of traits for TROP were 0.36, 0.50, 0.39, 0.30, 0.41, 0.32, 0.04, 0.09, 0.11, 0.25, and 0.11. The direction and magnitude of rg between traits were similar for TEMP and TROP, particularly between CWT, RBY, IMF, P8, and RIB. Genetic correlations of RBY were moderate and negative with all measures of fatness, including IMF (–0.38 TEMP and –0.43 TROP). Positive rg existed between all measures of fatness, with MARB and IMF close to unity. Negative rg was estimated between CWT and all fat measurements. Also negative were the rg and rp estimates between CWT and MEATC. For all traits in both TEMP and TROP, domestic weight carcasses exhibited lower additive variance than export market carcasses. However, genetic correlations between traits across market weight endpoints were positive and close to unity, with the exception of RBY for TROP. For TEMP breeds, genetic correlations between finishing regimes were close to unity. However, possible genotype by environment interactions were found for TROP for P8, MEATC, and MARB between finishing in different geographic regions, and between feedlot and pasture finished animals for RBY and MEATC. Genetic improvement of carcass traits is a possibility given the moderate heritabilities, moderate to strong genetic correlations, and little evidence of genotype by environment interactions.