Effect of DNA markers on the fertility traits of Japanese Black cattle for improving beef quantity and quality.

Carcass traits have been efficiently improved by recent selection using DNA markers in beef cattle. Additionally, DNA markers might have an effect on other traits such as fertility traits; therefore attention should also be paid to such pleiotropic effects. However, the effects of the markers on both carcass and fertility traits have never been evaluated in the same population, since they are generally measured in different populations. The objective in the current study was to discuss effectiveness of DNA markers developed for carcass traits through investigation of their effects on carcass and fertility traits in a population. We genotyped six markers SCD V293A, FASN g.841G>C, PLAG1 g.49066C>G, NCAPG I442M, DGAT1 K232A, and EDG1 g.1471620G>T in a Japanese Black cattle population ( n = 515 ). To investigate their effects on carcass and fertility traits, we performed statistical analysis (ANOVA and the Tukey-Kramer honestly significant difference (HSD) test). In the results, three of six markers, SCD V293A, NCAPG I442M, and EGD1 g.1471620G>T, were significantly associated with both carcass and fertility traits. Remarkably, the same allele for each marker had positive effects on both traits, suggesting that we would be able to simultaneously improve them using these markers in this population. However, previous studies reported that the effects of DNA markers could differ among populations. Therefore, it is necessary to confirm the effect of the marker in each population before it is used for improvement.

IARS mutation causes prenatal death in Japanese Black cattle.

Isoleucyl-tRNA synthetase (IARS) c.235G > C (p.V79L) is a causative mutation for a recessive disease called IARS disorder in Japanese black cattle. The disease is involved in weak calf syndrome and is characterized by low birth weight, weakness and poor suckling. The gestation period is often slightly extended, implying that intrauterine growth is retarded. In a previous analysis of 2597 artificial insemination (AI) procedures, we suggested that the IARS mutation might contribute toward an increase in the incidence of prenatal death. In this study, we extended this analysis to better clarify the association between the IARS mutation and prenatal death. The IARS genotypes of 92 animals resulting from crosses between carrier (G/C) × G/C were 27 normal (G/G), 55 G/C and 10 affected animals (C/C) (expected numbers: 23, 46 and 23, respectively). Compared to the expected numbers, there were significantly fewer affected animals in this population (P < 0.05), suggesting that more than half of the affected embryos died prenatally. When the number of AI procedures examined was increased to 11 580, the frequency of re-insemination after G/C × G/C insemination was significantly higher at 61-140 days (P < 0.001). The findings suggested that the homozygous IARS mutation not only causes calf death, but also embryonic or fetal death.

Allelic frequencies and association with carcass traits of six genes in local subpopulations of Japanese Black cattle.

Marker-assisted selection (MAS) is expected to accelerate the genetic improvement of Japanese Black cattle. However, verification of the effects of the genes for MAS in different subpopulations is required prior to the application of MAS. In this study, we investigated the allelic frequencies and genotypic effects for carcass traits of six genes, which can be used in MAS, in eight local subpopulations. These genes are SCD, FASN and SREBP1, which are associated with the fatty acid composition of meat, and NCAPG, MC1R and F11, which are associated with carcass weight, coat color and blood coagulation abnormality, respectively. The frequencies of desirable alleles of SCD and FASN were relatively high and that of NCAPG was relatively low, and NCAPG was significantly associated with several carcass traits, including carcass weight. The proportions of genotypic variance explained by NCAPG to phenotypic variance were 4.83 for carcass weight. We thus confirmed that NCAPG is a useful marker for selection of carcass traits in these subpopulations. In addition, we found that the desirable alleles of six genes showed no negative effects on carcass traits. Therefore, selection using these genes to improve target traits should not have negative impacts on carcass traits.

Genetic variants in the upstream region of activin receptor IIA are associated with female fertility in Japanese Black cattle.

