Associations between newly discovered polymorphisms in the Bos taurusgrowth hormone receptor gene and performance traits in Holstein–Friesian dairy cattle

Variations in the growth hormone receptor (GHR) gene sequence are associated with performance traits in cattle. For example, the single nucleotide polymorphism (SNP) F279Y in transmembrane exon 8 has a strong association with milk yield. In this study, 32 previously unreported, putative novel SNPs (31 in the 5′ non‐coding region) were identified by resequencing ～19 kb of the GHR gene in genomic DNA from 22 cattle of multiple breeds. A population of 848 Holstein–Friesian AI sires was subsequently genotyped for the 32 putative novel SNPs and seven published SNPs (including F279Y, one in exon 1A promoter and five in exon 10). Associations between each segregating SNP and genetic merit for performance were quantified in the 848 Holstein–Friesians using weighted animal linear mixed models. Six of the published SNPs and seven of the novel SNPs were associated with at least one of the traits – milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, calving interval, survival and growth and size traits. Even when the allelic substitution effect (P < 0.001) of F279Y was accounted for, the allelic substitution effect of one of the novel SNPs (GHR4.2) in the 5′ non‐coding region of GHR was associated with a lactation milk yield of 37.46 kg (P < 0.001). GHR4.2 and F279Y were not in linkage disequilibrium (r² = 0.00, D’ = 0.04) in the 848 Holstein–Friesians, indicating that their association with milk yield was independent.     

Single-nucleotide polymorphisms in the promoter of the growth hormone-releasing hormone receptor gene are associated with growth and reproduction traits in chickens

Growth hormone‐releasing hormone receptor (GHRHR) plays a critical role in growth hormone (GH) synthesis, release and regulation in animals. The objective of this study was to investigate variations of the chicken GHRHR gene and their associations with growth and reproduction traits in 768 Beijing You chickens. Results revealed three single nucleotide polymorphisms (SNPs) in the promoter region of the gene (g.‐1654A>G, g.‐1411A>G and g.‐142T>C). Association analysis revealed that the novel SNP g.‐1654A>G had significant effects on chicken body weight at 7, 9, 11, 13, 17 weeks of age and the age of first egg as well as egg number at 32, 36 and 40 weeks. Significant association was also observed between g.‐1411A>G and g.‐142T>C with EN24. Moreover, the age of first egg was distinctly related with g.‐142T>C (P < 0.05). Although significant statistical difference was not detected in GHRHR mRNA levels among genotypes of the SNPs (P > 0.05), strong expression variations of the gene were found between the ages 17 and 20 weeks in the population (P < 0.05). These results suggest that the three SNPs in the GHRHR promoter could be used as potential genetic markers to improve the growth and reproductive traits in chickens.     

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.     

A candidate gene association study for nine economically important traits in Italian Holstein cattle

We genotyped 58 single nucleotide polymorphisms (SNPs) in 25 candidate genes in about 800 Italian Holstein sires. Fifty‐six (minor allele frequency >0.02) were used to evaluate their association with single traits: milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP), milk protein percentage (PP), milk somatic cell count (MSCC); and complex indexes: longevity, fertility and productivity–functionality type (PFT), using deregressed proofs, after adjustment for familial relatedness. Thirty‐two SNPs were significantly associated (proportion of false positives <0.05) with different traits: 16 with MSCC, 15 with PY, 14 with MY, 12 with PFT, eight with longevity, eight with FY, eight with PP, five with FP and two with fertility. In particular, a SNP in the promoter region of the PRLR gene was associated with eight of nine traits. DGAT1 polymorphisms were highly associated with FP and FY. Casein gene markers were associated with several traits, confirming the role of the casein gene cluster in affecting milk yield, milk quality and health traits. Other SNPs in genes located on chromosome 6 were associated with PY, PP, PFT, MY (PPARGC1A) and MSCC (KIT). This latter association may suggest a biological link between the degree of piebaldism in Holstein and immunological functions affecting somatic cell count and mastitis resistance. Other significant SNPs were in the ACACA, CRH, CXCR1, FASN, GH1, LEP, LGB (also known as PAEP), MFGE8, SRC, TG, THRSP and TPH1 genes. These results provide information that can complement QTL mapping and genome‐wide association studies in Holstein.     

Genome‐wide association studies for milk production traits in Valle del Belice sheep using repeated measures

Genome‐wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test‐day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.     

CpG site DNA methylation of the CCAAT/enhancer‐binding protein,alpha promoter in chicken lines divergently selected for fatness

