Pool‐based genome‐wide association study identified novel candidate regions on BTA9 and 14 for oleic acid percentage in Japanese Black cattle

Fatty acid composition is an important indicator of beef quality. The objective of this study was to search the potential candidate region for fatty acid composition. We performed pool‐based genome‐wide association studies (GWAS) for oleic acid percentage (C18:1) in a Japanese Black cattle population from the Hyogo prefecture. GWAS analysis revealed two novel candidate regions on BTA9 and BTA14. The most significant single nucleotide polymorphisms (SNPs) in each region were genotyped in a population (n = 899) to verify their effect on C18:1. Statistical analysis revealed that both SNPs were significantly associated with C18:1 (p = .0080 and .0003), validating the quantitative trait loci (QTLs) detected in GWAS. We subsequently selected VNN1 and LYPLA1 genes as candidate genes from each region on BTA9 and BTA14, respectively. We sequenced full‐length coding sequence (CDS) of these genes in eight individuals and identified a nonsynonymous SNP T66M on VNN1 gene as a putative candidate polymorphism. The polymorphism was also significantly associated with C18:1, but the p value (p = .0162) was higher than the most significant SNP on BTA9, suggesting that it would not be responsible for the QTL. Although further investigation will be needed to determine the responsible gene and polymorphism, our findings would contribute to development of selective markers for fatty acid composition in the Japanese Black cattle of Hyogo.     

Genome‐wide association study for fatty acid composition in Japanese Black cattle

Fatty acid composition is one of the important traits in beef. The aim of this study was to identify candidate genomic regions for fatty acid composition by genome‐wide association study with 50 K single nucleotide polymorphism (SNP) array in Japanese Black cattle. A total of 461 individuals and 40 657 SNPs were used in this study. We applied genome‐wide rapid association using mixed model and regression (GRAMMAR) and genomic control approaches to estimate the associations between genotypes and fatty acid composition. In addition, two SNPs in fatty acid synthase (FASN) (T1952A) and stearoyl‐CoA desaturase (SCD) (V293A) genes were also genotyped. Association analysis revealed that 30 significant SNPs for several fatty acids (C14:0, C14:1, C16:1 and C18:1) were located in the BTA19 FASN gene located within this region but the FASN mutation had no significant effect on any traits. We also detected one significant SNP for C18:1 on BTA23 and two SNPs for C16:0 on BTA25. The region around 17 Mb on BTA26 harbored two significant SNPs for C14:1 and SNP in SCD in this region showed the strongest association with C14:1. This study demonstrated novel candidate regions in BTA19, 23 and 25 for fatty acid composition.     

The g.841G>C SNP of FASN gene is associated with fatty acid composition in beef cattle

The objective of the current study is to evaluate the association between fatty acid composition and fatty acid synthase gene polymorphisms as responsible mutations. For this purpose, we selected seven previously reported single nucleotide polymorphisms (SNPs) in FASN gene, including one within promoter region (g.841G>C) and six non‐synonymous SNPs (g.8805C>T, g.13126C>T, g.15532A>C, g.16024A>G, g.16039C>T, g.17924A>G), and genotyped them in Japanese Black cattle. Genotyping results revealed that g.8805 C>T and g.17924 A>G were monomorphic loci. Genome‐wide association analysis including the other five SNPs revealed that only g.841G>C showed significant associations with the percentages of C14:0, C14:1, C16:1 and C18:1 at 5% genome‐wide significance level. In order to further evaluate the effect, we genotyped g.841G>C using additional three populations, including two Japanese Black populations and a Holstein cattle population. g.16024A>G was also genotyped and included in the analysis because it has been reported to be associated with fatty acid composition in Japanese Black cattle. In the result of analysis of variance, g.841G>C showed stronger effects on fatty acid percentage than those of g.16024A>G in all populations. These results suggested that g.841G>C would be a responsible mutation for fatty acid composition and contribute to production of high‐grade beef as a selection marker in beef cattle.     

UTS2R gene polymorphisms are associated with fatty acid composition in Japanese beef cattle

Fatty acid composition of beef adipose tissue is one of its important traits because a high proportion of monounsaturated fatty acid is related to favorable beef flavor and tenderness. In this study, we searched polymorphisms in full length coding DNA sequence of urotensin 2 recepter and investigated the effects on fatty acid composition (C14:0, C14:1, C16:0, C16:1, C18:0, C18:1, C18:2, monounsaturated fatty acid, saturated fatty acid). Eight single nucleotide polymorphisms (SNP) were identified by sequence comparison among eight animals, including five Japanese Black and three Holstein cattle. One of these SNP (c.866C>T) was predicted to cause amino acid substitutions (P289L) and the other seven synonymous SNP, including c.267C>T, were presumed to be in linkage disequilibrium. Therefore we selected two SNP (c.267C>T and c.866C>T) for further analysis. We investigated associations between these genotypes and fatty acid composition in three Japanese Black populations (n = 560, 245 and 287) and a Holstein population (n = 202). Tukey‐Kramer's honestly significant difference test revealed that CC genotype in c.267C>T indicated lower C14:0 and higher C18:1 than the other genotypes in Japanese Black cattle and CC genotype in c.866C>T showed lower C16:1 than CT genotype in Holstein cattle (P < 0.05). These results suggested that these genotypes would contribute to production of high‐grade meat as selection markers in beef cattle.     

Identification of leptin gene polymorphisms associated with carcass traits and fatty acid composition in Japanese Black cattle

Previous studies have indicated that some leptin gene polymorphisms were associated with economically important traits in cattle breeds. However, polymorphisms in the leptin gene have not been reported thus far in Japanese Black cattle. Here, we aimed to identify the leptin gene polymorphisms which are associated with carcass traits and fatty acid composition in Japanese Black cattle. We sequenced the full‐length coding sequence of leptin gene for eight Japanese Black cattle. Sequence comparison revealed eight single nucleotide polymorphisms (SNPs). Three of these were predicted to cause amino acid substitutions: Y7F, R25C and A80V. Then, we genotyped these SNPs in two populations (JB1 with 560 animals and JB2 with 450 animals) and investigated the effects on the traits. Y7F in JB1 and A80V in JB2 were excluded from statistical analysis because the minor allele frequencies were low (< 0.1). Association analysis revealed that Y7F had a significant effect on the dressed carcass weight in JB2; R25C had a significant effect on C18:0 and C14:1 in JB1 and JB2, respectively; and A80V had a significant effect on C16:0, C16:1, C18:1, monounsaturated fatty acid and saturated fatty acid in JB1. The results suggested that these SNPs could be used as an effective marker for the improvement of Japanese Black cattle.     

Genome‐wide association studies reveal additional related loci for fatty acid composition in a Duroc pig multigenerational population

The aim of the present study was to detect quantitative trait loci affecting fatty acid composition in back fat and intramuscular fat in a Duroc pig population comprising seventh‐generation pedigrees using genome‐wide association studies (GWAS). In total, 305 animals were genotyped using single nucleotide polymorphisms (SNPs) array and five selected SNPs from regions containing known candidate genes related to fatty acid synthesis or metabolism. In total, 24 genome‐wide significant SNP regions were detected in 12 traits, and 76 genome‐wide suggestive SNP regions were detected in 33 traits. The Sus scrofa chromosome (SSC) 7 at 10.3 Mb was significantly associated with C17:0 in intramuscular fat, while the SSC9 at 13.6 Mb was significantly associated with C14:0 in intramuscular fat. The SSC12 at 1.0 Mb was significantly associated with C14:0 in back fat and the SSC14 at 121.0 Mb was significantly associated with C18:0 in intramuscular fat. These regions not only replicated previously reported loci containing some candidate genes involved in fatty acid composition (fatty acid synthase and stearoyl‐CoA desaturase) but also included several additional related loci.     

Detection of candidate polymorphisms around the QTL for fat area ratio to rib eye area on BTA7 using whole‐genome resequencing in Japanese Black cattle

In our previous study, we performed genome‐wide association study (GWAS) to identify the genomic region associated with Fat area ratio to rib eye area (FAR) and detected a candidate in BTA7 at 10–30 Mbp. The present study aims to comprehensively detect all polymorphisms in the candidate region using whole‐genome resequencing data. Based on whole‐genome resequencing of eight animals, we detected 127,090 polymorphisms within the region. Of these, 31,945 were located within the genes. We further narrowed the polymorphisms to 6,044 with more than five allele differences between the high and low FAR groups that were located within 179 genes. We subsequently investigated the functions of these genes and selected 170 polymorphisms in eight genes as possible candidate polymorphisms. We focused on SLC27A6 K81M as a putative candidate polymorphism. We genotyped the SNP in a Japanese Black population (n = 904) to investigate the effect on FAR. Analysis of variance revealed that SLC27A6 K81M had a lower p‐value (p = .0009) than the most significant SNP in GWAS (p = .0049). Although only SLC27A6 K81M was verified in the present study, subsequent verification of the remaining candidate genes and polymorphisms could lead to the identification of genes and polymorphisms responsible for FAR.     

