Mitochondrial DNA diversity patterns in Pakistani buffalo

The Indian subcontinent is considered to be the likely centre of river buffalo domestication, based on population dynamics, archaeological evidence and genetic diversity. Recent studies on mitochondrial DNA diversity have drawn useful conclusions about the domestication history of Bubalus bubalis. The conclusions of these studies are, however, incomplete, unless samples can also be analysed from Pakistan, which contains the second largest buffalo population of the world. Here, we report the results of the first study on mitochondrial D-loop sequence diversity in five breeds of Pakistani buffalo. Analysis of sequence variations in 503-bp of the D-loop region of 123 animals revealed 52 haplotypes, including 40 singletons. Multidimensional display of breed pairwise F(ST) values revealed no strong clustering of breeds. Bayesian, maximum parsimony, neighbour joining and UPGMA trees revealed a topology consistent with domestication as well as subsequent introgression of multiple maternal lineages from the wild stocks. Reduced median network analysis provided evidence of population expansion from more than one set of haplotypes. The study also confirmed that Pakistani buffalo are of the river type. The observed mitochondrial D-loop sequence diversity suggests that Pakistani areas bordering India might have contributed to the initiation of domestication of the present-day river buffalo.     

DNA fingerprinting of Pakistani buffalo breeds (Nili-Ravi, Kundi) using microsatellite and cytochrome b gene markers

Cytochrome b gene markers have been proved as an efficient and powerful tool for breed characterization and species identification of buffaloes. This study represents the substantial analysis of mitochondrial DNA variation in Pakistani buffalo breeds and provides information about their genetic diversity. In this study partial amplification of cytochrome b gene of 1,061 bp was done and sequencing results showed ten haplotypes. Comparing all fifty samples from two buffalo breeds of Pakistan, fifteen polymorphic sites were observed out of which, twelve codons 42, 71, 118, 120, 199, 235, 269, 297, 318, 327, 350, 355 of mitochondrial cytochrome b gene are monomorphic which translate same amino acids as in the reference protein sequence due to silent mutation while different in DNA sequence. Similarly three codons 163, 246, 337 of mitochondrial cytochrome b are polymorphic and different from the reference sequence with respect to DNA as well as protein sequence. For the further confirmation a panel of nine microsatellite markers was used with high polymorphism information content (PIC). The frequency distribution of these alleles varies from three to eight allele at locus CSSM66 and ILST029 respectively. The results obtained from this study may contribute to the establishment of routine genotyping service of buffalo breeds for buffalo farmers for animal forensic application in case of any dispute. Additionally this study may help for breed characterization and phylogeny of aforementioned breeds of buffalo.     

Genetic diversity analysis of an indigenous Chinese buffalo breed and hybrids based on microsatellite data

Chinese native buffaloes have faced the threat of extinction, along with an increase in crossbreeding with domesticated river buffaloes; consequently, conservation of local buffalo genetic resources has become a priority. A Chinese native breed, Jianghan, is often crossed intentionally and unintentionally with imported breeds from India and Pakistan, Murrah, and Nili-Ravi. A total of 128 buffaloes of the breeds Jianghan, Murrah, and Nili-Ravi and their presumed hybrid offspring were genotyped for 10 microsatellite markers. Heterozygosity and Wright's F-statistics were calculated to determine the genetic variation in those populations. The observed average heterozygosities ranged from 0.836 (Murrah) to 0.986 (Jianghan), higher than the expected heterozygosities and all the inbreeding values within the populations were negative. The genetic distances between the presumed hybrid buffaloes and the two imported river type dairy buffalo breeds (Murrah and Nili-Ravi) were lower than with the native Jianghan, indicating strong contributions of the imported breeds to this presumed hybrid buffalo population. This information will be useful for the development of rational breeding for the dairy buffalo industry and for conservation strategies for the Jianghan buffalo.     

Mitochondrial DNA Cleavage Patterns Distinguish Independent Origin of Chinese Domestic Geese and Western Domestic Geese

