Quantitative trait loci for live animal and carcass composition traits in Jersey and Limousin back-cross cattle finished on pasture or feedlot

Identification of novel candidate genes for age at first calving in Nellore cows using a SNP chip specifically developed for Bos taurus indicus cattle

The age at first calving has a great economic impact on the beef cattle system and calving at 24 months is an objective of selection for a more efficient herd. However, an age at first calving around 36 months has been observed for Nellore cattle in Brazil. Thus, a genome-wide association study (GWAS) was carried out with 8376 records of age at first calving and 3239 animals genotyped with the GGP-Indicus 35K, which has been developed specifically for Bos taurus indicus. The weighted single-step genomic best linear unbiased prediction method was used, with adjacent SNPs (single nucleotide polymorphisms) in genomic windows of 1.0 Mb. After quality control, 3239 (2161 males and 1078 females) animals genotyped for 30,519 SNPs were used in GWAS analysis. The average and standard deviation of age at first calving were 1041.7 and 140.6 days, respectively. The heritability estimate was 0.10 ± 0.02. The GWAS analysis found seven genomic regions in BTA1, 2, 5, 12, 18, 21, and 24, which explained a total of 11.24% of the additive genetic variance of age at first calving. In these regions were found 62 protein coding genes, and the genes HSD17B2, SERPINA14, SERPINA1, SERPINA5, STAT1, NFATC1, ATP9B, CTDP1, THPO, ECE2, PSMD2, EIF4G1, EIF2B2, DVL3, POLR2H, TMTC2, and GPC6 are possible candidates for age at first birth due their function. Moreover, two molecular functions ("serine-type endopeptidase inhibitor activity" and "negative regulation of endopeptidase activity") were significant, which depend on several serpin genes. The use of a SNP chip developed especially for Bos taurus indicus allowed to find genomic regions for age at first calving, which are close to QTLs previously reported for other reproduction-related traits. Future studies can reveal the causal variants and their effects on reproductive precocity of Nellore cows.

Including gene networks to predict calving difficulty in Holstein, Brown Swiss and Jersey cattle

Background

Calving difficulty or dystocia has a great economic impact in the US dairy industry. Reported risk factors associated with calving difficulty are feto-pelvic disproportion, gestation length and conformation. Different dairy cattle breeds have different incidence of calving difficulty, with Holstein having the highest dystocia rates and Jersey the lowest. Genomic selection becomes important especially for complex traits with low heritability, where the accuracy of conventional selection is lower. However, for complex traits where a large number of genes influence the phenotype, genome-wide association studies showed limitations. Biological networks could overcome some of these limitations and better capture the genetic architecture of complex traits. In this paper, we characterize Holstein, Brown Swiss and Jersey breed-specific dystocia networks and employ them in genomic predictions.

Results

Marker association analysis identified single nucleotide polymorphisms explaining the largest average proportion of genetic variance on BTA18 in Holstein, BTA25 in Brown Swiss, and BTA15 in Jersey. Gene networks derived from the genome-wide association included 1272 genes in Holstein, 1454 genes in Brown Swiss, and 1455 genes in Jersey. Furthermore, 256 genes in Holstein network, 275 genes in the Brown Swiss network, and 253 genes in the Jersey network were within previously reported dystocia quantitative trait loci. The across-breed network included 80 genes, with 9 genes being within previously reported dystocia quantitative trait loci. The gene-gene interactions in this network differed in the different breeds. Gene ontology enrichment analysis of genes in the networks showed Regulation of ARF GTPase was very significant (FDR ≤ 0.0098) on Holstein. Neuron morphogenesis and differentiation was the term most enriched (FDR ≤ 0.0539) on the across-breed network. Genomic prediction models enriched with network-derived relationship matrices did not outperform regular GBLUP models.

Conclusions

Regions identified in the genome were in the proximity of previously described quantitative trait loci that would most likely affect calving difficulty by altering the feto-pelvic proportion. Inclusion of identified networks did not increase prediction accuracy. The approach used in this paper could be extended to any instance with asymmetric distribution of phenotypes, for example, resistance to disease data.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0606-y) contains supplementary material, which is available to authorized users.

Polymorphisms and genes associated with puberty in heifers

Puberty onset is a multifactorial process influenced by genetic determinants and environmental conditions, especially nutritional status. Genes, genetic variations, and regulatory networks compose the molecular basis of achieving puberty. In this article, we reviewed the discovery of multiple polymorphisms and genes associated with heifer puberty phenotypes and discuss the opportunities to use this evolving knowledge of genetic determinants for breeding early pubertal Bos indicus-influenced cattle. The discovery of polymorphisms and genes was mainly achieved through candidate gene studies, quantitative trait loci analyses, genome-wide association studies, and recently, global gene expression studies (transcriptome). These studies are recapitulated and summarized in the current review.

