Selective and ATP-competitive kinesin KIF18A inhibitor suppresses the replication of influenza A virus

The influenza virus is one of the major public health threats. However, the development of efficient vaccines and therapeutic drugs to combat this virus is greatly limited by its frequent genetic mutations. Because of this, targeting the host factors required for influenza virus replication may be a more effective strategy for inhibiting a broader spectrum of variants. Here, we demonstrated that inhibition of a motor protein kinesin family member 18A (KIF18A) suppresses the replication of the influenza A virus (IAV). The expression of KIF18A in host cells was increased following IAV infection. Intriguingly, treatment with the selective and ATP-competitive mitotic kinesin KIF18A inhibitor BTB-1 substantially decreased the expression of viral RNAs and proteins, and the production of infectious viral particles, while overexpression of KIF18A enhanced the replication of IAV. Importantly, BTB-1 treatment attenuated the activation of AKT, p38 MAPK, SAPK and Ran-binding protein 3 (RanBP3), which led to the prevention of the nuclear export of viral ribonucleoprotein complexes. Notably, administration of BTB-1 greatly improved the viability of IAV-infected mice. Collectively, our results unveiled a beneficial role of KIF18A in IAV replication, and thus, KIF18A could be a potential therapeutic target for the control of IAV infection.     

Molecular characterization, polymorphism of bovine ZBTB38 gene and association with body measurement traits in native Chinese cattle breeds

Zinc finger and BTB domain containing 38 (ZBTB38), binding to and repressing methylated DNA, is an important candidate gene for selection of body measurement traits through marker-assisted selection (MAS). The expression of ZBTB38 is regulated in human and animal height as well as other stature indexes. Genomic structural analysis shows that bovine ZBTB38 shares much similarity with human ZBTB38. We discovered and evaluated the potential association of the single nucleotide polymorphism (SNP) of the bovine ZBTB38 gene with body measurement traits in 722 individuals. The latest findings demonstrate that the A841G SNP in exon 1 is significantly associated with Body Length (BL), Hip Height (HH) and Heart Girth (HG). Furthermore, the analysis of A841G SNP marker shows that there are significant effects on the BL (P = 0.0389) in 722 individuals, significant effects on the HH (P = 0.0173) and HG (P = 0.0147) in Qinchuan improvement steers (QI) population, as well as significant effects on the WH (P = 0.0094) in Xuelong (XL) population. These results clearly suggest that the ZBTB38 gene is among of target genes for body measurement traits in bovine reproduction and breeding, and thus provide data for establishment of an animal model using cattle to study big animal body type.     

No evidence for an association between common nonsynonymous polymorphisms in delta and bristle number variation in natural and laboratory populations of Drosophila melanogaster

We test the hypothesis that naturally occurring nonsynonymous variants in the Delta ligand of the Notch signaling pathway contribute to standing variation in sternopleural and/or abdominal bristle number in Drosophila melanogaster, for both a large cohort of wild-caught flies and previously described laboratory lines. We sequenced the transcribed region of Delta for 16 naturally occurring chromosomes and 65 SNPs, including 7 nonsynonymous SNPs (nsSNPs), were observed. Identified nsSNPs and 6 additional common SNPs, all located in exon 6 and the 3' UTR, were genotyped in 2060 wild-caught flies using an OLA-based methodology and genotyped in 38 additional natural chromosomes via DNA sequencing. None of the genotyped nsSNPs were significantly associated with natural variation in bristle number as assessed by a permutation test. A 95% upper bound on the additive genetic variance attributable to each genotyped SNP in the large natural cohort is <2% of the total phenotypic variation. Results suggest that two previously detected genotype/phenotype associations between bristle number and variants in the introns of Delta cannot be explained by linkage disequilibrium between these variants and nearby nonsynonymous variants. Unidentified regulatory variants more parsimoniously explain previous observations.     

Whole-exome sequencing identifies mutations of KIF22 in spondyloepimetaphyseal dysplasia with joint laxity, leptodactylic type

