Measuring conservation of contiguous sets of autosomal markers on bovine and porcine genomes in relation to the map of the human genome.

Based on published information, we have identified 991 genes and gene-family clusters for cattle and 764 for pigs that have orthologues in the human genome. The relative linear locations of these genes on human sequence maps were used as "rulers" to annotate bovine and porcine genomes based on a CSAM (contiguous sets of autosomal markers) approach. A CSAM is an uninterrupted set of markers in one genome (primary genome; the human genome in this study) that is syntenic in the other genome (secondary genome; the bovine and porcine genomes in this study). The analysis revealed 81 conserved syntenies and 161 CSAMs between human and bovine autosomes and 50 conserved syntenies and 95 CSAMs between human and porcine autosomes. Using the human sequence map as a reference, these 991 and 764 markers could correlate 72 and 74% of the human genome with the bovine and porcine genomes, respectively. Based on the number of contiguous markers in each CSAM, we classified these CSAMs into five size groups as follows: singletons (one marker only), small (2-4 markers), medium (5-10 markers), large (11-20 markers), and very large (> 20 markers). Several bovine and porcine chromosomes appear to be represented as di-CSAM repeats in a tandem or dispersed way on human chromosomes. The number of potential CSAMs for which no markers are currently available were estimated to be 63 between human and bovine genomes and 18 between human and porcine genomes. These results provide basic guidelines for further gene and QTL mapping of the bovine and porcine genomes, as well as insight into the evolution of mammalian genomes.     

The Mammalian 2′-5′ Oligoadenylate Synthetase Gene Family: Evidence for Concerted Evolution of Paralogous Oas1 Genes in Rodentia and Artiodactyla

Multiple 2′-5′ oligoadenylate (2-5A) synthetases are important components of innate immunity in mammals. Gene families encoding these proteins have previously been studied mainly in humans and mice. To reconstruct the evolution of this gene family in mammals, a search for additional 2-5A synthetase genes was performed in rat, cattle, pig, and dog. Twelve 2′-5′ oligoadenylate synthetase (Oas) genes were identified in the rat genome, including eight Oas1 genes, two Oas1 pseudogenes, single copies of Oas2 and Oas3, and two Oas-like genes, Oasl1 and Oasl2. Four OAS genes were detected in the pig genome and five OAS genes were found in both the cattle and dog genomes. An OAS3 gene was not found in either the cattle or the pig genome. While two tandemly duplicated OAS-like (OASL) genes were identified in the dog genome, only a single OASL orthologue was found in both the cattle and the pig genomes. The bovine and porcine OASL genes contain premature stop codons and encode truncated proteins, which lack the typical C-terminal double ubiquitin domains. The cDNA sequences of the rat, cattle, pig, and dog OAS genes were amplified, sequenced and compared with each other and with those in the human, mouse, horse, and chicken genomes. Evidence of concerted evolution of paralogous 2′-5′ oligoadenylate synthetase 1 genes was obtained in rodents (Rodentia) and even-toed ungulates (Artiodactyla). Calculations using the nonparametric Kolmogorov-Smirnov test suggested that the homogenization of paralogous OAS1 sequences was due to gene conversion rather than stabilizing selection. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. Martin Kreitman     

Multi-species comparative mapping in silico using the COMPASS strategy

MOTIVATION: The completion of human and mouse genome sequences provides a valuable resource for decoding other mammalian genomes. The comparative mapping by annotation and sequence similarity (COMPASS) strategy takes advantage of the resource and has been used in several genome-mapping projects. It uses existing comparative genome maps based on conserved regions to predict map locations of a sequence. An automated multiple-species COMPASS tool can facilitate in the genome sequencing effort and comparative genomics study of other mammalian species. RESULTS: The prerequisite of COMPASS is a comparative map table between the reference genome and the predicting genome. We have built and collected comparative maps among five species including human, cattle, pig, mouse and rat. Cattle-human and pig-human comparative maps were built based on the positions of orthologous markers and the conserved synteny groups between human and cattle and human and pig genomes, respectively. Mouse-human and rat-human comparative maps were based on the conserved sequence segments between the two genomes. With a match to human genome sequences, the approximate location of a query sequence can be predicted in cattle, pig, mouse and rat genomes based on the position of the match relatively to the orthologous markers or the conserved segments. AVAILABILITY: The COMPASS-tool and databases are available at http://titan.biotec.uiuc.edu/COMPASS/     

