Identification of contemporary selection signatures using composite log likelihood and their associations with marbling score in Korean cattle

Positive selection not only increases beneficial allele frequency but also causes augmentation of allele frequencies of sequence variants in close proximity. Signals for positive selection were detected by the statistical differences in subsequent allele frequencies. To identify selection signatures in Korean cattle, we applied a composite log‐likelihood (CLL)‐based method, which calculates a composite likelihood of the allelic frequencies observed across sliding windows of five adjacent loci and compares the value with the critical statistic estimated by 50 000 permutations. Data for a total of 11 799 nucleotide polymorphisms were used with 71 Korean cattle and 209 foreign beef cattle. As a result, 147 signals were identified for Korean cattle based on CLL estimates (P < 0.01). The signals might be candidate genetic factors for meat quality by which the Korean cattle have been selected. Further genetic association analysis with 41 intragenic variants in the selection signatures with the greatest CLL for each chromosome revealed that marbling score was associated with five variants. Intensive association studies with all the selection signatures identified in this study are required to exclude signals associated with other phenotypes or signals falsely detected and thus to identify genetic markers for meat quality.     

Genome-wide association study reveals five nucleotide sequence variants for carcass traits in beef cattle

Genetic associations of nucleotide sequence variants with carcass traits in beef cattle were investigated using a genome-wide single nucleotide polymorphism (SNP) assay. Three hundred and thirteen Korean cattle were genotyped with the Illumina BovineSNP50 BeadChip, and 39,129 SNPs from 311 animals were analysed for each carcass phenotype after filtering by quality assurance. Five sequence markers were associated with one of the meat quantity or quality traits; rs109593638 on chromosome 3 with marbling score, rs109821175 on chromosome 11 and rs110862496 on chromosome 13 with backfat thickness (BFT), and rs110228023 on chromosome 6 and rs110201414 on chromosome 16 with eye muscle area (EMA) (P < 1.27 × 10(-6) , Bonferonni P < 0.05). The ss96319521 SNP, located within a gene with functions of muscle development, dishevelled homolog 1 (DVL1), would be a desirable candidate marker. Individuals with genotype CC at this gene appeared to have increased both EMA and carcass weight. Fine-mapping would be required to refine each of the five association signals shown in the current study for future application in marker-assisted selection for genetic improvement of beef quality and quantity.     

Massively parallel sequencing of Chikso (Korean brindle cattle) to discover genome-wide SNPs and InDels

Since the completion of the bovine sequencing projects, a substantial number of genetic variations such as single nucleotide polymorphisms have become available across the cattle genome. Recently, cataloguing such genetic variations has been accelerated using massively parallel sequencing technology. However, most of the recent studies have been concentrated on European Bos taurus cattle breeds, resulting in a severe lack of knowledge for valuable native cattle genetic resources worldwide. Here, we present the first whole-genome sequencing results for an endangered Korean native cattle breed, Chikso, using the Illumina HiSeq 2,000 sequencing platform. The genome of a Chikso bull was sequenced to approximately 25.3-fold coverage with 98.8% of the bovine reference genome sequence (UMD 3.1) covered. In total, 5,874,026 single nucleotide polymorphisms and 551,363 insertion/deletions were identified across all 29 autosomes and the X-chromosome, of which 45% and 75% were previously unknown, respectively. Most of the variations (92.7% of single nucleotide polymorphisms and 92.9% of insertion/deletions) were located in intergenic and intron regions. A total of 16,273 single nucleotide polymorphisms causing missense mutations were detected in 7,111 genes throughout the genome, which could potentially contribute to variation in economically important traits in Chikso. This study provides a valuable resource for further investigations of the genetic mechanisms underlying traits of interest in cattle, and for the development of improved genomics-based breeding tools.     

Genome-Wide Detection of Selective Signature in Chinese Holstein