CCAAT/enhancer‐binding protein (C/EBP), alpha (CEBPA) is a master regulator of adipogenesis and, together with peroxisome proliferator‐activated receptor gamma (PPARG), plays a critical role in adipocyte differentiation. Previous studies have demonstrated that CEBPA is regulated by DNA methylation and involved in the osteogenesis and adipogenesis of mouse C3H10T1/2 and bone marrow stromal cells. However, it is unclear whether CEBPA is regulated by DNA methylation in adipose tissues. Thus, the objectives of the present study were to investigate CpG site methylation in a 357‐bp CEBPA promoter region and to assess the correlation between promoter CpG site methylation and CEBPA gene expression in the abdominal adipose tissues of Northeast Agricultural University broiler lines divergently selected for abdominal fat content. The results showed that the methylation percentage of the analyzed CEBPA promoter region was significantly higher in lean broilers than in fat broilers at 2 weeks (80.3% vs. 43.4%, P < 0.0001), 3 weeks (95.4% vs. 74.0%, P < 0.0001) and 7 weeks of age (82.6% vs. 57.2%, P < 0.0001). Real‐time quantitative RT‐PCR analysis showed that CEBPA expression was significantly higher in the fat vs. the lean line at 2 weeks of age (P = 0.0013) but not at 3 or 7 weeks of age. The correlation analysis showed that only at 2 weeks of age was the methylation percentage negatively correlated with CEBPA expression (Pearson's r = ?0.8312, P = 0.0029). Of all seven tested CpGs, only two, the CpGs at ?1494 and ?1478 bp, displayed a significantly negative correlation with CEBPA mRNA expression. These results suggest that the CEBPA is methylated in adipose tissue and may regulate chicken early adipose development.     

Single nucleotide polymorphisms in bovine <Emphasis Type="Italic">PRL</Emphasis> gene and their associations with milk production traits in Chinese Holsteins

Prolactin (PRL) plays a crucial role in the initiation and maintenance of lactation in mammals. In this study, seven PCR fragments representing most important functional domains of PRL gene were screened for single nucleotide polymorphisms (SNPs) in Chinese Holsteins by single-strand conformation polymorphisms and amplicons sequencing, and their genetic effects on milk production traits were evaluated. A total of four SNPs, including two in the promoter (−1043A>G and −402A>G), one in intron 1 (+2723C>T) and one in exon 4 (+8398G>A) were identified in PRL gene. Statistical results showed significant associations between the promoter genotypes and the milk performance traits in Chinese Holsteins. Cows with genotype P1-GG showed higher milk yields (P < 0.01), while cows with genotype P1-AA showed higher fat contents (P < 0.01). Haplotype analysis of two SNPs in promoter region revealed that the Hap(AG) was significantly associated with increased milk yields and Hap(AA) was associated with increased fat contents (P < 0.01). This is the second study reporting SNPs in the 5′-regulatory region of PRL gene, which interfere with milk production traits.     

Variation in IGF2BP2 Interacts With Adiposity to Alter Insulin Sensitivity in Mexican Americans

Genome‐wide association studies showed variation in insulin‐like growth factor‐2 binding protein 2 (IGF2BP2) to be associated with type 2 diabetes mellitus (T2DM). We examined a 20‐kb region of IGF2BP2 for association with T2DM‐related quantitative traits in Mexican American families of a proband with gestational diabetes mellitus (GDM) from the BetaGene study. We genotyped 14 single‐nucleotide polymorphisms (SNPs) in 717 individuals from 146 families phenotyped by oral glucose tolerance test (OGTT), intravenous glucose tolerance tests (IVGTTs) with minimal model analysis, and dual‐energy X‐ray absorptiometry scan for percent body fat. Three SNPs and one SNP combination that captured the majority of the variation in the region were tested for association with T2DM‐related quantitative traits using a variance components framework. After correction for multiple testing, rs11705701 showed association with percent body fat (P_ACT = 0.041) with body fat decreasing ～1.5–2% per copy of the A allele. We next tested whether the interaction between rs11705701 and body fat was associated with T2DM‐relative quantitative traits. rs11705701 was significantly associated with insulin sensitivity (Bonferroni P = 0.028) and marginally associated with OGTT 2‐h insulin (Bonferroni P = 0.066) and disposition index (DI) (Bonferroni P = 0.072). We conclude that rs11705701 in IGF2BP2 is associated with body fat and this effect on body fat influences insulin resistance which may contribute to T2DM risk.     

Genome‐wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome‐wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F₂ progeny were analyzed for each blood‐related trait after filtering for quality control. Data were analyzed by the genome‐wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10^?6) on SSC3, 6, 8, 13 and 17 were identified for blood‐related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in pigs.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

Genome‐wide association study reveals novel genes for the ear size in sheep (Ovis aries)