Spondyloepimetaphyseal dysplasia with joint laxity (SEMDJL), leptodactylic (lepto-SEMDJL) or Hall type, is an autosomal-dominant skeletal dysplasia manifesting with short stature, joint laxity with dislocation(s), limb malalignment, and spinal deformity. Its causative gene mutation has not yet been discovered. We captured and sequenced the exomes of eight affected individuals in six unrelated kindreds (three individuals in a family and five simplex individuals). Five novel sequence variants in KIF22, which encodes a member of the kinesin-like protein family, were identified in seven individuals. Sanger sequencing of KIF22 confirmed that c.443C>T (p.Pro148Ser) cosegregated with the phenotype in the affected individuals in the family; c.442C>T (p.Pro148Leu) or c.446G>A (p.Arg149Gln) was present in four of five simplex individuals, but was absent in unaffected individuals in their family and 505 normal cohorts. KIF22 mRNA was detected in human bone, cartilage, joint capsule, ligament, skin, and primary cultured chondrocytes. In silico analysis of KIF22 protein structure indicates that Pro148 and Arg149 are important in maintaining hydrogen bonds in the ATP binding and motor domains of KIF22. We conclude that these mutations in KIF22 cause lepto-SEMDJL.     

Exome sequencing reveals genetic architecture in patients with isolated or syndromic short stature

Short stature is among the most common endocrinological disease phenotypes of childhood and may occur as an isolated finding or in conjunction with other clinical manifestations. Although the diagnostic utility of clinical genetic testing in short stature has been implicated, the genetic architecture and the utility of genomic studies such as exome sequencing(ES) in a sizable cohort of patients with short stature have not been investigated systematically. In this study, we recruited 561 individuals with short stature from two centers in China during a 4-year period. We performed ES for all patients and available parents. All patients were retrospectively divided into two groups: an isolated short stature group(group I, n = 257) and an apparently syndromic short stature group(group II, n = 304). Causal variants were identified in 135 of 561(24.1%) patients. In group I, 29 of 257(11.3%) of the patients were solved by variants in 24 genes. In group II, 106 of 304(34.9%) patients were solved by variants in 57 genes. Genes involved in fundamental cellularprocess played an important role in the genetic architecture of syndromic short stature. Distinct genetic architectures and pathophysiological processes underlie isolated and syndromic short stature.     

FAM111A Mutations Result in Hypoparathyroidism and Impaired Skeletal Development

Kenny-Caffey syndrome (KCS) and the similar but more severe osteocraniostenosis (OCS) are genetic conditions characterized by impaired skeletal development with small and dense bones, short stature, and primary hypoparathyroidism with hypocalcemia. We studied five individuals with KCS and five with OCS and found that all of them had heterozygous mutations in FAM111A. One mutation was identified in four unrelated individuals with KCS, and another one was identified in two unrelated individuals with OCS; all occurred de novo. Thus, OCS and KCS are allelic disorders of different severity. FAM111A codes for a 611 amino acid protein with homology to trypsin-like peptidases. Although FAM111A has been found to bind to the large T-antigen of SV40 and restrict viral replication, its native function is unknown. Molecular modeling of FAM111A shows that residues affected by KCS and OCS mutations do not map close to the active site but are clustered on a segment of the protein and are at, or close to, its outer surface, suggesting that the pathogenesis involves the interaction with as yet unidentified partner proteins rather than impaired catalysis. FAM111A appears to be crucial to a pathway that governs parathyroid hormone production, calcium homeostasis, and skeletal development and growth.     

Genetic studies of blood markers in Przewalski's horses

Ninety-six Przewalski's horses (Equus przewalskii) were blood typed using systems of inherited blood variants known to be highly effective for parentage testing of domestic horses (E. caballus). Sixteen red cell antigenic factors detected using sera prepared by alloimmunization of domestic horses were shown to be inherited in six systems (A, C, D, P, Q, and U) and in the same patterns as domestic horses. Family data confirmed autosomal, codominant inheritance at five loci of serum protein variants (Al, Tf, Xk, Pi, and Es) and three loci of red cell proteins (PGM, PHI, and Hb). One serum protein locus (Gc) and two red cell protein loci (PGD and CA) appeared to be monomorphic. Despite the narrow genetic base and high inbreeding coefficients of captive Przewalski's horses, average heterozygosity calculated over 18 loci was estimated to be 0.320 +/- 0.05, which was similar to that found in five breeds of domestic horses.     