Selective signatures in whole genome can help us understand the mechanisms of selection and target causal variants for breeding program. In present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions harboring such signals in Chinese Holstein, and then verified the biological significance of these identified regions based on commonly-used bioinformatics analyses. Results showed a total of 125 significant regions in entire genome containing some of important functional genes such as LEP, ABCG2, CSN1S1, CSN3 and TNF based on the Gene Ontology database. Some of these annotated genes involved in the core regions overlapped with those identified in our previous GWAS as well as those involved in a recently constructed candidate gene database for cattle, further indicating these genes under positive selection maybe underlie milk production traits and other important traits in Chinese Holstein. Furthermore, in the enrichment analyses for the second level GO terms and pathways, we observed some significant terms over represented in these identified regions as compared to the entire bovine genome. This indicates that some functional genes associated with milk production traits, as reflected by GO terms, could be clustered in core regions, which provided promising evidence for the exploitability of the core regions identified by EHH tests. Findings in our study could help detect functional candidate genes under positive selection for further genetic and breeding research in Chinese Holstein.     

Genetic Polymorphisms of the CACNA2D1 Gene and Their Association with Carcass and Meat Quality Traits in Cattle

The objective of this study was to identify genetic polymorphisms of the CACNA2D1 gene and to analyze associations between SNPs and carcass and meat quality traits in cattle. Through PCR-RFLP and DNA sequencing methods, a new allelic variant corresponding to the A → G mutation (aspartic to glycine amino acid replacement) of the bovine CACNA2D1 gene was detected. Two alleles and three genotypes (AA, AG, and GG) were defined. Genetic character indicated that the A526745G locus showed moderate polymorphism and was in Hardy–Weinberg equilibrium. Gene-specific SNP marker association analysis showed that the A526745G mutant was significantly associated with carcass weight, dressing percentage, meat percentage, and backfat thickness. The results add new evidence that CACNA2D1 is an important candidate gene for the selection of carcass and meat quality traits in the cattle industry.     

FTO gene variants are associated with growth and carcass traits in cattle

An important aim in animal breeding is the improvement of growth and meat quality traits. Previous studies have demonstrated that genetic variants in the fat mass and obesity associated (FTO) gene have a relatively large effect on human obesity as well as on body composition in rodents and, more recently, in livestock. Here, we examined the effects of the FTO gene variants on growth and carcass traits in the Slovenian population of Simmental (SS) and Brown (SB) cattle. To validate and identify new polymorphisms, we used sequencing, PCR‐RFLP analysis and TaqMan assays in the SS breed and FTO gene variants data from the Illumina BovineSNP50 v1 array for the SB breed. Sequencing of the eight samples of progeny‐tested SS sires detected 108 single nucleotide polymorphisms (SNPs) in the bovine FTO gene. Statistical analyses between growth and carcass traits and 34 FTO polymorphisms revealed significant association of FTO variants with lean meat percentage in both breeds. Additionally, FTO SNPs analyzed in SS cattle were associated with fat percentage, bone weight and live weight at slaughter. The FTO gene can thus be regarded as a candidate gene for the marker‐assisted selection programs in our and possibly other populations of cattle. Future studies in cattle might reveal novel roles for the FTO gene in shaping carcass traits in livestock species as well as body composition control in other mammals.     

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

Background

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F_ST. The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

Results

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F_ST, potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F_ST signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F_ST values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

Conclusions

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-876) contains supplementary material, which is available to authorized users.     

The polymorphisms of LYRM1 gene and their association with body measurement and ultrasound traits of Qinchuan cattle

Body measurement and meat quality traits which play important roles in the assessment of productivity and economy in cattle were influenced by genes and environmental factors. Latest studies showed that LYR motif containing 1 (LYRM1) may be involved in influencing fatness deposition in animals. The objective of this study was to detect bovine LYRM1 gene polymorphism and analyze its association with body measurement and meat quality traits of cattle. Blood samples were taken from a total of 404 Qinchuan cattle aged from 18–24 months. Created restriction site-polymerase chain reaction-restriction fragment length polymorphism (CRS-PCR–RFLP) and DNA sequencing were used to find out LYRM1 single polymorphism nucleotide (SNPs). Sequence analysis of LYRM1 gene revealed two SNPs (g.165 C > A, g.193 A > G) in 3′ untranslated region (3′UTR) of exon 3. And g.165 C > A showed two genotypes namely AC and CC while g.193 A > G showed three genotypes: AA, AG and GG. Analysis results showed that there were significant associations between polymorphism of these two and body measurement and meat quality traits in Qinchuan cattle population. Based on the results obtained from this study, it is inferred that LYRM1 gene may have potential effects on body measurement and meat quality traits in Qinchuan cattle population and could be used for marker-assisted selection.     

