Haplotype structure and phylogenetic shadowing of a hypervariable region in the CAPN10 gene

It has been proposed that variation in calpain 10 (CAPN10) contributes to the risk of type 2 diabetes (T2D). A previous survey of CAPN10 in ethnically diverse populations revealed an intronic region with a significant excess of polymorphism levels relative to inter-species sequence divergence, suggesting that this region was the target of long-standing balancing selection. Based on the thrifty genotype hypothesis, variation that increases risk to T2D in contemporary humans at one time conferred a survival advantage in ancestral populations. Thus, the signature of positive natural selection in a T2D candidate gene could identify a genomic region containing variation that influences disease susceptibility. Here, we investigate this hypothesis by re-sequencing the CAPN10 region with unusual polymorphism levels in T2D cases and controls (n=91) from a Mexican American (MA) population, and by using networks to infer the evolutionary relationships between the major haplotypes. Haplotype tag SNPs (htSNPs) were then selected in each population sample and in MA cases and controls. By placing the htSNPs on the haplotype network, we investigate how cross-population differences in CAPN10 genetic architecture may affect the detection of the disease association. Interestingly, despite the small scale of our case-control study, we observe a nearly significant signal of association between T2D and variation in the putative target of balancing selection. Finally, we use phylogenetic shadowing across 10 primate species to search for conserved non-coding elements that may affect the expression and function of CAPN10. These elements are postulated to be the targets of long-standing balancing selection.     

<Emphasis Type="Italic">CAPN1</Emphasis> markers in three Argentinean cattle breeds: report of a new InDel polymorphism within intron 17

In this study, the genotype distribution and allelic frequencies of CAPN1 (Calcium activated neutral protease) single nucleotide polymorphisms (SNPs) were analyzed taking advantage of the different genetic backgrounds provided by Hereford, Brahman and Braford cattle. We report a new insertion/deletion (InDel) polymorphism, consisting of a change of seven nucleotides for only one nucleotide (TCTGGGT → C) within intron 17 of the CAPN1 gene. The segregation pattern of this polymorphism was analyzed together with the markers CAPN316, CAPN530 and CAPN4751 already described. The allele distribution of CAPN1 markers in the Braford crossbreed (3/8 Brahman 5/8 Hereford) is described for the first time. Four assays of allelic discrimination were designed: the tetra primer ARMS-PCR technique for genotyping the new InDel and the CAPN4751 marker, and a PCR-RFLP method for genotyping the markers CAPN316 and CAPN530. The genotypic and minor allele frequencies (MAFs) obtained showed that the InDel polymorphism does not provide redundant information to that already provided by the other CAPN1 markers and segregates differently between breeds, being a common SNP (MAF ≥ 0.05) in the herds with a high percentage of Bos indicus background. The high percentage of heterozygous individuals found in the Braford crossbreed for the markers assessed reveals enough genetic variation that could help to solve the tenderness problem of tropical-adapted cattle.     

Characteristics of myotonic dystrophy in Istria: molecular genetic approach. Part II: Analysis of genetic polymorphisms

One of the world highest prevalence estimates of myotonic dystrophy (DM) has been reported in the Croatian region Istria. To analyse the population genetic characteristics of DM locus in Istria, two intragenic and three extragenic polymorphic markers were tested. The Southern blot technique was used for D19S63 locus analysis, whereas PCR analysis was performed for CKMM, Alu polymorphism, DMPK (G/T) intron 9/HinfI polymorphism, and D19S207 genetic markers. The compound haplotypes segregating with DM were established. A complete association between the DM mutation and D19S63, D19S207, intron 9/HinfI polymorphism and Alu polymorphism markers were found. In all DM chromosomes: D19S63 and Alu markers had the allele 1 in common; D19S207 had the allele 3 in common, DMPK (G/T) intron 9/HinfI marker had the allele 2 in common. The analysis of CKMM polymorphism revealed genotype heterogeneity; in DM chromosomes either allele 2 or allele 4 were found. The haplotype analysis in the population of Croatian Istria supports the linkage disequilibrium between the DM mutation and Alu polymorphism, intron 9/HinfI polymorphism, D19S63 and D19S207 markers as reported worldwide. The results of the haplotype analysis suggest a common origin of the mutation in Istrian population.     

GPR30 gene polymorphisms are associated with progesterone receptor status and histopathological characteristics of breast cancer patients