Comparative chromosome painting defines the high rate of karyotype changes between pigs and bovids

Human and sheep chromosome-specific probes were used to construct comparative painting maps between the pig (Suiformes), cattle and sheep (Bovidae), and humans. Various yet unknown translocations were observed that would assist in a more complete reconstruction of homology maps of these species. The number of homologous segments that can be identified with sheep probes in the pig karyotype exceeds that described previously by chromosome painting between two non-primate mammals belonging to the same order. Sheep probes painted 62 segments on pig autosomes and delineated not only translocations, but also 9 inversions. All inversions were paracentric and indicate that these rearrangements may be characteristic for chromosomal changes in suiforms. Hybridizations of all sheep painting probes to cattle chromosomes confirmed the chromosome conservation in bovids. In addition, we observed a small translocation that was previously postulated from linkage mapping data, but was not yet described by physical mapping. The chromosome painting data are complemented with a map of available comparative gene mapping data between pig and sheep genomes. A detailed table listing the comparative gene mapping data between pig and cattle genomes is provided. The reanalysis of the pig karyotype with a new generation of human paint probes provides an update of the human/pig comparative genome map and demonstrates two new chromosome homologies. Seven conserved segments not yet identified by chromosome painting are also reported. Received: 2 October 2000 / Accepted: 15 January 2001     

Mapping, Phylogenetic and Expression Analysis of the RNase (RNaseA) Locus in Cattle

The mammalian secreted ribonucleases (RNases) comprise a large family of structurally related proteins displaying considerable sequence variation, and have been used in evolutionary studies. RNase 1 (RNase A) has been assumed to play a role in digestion, while other members have been suggested to contribute to host defence. Using the recently assembled bovine genome sequence, we characterised the complete repertoire of genes present in the RNaseA family locus in cattle, and compared this with the equivalent locus in the human and mouse genomes. Several additions and corrections to the earlier analysis of the RNase locus in the mouse genome are presented. The bovine locus encodes 19 RNases, of which only six have unambiguous equivalent genes in the other two species. Chromosomal mapping and phylogenetic analysis indicate that a number of distinct gene duplication events have occurred in the cattle lineage since divergence from the human and mouse lineages. Substitution analysis suggests that some of these duplicated genes are under evolutionary pressure for purifying selection and may therefore be important to the physiology of cattle. Expression analysis revealed that individual RNases have a wide pattern of expression, including diverse mucosal epithelia and immune-related cells and tissues. These data clarify the full repertoire of bovine RNases and their relationships to those in humans and mice. They also suggest that RNase gene duplication within the bovine lineage accompanied by altered tissue-specific expression has contributed a survival advantage.     

Comparative human-mouse-rat sequence analysis of the ICAM gene cluster on HSA 19p13.2 and a 185-kb porcine region from SSC 2q

The human intercellular adhesion molecule gene (ICAM) cluster is located in a GC-rich and gene-rich region on HSA 19p13.2. We determined the complete DNA sequence of a 185-kb porcine bacterial artificial chromosome (BAC) clone containing parts of the ICAM gene cluster. We used the porcine sequence for a detailed comparative analysis between human, pig, mouse and rat. The 185 kb of porcine sequence covered 220 kb of homologous sequence in the human genome, which adds to the growing evidence that the porcine genome is somewhat smaller than the human genome. The genomic sequences of the four species showed a high level of conserved synteny and no rearrangements in gene order were observed. During evolution, the ICAM3 gene was inactivated by mutation in the mouse and rat genome, whereas it is still present in the human and pig genome. The loss of Icam3 in rodent genomes might be relevant for rodent-specific properties of the T-cell-mediated immune response. All the other investigated genes are conserved across all four investigated sequences.     

