Characterization of hARD2, a processed hARD1 gene duplicate, encoding a human protein N-α-acetyltransferase
© Arnesen et al; licensee BioMed Central Ltd. 2006
Received: 21 November 2005
Accepted: 25 April 2006
Published: 25 April 2006
Protein acetylation is increasingly recognized as an important mechanism regulating a variety of cellular functions. Several human protein acetyltransferases have been characterized, most of them catalyzing ε-acetylation of histones and transcription factors. We recently described the human protein acetyltransferase hARD1 (h uman Ar rest D efective 1). hARD1 interacts with NATH (N-A cetyl T ransferase H uman) forming a complex expressing protein N-terminal α-acetylation activity.
We here describe a human protein, hARD2, with 81 % sequence identity to hARD1. The gene encoding hARD2 most likely originates from a eutherian mammal specific retrotransposition event. hARD2 mRNA and protein are expressed in several human cell lines. Immunoprecipitation experiments show that hARD2 protein potentially interacts with NATH, suggesting that hARD2-NATH complexes may be responsible for protein N-α-acetylation in human cells. In NB4 cells undergoing retinoic acid mediated differentiation, the level of endogenous hARD1 and NATH protein decreases while the level of hARD2 protein is stable.
A human protein N-α-acetyltransferase is herein described. ARD2 potentially complements the functions of ARD1, adding more flexibility and complexity to protein N-α-acetylation in human cells as compared to lower organisms which only have one ARD.
Protein acetylation is a very common modification with a significant impact on several cellular processes. Acetylation occurs both at lysine residues within proteins (Nε-acetylation) and at the N-terminus of proteins (Nα-acetylation). In yeast, N-acetyltransferase 1 (Nat1p) complexes with Arrest defective 1 (Ard1p) to generate a functional NatA protein Nα-acetyltransferase , Ard1p being the catalytic subunit. Proteins with Ser-, Thr-, Gly-, or Ala- N-termini are described to be substrates of NatA after methionine cleavage . The yeast NatB and NatC complexes acetylates different subsets of methionine N-termini [2–4]. Almost all known N-terminally acetylated yeast proteins are products of one of these Nat complexes. Protein N-terminal acetylation is generally believed to be a cotranslational process linked to the ribosome [6–10]. hARD1, the human protein with highest sequence similarity to yeast ARD1, has been described on the genomic (TE2, GenBank [NM_003491]) , mRNA , protein, and enzyme activity levels . Endogenous hARD1 was demonstrated to interact with NATH and express protein Nα-acetyltransferase activity. The complex was found to interact with ribosomal subunits supporting its function in cotranslational acetylation . In vitro translated mouse homologues, mNAT1 and mARD1, have also been shown to interact and express N-acetyltransferase activity . In S. cerevisiae and D. melanogaster, a third subunit of the NatA complex has been described and named Nat5p and San, respectively [8, 14]. The function of this subunit is unknown, but sequence analysis suggests that Nat5p/San is an acetyltransferase. The human orthologue, hNAT5, was also recently demonstrated to be a part of the human NatA complex .
Even though 80–90 % of all mammalian proteins and 50 % of yeast proteins are estimated to be cotranslationally Nα-acetylated [4, 16–20], only a few examples exist describing the functional importance of proper Nα-acetylation. For instance, the function of the yeast proteins Orc1p and Sir3p in telomeric silencing is dependent on proper NatA-mediated Nα-acetylation of these proteins [21, 22].
Using yeast null strains, NatA activity has been demonstrated to be associated with Go entry, cell growth, and the ability to sporulate [23–26]. The importance of protein Nα-acetylation has also been described in C. elegans, where knockdown of either the ard1 or nat1 homologues resulted in embryonal lethality . The human NatA complex has also recently been demonstrated to be essential for normal cellular viability. RNA interference mediated knockdown of NATH or hARD1 induced apoptosis in HeLa cells .
Mouse ARD1 was also reported to be implicated in the acetylation of lysine 532 of HIF-1α, contributing to its degradation in normoxia . However, several independent investigations have reported that at least the wildtype hARD1 protein does not mediate Nε-acetylation of the lysine residue 532 of HIF-1α [29–32].
The hARD1 gene is located on chromosome X (Xq28). Database searches revealed the presence on chromosome 4 (4q21.23) of a putative human paralogue of the previously published hARD1 gene (GeneID:84779, hypothetical protein [MGC10646]). We named this hypothetical human ARD, hARD2.
Here we describe the cloning and expression of hARD2. The entire ORF of hARD2 is intronless, resembling a gene duplicate. Many gene duplicates are non-functional pseudogenes but some, including hARD2, are active genes producing mRNAs and proteins [33–35]. Similar to hARD1, hARD2 interacts with NATH and expresses N-α-acetyltransferase activity.
hARD2 cloning and expression
Evolution of the hARD2 gene
Structural comparison between hARD1 and hARD2
Overexpressed hARD2 co-immunoprecipitate NATH and HIF-1α
The interaction between hARD1 and HIF-1α has been described . To investigate whether also hARD2 is capable of interacting with HIF-1α, we cotransfected MCF-7 cells with plasmids encoding Xpress-hARD2 and HA-HIF-1α. HA-HIF-1α will then accumulate under these normoxic conditions probably due to saturation of VHL-mediated degradation of HIF-1α. The resulting cell lysates were subjected to immunoprecipitation with the anti-Xpress antibody. Western blotting analysis using anti-HA demonstrated that a fraction of HA-HIF-1α co-immunoprecipitates with Xpress-hARD2 (Figure 5B).
Subcellular localization of hARD2
N-α-acetyltransferase activity of hARD2
hARD1 and hARD2 are differently regulated during NB4 cell differentiation
We here report on the presence of a human gene encoding a protein acetyltransferase at the mRNA, protein and enzyme activity levels. The ARD2 acetyltransferase has originated through a retrotransposition event and represents an additional ARD gene that can evolve independently from ARD1. Our results suggest that the duplication event is eutherian mammal specific. Amino acid sequence differences between hARD1 and hARD2 demonstrate that changes already have taken place, particularly in the C-terminal region of hARD2. As proposed for hARD1, hARD2 could have several different roles in the cell. The interaction with NATH indicates a role in the cytoplasm, possibly cotranslational N-acetylation. It is likely that hARD1 and hARD2 bind the same region within NATH and therefore compete for NATH binding. The difference in activity between hARD1 and hARD2 could suggest that the ACTH peptide used in this assay is not an optimal substrate for hARD2. It is not unlikely that the substrate preference is different between these two enzymes. A putative HIF-1α associated function is supported by the immunoprecipitation results demonstrating the potential interaction between hARD2 and HIF-1α (Figure 5B). Whether or not HIF-1α is a direct substrate for hARD2 mediated acetylation has not been investigated in this study and recent studies question the role of hARD1 in destabilizing HIF-1α [29–32]. Thus, the potential link between HIF-1α and hARD1/hARD2 remains an unresolved issue.
Unfortunately, the multiple sequence alignment outside of the acetyl transferase domain used for phylogenetic studies is not good enough for evolutionary studies involving calculation of substitution rates at synonymous and nonsynonymous sites or other likelihood ratio tests. Future sequencing from additional mammals should improve alignment quality to enable such studies. The ultimate question regarding differences in function between hARD1 and hARD2 is one of selective pressures and neofunctionalization vs. subfunctionalization (see for example, [44, 45]), which requires knowledge of the expression patterns and functions of the ancestral state. The metatherian, frog, and fish sequences, which are currently available, can serve as proxies for the ancestral state as they appear not to have duplicate ARD copies. In an attempt to address the function of zebrafish ARD, an antibody to hARD1 was directed against zebrafish tissue, but no binding was detected, probably due to sequence divergence (data not shown). Future work in zebrafish or other non-eutherian mammal species closer to human can address the function of the ancestral state and the relative selective pressures on hARD1 and hARD2.
There are several indications that hARD1 and NATH may be linked to differentiation. As mentioned, mNAT1 and mARD1 mRNAs are downregulated during neuronal differentiation in mouse . In the present study, we show that hARD1 and NATH proteins were both downregulated after induction of differentiation in promyelocytic leukaemia cells. This is the first description of NAT-ARD1 downregulation induced by differentiation in human cells and also the first verification of downregulation at the protein level. Interestingly, the hARD1 gene was one of twelve genes identified to be elevated in dedifferentiated hepatocellular carcinomas . It should also be noted that the yeast ARD1 gene originally was implicated in controlling the switch between the mitotic cell cycle and developmental pathways . Whether or not hARD1 and NATH influences differentiation per se, awaits further studies, but our present results add support to a link between the NATH-hARD1 complex and differentiation. The lack of hARD2 downregulation during granulocytic differentiation could suggest that a specific subset of proteins then is acetylated at a stable level. In contrast, the hARD1 specific substrates would be less acetylated during differentiation. This balance between acetylation of different subsets of proteins could have an impact on the differentiation process itself or alternatively facilitate cellular adaptations associated with the process.
In summary, we have identified and characterized a human protein N-α-acetyltransferase which probably originates from a recent gene duplication event. The hARD2 protein displays similar properties as hARD1 in terms of subcellular localization and potential interactions with NATH and HIF-1α. Further expression and enzymatic studies are required to assess the overall functional contribution of hARD2 in cellular processes.
hARD2 cloning and expression
NCBI BLAST database was employed in the search of human homologues of yeast ARD1. The second best hit after hARD1 was the protein MGC10646, termed hARD2. Primers pr1 and pr3 were designed to clone the hARD2 gene and to detect the presence of hARD2 mRNA. Plasmid encoding V5- and Xpress-tagged hARD2 was constructed from a gene-specific PCR of cDNA made from total RNA isolated from human ARO cells. PCR product was inserted into the TOPO TA vectors pcDNA3.1/V5His and pcDNA 4/HisMax (Invitrogen). Simultaneously a vector encoding wildtype hARD2 was made using a reverse primer including the hARD2 stop codon. cDNA was made as previously described . The primers for amplifying the hARD2 gene (Figure 1) were as follows: pr1 (hARD2 forward), 5'-ATG AAC ATC CGC AAC GCT CAG-3'; pr3 (hARD2 rev), 5'- GGA GGT GGA ATC CGA GCT TTC-3'; pr2 (hARD2 545forward), 5'-CAG CAC ACT TTC TGA TTC TGA AG-3'; pr4 (hARD2 926reverse), 5'-GTA ATG GCA GGT CTC AAA GTTC-3'. β-Actin primers: Actin-F, 5'- GGC ACC ACA CCT TCT ACA 3'; Actin-R, 5'- AGG AAG GCT GGA AGA GTG 3'. Primers for amplifying genomic M. giganteus ARD: For: 5'-gtg aar cgy tcn cac cgg cgc cty ggy ctg-3'; Rev: 5'-ctc ttc ctg rcr tgy agr gas acr tay ttk gc-3'. Primers for amplifying M. giganteus ARD from cDNA: kARD1/2 ORFf, 5'-ATG AAC ATC CGC AAY GCK MRG CCA GAS GAC C; kARD1/2 ORFr, 5'-CTA GGA GGY KGA RTC SGA GSY YTC TGA GCT GTC C (r = a + g, y = c + t, n = a + c + t + g, s = c + g, k = g + t).
Cell culture and transfection
Cells were cultured at 37°C, 5% CO2 in DMEM (HEK293, embryonal kidney, ATCC CRL-1573 and HeLa, epithelial cervix adenocarcinoma, ATCC CCL-2), RPMI 1640 (NB4, acute promyelocytic leukemia, DSMZ ACC 207 and MCF-7, epithelial mammary gland, breast adenocarcinoma, ATCC HTB-22) or EMEM (SK-MEL2, malignant melanoma, ATCC HTB-68) supplemented with 10% FBS and 3% l-glutamine. Transfections were performed using Fugene6 (Roche) according to the instruction manual. The plasmid pHA-HIF-1α has been described . A hARD2 specific rabbit antibody was generated by Biogenes GmbH using a peptide corresponding to amino acids 192–206 of hARD2. Western blotting was performed as described . Dilutions: anti-hARD2 1:500; anti-hARD1  1:500; anti-NATH  1:500; anti-V5 (Invitrogen) 1:1000; anti-β-tubulin (Sigma) 1:1000; anti-HA (Sigma) 1:1000.
Immunofluorescence and immunoprecipitation
HeLa cells were transfected using Fugene6 and grown 24 h on coverslips. Then cells were prepared for immunofluorescence as described . HEK293 or MCF-7 cells (~2 × 106) were transfected using Fugene6 and incubated 48 hours before harvesting and lysis in 300 μl lysis buffer. Immunoprecipitation was performed as previously described .
Immunoprecipitation of X-press-lacZ, X-press-hARD1 or X-press-hARD2 was performed as described above. Pellets of Protein A/G-Agarose bound X-press-hARD2 was added 10 μl ACTH (0.5 mM, human corticotropin fragment 1–24, Calbiochem), 4 μl [3H]Acetyl-CoA (1 μCi, 107 GBq/mmol, Amersham Biosciences) and 136 μl 0.2 M K2HPO4 (pH 8.1). The mixture was incubated for 2 hours at 37°C. After centrifugation the supernatant was added to 150 μl SP Sepharose (50% slurry in 0.5 M Acetic acid, Sigma) and incubated on a rotor for 5 min. The mixture was centrifuged and the pellet was washed three times with 0.5 M acetic acid and finally with methanol. Radioactivity in the ACTH containing pellet was determined by scintillation counting.
Alignment and tree building
Peptide sequence alignments were made using T-Coffee  with the default settings. Coding sequence alignments were produced by aligning the coding sequences with reference to the alignment of the corresponding peptide sequences. The phylogenetic tree from the initial peptide sequence alignment was built using MrBayes  with the Jones matrix of substitution run on four MCMC chains for 750000 generations, to generate a majority rule consensus tree of all compatible partitions of the final 500000 generations (sampled rate 100). MrBayes was also used to produce the coding sequence alignment based tree, with the same settings as before except for a nucleotide model of substitution allowing different rates for transitions and transversions and a gamma distribution of rates across sites. Gene trees were rooted by mapping them on to the NCBI tree of life whilst minimizing gene duplication and loss events and allowing poorly supported branches (those with posterior probabilities less than 0.7) to be rearranged according to the NCBI taxonomy as a reference species tree (Berglund, Steffansson, Betts and Liberles, Manuscript submitted).
hypoxia inducible factor-1α
polymerase chain reaction.
We thank C. Baldersheim, K. Jacobsen and C. Hoff for technical assistance and S. Knappskog for advice. R. McKinnon, A. El-Merhibi, G.H. Robstad, Kristiansand Zoo and C. Krakstad are thanked for providing kangaroo and rodent specimens. This work was supported by The Norwegian Cancer Society (Grants to TA, JEV, JRL), The Locus of Experimental Cancer Research (University of Bergen), The Meltzer Foundation (Grant to TA), FUGE and the International Agency for Research on Cancer (grant to FP).
- Park EC, Szostak JW: Ard1 and Nat1 Proteins Form A Complex That Has N-Terminal Acetyltransferase Activity. EMBO J. 1992, 11 (6): 2087-93.PubMed CentralPubMed
- Polevoda B, Norbeck J, Takakura H, Blomberg A, Sherman F: Identification and specificities of N-terminal acetyltransferases from Saccharomyces cerevisiae. EMBO J. 1999, 18: 6155-6168. 10.1093/emboj/18.21.6155.PubMed CentralPubMedView Article
- Arnold RJ, Polevoda B, Reilly JP, Sherman F: The action of N-terminal acetyltransferases on yeast ribosomal proteins. J Biol Chem. 1999, 274: 37035-37040. 10.1074/jbc.274.52.37035.PubMedView Article
- Polevoda B, Sherman F: N-terminal acetyltransferases and sequence requirements for N-terminal acetylation of eukaryotic proteins. J Mol Biol. 2003, 325: 595-622. 10.1016/S0022-2836(02)01269-X.PubMedView Article
- Polevoda B, Sherman F: Composition and function of the eukaryotic N-terminal acetyltransferase subunits. Biochem Biophys Res Commun. 2003, 308: 1-11. 10.1016/S0006-291X(03)01316-0.PubMedView Article
- Arnesen T, Anderson D, Baldersheim C, Lanotte M, Varhaug JE, Lillehaug JR: Identification and characterization of the human ARD1-NATH protein acetyltransferase complex. Biochem J. 2005, 386: 433-443. 10.1042/BJ20041071.PubMed CentralPubMedView Article
- Driessen HP, de Jong WW, Tesser GI, Bloemendal H: The mechanism of N-terminal acetylation of proteins. CRC Crit Rev Biochem. 1985, 18: 281-325.PubMedView Article
- Gautschi M, Just S, Mun A, Ross S, Rucknagel P, Dubaquie Y, Ehrenhofer-Murray A, Rospert S: The yeast N(alpha)-acetyltransferase NatA is quantitatively anchored to the ribosome and interacts with nascent polypeptides. Mol Cell Biol. 2003, 23: 7403-7414. 10.1128/MCB.23.20.7403-7414.2003.PubMed CentralPubMedView Article
- Pestana A, Pitot HC: Acetylation of nascent polypeptide chains on rat liver polyribosomes in vivo and in vitro. Biochemistry. 1975, 14: 1404-1412. 10.1021/bi00678a010.PubMedView Article
- Yamada R, Bradshaw RA: Rat liver polysome N alpha-acetyltransferase: isolation and characterization. Biochemistry. 1991, 30: 1010-1016. 10.1021/bi00218a018.PubMedView Article
- Tribioli C, Mancini M, Plassart E, Bione S, Rivella S, Sala C, Torri G, Toniolo D: Isolation of new genes in distal Xq28: transcriptional map and identification of a human homologue of the ARD1 N-acetyl transferase of Saccharomyces cerevisiae. Hum Mol Genet. 1994, 3: 1061-1067.PubMedView Article
- Jeong JW, Bae MK, Ahn MY, Kim SH, Sohn TK, Bae MH, Yoo MA, Song EJ, Lee KJ, Kim KW: Regulation and destabilization of HIF-1alpha by ARD1-mediated acetylation. Cell. 2002, 111: 709-720. 10.1016/S0092-8674(02)01085-1.PubMedView Article
- Sugiura N, Adams SM, Corriveau RA: An evolutionarily conserved N-terminal acetyltransferase complex associated with neuronal development. J Biol Chem. 2003, 278: 40113-40120. 10.1074/jbc.M301218200.PubMedView Article
- Williams BC, Garrett-Engele CM, Li ZX, Williams EV, Rosenman ED, Goldberg ML: Two putative acetyltransferases, San and deco, are required for establishing sister chromatid cohesion in drosophila. Curr Biol. 2003, 13: 2025-2036. 10.1016/j.cub.2003.11.018.PubMedView Article
- Arnesen T, Anderson D, Torsvik J, Halseth HB, Varhaug JE, Lillehaug JR: Cloning and characterization of hNAT5/hSAN: An evolutionarily conserved component of the NatA protein N-alpha-acetyltransferase complex. Gene. 2006, 371: 291-295. 10.1016/j.gene.2005.12.008.PubMedView Article
- Brown JL, Roberts WK: Evidence That Approximately 80 Per Cent of Soluble-Proteins from Ehrlich Ascites-Cells Are N-Alpha-Acetylated. J Biol Chem. 1976, 251: 1009-1014.PubMed
- Brown JL: A comparison of the turnover of alpha-N-acetylated and nonacetylated mouse L-cell proteins. J Biol Chem. 1979, 254: 1447-1449.PubMed
- Jõnvall H: Acetylation of protein N-terminal amino groups: structural observations on alpha-amino acetylated proteins. J Theor Biol. 1975, 55 (1): 1-12.View Article
- Lee FJ, Lin LW, Smith JA: N alpha-acetyltransferase deficiency alters protein synthesis in Saccharomyces cerevisiae. FEBS Lett. 1989, 256: 139-142. 10.1016/0014-5793(89)81734-X.PubMedView Article
- Persson B, Flinta C, Vonheijne G, Jornvall H: Structures of N-Terminally Acetylated Proteins. Eur J Biochem. 1985, 152: 523-527. 10.1111/j.1432-1033.1985.tb09227.x.PubMedView Article
- Geissenhoner A, Weise C, Ehrenhofer-Murray AE: Dependence of ORC silencing function on NatA-mediated Nalpha acetylation in Saccharomyces cerevisiae. Mol Cell Biol. 2004, 24: 10300-10312. 10.1128/MCB.24.23.10300-10312.2004.PubMed CentralPubMedView Article
- Wang X, Connelly JJ, Wang CL, Sternglanz R: Importance of the Sir3 N terminus and its acetylation for yeast transcriptional silencing. Genetics. 2004, 168: 547-551. 10.1534/genetics.104.028803.PubMed CentralPubMedView Article
- Lee FJS, Lin LW, Smith JA: N-Alpha Acetylation Is Required for Normal Growth and Mating of Saccharomyces-Cerevisiae. J Bact. 1989, 171: 5795-5802.PubMed CentralPubMed
- Mullen JR, Kayne PS, Moerschell RP, Tsunasawa S, Gribskov M, Colavitoshepanski M, Grunstein M, Sherman F, Sternglanz R: Identification and Characterization of Genes and Mutants for An N-Terminal Acetyltransferase from Yeast. EMBO J. 1989, 8: 2067-2075.PubMed CentralPubMed
- Whiteway M, Szostak JW: The Ard1 Gene of Yeast Functions in the Switch Between the Mitotic Cell-Cycle and Alternative Developmental Pathways. Cell. 1985, 43: 483-492. 10.1016/0092-8674(85)90178-3.PubMedView Article
- Whiteway M, Freedman R, Vanarsdell S, Szostak JW, Thorner J: The Yeast Ard1 Gene-Product Is Required for Repression of Cryptic Mating-Type Information at the Hml Locus. Mol Cell Biol. 1987, 7: 3713-3722.PubMed CentralPubMedView Article
- Sonnichsen B, Koski LB, Walsh A, Marschall P, Neumann B, Brehm M, Alleaume AM, Artelt J, Bettencourt P, Cassin E, Hewitson M, Holz C, Khan M, Lazik S, Martin C, Nitzsche B, Ruer M, Stamford J, Winzi M, Heinkel R, Roder M, Finell J, Hantsch H, Jones SJM, Jones M, Piano F, Gunsalus KC, Oegema K, Gonczy P, Coulson A, Hyman AA, Echeverri CJ: Full-genome RNAi profiling of early embryogenesis in Caenorhabditis elegans. Nature. 2005, 434: 462-469. 10.1038/nature03353.PubMedView Article
- Arnesen T, Gromyko D, Pendino F, Ryningen A, Varhaug JE, Lillehaug JR: Induction of apoptosis in human cells by RNAi-mediated knockdown of hARD1 and NATH, components of the protein N-alpha-acetyltransferase complex. Oncogene. 2006, PMID: 16518407
- Arnesen T, Kong X, Evjenth R, Gromyko D, Varhaug JE, Lin Z, Sang N, Caro J, Lillehaug JR: Interaction between HIF-1alpha (ODD) and hARD1 does not induce acetylation and destabilization of HIF-1alpha. FEBS Lett. 2005, 579: 6428-6432. 10.1016/j.febslet.2005.10.036.PubMed CentralPubMedView Article
- Bilton R, Mazure N, Trottier E, Hattab M, Dery MA, Richard DE, Pouyssegur J, Brahimi-Horn C: Arrest-defective-1 protein, an acetyltransferase, does not alter stability of hypoxia-inducible factor (HIF)-1alpha and is not induced by hypoxia or HIF. J Biol Chem. 2005, 280: 31132-31140. 10.1074/jbc.M504482200.PubMedView Article
- Fisher TS, Des Etages S, Hayes L, Crimin K, Li BY: Analysis of ARD1 function in hypoxia response using retroviral RNA interference. J Biol Chem. 2005, 280: 17749-17757. 10.1074/jbc.M412055200.PubMedView Article
- Arnesen T: HIF1 alpha and ARD1: enemies, friends or neither?. Nat Rev Cancer. 2006, 6 (3): doi:10.1038/nrc1779-c1 MAR 2006
- Bristow J, Gitelman SE, Tee MK, Staels B, Miller WL: Abundant Adrenal-Specific Transcription of the Human P450C21A Pseudogene. J Biol Chem. 1993, 268: 12919-12924.PubMed
- Kandouz M, Bier A, Carystinos GD, Alaoui-Jamali MA, Batist G: Connexin43 pseudogene is expressed in tumor cells and inhibits growth. Oncogene. 2004, 23: 4763-4770. 10.1038/sj.onc.1207506.PubMedView Article
- Mccarrey JR, Thomas K: Human Testis-Specific Pgk Gene Lacks Introns and Possesses Characteristics of A Processed Gene. Nature. 1987, 326: 501-505. 10.1038/326501a0.PubMedView Article
- Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28: 4364-4375. 10.1093/nar/28.21.4364.PubMed CentralPubMedView Article
- Roth C, Betts MJ, Steffansson P, Saelensminde G, Liberles DA: The Adaptive Evolution Database (TAED): a phylogeny based tool for comparative genomics. Nucleic Acids Res. 2005, D495-D497. 33 Database
- Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, D138-D141. 10.1093/nar/gkh121. 32 Database
- Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP, Bork P: SMART 4.0: towards genomic data integration. Nucleic Acids Res. 2004, 32: D142-D144. 10.1093/nar/gkh088.PubMed CentralPubMedView Article
- Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.PubMedView Article
- Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003, 31: 3381-3385. 10.1093/nar/gkg520.PubMed CentralPubMedView Article
- Kraulis PJ: Molscript – A Program to Produce Both Detailed and Schematic Plots of Protein Structures. J Appl Cryst. 1991, 24: 946-950. 10.1107/S0021889891004399.View Article
- Linding R, Russell RB, Neduva V, Gibson TJ: GlobPlot: exploring protein sequences for globularity and disorder. Nucleic Acids Res. 2003, 31: 3701-3708. 10.1093/nar/gkg519.PubMed CentralPubMedView Article
- He XL, Zhang JZ: Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics. 2005, 169: 1157-1164. 10.1534/genetics.104.037051.PubMed CentralPubMedView Article
- Rastogi S, Liberles DA: Subfunctionalization of duplicated genes as a transition state to neofunctionalization. BMC Evol Biol. 2005, 5: 28-10.1186/1471-2148-5-28.PubMed CentralPubMedView Article
- Midorikawa Y, Tsutsumi S, Taniguchi H, Ishii M, Kobune Y, Kodama T, Makuuchi M, Aburatani H: Identification of genes associated with dedifferentiation of hepatocellular carcinoma with expression profiling analysis. Jpn J Canc Res. 2002, 93 (6): 636-43.View Article
- Arnesen T, Gromyko D, Horvli O, Fluge O, Lillehaug J, Varhaug JE: Expression of N-acetyl transferase human and human arrest defective 1 proteins in thyroid neoplasms. Thyroid. 2005, 15: 1131-1136. 10.1089/thy.2005.15.1131.PubMedView Article
- Yu F, White SB, Zhao Q, Lee FS: HIF-1 alpha binding to VHL is regulated by stimulus-sensitive proline hydroxylation. Proc Natl Acad Sci U S A. 2001, 98 (17): 9630-5. 10.1073/pnas.181341498.PubMed CentralPubMedView Article
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-217. 10.1006/jmbi.2000.4042.PubMedView Article
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.