Multiple Linkage Disequilibrium Mapping Methods to Validate Additive Quantitative Trait Loci in Korean Native Cattle (Hanwoo)

The efficiency of genome-wide association analysis (GWAS) depends on power of detection for quantitative trait loci (QTL) and precision for QTL mapping. In this study, three different strategies for GWAS were applied to detect QTL for carcass quality traits in the Korean cattle, Hanwoo; a linkage disequilibrium single locus regression method (LDRM), a combined linkage and linkage disequilibrium analysis (LDLA) and a BayesCπ approach. The phenotypes of 486 steers were collected for weaning weight (WWT), yearling weight (YWT), carcass weight (CWT), backfat thickness (BFT), longissimus dorsi muscle area, and marbling score (Marb). Also the genotype data for the steers and their sires were scored with the Illumina bovine 50K single nucleotide polymorphism (SNP) chips. For the two former GWAS methods, threshold values were set at false discovery rate <0.01 on a chromosome-wide level, while a cut-off threshold value was set in the latter model, such that the top five windows, each of which comprised 10 adjacent SNPs, were chosen with significant variation for the phenotype. Four major additive QTL from these three methods had high concordance found in 64.1 to 64.9Mb for Bos taurus autosome (BTA) 7 for WWT, 24.3 to 25.4Mb for BTA14 for CWT, 0.5 to 1.5Mb for BTA6 for BFT and 26.3 to 33.4Mb for BTA29 for BFT. Several candidate genes (i.e. glutamate receptor, ionotropic, ampa 1 [GRIA1], family with sequence similarity 110, member B [FAM110B], and thymocyte selection-associated high mobility group box [TOX]) may be identified close to these QTL. Our result suggests that the use of different linkage disequilibrium mapping approaches can provide more reliable chromosome regions to further pinpoint DNA makers or causative genes in these regions.

Modeling tenderness for genetic and quantitative trait loci analyses

Tenderness is the one of the most important attributes of beef with regards to the intent of consumers to repurchase. Beef tenderness is frequently measured objectively as Warner-Bratzler shear force. However, to better understand the genetics underlying tenderness, better descriptors were derived herein. Two traits, adjusted shear force and aging rate, were calculated to incorporate measures of Warner-Bratzler shear force repeated over aging time of 2 muscles, longissimus dorsi (LD) and semitendinosus (ST). Analysis of these traits indicated that the shear forces of the 2 muscles are lowly correlated and aging rate, as calculated, was not correlated with adjusted shear force. Breed, SNP, and QTL effects were analyzed for the 4 traits (2 traits × 2 muscles). For adjusted shear force, the QTL were located on cattle chromosomes BTA 5, 18, 25, and 29. With the exception of the QTL on BTA 5, the adjusted shear force QTL were observed for either the LD or the ST but not both muscles. The QTL for aging rate were located on BTA 1, 4, 7, 11, 13, 19, and 20.

QTL and association analysis for skin and fibre pigmentation in sheep provides evidence of a major causative mutation and epistatic effects

The pursuits of white features and white fleeces free of pigmented fibre have been important selection objectives for many sheep breeds. The cause and inheritance of non-white colour patterns in sheep has been studied since the early 19th century. Discovery of genetic causes, especially those which predispose pigmentation in white sheep, may lead to more accurate selection tools for improved apparel wool. This article describes an extended QTL study for 13 skin and fibre pigmentation traits in sheep. A total of 19 highly significant, 10 significant and seven suggestive QTL were identified in a QTL mapping experiment using an Awassi × Merino × Merino backcross sheep population. All QTL on chromosome 2 exceeded a LOD score of greater than 4 (range 4.4-30.1), giving very strong support for a major gene for pigmentation on this chromosome. Evidence of epistatic interactions was found for QTL for four traits on chromosomes 2 and 19. The ovine TYRP1 gene on OAR 2 was sequenced as a strong positional candidate gene. A highly significant association (P < 0.01) of grandparental haplotypes across nine segregating SNP/microsatellite markers including one non-synonymous SNP with pigmentation traits could be shown. Up to 47% of the observed variation in pigmentation was accounted for by models using TYRP1 haplotypes and 83% for models with interactions between two QTL probabilities, offering scope for marker-assisted selection for these traits.

A chromosome-wide QTL study on BTA29 affecting temperament traits in German Angus beef cattle and mapping of DRD4

