Molecular analysis of alternative transcripts of the equine cordon-bleu WH2 repeat protein-like 1 (COBLL1) gene.
The Thoroughbred is the most famous breed for horse racing because of its speed and agility. The Thoroughbred industry as a large agribusiness creates annual revenue of around $34 billion in the United States alone, and provides about 470,000 jobs related to the network of farms, training centers and race tracks. Genes related to their athletic phenotypes have not been identified, although physical and physiological adaptations related to the potentials of their superior athleticism are well characterized (Gu et al., 2009).
In a previous study, the whole transcriptome for blood and muscle tissues of six Thoroughbred horses before and after 30 min of exercises were analyzed using RNA-sequencing. From the analysis, 32,361 unigene clusters were identified, of which 20,428 have been newly reported as belonging to the equine gene model. Differentially expressed genes were also found in blood (142) and muscle (1,163). In addition, alternative splicing forms of cytoplasmic dynein 1 (DYNCl), AXL receptor tyrosine kinase (AXL), pleckstrin homology domain containing, family G (with RhoGef Domain) member 1 (PLEKHG1), and the cordon-bleu WH2 repeat protein-like 1 (COBLL1) showed reversed expression patterns before and after exercise in skeletal muscle (Park et al., 2012).
The COBLL1 gene was originally cloned from the human brain cDNA library and is highly expressed in most human tissues including lung, liver, kidney, pancreas, ovary, spinal cord and brain (Nagase et al., 1999). COBLL1 was identified as another class of Cobl (Carroll et al., 2003). COBLL1 has been reported as a negative regulator of apoptosis in tumor cells (Gordon et al., 2011) and has been associated with metabolic diseases related to insulin resistance (Mancina et al., 2013). Despite various studies on COBLL1, the specific roles of both genes are still unclear, especially in horses. In this study, we focused on the analysis of the COBLL1 gene which was one of the alternatively spliced genes during exercise. We confirmed alternative splicing forms of COBLL1 in various horse tissues and examined the expression pattern of each alternative splicing form in response to exercise.
MATERIALS AND METHODS
Tissue sampling in horse
Three male Thoroughbred horses aged 5, 9, and 10 years, which were maintained at Ham-an Racing Horse Resort and Training Center, were used to obtain the blood and skeletal muscle samples before and after exercise. Exercise was performed by trotting at the speed of 13 km/h for 30 min. The National Institute of Subtropical Agriculture, Rural Development Administration, provided three Jeju horses which were used for collecting the tissue samples such as skeletal muscle, kidney, thyroid, lung, appendix, colon, spinal cord and heart. All tissue samples were kept in liquid nitrogen tank until extraction of RNA was done. All procedures were conducted by following the protocol that had been reviewed and approved by the Institutional Animal Care and Use Committee at Pusan National University (protocol numbers: PNU-2013-0417, PNU-2013-0411).
RNA extraction and cDNA synthesis
Trizol reagent (Invitrogen, Carlsbad, CA, USA) was used to extract total RNA from tissue samples (skeletal muscle, kidney, thyroid, lung, colon, spinal cord, and heart) and leukocytes after exercise in horse, according to the Invitrogen manual. In order to prevent contamination of genomic DNA, RNase-free DNase kit (Qiagen, Venlo, Netherlands) was used according to the manufacturer's operating manual. Total RNA quantification was performed by using a NanoDrop ND-1000 Spectrophotometer. The cDNAs were synthesized in a reaction with oligo-dT primers, moloney-murine leukemia virus reverse transcriptase (Promega, Madison, WI, USA), RNase inhibitor (Promega, USA) and RNase-free dd[H.sub.2]O, which were incubated at 37[degrees]C for 4 h.
Reverse transcription polymerase chain reaction amplification and sequencing
The horse COBLL1 transcripts were analyzed by reverse transcription polymerase chain reaction (RT-PCR) amplification. The original horse COBLL1 transcripts were amplified by the primer pairs S (5'-ATTCCCTGGAGTCAGTCCATG CAG-3') and AS (5'-GGGAAGTGCCTAAAGATGAGG CTG-3'). These primers were designed using the PRIMER3 program (http://frodo.wi.mit.edu/primer3/) and sequence information was obtained from the UCSC genome browser (Accession number: NM-001081812.1). The expected sizes of PCR products were 276 bp and 156 bp for COBLL1a and COBLL1b, respectively. The RT-PCR conditions were as follows: an initial step of 94[degrees]C for 4 min, 35 cycles of 94[degrees]C for 40 s, 56[degrees]C for 1 min, and 72[degrees]C for 2 min were performed. The RT-PCR products were analyzed after electrophoresis on a 1.5% agarose gel. COBLL1a and COBLL1b were cloned into pGEM-T Easy-cloning vector (Promega, USA) and their nucleotide sequences were determined by Sanger sequencing. The sequenced nucleotides were aligned by Blast (National Center for Biotechnology Information; NCBI, Bethesda, MD, USA).
Real-time quantitative polymerase chain reaction analysis
To analyze the expression level of COBLL1 alternative splicing isoforms in muscle before and after exercise, real-time quantitative polymerase chain reaction (RT-qPCR) was conducted by using the BioRad CFX-96 machine (BioRad, Hercules, CA, USA). Primer sequences were designed to detect each type of alternative splicing forms. Two primer pairs, Cobll1a-S (5'-ATCAGTAGACCCGCACGGAA-3') and Cobll1a-AS (5'-ACACTGACGGAGGAAACCTC-3') for COBLL1a, and Cobll1b-S (5'-AGGGTCTGTCCAGCC CAGCT-3') and Cobll1b-AS (5'-CCATAACGTGTGGTC CCTGTCC-3') for COBLL1b, were used to specifically detect each isoform. Each reaction was executed in a total 25 [micro]L of mixture containing 14 [micro]L of SYBR green master mix, 2 [micro]L of forward primer (5 pmol), 2 [micro]L of reverse primer (5 pmol), 5 [micro]L of distilled water, and 2 [micro]L (50 ng/[micro]L) of cDNA. The PCR conditions were at 94[degrees]C for 5 min of pre-denaturation step, 39 cycles of 94[degrees]C for 20 s, 56[degrees]C for 20 s and 72[degrees]C for 30 s, and followed by 72[degrees]C for 10 min as a final step. All samples were measured in triplicate to ensure reproducibility, and [C.sub.t] value was used to calculate the fold change by using the [2.sup.-[DELTA][DELTA]Ct] method (Livak and Schmittgen, 2001). The glyceraldehyde 3-phosphate dehydrogenase gene was used for reference.
Phylogenetic analysis of equine COBLL1
The amino acid sequences of COBLL1 of various species were retrieved from the Gene database of NCBI (wild horse, XP_008531702.1; human, NP_055715.3; chimpanzee, XP_003309326.1; mouse, NP_795999.2; rat, XP_006234363.1; cow, XP_002685420.3; dog, XP_005640312.1; chicken, XP_004942872.1; frog, XP_002939319.2). The amino acid sequences were aligned and used for generating phylogenetic tree by using the MUSCLE and Neighbor-Joining methods, respectively, in MEGA6 program (Tamura et al., 2013). The web-based protein prediction program, PredictProtein (www.predictprotein.org), was used to predict the secondary and tertiary structure (Rost et al., 2004).
RESULTS AND DISCUSSION
Comparison of COBLL1 amino acid sequences among various species and structural analysis of equine COBLL1 gene according to alternative splicing
The horse COBLL1 gene has 14 exons and 13 introns spanning about 60.5 kb on the chromosome 18 and six transcriptional variants had been predicted in NCBI Gene database. When the amino acid sequences of COBLL1 were compared among various species, COBLL1 showed higher identity at the N-terminal sequences containing a Cordonbleu ubiquitin-like domain than at the C-terminal (Figure 1A, solid box). Interestingly, one of the alternative splicing forms of COBLL1 gene had been observed within the Cordon-bleu ubiquitin-like domain of various species including horse, wild horse, human, dog, chimpanzee, and chicken. However, this variant did not seem to affect the function of the domain as it was predicted as the Cordonbleu ubiquitin-like domain. Another alternative splicing was observed between exon 8 and 10 (with skipping of exon 9) in many species including horse, wild horse, human, mouse, rat, dog, and chicken (Figure 1A, dashed box). The similarity of amino acid sequences around these sequences was relatively lower than the other N-terminal region shown in Figure 1A. Horse COBLL1 was evolutionary close to cow COBLL1 while largely divergent from frog (Figure 1B). COBLL1 has been known to regulate apoptosis or to be associated with metabolic diseases related to insulin resistance. It has been reported that alternative splicing forms of a gene may be differentially expressed under conditions of wound response, tumorogenesis, embryonic development and response to external stimuli (Vitale et al., 1994; Shin and Manley, 2004; Huang et al., 2005; Gardina et al., 2006; Lynch, 2007). Therefore, we assumed that COBLL1 was evolutionally conserved through vertebrates, and the expression patterns of the transcript variants were different according to conditions during stimulation.
Validation and expression pattern of splicing variants of equine COBLL1
In a previous study, two alternative splicing forms of COBLL1 transcripts were identified in horses (Park et al., 2012). And they are found in both NCBI (XM_005601553.1 and XM_005601555.1) and Ensembl (ENSECAT00000027104 and ENSECAT00000027105) data base and designated to COBLL1a and COBLL1b, respectively (Figure 2A). Specific primers were designed for binding to exon 7 and 8 forward primers and exon 10 reverse primer to confirm the expression of COBLL1a and COBLL1b in horse. When PCR for horse COBLL1 transcripts in blood was performed, two transcriptional variants with different sizes (277 bp and 157 bp) were detected (Figure 2B). Each amplicon for COBLL1a and COBLL1b transcript was detected with 120 bp difference in size by gel electrophoresis. By sequencing analysis with the amplicons, the deletion of exon 9 was observed in COBLL1b as expected (Figure 2C).
Expression patterns of COBLL1a and COBLL1b in various tissues and response to exercise
To investigate the expression patterns of COBLL1 variants in horse tissues, RT-PCR was conducted with cDNAs from the various tissues including kidney, spinal cord, heart, thyroid, skeletal muscle, colon and lung (Figure 3A). The COBLL1a showed strong expression in the kidney, spinal cord and lung, moderate expression in the heart and skeletal muscle, and low expression in the thyroid and colon. Unlike COBLL1a expression, COBLL1b was expressed evenly in all tissues. The expression patterns analysis of COBLL1a and COBLL1b in muscle tissues before and after exercise revealed that the expression of both COBLL1a and COBLL1b decreased after exercise, but only COBLL1b showed significant decrease (Figure 3B). Thus, it is speculated that COBLL1a and COBLL1b responded similarly to the effects of exercise in skeletal muscle.
According to a recent genome-wide association study, the human COBLL1 gene was found to be located at the genomic region associated with a metabolic traits such as fasting insulin level, type 2 diabetes, and triacylglycerol metabolism (Manning et al., 2012; Albrechtsen et al., 2013; Desmarchelier et al., 2014). In addition, the C allele of COBLL1 is functionally linked to lower serum insulin levels and insulin resistance in humans (Mancina et al., 2013). These two reports suggest that COBLL1 may play an important role in metabolic pathways, especially in the regulation of glucose uptake by insulin. Generally, exercise increases the uptake of glucose mediated by up-regulation of GLUT4 protein content and leads to enhanced insulin resistance in muscle (Goodyear and Kahn, 1998; Borghouts and Keizer, 2000). Given the importance of glucose metabolism in exercise, it is of interest that the regulation of COBLL1 expression maybe important for regulating glucose level or switching of energy source, possibly through an insulin signaling pathway, in muscle after exercise. Taken together, this information may be helpful to investigate the role of equine COBLL1 on the metabolic changes related to the athletic performance of racehorses and to develop COBLL1 as a useful biomarker for athletic performance.
Protein prediction of COBLL1a and COBLL1b
In analyzing the amino acid sequences of COBLL1 splicing forms, COBLL1b has a deletion of 40 amino acids due to the exclusion of exon 9 by alternative splicing compared to COBLL1a (Figure 4A). These sequences were predicted to include three putative phosphorylation sites by casein kinase 1, p38 mitogen-activated protein kinase, and protein kinase C (Figure 4B). In many cases, protein phosphorylation is closely related to the control of protein stabilization, serving as a marker that triggers subsequent ubiquitination, in particular where ubiquitination leads to degradation (Treier et al., 1994; Fuchs et al., 1996; Magnani et al., 2000). Therefore, COBLL1a could be elaborately controlled in each tissue, whereas COBLL1b could be stable in various tissues due to lack of these phosphorylation sites by alternative splicing (Figure 3A).
This work was supported by a grant from the Next-Generation BioGreen 21 Program (No. PJ011173, PJ011044), Rural Development Administration, Republic of Korea. Authors thank Ms S.H. Moon, MSc, for editorial help.
Albrechtsen, A., N. Grarup, Y Li, T. Sparso, G Tian, H. Cao, T. Jiang, S. Y Kim, T. Korneliussen, and Q. Li et al. 2013. Exome sequencing-driven discovery of coding polymorphisms associated with common metabolic phenotypes. Diabetologia 56:298-310.
Borghouts, L. B. and H. A. Keizer. 2000. Exercise and insulin sensitivity: A review. Int. J. Sports Med. 21:1-12.
Carroll, E. A., D. Gerrelli, S. Gasca, E. Berg, D. R. Beier, A. J. Copp, and J. Klingensmith. 2003. Cordon-bleu is a conserved gene involved in neural tube formation. Dev. Biol. 262:16-31.
Desmarchelier, C., J. C. Martin, R. Planells, M. Gastaldi, M. Nowicki, A. Goncalves, R. Valero, D. Lairon, and P. Borel. 2014. The postprandial chylomicron triacylglycerol response to dietary fat in healthy male adults is significantly explained by a combination of single nucleotide polymorphisms in genes involved in triacylglycerol metabolism. J. Clin. Endocrinol. Metab. 99:E484-488.
Fuchs, S. Y, L. Dolan, R. J. Davis, and Z. Ronai. 1996. Phosphorylation-dependent targeting of c-Jun ubiquitination by Jun N-kinase. Oncogene 13:1531-1535.
Gardina, P. J., T. A. Clark, B. Shimada, M. K. Staples, Q. Yang, J. Veitch, A. Schweitzer, T. Awad, C. Sugnet, S. Dee, C. Davies, A. Williams, and Y Turpaz. 2006. Alternative splicing and differential gene expression in colon cancer detected by a whole genome exon array. BMC Genomics 7:325.
Goodyear, L. J. and B. B. Kahn. 1998. Exercise, glucose transport, and insulin sensitivity. Annu. Rev. Med. 49:235-261.
Gordon, G J., R. Bueno, and D. J. Sugarbaker. 2011. Genes associated with prognosis after surgery for malignant pleural mesothelioma promote tumor cell survival in vitro. BMC Cancer 11:169.
Gu, J., N. Orr, S. D. Park, L. M. Katz, G Sulimova, D. E. MacHugh, and E. W. Hill. 2009. A genome scan for positive selection in thoroughbred horses. PLoS One 4(6):e5767.
Huang, X., Y. Guo, Y. Shui, S. Gao, H. Yu, H. Cheng, and R. Zhou. 2005. Multiple alternative splicing and differential expression of dmrtl during gonad transformation of the rice field eel. Biol. Reprod. 73:1017-1024.
Livak, K. J. and T. D. Schmittgen. 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25:402-408.
Lynch, K. W. 2007. Regulation of alternative splicing by signal transduction pathways. Adv. Exp. Med. Biol. 623:161-174. Magnani, M., R. Crinelli, M. Bianchi, and A. Antonelli. 2000. The ubiquitin-dependent proteolytic system and other potential targets for the modulation of nuclear factor-kB (NF-kB). Curr. Drug Targets 1:387-399.
Mancina, R. M., M. A. Burza, C. Maglio, C. Pirazzi, F. Sentinelli, M. Incani, T. Montalcini, A. Pujia, T. Congiu, S. Loche, S. Pilia, O. Wiklund, J. Boren, S. Romeo, and M. G. Baroni. 2013. The COBLL1 C allele is associated with lower serum insulin levels and lower insulin resistance in overweight and obese children. Diabetes Metab. Res. Rev. 29:413-416.
Manning, A. K., M. F. Hivert, R. A. Scott, J. L. Grimsby, N. Bouatia-Naji, H. Chen, D. Rybin, C. T. Liu, L. F. Bielak, and I. Prokopenko et al. 2012. A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance. Nat. Genet. 44:659-669.
Nagase, T., K. Ishikawa, M. Suyama, R. Kikuno, M. Hirosawa, N. Miyajima, A. Tanaka, H. Kotani, N. Nomura, and O. Ohara. 1999. Prediction of the coding sequences of unidentified human genes. XIII. The complete sequences of 100 new cDNA clones from brain which code for large proteins in vitro. DNA Res. 6:63-70.
Park, K. D., J. Park, J. Ko, B. C. Kim, H. S. Kim, K. Ahn, K. T. Do, H. Choi, H. M. Kim, S. Song, S. Lee, S. Jho, H. S. Kong, Y. M. Yang, B. H. Jhun, C. Kim, T. H. Kim, S. Hwang, J. Bhak, H. K. Lee, and B. W. Cho. 2012. Whole transcriptome analyses of six thoroughbred horses before and after exercise using RNA-Seq. BMC Genomics 13:473.
Rost, B., G. Yachdav, and J. Liu. 2004. The PredictProtein server. Nucl. Acids Res. 32:W321-326.
Shin, C. and J. L. Manley. 2004. Cell signalling and the control of pre-mRNA splicing. Nat. Rev. Mol. Cell Biol 5:727-738.
Tamura, K., G. Stecher, D. Peterson, A. Filipski, and S. Kumar. 2013. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30:2725-2729.
Treier, M., L. M. Staszewski, and D. Bohmann. 1994. Ubiquitin-dependent c-Jun degradation in vivo is mediated by the delta domain. Cell 78:787-798.
Vitale, A. T., M. Pedroza-Seres, V. Arrunategui-Correa, S. J. Lee, S. DiMeo, C. S. Foster, and R. B. Colvin. 1994. Differential expression of alternatively spliced fibronectin in normal and wounded rat corneal stroma versus epithelium. Invest. Ophthalmol. Vis. Sci. 35:3664-3672.
Submitted Sept. 16, 2014; Revised Nov. 13, 2014; Accepted Jan. 26, 2015
Jeong-Woong Park (1, a), Hyun-Jun Jang (2, 3), Sangsu Shin (1, 3), Hyun-Woo Cho (1), Jae-Young Choi (1), Nam-Young Kim (4), Hak-Kyo Lee (5), Kyong-Tak Do (6), Ki-Duk Song (5), *, and Byung-Wook Cho (1, 3), *
(1) Department of Animal Science, College of Natural Resources and Life Science, Pusan National University, Miryang 627-706, Korea
* Corresponding Author: Ki-Duk Song. Tel: +82-31-670-5656, Fax: +82-31-670-5661, E-mail: firstname.lastname@example.org / Byung-Wook Cho. Tel: +82-55-350-5515, Fax: +82-55-581-2962, E-mail: bwcho@pu san.ac.kr
(2) College of Pharmacy, Dankook University, Cheonan 330-714, Korea.
(3) Life and Industry Convergence Research Institute, College of Natural Resources and Life Science, Pusan National University, Miryang 627-706, Korea.
(4) Subtropical Animal Experiment Station, National Institute of Animal Science, Rural Development Administration, Jeju 690-150, Korea.
(5) Genomic Informatics Center, Hankyong National University, Anseong 456-749, Korea.
(6) Department of Equine Sciences, Sorabol College, Gyeongju 780711, Korea.
(a) These authors contributed equally to this work.
|Printer friendly Cite/link Email Feedback|
|Author:||Park, Jeong-Woong; Jang, Hyun-Jun; Shin, Sangsu; Cho, Hyun-Woo; Choi, Jae-Young; Kim, Nam-Young; Lee|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Jun 1, 2015|
|Previous Article:||Estimation of pork quality traits using exsanguination blood and postmortem muscle metabolites.|
|Next Article:||Biocomputational characterization and evolutionary analysis of bubaline Dicer1 enzyme.|