Signatures of selection in Charolais beef cattle identified by genome‐wide analysis

Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8–59.3 Mb on BTA 4; 67.5–69.3 on BTA 5; 39.5–41.0 Mb on BTA 6; and 26.4–29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.     

Multiple marker effects of single nucleotide polymorphisms in three genes,AKIRIN2, EDG1 and RPL27A,for marbling development in Japanese Black cattle

Marbling in beef, measured by Beef Marbling Standard (BMS) number, is an economically important trait for beef cattle breeding and markets in Japan. We previously detected three single nucleotide polymorphisms (SNPs) associated with BMS number of Japanese Black in Oita prefecture: c.*188G>A in AKIRIN2, g.1471620G>T in EDG1 and g.3109537C>T in RPL27A. Here, we carried out single and multiple marker association analyses for the three SNPs in a different commercial Japanese Black population of 892 genotyped animals. The single marker analyses with the model including a single SNP showed significant associations of all SNPs with BMS number. The multiple marker analysis with the model including the main effects of the three SNPs and their interactions detected only significant main effects of g.1471620G>T and g.3109537C>T and a significant interaction between c.*188G>A and g.1471620G>T. These findings suggest the presence of inter‐allelic interactions among genes affecting the development of beef marbling. For effective marker‐assisted selection for BMS number, interactions among these markers need to be considered.     

MC1R c.310G>‐ and c.871G > A determine the coat color of Kumamoto sub‐breed of Japanese Brown cattle

Coat color is one of the important factors characterizing breeds for domestic animals. Melanocortin 1 receptor (MC1R) is a representative responsible gene for this phenotype. Two single‐nucleotide polymorphisms (SNPs) in bovine MC1R gene, c.296T > C and c.310G>‐, have been well characterized, but these SNPs are not enough to explain cattle coat color. As far as we know, MC1R genotypes of Kumamoto sub‐breed of Japanese Brown cattle have not been analyzed. In the current study, genotyping for c.296T > C and c.310G>‐ was performed to elucidate the role of MC1R in determining the coat color of this sub‐breed. As a result, most animals were e/e genotype, suggesting the coat color of this sub‐breed is derived from the e allele of MC1R gene. However, we found six animals with E/e genotype, which coat color would be black theoretically. Subsequently, sequence comparison was performed with these animals to identify other polymorphisms affecting coat color, elucidating that these animals possessed the A allele of c.871G > A commonly. c.871G > A was a non‐synonymous mutation in the seventh transmembrane domain, suggesting alteration of the function and/or the structure of MC1R protein. Our data indicated that the A allele of c.871G > A might be a loss‐of‐function mutation.     

Polymorphism of glutathione S-transferase P1 of dogs with mammary tumours

Mammary tumours constitute more than half of neoplasms in female dogs from different countries. Genome sequences are associated with cancer susceptibility but there is little information available about genetic polymorphisms of glutathione S-transferase P1 (GSTP1) in canine cancers. The aim of this study was to find single nucleotide polymorphisms (SNPs) in GSTP1 of dogs (Canis lupus familiaris) with mammary tumours compared to healthy dogs and to determine the association between GSTP1 polymorphisms and the occurrence of these tumours. The study population included 36 client-owned female dogs with mammary tumours and 12 healthy female dogs, with no previous diagnosis of cancer. DNA was extracted from blood and amplified by PCR assay. PCR-products were sequenced by Sanger method and analysed manually. The 33 polymorphisms were found in GSTP1: 1 coding SNP (exon 4), 24 non-coding SNPs (9 in exon 1), 7 deletions and 1 insertion. The 17 polymorphisms have been found in introns 1, 4, 5 and 6. The dogs with mammary tumours have significant difference from healthy in SNPs I4 c.1018 + 123 T > C (OR 13.412, 95%CI 1.574–114.267, P = .001), I5 c.1487 + 27 T > C (OR 10.737, 95%CI 1.260–91.477, P = .004), I5 c.1487 + 842 G > C (OR 4.714, 95% CI 1.086–20.472, P = .046) and I6 c.2481 + 50 A > G (OR 12.000, 95% CI 1.409–102.207, P = .002). SNP E5 c.1487 T > C and I5 c.1487 + 829 delG also differed significantly (P = .03) but not to the confidence interval. The study, for the first time, showed a positive association of SNPs in GSTP1 with mammary tumours of dogs, that can possibly be used to predict the occurrence of this pathology.     