A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between Type 2 Diabetes and Prostate Cancer

There is increasing evidence that pleiotropy, the association of multiple traits with the same genetic variants/loci, is a very common phenomenon. Cross-phenotype association tests are often used to jointly analyze multiple traits from a genome-wide association study (GWAS). The underlying methods, however, are often designed to test the global null hypothesis that there is no association of a genetic variant with any of the traits, the rejection of which does not implicate pleiotropy. In this article, we propose a new statistical approach, PLACO, for specifically detecting pleiotropic loci between two traits by considering an underlying composite null hypothesis that a variant is associated with none or only one of the traits. We propose testing the null hypothesis based on the product of the Z-statistics of the genetic variants across two studies and derive a null distribution of the test statistic in the form of a mixture distribution that allows for fractions of variants to be associated with none or only one of the traits. We borrow approaches from the statistical literature on mediation analysis that allow asymptotic approximation of the null distribution avoiding estimation of nuisance parameters related to mixture proportions and variance components. Simulation studies demonstrate that the proposed method can maintain type I error and can achieve major power gain over alternative simpler methods that are typically used for testing pleiotropy. PLACO allows correlation in summary statistics between studies that may arise due to sharing of controls between disease traits. Application of PLACO to publicly available summary data from two large case-control GWAS of Type 2 Diabetes and of Prostate Cancer implicated a number of novel shared genetic regions: 3q23 (ZBTB38), 6q25.3 (RGS17), 9p22.1 (HAUS6), 9p13.3 (UBAP2), 11p11.2 (RAPSN), 14q12 (AKAP6), 15q15 (KNL1) and 18q23 (ZNF236).     

Analysing the Effect of Mutation on Protein Function and Discovering Potential Inhibitors of CDK4: Molecular Modelling and Dynamics Studies

The cyclin-dependent kinase 4 (CDK4)-cyclin D1 complex plays a crucial role in the transition from the G1 phase to S phase of the cell cycle. Among the CDKs, CDK4 is one of the genes most frequently affected by somatic genetic variations that are associated with various forms of cancer. Thus, because the abnormal function of the CDK4-cyclin D1 protein complex might play a vital role in causing cancer, CDK4 can be considered a genetically validated therapeutic target. In this study, we used a systematic, integrated computational approach to identify deleterious nsSNPs and predict their effects on protein-protein (CDK4-cyclin D1) and protein-ligand (CDK4-flavopiridol) interactions. This analysis resulted in the identification of possible inhibitors of mutant CDK4 proteins that bind the conformations induced by deleterious nsSNPs. Using computational prediction methods, we identified five nsSNPs as highly deleterious: R24C, Y180H, A205T, R210P, and R246C. From molecular docking and molecular dynamic studies, we observed that these deleterious nsSNPs affected CDK4-cyclin D1 and CDK4-flavopiridol interactions. Furthermore, in a virtual screening approach, the drug 5_7_DIHYDROXY_ 2_ (3_4_5_TRI HYDROXYPHENYL) _4H_CHROMEN_ 4_ONE displayed good binding affinity for proteins with the mutations R24C or R246C, the drug diosmin displayed good binding affinity for the protein with the mutation Y180H, and the drug rutin displayed good binding affinity for proteins with the mutations A205T and R210P. Overall, this computational investigation of the CDK4 gene highlights the link between genetic variation and biological phenomena in human cancer and aids in the discovery of molecularly targeted therapies for personalized treatment.     

Refining Susceptibility Loci of Chronic Obstructive Pulmonary Disease with Lung eqtls

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD.     

Linkage relationships of seven enzyme and two pigmentation loci in the snail Biomphalaria glabrata