Combinations of Susceptibility Genes Are Associated with Higher Risk for Multiple Sclerosis and Imply Disease Course Specificity

Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system that predominantly affects young adults. The genetic contributions to this multifactorial disease were underscored by a genome wide association study (GWAS) conducted by the International Multiple Sclerosis Genetic Consortium in a multinational cohort prompting the discovery of 57 non-MHC MS-associated common genetic variants. Hitherto, few of these newly reported variants have been replicated in larger independent patient cohorts. We genotyped a cohort of 1033 MS patients and 644 healthy controls with a consistent genetic background for the 57 non-MHC variants reported to be associated with MS by the first large GWAS as well as the HLA DRB1*1501 tagging SNP rs3135388. We robustly replicated three of the 57 non-MHC reported MS-associated single nucleotide polymorphisms (SNPs). In addition, our study revealed several genotype-genotype combinations with an evidently higher degree of disease association than the genotypes of the single SNPs. We further correlated well-defined clinical phenotypes, i.e. ataxia, visual impairment due to optic neuritis and paresis with single SNPs and genotype combinations, and identified several associations. The results may open new avenues for clinical implications of the MS associated genetic variants reported from large GWAS.     

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.     

Lack of association between AQP4 polymorphisms and risk of inflammatory demyelinating disease in a Korean population

Multiple sclerosis (MS) and neuromyelitis optica (NMO) are demyelinating autoimmune inflammatory diseases that affect the central nervous system (CNS). Previous genome-wide or candidate gene studies have suggested that genetic variants might be associated with the risk of MS or NMO. Aquaporin 4 (AQP4) is a commonly distributed water channel in astrocytes of the CNS, and its expression is decreased in NMO lesions due to astrocyte cytotoxicity. Previous studies have suggested the associations of AQP4 single nucleotide polymorphisms (SNPs) with MS and/or NMO. However, there have been few replication studies in various ethnic populations. This study, as the first of its kind performed in an Asian population, investigated associations of AQP4 SNPs with the risk of inflammatory demyelinating disease (IDD), including MS and NMO, in a Korean population. A total of seven common AQP4 SNPs were selected based on status of linkage disequilibrium (LD), and then genotyped in 178 IDD cases (79 MS and 99 NMO patients) and 237 normal controls. Statistical analyses showed no significant associations between AQP4 SNPs/haplotypes and development of IDD, including MS and NMO (P > 0.05). Further replications in larger cohorts and other ethnic groups are needed.     

mRNA expression pattern and association study with growth traits of bovine vaspin gene

Visceral adipose tissue-derived serine protease inhibitor (vaspin) is an interesting novel adipocytokine with insulin-sensitizing effects. Some studies have suggested that vaspin could play an important role in the development of obesity and metabolic disorders. However, the tissue expression patterns in cattle and impact of vaspin gene variants on the growth traits has not been determined yet. Herein, we firstly investigated the tissue expression patterns of vaspin gene in new born and adult cattle. The results showed that vaspin was ubiquitously expressed in most tissues and strongly expressed in the heart, skeletal muscle and fat. Then, genetic variants within bovine vaspin gene were screened in 1235 individuals from five Chinese indigenous cattle breeds. Two novel mutations in coding region (NW_001494061: g.1124477 G>A and g.1118561 T>C) of bovine vaspin gene were identified using MspI PCR–RFLP and HhaI ACRS PCR–RFLP detection. Association analysis revealed both two mutations were significantly associated with bodyweight and chest girth at 24 months in cattle (P < 0.05). Therefore, the MspI and HhaI genetic variants of bovine vaspin gene were recommended as DNA markers related to growth traits through marker-assisted selection for genetics and breeding in cattle.     

