Effects of SNPs and alternative splicing within HGF gene on its expression patterns in Qinchuan cattle
© Cai et al. 2015
Received: 7 April 2015
Accepted: 9 December 2015
Published: 22 December 2015
Identification of genetic variants, including SNPs (Single Nucleotide Polymorphisms), CNVs (Copy Number Variations) and alternative splicing, within functional genes has received increasing attention in animal science research. HGF (Hepatocyte Growth Factor) is a very important growth factor that works as a mitogen or a morphogen during tissue growth, development and regeneration. However, to date, the functions of genetic variants within the bovine HGF gene, particularly their effects on mRNA expression, have not been determined well.
The present study aimed to perform association analysis between genetic variants and mRNA expression for the bovine HGF gene in Qinchuan cattle using various strategies, including PCR-RFLP (Restriction Fragment Length Polymorphism), qPCR (Quantitative Real-time quantitative PCR), TA cloning, DNA sequencing and bioinformatics analysis. A total of five SNPs were identified and only SV1 locus significantly affected HGF mRNA expression in fetal skeletal muscle (P < 0.05). Heterozygous genotype individuals showed significantly higher HGF expression (P < 0.05), which was significantly greater in the “CTCCAGGGTT” combined genotype than that in the “CCCCGGGGTT” combined genotype (P < 0.05). In addition, two alternative splicing variations, HGF-W and HGF-M, were identified, which resulted from alternative 3′ splice sites of exon 5, and HGF-W showed higher mRNA levels than HGF-M in all tissues.
In summary, genetic variations within the HGF gene affected mRNA expression. These findings provide new insight into the molecular characteristics and functions of bovine HGF.
KeywordsAlternative splicing Expression HGF Qinchuan cattle SNPs
HGF, firstly described as a potent mitogen for hepatocytes, has mitogenic , motogenic, morphogenic , anti-apoptotic  and anti-fibrotic  activities in various cell types, via c-Met, the specific receptor of HGF . Three signaling pathways, including PI3K (phosphatidy linositol 3-kinase), STAT (signal transducer and activator of transcription) and MAPK (mitogen activated protein kinase) [6–8], regulating growth and development in mammals and are related to the HGF-Met pathway. However, the functions of the bovine HGF remain unclear.
SNPs are commonly used as molecular markers. Exonic SNPs have a direct effects on the properties of proteins, while SNPs within introns and untranslated region can affect the expression and splicing of mRNA [9, 10]. Numerous studies have shown that SNPs in the HGF gene are associated with normal growth and development [11–13]. However, it remains unclear that how these mutations affect gene expression, including the bovine HGF gene examined in this study.
Alternative splicing, sequences of pre-mRNA could be alternatively included into the mature mRNA or removed, has been widely examined in genomes studies. Increasing evidences have demonstrated that alternative splicing contributes to protein localization, enzymatic properties, and protein binding domains [14–16], and so on. Three transcript variants, X1 (XM_005205314), X2 (XM_005205315) and X3 (XM_005205316), were predicted to be present in Hereford cattle. Moreover, biologically active rat HGF protein contains 697 or 692 amino acids , indicating that there are two alternative splicing transcripts of HGF in mice. However, there have been no reports regarding HGF splicing variations in Chinese cattle.
In this study, firstly, the effects of SNPs in the bovine HGF gene on mRNA expression were conformed and found that variants within bovine HGF gene were associated to its mRNA levels. Secondly, splicing variants and their tissue expression were investigated. These results provide a foundation for further studies on the bovine HGF.
Sample collection and genomic DNA/total RNA isolation
Skeletal muscle samples from 48 male and 37 female fetal Qinchuan cattle were obtained for genomic DNA and total RNA isolation to detect the correlation between mutations and relative mRNA expression levels. Different tissue samples, including skeletal muscle, liver, heart, kidney, lung and spleen from fetal (5 months of gestation) and adult (24 months old) of Qinchuan cattle (n = 5, respectively) were collected for total RNA isolation to evaluate tissue expression.
Total RNA was extracted from different tissues using the Trizol kit (Takara, Japan). cDNA was synthesized according to the PrimeScript RT Reagent Kit (Perfect Real Time) (Takara, Japan). Genomic DNA from skeletal muscle was extracted by proteinase K digestion, chloroform extraction and absolute ethanol precipitation.
Primers used for mutation detecting, cloning and expression survey
Identification of bovine HGF splicing variants and bioinformatics analysis
Sixty cDNA samples were collected from 6 tissues of 10 cattle individuals (5 fetal and 5 adults) to prepare a cDNA pool, which was used as a template for amplification using the primer P-HGF-CDS (Table 1), which covers the whole coding region of bovine HGF gene based on the reference sequence (NM_001031751.1). The product was purified and ligated into the pGEM-T easy vector (Promega, USA) and surveyed by sequencing. Sequences were compared with the reference sequence by BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). SMART (http://smart.embl-heidelberg.de) was used to predict the protein domains of bovine HGF.
Quantitative Real-time quantitative PCR (qPCR)
qPCR was performed using gene specific primers designed by Primer 5.0, with glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as an endogenous control (Table 1). The PCR specificity of these primers was evaluated using Primer-Blast (http://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome). The reaction was carried out with SYBR premix ExTaq II (Takara, Japan) and a CFX 96™ Real Time Detection System (Bio-Rad, USA) according to the manufacturer’s instruction. The qPCR was performed in 20 μL of reaction mixtures consisting of 50 ng genomic DNA, 10 μL SYBR Premix ExTaq II (Takara, Japan), and 20 pmol each primer. Thermal cycling conditions were: 95 °C for 1 min, followed by 40 cycles of 95 °C for 10 s, 60 °C for 30 s. Melting curve analysis was performed to confirm specific amplification. Each PCR was performed in triplicate.
The relative expression ratios were calculated using the following formula as described by Schmittgen and Livak the 2-ΔΔCt method , in which ΔΔCt = ΔCt(test) ‐ ΔCt(Calibrator), ΔCt(test) = Ct(target gene) ‐ Ct(reference gene), ΔCt(Calibrator) = Ct(target Calibrator) ‐ Ct(reference Calibrator). The data were expressed as the mean ± SE (Standard Error).
After genotyping, SPSS 18.0 software was applied to analyze the association betwween genotypes and genotypic combinations with expression levels in Qinchuan cattle using one-way analysis of variance. The following adjusted linear model as follow was applied: Yil = μ + Gi + Sl + GiSl + eil, where Yil is the expression; μ is the overall mean; Gi is the effect of genotypes when this model was used to analyze the association between a single locus and mRNA expression level, and Gi is the effect of genotypic combinations when this model was used for the association analysis between combined genotypes and expression level; Sl is the effect of sex; and eil is the residual effect. The differences in expression levels between two alternative splicings in every tissue was analyzed using Student t test. The threshold of significance was P < 0.05.
HGF mRNA expression based on genotypes in skeletal muscle of fetal bovine
Identification and verification of alternative splicing in bovine HGF
SMART analysis showed no differences in protein domain types between HGF-W and HGF-M containing an N-terminal transmembrane segment, followed successively by a PAN_AP domain, four tandem Kringle domains, and a C- terminal Tryp_SPc domain, which was similar to the reported HGF structures in other species . The bovine HGF-W sequence was used as a query sequence to search for similar sequences using BLAST. The amino acid sequence of Qinchuan cattle HGF shared 99 % identity with all tested species, indicating that HGF protein was highly conserved during evolution. Furthermore, splicing events are also observed in human (NP_001010932.1 and NP_000592.3), dog (XP_005630939.1 and NP_001002964.1) and other species.
Quantification of alternative splicing in bovine HGF
To evaluate the tissue expression of HGF-W and HGF-M, two pairs of specific primers (P-qHGF-W and P-qHGF-M) were designed (Fig. 3b and Table 1). As shown in Fig. 2c–d, these two splicing variants were measured in six tested tissues of fetal and adult cattle. Interestingly, the expression of HGF-W was higher than that of HGF-M in all tested tissues, suggesting that HGF-W was the domain alternative splicing variant. In particular, the statistical differences were significant in the heart (P = 0.046), liver (P = 0.006), spleen (P = 0.025), lung (P = 0.030) in fetal bovine. However, in adult bovine, a statistical difference was observed only in the spleen (P = 0.015). The expression level of HGF-W in the liver was the highest among the expression levels of tissues in both fetal and adult cattle. HGF-M expression was the highest in the spleen of fetal cattle, but in the lung of adult cattle. In contrast, the lowest expression was observed in the heart in both stages. These two alternative splicing events within the bovine HGF gene resulted in varying expression levels among different tissues.
Bioinformatics analysis showed that HGF was highly conserved across species, and thus it was hypothesized that this protein plays a vital role in the growth and development of Qinchuan cattle. Given the emerging roles of HGF in normal growth and development, the identification of genetic variations, transcript variations and expression of HGF appears an essential step for further studies on HGF in bovine growth and development.
A large number of studies have suggested that gene polymorphisms are associated with expression, which was also confirmed in this study. Five SNPs in the bovine HGF gene were identified, only one SNP in intron 1 was significantly associated with HGF expression in Qinchuan cattle. In some genes, intron 1 is in the range of promoter and cis-regulating elements, which may play an essential role in controlling transcription and expression . Thus, SV1 may affect HGF mRNA expression.
It is well-known that all muscle fibers, the structural units of skeletal muscle, are formed in the prenatal stage . During this period, HGF helps myogenic precursor cells to move from somite to limb buds. Postnatal muscle growth invovles muscle fiber hypertrophy, hyperplasia  and regeneration . During regeneration, satellite cells may be activated, resulting in the expression of HGF and its receptor . Skeletal muscle is a major part of the animal body and is related to animal health, growth traits and animal products . Thus, the fetal skeletal muscle was used to estimate the impact of genotypes on mRNA expression. The result showed that SV1 was significantly related to the expression level of the HGF gene. Additionally, SV1 was significantly associated with growth traits in Chinese cattle . These resluts suggested that SV1 affects HGF mRNA expression in fetal skeletal muscle, and then influences the function of bovine HGF and growth traits. These findings provide insight into the molecular mechanism of HGF post-transcriptional regulation.
Alternative splicing is predicted to exist in more than 90 % human genes . Since alternative splicing is a primary source of biological complexity in mammals , two transcript variants and their expression of HGF was identified. Various mechanisms are involved in alternative splicing: alternative 5′ splice selection, alternative 3′ splicing selection, exon skipping and intron retention . Sequence alignment indicated that HGF-M is subject to exon skipping and results from the alternative 3′ splicing site selection of exon 5. According to previous studies, two main mechanisms create this 3′ splice site: random mutations in intronic sequences and RNA editing, which can lead to spliceosome recognition and exonization . However, no other mutations were found in the splicing region according to the DNA pool sequencing result  and resequencing analysis [unpublished observations, Xu et al.]. In contrast, AG was located in the 3′ end of splicing region (Fig. 3), which was the same as the 3′ splice acceptor site in exon-intron boundaries of the bovine HGF gene. Therefore, this 3′ splice site in bovine HGF may have been induced by RNA editing.
qPCR showed that two splicing variants were ubiquitously expressed and that the expression of HGF-W was higher than HGF-M in all examined tissues. According to an NCBI search, the two alternative splicing variations exist in many species, including NP_001010932.1 and NP_000592.3 in humans. In SMART analysis, the deduced protein domains of HGF-W were compared to that of HGF-M. The result showed that the splicing fragment was located in 4 Kringle domains, which were connected via the HGF receptor, Met . Thus, alternative splicing of HGF-M may influence the HGF-Met pathway by changing the combination of HGF to Met. This hypothesis is consistent with the results of previous studies examining alternative splicing altering the binding to other proteins.
In summary, genetic variations within the bovine HGF gene significantly affect gene expression and then works on gene final function. This is the first report describing two bovine HGF transcript variants in Qinchuan cattle, which were confirmed without tissue specificity and species specificity. These results increasing the understanding of molecular mechanisms of bovine HGF.
one-way analysis of variance
Copy number variations
Hepatocyte Growth Factor
mitogen activated protein kinase
phosphatidy linositol 3-kinase
Quantitative Real-time quantitative PCR
Single Nucleotide Polymorphisms
signal transducer and activator of transcription
This study was supported by the National Natural Science Foundation of China (No. 31272408), Program of National Beef Cattle Industrial Technology System (NO. CARS-38) and National 863 Program of China (Grant No. 2013AA102505).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Nakamura T, Nawa K, Ichihara A. Partial purification and characterization of hepatocyte growth factor from serum of hepatectomized rats. Biochem Biophys Res Commun. 1984;122:1450–9.PubMedView ArticleGoogle Scholar
- Birchmeier C, Gherardi E. Developmental roles of HGF/SF and its receptor, the c-Met tyrosine kinase. Trends Cell Biol. 1998;8:404–10.PubMedView ArticleGoogle Scholar
- Kosai KI, Matsumot K, Nagata S, Tsujimoto Y, Nakamura T. Abrogation of Fas-induced fulminant hepatic failure in mice by hepatocyte growth factor. Biochem Biophys Res Commun. 1998;244:683–90.PubMedView ArticleGoogle Scholar
- Giebeler A, Boekschoten MV, Klein C, Borowiak M, Birchmeier C, Gassler N, et al. c-Met confers protection against chronic liver tissue damage and fibrosis progression after bile duct ligation in mice. Gastroenterol. 2009;137:297–308.View ArticleGoogle Scholar
- Naldini L, Vigna E, Narsimhan RP, Gaudino G, Zarnegar R, Michalopoulos GK, et al. Hepatocyte growth factor (HGF) stimulates the tyrosine kinase activity of the receptor encoded by the proto-oncogene c-MET. Oncogene. 1991;6:501–4.PubMedGoogle Scholar
- Ramos-Nino ME, Blumen SR, Sabo-Attwood T, Pass H, Carbone M, Testa JR, et al. HGF Mediates Cell Proliferation of Human Mesothelioma Cells through a PI3K/MEK5/Fra-1 Pathway. Am J Respir Cell Mol Biol. 2008;38:209–17.PubMedPubMed CentralView ArticleGoogle Scholar
- Takeda K, Akira S. STAT family of transcription factors in cytokine-mediated biological responses. Cytokine Growth Factor Rev. 2000;11:199–207.PubMedView ArticleGoogle Scholar
- Grotegut S, Schweintz D, Christofori G, Lehembre F. Hepatocyte growth factor induces cell scattering through MAPK/Egr-1-mediated upregulation of Snail. EMBO J. 2006;25:3534–45.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang X, Li T, Zhao HB, Khatib H. A mutation in the 3′ untranslated region diminishes microRNA binding and alters expression of the OLR1 gene. J Dairy Sci. 2013;96:6525–8.PubMedView ArticleGoogle Scholar
- Wang X, Zhong J, Gao Y, Ju Z, Huang J. A SNP in intron 8 of CD46 causes a novel transcript associated with mastitis in Holsteins. BMC Genomics. 2014;15:630.PubMedPubMed CentralView ArticleGoogle Scholar
- Hoshino K, Satoh T, Kawaguchi Y, Kuwana M. Association of hepatocyte growth factor promoter polymorphism with severity of interstitial lung disease in Japanese patients with systemic sclerosis. Arthritis Rheumatol. 2011;63:2465–72.View ArticleGoogle Scholar
- Kaya M, Çokakli M, Berk AT, Yaman A, Yesilirmak D, Kumral A, et al. Associations of VEGF/VEGF-Receptor and HGF/c-Met promoter polymorphisms with progression/regression of retinopathy of prematurity. Curr Eye Res. 2013;38:137–42.PubMedView ArticleGoogle Scholar
- Cai H, Lan X, Li A, Zhou Y, Sun J, Lei CZ, et al. SNPs of bovine HGF gene and their association with growth traits in Nanyang cattle. Res Vet Sci. 2013;95:483–8.PubMedView ArticleGoogle Scholar
- Unoki M, Shen JC, Zheng ZM, Harris CC. Novel splice variants of ING4 and their possible roles in the regulation of cell growth and motility. J Biol Chem. 2006;281:134677–86.View ArticleGoogle Scholar
- Lira ME, Loomis AK, Paciga SA, Lloyd DB, Thompson JF. Expression of CETP and of splice variants induces the same level of ER stress despite secretion efficiency differences. J Lipid Res. 2008;49:1955–62.PubMedView ArticleGoogle Scholar
- Vendel AC, Terry MD, Striegel AR, Iverson NM, Leuranguer V, Rithner CD, et al. Alternative splicing of the voltage-gated Ca2+ channel β4 subunit creates a uniquely folded N-terminal protein binding domain with cell-specific expression in the cerebellar cortex. J Neurosci. 2006;26:2635–44.PubMedView ArticleGoogle Scholar
- Nakamura T, Nishizawa T, Hagiya M, Seki T, Shimonishi M, Sugimura A, et al. Molecular cloning and expression of human hepatocyte growth factor. Nature. 1989;342:440–3.PubMedView ArticleGoogle Scholar
- Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.PubMedView ArticleGoogle Scholar
- Scohy S, Gabant P, Szpirer C, Szpirer J. Identification of an enhancer and alternative promoter in the first intron of α-fetoprotein gene. Nucleic Acids Res. 2000;28:3743–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Stickland NC, Handel SE. The numbers and types of muscle fibres in large and small breeds of pigs. J Anat. 1986;147:181–9.PubMedPubMed CentralGoogle Scholar
- Kelley G. Mechanical overload and skeletal muscle fiber hyperplasia: a meta-analysis. J Appl Physiol. 1996;81:1584–8.PubMedGoogle Scholar
- Tajbakhsh S. Skeletal musele stem cells in developmental versus regenerative myogenesis. J Intern Med. 2009;266:372–89.PubMedView ArticleGoogle Scholar
- Yamada M, Tatsumi R, Yamanouchi K, Hosoyama T, Shiratsuchi SI, Sato A, et al. High concentrations of HGF inhibit skeletal muscle satellite cell proliferation in vitro by inducing expression of myostatin: a possible mechanism for reestablishing satellite cell quiescence in vivo. Am J Physiol Cell Physiol. 2010;298:C465–76.PubMedPubMed CentralView ArticleGoogle Scholar
- Picard B, Lefaucheur L, Berri C, Duclos MJ. Muscle fibre ontogenesis in farm animal species. Reprod Nutr Dev. 2002;42:415–32.PubMedView ArticleGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456:470–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim E, Magen A, Ast G. Different levels of alternative splicing among eukaryotes. Nucleic Acids Res. 2007;35:125–31.PubMedPubMed CentralView ArticleGoogle Scholar
- Keren H, Lev-Maor G, Ast G. Alternative splicing and evolution: diversification, exon definition and function. Nat Rev Genet. 2010;11:345–55.PubMedView ArticleGoogle Scholar
- Matsumoto K, Kataoka H, Date K, Nakamura T. Cooperative interaction between α-and β-chains of hepatocyte growth factor on c-Met receptor confers ligand-induced receptor tyrosine phosphorylation and multiple biological responses. J Biol Chem. 1998;273:22913–20.PubMedView ArticleGoogle Scholar