It has generally been assumed, based on morphology, that Chinese domestic goose breeds were derived from the swan goose (Anser cygnoides) and that European and American breeds were derived from the graylag goose (Anser anser). To test the validity of this assumption, we investigated the mtDNA cleavage patterns of 16 Chinese breeds and 2 European breeds as well as hybrids produced between a Chinese breed and a European breed. After 224 mtDNAs, isolated from the Chinese and European breeds, were digested by 19 restriction endonucleases, variations of the cleavage patterns were observed for four enzymes (EcoRV, HaeII, HincII, and KpnI). All Chinese breeds and their maternal hybrids except the Yili breed showed an identical haplotype, named haplotype I or the Chinese haplotype; the European breeds and the Yili breed showed another haplotype, named haplotype II or the western haplotype. None of the haplotype found in the Chinese type was detectable in the western type and vice versa. The two haplotypes were found to differ from each other at 8.0% of the sites surveyed and with a 0.72% sequence divergence. Using 2% substitution per million years calibrated from the genera Anser and Branta, the two domestic geese haplotypes were estimated to have diverged approximately 360,000 years ago, well outside the 3000-6000 years in domestic history. Our findings provide the first molecular genetic evidence to support the dual origin assumption of domestic geese in the world. Meanwhile, the four mtDNA restriction fragment length polymorphisms can be used as maternal genetic markers to distinguish the two types of domestic geese.     

Identification of X- and Y-chromosome bearing buffalo (Bubalus bubalis) sperm

The aim of this study was to identify the difference in DNA content characterizing the X- and Y-chromosome bearing sperm of buffalo. Sperm from six Murrah buffaloes and six Nili-Ravi buffaloes were collected and stained with Hoechst 33342 followed by flow cytometry analysis of the DNA content. Two symmetrical, separate but overlapping peaks presumed to be X- and Y-chromosome bearing sperm were detected. The difference in fluorescence intensity, which related to the DNA content, between the X- and Y-sperm was 3.59+/-0.11% for Murrah buffalo and 3.55+/-0.14% for Nili-Ravi buffalo, respectively. Significant differences were observed among males within each breed, but there were no differences between the averages of the two breeds. The results indicate that flow cytometric sorting of X- and Y-sperm of buffalo is feasible.     

Prion protein gene polymorphisms in four goat breeds of Pakistan

Four different goat breeds (Pak-Angora, Dera Din Panah, Naachi and Teddy) of Pakistan were selected to investigate polymorphism in the prion protein gene (PrP gene) responsible for scrapie disease resistance in goats. Initially, genotyping of 187 animals of these four breeds by restriction fragment length polymorphism (RFLP) was done to see the genotype for codon 136 and 154. All the animals were monomorphic with a genotype of AARR except one animal of Teddy breed having the genotype of AARH. Sequencing of PrP gene of twenty animals representing these four goat breeds revealed two genotypes PPSSSS and PPSSPS with haplotypes PSS and PSP of PrP gene at the codon numbers 42, 138, and 240. All four breeds showed both wild type monomorphic sequence and mutant polymorphic sequences of these codons. The mutants of 42 and 138 codons translate the same amino acids as with the wild type sequences, while the mutant of codon 240 is responsible for a different amino acid translation i.e., serine to proline. In short, this study provides preliminary information about alleles and genotypes of PrP gene in four goat breeds of Pakistan.     

Mitochondrial haplotypes reveal a strong genetic structure for three Indian sheep breeds

This survey represents the first characterization of mitochondrial DNA diversity within three breeds of Indian sheep (two strains of the Deccani breed, as well as the Bannur and Garole breeds) from different geographic regions and with divergent phenotypic characteristics. A 1061-bp fragment of the mitochondrial genome spanning the control region, a portion of the 12S rRNA gene and the complete phenyl tRNA gene, was sequenced from 73 animals and compared with the corresponding published sequence from European and Asian breeds and the European Mouflon (Ovis musimon). Analysis of all 156 sequences revealed 73 haplotypes, 52 of which belonged to the Indian breeds. The three Indian breeds had no haplotypes in common, but one Indian haplotype was shared with European and other Asian breeds. The highest nucleotide and haplotype diversity was observed in the Bannur breed (0.00355 and 0.981 respectively), while the minimum was in the Sangamneri strain of the Deccani breed (0.00167 and 0.882 respectively). All 52 Indian haplotypes belonged to mitochondrial lineage A. Therefore, these Indian sheep are distinct from other Asian and European breeds studied so far. The relationships among the haplotypes showed strong breed structure and almost no introgression among these Indian breeds, consistent with Indian sheep husbandry, which discourages genetic exchange between breeds. These results have implications for the conservation of India's ovine biodiversity and suggest a common origin for the breeds investigated.     

Novel SNPs of the mtDNA ND5 Gene and Their Associations with Several Growth Traits in the Nanyang Cattle Breed

Polymorphism of the mtDNA ND5 gene was detected by PCR-SSCP and DNA sequencing methods in 714 individuals from six Chinese cattle breeds. The results showed that at the P2 locus, two kinds of haplotypes, named A and B, with three SNPs (T12900C, A12923T, C12924T) were detected. The frequencies of haplotype A in six breeds were 0.041–0.167. Polymorphism of the ND5 gene was shown to be associated with growth in the Nanyang breed. Individuals with haplotype B had greater hucklebone width than those with haplotype A when 6 months old (P < 0.01), as well as better body height, body length, body weight, and average daily gain at 6 months (P < 0.05). Therefore, haplotype B of the ND5 gene was likely to have a positive effect on growth traits at 6 months in the Nanyang breed.     

East Asian contributions to Dutch traditional and western commercial chickens inferred from mtDNA analysis

Understanding the complex origin of domesticated populations is of vital importance for understanding, preserving and exploiting breed genetic diversity. Here, we aim to assess Asian contributions to European traditional breeds and western commercial chickens for mitochondrial genetic diversity. To this end, a 365‐bp fragment of the chicken mtDNA D‐loop region of 16 Dutch fancy breeds (113 individuals) was surveyed, comprising almost the entire breed diversity of The Netherlands. We also sequenced the same fragment for 160 commercial birds representing all important commercial types from multiple commercial companies that together represent more than 50% of the worldwide commercial value. We identified 20 different haplotypes. The haplotypes clustered into five clades. The commonest clade (E‐clade) supposedly originates from the Indian subcontinent. In addition, both in commercial chicken and Dutch fancy breeds, many haplotypes were found with a clear East Asian origin. However, the erratic occurrence of many different East Asian mitochondrial clades indicates that there were many independent instances where breeders used imported exotic chickens for enhancing local breeds. Nucleotide diversity and haplotype diversity analyses showed the influence of the introgression of East Asian chicken on genetic diversity. All populations that had haplotypes of multiple origin displayed high inferred diversity, as opposed to most populations that had only a single mitochondrial haplotype signature. Most fancy breeds were found to have a much lower within‐population diversity compared to broilers and layers, although this is not the case for mitochondrial estimates in fancy breeds that have multiple origin haplotypes.     

中国部分地方黄牛品种线粒体D-loop区的遗传变异与进化分析

测定了13个黄牛品种125个个体的线粒体D-loop区段的全序列,包括12个中国地方黄牛品种的123个个体和德国黄牛2个个体,并进行了分析。结果显示,共检测到93个变异位点,57个单倍型,平均核苷酸差异(average number ofnucleotide differences,k)为22.708,核苷酸多样度(nucleotide diversity,π)为0.0251±0.00479,单倍型多样度(haplotypediversity,Hd)为0.888±0.026,表明我国黄牛品种遗传多样性非常丰富。构建的Neighbor-Joining进化树显示这13个品种主要分成两大类型:普通牛和瘤牛;新发现的特殊类型Ⅲ只有一个西藏阿沛甲咂牛的个体,它与牦牛D-loop序列最相近,证明西藏地区的黄牛与牦牛之间存在基因渗入现象。普通牛和瘤牛在日喀则驼峰牛中占的比例分别是64.3%和35.7%,在阿沛甲咂牛中占的比例分别是50.0%和50.0%,证明了西藏的黄牛也有瘤牛类型。云南牛品种的单倍型非常丰富证明了云南在中国黄牛起源上的重要地位;在27个中国黄牛品种中(本研究11个品种以及GenBank上的16个品种)找到了中国瘤牛的核心单倍型i1,并且对它进行了讨论。同时证明了西藏瘤牛独立于中国瘤牛核心类群的特殊性。     

Marker assisted evaluation of morphological and genetic attributes of sub-populations of Nili-Ravi buffalo: A vulnerable dairy type riverine breed of India

