Accuracy of imputation of single nucleotide polymorphism marker genotypes from low-density panels in Japanese Black cattle.

Using target and reference fattened steer populations, the performance of genotype imputation using lower-density marker panels in Japanese Black cattle was evaluated. Population imputation was performed using BEAGLE software. Genotype information for approximately 40,000 single nucleotide polymorphism (SNP) markers by Illumina BovineSNP50 BeadChip was available, and imputation accuracy was assessed based on the average concordance rates of the genotypes, varying equally spaced SNP densities, and the number of individuals in the reference population. Two additional statistics were also calculated as indicators of imputation performance. The concordance rates tended to be lower for SNPs with greater minor allele frequencies, or those located near the ends of the chromosomes. Longer autosomes yielded greater imputation accuracies than shorter ones. When SNPs were selected based on linkage disequilibrium information, relative imputation accuracy was slightly improved. When 3000 and 10,000 equally spaced SNPs were used, the imputation accuracies were greater than 90% and approximately 97%, respectively. These results indicate that combining genotyping using a lower-density SNP chip with genotype imputation based on a population of individuals genotyped using a higher-density SNP chip is a cost-effective and valid approach for genomic prediction.

PLAG1 and NCAPG-LCORL in livestock.

A recent progress on stature genetics has revealed simple genetic architecture in livestock animals in contrast to that in humans. PLAG1 and/or NCAPG-LCORL, both of which are known as a locus for adult human height, have been detected for association with body weight/height in cattle and horses, and for selective sweep in dogs and pigs. The findings indicate a significant impact of these loci on mammalian growth or body size and usefulness of the natural variants for selective breeding. However, association with an unfavorable trait, such as late puberty or risk for a neuropathic disease, was also reported for the respective loci, indicating an importance to discriminate between causality and association. Here I review the recent findings on quantitative trait loci (QTL) for stature in livestock animals, mainly focusing on the PLAG1 and NCAPG-LCORL loci. I also describe our recent efforts to identify the causative variation for the third major locus for carcass weight in Japanese Black cattle.

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.

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.

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.