BACKGROUND: Female fertility, a fundamental trait required for animal reproduction, has gradually declined in the last 2 decades in Japanese Black cattle. To identify associated genetic variants in Japanese Black cattle, we evaluated female fertility as a metric to describe the average inverse of the number of artificial inseminations required for conception from the first through the fourth parity (ANAI4) and conducted a genome-wide association study (GWAS) using 430 animals with extreme ANAI4 values from 10,399 animals. RESULTS: We found that 2 variants, namely a single-nucleotide polymorphisms (SNP; g.48476925C > T) and a 3-bp indel (g.48476943_48476946insGGC), in the upstream region of the activin receptor IIA gene (ACVR2A) were associated with ANAI4. ACVR2A transcripts from Japanese Black cattle of the Q haplotype, defined by the SNP and the 3-bp indel, with increased ANAI4 were 1.29-1.32-fold more abundant than q-derived transcripts. In agreement, reporter assay results revealed that the activity of the ACVR2A promoter was higher in reporter constructs with the Q haplotype than in those with the q haplotype by approximately 1.2 fold. Expression of exogenous ACVR2A induced dose-dependent increases of reporter activity from the follicle-stimulating hormone, beta polypeptide (FSHB) promoter in response to activin A in a pituitary gonadotrophic cell line. The findings suggested that sequence variations in the upstream region of ACVR2A with the Q haplotype increased ACVR2A transcription, which in turn induced FSHB expression. This association was replicated using a sample population size of 1,433 animals; the frequency of the Q haplotype was 0.39, and Q-to-q haplotype substitution resulted in an increase of 0.02 in terms of ANAI4. CONCLUSIONS: This GWAS identified variants in the upstream region of ACVR2A, which were associated with female fertility in Japanese Black cattle. The variants affected the level of ACVR2A mRNA expression, which could lead to an allelic imbalance. This association was replicated with a sample population of 1,433 animals. Thus, the results suggest that the Q haplotype could serve as a useful marker to select Japanese Black cattle with superior female fertility.

Non-synonymous FGD3 Variant as Positional Candidate for Disproportional Tall Stature Accounting for a Carcass Weight QTL (CW-3) and Skeletal Dysplasia in Japanese Black Cattle.

Recessive skeletal dysplasia, characterized by joint- and/or hip bone-enlargement, was mapped within the critical region for a major quantitative trait locus (QTL) influencing carcass weight; previously named CW-3 in Japanese Black cattle. The risk allele was on the same chromosome as the Q allele that increases carcass weight. Phenotypic characterization revealed that the risk allele causes disproportional tall stature and bone size that increases carcass weight in heterozygous individuals but causes disproportionately narrow chest width in homozygotes. A non-synonymous variant of FGD3 was identified as a positional candidate quantitative trait nucleotide (QTN) and the corresponding mutant protein showed reduced activity as a guanine nucleotide exchange factor for Cdc42. FGD3 is expressed in the growth plate cartilage of femurs from bovine and mouse. Thus, loss of FDG3 activity may lead to subsequent loss of Cdc42 function. This would be consistent with the columnar disorganization of proliferating chondrocytes in chondrocyte-specific inactivated Cdc42 mutant mice. This is the first report showing association of FGD3 with skeletal dysplasia.

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.

Variants in the 3' UTR of general Transcription factor IIF, polypeptide 2 affect female calving efficiency in Japanese Black cattle.