G-protein coupled receptor GPR30 has been demonstrated to mediate estrogenic effects on essential features of human breast cancer cells. Polymorphisms in GPR30 gene might therefore affect breast cancer susceptibility or tumor characteristics. This is the first study examining allele and genotype frequencies of GPR30 single nucleotide polymorphisms (SNPs) in breast cancer patients. A total of 257 sporadic breast cancer cases and 247 age-matched controls were genotyped for three GPR30 polymorphisms by means of allele-specific tetra-primer PCR. Comparison of the breast cancer case and the control group with regard to the SNP allele, genotype and haplotype frequencies did not show significant differences. In contrast, the GPR30 SNPs tested were significantly associated with tumor size, histological grading, nodal status and progesterone receptor (PR) status. The A allele of SNP rs3808351 was significantly less frequent in patients with large or G3 tumors, T allele of SNP rs11544331 less frequently occurred in patients with positive nodal status, suggesting that both SNPs might exert protective effects regarding aggressive breast cancer entities. Both homozygous GG genotype of promoter SNP rs3808350 and T allele of missense SNP rs11544331 were inversely associated with PR-negativity, suggesting that they might exert protective effects regarding development of PR-negative cancer. In conclusion, the results of this study support the important role of GPR30 in breast cancer and encourage functional studies on the molecular mechanisms underlying the association of GPR30 polymorphisms with PR status and tumor growth.     

Adiponectin gene haplotype is associated with preeclampsia

We determined whether the polymorphism of the gene encoding adiponectin contributes to susceptibility to preeclampsia. The study involved 133 Finnish women with preeclampsia and 245 healthy control subjects. All women were genotyped for two single nucleotide polymorphisms (SNPs), SNP45 in exon 2 and SNP276 in intron 2, in the adiponectin gene. Chi2 analysis was used to assess genotype and allele frequency differences between the preeclamptic and control groups. In addition, the pair of loci haplotype analysis, using the expectation-maximization (EM) algorithm, was used to examine the estimated haplotype frequencies of the two SNPs, among the two groups. The TT genotype versus the pooled G genotypes in SNP276 was associated with protection against preeclampsia (p = 0.012) at an odds ratio of 0.27 (95% confidence interval [CI]: 0.09-0.80). Also the genotype and allele frequency distributions of SNP276 differed significantly between the preeclampsia group and the control group (p = 0.035 and p = 0.043, respectively). Single-point genotype and allele distributions in SNP45 of the adiponectin gene were not statistically different between the groups. In the haplotype estimation analysis, the pooled G haplotypes versus the TT haplotype were significantly overrepresented in the preeclampsia group (p = 0.042 +/- 0.005). Polymorphisms of the adiponectin gene show a weak, but statistically significant, haplotype association with susceptibility to preeclampsia in Finnish women.     

Are variants in the CAPN10 gene related to risk of type 2 diabetes? A quantitative assessment of population and family-based association studies

The calpain-10 gene (CAPN10) on chromosome 2q37.3 was the first candidate gene for type 2 diabetes (T2D) identified through a genomewide screen and positional cloning. One polymorphism (UCSNP-43: G-->A) and a specific haplotype combination defined by three polymorphisms (UCSNP-43, -19, and -63) were linked to an increased risk of T2D in several populations. To quantitatively assess the collective evidence for the effects of CAPN10 on risk of T2D, we conducted a meta-analysis of both population-based and family-based association studies. We retrieved data from the MEDLINE, PubMed, and Online Mendelian Inheritance in Man databases, as well as from other relevant reports and abstracts published up to July 2003. From a total of 26 studies with primary data (21 population-based studies: 5,013 cases and 5,876 controls; 5 family-based studies: 487 parent-offspring trios), we developed a summary database that contains variables of study design, study population/ethnicity, specific polymorphisms and haplotype combinations in CAPN10, and diabetes-related metabolic phenotypes. For population-based studies, we used both fixed-effects and random-effects models to calculate the pooled odds ratio (OR) and 95% confidence interval (CI) for the associations of CAPN10 genotypes with the risk of T2D. We also calculated weighted mean differences for the associations between CAPN10 and diabetes-related quantitative traits. Under either an additive or a dominant effect model, we found no statistically significant relation between CAPN10 genotypes in the UCSNP-43 locus and T2D risk. However, under a recessive model, individuals homozygous for the common G allele had a statistically significant 19% higher risk of T2D than carriers of the A allele (OR 1.19; 95% CI 1.07-1.33). The association between the 112/121 haplotype combination and T2D risk appeared to be overestimated by several initial small studies with positive findings (OR 1.38; 95% CI 1.04-1.84). After we removed these initial studies, this association became nonsignificant (OR 1.11; 95% CI 0.91-1.35). Moreover, we found no evidence for the associations between the UCSNP-43 G/G genotype and the 112/121 haplotype combination and metabolic phenotypes. Our meta-analysis of family-based studies showed only an overtransmission of the rare allele C in UCSNP-44 from heterozygous parents to their affected offspring with T2D. Our analysis indicates that inadequate statistical power, racial/ethnic differences in frequencies of alleles, haplotypes and haplotype combinations, potential gene-gene or gene-environment interactions, publication bias, and multiple hypothesis testing may contribute to the significant heterogeneity in previous studies of CAPN10 and T2D. Our findings also suggest that both large-scale, well-designed association studies and functional studies are warranted to either reliably confirm or conclusively refute the initial hypothesis regarding the role of CAPN10 in T2D risk.     