Phylogenetic analysis and identification of pseudogenes reveal a progressive loss of zona pellucida genes during evolution of vertebrates

Vertebrate eggs are surrounded by an extracellular matrix with similar functions and conserved individual components: the zona pellucida (ZP) glycoproteins. In mammals, chickens, frogs, and some fish species, we established an updated list of the ZP genes, studied the relationships within the ZP gene family using phylogenetic analysis, and identified ZP pseudogenes. Our study confirmed the classification of ZP genes in six subfamilies: ZPA/ZP2, ZPB/ZP4, ZPC/ZP3, ZP1, ZPAX, and ZPD. The identification of a Zpb pseudogene in the mouse genome, Zp1 pseudogenes in the dog and bovine genomes, and Zpax pseudogenes in the human, chimpanzee, macaque, and bovine genomes showed that the evolution of ZP genes mainly occurs by death of genes. Our study revealed that the extracellular matrix surrounding vertebrate eggs contains three to at least six ZP glycoproteins. Mammals can be classified in three categories. In the mouse, the ZP is composed of three ZP proteins (ZPA/ZP2, ZPC/ZP3, and ZP1). In dog, cattle and, putatively, pig, cat, and rabbit, the zona is composed of three ZP proteins (ZPA/ZP2, ZPB/ZP4, and ZPC/ZP3). In human, chimpanzee, macaque, and rat, the ZP is composed of four ZP proteins (ZPA/ZP2, ZPB/ZP4, ZPC/ZP3, and ZP1). Our review provides new directions to investigate the molecular basis of sperm-egg recognition, a mechanism which is not yet elucidated.     

A comprehensive radiation hybrid map of bovine Chromosome 26 (BTA26): comparative chromosomal organization between HSA10q and BTA26 and BTA28

In this study, we present a comprehensive 3000-rad radiation hybrid (RH) map of bovine Chromosome (Chr) 26 (BTA26) with 80 markers including 50 genes or ESTs: 44 have an ortholog mapping to human Chr 10 (HSA10) and 29 to mouse Chr (MMU) 7, 10, and 19. Moreover, 12 other HSA10 genes were integrated in a newly developed RH map of BTA28 (seven represent new assignments). The available draft of the mouse genome allowed us to present a detailed picture of the distribution of conserved synteny segments among the three species (human, cattle, and mouse) and to propose a simple model of the comparative chromosomal organization between the long arm of HSA10 and BTA26 and 28. Finally, the INRA bovine BAC library was screened for most of the BTA26 markers considered in this study to provide anchors for the bovine physical map.     

Integration of hybridization-based markers (overgos) into physical maps for comparative and evolutionary explorations in the genus Oryza and in Sorghum

Background

On porcine chromosome 7, the region surrounding the Major Histocompatibility Complex (MHC) contains several Quantitative Trait Loci (QTL) influencing many traits including growth, back fat thickness and carcass composition. Previous studies highlighted that a fragment of ~3.7 Mb is located within the Swine Leucocyte Antigen (SLA) complex. Internal rearrangements of this fragment were suggested, and partial contigs had been built, but further characterization of this region and identification of all human chromosomal fragments orthologous to this porcine fragment had to be carried out.

Results