A Non-Coding Genomic Duplication at the HMX1 Locus Is Associated with Crop Ears in Highland Cattle

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.     

Deciphering signature of selection affecting beef quality traits in Angus cattle

Artificial selection towards a desired phenotype/trait has modified the genomes of livestock dramatically that generated breeds that greatly differ in morphology, production and environmental adaptation traits. Angus cattle are among the famous cattle breeds developed for superior beef quality. This paper aimed at exploring genomic regions under selection in Angus cattle that are associated with meat quality traits and other associated phenotypes. The whole genome of 10 Angus cattle was compared with 11 Hanwoo (A-H) and 9 Jersey (A-J) cattle breeds using a cross-population composite likelihood ratio (XP-CLR) statistical method. The top 1% of the empirical distribution was taken as significant and annotated using UMD3.1. As a result, 255 and 210 genes were revealed under selection from A–H and A–J comparisons, respectively. The WebGestalt gene ontology analysis resulted in sixteen (A–H) and five (A–J) significantly enriched KEGG pathways. Several pathways associated with meat quality traits (insulin signaling, type II diabetes mellitus pathway, focal adhesion pathway, and ECM-receptor interaction), and feeding efficiency (olfactory transduction, tight junction, and metabolic pathways) were enriched. Genes affecting beef quality traits (e.g., FABP3, FTO, DGAT2, ACS, ACAA2, CPE, TNNI1), stature and body size (e.g., PLAG1, LYN, CHCHD7, RPS20), fertility and dystocia (e.g., ESR1, RPS20, PPP2R1A, GHRL, PLAG1), feeding efficiency (e.g., PIK3CD, DNAJC28, DNAJC3, GHRL, PLAG1), coat color (e.g., MC1-R) and genetic disorders (e.g., ITGB6, PLAG1) were found to be under positive selection in Angus cattle. The study identified genes and pathways that are related to meat quality traits and other phenotypes of Angus cattle. The findings in this study, after validation using additional or independent dataset, will provide useful information for the study of Angus cattle in particular and beef cattle in general.     

Polymorphisms of the bovine DKK2 and their associations with body measurement traits and meat quality traits in Qinchuan cattle

The objective of this research were to detect bovine Dickkopf 2 (DKK2) gene polymorphism and analyze their associations with body measurement traits (BMT) and meat quality traits (MQT) of animals. Blood samples were taken from a total of 541 Qinchuan cattle aged from 18 to 24 months. Polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) was employed to find out DKK2 single-polymorphism nucleotide (SNPs) and to explore their possible association with BMT and MQT. Sequence analysis of DKK2 gene revealed 2 SNPs (C29 T and A169C) in 5′ untranslated region (5′UTR) of exon 1.C29T and A164T SNPs are both synonymous mutation, which showed 2 genotypes namely (CC, CT) and (AA and AC), respectively. Association analysis of polymorphism with body measurement and meat quality traits at the two locus showed that there were significant effects on CT, BL, RL, PBW, BFT, LMA, and IFC. These results suggest that the DKK2 gene might have potential effects on BMT and MQT in Qinchuan cattle population and could be used for marker-assisted selection.     

A genome‐wide scan for selection signatures in Nellore cattle

Brazilian Nellore cattle (Bos indicus) have been selected for growth traits for over more than four decades. In recent years, reproductive and meat quality traits have become more important because of increasing consumption, exports and consumer demand. The identification of genome regions altered by artificial selection can potentially permit a better understanding of the biology of specific phenotypes that are useful for the development of tools designed to increase selection efficiency. Therefore, the aims of this study were to detect evidence of recent selection signatures in Nellore cattle using extended haplotype homozygosity methodology and BovineHD marker genotypes (>777 000 single nucleotide polymorphisms) as well as to identify corresponding genes underlying these signals. Thirty‐one significant regions (P < 0.0001) of possible recent selection signatures were detected, and 19 of these overlapped quantitative trait loci related to reproductive traits, growth, feed efficiency, meat quality, fatty acid profiles and immunity. In addition, 545 genes were identified in regions harboring selection signatures. Within this group, 58 genes were associated with growth, muscle and adipose tissue metabolism, reproductive traits or the immune system. Using relative extended haplotype homozygosity to analyze high‐density single nucleotide polymorphism marker data allowed for the identification of regions potentially under artificial selection pressure in the Nellore genome, which might be used to better understand autozygosity and the effects of selection on the Nellore genome.     