BACKGROUND: Calving efficiency can be described as the measure of a cow's ability to produce viable offspring within a specific period of time. This trait is crucial in beef cattle because calves are necessary both for the production of beef and for heifer replacements. Recently, the number of calves produced at 4 years of age (NCP4) has been used to evaluate the calving efficiency of Japanese Black cattle. To identify variants associated with calving efficiency in Japanese Black cattle, we conducted a genome-wide association study (GWAS) using 688 animals with extreme NCP4 values selected from 15,225 animals. RESULTS: We identified genetic variants on bovine chromosome 12 (BTA12) that were associated with NCP4. The General Transcription Factor IIF, polypeptide 2 (GTF2F2), located in the 132 kbp-associated region, proved to be in strong linkage disequilibrium. We found 15 associated variants in the promoter and the 3' UTR regions. Consistent with this finding, transcripts of GTF2F2 derived from the haplotype (Q) with the increased number of calves were 1.33-fold more abundant than q-derived transcripts. Furthermore, luciferase assays revealed that the activity of the 3' UTR, a region that includes nine SNPs, was higher in constructs with the Q haplotype than in those with the q haplotype by approximately 1.35-fold. In contrast, the activity of the promoter region did not differ between haplotypes. The association was replicated in an independent sample of 827 animals that were randomly selected from the remainder of the cohort from the same farms used in the GWAS. In the replicated population, the frequency of the Q haplotype is 0.313, and this haplotype accounts for 2.69% of the total phenotypic variance. The effect of the Q to q haplotype substitution on NCP4 was 0.054 calves. These findings suggest that variants in the 3' UTR of GTF2F2 affect the level of GTF2F2 mRNA, which is associated with calving efficiency. CONCLUSIONS: This GWAS has identified variants in the 3' UTR of GTF2F2 that were associated with the NCP4 of Japanese Black cattle, and this association was validated in an independent sample. The Q haplotype will be immediately useful in improving the calving efficiency of Japanese Black cattle.

Mapping and exome sequencing identifies a mutation in the IARS gene as the cause of hereditary perinatal weak calf syndrome.

We identified an IARS (isoleucyl-tRNA synthetase) c.235G>C (p.Val79Leu) substitution as the causative mutation for neonatal weakness with intrauterine growth retardation (perinatal weak calf syndrome). In Japanese Black cattle, the syndrome was frequently found in calves sired by Bull A. Hence, we employed homozygosity mapping and linkage analysis. In order to identify the perinatal weak calf syndrome locus in a 4.04-Mb region of BTA 8, we analysed a paternal half-sibling family with a BovineSNP50 BeadChip and microsatellites. In this critical region, we performed exome sequencing to identify a causative mutation. Three variants were detected as possible candidates for causative mutations that were predicted to disrupt the protein function, including a G>C (p.Val79Leu) mutation in IARS c.235. The IARS c.235G>C mutation was not a homozygous risk allele in the 36 healthy offspring of Bull A. Moreover, the IARS Val79 residue and its flanking regions were evolutionarily and highly conserved. The IARS mutant (Leu79) had decreased aminoacylation activity. Additionally, the homozygous mutation was not found in any of 1526 healthy cattle. Therefore, we concluded that the IARS c.235G>C mutation was the cause of hereditary perinatal weak calf syndrome.

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 assessment of genetic diversity within and among the eight subpopulations of Japanese Black cattle using 52 microsatellite markers.

Japanese Black cattle are at risk for genetic homogeneity due to intensive use of a few sires. Therefore, assessment of the actual genetic diversity of this breed is important for future breeding plans. In the present study, we investigated the genetic diversity within and among eight subpopulations of Japanese Black cattle using 52 microsatellite markers. The parameters for genetic diversity of Japanese Black cattle were comparable to those of other cattle breeds, suggesting that the relatively high genetic diversity of the breed. However, upon comparison among the eight subpopulations, the Hyogo subpopulation showed markedly low genetic diversity. The results of the pairwise FST values, phylogenetic network and structure analysis indicated that the Hyogo population has remarkably high level of genetic differentiation from other populations, while Yamagata, Niigata, Hiroshima and Kagawa populations have low levels of genetic differentiation. Furthermore, multidimensional scaling plots indicated that individuals in some subpopulations were separated from individuals in the other subpopulations. We conclude that while the overall genetic diversity of Japanese Black cattle is still maintained at a relatively high level, that of a particular subpopulation is significantly reduced, and therefore the effective population size of the breed needs to be controlled by correct mating strategies.

BLV-CoCoMo-qPCR: a useful tool for evaluating bovine leukemia virus infection status.