A whole physical map of the region was constructed by integrating Radiation Hybrid (RH) mapping, BAC fingerprinting data of the INRA BAC library and anchoring BAC end sequences on the human genome. 17 genes and 2 reference microsatellites were ordered on the high resolution IMNpRH2_12000rad Radiation Hybrid panel. A 1000:1 framework map covering 550 cR₁₂₀₀₀ was established and a complete contig of the region was developed. New micro rearrangements were highlighted between the porcine and human genomes. A bovine RH map was also developed in this region by mapping 16 genes. Comparison of the organization of this region in pig, cattle, human, mouse, dog and chicken genomes revealed that 1) the translocation of the fragment described previously is observed only on the bovine and porcine genomes and 2) the new internal micro rearrangements are specific of the porcine genome.

Conclusion

We estimate that the region contains several rearrangements and covers 5.2 Mb of the porcine genome. The study of this complete BAC contig showed that human chromosomal fragments homologs of this heavily rearranged QTL region are all located in the region of HSA6 that surrounds the centromere. This work allows us to define a list of all candidate genes that could explain these QTL effects.     

A high utility integrated map of the pig genome

Background

The domestic pig is being increasingly exploited as a system for modeling human disease. It also has substantial economic importance for meat-based protein production. Physical clone maps have underpinned large-scale genomic sequencing and enabled focused cloning efforts for many genomes. Comparative genetic maps indicate that there is more structural similarity between pig and human than, for example, mouse and human, and we have used this close relationship between human and pig as a way of facilitating map construction.

Results

Here we report the construction of the most highly continuous bacterial artificial chromosome (BAC) map of any mammalian genome, for the pig (Sus scrofa domestica) genome. The map provides a template for the generation and assembly of high-quality anchored sequence across the genome. The physical map integrates previous landmark maps with restriction fingerprints and BAC end sequences from over 260,000 BACs derived from 4 BAC libraries and takes advantage of alignments to the human genome to improve the continuity and local ordering of the clone contigs. We estimate that over 98% of the euchromatin of the 18 pig autosomes and the X chromosome along with localized coverage on Y is represented in 172 contigs, with chromosome 13 (218 Mb) represented by a single contig. The map is accessible through pre-Ensembl, where links to marker and sequence data can be found.

Conclusion

The map will enable immediate electronic positional cloning of genes, benefiting the pig research community and further facilitating use of the pig as an alternative animal model for human disease. The clone map and BAC end sequence data can also help to support the assembly of maps and genome sequences of other artiodactyls.     

In silico identification and Bayesian phylogenetic analysis of multiple new mammalian kallikrein gene families

Kallikrein gene families have been identified previously in genomes of the human, the mouse, and the rat, and individual kallikrein-like genes have been found in many more species. This study presents the in silico identification of kallikrein gene families in the recently sequenced genomes of four additional mammalian species, the chimpanzee, the dog, the pig, and the opossum. Phylogenies were constructed with gene sequences from all seven mammalian families, using Bayesian analysis, which clarified the evolutionary relationships between these genes. Individual gene sequences, as well as concatenated constructs of multiple sequences, were used. Fifteen kallikrein genes were located in the chimpanzee (Pan troglodytes) genome, while only 14 were identified in the canine (Canis familiaris) genome as no orthologue to human KLK3 was found. Thirteen genes were identified from the pig (Sus scrofa) genome, which lacked homologues to KLK2 and KLK3, and 11 genes, orthologous to human KLK5 through KLK15, were found in the opossum (Monodelphis domestica) genome. No kallikrein genes were identified from the available genome sequences of the chicken (Gallus gallus) or African clawed frog (Xenopus tropicalis). Within the family of kallikreins several subfamilies were suggested by phylogenetic analysis. One consisted of KLK4, KLK5, and KLK14; another of KLK9, KLK11, and KLK15; a third of KLK10 and KLK12; a fourth of KLK6 and KLK13; and finally one of KLK8 and the classical kallikreins (KLK1, KLK2, and KLK3).     

Patterns of transitional mutation biases within and among mammalian genomes