Copy number variation detection in Chinese indigenous cattle by whole genome sequencing

Copy number variation (CNV) refers to a kind of structural variation, having functional and evolutionary effects on phenotypes. Thus far, further elucidation of the CNVs in different Chinese indigenous cattle breeds by whole genome sequencing have yet not been done. In this study, a comprehensive genomic analysis was performed on 75 cattle individuals including six Chinese indigenous cattle breeds and two non-native specialized beef cattle breeds. Based on the 11,486 CNVRs discovered, population analysis was performed, showed that all the cattle breeds clustered in to three clades, consistent with their lineages Bos taurus, Bos taurus × Bos indicus and Bos indicus. Importantly, a set of CNVRs related genes were found to be associated with the traits of interest, which include meat production or quality (CAST, ACTC1, etc.), adaption (BLA-DQB, EGLN2, etc.) and coat color (KIT, MITF, etc.). These results provide valuable full genome variation resources for Chinese bovine genome research and would be helpful for cattle breeding and selection programs in the future.     

Novel polymorphisms of SIX4 gene and their association with body measurement traits in Qinchuan cattle

Sine oculis homeobox homolog 4 (SIX4) gene belongs to the sine oculis/SIX gene family, which includes six members in vertebrates. SIX4 gene plays a crucial role in skeletal myogenesis, and its genetic variations or deficiency may cause hypopituitarism, suggesting that SIX4 gene is a potential candidate gene affecting body measurement traits (BMTs) in animals. Herein, the objectives of this study were to identify genetic polymorphisms of bovine SIX4 gene and to analyze potential association between single nucleotide polymorphisms (SNPs) and body measurement traits in Qinchuan cattle. In the present study, we investigated polymorphisms of SIX4 gene in 426 Qinchuan cattle using DNA sequencing and polymerase chain reaction–restriction fragment length polymorphisms. Three novel SNPs were identified within bovine SIX4 gene. Associations between body measurement traits and SIX4 gene polymorphisms were investigated, and significant statistical associations were found between polymorphisms of these three SNPs and body measurement traits (P < 0.05). Hence, based on results obtained from this study, we conjectured that SIX4 gene may have potential effects on body measurement traits in Qinchuan cattle population and could be used for marker-assisted selection.     

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.     

Rice-Map: a new-generation rice genome browser

Background

Because the Japanese native cattle Kuchinoshima-Ushi have been isolated in a small island and their lineage has been intensely protected, it has been assumed to date that numerous and valuable genomic variations are conserved in this cattle breed.

Results

In this study, we evaluated genetic features of this breed, including single nucleotide polymorphism (SNP) information, by whole-genome sequencing using a Genome Analyzer II. A total of 64.2 Gb of sequence was generated, of which 86% of the obtained reads were successfully mapped to the reference sequence (Btau 4.0) with BWA. On an average, 93% of the genome was covered by the reads and the number of mapped reads corresponded to 15.8-fold coverage across the covered region. From these data, we identified 6.3 million SNPs, of which more than 5.5 million (87%) were found to be new. Out of the SNPs annotated in the bovine sequence assembly, 20,432 were found in protein-coding regions containing 11,713 nonsynonymous SNPs in 4,643 genes. Furthermore, phylogenetic analysis using sequence data from 10 genes (more than 10 kbp) showed that Kuchinoshima-Ushi is clearly distinct from European domestic breeds of cattle.

Conclusions

These results provide a framework for further genetic studies in the Kuchinoshima-Ushi population and research on functions of SNP-containing genes, which would aid in understanding the molecular basis underlying phenotypic variation of economically important traits in cattle and in improving intrinsic defects in domestic cattle breeds.