BACKGROUND: Bovine leukemia virus (BLV) is associated with enzootic bovine leukosis, which is the most common neoplastic disease of cattle. BLV infects cattle worldwide, imposing a severe economic impact on the dairy cattle industry. Recently, we developed a new quantitative real-time polymerase chain reaction (PCR) method using Coordination of Common Motifs (CoCoMo) primers to measure the proviral load of known and novel BLV variants in BLV-infected animals. Indeed, the assay was highly effective in detecting BLV in cattle from a range of international locations. This assay enabled us to demonstrate that proviral load correlates not only with BLV infection capacity as assessed by syncytium formation, but also with BLV disease progression. In this study, we compared the sensitivity of our BLV-CoCoMo-qPCR method for detecting BLV proviruses with the sensitivities of two real-time PCR systems, and also determined the differences of proviral load with serotests. RESULTS: BLV-CoCoMo-qPCR was found to be highly sensitive when compared with the real-time PCR-based TaqMan MGB assay developed by Lew et al. and the commercial TaKaRa cycleave PCR system. The BLV copy number determined by BLV-CoCoMo-qPCR was only partially correlated with the positive rate for anti-BLV antibody as determined by the enzyme-linked immunosorbent assay, passive hemagglutination reaction, or agar gel immunodiffusion. This result indicates that, although serotests are widely used for the diagnosis of BLV infection, it is difficult to detect BLV infection with confidence by using serological tests alone. Two cattle were experimentally infected with BLV. The kinetics of the provirus did not precisely correlate with the change in anti-BLV antibody production. Moreover, both reactions were different in cattle that carried different bovine leukocyte antigen (BoLA)-DRB3 genotypes. CONCLUSIONS: Our results suggest that the quantitative measurement of proviral load by BLV-CoCoMo-qPCR is useful tool for evaluating the progression of BLV-induced disease. BLV-CoCoMo-qPCR allows us to monitor the spread of BLV infection in different viewpoint compared with classical serotest.

The diversity of bovine MHC class II DRB3 and DQA1 alleles in different herds of Japanese Black and Holstein cattle in Japan.

In cattle, bovine leukocyte antigens (BoLAs) have been extensively used as markers for bovine diseases and immunological traits. In this study, we sequenced alleles of the BoLA class II loci, BoLA-DRB3 and BoLA-DQA1, from 650 Japanese cattle from six herds [three herds (507 animals) of Japanese Black cattle and three herds (143 animals) of Holstein cattle] using polymerase chain reaction-sequence-based typing (PCR-SBT) methods. We identified 26 previously reported distinct DRB3 alleles in the two populations: 22 in Japanese Black and 17 in Holstein. The number of DRB3 alleles detected in each herd ranged from 9 to 20. Next, we identified 15 previously reported distinct DQA1 alleles: 13 in Japanese Black and 10 in Holstein. The number of alleles in each herd ranged from 6 to 10. Thus, allelic divergence is significantly greater for DRB3 than for DQA1. A population tree on the basis of the frequencies of the DRB3 and DQA1 alleles showed that, although the genetic distance differed significantly between the two cattle breeds, it was closely related within the three herds of each breed. In addition, Wu-Kabat variability analysis indicated that the DRB3 gene was more polymorphic than the DQA1 gene in both breeds and in all herds, and that the majority of the hypervariable positions within both loci corresponded to pocket-forming residues. The DRB3 and DQA1 heterozygosity for both breeds within each herd were calculated based on the Hardy-Weinberg equilibrium. Only one Japanese Black herd showed a significant difference between the expected and observed heterozygosity at both loci. This is the first report presenting a detailed study of the allelic distribution of BoLA-DRB3 and -DQA1 genes in Japanese Black and Holstein cattle from different farms in Japan. These results may help to develop improved livestock breeding strategies in the future.

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method for mtDNA typing in hokkaido brown bear (Ursus arctos yesoensis).

To develop an easy method of typing of mitochondrial DNA (mtDNA) in Hokkaido brown bears (Ursus arctos yesoensis), the PCR-RFLP technique was improved using four restriction enzymes: Mbo 1, Cfr 13 I, TspE 1, and Fok 1. This approach identified seven groups of mtDNA haplotypes, HB1/2/5-7, HB 3, HB4, HB8/9, HB10/11, HB12 and HB13 from 102 brown bears of northern, central and eastern Hokkaido.