Variations in ear size can be observed in livestock such as sheep; however, the genetic basis of variable ear size in sheep is still poorly understood. To investigate causative genes associated with ear size in sheep, a genome‐wide association study was performed in 115 adult Duolang sheep with different‐sized floppy ears using the Ovine Infinium HD BeadChip. We found 38 significant SNPs at the genome‐wide or chromosome‐wise 5% significance level after Bonferroni correction. The most significant association (P = 1.61 × 10^?6) was found at SNP rs402740419, located in the DCC gene, which plays a critical role in ear development. Also, we observed two additional significant SNPs, rs407891215 in PTPRD and rs407769095 in SOX5, both of which are functionally associated with ear developmental processes. Our results are useful for future sheep breeding and provide insights into the genetic basis of ear size development in sheep and other livestock.     

Targeted genotyping to identify potential functional variants associated with cholesterol content in bovine milk

High blood cholesterol concentration, mainly caused by high dietary cholesterol, is a potential risk factor for human health. Dairy products are important sources of human dietary cholesterol intake. Therefore, monitoring bovine milk cholesterol concentration is important for human health benefit. Genetic selection for improvement of cow milk cholesterol content requires understanding of the genetics of milk cholesterol. For this purpose, we performed analyses of additive and dominance effects of 126 potentially functional SNPs within 43 candidate genes with milk cholesterol content [expressed as mg of cholesterol in 100 g of fat (CHL_fat) or in 100 mg of milk (CHL_milk)]. The additive and dominance effects of SNPs rs380643365 in AGPAT1 (P = 0.04) and rs134357240 in SOAT1 (P = 0.035) genes associated significantly with CHL_fat. Moreover, five (rs109326954 and rs523413537 in DGAT1, rs109376747 in LDLR, rs42781651 in FAM198B and rs109967779 in ACAT2) and four (rs137347384 in RBM19, rs109376747 in LDLR, rs42016945 in PPARG and rs110862179 in SCAP) SNPs were significantly associated with CHL_milk (P < 0.05) based on additive and dominance effect analyses respectively. Rs109326954 and rs523413537 in DGAT1 explained a considerable portion of the phenotypic variance of CHL_milk (7.54 and 6.84% respectively), and might be useful in selection programs for reduced milk cholesterol content. Several significantly associated SNPs were in genes (such as ACAT2 and LDLR) involved in cholesterol metabolism in the liver or cholesterol transport, suggesting multiple mechanisms regulating milk cholesterol content. Nine and seven SNPs identified by additive or dominance effect analyses associated significantly with milk yield and fat yield respectively. Further analyses are required to better understand the consequences of these variants and their potential use in genomic selection of the studied traits.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Genetic markers on BTA14 predictive for residual feed intake in beef steers and their effects on carcass and meat quality traits

With the high cost of feed for animal production, genetic selection for animals that metabolize feed more efficiently could result in substantial cost savings for cattle producers. The purpose of this study was to identify DNA markers predictive for differences among cattle for traits associated with feed efficiency. Crossbred steers were fed a high‐corn diet for 140 days and average daily feed intake (ADFI), average daily gain (ADG), and residual feed intake (RFI) phenotypes were obtained. A region on chromosome 14 was previously associated with RFI in this population of animals. To develop markers with the highest utility for predicting an animal's genetic potential for RFI, we genotyped additional markers within this chromosomal region. These polymorphisms were genotyped on the same animals (n = 1066) and tested for association with ADFI, ADG and RFI. Six markers within this region were associated with RFI (P ≤ 0.05). After conservative correction for multiple testing, one marker at 25.09 Mb remained significant (P = 0.02) and is responsible for 3.6% of the RFI phenotypic variation in this population of animals. Several of these markers were also significant for ADG, although none were significant after correction. Marker alleles with positive effects on ADG corresponded to lower RFI, suggesting an effect increasing growth without increasing feed intake. All markers were also assessed for their effects on meat quality and carcass traits. All of the markers associated with RFI were associated with adjusted fat thickness (AFT, P ≤ 0.009) and three were also associated with hot carcass weight (HCW, P ≤ 0.003). Marker alleles associated with lower RFI were also associated with reduced AFT, and if they were associated for HCW, the effect was an increase in weight. These markers may be useful as prediction tools for animals that utilize feed more efficiently; however, validation with additional populations of cattle is required.     

Genome‐wide association analysis identifies quantitative trait loci for growth in a Landrace purebred population

Growth‐related traits are complex and economically important in the livestock industry. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting growth in pigs. A genome‐wide association study (GWAS) was performed using the porcine single‐nucleotide polymorphism (SNP) 60K bead chip. A mixed‐effects model and linear regression approach were used for the GWAS. The data used in the study included 490 purebred Landrace pigs. All experimental animals were genotyped with 39 438 SNPs located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 16 and body weight at 71 days of age (ALGA0092396, P = 5.35 × 10^?9, Bonferroni adjusted P < 0.05). The SNP marker was located near the genomic region containing IRX4, which encodes iroquois homeobox 4. This SNP marker could be useful in the selective breeding program after validating its effect on other populations.