Significant transition/transversion mutation bias is a well-appreciated aspect of mammalian nuclear genomes; however, patterns of bias among genes within a genome and among species remain largely uncharacterized. Understanding these patterns is important for understanding similarities and differences in mutational patterns among genomes and genomic regions. Therefore, we have conducted an analysis of 7,587 pairs of sequences of 4,347 mammalian protein-coding genes from seven species (human, mouse, rat, cow, sheep, pig, and macaque) and from the introns of 51 gene pairs and multiple intergenic regions (37 kbp, 52 kbp and 65 kbp) from the human, chimpanzee, and baboon genomes. Our analyses show that genes and regions with widely varying base composition exhibit uniformity of transition mutation rate both within and among mammalian lineages, as long as the transitional mutations caused by CpG hypermutability are excluded. The estimates show no relationship to potential intrachromosomal or interchromosomal effects. This uniformity points to similarity in point mutation processes in genomic regions with substantially different GC-content biases.     

Features and trend of loss of promoter-associated CpG islands in the human and mouse genomes

CpG islands (CGIs) are often considered as gene markers, but the number of CGIs varies among mammalian genomes that have similar numbers of genes. In this study, we investigated the distribution of CGIs in the promoter regions of 3,197 human-mouse orthologous gene pairs and found that the mouse genome has notably fewer CGIs in the promoter regions and less pronounced CGI characteristics than does the human genome. We further inferred CGI's ancestral state using the dog genome as a reference and examined the nucleotide substitution pattern and the mutational direction in the conserved regions of human and mouse CGIs. The results reveal many losses of CGIs in both genomes but the loss rate in the mouse lineage is two to four times the rate in the human lineage. We found an intriguing feature of CGI loss, namely that the loss of a CGI usually starts from erosion at the both edges and gradually moves towards the center. We found functional bias in the genes that have lost promoter-associated CGIs in the human or mouse lineage. Finally, our analysis indicates that the association of CGIs with housekeeping genes is not as strong as previously estimated. Our study provides a detailed view of the evolution of promoter-associated CGIs in the human and mouse genomes and our findings are helpful for understanding the evolution of mammalian genomes and the role of CGIs in gene function.     

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.     

A set of 99 cattle microsatellites: characterization,synteny mapping,and polymorphism

Cattle microsatellite clones (136) were isolated from cosmid (10) and plasmid (126) libraries and sequenced. The dinucleotide repeats were studied in each of these sequences and compared with dinucleotide repeats found in other vertebrate species where information was available. The distribution in cattle was similar to that described for other mammals, such as rat, mouse, pig, or human. A major difference resides in the number of sequences present in the bovine genome, which seemed at best one-third as large as in other species. Oligonucleotide primers (117 pairs) were synthesized, and a PCR product of expected size was obtained for 88 microsatellite sequences (75%). Synteny or chromosome assignment was searched for each locus with PCR amplification on a panel of 36 hamster/bovine somatic cell hybrids. Of our bovine microsatellites, eighty-six could be assigned to synteny groups of chromosomes. In addition, 10 other microsatellites—HEL 5, 6, 9, 11, 12, 13 (Kaukinen and Varvio 1993), HEL 4, 7, 14, 15—as well as the microsatellite found in the -casein gene (Fries et al. 1990) were mapped on the hybrids. Microsatellite polymorphism was checked on at leat 30 unrelated animals of different breeds. Almost all the autosomal and X Chr microsatellites displayed polymorphism, with the number of alleles varying between two and 44. We assume that these microsatellites could be very helpful in the construction of a primary public linkage map of the bovine genome, with an aim of finding markers for Economic Trait Loci (ETL) in cattle.     

Arachnomelia in Brown Swiss cattle maps to chromosome 5