Crossing experiments with inbred stocks of the snail (Biomphalaria glabrata) demonstrated that variants at two loci determining pigmentation and seven enzyme-determining loci exhibited normal Mendelian segregation ratios in F2 progeny. Among 39 pairwise comparisons for joint segregation, there was evidence of genetic linkage between a locus controlling mantle pigmentation (S) and 6-phosphogluconate dehydrogenase (Pgd) and confirmation of a previously described linkage between esterase-2 (Est-2) and catalase (Cat). Recombination fractions were estimated to be 17 +/- 4 for S-Pgd and 33 +/- 5 for Est-2-Cat. The remaining five loci--Acon-1, Pgm-1, Lap-1, Lap-2, and Pgd--assorted independently. This brings to 17 the number of loci examined for segregation and assortment in this medically important species. As Biomphalaria has a chromosome number n = 18, markers should soon be available for most or all of the linkage groups.     

In-Silico Analyses of Nonsynonymous Variants in the BRCA1 Gene

BReast CAncer gene 1 (BRCA1)—a tumor suppressor gene plays an important role in the DNA repair mechanism. Several BRCA1 variants perturb its structure and function, including synonymous and nonsynonymous single nucleotide polymorphisms (SNPs). In the present study, we performed in-silico analyses of nonsynonymous SNPs (nsSNPs) of the BRCA1 gene. In total, 122 nsSNPs were retrieved from the NCBI SNP database and in-silico analyses were performed using computational prediction tools: SIFT, PROVEAN, Mutation Taster, PolyPhen-2, MutPred, and ConSurf. Of these tools, SIFT, PROVEAN, and Mutation Taster predicted 61 out of 122 nsSNPs as “damaging”, based on structural homology analysis. PolyPhen-2 classified 22 nsSNPs as “probably damaging”. These nsSNPs were further analyzed by MutPred to predict basic molecular mechanisms of amino acid alteration. ConSurf analysis predicted eleven conserved amino acid residues with structural and functional consequences. We identified five amino acid residues in the RING finger domain (L22, C39, H41, C44, and C47) and two in the BRCT domain (P1771 and I1707) with the potential to deter the BRCA1 protein function. This study provides insights into the effect of nsSNPs and amino acid substitutions in BRCA1.

     

Consistent interethnic differences in the distribution of clinically relevant endothelial nitric oxide synthase genetic polymorphisms.

A maldistribution of endothelial nitric oxide synthase (eNOS) genetic variants may explain differences in NO-mediated effects and response to drugs among black and white subjects. While interethnic differences in the distribution of eNOS genetic variants exist in the American population, it is not known whether such interethnic differences exist in other populations. To test this possibility, we examined the distribution of genetic variants of three clinically relevant eNOS polymorphisms (T(-786)C in the promoter, the VNTR in intron 4, and the Glu298Asp variant in exon 7) in 136 black and 154 white subjects from a Brazilian population, which is very heterogeneous. We also estimated the haplotype frequency and evaluated associations between these variants. The Asp298 variant was more common in whites (32.8%) than in blacks (15.1%) (P < 0.004). Similarly, the C(-786) variant was more common in whites (41.9%) than in blacks (19.5%) (P < 0.0004). However, the 4a variant was more common in blacks (32.0%) than in whites (17.9%) (P < 0.003). The most common predicted haplotype in both ethnic groups combined only wild-type variants. While the second most common haplotype in blacks includes the variant 4a and the wild-type variants for the remaining polymorphisms, the second most common haplotype in whites includes the variants Asp298 and C(-786) and the wild-type variant for polymorphism in intron 4. The marked interethnic differences that we found in Brazilians are very similar to those previously reported in Americans. These findings strongly suggest a consistent difference in the distribution of eNOS genetic variants in blacks compared with whites and indicate that the interethnic differences do not vary with geographic origin.     

Prediction of the Damage-Associated Non-Synonymous Single Nucleotide Polymorphisms in the Human MC1R Gene

