Whole‐genome resequencing to identify candidate genes for the QTL for oleic acid percentage in Japanese Black cattle

In our previous study, we detected a QTL for the oleic acid percentage (C18:1) on BTA9 in Japanese Black cattle through a genome‐wide association study (GWAS). In this study, we performed whole‐genome resequencing on eight animals with higher and lower C18:1 to identify candidate polymorphisms for the QTL. A total of 39,658 polymorphisms were detected in the candidate region, which were narrowed to 1993 polymorphisms within 23 genes based on allele differences between the high and low C18:1 groups. We subsequently selected three candidate genes, that is, CYB5R4, MED23, and VNN1, among the 23 genes based on their function in fatty acid metabolism. In each candidate gene, three SNPs, that is, CYB5R4 c.*349G > T, MED23 c.3700G > A, and VNN1 c.197C > T, were selected as candidate SNPs to verify their effect on C18:1 in a Japanese Black cattle population (n = 889). The statistical analysis showed that these SNPs were significantly associated with C18:1 (p < 0.05), suggesting that they were candidates for the QTL. In conclusion, we successfully narrowed the candidates for the QTL by detecting possible polymorphisms located within the candidate region. It is expected that the responsible polymorphism can be identified by demonstrating their effect on the gene's function.     

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.     

A genome‐wide association study for fat‐related traits computed by image analysis in Japanese Black cattle

The objective of this study was to identify genomic regions associated with fat‐related traits using a Japanese Black cattle population in Hyogo. From 1836 animals, those with high or low values were selected on the basis of corrected phenotype and then pooled into high and low groups (n = 100 each), respectively. DNA pool‐based genome‐wide association study (GWAS) was performed using Illumina BovineSNP50 BeadChip v2 with three replicate assays for each pooled sample. GWAS detected that two single nucleotide polymorphisms (SNPs) on BTA7 (ARS‐BFGL‐NGS‐35463 and Hapmap23838‐BTA‐163815) and one SNP on BTA12 (ARS‐BFGL‐NGS‐2915) significantly affected fat percentage (FAR). The significance of ARS‐BFGL‐NGS‐35463 on BTA7 was confirmed by individual genotyping in all pooled samples. Moreover, association analysis between SNP and FAR in 803 Japanese Black cattle revealed a significant effect of SNP on FAR. Thus, further investigation of these regions is required to identify FAR‐associated genes and mutations, which can lead to the development of DNA markers for marker‐assisted selection for the genetic improvement of beef quality.     

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.     

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.     

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.     

An imputation‐based genome‐wide association study on traits related to male reproduction in a White Duroc × Erhualian F2 population

Boar reproductive traits are economically important for the pig industry. Here we conducted a genome‐wide association study (GWAS) for 13 reproductive traits measured on 205 F₂ boars at day 300 using 60 K single nucleotide polymorphism (SNP) data imputed from a reference panel of 1200 pigs in a White Duroc × Erhualian F₂ intercross population. We identified 10 significant loci for seven traits on eight pig chromosomes (SSC). Two loci surpassed the genome‐wide significance level, including one for epididymal weight around 60.25 Mb on SSC7 and one for semen temperature around 43.69 Mb on SSC4. Four of the 10 significant loci that we identified were consistent with previously reported quantitative trait loci for boar reproduction traits. We highlighted several interesting candidate genes at these loci, including APN, TEP1, PARP2, SPINK1 and PDE1C. To evaluate the imputation accuracy, we further genotyped nine GWAS top SNPs using PCR restriction fragment length polymorphism or Sanger sequencing. We found an average of 91.44% of genotype concordance, 95.36% of allelic concordance and 0.85 of r² correlation between imputed and real genotype data. This indicates that our GWAS mapping results based on imputed SNP data are reliable, providing insights into the genetic basis of boar reproductive traits.     

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.     

Genome‐wide association study for carcass traits,fatty acid composition,chemical composition,sugar, and the effects of related candidate genes in Japanese Black cattle