Arachnomelia in Brown Swiss cattle is a monogenic autosomal recessive inherited congenital disorder of the skeletal system giving affected calves a spidery look (OMIA ID 000059). Over a period of 20 years 15 cases were sampled in the Swiss and Italian Brown cattle population. Pedigree data revealed that all affected individuals trace back to a single acknowledged carrier founder sire. A genome scan using 240 microsatellites spanning the 29 bovine autosomes showed homozygosity at three adjacent microsatellite markers on bovine Chr 5 in all cases. Linkage analysis confirmed the localization of the arachnomelia mutation in the region of the marker ETH10. Fine-mapping and haplotype analysis using a total of 34 markers in this region refined the critical region of the arachnomelia locus to a 7.19-Mb interval on bovine Chr 5. The disease-associated IBD haplotype was shared by 36 proven carrier animals and allows marker-assisted selection. As the corresponding human and mouse chromosome segments do not contain any clear functional candidate genes for this disorder, the mutation causing arachnomelia in the Brown Swiss cattle might help to identify an unknown gene in bone development. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Identification, evolution, and regulation of expression of Guinea pig trappin with an unusually long transglutaminase substrate domain

Trappins are found in human, bovine, hippopotamus, and members of the pig family, but not in rat and mouse. To clarify the evolution of the trappin genes and the functional significance of their products, we isolated the trappin gene in guinea pig, a species belonging to a rodent family distinct from rat and mouse. Guinea pig trappin was confirmed to encode the same domain structure as trappin, consisting of a signal sequence, an extra large transglutaminase substrate domain, and a whey acidic protein motif. Northern blot analysis and in situ hybridization histochemistry as well as immunohistochemistry demonstrated that guinea pig trappin is expressed solely in the secretory epithelium of the seminal vesicle and that its expression is androgen-dependent. We confirmed that guinea pig trappin is cross-linked by prostate transglutaminase and that the whey acidic protein motif derived from guinea pig trappin has an inhibitory activity against leukocyte elastase. Genome sequence analysis showed that guinea pig trappin belongs to the family of REST (rapidly evolving seminal vesicle transcribed) genes.     

Comparative analysis of genome maintenance genes in naked mole rat,mouse, and human

Genome maintenance (GM) is an essential defense system against aging and cancer, as both are characterized by increased genome instability. Here, we compared the copy number variation and mutation rate of 518 GM‐associated genes in the naked mole rat (NMR), mouse, and human genomes. GM genes appeared to be strongly conserved, with copy number variation in only four genes. Interestingly, we found NMR to have a higher copy number of CEBPG, a regulator of DNA repair, and TINF2, a protector of telomere integrity. NMR, as well as human, was also found to have a lower rate of germline nucleotide substitution than the mouse. Together, the data suggest that the long‐lived NMR, as well as human, has more robust GM than mouse and identifies new targets for the analysis of the exceptional longevity of the NMR.     

Relationship between the major protein kinase C substrates acidic 80-kDa protein-kinase-C substrate (80K) and myristoylated alanine-rich C-kinase substrate (MARCKS). Members of a gene family or equivalent genes in different species.

Two major protein-kinase-C (PKC) substrates have been described in the literature; an 87-kDa bovine and human PKC substrate, called MARCKS, and an acidic 80-kDa PKC substrate, isolated from rat brain and Swiss 3T3 cells, termed 80K. Since there is only 66-74% sequence similarity between MARCKS and 80K, we have further investigated their relationship in this study. Southern-blot experiments with gene-specific probes demonstrated the presence of the 80K, but not MARCKS, gene in the mouse genome. Furthermore, polymerase-chain-reaction (PCR) analyses using three pairs of primers that specifically recognise either 80K, MARCKS or conserved sequences of both genes, revealed the presence of only the 80K gene in the mouse and rat genomes and only the MARCKS gene in the bovine and human genomes with mRNA expression in the corresponding brain tissues. Northern-blot analysis of a variety of tissues indicated that both 80K and MARCKS have similar patterns of expression. Most components of signal-transduction pathways are present in multiple molecular isoforms as members of a gene family. In contrast, the findings presented in this study indicate that rodent 80K and bovine and human MARCKS are not distinct members of a gene family, but represent the equivalent substrates in different species.     

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.