The melanocortin 1 receptor (MC1R) is involved in the control of melanogenesis. Polymorphisms in this gene have been associated with variation in skin and hair color and with elevated risk for the development of melanoma. Here we used 11 computational tools based on different approaches to predict the damage-associated non-synonymous single nucleotide polymorphisms (nsSNPs) in the coding region of the human MC1R gene. Among the 92 nsSNPs arranged according to the predictions 62% were classified as damaging in more than five tools. The classification was significantly correlated with the scores of two consensus programs. Alleles associated with the red hair color (RHC) phenotype and with the risk of melanoma were examined. The R variants D84E, R142H, R151C, I155T, R160W and D294H were classified as damaging by the majority of the tools while the r variants V60L, V92M and R163Q have been predicted as neutral in most of the programs The combination of the prediction tools results in 14 nsSNPs indicated as the most damaging mutations in MC1R (L48P, R67W, H70Y, P72L, S83P, R151H, S172I, L206P, T242I, G255R, P256S, C273Y, C289R and R306H); C273Y showed to be highly damaging in SIFT, Polyphen-2, MutPred, PANTHER and PROVEAN scores. The computational analysis proved capable of identifying the potentially damaging nsSNPs in MC1R, which are candidates for further laboratory studies of the functional and pharmacological significance of the alterations in the receptor and the phenotypic outcomes.     

Computational Identification of Pathogenic Associated nsSNPs and its Structural Impact in UROD Gene: A Molecular Dynamics Approach

Uroporphyrinogen decarboxylase is a cytosolic enzyme involved in the biosynthetic pathway of heme production. Decreased activity of this enzyme results in porphyria cutanea tarda and hepato erythropoietic porphyria. Nonsynonymous single nucleotide polymorphisms (nsSNPs) alter protein sequence and can cause disease. Identifying the deleterious nsSNPs that contribute to disease is an important task. We used five different in silico tools namely SIFT, PANTHER, PolyPhen2, SNPs&GO, and I-mutant3 to identify deleterious nsSNPs in UROD gene. Further, we used molecular dynamic (MD) approach to evaluate the impact of deleterious mutations on UROD protein structure. By comparing the results of all the five prediction results, we screened 35 (51.47 %) nsSNPs as highly deleterious. MD analysis results show that all the three L161Q, L282R, and I334T deleterious variants were affecting the UROD protein structural stability and flexibility. Our findings provide strong evidence on the effect of deleterious nsSNPs in UROD gene. A detailed MD study provides a new insight in the conformational changes occurred in the mutant structures of UROD protein.     

The 18p11.22 locus is associated with never smoker non-small cell lung cancer susceptibility in Korean populations

The proportion of never smoker non-small cell lung cancer (NSCLC) in Asia is about 30-40%. Despite the striking demographics and high prevalence of never smoker NSCLC, the exact causes still remain undetermined. Although several genome wide association (GWA) studies were conducted to find susceptibility loci for lung cancer in never smokers, no regions were replicated except for 5p15.33, suggesting locus heterogeneity and different environmental toxic effects. To identify genetic loci associated with susceptibility of lung cancer in never smokers, we performed a GWA analysis using the Affymetrix 6.0 SNP array. For discovery GWA set, we recruited 446 never smoking Korean patients with NSCLC and 497 normal subjects. We tested association of SNPs with lung cancer susceptibility using the Cochran-Armitage trend test. For validation, 39 SNPs were selected from the top 50 SNPs and five additional SNPs were selected in the DAB1 gene region which showed significant associations in the GWA analysis. The validation SNPs were genotyped in an independent sample including 434 patients and 1,000 controls. Among the 44 validation SNPs, two SNPs (rs11080466 and rs11663246) near the APCDD1, NAPG and FAM38B genes in the 18p11.22 region were replicated. P value of rs11080466 was 1.08 × 10(-6) in the combined sets (2.68 × 10(-5) in the discovery set and 2.60 × 10(-3) in the validation set) and odds ratio was 0.68 (0.58-0.79). We observed similar association for rs11663246. Our result suggests the 18p11.22 region as a novel lung cancer susceptibility locus in never smokers.     

Mechanistic Phenotypes: An Aggregative Phenotyping Strategy to Identify Disease Mechanisms Using GWAS Data