Three SNPs within exons of INHA and ACVR2B genes are significantly associated with litter size in Dazu black goats

The aim of this study was to analyse the association between single-nucleotide polymorphisms within INHA and ACVR2B and litter size in Dazu black goats. In total, twenty-two SNPs were genotyped in 190 individuals by SNaPshot and resequencing. The results showed that three SNPs (SNP_1, SNP_12 and SNP_13 in this study) were detected to have significant additive genetic effect on the recorded goat litter size (p < .05). The SNP_1 (NC_030809.1), a non-synonymous substitution of G for T at chr2-g. 28314990 in the exon 2 of INHA gene (NM_001285606.1), resulted in homozygote 2 (HOM2) contributed 0.25 and heterozygote (HET) contributed 0.12 larger litter than homozygote 1 (HOM1). Meanwhile, SNP_12 (Chr22-g. 11721225 A > T) and SNP_13 (Chr22-g. 11721227 A > C) (NC_030829.1) simultaneously mutated at the first and third position of a triplet AAA (lysine, K) in the exon 4 of ACVR2B gene (XM_018066623.1) had estimated genetic effects of HOM1 (0.00) and HOM2 (0.03) larger than HET (−0.12). In conclusion, one SNPs (chr2-g. 28314990 T > G) within the exon 2 of INHA and two SNPs (Chr22-g. 11721225 A > T and Chr22-g. 11721227 A > C) i n the exon 4 of ACVR2B gene were highly recommended as candidate markers of litter size in Dazu black goats. A large-scale association study to assess the impact of these variants on litter size is still necessary.     

A genomewide association mapping study using ultrasound‐scanned information identifies potential genomic regions and candidate genes affecting carcass traits in Nellore cattle

The aim of this study was to identify candidate genes and genomic regions associated with ultrasound‐derived measurements of the rib‐eye area (REA), backfat thickness (BFT) and rumpfat thickness (RFT) in Nellore cattle. Data from 640 Nellore steers and young bulls with genotypes for 290 863 single nucleotide polymorphisms (SNPs) were used for genomewide association mapping. Significant SNP associations were explored to find possible candidate genes related to physiological processes. Several of the significant markers detected were mapped onto functional candidate genes including ARFGAP3, CLSTN2 and DPYD for REA; OSBPL3 and SUDS3 for BFT; and RARRES1 and VEPH1 for RFT. The physiological pathway related to lipid metabolism (CLSTN2, OSBPL3, RARRES1 and VEPH1) was identified. The significant markers within previously reported QTLs reinforce the importance of the genomic regions, and the other loci offer candidate genes that have not been related to carcass traits in previous investigations.     

Genetic parameter estimation and gene association analyses for meat quality traits in open-air free-range Iberian pigs

Meat quality of Iberian pigs is defined by the combination of their genetic characteristics and the particular production system. To carry out a genetic analysis of the main meat quality traits, we estimated their heritabilities, genetic correlations and the association effects of 32 selected SNPs of 12 candidate genes. A total of ten traits were measured in longissimus dorsi samples from 1,199 Iberian pigs fattened in the traditional free-range system: water holding capacity (thawing, cooking and centrifuge force water losses), instrumental colour (lightness L*, redness a* and yellowness b*), myoglobin content, shear force on cooked meat, and shear force and maximum compression force on dry-cured loin. Estimated heritability values were low to moderate (0.01 to 0.43) being the lowest for L* and the highest for cooking loss. Strong genetic correlations between water holding capacity traits (0.93 to 0.96) and between myoglobin content and a* (0.94) were observed. The association analyses revealed 19 SNPs significantly associated with different traits. Consistent and strong effects were observed between PRKAG3 SNPs (rs319678464G > C and rs330427832C > T), MYH3_rs81437544T > C, CASP3_rs319658214G > T and CTSL_rs332171512A > G and water losses. Also for CAPN1_rs81358667G > A and CASP3_rs319658214G > T and shear force. The SNPs mapping on PRKAG3 showed the highest effects on Minolta colour traits. Genotyping of these SNPs could be useful for the selection of Iberian young boars with similar estimated breeding values for productive traits.     

Association of SNPs in GnRH gene with sperm quality traits of Chinese water buffalo