We performed a genome‐wide association study (GWAS) and candidate gene analysis to: (i) evaluate the effectiveness of the GWAS in our small population by performing GWAS for carcass weight (CW) and fatty acid composition; (ii) detect novel candidate regions affecting non‐CW carcass traits, chemical composition and sugar; and (iii) evaluate the association of the candidate genes previously detected in CW and fatty acid composition with other economically important traits. A total of 574 Japanese Black cattle and 40 657 Single nucleotide polymorphisms were used. In addition, candidate gene analyses were performed to evaluate the association of three CW‐related genes and two fatty acid‐related genes with carcass traits, fatty acid composition, chemical composition and sugar. The significant regions with the candidate genes were detected for CW and fatty acid composition, and these results showed that a significant region would be detectable despite the small sample size. The novel candidate regions were detected on BTA23 for crude protein and on BTA19 for fructose. CW‐related genes associated with the rib‐eye area and fatty acid composition were identified, and fatty acid‐related genes had no relationship with other traits. Moreover, the favorable allele of CW‐related genes had an unfavorable effect on fatty acid composition.     

Genome‐wide association study of cell‐mediated immune response in chicken

Cell‐mediated immunity (CMI) causes the intracellular destruction of the antigen or elimination of the host cell to make animals resistant against exogenous antigens and cancers. In this study, a genome‐wide association study (GWAS) was carried out to identify genomic regions associated with CMI in chicken using chicken 60k high‐density single nucleotide polymorphism (SNP) array. Genomic relationships were taken into account to adjust for population structure. In order to account for multiple testing, chromosome‐wise false discovery rate was controlled at 5% and 10% levels. Moreover, a comparison of the power of fixed and mixed linear models based on genomic inflation factor was carried out. Mixed linear model (MLM) had better inflation rate, and therefore the results from MLM were used for subsequent analysis. Three significantly associated SNPs (FDR < 0.05) on chromosome 24 and linkage group E22C19W28_E50C23, and three suggestively associated SNPs (FDR < 0.1) on chromosome 1, 5 and 16 were identified. Pathway analysis showed that two biological pathways, which are related to immune response, were strongly associated with the candidate genes surrounding identified SNPs, and their influences were mostly on antigen processing and presentation, and cellular structure.     

Application of highly differentiated SNPs between Japanese Black and Holstein to a breed assignment test between Japanese Black and F1 (Japanese Black x Holstein) and Holstein

Two taurine breeds, Japanese Black and Holstein, established from geographically distant origins and selected for different uses, beef and dairy, were extensively genotyped using a genome‐wide single nucleotide polymorphism (SNP) chip with more than 1000 animals of each breed. The genetic structure was examined by principal component analysis, in which the first principal component clearly separated the two breeds and explained more than 15% of the variance. Highly differentiated SNPs were detected throughout the genome, some of which were clustered within small regions on BTA4 (79.2–79.7 Mb, Btau4.0) and BTA26 (22.2–23.6 Mb). A breed assignment test was developed using 18 highly differentiated SNPs to distinguish Japanese Black from F₁ (Japanese Black × Holstein) and Holstein. The error rate that an F₁ or Holstein animal is misjudged as Japanese Black was expected to be < 0.8%, while the error rate that a Japanese Black animal is misjudged as F₁ or Holstein was expected to be < 0.001%. This test provides a reliable and powerful method to detect breed label falsification in retail beef.     

Genome‐wide association study for feedlot average daily gain in Nellore cattle (Bos indicus)

The genome‐wide association study (GWAS) results are presented for average daily gain (ADG) in Nellore cattle. Phenotype of 720 male Bos indicus animals with information of ADG in feedlots and 354 147 single‐nucleotide polymorphisms (SNPs) obtained from a database added by information from Illumina Bovine HD (777 962 SNPs) and Illumina BovineSNP50 (54 609) by imputation were used. After quality control and imputation, 290 620 SNPs remained in the association analysis, using R package Genome‐wide Rapid Association using Mixed Model and Regression method GRAMMAR‐Gamma. A genomic region with six significant SNPs, at Bonferroni‐corrected significance, was found on chromosome 3. The most significant SNP (rs42518459, BTA3: 85849977, p = 9.49 × 10^?8) explained 5.62% of the phenotypic variance and had the allele substitution effect of ?0.269 kg/day. Important genes such as PDE4B, LEPR, CYP2J2 and FGGY are located near this region, which is overlapped by 12 quantitative trait locus (QTLs) described for several production traits. Other regions with markers with suggestive effects were identified in BTA6 and BTA10. This study showed regions with major effects on ADG in Bos indicus in feedlots. This information may be useful to increase the efficiency of selecting this trait and to understand the physiological processes involved in its regulation.     