In the present study, we report the distribution of true to type and atypical Nili-Ravi buffalo, a vulnerable dairy type riverine breed of North India and its underlying genetic structure. Out of total investigated buffaloes 73.5% had bilateral wall eyes while 5.4% had unilateral wall eyes and 21.1% had no wall eyes. 41.15% of Nili-Ravi buffaloes maintained in the breeding farm were having typical true to the type characteristics (both eyes walled, white markings in forehead, muzzle/chin, all the four legs and tail) while only 28.5% of Nili-Ravi buffaloes were true to the type under field conditions. Genotypic data were generated in four groups of Nili-Ravi buffalo (FMTNR—Typical Nili-Ravi from farm; FMANR—Atypical Nili-Ravi from farm; FDTNR—Typical Nili-Ravi from field; FDANR—Atypical Nili-Ravi from field) at 16 microsatellite loci. Comparative genetic analysis of various groups of Nili-Ravi buffaloes with Murrah revealed significant between group differences with an estimated global F ST of 0.063. Pair-wise F ST values ranged from 0.003 (between FDTNR and FDANR) to 0.112 (between FMTNR and FDTNR). Phylogenetic analysis and multi-dimensional scaling revealed clustering of FDTNR and FDANR together while FMTNR and FMANR clustered separately with Murrah in between farm and field Nili-Ravi buffaloes. Based on the results, the paper also proposes three pronged strategy for conservation and sustainable genetic improvement of Nili-Ravi buffalo in India.     

PRNP gene variation in Pakistani cattle and buffaloes

Bovine spongiform encephalopathy (BSE) is a neurodegenerative prion protein misfolding disorder of cattle. BSE is of two types, classical BSE and atypical BSE which in turn is of two types, H-type BSE and L-type BSE. Both H-type BSE and L-type BSE are primarily sporadic prion disorders. However, one case of H-type BSE has recently been associated with E211K polymorphism in the prion protein gene (PRNP). Two polymorphisms in the bovine PRNP are also associated with susceptibility to classical BSE: a 23 bp insertion/deletion (indel) in the PRNP promoter region and a 12 bp indel in the first intron. No information regarding BSE susceptibility in Pakistani cattle is available. The present study aimed at achieving this information. A total of 236 cattle from 7 breeds and 281 buffaloes from 5 breeds were screened for E211K polymorphism and 23 bp and 12 bp indels employing triplex PCR. The E211K polymorphism was not detected in any of the animals studied. The 23 bp insertion allele was underrepresented in studied cattle breeds while the 12 bp insertion allele was overrepresented. Both 23 bp and 12 bp insertion alleles were overrepresented in studied buffalo breeds. Almost 90% of alleles were insertion alleles across all studied buffalo breeds. The average frequency of 23 bp and 12 bp insertion alleles across all studied cattle breeds was found to be 0.1822 and 0.9407, respectively. There were significant differences between Pakistani and worldwide cattle in terms of allele, genotype and haplotype frequencies of 23 bp and 12 bp indels. The higher observed frequency of 12 bp insertion allele suggests that Pakistani cattle are relatively more resistant to classical BSE than European cattle. However, the key risk factor for classical BSE is the dietary exposure of cattle to contaminated feedstuffs.     

Genetic polymorphism of alpha-lactalbumin gene in riverine buffalo.

Alpha-lactalbumin (alpha-LA) is a major whey protein found in milk. Polymorphs of alpha-LA gene are reported to be significantly associated with milk production and constituent traits. Therefore, the present study was undertaken to detect polymorphism in alpha-LA at the genic level and to explore allelic variability at this locus. A total of 196 animals, belonging to four breeds of riverine buffalo viz. Bhadwari, Mehsana, Surti and Murrah were included under the present investigation. Two fragments i.e. 133 bp (Exon 1) and 159 bp (Exon 2) of alpha-LA gene were amplified by polymerase chain reaction and subsequently, single strand confirmation polymorphism (SSCP) study was carried out to identify different allelic pattern and genotypes of the animal included in the study. Both fragment of alpha-LA gene was found to be polymorphic in all the four breeds of riverine buffalo. Number of genotypes and allele varied breed to breed for both the fragments. In case of 133 bp fragment, four alleles A, B, C and D were found among different breeds of buffalo whereas in 159 bp fragment, five alleles namely A, B, C, D and E was found in different breeds. Nucleotide sequence data of different alleles showed the presence of both silent as well as functional mutation leading to variability in polypeptide chain.     

Use of linked loci as individuals or haplotypes for marker-assisted breed assignment