High frequency of polymorphism but no mutations found in the GLUT1 glucose transporter gene in NIDDM and familial obesity by SSCP analysis

To evaluate whether a structural defect in the human glucose transporter gene GLUT1 could be involved in the aetiology of insulin resistance, a key factor of non-insulin-dependent diabetes mellitus (NIDDM) and obesity, we performed single-strand conformation polymorphism (SSCP) analysis in 40 subjects (20 NIDDM patients and 20 subjects with familial obesity). The GLUT1 gene, which is involved in basal glucose transport in most tissues, consists of ten exons and encodes a 492 amino acid protein. Population studies have shown a strong association between the X1 allele of an XbaI restriction fragment length polymorphism of the GLUT1 gene and NIDDM. We therefore performed SSCP analysis in NIDDM subjects known to carry at least one X1 allele. Variant SSCP patterns were detected in exons 2, 4, 5 and 9. Sequence analysis of the SSCP variants revealed the presence, in all exons examined, of silent mutations consisting of single-nucleotide substitutions with no amino acid changes. Both NIDDM and obese patients showed a high frequency of polymorphism in the sequence (50% and 35%, respectively). We conclude that the GLUT1 gene is unlikely to play a role in the aetiology of NIDDM and obesity. However, the strong association between the GLUT1 gene and NIDDM, together with the recent family studies showing linkage between chromosome 1p and NIDDM warrant further studies on this chromosomal region. Received: 18 August 1997 / Accepted: 10 December 1997     

Selection of genetic markers for association analyses,using linkage disequilibrium and haplotypes

The genotyping of closely spaced single-nucleotide polymorphism (SNP) markers frequently yields highly correlated data, owing to extensive linkage disequilibrium (LD) between markers. The extent of LD varies widely across the genome and drives the number of frequent haplotypes observed in small regions. Several studies have illustrated the possibility that LD or haplotype data could be used to select a subset of SNPs that optimize the information retained in a genomic region while reducing the genotyping effort and simplifying the analysis. We propose a method based on the spectral decomposition of the matrices of pairwise LD between markers, and we select markers on the basis of their contributions to the total genetic variation. We also modify Clayton's "haplotype tagging SNP" selection method, which utilizes haplotype information. For both methods, we propose sliding window-based algorithms that allow the methods to be applied to large chromosomal regions. Our procedures require genotype information about a small number of individuals for an initial set of SNPs and selection of an optimum subset of SNPs that could be efficiently genotyped on larger numbers of samples while retaining most of the genetic variation in samples. We identify suitable parameter combinations for the procedures, and we show that a sample size of 50-100 individuals achieves consistent results in studies of simulated data sets in linkage equilibrium and LD. When applied to experimental data sets, both procedures were similarly effective at reducing the genotyping requirement while maintaining the genetic information content throughout the regions. We also show that haplotype-association results that Hosking et al. obtained near CYP2D6 were almost identical before and after marker selection.     

Effects of genetic variants and mapping assignments of the ovine calpain regulatory subunit gene on chromosome 14

