Quantitative genomics of 30 complex phenotypes in Wagyu x Angus F1 progeny.

In the present study, a total of 91 genes involved in various pathways were investigated for their associations with six carcass traits and twenty-four fatty acid composition phenotypes in a Wagyu×Angus reference population, including 43 Wagyu bulls and their potential 791 F(1) progeny. Of the 182 SNPs evaluated, 102 SNPs that were in Hardy-Weinberg equilibrium with minor allele frequencies (MAF>0.15) were selected for parentage assignment and association studies with these quantitative traits. The parentage assignment revealed that 40 of 43 Wagyu sires produced over 96.71% of the calves in the population. Linkage disequilibrium analysis identified 75 of 102 SNPs derived from 54 genes as tagged SNPs. After Bonferroni correction, single-marker analysis revealed a total of 113 significant associations between 44 genes and 29 phenotypes (adjusted P<0.05). Multiple-marker analysis confirmed single-gene associations for 10 traits, but revealed two-gene networks for 9 traits and three-gene networks for 8 traits. Particularly, we observed that TNF (tumor necrosis factor) gene is significantly associated with both beef marbling score (P=0.0016) and palmitic acid (C16:0) (P=0.0043), RCAN1 (regulator of calcineurin 1) with rib-eye area (P=0.0103), ASB3 (ankyrin repeat and SOCS box-containing 3) with backfat (P=0.0392), ABCA1 (ATP-binding cassette A1) with both palmitic acid (C16:0) (P=0.0025) and oleic acid (C18:1n9) (P=0.0114), SLC27A1(solute carrier family 27 A1) with oleic acid (C18:1n9) (P=0.0155), CRH (corticotropin releasing hormone) with both linolenic acid (OMEGA-3) (P=0.0200) and OMEGA 6:3 RATIO (P=0.0054), SLC27A2 (solute carrier family 27 A2) with both linoleic acid (OMEGA-6) (P=0.0121) and FAT (P=0.0333), GNG3 (guanine nucleotide binding protein gamma 3 with desaturase 9 (P=0.0115), and EFEMP1 (EGF containing fibulin-like extracellular matrix protein 1), PLTP (phospholipid transfer protein) and DSEL (dermatan sulfate epimerase-like) with conjugated linoleic acid (P=0.0042-0.0044), respectively, in the Wagyu x Angus F(1) population. In addition, we observed an interesting phenomenon that crossbreeding of different breeds might change gene actions to dominant and overdominant modes, thus explaining the origin of heterosis. The present study confirmed that these important families or pathway-based genes are useful targets for improving meat quality traits and healthful beef products in cattle.

Discovery of novel genetic networks associated with 19 economically important traits in beef cattle.

Quantitative or complex traits are determined by the combined effects of many loci, and are affected by genetic networks or molecular pathways. In the present study, we genotyped a total of 138 mutations, mainly single nucleotide polymorphisms derived from 71 functional genes on a Wagyu x Limousin reference population. Two hundred forty six F(2) animals were measured for 5 carcass, 6 eating quality and 8 fatty acid composition traits. A total of 2,280 single marker-trait association runs with 120 tagged mutations selected based on the HAPLOVIEW analysis revealed 144 significant associations (P < 0.05), but 50 of them were removed from the analysis due to the small number of animals (< or = 9) in one genotype group or absence of one genotype among three genotypes. The remaining 94 single-trait associations were then placed into three groups of quantitative trait modes (QTMs) with additive, dominant and overdominant effects. All significant markers and their QTMs associated with each of these 19 traits were involved in a linear regression model analysis, which confirmed single-gene associations for 4 traits, but revealed two-gene networks for 8 traits and three-gene networks for 5 traits. Such genetic networks involving both genotypes and QTMs resulted in high correlations between predicted and actual values of performance, thus providing evidence that the classical Mendelian principles of inheritance can be applied in understanding genetic complexity of complex phenotypes. Our present study also indicated that carcass, eating quality and fatty acid composition traits rarely share genetic networks. Therefore, marker-assisted selection for improvement of one category of these traits would not interfere with improvement of another.

Genome assembly anchored QTL map of bovine chromosome 14.

Bovine chromosome 14 (BTA14) has been widely explored for quantitative trait loci (QTL) and genes related to economically important traits in both dairy and beef cattle. We reviewed more than 40 investigations and anchored 126 QTL to the current genome assembly (Btau 4_0). Using this anchored QTL map, we observed that, in dairy cattle, the region spanning 0 - 10 Mb on BTA14 has the highest density QTL map with a total of 56 QTL, mainly for milk production traits. It is very likely that both somatic cell score (SCS) and clinical mastitis share some common QTL in two regions: 61.48 Mb - 73.84 Mb and 7.86 Mb - 39.55 Mb, respectively. As well, both ovulation rate and twinning rate might share a common QTL region from 34.16 Mb to 65.38 Mb. However, there are no common QTL locations in three pregnancy related phenotypes: non-return rate, pregnancy rate and daughter pregnancy rate. In beef cattle, the majority of QTL are located in a broad region of 15 Mb - 45 Mb on the chromosome. Functional genes, such as CRH, CYP11B1, DGAT1, FABP4 and TG, as potential candidates for some of these QTL, were also reviewed. Therefore, our review provides a standardized QTL map anchored within the current genome assembly, which would enhance the process of selecting positional and physiological candidate genes for many important traits in cattle.

Evaluation of candidate gene effects for beef backfat via Bayesian model selection.

Candidate gene approaches provide tools for exploring and localizing causative genes affecting quantitative traits and the underlying variation may be better understood by determining the relative magnitudes of effects of their polymorphisms. Diacyglycerol O-acyltransferase 1 (DGAT1), fatty acid binding protein (heart) 3 (FABP3), growth hormone 1 (GH1), leptin (LEP) and thyroglobulin (TG) have been previously identified as genes contributing to genetic control of subcutaneous fat thickness (SFT) in beef cattle. In the present research, Bayesian model selection was used to evaluate effects of these five candidate genes by comparing competing non-nested models and treating candidate gene effects as either random or fixed. The analyses were implemented in SAS to simplify the programming and computation. Phenotypic data were gathered from a F(2) population of Wagyu x Limousin cattle. The five candidate genes had significant but varied effects on SFT in this population. Bayesian model selection identified the DGAT1 model as the one with the greatest model probability, whether candidate gene effects were considered random or fixed, and DGAT1 had the greatest additive effect on SFT. The SAS codes developed in the study are freely available and can be downloaded at: http://www.ansci.wsu.edu/programs/.

Significant associations of the mitochondrial transcription factor A promoter polymorphisms with marbling and subcutaneous fat depth in Wagyu x Limousin F2 crosses.

Mitochondrial transcription factor A (TFAM), a nucleus-encoded protein, regulates the initiation of transcription and replication of mitochondrial DNA (mtDNA). Decreased expression of nuclear-encoded mitochondrial genes has been associated with onset of obesity in mice. Therefore, we hypothesized genetic variants in TFAM gene influence mitochondrial biogenesis consequently affecting body fat deposition and energy metabolism. In the present study, both cDNA (2259 bp) and genomic DNA (16,666 bp) sequences were generated for the bovine TFAM gene using a combination of in silico cloning with targeted region PCR amplification. Alignment of both cDNA and genomic sequences led to the determination of genomic organization and characterization of the promoter region of the bovine TFAM gene. Two closely linked A/C and C/T single nucleotide polymorphisms (SNPs) were found in the bovine TFAM promoter and then genotyped on 237 Wagyu x Limousin F(2) animals with recorded phenotypes for marbling and subcutaneous fat depth (SFD). Statistical analysis demonstrated that both SNPs and their haplotypes were associated with marbling (P=0.0153 for A/C, P=0.0026 for C/T, and P=0.0004 for haplotype) and SFD (P=0.0200 for A/C, P=0.0039 for C/T, and P=0.0029 for haplotype), respectively. A search for transcriptional regulatory elements using MatInspector indicated that both SNPs lead to a gain/loss of six putative-binding sites for transcription factors relevant to fat deposition and energy metabolism. Our results suggest for the first time that TFAM gene plays an important role in lipid metabolism and may be a strong candidate gene for obesity in mammals.

Comparative gene-based in silico analysis of transcriptomes in different bovine tissues and (or) organs.

A gene-based approach was used to annotate 322,168 cattle expressed sequence tags (ESTs) based on human genes in order to census the transcriptomes, analyze their expression similarities, and identify genes preferentially expressed in different bovine tissues and (or) organs. Of the 34,157 human coding genes used in a standalone BLAST search, 14,928 could be matched with provisional orthologous sequences in a total of 230,135 bovine ESTs. The remaining 92,033 bovine ESTs were estimated to represent an additional 5970 genes in cattle. On average, approximately 8600 genes were estimated to be expressed in a single tissue and (or) organ and 13,000 in a pooled tissue library. On the basis of the estimated numbers of genes, no more than 3% of genes would be missed when approximately 34,000 ESTs were sequenced from a single tissue and (or) organ library and approximately 40,000 ESTs from a pooled source, respectively. Cluster analyses of the gene expression patterns among 12 single tissues and (or) organs in cattle revealed that their expression similarities would depend on physiological functions. In addition, a total of 1502 genes were identified as preferentially expressed genes in these 12 single tissues and (or) organs with LOD (logarithm of the odds, base 10) > or = 3.0. Therefore, our study provides some insights for further investigating the developmental and functional relations of various tissues and organs in mammals.