The behaviour of beef cattle is important for the safety and welfare of stockmen and animals. Ten microsatellites spanning BTA29 and, in addition, the candidate gene, dopamine receptor D4 gene, were analysed in 545 German Angus calves of six sires and included in a quantitative trait locus (QTL) study on the basis of three different behaviour tests. A putative QTL for the score while entering the scale (ScE) was detected at BMS764. The DRD4 fragment was mapped in the distal region of BTA29 15.3 cM distal of ILSTS081. The results clearly indicate that BTA29 with a putative QTL in the proximal part and the candidate gene, DRD4, in the distal part plays an important role in the regulation of temperament. During the study one of the sires was detected to be a blood chimera.

Association between GNRHR, LHR and IGF1 polymorphisms and timing of puberty in male Angus cattle

Background

In bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Although the timing of puberty is likely to be a multigenic trait, previous studies indicate that there may also be single genes that exert major effects on the timing of puberty within the general population. Despite its economic importance, there are not many SNPs or genetic markers associated with the age of puberty in male cattle. In the present work, we selected three candidate genes, GNRHR, LHR and IGF1, and associated their polymorphisms with the age of puberty in Angus male cattle.

Results

After weaning, 276 Angus males were measured every month for weight (W), scrotal circumference (SC), sperm concentration (C) and percentage of motility (M). A total of 4 SNPs, two within GNRHR, one in LHR and one in IGF1 were genotyped using the pyrosequencing technique. IGF1-SnaBI SNP was significant associated (P < 0.01) with age at SC 28 cm, but it were not associated with age at M 10% and C 50 million. Genotype CC exhibited an average age at SC 28 cm of 7 and 11 days higher than CT (p = 0.037) and TT (p = 0.012), respectively. This SNP explained 1.5% of the genetic variance of age of puberty at SC28. LHR-I499L, GNRHR-SNP5 and GNRHR-SNP6 were not associated with any of the measurements. However, GNRHR haplotypes showed a suggestive association with age at SC 28 cm.

Conclusions

The findings presented here could support the hypothesis that IGF1 is a regulator of the arrival to puberty in male calves and is involved in the events that precede and initiate puberty in bull calves. Given that most studies in cattle, as well as in other mammals, were done in female, the present results are the first evidence of markers associated with age at puberty in male cattle.

Quantitative trait loci for organ weights and adipose fat composition in Jersey and Limousin back-cross cattle finished on pasture or feedlot

A QTL study of live animal and carcass traits in beef cattle was carried out in New Zealand and Australia. Back-cross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. This paper reports on weights of eight organs (heart, liver, lungs, kidneys, spleen, gastro-intestinal tract, fat, and rumen contents) and 12 fat composition traits (fatty acid (FA) percentages, saturated and monounsaturated FA subtotals, and fat melting point). The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. For organ weights and fat composition traits, 10 and 12 significant QTL locations (P<0.05), respectively, were detected on a genome-wide basis, in combined-sire or within-sire analyses. Seven QTL significant for organ weights were found at the proximal end of chromosome 2. This chromosome carries a variant myostatin allele (F94L), segregating from the Limousin ancestry, and this is a positional candidate for the QTL. Ten significant QTL for fat composition were found on chromosomes 19 and 26. Fatty acid synthase and stearoyl-CoA desaturase (SCD1), respectively, are positional candidate genes for these QTL. Two FA QTL found to be common to sire groups in both populations were for percentages of C14:0 and C14:1 (relative to all FAs) on chromosome 26, near the SCD1 candidate gene.

Genetic mapping of quantitative trait loci for meat quality and muscle metabolic traits in cattle

A whole-genome scan was carried out in New Zealand and Australia to detect quantitative trait loci (QTL) for live animal and carcass composition traits and meat quality attributes in cattle. Backcross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. This paper reports on meat quality traits (tenderness measured as shear force at 4-5 ages on two muscles as well as associated traits of meat colour, pH and cooking loss) and a number of metabolic traits. For meat quality traits, 18 significant QTL (P < 0.05), located in nine linkage groups, were detected on a genome-wise basis, in combined-sire (seven QTL) or within-sire analyses (11 QTL). For metabolic traits, 11 significant QTL (P < 0.05), located in eight linkage groups, were detected on a genome-wise basis, in combined-sire (five QTL) or within-sire analyses (six QTL). BTA2 and BTA3 had QTL for both metabolic traits and meat quality traits. Six significant QTL for meat quality and metabolic traits were found at the proximal end of chromosome 2. BTA2 and BTA29 were the most common chromosomes harbouring QTL for meat quality traits; QTL for improved tenderness were associated with Limousin-derived and Jersey-derived alleles on these two chromosomes, respectively.