Hypothalamic gonadotropin-releasing hormone (GnRH) controls the activity of hypothalamic–pituitary–gonadal axis and plays a key role in the reproductive performance of animals. In this study, five single nucleotide polymorphisms (SNPs), namely g.991T > C, g.1041T > C g.3424T > C, g.3462C > A and g.3463Inde A, were detected in the GnRH gene of 162 water buffaloes by Sanger sequencing. Each SNP was associated with more than two sperm quality traits of ejaculate volume, sperm concentration, post-thaw sperm motility and sperm abnormality. g.3424T > C and g.3462C > A were related to these four traits and had a remarkable effect on ejaculate volume. The three other SNPs were related to sperm concentration, post-thaw sperm motility and sperm abnormality. Moreover, six haplotypes (H1: TCCAI, H2: CTTC-, H3: TCCCI, H4: CTTA-, H5: CCTA- and H6: CTCC-) composed of five SNPs comprising seven different combined genotypes were generated by linkage disequilibrium analysis. Statistics followed by one-way ANOVA indicated that water buffaloes with the haplotype combination H1H1 had the highest genotypic frequency, and those with the H4H4 haplotype combination had the highest ejaculate volume. The sperm concentration of those with haplotype combination H1H5 was higher than that of the other genotypes. In summary, our study showed a remarkable association between the SNPs of GnRH and sperm quality traits of Chinese water buffalo.     

Comparison of the effects explained by variations in the bovine PLAG1 and NCAPG genes on daily body weight gain,linear skeletal measurements and carcass traits in Japanese Black steers from a progeny testing program

The c.1326T>G single nucleotide polymorphism (SNP) in the NCAPG gene, which leads to an amino acid change of Ile442 to Met442, was previously identified as a candidate causative variation for a bovine carcass weight quantitative trait loci (QTL) on chromosome 6, which was associated with linear skeletal measurement gains and daily body weight gain at puberty. Recently, we identified the stature quantitative trait nucleotides (QTNs) in the PLAG1‐CHCHD7 intergenic region as the causative variations for another carcass weight QTL on chromosome 14. This study aimed to compare the effects of the two QTL on growth and carcass traits using 768 Japanese Black steers from a progeny testing program and to determine whether a genetic interaction was present between them. The FJX_250879 SNP representing the stature QTL was associated with linear skeletal measurements and average daily body weight gain at early and late periods during adolescence. A genetic interaction between FJX_250879 and NCAPG c.1326T>G was detected only for body and rump lengths. Both were associated with increased carcass weight and Longissimus muscle area, and NCAPG c.1326T>G was also associated with reduced subcutaneous fat thickness and increased carcass yield estimate. These results will provide useful information to improve carcass weight in Japanese Black cattle.     

Single nucleotide polymorphisms associated with thermoregulation in lactating dairy cows exposed to heat stress

Dairy cows with increased rectal temperature experience lower milk yield and fertility. Rectal temperature during heat stress is heritable, so genetic selection for body temperature regulation could reduce effects of heat stress on production. One aim of the study was to validate the relationship between genotype and heat tolerance for single nucleotide polymorphisms (SNPs) previously associated with resistance to heat stress. A second aim was to identify new SNPs associated with heat stress resistance. Thermotolerance was assessed in lactating Holsteins during the summer by measuring rectal temperature (a direct measurement of body temperature regulation; n = 435), respiration rate (an indirect measurement of body temperature regulation, n = 450) and sweating rate (the major evaporative cooling mechanism in cattle, n = 455). The association between genotype and thermotolerance was evaluated for 19 SNPs previously associated with rectal temperature from a genomewide analysis study (GWAS), four SNPs previously associated with change in milk yield during heat stress from GWAS, 2 candidate gene SNPs previously associated with rectal temperature and respiration rate during heat stress (ATPA1A and HSP70A) and 66 SNPs in genes previously shown to be associated with reproduction, production or health traits in Holsteins. For SNPs previously associated with heat tolerance, regions of BTA4, BTA6 and BTA24 were associated with rectal temperature; regions of BTA6 and BTA24 were associated with respiration rate; and regions of BTA5, BTA26 and BTA29 were associated with sweating rate. New SNPs were identified for rectal temperature (n = 12), respiration rate (n = 8) and sweating rate (n = 3) from among those previously associated with production, reproduction or health traits. The SNP that explained the most variation were PGR and ASL for rectal temperature, ACAT2 and HSD17B7 for respiration rate, and ARL6IP1 and SERPINE2 for sweating rate. ARL6IP1 was associated with all three thermotolerance traits. In conclusion, specific genetic markers responsible for genetic variation in thermoregulation during heat stress in Holsteins were identified. These markers may prove useful in genetic selection for heat tolerance in Holstein cattle.