The aim of this study was to search genetic variants, to map the ovine calpain regulatory subunit gene (CAPN4S) on chromosomes, and to find associations between single nucleotide (nt) polymorphisms (SNPs) and weight traits. Genotyping was performed with 226 purebred, 133 crossbred and 122 animals from eight sheep reference families. The polymerase chain reaction–single strand conformational polymorphism analysis detected genetic variants in the genomic regions between exons 4 and 8. Nt polymorphisms were confirmed at positions 115 (G/C) with CAPN445, 71 (A/G) with CAPN456, and 25 (A/G) and 44 (A/T) with the CAPN478 segments. The Polypay and Targhee purebreds showed significantly different allele frequencies while the mixed population was placed between pure breeds. The two-point linkage analysis using the CRI-MAP version 2.4 identified locations of the CAPN4S segments, showing that CAPN4S is closely linked with the LSCV30 microsatellite marker on chromosome 14. The genotypes of g.45923178G>C had significant mean differences for BWT and WWT, meaning that the identified SNP may be used in marker assisted selection programs. A significant haplo-block was constructed between g.45924950T>C and g.45924969A>T, and major frequencies were observed as AA (0.727) and GT (0.267).     

Association between insertion-deletion polymorphism of the angiotensin-converting enzyme gene and development of angiopathies in patients with non-insulin dependent diabetes mellitus from the Chuvash Republic

Insertion/deletion polymorphism of the angiotensin-converting enzyme (ACE) gene was analyzed in patients with non-insulin-dependent diabetes mellitus (NIDDM) and in the control group consisting of healthy subjects. The insertion allele (I) and genotype II were found to be associated with NIDDM. The frequencies of diabetic retinopathy and nephropathy in NIDDM patients were not associated with this polymorphism. However, an association was found between the DD genotype of the ACE gene and diabetic angiopathy in lower extremities.     

APOC3基因SstⅠ单核苷酸多态性与冠心病、Ⅱ型糖尿病合并高甘油三酯血症的关联性研究

对许多人群研究表明 ,位于APOA1/C3/A4 /A5基因簇上的载脂蛋白C3基因 (APOC3)SstⅠ多态性与高甘油三酯血症 (Hypertriglyceridaemia ,HTG)密切相关 ,高甘油三酯是冠心病和糖尿病的独立危险因素。为探讨中国人群APOC3基因SstⅠ单核苷酸多态性与冠状动脉粥样硬化性心脏病 (coronaryatheroscleroticheartdisease,CHD)合并高甘油三酯血症 (HTG)、Ⅱ型糖尿病 (non insulin dependentdiabetesmellitus,NIDDM)合并高甘油三酯血症 (HTG)患者的相关性 ,应用聚合酶链反应 限制性片段长度多态性 (PCR RFLP)的方法 ,分析了 2 6 7例CHD患者、2 4 6例NIDDM患者及 4 91例健康对照APOC3基因SstⅠ位点 (S1/S2 )多态性。CHD组、NIDDM组和对照组的APOC3基因SstⅠ多态位点S2等位基因频率分别为 0 30 1、0 30 7和 0 2 86 ,其基因型频率和等位基因频率分布与对照组比较均无显著性差异 (P >0 0 5 )。以TG >1 90mmol/L为标准将CHD组、NIDDM组分为正常甘油三酯组 (NTG)和高甘油三酯组(HTG)发现 ,在CHD患者 ,HTG亚组S1S2基因型频率显著高于NTG亚组 (0 5 4 2 >0 35 7,χ2 =8 77,P =0 0 12 4 ) ;在NIDDM患者 ,HTG亚组S2S2基因型频率显著高于NTG亚组 (0 2 0 0 >0 0 5 5 ,χ2 =2 0 2 1,P =0 0 0 0 0 ) ,两亚组间等位基因频     

Contrasting evolutionary histories of two introns of the duchenne muscular dystrophy gene, Dmd, in humans