The objective of this study was to use simulation to evaluate the benefits of considering haplotypes of loci when linked single nucleotide polymorphisms are used for breed assignment. Three breeds of 10 000 females each were simulated under eight scenarios that differed according to the number of generations separating the breeds, size of breed founder populations and recombination rate between linked loci. Molecular genotypes consisted of 20 groups of three linked loci each. Breed assignment was performed in the final generation and was based on the frequency method. Haplotypes were reconstructed using the expectation–maximization algorithm. Accuracy of breed assignment was based on the frequency of correct breed assignment. Assignment accuracy increased as more genotypes (loci or haplotypes) were considered and more animals were used to estimate genotypic frequencies within breed. For most scenarios, use of haplotypes yielded equal or greater accuracies than when loci were considered independent. The advantage of haplotypes tended to increase as linkage disequilibrium between adjacent loci increased. The greatest advantage for using haplotypes was observed when recombination rate was low (0.001), breeds were separated by few generations (100), and a relatively large number of founder animals (110) was used to form new breeds. In this situation, 90% accuracy of breed assignment was achieved using nine to 14 haplotypes (i.e. 27–42 loci) depending on breed, vs. 39–57 individual loci.     

Mitochondrial DNA variation in cattle of South China: origin and introgression

Ten restriction endonucleases were used to investigate the mitochondrial DNA restriction fragment length polymorphism (mtDNA RFLP) of 11 native cattle breeds and one cultivated cattle breed in South China. Twenty-three restriction morphs were detected, which can be sorted into five haplotypes. A phylogenetic tree of the haplotypes was constructed by using the 'upgMa' method. Our study showed that haplotype I and II are identical to the zebu (Bos indicus) and taurine (Bos taurus) haplotypes, respectively. Zebu and taurine were the two major origins of cattle populations in South China, and the zebu probably had more influence on the native cattle population than taurine did. Haplotype III is identical to haplotype I of yak (Bos grunniens), which was only detected in the Diqing cattle breed. Haplotype IV was detected for the first time. This haplotype, found only in Dehong cattle, might be from an independent domestication event, probably from another Bos indicus population. Divergence of haplotypes I and IV occurred about 268,000-535,000 years ago, much earlier than the 10,000-year history of cattle husbandry. Our results also suggest a secondary introgession of mtDNA from yak to Diqing cattle.     

Ancestral matrilineages and mitochondrial DNA diversity of the Lidia cattle breed

To clarify the genetic ancestry and the mitochondrial DNA (mtDNA) diversity of the Lidia cattle breed, a 521-bp D-loop fragment was sequenced in 527 animals belonging to 70 herds distributed across 29 lineages. The mtDNA diversity recorded was similar to that seen for Middle Eastern breeds and greater than that recorded for the majority of European breeds. Haplotype T3 was the most common (81%), followed by the African T1 haplotype (17%); very low frequencies were recorded for haplotypes T and T2. The results agree with there being two major ancestral lines for the Lidia breed, European and African, similar to that seen for other Mediterranean breeds. A wide range of variation in haplotype frequencies was seen between the examined lineages. Haplotype T3 was present in all those analysed; in five it was the only one present, and in only one lineage (Miura) was its frequency lower than that of T1. T1*, a haplotype reported in Criollo breeds and to date in only a single European breed (the Retinta breed from Spain), was found in a single animal belonging to the Concha y Sierra lineage. Network analysis of the Lidia breed revealed the presence of two major haplotypes: T3 and T1. The Lidia breed appears to be more closely related to prehistoric Iberian and Italian than to British aurochs.     

Recovery of Native Genetic Background in Admixed Populations Using Haplotypes,Phenotypes, and Pedigree Information – Using Cika Cattle as a Case Breed