Metabolomic profiles indicate distinct physiological pathways affected by two loci with major divergent effect on Bos taurus growth and lipid deposition

Identifying trait-associated genetic variation offers new prospects to reveal novel physiological pathways modulating complex traits. Taking advantage of a unique animal model, we identified the I442M mutation in the non-SMC condensin I complex, subunit G ( NCAPG) gene and the Q204X mutation in the growth differentiation factor 8 ( GDF8) gene as substantial modulators of pre- and/or postnatal growth in cattle. In a combined metabolomic and genotype association approach, which is the first respective study in livestock, we surveyed the specific physiological background of the effects of both loci on body-mass gain and lipid deposition. Our data provided confirming evidence from two historically and geographically distant cattle populations that the onset of puberty is the key interval of divergent growth. The locus-specific metabolic patterns obtained from monitoring 201 plasma metabolites at puberty mirror the particular NCAPG I442M and GDF8 Q204X effects and represent biosignatures of divergent physiological pathways potentially modulating effects on proportional and disproportional growth, respectively. While the NCAPG I442M mutation affected the arginine metabolism, the 204X allele in the GDF8 gene predominantly raised the carnitine level and had concordant effects on glycerophosphatidylcholines and sphingomyelins. Our study provides a conclusive link between the well-described growth-regulating functions of arginine metabolism and the previously unknown specific physiological role of the NCAPG protein in mammalian metabolism. Owing to the confirmed effect of the NCAPG/LCORL locus on human height in genome-wide association studies, the results obtained for bovine NCAPG might add valuable, comparative information on the physiological background of genetically determined divergent mammalian growth.

Mapping Quantitative Trait Loci (QTL) in sheep. III. QTL for carcass composition traits derived from CT scans and aligned with a meta-assembly for sheep and cattle carcass QTL

An (Awassi × Merino) × Merino single-sire backcross family with 165 male offspring was used to map quantitative trait loci (QTL) for body composition traits on a framework map of 189 microsatellite loci across all autosomes. Two cohorts were created from the experimental progeny to represent alternative maturity classes for body composition assessment. Animals were raised under paddock conditions prior to entering the feedlot for a 90-day fattening phase. Body composition traits were derived in vivo at the end of the experiment prior to slaughter at 2 (cohort 1) and 3.5 (cohort 2) years of age, using computed tomography. Image analysis was used to gain accurate predictions for 13 traits describing major fat depots, lean muscle, bone, body proportions and body weight which were used for single- and two-QTL mapping analysis. Using a maximum-likelihood approach, three highly significant (LOD ≥ 3), 15 significant (LOD ≥ 2), and 11 suggestive QTL (1.7 ≤ LOD < 2) were detected on eleven chromosomes. Regression analysis confirmed 28 of these QTL and an additional 17 suggestive (P < 0.1) and two significant (P < 0.05) QTL were identified using this method. QTL with pleiotropic effects for two or more tissues were identified on chromosomes 1, 6, 10, 14, 16 and 23. No tissue-specific QTL were identified.A meta-assembly of ovine QTL for carcass traits from this study and public domain sources was performed and compared with a corresponding bovine meta-assembly. The assembly demonstrated QTL with effects on carcass composition in homologous regions on OAR1, 2, 6 and 21.

Genetic variation in the beta, beta-carotene-9', 10'-dioxygenase gene and association with fat colour in bovine adipose tissue and milk

beta, beta-carotene-9', 10'-dioxygenase (BCO2) plays a role in cleaving beta-carotene eccentrically, and may be involved in the control of adipose and milk colour in cattle. The bovine BCO2 gene was sequenced as a potential candidate gene for a beef fat colour QTL on chromosome (BTA) 15. A single nucleotide base change located in exon 3 causes the substitution of a stop codon (encoded by the A allele) for tryptophan(80) (encoded by the G allele) (c. 240G>A, p.Trp80stop, referred to herein as SNP W80X). Association analysis showed significant differences in subcutaneous fat colour and beta-carotene concentration amongst cattle with different BCO2 genotypes. Animals with the BCO2 AA genotype had more yellow beef fat and a higher beta-carotene concentration in adipose tissues than those with the GA or GG genotype. QTL mapping analysis with the BCO2 SNP W80X fitted as a fixed effect confirmed that this SNP is likely to represent the quantitative trait nucleotide (QTN) for the fat colour-related traits on BTA 15. Moreover, animals with the AA genotype had yellower milk colour and a higher concentration of beta-carotene in the milk.