The Duchenne muscular dystrophy (Dmd) locus lies in a region of the X chromosome that experiences a high rate of recombination and is thus expected to be relatively unaffected by the effects of selection on nearby genes. To provide a picture of nucleotide variability at a high-recombination locus in humans, we sequenced 5. 4 kb from two introns of Dmd in a worldwide sample of 41 alleles from Africa, Asia, Europe, and the Americas. These same regions were also sequenced in one common chimpanzee and one orangutan. Dramatically different patterns of genetic variation were observed at these two introns, which are separated by >500 kb of DNA. Nucleotide diversity at intron 44 pi = 0.141% was more than four times higher than nucleotide diversity at intron 7 pi = 0.034% despite similar levels of divergence for these two regions. Intron 7 exhibited significant linkage disequilibrium extending over 10 kb and also showed a significant excess of rare polymorphisms. In contrast, intron 44 exhibited little linkage disequilibrium and no skew in the frequency distribution of segregating sites. Intron 7 was much more variable in Africa than in other continents, while intron 44 displayed similar levels of variability in different geographic regions. Comparison of intraspecific polymorphism to interspecific divergence using the HKA test revealed a significant reduction in variability at intron 7 relative to intron 44, and this effect was most pronounced in the non-African samples. These results are best explained by positive directional selection acting at or near intron 7 and demonstrate that even genes in regions of high recombination may be influenced by selection at linked sites.     

Type 2 Diabetes and Three Calpain-10 Gene Polymorphisms in Samoans: No Evidence of Association

Although genomewide scans have identified several potential chromosomal susceptibility regions in several human populations, finding a causative gene for type 2 diabetes has remained elusive. Others have reported a novel gene, calpain-10 (CAPN10), located in a previously identified region on chromosome 2q37.3, as a putative susceptibility gene for type 2 diabetes. Three single-nucleotide polymorphisms (SNPs) (UCSNP43, UCSNP19, and UCSNP63) were shown to be involved in increased risk of the disease among Mexican Americans. We have tested the association of these three SNPs with type 2 diabetes among the Samoans of Polynesia, who have a very high prevalence of the disease. In the U.S. territory of American Samoa, prevalence is 25% and 15% in men and women, respectively, whereas, in the independent nation of Samoa, prevalence is 3% and 5% in men and women, respectively. In our study sample, which consisted of 172 unrelated affected case subjects and 96 control subjects, we failed to detect any association between case subjects and control subjects in allele frequencies, haplotype frequencies, or haplotype combinations of UCSNP43, -19, and -63. Also, our data showed no evidence of linkage, among 201 affected sib pairs, in the region of chromosome 2 that contains these SNPs. Three plausible scenarios could explain these observations. (1) CAPN10 is a susceptibility gene only in particular ethnic groups; (2) our study lacks power to detect the effects of CAPN10 polymorphisms (but our sample size is comparable to that of earlier reports); or (3) the underlying biological mechanism is too complex and requires further research.     

Association between the Insertion/Deletion Polymorphism of the Angiotensin-Converting Enzyme Gene and Angiopathy in Patients with Non-Insulin Dependent Diabetes Mellitus in Chuvash Republic

Insertion/deletion polymorphism of the angiotensin-converting enzyme (ACE) gene was analyzed in patients with non-insulin-dependent diabetes mellitus (NIDDM) and in the control group consisting of healthy subjects. The insertion allele (I) and genotype II were found to be associated with NIDDM. The frequencies of diabetic retinopathy and nephropathy in NIDDM patients were not associated with this polymorphism. However, an association was found between the DD genotype of the ACE gene and diabetic angiopathy in lower extremities.     

Calpain-10 variants and haplotypes are associated with polycystic ovary syndrome in Caucasians

PCOS is known to be associated with an increased risk of T2DM and has been proposed to share a common genetic background with T2DM. Recent studies suggest that the Calpain-10 gene (CAPN10) is an interesting candidate gene for PCOS susceptibility. However, contradictory results were reported concerning the contribution of certain CAPN10 variants, especially of UCSNP-44, to genetic predisposition to T2DM, hirsutism, and PCOS. By means of MALDI-TOF MS technique, we genotyped an expanded single nucleotide polymorphism panel, including the CAPN10 UCSNP-44, -43, -56, ins/del-19, -110, -58, -63, and -22 in a sample of 146 German PCOS women and 606 population-based controls. Statistical analysis revealed an association between UCSNP-56 and susceptibility to PCOS with an odds ratio (OR) of 2.91 (95% CI=1.51-5.61) for women carrying an AA genotype compared with GG. As expected, the 22-genotype of the ins/del-19 variant, which is in high linkage disequilibrium (r2=0.98) with UCSNP-56, was also significantly associated (OR=2.98, 95% CI=1.55-5.73). None of the additionally tested variants alone showed any significant association with PCOS. A meta-analysis including our study (altogether 623 PCOS cases and 1,224 controls) also showed significant association only with ins/del-19. The most common haplotype TGG3AGCA was significantly associated with a lower risk for PCOS (OR=0.487, P=0.0057). In contrast, the TGA2AGCA haplotype was associated with an increased risk for PCOS (OR=3.557, P=0.0011). By investigating a broad panel of CAPN10 variants, our results pointed to an allele dose-dependent association of UCSNP-56 and ins/del-19 with PCOS.     