The aim of this study was to obtain unbiased estimates of the diversity parameters, the population history, and the degree of admixture in Cika cattle which represents the local admixed breeds at risk of extinction undergoing challenging conservation programs. Genetic analyses were performed on the genome-wide Single Nucleotide Polymorphism (SNP) Illumina Bovine SNP50 array data of 76 Cika animals and 531 animals from 14 reference populations. To obtain unbiased estimates we used short haplotypes spanning four markers instead of single SNPs to avoid an ascertainment bias of the BovineSNP50 array. Genome-wide haplotypes combined with partial pedigree and type trait classification show the potential to improve identification of purebred animals with a low degree of admixture. Phylogenetic analyses demonstrated unique genetic identity of Cika animals. Genetic distance matrix presented by rooted Neighbour-Net suggested long and broad phylogenetic connection between Cika and Pinzgauer. Unsupervised clustering performed by the admixture analysis and two-dimensional presentation of the genetic distances between individuals also suggest Cika is a distinct breed despite being similar in appearance to Pinzgauer. Animals identified as the most purebred could be used as a nucleus for a recovery of the native genetic background in the current admixed population. The results show that local well-adapted strains, which have never been intensively managed and differentiated into specific breeds, exhibit large haplotype diversity. They suggest a conservation and recovery approach that does not rely exclusively on the search for the original native genetic background but rather on the identification and removal of common introgressed haplotypes would be more powerful. Successful implementation of such an approach should be based on combining phenotype, pedigree, and genome-wide haplotype data of the breed of interest and a spectrum of reference breeds which potentially have had direct or indirect historical contribution to the genetic makeup of the breed of interest.     

Iberian origin of Brazilian local pig breeds based on Cytochrome b (MT-CYB) sequence

The aim of this work was to investigate the possible origin of local Brazilian pig breeds through Cytochrome b ( MT-CYB ) mitochondrial analyses. The results indicated that the main local pig breeds descended from two different European maternal lineages, both Iberian varieties. The haplotype relationship analysis showed that Monteiro, Nilo, Piau and Tatu breeds share haplotypes only with Iberian varieties, while the Moura breed presented a different maternal lineage. The Moura appears to share a high frequency of haplotypes with the Black Hairy Iberian variety and Hungarian Mangalica breed.     

Biometry of frozen-thawed sperm from eight breeds of Indian buffaloes (Bubalus bubalis)

Sperm morphometry, in combination with other objective traits, can be useful for developing a fertility index. The objective of the present study was to measure various biometric end points of frozen-thawed sperm from eight breeds of Indian buffaloes (Murrah, Surti, Tarai, Mehsana, Jaffrabadi, Bhadawari, Pandharpuri and Nili-Ravi). The sperm head of Pandharpuri buffaloes had the greatest length (10.21 microm), width (6.05 microm), area (52.31 microm(2)) and perimeter (31.86 microm). The ratio of sperm width to length was also greatest (0.61) in Pandharpuri as well as in two other breeds, viz. Nili-Ravi and Jaffrabadi. Murrah had the smallest sperm head width (4.75 microm), area (41.65 microm(2)) and perimeter (29.17 microm), but its sperm tail was longest (57.02 microm), along with that of Jaffrabadi buffaloes (56.96 microm). Based on mean values of sperm tail length, mid piece length and its width the eight buffalo breeds were categorized into three, four and five groups, respectively. Multivariate analysis and clustering put six breeds (Surti, Tarai, Mehsana, Jaffrabadi, Bhadawari and Nili-Ravi) in one cluster, whereas Murrah and Pandharpuri appeared as separate entities.     

Genetic Variability and Linkage Disequilibrium Patterns in the Bovine <Emphasis Type="Italic">DNAJA1</Emphasis> Gene

Correlation between expression level of the bovine DNAJA1 gene and meat tenderness was recently found in Charolais longissimus thoracis muscle samples, suggesting that this gene could play an important role in meat tenderness. Here, we report the validation of polymorphisms within the bovine DNAJA1 gene, and the haplotype variability and extent of linkage disequilibrium in the three main French beef breeds (Blonde d’Aquitaine, Charolais, Limousin). Genotyping 18 putative SNPs revealed that 16 SNPs were polymorphic within the breeds tested. Two SNPs were removed from further analyses as one SNP had a low genotyping call rate, while the other SNP was not in Hardy–Weinberg equilibrium. The degree of heterozygosity observed for the remaining 14 SNPs varied between breeds, with Charolais being the breed with the highest genetic variation and Blonde d’Aquitaine the lowest. Linkage disequilibrium and haplotype structure of DNAJA1 were different between breeds. Eighteen different haplotypes, including three shared by all breeds, were discovered, and two to three tag SNPs (depending on the breed) are sufficient to capture all the genetic variability seen in these haplotypes. The results of this study will facilitate the design of optimal future association studies evaluating the role of the DNAJA1 gene in meat tenderness.