A single mutation can alter cellular and global homeostatic mechanisms and give rise to multiple clinical diseases. We hypothesized that these disease mechanisms could be identified using low minor allele frequency (MAF<0.1) non-synonymous SNPs (nsSNPs) associated with “mechanistic phenotypes”, comprised of collections of related diagnoses. We studied two mechanistic phenotypes: (1) thrombosis, evaluated in a population of 1,655 African Americans; and (2) four groupings of cancer diagnoses, evaluated in 3,009 white European Americans. We tested associations between nsSNPs represented on GWAS platforms and mechanistic phenotypes ascertained from electronic medical records (EMRs), and sought enrichment in functional ontologies across the top-ranked associations. We used a two-step analytic approach whereby nsSNPs were first sorted by the strength of their association with a phenotype. We tested associations using two reverse genetic models and standard additive and recessive models. In the second step, we employed a hypothesis-free ontological enrichment analysis using the sorted nsSNPs to identify functional mechanisms underlying the diagnoses comprising the mechanistic phenotypes. The thrombosis phenotype was solely associated with ontologies related to blood coagulation (Fisher''s p = 0.0001, FDR p = 0.03), driven by the F5, P2RY12 and F2RL2 genes. For the cancer phenotypes, the reverse genetics models were enriched in DNA repair functions (p = 2×10−5, FDR p = 0.03) (POLG/FANCI, SLX4/FANCP, XRCC1, BRCA1, FANCA, CHD1L) while the additive model showed enrichment related to chromatid segregation (p = 4×10−6, FDR p = 0.005) (KIF25, PINX1). We were able to replicate nsSNP associations for POLG/FANCI, BRCA1, FANCA and CHD1L in independent data sets. Mechanism-oriented phenotyping using collections of EMR-derived diagnoses can elucidate fundamental disease mechanisms.     

ZBTB38 suppresses prostate cancer cell proliferation and migration via directly promoting DKK1 expression

Prostate cancer is still one of the most common malignancies in men all around the world. The mechanism of how prostate cancer initiates and develops is still not clear. Here in this study, we show that tumor suppressor ZBTB38 could suppress the migration and proliferation of prostate cancer cells. We find lower ZBTB38 expression in prostate cancer tissues, which also strongly predicts a poorer prognosis of prostate cancer. ZBTB38 binds DKK1 (Dickkopf WNT signaling pathway inhibitor 1) locus and promotes DKK1 expression in prostate cancer cell lines. Consistently, reduction of DKK1 expression significantly restores ZBTB38-mediated suppression of migration and proliferation of prostate cancer cell lines. Mechanistically, we find that ZBTB38 primarily binds the promoters of target genes, and differentially regulates the expression of 1818 genes. We also identify PRKDC (protein kinase, DNA-activated, catalytic subunit) as a ZBTB38-interacting protein that could repress the function of ZBTB38 in suppressing migration and proliferation of prostate cancer cells. Taken together, our results indicate that ZBTB38 could repress cell migration and proliferation in prostate cancer via promoting DKK1 expression, and also provide evidence supporting ZBTB38 as a potential prognosis marker for prostate cancer.Subject terms: Tumour-suppressor proteins, Prostate cancer     

ZBTB38 gene polymorphism associated with body measurement traits in native Chinese cattle breeds

Body measurement traits, influenced by genes and environmental factors, play numerous important roles in the value assessment of productivity and economy. In this study, we investigated the association between genetic polymorphisms of the zinc finger and BTB domain containing 38 gene (ZBTB38) and body measurement traits in native Chinese cattle. Using direct DNA sequencing in 423 individuals of 8 different cattle subpopulations, 9 novel polymorphisms were identified for genotyping within 647 bp region of exon 1 of ZBTB38. Linkage disequilibrium and association analysis revealed that two coding exon polymorphisms (2323 G>A and 2325 C>T polymorphisms), missense mutations valine GTC(T) to isoleucine ATC(T), were associated with body length (BL), withers height (WH) and rump length (RL). Furthermore, the analysis of 2323 G>A and 2325 C>T SNP markers shows that there are significant effects on the BL (P = 0.0488), WH (P = 0.0044) and RL (P = 0.0314) in the total population. These results clearly suggest that the ZBTB38 gene is among the target genes for body measurement traits in bovine breeding, and provide data for establishment of an animal model using cattle to study big animal body type.