Genome‐wide association studies identified variants for taurine concentration in Japanese Black beef

We performed genome‐wide association studies (GWAS) using the BovineSNP50 array to detect significant single nucleotide polymorphisms (SNPs) that may affect the concentration of 22 free amino acids and three peptides in Japanese Black beef cattle. A total of 574 Japanese Black cattle and 40,657 SNPs from the array were used for this study. Genome‐wide significant SNPs were detected for β‐alanine (three SNPs on chromosomes 22 and 29) and taurine (26 SNPs on chromosome 22). Importantly, the top two SNPs for taurine were highly significant (p = 6.2 × 10^?21), and the frequency of the increase‐concentration allele (Q) for taurine was found to be 0.73. The Q allele frequency of this population was similar to that of the other unrelated Japanese Black cattle, but different from that of the other breeds. In addition, the significant SNPs were not associated with carcass traits or fatty acid compositions. Interestingly, the top three of the four most significant SNPs for taurine were located near solute carrier family 6, member 6 (SLC6A6), which is a membrane transporter for taurine. We also found two associated variants in the 5′‐upstream region of SLC6A6; however, they were less significantly associated than the SNPs from the BovineSNP50 array.     

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.     

Genome‐wide association studies for two exterior traits in Chinese Dongxiang spotted pigs

Pigs show extensive variation in exterior appearance. This variation has been explored as one of the selection targets to form breeding features in the pig industry. In this study, a customized Affymetrix Axiom Pig1.4M array plate was used to conduct genome‐wide association studies (GWAS) for two exterior traits—spotted coat color (SC) and facial type (FT)—in Chinese Dongxiang spotted pigs. Two and 12 single nucleotide polymorphisms (SNPs) were identified to be associated with SC and FT at the 5% genome‐wide significant level, respectively. Of these SNPs, two associated with SC were located around the DACH1 gene on Sus Scrofa chromosome (SSC)11. Eleven SNPs associated with FT were located within a 7.68 Mb region (29.03–36.71 Mb) on SSC7. Based on the GWAS results and the biological functions of genes, we highlight EDNRB as a candidate gene for SC and HMGA1 and RPS10 as the potential genes affecting facial variation. The findings contribute to the final characterization of causative genes and mutations underlying the effects of the SSC7 and SSC11 loci and improve our understanding of the genetic basis of phenotypic variation in Chinese indigenous pigs.     

Genetic differences in oestrous signs and oestrogen metabolism‐related genes between Chinese Mi and European Landrace‐Large White pigs

Oestrous signs affect timely mating and reproductive efficiency in swine breeding herds. To study the genetic difference of oestrous signs between Chinese and European pigs, 100 Landrace‐Large White (LLW) cross gilts and 50 Chinese Mi gilts were assessed for oestrous signs and the concentrations of serum estradiol‐17β and progesterone were determined. The genotype of 39 single nucleotide polymorphisms (SNPs) in 11 oestrogen metabolism and function‐related genes was determined by Sequenom iPLEX platform. Compared with LLW gilts, Mi gilts had longer time of standing reflex (p < .001), higher scores of vulva reddening (p = .001) and greater serum estradiol‐17β concentration (p < .01). Gilts with greater serum estradiol‐17β concentrations also had greater (p < .05) scores for oestrous signs. Genetic polymorphisms of nine genes in oestrogen metabolism pathways had significant differences (p < .05) between LLW and Mi gilts. There were three and six haploblocks of SNPs in LLW and Mi, respectively. Compared with LLW, the distribution of haplotypes was more centralized in Mi pigs. Genetic polymorphisms of oestrogen metabolism‐related genes have considerable differences between Chinese Mi and European LLW pigs. Because of the important roles of oestrogen during the oestrus, some genes of oestrogen metabolism pathway could be considered as candidate genes for oestrous signs.