Application of highly differentiated SNPs between Japanese Black and Holstein to a breed assignment test between Japanese Black and F(1) (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(1) (Japanese Black × Holstein) and Holstein. The error rate that an F(1) 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(1) 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 identified three major QTL for carcass weight including the PLAG1-CHCHD7 QTN for stature in Japanese Black cattle.

BACKGROUND: Significant quantitative trait loci (QTL) for carcass weight were previously mapped on several chromosomes in Japanese Black half-sib families. Two QTL, CW-1 and CW-2, were narrowed down to 1.1-Mb and 591-kb regions, respectively. Recent advances in genomic tools allowed us to perform a genome-wide association study (GWAS) in cattle to detect associations in a general population and estimate their effect size. Here, we performed a GWAS for carcass weight using 1156 Japanese Black steers. RESULTS: Bonferroni-corrected genome-wide significant associations were detected in three chromosomal regions on bovine chromosomes (BTA) 6, 8, and 14. The associated single nucleotide polymorphisms (SNP) on BTA 6 were in linkage disequilibrium with the SNP encoding NCAPG Ile442Met, which was previously identified as a candidate quantitative trait nucleotide for CW-2. In contrast, the most highly associated SNP on BTA 14 was located 2.3-Mb centromeric from the previously identified CW-1 region. Linkage disequilibrium mapping led to a revision of the CW-1 region within a 0.9-Mb interval around the associated SNP, and targeted resequencing followed by association analysis highlighted the quantitative trait nucleotides for bovine stature in the PLAG1-CHCHD7 intergenic region. The association on BTA 8 was accounted for by two SNP on the BovineSNP50 BeadChip and corresponded to CW-3, which was simultaneously detected by linkage analyses using half-sib families. The allele substitution effects of CW-1, CW-2, and CW-3 were 28.4, 35.3, and 35.0 kg per allele, respectively. CONCLUSION: The GWAS revealed the genetic architecture underlying carcass weight variation in Japanese Black cattle in which three major QTL accounted for approximately one-third of the genetic variance.

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

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

Dissection of genetic factors modulating fetal growth in cattle indicates a substantial role of the non-SMC condensin I complex, subunit G (NCAPG) gene.

The increasing evidence of fetal developmental effects on postnatal life, the still unknown fetal growth mechanisms impairing offspring generated by somatic nuclear transfer techniques, and the impact on stillbirth and dystocia in conventional reproduction have generated increasing attention toward mammalian fetal growth. We identified a highly significant quantitative trait locus (QTL) affecting fetal growth on bovine chromosome 6 in a specific resource population, which was set up by consistent use of embryo transfer and foster mothers and, thus, enabled dissection of fetal-specific genetic components of fetal growth. Merging our data with results from other cattle populations differing in historical and geographical origin and with comparative data from human whole-genome association mapping suggests that a nonsynonymous polymorphism in the non-SMC condensin I complex, subunit G (NCAPG) gene, NCAPG c.1326T>G, is the potential cause of the identified QTL resulting in divergent bovine fetal growth. NCAPG gene expression data in fetal placentomes with different NCAPG c.1326T>G genotypes, which are in line with recent results about differential NCAPG expression in placentomes from studies on assisted reproduction techniques, indicate that the NCAPG locus may give valuable information on the specific mechanisms regulating fetal growth in mammals.

Cross-breed comparisons identified a critical 591-kb region for bovine carcass weight QTL (CW-2) on chromosome 6 and the Ile-442-Met substitution in NCAPG as a positional candidate.

BACKGROUND: Growth-related traits have been mapped on bovine chromosome 6 (BTA 6) in various bovine breed populations. We previously mapped a significant quantitative trait locus (QTL) for carcass and body weight (CW-2) between 38 and 55 cM on BTA 6 using a Japanese Black half-sib family. Additional QTL mapping studies detected four QTL for body or carcass weight that overlapped with CW-2 in Japanese Black and Japanese Brown half-sib families. To map the region in greater detail, we applied cross-breed comparisons of haplotypes that have been shown to be powerful in canine. RESULTS: We used 38 microsatellite markers to search for a shared Q (increasing carcass and/or body weight) haplotype within the 17-cM CW-2 region among five sires. Linkage disequilibrium mapping using maternal alleles of the offspring showed that an 815-kb shared Q haplotype was associated with body or carcass weight in both breeds. The addition of 43 single nucleotide polymorphism (SNP) markers narrowed the region to 591 kb containing 4 genes. The SNP changing Ile-442 to Met in NCAPG (chromosome condensation protein G) was significantly associated with carcass weight (p < 1.2 x 10-11) in a large Japanese Black population as well as in the five families. The Q allele of the SNP was also associated with a larger longissimus muscle area and thinner subcutaneous fat thickness in steers of all five families, indicating that the CW-2 locus is pleiotropic and favorable for marker-assisted selection of beef cattle. CONCLUSION: A 591-kb critical region for CW-2 was identified. The SNP changing Ile-442 to Met in NCAPG (chromosome condensation protein G) can be used as a positional candidate of CW-2 for marker-assisted selection.

Structural and functional analysis of the bovine Mx1 promoter.

The bovine Mx1 promoter region was found to contain 4 IFN-stimulated response elements (ISREs), 7 GC boxes, 2 IL-6 responsive elements, 2 NFκB-binding sites and 2 AP-1-binding sites. Among Holstein, Charolai, and Brahman breeds, 5 nucleotide substitutions were detected in the promoter region. After the Mx1 promoter region from Holstein had been constructed with pGL-basic expression vector, the transfected cells showed promoter activity after IFN induction. Several artificial deletion mutants were prepared to determine the important regulatory elements responsible for the promoter activity, and it was found that ISRE has a key function in IFN response. The proximal ISRE1 showed potential induction by IFN. Furthermore, the proximal GC boxes were found to be essential for IFN response in the bovine Mx1 promoter with the deletion mutants. In this case, the 2 GC boxes exhibited a synergistic activation in the IFN response. Mithramycin A, an agent that inhibits gene expression selectively by coating GC boxes, was used, and Mx mRNA expression in MDBK cells was suppressed by this chemical. Therefore, GC boxes were also shown to be essential for IFN response in the bovine Mx1 gene.

Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping.

To map quantitative trait loci (QTL) for growth and carcass traits in a purebred Japanese Black cattle population, we conducted multiple QTL analyses using 15 paternal half-sib families comprising 7860 offspring. We identified 40 QTL with significant linkages at false discovery rates of less than 0.1, which included 12 for intramuscular fat deposition called marbling and 12 for cold carcass weight or body weight. The QTL each explained 2%-13% of the phenotypic variance. These QTL included many replications and shared hypothetical identical-by-descent (IBD) alleles. The QTL for CW on BTA14 was replicated in five families with significant linkages and in two families with a 1% chromosome-wise significance level. The seven sires shared a 1.1-Mb superior Q haplotype as a hypothetical IBD allele that corresponds to the critical region previously refined by linkage disequilibrium mapping. The QTL for marbling on BTA4 was replicated in two families with significant linkages. The QTL for marbling on BTA6, 7, 9, 10, 20, and 21 and the QTL for body weight on BTA6 were replicated with 1% and/or 5% chromosome-wise significance levels. There were shared IBD Q or q haplotypes in the marbling QTL on BTA4, 6, and 10. The allele substitution effect of these haplotypes ranged from 0.7 to 1.2, and an additive effect between the marbling QTL on BTA6 and 10 was observed in the family examined. The abundant and replicated QTL information will enhance the opportunities for positional cloning of causative genes for the quantitative traits and efficient breeding using marker-assisted selection.

Identification of a 1.1-Mb region for a carcass weight QTL onbovine Chromosome 14.

We previously mapped a quantitative trait locus for carcass weight, designated Carcass Weight-1 (CW-1), to bovine Chromosome 14 using a purebred Wagyu pedigree based on progeny design analysis. To refine the critical region within 8.1 cM flanked by microsatellites BMS1941 and INRA094, we constructed a bacterial artificial chromosome (BAC) contig composed of 60 tiled BAC clones and prepared a high-density physical map including 80 microsatellites, of which 55 were developed in this study. We conducted linkage disequilibrium (LD) mapping in the CW-1 region with 47 microsatellites using paternal half-sib pedigrees whose sires exhibited homozygous CW-1 Q alleles in the region. The LD mapping study significantly narrowed the CW-1 locus to the 1.1-Mb region between microsatellites DIK7012 and DIK7020. Finally, we surveyed the 1.1-Mb-region genotypes of 1700 steers from 11 bulls having -/-, Q/-, or Q/Q alleles in the region, and we examined the effect of the CW-1 Q allele on carcass weight. The presence of the first Q increased carcass weight by 23.6 kg (95% confidence interval [CI], 17.6-29.5 kg), and the second Q increased carcass weight an additional 15.2 kg (95% CI, 10.7-19.7 kg). These results indicate the presence of a gene responsible for carcass weight within the 1.1-Mb region.

A comprehensive radiation hybrid map of the bovine genome comprising 5593 loci.

A bovine whole genome 7000-rad radiation hybrid (RH) panel, SUNbRH(7000-rad), was constructed to build a high-resolution RH map. The Shirakawa-USDA linkage map served as a scaffold to construct a framework map of 3216 microsatellites on which 2377 ESTs were ordered. The resulting RH map provided essentially complete coverage across the genome, with 1 cR7000 corresponding to 114 kb, and a cattle-human comparative map of 1716 bovine genes and sequences annotated in the human genome, which covered 79 and 72% of the bovine and human genomes, respectively. We then integrated the bovine RH and comparative maps with BAC fingerprint information in to construct a detailed, BAC-based physical map covering a reported 40-cM quantitative trait locus region for intramuscular fat or "marbling" on BTA 4. In summary, the new, high-resolution SUNbRH7000-rad, comparative, Shirakawa-USDA linkage, and BAC fingerprint maps provide a set of genomic tools for fine mapping regions of interest in cattle.

Abnormal development of nephrons in claudin-16-defective Japanese black cattle.

The kidneys of 37 Japanese Black calves aged 2 to 65 months diagnosed with Claudin 16 (CL-16) defect by the DNA-based test were examined pathologically. The animals exhibited clinical symptoms such as growth impairment, renal failure, overgrowth of hooves, and anemia at a young age. There was no correlation between the time of onset and age. Kidney weights relative to body weight were similar to those in normal animals, but both kidney net weights and size were reduced due to atrophy in animals that showed severe renal dysfunction. Histopathological examination of the kidneys showed reduction in the number of glomeruli, compensatory hypertrophy of glomeruli and tubules, and glomerular and tubular atrophy accompanied by interstitial fibrosis and lymphocytic infiltration. Glomeruli were clearly less in number in the kidneys of CL-16-defective animals than those of normal animals even in the cases with mild lesions. A small number of immature glomeruli and tubules were also detected, suggesting that there were fewer nephrons developed at birth in CL-16-defective animals. It was suggested that a defect of the CL-16 gene is involved in the "abnormal development of nephrons". Immunohistopathological examination of the kidneys showed that the epithelium of thick ascending limb of Henle was stained with anti-CL-16 antibody in the control animals, but not in the affected animals, suggesting a defect of CL-16 in the epithelium of renal tubules in the affected animals.

Integrating linkage and radiation hybrid mapping data for bovine chromosome 15.

BACKGROUND: Bovine chromosome (BTA) 15 contains a quantitative trait loci (QTL) for meat tenderness, as well as several breaks in synteny with human chromosome (HSA) 11. Both linkage and radiation hybrid (RH) maps of BTA 15 are available, but the linkage map lacks gene-specific markers needed to identify genes underlying the QTL, and the gene-rich RH map lacks associations with marker genotypes needed to define the QTL. Integrating the maps will provide information to further explore the QTL as well as refine the comparative map between BTA 15 and HSA 11. A recently developed approach to integrating linkage and RH maps uses both linkage and RH data to resolve a consensus marker order, rather than aligning independently constructed maps. Automated map construction procedures employing this maximum-likelihood approach were developed to integrate BTA RH and linkage data, and establish comparative positions of BTA 15 markers with HSA 11 homologs. RESULTS: The integrated BTA 15 map represents 145 markers; 42 shared by both data sets, 36 unique to the linkage data and 67 unique to RH data. Sequence alignment yielded comparative positions for 77 bovine markers with homologs on HSA 11. The map covers approximately 32% of HSA 11 sequence in five segments of conserved synteny, another 15% of HSA 11 is shared with BTA 29. Bovine and human order are consistent in portions of the syntenic segments, but some rearrangement is apparent. Comparative positions of gene markers near the meat tenderness QTL indicate the region includes separate segments of HSA 11. The two microsatellite markers flanking the QTL peak are between defined syntenic segments. CONCLUSIONS: Combining data to construct an integrated map not only consolidates information from different sources onto a single map, but information contributed from each data set increases the accuracy of the map. Comparison of bovine maps with well annotated human sequence can provide useful information about genes near mapped bovine markers, but bovine gene order may be different than human. Procedures to connect genetic and physical mapping data, build integrated maps for livestock species, and connect those maps to more fully annotated sequence can be automated, facilitating the maintenance of up-to-date maps, and providing a valuable tool to further explore genetic variation in livestock.

A comprehensive genetic map of the cattle genome based on 3802 microsatellites.

A microsatellite-based high-density genetic map facilitates for fine mapping of hereditary traits of interest, characterization of meiosis, and providing a foundation for physical map construction. Here, we developed a comprehensive genetic map on the basis of >880,000 genotypes across the USDA MARC cattle reference families with a potential genetic resolution of 0.8 cM at the 95% confidence level ( approximately 800 kb in the bovine genome). We incorporated 2325 microsatellites into the second-generation genetic map by linkage analysis based on sex-averaged two-point LOD scores (>3.0), of which 2293 were fine-mapped by multipoint linkage analysis. The new 3160-cM map comprised of 29 sex-averaged autosomal linkage groups and a sex-specific X-chromosome linkage group includes 3960 markers with 2389 positions, resulting in an average interval size of 1.4 cM. More than half (51%) of the total length of the map is covered with intervals of 2.0 cM or less, and the largest gap is a 10.2-cM interval on the X-linkage group. The new map should accelerate fine mapping and positional cloning of genes for genetic diseases and economically important traits in cattle, as well as related livestock species, such as sheep and goat.

Cloning and expression of type XII collagen isoforms during bovine adipogenesis.

In order to isolate candidate genes involved in bovine adipocyte differentiation, we have constructed a subtraction library from a clonal bovine intra-muscular pre-adipocyte (BIP) cell line using the suppression subtractive hybridization method. We have isolated a set of subtracted cDNA fragments whose respective mRNA levels are up-regulated during the adipogenic differentiation of BIP cells, and cloned cDNAs from a differentiated BIP-lambda ZAP II cDNA library. Two cDNA clones were highly homologous to the sequence of mouse and human type XII collagen alpha-1, determined by a BLAST homology search. As type XII collagen has been reported to have four types of splicing isoform, two clones were determined to be XII-1 and XII-2 splicing isoforms, respectively, because of a difference in the C-terminal NC1 domain. From the expression analysis of type XII collagen, the XIIA-2 isoform was mainly expressed in differentiated BIP cells and adipose tissues. Although the function of type XII collagen has not been established as yet, these results suggest that type XII collagen may be associated with adipocyte differentiation and adipose formation in cattle and is a potentially useful marker for adipogenesis.

Differential expression patterns of claudins, tight junction membrane proteins, in mouse nephron segments.

As the first step in understanding the physiologic functions of claudins (tight junction integral membrane proteins) in nephrons, the expression of claudin-1 to -16 in mouse kidneys was examined by Northern blotting. Among these claudins, only claudin-6, -9, -13, and -14 were not detectable. Claudin-5 and -15 were detected only in endothelial cells. Polyclonal antibodies specific for claudin-7 and -12 were not available. Therefore, the distributions of claudin-1, -2, -3, -4, -8, -10, -11, and -16 in nephron segments were examined with immunofluorescence microscopy. For identification of individual segments, antibodies specific for segment markers were used. Immunofluorescence microscopic analyses of serial frozen sections of mouse kidneys with polyclonal antibodies for claudins and segment markers revealed that claudins demonstrated very complicated, segment-specific, expression patterns in nephrons, i.e., claudin-1 and -2 in Bowman's capsule, claudin-2, -10, and -11 in the proximal tubule, claudin-2 in the thin descending limb of Henle, claudin-3, -4, and -8 in the thin ascending limb of Henle, claudin-3, -10, -11, and -16 in the thick ascending limb of Henle, claudin-3 and -8 in the distal tubule, and claudin-3, -4, and -8 in the collecting duct. These segment-specific expression patterns of claudins are discussed, with special reference to the physiologic functions of tight junctions in nephrons.