|
|
||||||||
RAPID COMMUNICATIONS |
Growth Biology Laboratory, Animal and Natural Resources Institute, USDA-ARS, Beltsville, MD 20705
2 Corresponding author: monika{at}anri.barc.usda.gov
| ABSTRACT |
|---|
|
|
|---|
) and 2 regulatory (ß and
) subunits. Mutations of the
subunit genes are known to affect AMPK functioning. In this study, we characterized the genomic organization and expression of 3 chicken AMPK
subunit genes (cPRKAG). Alternative splicing of the second exon of the cPRKAG1 gene resulted in 2 transcript variants that code for predicted proteins of 298 and 276 amino acids. Use of an alternate promoter and alternative splicing of the cPRKAG2 gene resulted in 4 transcript variants that code for predicted proteins of 567, 452, 328, and 158 amino acids. Alternative splicing of exon 3 of the cPRKAG3 gene resulted in the production of "long" and "short" transcript variants that code for predicted proteins of 382 and 378 amino acids, respectively. We found evidence for differential expression of individual
subunit gene transcript variants and, in some cases, tissue-specific expression was observed. The cPRKAG subunit genes displayed similar structural features and high sequence homology compared with corresponding mammalian
subunit gene homologues.
Key Words: 5'-AMP-activated protein kinase gamma subunit alternative splicing alternate promoter usage energy homeostasis
| INTRODUCTION |
|---|
|
|
|---|
) and 2 regulatory (ß and
) subunits (Mitchelhill et al., 1994; Stapleton et al., 1994). In mammals, there are 2
subunit isoforms, designated
-1 and
-2 (Stapleton et al., 1996); 2 ß subunit isoforms, ß-1 and ß-2 (Stapleton et al., 1997); and 3
subunit isoforms,
-1,
-2, and
-3 (Cheung et al., 2000). Each of these subunit isoforms is encoded by a separate gene (Kahn et al., 2005).
The AMPK is activated in response to metabolic and nutritional stresses that cause a depletion of cellular ATP and increase the intracellular AMP:ATP ratio (Hardie et al., 1998). Activation of AMPK requires phosphorylation of a critical Thr residue (Thr 172) located in the activation loop of the N-terminal kinase domain of the
catalytic subunit by an upstream protein kinase such as the tumor suppressor LKB1 or Ca2+/calmodulin dependent protein luinase kinase (Hawley et al., 1996, 2003; Hurley et al., 2005). Once activated, AMPK phosphorylates a variety of downstream protein targets that affect carbohydrate, protein, and lipid metabolism (Kahn et al., 2005).
The
subunit constitutes the energy-sensing module within the AMPK complex by binding AMP or ATP in a mutually exclusive manner (Scott et al., 2004). Gamma subunit proteins contain 4 repeats of a structural motif consisting of approximately 60 amino acids and classified as a cystathionine ß-synthase (CBS) domain (Bateman, 1997). Two CBS domains form a functional unit recently designated as a "Bateman domain" (Kemp, 2004), which serves as a cooperative regulatory AMP- and ATP-binding site within the AMPK complex (Scott et al., 2004). Several mutations have been identified in human
-2 and pig
-3 subunit genes that occur within or close to a CBS domain. Six different point mutations in the human
-2 gene are associated with certain disease states (e.g., Wolf-Parkinson-White syndrome) characterized by hypertrophic cardiomyopathy, abnormalities in electrical conductance, glycogen overload, and heart failure (Daniel and Carling, 2002). A nonconservative missense mutation (R225Q) in the porcine
-3 gene causes a dominant phenotype characterized by abnormally high skeletal muscle glycogen content and significant effects on meat quality (Milan et al., 2000).
Little is currently known about the structure of chicken AMPK subunit genes. Recently, we reported evidence for 7 chicken AMPK subunit gene homologues located on different chromosomes (Proszkowiec-Weglarz et al., 2006). However, the complete sequence and structural features of individual AMPK subunit genes in chickens, especially those encoding the
subunit isoforms, remain largely unknown. Because AMPK plays a central role in the regulation of energy balance, AMPK subunit genes are important potential candidates to be used as selection markers for feed efficiency and meat quality traits in commercial livestock populations (Benkel et al., 2005). Knowledge about
subunit gene structure and expression is also essential to understanding how AMPK influences physiological function. Therefore, the purpose of this work was to clone, sequence, and characterize transcripts from the 3 chicken AMPK
subunit genes (cPRKAG) and to study their expression in different tissues.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Nucleotide Sequencing
Total RNA from liver, heart, and skeletal muscle and a primer-directed reverse transcription PCR (RT-PCR)-based strategy were used to clone portions of the chicken AMPK
-1,
-2, and
-3 genes (cPRKAG1, cPRKAG2, and cPRKAG3, respectively). Nucleotide sequences for cPRKAG1, cPRKAG2, and cPRKAG3 cDNA, including the complete coding region and portions of the 5'- and 3'-untranslated regions (UTR), were derived using primer sets initially based on predicted genomic sequences, on available EST sequences, or both. With these primers, a series of overlapping PCR products was generated. The overlapping PCR products for each gene were then assembled into a single fragment of contiguous sequence. The presence and identity of unique cPRKAG gene transcript variants were further verified by RT-PCR using a unique 5'-end primer (when possible) and a common downstream 3'-end primer to amplify the entire coding region and portions of the 5'- and 3'-UTR to confirm that all cloned transcript variants corresponded to actual mRNA transcripts. Each PCR product was subjected to bidirectional automated fluorescent DNA sequencing using a Beckman Coulter CEQ 8000XL Genetic Analysis System with the dye terminator cycle sequencing method (Quick Start Kit, Beckman Coulter Inc., Fullerton, CA).
Rapid Amplification of cDNA Ends
Rapid amplification of cDNA ends (RACE) was successfully used to characterize the 5'-ends of cPRKAG2 and cPRKAG3 cDNA. Total RNA (1.0 µg) was used to prepare 5'-RACE-ready cDNA using the SMART RACE cDNA Amplification Kit (BD Biosciences Clontech, Palo Alto, CA). Polymerase chain reaction was performed using Platinum Taq DNA polymerase with 3.5 mM Mg+2 (Invitrogen), touchdown PCR, and cPRKAG2 or cPRKAG3 gene-specific 5'-RACE primers based on sequences obtained by prior sequencing of portions of the coding regions. We attempted to use 5'-RACE to determine 5'-end sequence for cPRKAG1, but we were unsuccessful due to high GC content of this region. In addition, 3'-RACE was not successful in generating 3'-end products for any of the cPRKAG transcripts. Instead, we utilized existing EST and genomic sequence in conjunction with a primer-walking strategy to sequence and verify these regions.
Cloning
Amplified cDNA fragments and RACE products containing sequences corresponding to the 5'-end were either subjected to direct sequencing or were subcloned into the pCR2.1-TOPO vector using the TOPO TA Cloning Kit (Invitrogen) and sequenced using M13 forward and reverse primers.
RT-PCR
Reverse transcription reactions (20 µL) consisted of the following: 1.0 µg of total RNA, 50 units of SuperScript III Reverse Transcriptase (Invitrogen), 40 units of an RNase inhibitor (Invitrogen), 0.5 mM dNTPs, and 100 ng of random hexamer primers. Polymerase chain reaction was performed in 25-µL reactions containing the following: 20 mM Tris-HCl, pH 8.4; 50 mM KCl; 1.0 unit of Platinum Taq DNA Polymerase (Invitrogen); 0.2 mM dNTPs; 2.0 mM Mg2+; 10 pmol of each gene-specific primer (Table 1
); 5 pmol each of an appropriate mixture of primers:competimers specific for 18S rRNA (QuantumRNA Universal 18S Standards kit, Ambion Inc., Austin, TX); and 1 µL of the reverse transcription reaction. Thermal cycling parameters were as follows: 1 cycle at 94°C for 2 min, followed by 30 cycles of 94°C for 30 s, 60°C for 30 s, 72°C for 1 min, with a final extension at 72°C for 8 min. For cPRKAG1 (transcript variants 1 and 2 combined), touch-down PCR was used as follows: 1 cycle at 94°C for 2 min, followed by 15 cycles, 94°C for 30 s, 72°C for 30 s, 72°C for 1.5 min (after each cycle, the temperature was decreased 0.5°C from 72°C to 65°C), followed by 25 cycles, 94°C for 30 s, 65°C for 30 s, 72°C for 1.5 min, with a final extension at 72°C for 8 min. The cPRKAG subunit gene transcript variants were coamplified along with 18S rRNA in a multiplex PCR format for relative quantitative assays. Negative controls were run to ensure PCR accuracy and specificity.
|
Statistical Analysis
Gene expression data were subjected to ANOVA using the GLM procedure of SAS software (SAS System for Windows, Version 8.2, SAS Institute Inc., Cary, NC). The Duncans multiple range test option of the GLM procedure for SAS was used to determine significance of mean differences. Statistical significance was set at P < 0.05.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
-1 subunit protein.
|
|
|
|
-2 subunit proteins ranged from 51 to 95% amino acid identity, whereas the identity specifically within the CBS domain-containing portion of the subunits ranged from 84 to 100%. The predicted proteins from transcript variants 1 to 3 each contain 4 consecutive CBS domains. Use of an internal translation initiation codon (ATGT4 in exon 11) in transcript variant 4 would produce a predicted protein (158 amino acids) that contains only the last 2 CBS domains due to truncation at its N-terminal end. However, it must be noted that transcript 4 possesses 2 potential translation initiation sites (ATGT3,T4). The use of the upstream translation initiation codon (ATGT3 located in exon 6) would produce a more severely truncated
-2 subunit protein (78 amino acids) due to the introduction of a premature termination codon. This truncated protein is predicted to contain some N-terminal sequence and only a portion of the first CBS domain. It is not possible to determine whether the upstream, the downstream, or both translation initiation codons would be utilized for transcript variant 4. In eukaryotic species, a scanning mechanism that recognizes sequence flanking the AUG codon determines the site(s) of translation initiation, and this process does permit the production of multiple proteins from a single mRNA (Kozak, 2001). The production of 2 different size peptides from a single mRNA by 2 initiation sites has recently been reported for the human neuropeptide Y gene (Kaipio et al., 2005). Moreover, the use of alternative initiation of translation in humans has been reported to occur most frequently in gene transcripts coding for regulatory proteins (Prats and Prats, 2002). However, it is not known to what extent this applies to cPRKAG2 transcript variant 4, because it is generally not highly expressed (see below) and both predicted proteins are severely truncated and presumably functionally impaired.
Individual cPRKAG2 gene transcript variants were expressed in all tissues, and some exhibited a tissue-specific expression pattern (Figure 5
, panels A to F). Transcript variant 1 was preferentially expressed in liver (Figure 5
, panel C). The highest level of expression of transcript variant 2 was observed in liver and heart tissue (Figure 5
, panel D). Brain, heart, kidney, skeletal muscle, abdominal fat, and hypothalamus preferentially expressed transcript variant 3 (Figure 5
, panel E), whereas, in spleen and duodenum, transcript variants 1 to 3 were expressed at similar levels (Figure 5
, panels B to E). The lowest expression of all 4 transcripts for cPRKAG2 was observed in pancreas (Figure 5
, panels B to F). Transcript variant 4 was most highly expressed in skeletal muscle tissue (Figure 5
, panel F), albeit at a reduced level compared with transcript variants 1 to 3 (Figure 5
, panel B).
|
subunit genes. It consists of 12 coding exons and 1 noncoding exon, ranging in size from 38 to >577 bp, and 12 introns, ranging in size from 73 to 759 bp, with the first intron being the largest (Figure 6
-3 subunit protein as compared with mammalian proteins. The core group of 8 exons (4 to 11) that encode a significant portion of the CBS domain region of the
-3 isoforms is totally conserved with respect to size and order when compared with mammalian PRKAG3 counterpart genes, as well as with cPRKAG1 and cPRKAG2 genes. Amarger et al. (2003) compared the genomic organization of PRKAG3 genes of the pig, human, and mouse and observed that exon size was perfectly conserved, with the exception of the first exon. Similar findings were reported for the equine and zebrafish PRKAG3 genes (Park et al., 2003). The cPRKAG3 gene, located on chromosome 7, spans 4.1 kb, making it much more compact than cPRKAG2. This organization is similar to the human PRKAG3 gene, which is located on chromosome 2 (2q35). The cPRKAG3 gene exhibits syntenic conservation as compared with its human counterpart, based on a consensus linkage map of the chicken genome (Groenen et al., 2000). Park et al. (2003) reported the presence of a conserved block of synteny involving the linkage of PRKAG3 with the KIAA0173 and serine/ threonine kinase genes in mammals and 2 fish species. We have found the same linkages to occur on chromosome 7 in the chicken as well.
|
-3 subunit proteins ranged from 65 to 67%, whereas the sequence identity specifically within the CBS domain-containing portion of the subunits ranged from 68 to 88%.
The cPRKAG3 gene is expressed exclusively in heart and skeletal muscle. Both of these tissues demonstrated a predominance of transcript variant 2 ("short" form), which is expressed at a significantly higher level (1.5-fold) compared with transcript variant 1 ("long" form). Moreover, the expression of both transcripts was higher (2-fold) in skeletal muscle as compared with heart tissue (Figure 7
). In contrast, significant expression of PRKAG3 in mammals has only been observed in skeletal muscle, but not in heart (Cheung et al., 2000). Moreover, Roux et al. (2006) observed preferential expression of the "short" transcript of the bovine PRKAG3 gene.
|
subunit genes in the chicken. The combined effects of alternative splicing and alternate promoter usage resulted in the production of multiple transcripts for each of the cPRKAG genes and, in some cases, tissue-specific expression indicative of additional levels of complexity in the regulation of these genes. Lareau et al. (2004) have suggested that alternative splicing, common in many eukaryotes, may play an important role in regulating protein expression, thus contributing substantially to biological complexity. Although multiple cPRKAG subunit gene transcripts are expressed, it remains to be shown that the proteins encoded by them are actually produced and are indeed functional, especially those that contain missing or truncated CBS domains. This will require further study with the production and use of specific antibody reagents that recognize each of the
subunit proteins. Each of the cPRKAG genes showed a high degree of similarity to mammalian counterparts, indicative of a conserved function for the
subunit proteins and possibly reflecting their derivation by duplication of a common ancestral gene, as has been suggested by Park et al. (2003) for mammalian PRKAG genes. Our findings offer new insights into chicken AMPK subunit gene structure and expression.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Received for publication April 19, 2006. Accepted for publication June 12, 2006.
| REFERENCES |
|---|
|
|
|---|
Bateman, A. 1997. The structure of a domain common to archae-bacteria and the homocystinuria disease protein. Trends Biochem. Sci. 22:1213.[ISI][Medline]
Benkel, B., S. Kollers, R. Fries, A. Sazanov, E. Yoshida, E. Valle, J. Davoren, and D. Hickey. 2005. Characterization of the bovine ampk
1 gene. Mamm. Genome 16:194200.[ISI][Medline]
Carling, D. 2005. AMP-activated protein kinase: Balancing the scales. Biochimie 87:8791.[Medline]
Cheung, P. C., I. P. Salt, S. P. Davies, D. G. Hardie, and D. Carling. 2000. Characterization of AMP-activated protein kinase
-subunit isoforms and their role in AMP binding. Biochem. J. 346:659669.
Daniel, T., and D. Carling. 2002. Functional analysis of mutations in the
2 subunit of AMP-activated protein kinase associated with cardiac hypertrophy and Wolff-Parkinson-White syndrome. J. Biol. Chem. 277:5101751024.
Groenen, M. A. M., H. H. Cheng, N. Bumstead, B. F. Benkel, W. E. Briles, T. Burke, D. W. Burt, L. B. Crittenden, J. Dodgson, J. Hillel, S. Lamont, A. Ponce de Leon, M. Soller, H. Takahashi, and A. Vignal. 2000. A consensus linkage map of the chicken genome. Genome Res. 10:137147.
Hardie, D. G., D. Carling, and M. Carlson. 1998. The AMP-activated/SNF1 protein kinase subfamily: Metabolic sensors of the eukaryotic cell? Annu. Rev. Biochem. 67:821855.[ISI][Medline]
Hawley, S. A., J. Boudeau, J. L. Reid, K. J. Mustard, L. Udd, T. P. Makela, D. R. Alessi, and D. G. Hardie. 2003. Complexes between the LKB1 tumor suppressor, STRAD
/ß and MO25
/ß are upstream kinases in the AMP-activated protein kinase cascade. J. Biol. 2:28.[Medline]
Hawley, S. A., M. Davison, A. Woods, S. P. Davies, R. K. Beri, D. Carling, and D. G. Hardie. 1996. Characterization of the AMP-activated protein kinase from rat liver and identification of threonine-172 as the major site at which it phosphorylates AMP-activated protein kinase. J. Biol. Chem. 271:2787927887.
Hurley, R. L., K. A. Anderson, J. M. Franzone, B. E. Kemp, A. R. Means, and L. A. Witters. 2005. The Ca2+/calmodulin-dependent protein kinase kinases are AMP-activated protein kinase kinases. J. Biol. Chem. 280:2906029066.
Kahn, B. B., T. Alquier, D. Carling, and D. G. Hardie. 2005. AMP-activated protein kinase: Ancient energy gauge provides clues to modern understanding of metabolism. Cell Metab. 1:1525.[ISI][Medline]
Kaipio, K., J. Kallio, and U. Pesonen. 2005. Mitochondrial targeting signal in human neuropeptide Y gene. Biochem. Biophys. Res. Commun. 337:633640.[ISI][Medline]
Kemp, B. E. 2004. Bateman domains and adenosine derivatives form a binding contract. J. Clin. Invest. 113:182184.[ISI][Medline]
Kozak, M. 2001. Constraints on reinitiation of translation in mammals. Nucleic Acids Res. 29:52265232.
Lang, T., L. Yu, Q. Tu, J. Jiang, Z. Chen, Y. Xin, G. Liu, and S. Zhao. 2000. Molecular cloning, genomic organization, and mapping of PRKAG2, a heart abundant
2 subunit of 5'-AMP-activated protein kinase, to human chromosome 7q36. Genomics 70:258263.[ISI][Medline]
Lareau, L. F., R. E. Green, R. S. Bhatnagar, and S. E. Brenner. 2004. The evolving roles of alternative splicing. Curr. Opin. Struct. Biol. 14:273282.[ISI][Medline]
Milan, D., J.-T. Jeon, C. Looft, V. Amarger, A. Robic, M. The-lander, C. Rogel-Gaillard, S. Paul, N. Iannuccelli, L. Rask, H. Ronne, K. Lundström, N. Reinsch, J. Gellin, E. Kalm, P. Le Roy, P. Chardon, and L. Andersson. 2000. A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle. Science 288:12481251.
Mitchelhill, K. I., D. Stapleton, G. Gao, C. House, B. Michell, F. Katsis, L. A. Witters, and B. E. Kemp. 1994. Mammalian AMP-activated protein kinase shares structural and functional homology with the catalytic domain of yeast Snf1 protein kinase. J. Biol. Chem. 269:23612364.
Park, H. B., S. Marklund, J. T. Jeon, J. R. Mickelson, S. J. Valberg, K. Sandberg, and L. Andersson. 2003. Molecular characterization and mutational screening of the PRKAG3 gene in the horse. Cytogenet. Genome Res. 102:211216.[ISI][Medline]
Prats, A. C., and H. Prats. 2002. Translational control of gene expression: Role of IRESs and consequences for cell transformation and angiogenesis. Prog. Nucleic Acid Res. Mol. Biol. 72:367413.[ISI][Medline]
Proszkowiec-Weglarz, M., M. P. Richards, R. Ramachandran, and J. P. McMurtry. 2006. Characterization of the AMP-activated protein kinase pathway in chickens. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 143:92106.[Medline]
Richards, M. P., and S. M. Poch. 2002. Quantitative analysis of gene expression by reverse transcription polymerase chain reaction and capillary electrophoresis with laser-induced fluorescence detection. Mol. Biotechnol. 21:1937.[ISI][Medline]
Roux, M., A. Nizou, L. Forestier, A. Ouali, H. Leveziel, and V. Amarger. 2006. Characterization of the bovine PRKAG3 gene: Structure, polymorphism, and alternative transcripts. Mamm. Genome 17:8392.[ISI][Medline]
Scott, J. W., S. A. Hawley, K. A. Green, M. Anis, G. Stewart, G. A. Scullion, D. G. Norman, and D. G. Hardie. 2004. CBS domains form energy-sensing modules whose binding of adenosine ligands is disrupted by disease mutations. J. Clin. Invest. 113:274284.[ISI][Medline]
Shamsadin, R., K. Jantsan, I. Adham, and W. Engel. 2001. Cloning, organisation, chromosomal localization and expression analysis of the mouse PRKAG1 gene. Cytogenet. Cell Genet. 92:134138.[ISI][Medline]
Stapleton, D., G. Gao, B. J. Michell, J. Widmer, K. Mitchelhill, T. The, C. M. House, L. A. Witters, and B. E. Kemp. 1994. Mammalian 5'-AMP-activated protein kinase non-catalytic subunits are homologues of proteins that interact with yeast Snf1 protein kinase. J. Biol. Chem. 269:2934329346.
Stapleton, D., K. I. Mitchelhill, G. Gao, J. Widmer, B. J. Michell, T. The, C. M. House, C. S. Fernandez, T. Cox, L. A. Witters, and B. E. Kemp. 1996. Mammalian AMP-activated protein kinase subfamily. J. Biol. Chem. 271:611614.
Stapleton, D., E. Woollatt, K. I. Mitchelhill, J. K. Nicholl, C. S. Fernandez, B. J. Michell, L. A. Witters, D. A. Power, G. R. Sutherland, and B. E. Kemp. 1997. AMP-activated protein kinase isoenzyme family: Subunit structure and chromosomal location. FEBS Lett. 409:452456.[ISI][Medline]
This article has been cited by other articles:
![]() |
M. P. Richards and M. Proszkowiec-Weglarz Mechanisms Regulating Feed Intake, Energy Expenditure, and Body Weight in Poultry Poult. Sci., July 1, 2007; 86(7): 1478 - 1490. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |