MBE Advance Access originally published online on April 16, 2008
Molecular Biology and Evolution 2008 25(7):1482-1487; doi:10.1093/molbev/msn096
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
Molecular Cloning and Characterization of a Moss (Ceratodon purpureus) Nonsymbiotic Hemoglobin Provides Insight into the Early Evolution of Plant Nonsymbiotic Hemoglobins
Laboratorio de Biofísica y Biología Molecular, Facultad de Ciencias, Universidad Autónoma del Estado de Morelos, Cuernavaca, Morelos, México
E-mail: ra{at}buzon.uaem.mx.
| Abstract |
|---|
|
|
|---|
Nonsymbiotic hemoglobins (nsHbs) are widespread in plants including bryophytes. Bryophytes (such as mosses) are among the oldest land plants, thus an analysis of a bryophyte nsHb is of interest from an evolutionary perspective. However, very little is known about bryophyte nsHbs. Here, we report the cloning and characterization of an nshb gene (cerhb) from the moss Ceratodon purpureus. Sequence analysis showed that cerhb is interrupted by 3 introns in identical position as all known plant nshb genes, which suggests that the ancestral nshb gene was interrupted by 3 introns. Expression analysis showed that cerhb expresses in protonemas and gametophytes growing in normal conditions and that it overexpresses in protonemas subjected to osmotic (sucrose), heat-shock, cold-, and nitrate-stress conditions. Also, modeling of the Ceratodon nsHb (CerHb) tertiary structure suggests that CerHb is hexacoordinate and that it binds O2 with high affinity. Comparative analysis of the predicted CerHb with native rice Hb1 and soybean leghemoglobin a structures revealed that the major evolutionary changes that probably occurred during the evolution of plant Hbs were 1) a hexacoordinate to pentacoordinate transition at the heme prosthetic group, 2) a length decrease at the CD-loop and N- and C-termini regions, and 3) the compaction of the protein into a globular structure.
Key Words: Ceratodon evolution gene expression globin molecular modeling nonsymbiotic
| Introduction |
|---|
|
|
|---|
Nonsymbiotic hemoglobins (nsHbs) are plant proteins that reversibly bind O2 and other gaseous ligands. These proteins have been detected in seeds, roots, stems, leaves, and flowers from plants growing in normal conditions (Taylor et al. 1994
The nsHbs have been identified in bryophyte, gymnosperm, and angiosperm species and are therefore widespread in land plants (Garrocho-Villegas et al. 2007
). Gene and protein sequence comparison suggests that nsHbs-1 and nsHbs-2 evolved from a common ancestor; however, it is still unclear if an nsHb-1–like or nsHb-2–like protein was the ancestor to plant nsHbs. Full bryophyte cDNA sequences coding for moss (Physcomitrella patens and Ceratodon purpureus) nsHbs and a partial gene sequence coding for a liverwort (Marchanthia polymorpha) nsHb were deposited in the GenBank database (with accession numbers AF218049
[GenBank]
, AF309562
[GenBank]
, and AY026341
[GenBank]
, respectively). The phenetic analysis of plant Hb sequences, including those from moss nsHbs, showed that bryophyte nsHbs are ancestral but intermediate to nsHbs-1 and nsHbs-2 (Garrocho-Villegas et al. 2007
; Gopalasubramaniam, Kovacs, et al. 2008
). Thus, the apparent oldest land plant nsHbs are bryophyte nsHbs. However, despite limited reports (Arredondo-Peter et al. 2000
; Trent et al. 2001
), bryophyte nsHbs remain uncharacterized and thus very little is known about the properties of primeval plant nsHbs. Here, we report the cloning and molecular characterization of a moss (C. purpureus) nsHb, including an analysis of gene expression in plant organs and modeling of the protein tertiary structure.
| Materials and Methods |
|---|
|
|
|---|
Plant Material and Growing Conditions
Ceratodon purpureus protonemas were kindly provided by Dr Ralf Reski (Albert-Ludwigs-Universität, Freiburg, Germany) as an axenic culture. Protonemas were cultured in Knop medium (Reski and Abel 1985
The (nshb) gene expression in plants subjected to stress conditions was evaluated in 2-week-old protonemas by duplicate. Osmotic-, salt-, and nitrate-stress conditions were generated by adding 10% sucrose (Umeda et al. 1994
), 250 mM NaCl (Umeda et al. 1994
), and 85 mM Ca(NO3)2 (Ohwaki et al. 2005
) to the growth medium, respectively. Cold-stress conditions were generated by incubating protonemas at 4 °C for 24 h (Takahashi et al. 1994
). Dark-stress conditions were generated by incubating the protonemas in darkness for 24 h inside a plant growth chamber (Lira-Ruan et al. 2001
). Flooding-stress conditions were generated by fully covering the protonemas with distilled water for 24 h (Taylor et al. 1994
; Lira-Ruan et al. 2001
). Heat-shock stress conditions were generated by incubating the protonemas at 42 °C for 2 h (Higo K and Higo H 1993
). Control protonemas (i.e., those growing under normal conditions) were cultivated for 2 weeks, and control and stressed protonemas were harvested at the same time.
Total DNA and Poly(A+) RNA Isolation
Total DNA was isolated from
100 mg of C. purpureus protonemas using the cetyltrimethylammoniun bromide (CTAB) method (Crose and Amorese 1978
). Poly(A+) RNA was isolated from plant tissues using the QuickPrep Micro mRNA Purification kit (Amersham Biosciences, Little Chalfont Buckinghamshire, United Kingdom) following the manufacturer's instructions. Total DNA and poly(A+) RNA were quantitated by spectrophotometry assuming 1 A260 = 50 or 40 µg/ml for DNA or RNA, respectively (Ausubel et al. 1995
).
Polymerase Chain Reaction Amplification, Gene Cloning, and DNA Sequencing
Primers were designed for polymerase chain reaction (PCR) to amplify the C. purpureus nshb gene using sequences at the start and stop codons of the C. purpureus nsHb cDNA (GenBank accession number AF309562). The primer sequences (degenerated for the NcoI and EcoRI restriction sites [underlined]) were (sense) primer Ceratodon nsHb (CerHb)/ATG-NcoI: 5'-CCATGGCACCACCGACAGTCG-3' and (antisense) primer CerHb/TAA-EcoRI: 5'-GAATTCTTACTGAGCAGCAGCCCTC-3'. Total C. purpureus DNA (
1 µg) was used as a template for PCR amplification. PCR components and concentrations were 1 µM of each sense and antisense primer, 400 µM of each deoxynucleotide triphosphate, and 1 U of Taq DNA polymerase (Invitrogen, Carlsbad, CA) in 1X PCR buffer containing 4 mM MgCl2 in a final volume of 25 µl. PCR amplification was carried out for 35 cycles at 55 °C for annealing using a thermalcycler (Minicycler, MJ Research, Watertown, MA). PCR products were detected in a 1.2% agarose gel after staining with ethidium bromide, isolated from the gel using the GeneClean kit (Q-BIOgene, Carlsbad, CA), and cloned into the pCR2.1 cloning vector (Invitrogen) following the manufacturer's instructions. Insert DNA was fully sequenced in both orientations at the Molecular Biology Facility of the Cell Physiology Institute of the National Autonomous University of México.
Southern Blot Analysis
Southern blot analysis of the total C. purpureus DNA was done as described by Arredondo-Peter et al. (1997)
using the C. purpureus nshb gene as a probe. Membranes were prehybridized and hybridized at 55 °C and washed at high stringency (60 °C) in 2x standard saline citrate (SSC)/0.1% sodium dodecyl sulfate (SDS) twice for 5 min each and 0.5x SSC/0.1% SDS twice for 15 min each.
Gene Expression Analysis by Reverse Transcriptase–Polymerase Chain Reaction
Expression of the C. purpureus nshb gene was examined in protonemas and gametophytes (4 and 10 weeks old, respectively) and in protonemas grown under stress conditions (see above) by reverse transcriptase–polymerase chain reaction (RT-PCR) using the GeneAmp RNA PCR kit (Applied Biosystems, Branchburg, NJ) as described by Arredondo-Peter et al. (1997)
. Both reverse transcription and PCR amplification were performed using the CerHb/ATG-NcoI and CerHb/TAA-EcoRI primers and 100 ng of template poly(A+) RNA. Amplification conditions were identical to those used for the amplification of the C. purpureus nshb gene (see above). After gel electrophoresis, PCR fragments were stained with ethidium bromide, and the gel was photodocumented, scanned, and analyzed using the ImageJ 1.36b software (http://rsb.info.nih.gov/ij/). The relative abundance of the CerHb transcripts in control and differentially treated protonemas was calculated from the intensity of each band using as reference the intensity detected in nsHb amplicon observed in protonemas grown under normal conditions.
In Silico Analysis and Molecular Modeling
Multiple sequence alignment and cluster analysis of the C. purpureus nsHb and selected plant Hbs were performed using the Neighbor-Joining method of the ClustalX program (Thompson et al. 1997
). The tertiary structure of the C. purpureus nsHb was predicted by homology modeling using the crystal structure of the rice Hb1 (Hargrove et al. 2000
) (PDB ID 1D8U
[PDB]
) as a template as described by Gopalasubramaniam, Garrocho-Villegas, et al. (2008)
.
| Results and Discussion |
|---|
|
|
|---|
Cloning of a C. purpureus nshb Gene
The nsHbs evolved alongside land plants after bryophyte-like organisms colonized land about 450 MYA. Thus, it is likely that these proteins played roles in plant adaptation to the terrestrial environment (Ross et al. 2002
Characterization of the CerHb Protein Sequence
The cerhb gene codes for a predicted protein of 177 amino acid residues in length (supplementary fig. S2, Supplementary Material online) with a calculated molecular mass of 19,535 Da. Sequence analysis showed that CerHb contains distal (H86) and proximal (H121) His and Phe B10 (F52), which are highly conserved in plant Hbs. However, in Ceratodon and Physcomitrella nsHbs (see below), the CD1 position is occupied by Tyr (Y66) instead of Phe, which is conserved in evolved nsHbs (supplementary fig. S3, Supplementary Material online). This observation shows that Tyr in ancestral nsHbs occupied the CD1 position and was replaced by Phe during the evolution of nsHbs.
Arredondo-Peter et al. (2000)
identified and partially characterized an nsHb sequence from the moss P. patens. Multiple sequence alignment showed that moss (C. purpureus and P. patens) nsHbs are the longest known plant Hbs and that the size of plant nsHbs decreased over time, mostly due to a loss of amino acid residues at the N-terminal (prehelix A) region (supplementary fig. S3, Supplementary Material online) (Ross et al. 2002
). Also, sequence alignment revealed that the C. purpureus and P. patens nsHbs are 80% identical and that the highest sequence variability occurs at the N-termini regions (i.e., within the first 35 amino acid residues). The identity value of moss nsHbs is similar to that observed between evolved nsHbs. For example, the identity value between rice and maize nsHb1s is approximately 80%. This observation suggests that nsHbs evolved at approximately the same rate and that major changes during plant nsHbs evolution occurred at the N-terminal region (see below).
Identification of the cerhb Gene Copy Number
The number of nshb gene copies varies with each species in higher plants. For example, single and multiple nshb gene copies exist in barley (Taylor et al. 1994
) and rice (Arredondo-Peter et al. 1997
; Lira-Ruan et al. 2002
), respectively. Also, the search for homologs to the globin fold in the SUPERFAMILY database (Gough et al. 2001
) identified several Hb sequences in the P. patens genome (data not shown). When total C. purpureus DNA was subjected to Southern blotting with the cerhb probe, 2 hybridizing fragments of 4 and 6–7 kb in length were detected in the DNA digested with EcoRI (fig. 1). This observation suggests that at least 2 copies of the nshb gene exist in the C. purpureus genome. The existence of cDNAs coding for CerHb indicates that cerhb is functional in C. purpureus.
|
Expression of the cerhb Gene in C. purpureus Grown under Normal and Stress Conditions
The nsHb transcripts and proteins have been reported to exist in several organs from plants growing under normal and stress conditions (Taylor et al. 1994
|
Predicted Tertiary Structure of CerHb
Nothing is known about the tertiary structure of bryophyte nsHbs. Understanding the structural properties of nsHbs is important to clarifying the role these proteins play in plant organs. In silico methods exist to predict with high reliability the tertiary structure of proteins from template structures (Sáenz-Rivera et al. 2004
-helix exists within the prehelix A region of predicted CerHb, which may facilitate protein translocation. Ross et al. (2002)
-helix at the prehelix A region in predicted CerHb (fig. 3a) is consistent with the possibility that CerHb (and other bryophyte nsHbs) may function inside cell organelles. If this is correct, it is possible that nsHbs became cytoplasmic during evolution (Ross et al. 2002
|
The Structural Evolution of CerHb and Plant nsHbs
Comparative analysis of the CerHb, rice Hb1, and soybean leghemoglobin a (Lba) tertiary structures (i.e., primitive to evolved plant Hbs) reveals major evolutionary changes that probably occurred during plant Hbs evolution (fig. 4). These changes include 1) a hexacoordinate to pentacoordinate transition at the heme prosthetic group (Hoy et al. 2007
|
| Supplementary Material |
|---|
|
|
|---|
Supplementary figures S1–S3 are available at Molecular Biology and Evolutiononline (http://www.mbe.oxfordjournals.org/).
| Acknowledgements |
|---|
|
|
|---|
The authors wish to express their gratitude to Dr Ralf Reski (Albert-Ludwigs-Universität, Freiburg, Germany) for kindly providing the C. purpureus protonemas, to an anonymous reviewer for useful suggestions to improve the contents of this article, and to Miss Gillian Klucas for corrections made to improve the English language. This work was funded by Consejo Nacional de Ciencia y Tecnología (CoNaCyT project no. 42873Q), México. V.G.-V. is a postdoctoral fellow financed by CoNaCyT (IdAP 9272).
| Footnotes |
|---|
Geoffrey McFadden, Associate Editor
| References |
|---|
|
|
|---|
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. (1990) 215:403–410.[CrossRef][Web of Science][Medline]
Andersson CR, Jensen EO, Llewellyn DJ, Dennis ES, Peacock WJ. A new hemoglobin gene from soybean: a role for hemoglobin in all plants. Proc Natl Acad Sci USA (1996) 93:5682–5687.
Arredondo-Peter R, Hargrove MS, Moran JF, Sarath G, Klucas RV. Plant hemoglobins. Plant Physiol. (1998) 118:1121–1126.
Arredondo-Peter R, Hargrove MS, Sarath G, Moran JF, Lohrman J, Olson JS, Klucas RV. Rice hemoglobins: gene cloning, analysis and oxygen-binding kinetics of a recombinant protein synthesized in Escherichia coli. Plant Physiol. (1997) 115:1259–1266.[Abstract]
Arredondo-Peter R, Ramírez M, Sarath G, Klucas RV. Sequence analysis of an ancient hemoglobin cDNA isolated from the moss Physcomitrella patens (accession no. AF218049). Plant Physiol. (2000) 122:1457.
Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA, Struhl K. Current protocols in molecular biology (1995) New York: Wiley.
Crose J, Amorese B. Isolation of plant DNA from fresh tissue. Focus (1978) 9:3–6.[Medline]
Duff SMG, Wittenberg JB, Hill RD. Expression, purification and properties of recombinant barley (Hordeum sp.) hemoglobin: optical spectra and reactions with gaseous ligands. J Biol Chem. (1997) 272:16746–16752.
Garrocho-Villegas V, Gopalasubramaniam SK, Arredondo-Peter R. Plant hemoglobins: what we know six decades after their discovery. Gene (2007) 398:78–85.[CrossRef][Web of Science][Medline]
Gopalasubramaniam SK, Garrocho-Villegas V, Bustos G, Pastor N, Arredondo-Peter R. Use of in silico (computer) methods to predict and analyze the tertiary structure of plant hemoglobins. Methods Enzymol. (2008) 436:393–410.[Web of Science][Medline]
Gopalasubramaniam SK, Kovacs F, Violante-Mota F, Twigg P, Arredondo-Peter R, Sarath G. Cloning and characterization of a caesalpinoid (Chamaecrista fasciculata) hemoglobin: the structural transition from a nonsymbiotic hemoglobin to a leghemoglobin. Proteins (2008) doi: 10.1002/prot.21917.
Gough J, Karplus K, Hughey R, Chothia C. Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J Mol Biol. (2001) 313:903–919.[CrossRef][Web of Science][Medline]
Hargrove M, Brucker EA, Stec B, Sarath G, Arredondo-Peter R, Klucas RV, Olson JS, Philips GN Jr. Crystal structure of a non-symbiotic hemoglobin. Structure (2000) 8:1005–1014.[Medline]
Hebelstrup KH, Igamberdiev AU, Hill RD. Metabolic effects of hemoglobin gene expression in plants. Gene (2007) 398:86–93.[CrossRef][Web of Science][Medline]
Higo K, Higo H. Changes in translatable mRNA populations induced in rice seedlings by exposure to freeze-thaw stress. Biochem Mol Biol Int. (1993) 29:403–410.[Web of Science][Medline]
Hoy JA, Robinson H, Trent JT III, Kakar S, Smagghe BJ, Hargrove MS. Plant hemoglobins: a molecular fossil record for the evolution of oxygen transport. J Mol Biol. (2007) 371:168–179.[CrossRef][Web of Science][Medline]
Igamberdiev AU, Baron K, Manac'h-Little N, Stoimenova M, Hill RD. The haemoglobin/nitric oxide cycle: involvement in flooding stress and effects on hormone signaling. Ann Bot (Lond). (2005) 96:557–564.
Lira-Ruan V, Ross E, Sarath G, Klucas RV, Arredondo-Peter R. Mapping and analysis of a hemoglobin gene family from rice (Oryza sativa). Plant Physiol Biochem. (2002) 40:199–202.[CrossRef][Web of Science]
Lira-Ruan V, Sarath G, Klucas RV, Arredondo-Peter R. Synthesis of hemoglobins in rice (Oryza sativa var. Jackson) plants growing in normal and stress conditions. Plant Sci. (2001) 161:279–287.[Medline]
Ohwaki Y, Kawagishi-Kobayashi M, Wakasa K, Fujihara S, Yoneyama T. Induction of class-1 non-symbiotic hemoglobin genes by nitrate, nitrite and nitric oxide in cultured rice cells. Plant Cell Physiol. (2005) 46:324–331.
Rensing SA, Lang D, Zimmer AD, et al, (67 co-authors). The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science (2008) 319:64–69.
Reski R, Abel WO. Induction of budding on chloronemata and caulonemanta of the moss Physcomitrella patens using isopentenyladenine. Planta (1985) 206:354–358.
Ross EJH, Lira-Ruan V, Arredondo-Peter R, Klucas RV, Sarath G. Recent insights into plant hemoglobins. Rev Plant Biochem Biotechnol. (2002) 1:173–189.
Ross EJH, Shearman L, Mathiesen M, Zhou J, Arredondo-Peter R, Sarath G, Klucas RV. Non-symbiotic hemoglobins are synthesized during germination and in differentiating cell types. Protoplasma (2001) 218:125–133.[CrossRef][Web of Science][Medline]
Ross EJH, Stone JM, Elowsky CG, Arredondo-Peter R, Klucas RV, Sarath G. Activation of the Oryza sativa non-symbiotic haemoglobin-2 promoter by the cytokinin-regulated transcription factor, ARR1. J Exp Bot. (2004) 55:1721–1731.
Sáenz-Rivera J, Sarath G, Arredondo-Peter R. Modeling the tertiary structure of a maize (Zea mays ssp. mays) non-symbiotic hemoglobin. Plant Physiol Biochem. (2004) 42:891–897.[CrossRef][Web of Science][Medline]
Takahashi R, Joshee N, Kitagawa Y. Induction of chilling resistance by water stress, and cDNA sequence analysis and expression of water stress-regulated genes in rice. Plant Mol Biol. (1994) 26:339–352.[CrossRef][Web of Science][Medline]
Taylor ER, Nie XZ, MacGregor AW, Hill RD. A cereal haemoglobin gene is expressed in seed and root tissues under anaerobic conditions. Plant Mol Biol. (1994) 24:853–862.[CrossRef][Web of Science][Medline]
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. (1997) 24:4876–4882.
Trent JT III, Hvitved AN, Hargrove MS. A model for ligand binding to hexacoordinate hemoglobins. Biochemistry (2001) 40:6155–6163.[CrossRef][Web of Science][Medline]
Trevaskis B, Watts RA, Andersson SR, Llewellyn DJ, Hargrove MS, Olson JS, Dennis ES, Peacock WJ. Two hemoglobin genes in Arabidopsis thaliana: the evolutionary origins of leghemoglobins. Proc Natl Acad Sci USA (1997) 94:12230–12234.
Umeda M, Hara C, Matsubayashi Y, Li H, Liu Q, Tadokoro F, Aotsuka S, Uchimiya H. Expressed sequence tags from cultured cells of rice (Oryza sativa L.) under stressed conditions: analysis of transcripts of genes engaged in ATP-generating pathways. Plant Mol Biol. (1994) 25:469–478.[CrossRef][Web of Science][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