SimPed: a simulation program to generate haplotype and genotype data for pedigree structures

With the widespread availability of SNP genotype data, there is great interest in analyzing pedigree haplotype data. Intermarker linkage disequilibrium for microsatellite markers is usually low due to their physical distance; however, for dense maps of SNP markers, there can be strong linkage disequilibrium between marker loci. Linkage analysis (parametric and nonparametric) and family-based association studies are currently being carried out using dense maps of SNP marker loci. Monte Carlo methods are often used for both linkage and association studies; however, to date there are no programs available which can generate haplotype and/or genotype data consisting of a large number of loci for pedigree structures. SimPed is a program that quickly generates haplotype and/or genotype data for pedigrees of virtually any size and complexity. Marker data either in linkage disequilibrium or equilibrium can be generated for greater than 20,000 diallelic or multiallelic marker loci. Haplotypes and/or genotypes are generated for pedigree structures using specified genetic map distances and haplotype and/or allele frequencies. The simulated data generated by SimPed is useful for a variety of purposes, including evaluating methods that estimate haplotype frequencies for pedigree data, evaluating type I error due to intermarker linkage disequilibrium and estimating empirical p values for linkage and family-based association studies.     

Identification of copy number variation of CAPN10 in Thais with type 2 diabetes by multiplex PCR and denaturing high performance liquid chromatography (DHPLC)

Copy number variations (CNVs) have been shown to be associated with several diseases. They can cause deviation of genotypes from Hardy-Weinberg Equilibrium (HWE). Genetic case-control association studies in Thais revealed that genotype distribution of CAPN10 Indel19 was deviated from HWE after correction of genotyping error. Therefore, we aim to identify CNVs within CAPN10 Indel19 region. The semi-quantitative denaturating high performance liquid chromatography (DHPLC) method was used to detect CNVs in the region of CAPN10 Indel19 marker in cohort of 305 patients with type 2 diabetes and 250 control subjects without diabetes. CNVs in the region of CAPN10 Indel19 was successfully detected by DHPLC. After correction of genotype calling based on the status of identified CNVs, CAPN10 Indel19 genotypes were well-fitted for HWE (p>0.05). However, we did not find association between CNV genotypes and risk of type 2 diabetes in our population. CNVs in CAPN10 have been identified in Thais. These CNVs lead to deviation from HWE of CAPN10 Indel19 genotypes. After excluding identified CNVs from the analysis, CAPN10 Indel19 was associated with type 2 diabetes. The information obtained from our study would be helpful for genotyping accuracies of SNPs residing in the CNVs region.     

Genetic Variant of MYLK4 Gene and its Association with Growth Traits in Chinese Cattle

As a member of MYLK family, MYLK4 gene may play a vital role in muscle development. In this study, one novel single-nucleotide polymorphism (SNP) was identified the bovine MYLK4 by sequencing pooled DNA samples (pool-Seq) and forced polymerase chain reaction-restriction fragment length polymorphism (forced PCR-RFLP) methods. Overall, we reported one mutation (SNP1) in the intron 10 region within the bovine MYLK4 gene in 559 individuals representing five main cattle breeds from China (Nanyang, NY; Qinchuan; Jiaxian, JX; Pinan cattle; and Caidamu cattle, CDM). Genotype AA and allele A were predominant in the QC, PN, and XN populations. Association analysis with growth traits in the QC breed showed that the animals with genotype GG had significantly greater chest breadth and hip width (P?<?0.05). Meanwhile, the genotype GG was strongly associated with withers height and body length than those with genotype AA (P?<?0.01 or P?<?0.05) at 12 months in the NY breed. These statistical results exhibited that the MYLK4 gene might be a potential candidate gene to improve cattle’s growth traits, and the SNP could be used as molecular markers in early marker-assisted selection (MAS) in beef cattle breeding program.