<title>Background</title> <p>β-Galactosidases (Bgals) (EC 3.2.1.23) catalyze the hydrolysis of the glycosidic bonds of terminal non-reducing β-D-galactosyl residues of oligosaccharides and β-D-galactopyranosides. This group of enzymes has a broad distribution, which encompasses plants, animals and microorganisms. The biological functions of these enzymes include degradation of structural polysaccharides in plant cell walls [##REF##10792826##1##,##REF##10852969##2##], thereby acting to control fruit softening during ripening [##REF##12867545##3##], hydrolysis of dietary lactose [##REF##5953213##4##,##REF##14246671##5##] and degradation of glycolipids and proteoglycans in mammals [##UREF##0##6##,##REF##3310533##7##], and metabolism of lactose and other galactosides in microorganisms [##REF##4124306##8##,##REF##14118##9##]. The extensive diversity of Bgals raises the question of what is the basis of their different substrate specificities, which help determine their functions.</p> <p>Glycosyl hydrolases (GH) have been grouped and classified as families on the basis of structural similarity [##REF##1747104##10##] and Bgals fall into glycosyl hydrolase (GH) families 1, 2, 35, and 42, which are part of superfamily A (or Clan A). Based on amino acid sequence similarities [##REF##8352747##11##], plant Bgals that have been described belong to GH family 35 (GH35). Like other glycosidase families, GH35 includes multiple genes in various plant species, including <italic>Arabidopsis </italic>[##REF##17466346##12##], tomato [##REF##10889266##13##], papaya [##REF##15694277##14##], apple [##REF##7991682##15##], <italic>Vigna unguiculata </italic>[##UREF##1##16##] and barley [##UREF##2##17##], suggesting that GH35 gene multiplicity is ubiquitous in plants.</p> <p>Analysis of plant Bgal activities led to the proposal of two classes, I and II [##REF##15980190##18##]. Class I is made up of enzymes with well characterized exo-β-(1→4)-galactanase activities that can specifically act on pectic β-(1→4)-galactan. Class II has hydrolytic activity on the β-(1→3)- and β-(1→6)-galactosyl linkages of arabinogalactan proteins (AGPs), but lacks activity toward pectic β-(1→4)-galactan [##REF##16666809##19##,##UREF##3##20##], as reported in spinach leaves [##UREF##3##20##] and hypocotyls and young leaves of radish [##REF##15980190##18##]. Class I OsBgals can specifically act on β-(1→4)-galactosyl residues in pectin and xyloglucan. Therefore, they could play an important role in modification of the architecture of the cell wall and intercellular attachment [##REF##10792826##1##,##REF##10852969##2##]. The apparent involvement of Bgals in pectin disassembly during fruit ripening has been reported in various plant species, including tomato [##REF##12867545##3##,##REF##16662771##21##,##REF##7630937##22##], muskmelon [##REF##16653123##23##], kiwifruit [##UREF##4##24##], persimmon [##REF##8058842##25##], sweet cherry [##UREF##5##26##], mango [##REF##7766393##27##], and peach [##REF##8058842##25##]. Bgals expression also accompanies many stages of plant development in other tissues, for example, spinach leaves [##UREF##3##20##], mungbean seedlings [##REF##11393513##29##], radish hypocotyls and young leaves [##REF##15980190##18##], and meristem zones of roots, cotyledons, vascular tissues, trichomes, and pollen of tobacco [##REF##16704113##30##,##REF##15517348##31##]. Bgals that appear to catalyze galactose removal from xyloglucans during their disassembly were observed in cotyledons of nasturtium (<italic>Tropueolum mujus </italic>L.) seed [##REF##3126187##32##], <italic>Copaifera langsdorffii </italic>seed [##UREF##6##33##], and <italic>Hymenaea courbaril </italic>seed, where they were reported to act in cooperation with α-xylosidase, β-glucosidase, and other enzymes to achieve xyloglucan degradation [##REF##10729610##34##,##UREF##7##35##]. Therefore, Bgals appear to play a role in cell wall remodeling in many plant processes.</p> <p>In rice, GH35 has not yet been well characterized, and the biological function of rice Bgals remain mysterious. Konno and Tsumuki [##UREF##8##36##] identified both soluble and cell-wall-bound Bgal, of which a 42 kDa soluble Bgal was purified and shown to release galactose from larch wood and rice cell wall arabinogalactans. Likewise, Kaneko and Kobayashi [##REF##12723614##37##] isolated a Bgal with 40 and 47 kDa subunits from the medium of rice suspension cells. Recently, Chantarangsee et al. [##UREF##9##38##] characterized two recombinant Bgal isozymes, including the 90 kDa OsBgal1 and 72 kDa OsBgal2, which had different expression patterns, though both are found throughout plant growth. The complete rice genome sequence allows extensive study of rice GH35, so the gene structures, encoded protein sequences and phylogenetic relationships with other Bgals from rice and other organisms, which may provide clues to their evolution and possible functions, were determined. In addition, the transcript expression patterns of all rice GH35 <italic>Bgals </italic>have been determined and OsBgal13 was expressed in <italic>Escherichia coli </italic>to gain a clue to OsBgal physiological functions.</p>
<title>Methods</title> <title>Sequence data and database search</title> <p>To find all GH35 genes in rice, tBLASTn searches [##REF##9254694##57##] were performed in the National Center for Biotechnology Information (NCBI) Genbank nr, indica rice genome, and expressed sequence tags (dbEST) databases with <italic>OsBgal1 </italic>(AK102192) as the query. The retrieved <italic>OsBgal </italic>genes were used to identify their Unigene cluster, Gene Locus, conserved domains, gene position on the 12 <italic>Oryza sativa </italic>(rice) chromosomes, homolog genes, EST-based expression profiles and GEO profile in the NCBI databases (please see Availability &amp; requirements for more information). The number and positions of exons and introns for each individual gene were determined by manually comparing the cDNA and predicted cDNA sequences with their corresponding genomic DNA sequences. Homologous proteins from other organisms were retrieved by links from the CAZY website (please see Availability &amp; requirements for more information) or by BLASTp searches at NCBI. Putative signal peptides were predicted with the SignalP program (please see Availability &amp; requirements for more information) [##REF##15223320##58##];  and putative N-glycosylation sites were identified with the NetNGlyc [##UREF##13##59##]; (please see Availability &amp; requirements for more information) and manually inspected to remove NPS/T sites.</p> <p>Chromosomal locations of genes were identified on and the map drawn with the NCBI map viewer (please see Availability &amp; requirements for more information). Segmental duplication analysis was done with DAGchainer [##REF##15247098##60##] and the TIGR rice segmental duplication database (please see Availability &amp; requirements for more information) with the maximum length distance permitted between collinear gene pairs set to be 500 kb. LTR-retrotransposon elements, interspersed repeats and low complexity DNA sequences were identified with RepeatMasker (please see Availability &amp; requirements for more information).</p> <title>Construction of phylogenetic trees</title> <p>The multiple sequence alignment of β-galactosidase protein sequences from rice and other organisms was made with ClustalW [##REF##7984417##61##] and manually adjusted and edited to remove unconserved N- and C-terminal regions with Genedoc (please see Availability &amp; requirements for more information) [##UREF##14##62##]. Phylogenetic trees were constructed by the neighbor-joining facility in ClustalX 1.83 [##REF##9810230##63##] and the branching order verified for maximum parsimony with the Protpars program in the PHYLIP software suite [##UREF##15##64##].</p> <title>Plant material</title> <p>Indica rice (<italic>Oryza sativa </italic>L., cv. KDML 105) seeds were soaked under sterile conditions on tissue paper moistened with sterile distilled water at 28°C until germination. Germinated rice was moved to soil-filled plastic pots and grown until 15 days, then moved to large clay pots. Samples were collected every 3 days from seeding until 15 days after seeding, then sampling was continued every month. To attain more materials at the flowering stage, samples were collected from KDML105 fields at the Pathum thani Rice Research Center in November, 2005. The samples were immediately frozen in liquid nitrogen and stored at -80°C until use.</p> <title>Semi-quantitative reverse transcription-polymerase chain reaction (qRT-PCR)</title> <p>Total RNA was isolated from various tissues: germinated seed, root, shoot, leaf blade, leaf sheath, node, internode, initiating panicle, developing panicle, emerging panicle, anther, flower, milk grain, and grain during dry down by the following procedure. Tissues (100 mg) were ground to powder under liquid nitrogen, then RNA was extracted using RNeasy Plant Mini (QIAGEN GmbH, Hilden, Germany) or Sigma Spectrum RNA extraction (Sigma-Aldrich, St Louis, USA) kits. For starchy tissues: germinated seed and grain during dry down, RNA extractions were first done with Trizol reagent (Invitrogen, Carlsbad, CA, USA [##REF##2440339##65##]), and the RNA further purified from the extracts by the spin column procedure of the above RNA extraction kits. The RNA was quantified in a spectrophotometer at 260 nm.</p> <p>For RT-PCR, approximately 5 μg of total RNA was treated with RQ1 RNase-free DNase I (Promega Corporation, Madison, WI, USA). First stand cDNA was synthesized from RNA template primed with Oligo (dT)<sub>20 </sub>with the SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen). The gene-specific primers used for semi-quantitative PCR were designed from the 3'UTR or 3' coding sequence of each β-galactosidase gene (Table ##TAB##4##5##). The constitutively expressed <italic>β-Actin </italic>and <italic>Ubiquitin-6 </italic>genes were used for normalization. The PCR (10 μl total volume) was done with a 10 ng or 100 ng aliquot of the first stand cDNA as template and 0.5 units of <italic>Taq </italic>DNA polymerase (Promega). The relative expression abundances were obtained by dividing the total densimetric intensities measured on a Fluor-S gel documentation system with Quantity 1 software (Bio-RAD, Hercules, CA, USA) for each gene by those for the control gene. All samples were assayed 3 times in separate reactions to give means and standard deviations for the relative abundances.</p> <title>Amplification of cDNA clones and sequencing</title> <p>The first strand cDNA reverse-transcribed from RNA of developing and emerging panicle was used as a template to amplify <italic>OsBgal6</italic>, <italic>OsBgal11</italic>, and <italic>OsBgal15 </italic>cDNA. The set of primers for amplification of partial or full-length cDNA of each gene was designed from its genomic sequence. The PCR was catalyzed by <italic>Pfu </italic>polymerase (Promega) with heating at 94°C for 5 min, followed by 30 cycles of 1 min at 94°C, 30 s at the appropriate annealing temperature, (58°C for most primers) and 5 min at 72°C, and a final extension for 7 min at 72°C. The final product was gel purified with the QIAquick Gel Extraction Kit (QIAGEN), then sequenced by automated sequencing at Macrogen Inc. (Seoul, Korea).</p> <p>Two pairs of primers were used to amplify <italic>OsBgal6 </italic>cDNA: (1), bgal6_startF (TCAGTCAGTAGTCAGACATG) and bgal6_QIENEY_R (AATGCAGGCTCAATCATCAG); and (2), bgal6_QIENEY_F (GATGATTGAGCCTGCATTTG) and bgal6_stopR (AGTTTCCTGTGTTGCATCAC). The full-length sequence of <italic>OsBgal6 </italic>was determined with the above primer sets and additional primers, including bgal6_1163r (GGTGTGTTATGCTGATCGAAG), bgal6_1650f (GGATTCTGGCGCCTACATG), bgal6_1143F (CTTCGATCAGCATAACACACC), and bgal6_1671r (CCATGTAGGCG CCAGAATC).</p> <p>The oligonucleotide primers used for amplification of a partial <italic>OsBgal11 </italic>cDNA were: bgal11_Seq_r1 (TCATCGCATGTGCAGTG) and bgal11_STOP_R (CCTTCTTCCTAAGCCGCCTG).</p> <p>The <italic>Osbgal15 </italic>cDNA was amplified with 2 pairs of primers. First, the bgal15_START_F (CGCGTGCCGGCGATGAAG) and bgal15_QIENE_Rev (CTCGTTTTCAATCTGTGCCAG) primers were used to amplify a 5' cDNA fragment, and the bgal15_QIENEY_F (CTGGCACAGATTGAAAACGAG) and bgal15_STOP_R (TATCAACATGAAGCCTGAACGGTG) primers were used to amplify an overlapping 3' cDNA fragment. The products of these first two PCR were mixed and a sequential PCR was performed with the bgal15_START_F and bgal15_STOP_R primers. The full-length CDS cDNA was cloned into the <italic>EcoR</italic>Isite of pBluescript KS(+) by standard methods [##UREF##16##66##]. The recombinant plasmid was sequenced with the T7, M13, bgal15QIENE_F, bgal15QIENE_REV, bgal15_1138_F (ACTCATCTTCTGCCTGCTTG), and Bgal15_1797_R (TACGGTGCCATTGTTG TTG) primers.</p> <title>Expression of OsBgal13</title> <p>A cDNA encoding the predicted mature OsBgal13 protein was amplified from the full-length cDNA clone of Genbank Accession <ext-link ext-link-type="gen" xlink:href="AK065546">AK065546</ext-link>[##REF##12869764##67##], obtained from the Rice Genome Resource Center, Tsukuba, Japan. The plasmid clone was used as template for amplification with the bgal13matN (CACCGCCGCCGCGTATGC) and bgal13stopR (GATCACATCTCACCGCGAGGCTC) primers and <italic>Pfu </italic>polymerase (Promega), according the supplier's instructions. The PCR product was gel-purified and cloned into the pENTR/D-TOPO plasmid (Invitrogen), according to the supplier's instructions, and sequenced. The cDNA was transferred into the pET32/DEST vector [##REF##17196101##68##] by an LR clonase reaction (Invitrogen). The thioredoxin fusion protein was expressed from this plasmid and expressed as described for OsBgal1 by Chantarangsee et al. [##UREF##9##38##]. The protein was extracted and purified on Ni-NTA IMAC resin (QIAGEN), as with OsBgal1 [##UREF##9##38##], except that 0.1 mg/mL soybean trypsin inhibitor was added to the extraction buffer, and the protein was eluted from the IMAC column with 5–20 mM imidazole. The OsBgal13 protein was further purified by gel-filtration on a Superdex S200 column. The purified OsBgal13 was tested for glycone specificity by incubation of enzyme with 1 mM <italic>p</italic>-nitrophenyl (<italic>p</italic>NP) glycosides for 30 min at 30°C, followed by stopping the reaction with 2 volumes 0.4 M Na<sub>2</sub>CO<sub>3 </sub>and measuring the 405 nm absorbance on a microtiter plate. Action on galactose-containing oligosaccharides and polysaccharides was determined by incubating with 1 mM oligosaccharide or 0.5% polysaccharide at 30°C 24 h, then separating the products on TLC and detecting carbohydrates, as previously described [##UREF##9##38##].</p>
<title>Results and Discussion</title> <title>Identification of rice GH35 genes and their protein products</title> <p>A total of 15 GH35 genes (defined <italic>OsBgal</italic>) were identified in rice genome databases (Table ##TAB##0##1##), each of which encodes a protein that contains the GH35 β-galactosidase active site consensus sequence G-G-P- [LIVM]-x-Q-x-E-N-E- [FY] [##REF##9649738##39##]. Most of the encoded proteins contain a GH35 domain at the N-terminus and a Gal-lectin-like domain at the C-terminus. Although the speculated carbohydrate-binding function of these Gal-lectin domains is not yet proven, their existence and conserved nature has led to the suggestion that they may increase the catalytic efficiency on polysaccharide substrates [##REF##17466346##12##]. The proteins from <italic>OsBgal</italic>2, <italic>OsBgal</italic>7 and <italic>OsBgal9 </italic>lack the Gal-lectin-like domain, suggesting that it is not necessary for the functions of these rice GH35 members.</p> <p>The deduced amino acid sequences of the OsBgal proteins were used to predict their putative signal peptides, protein lengths, molecular masses, pI values, possible N-glycosylation sites, and cellular destinations (Table ##TAB##1##2##). All OsBgals contain putative signal peptides, which range in length from 20 to 36 amino acids, except for that of OsBgal8 which is predicted to have 62 amino acids. The mature OsBgal proteins were predicted to contain from 653 (OsBgal9) to 894 (OsBgal8) amino acids, corresponding to molecular masses of 73.5 to 97.9 kDa. Most of the proteins were predicted to have pI in the acidic range (5.56–6.76), except those of OsBgal2, OsBgal7 and OsBgal10 are in the basic range (7.65–9.1) (Table ##TAB##1##2##). Only OsBgal2 and OsBgal7 lacked putative N-glycosylation sites, while the other enzymes were predicted to contain from 2 to 10 sites. Most OsBgal members are predicted to localize to the organelles of the secretory pathway, for instance, the Golgi apparatus, endoplasmic reticulum, or vacuole, or to be secreted. Only 3 OsBgals, including OsBgal9, OsBgal10 and OsBgal11, were predicted to be localized in the lysosome-like vacuole, which correlates to the location of mammalian β-galactosidases in the lysosome (Table ##TAB##1##2##). Although OsBgal10 and OsBgal11 are very similar in predicted mature protein length and possible cellular destinations, their predicted pI values are quite different. OsBgal5, OsBgal12, OsBgal14 and OsBgal15 all have high numbers of glycosylation sites, similar MW, pI and predicted possible destinations, which suggests that they may have redundant or overlapping functions.</p> <title>Phylogenetic analysis</title> <p>The multiple alignment of full-length protein sequences was used to construct an unrooted phylogenetic tree. The phylogenetic tree, which includes all the GH35 genes identified in the rice (monocot), <italic>A. thaliana </italic>(dicot) and <italic>Physcomitrella patens </italic>(a bryophyte moss) genomes, as well as representatives of animals, fungi, protists, archaea, and eubacteria, has three major branches, one of which is nearly plant specific (except for one <italic>Dictyostelium discoideum </italic>gene product), one of which contains representatives of animals, eubacteria and plants, and one of which contains archaeal, eubacterial and fungal enzymes (Figure ##FIG##0##1##). Almost all the rice, <italic>A. thaliana</italic>, and <italic>P. patens </italic>enzymes fall within the plant-specific cluster, except <italic>OsBgal9</italic>, <italic>AtBGAL</italic>17 and <italic>P. patens </italic>EDQ62875, which fall in the cluster with animal β-galactosidases. Previously, it was suggested that the plant enzymes that cluster with animals might have been transferred to plants by horizontal gene transfer [##REF##17466346##12##], but the fact that the bryophyte, which is thought to have diverged from the vascular plants early in plant evolution, has this type of β-galactosidase and the broad distribution of organisms with this type of β-galactosidase suggest that plants maintained a copy of this gene when plants and animals diverged. Perhaps the more relevant question is how the plant specific Bgal ancestor came to be, and how the C-terminal Gal-lectin-like domain was acquired by this lineage. It appears that the fungal-type and plant-type lineages may have diverged from the animal-type β-galactosidases early in GH35 evolution, but plants retained the animal-type β-galactosidase, as well.</p> <p>Within the plant-type β-galactosidases, OsBgals can be divided into three distinct groups (a, b and c in Figure ##FIG##0##1##). The major group (group a) contains 5 clusters with total 7 OsBgal members. The cluster a1 contains 4 rice genes, including OsBgal1, OsBgal2, OsBgal3, and OsBgal7, along with 6 AtBGALs, the <italic>Asparagus officinalis </italic>Bgal, and 4 tomato and 5 chickpea Bgals that are not shown, some of which contain C-terminal galactose-binding lectin-like domains and some of which do not. Another cluster, a2, contains OsBgal13, AtBGAL9 and two putative <italic>P. patens </italic>Bgals (EDQ59880 and EDQ81397), which suggests this cluster is of ancient origin. Two other putative <italic>P. patens </italic>Bgals fall in cluster a3, which appears to be closely related to cluster a2. OsBgal8 groups with AtBGAL8 and RsBgal1 of radish, which specifically hydrolyzes the β-(1→3)- and β-(1→6)-linked oligogalactans of AGPs [##REF##15980190##18##], in subsgroup a4. OsBgal4 groups with AtBGAL10 in cluster a5. The second biggest group, b, contains 4 rice Bgals with high amino acid sequence similarity, namely OsBgal5, OsBgal12, OsBgal14, and OsBgal15. These isozymes have nearly the same protein lengths and gene structures, as described below. Within group c, OsBgal10 and OsBgal11 (cluster c1) are closely related, while OsBgal6 (cluster c2) is less clearly associated, though both distance-based and maximum parsimony trees back this model.</p> <title>Chromosomal locations of OsBgals</title> <p>As shown in the map of their chromosomal locations and directions of transcription in Figure ##FIG##1##2##, the <italic>OsBgals </italic>are distributed over all chromosomes, except chromosomes 4, 7 and 11. Three <italic>OsBgals </italic>are located on chromosome 1, while 2 <italic>OsBgal </italic>each are found on chromosomes 3, 5, 6, and 10, and one <italic>OsBgal </italic>each is present on chromosomes 2, 8, 9, and 12. Similar to other gene families, <italic>OsBgals </italic>appear to have undergone gene duplication, as the rice genome has undergone genome-wide duplication events, including polyploidy, which promote the amplification of gene family members. Segmental duplication analysis identified four <italic>OsBgals </italic>located on the duplicated segmental regions in rice chromosomes. <italic>OsBgal7 </italic>(Rice genome project locus Os02g0219200) and <italic>OsBgal2 </italic>(Os06g0573600) are located on a chromosomal segment duplicated between chromosomes 2 and 6, while <italic>OsBgal10 </italic>(Os08g0549200) and <italic>OsBgal11 </italic>(Os09g0539200) are located on segments duplicated between chromosomes 8 and 9. These duplications are supported by the close phylogenetic relationships of these genes. Therefore, chromosomal segment duplication appears to have played a role in multiplication of <italic>OsBgals </italic>within the rice genome. In contrast, <italic>OsBgal12 </italic>and <italic>OsBgal14</italic>, which have high sequence similarity, are located close together on the same chromosome. Thus, these two genes may have diverged by tandem duplication.</p> <title>Exon-intron organization and primary gene structure analysis</title> <p>Automated annotations for all <italic>OsBgal </italic>genes were available in the public databases, except for <italic>OsBgal6</italic>, <italic>OsBgal14 </italic>and <italic>OsBgal15</italic>. However, we confirmed and corrected exon-intron organization manually by comparing the corresponding full length cDNA and EST sequences, which were either already available in publicly accessible databases or were obtained during this study, with genomic sequences.</p> <p>Gene structure comparisons showed eight different splicing patterns. The sizes of the majority of the coding exons are conserved, but some appear to have had intron-loss events (Figure ##FIG##2##3##). The pattern with the highest number of exons, found in <italic>OsBgal1, OsBgal3, OsBgal4</italic>, <italic>OsBgal6 and OsBgal13</italic>, contains 19 exons with 18 introns, with exons 1 through 10 encoding the GH35 domain, and exons 18 and 19 encoding the Gal-lectin-like domain. Ahn et al [##REF##17466346##12##] observed this same pattern in <italic>A. thaliana Bgals </italic>and surmised from parsimony that it is likely to be the ancestral gene pattern and other patterns were derived from it, mainly by intron loss. The presence of this pattern in diverse rice genes and in a <italic>P. patens </italic>gene (EDQ59880) supports this conclusion. The second most similar pattern is 18 exons separated by 17 introns with the loss of intron 1, which occurred in <italic>OsBgal11</italic>. A pattern of two-intron loss was found in <italic>OsBgal5, OsBgal12</italic>, <italic>OsBgal14</italic>, and <italic>OsBgal15</italic>, which contain 17 exons with the loss of introns 17 and 18.</p> <p>The <italic>OsBgal2</italic>, <italic>OsBgal7 </italic>and <italic>OsBgal9 </italic>genes all lack the galactose-binding, lectin-like domain at the C-terminus, however, their splicing patterns were different. <italic>OsBgal2 </italic>and <italic>OsBgal7 </italic>have patterns similar to other rice GH35 genes and retain the 5' end of exon 18, but <italic>OsBgal2 </italic>contains 16 exons with the loss of introns 11 and 16, while intron 15 has been lost in <italic>OsBgal7</italic>, leaving 17 exons (Figure ##FIG##2##3##). In contrast, <italic>OsBgal9 </italic>appeared to be distantly diverged from other rice GH35 members since its exon/intron organization does not correspond to the other <italic>OsBgals</italic>, except that an intron is found in the same position as intron 1 in the other genes.</p> <p>Although the pattern and lengths of exons in rice GH35 genes are similar, <italic>OsBgal13 </italic>spans approximately 20.8 kb, because it contains over three-fold longer introns than the other genes in the family. Long terminal repeat (LTR)-retrotransposons elements and other repeated sequences were found within <italic>OsBgal13 </italic>gene introns 1, 2, 8, 9, 10, 13, 14 and 16, which contain one GC rice region, three AT rich regions and eight transposable elements (Table ##TAB##2##3##). Likewise, <italic>OsBgal6</italic>, the second largest <italic>OsBgal </italic>gene due to its long (4 kb) intron 1, had one AT rich region and seven transposable elements, all of which were located within intron 1. In <italic>OsBgal </italic>genes of normal length, such as <italic>OsBgal1</italic>, no transposable elements were observed. These results may reflect the observation that rice genome expansion is forced by transposable element amplification [##REF##16819716##40##].</p> <p>Comparison of gene splice patterns and phylogenetic relationships between plant-type Bgals from rice, <italic>A. thaliana</italic>, and <italic>P. patens </italic>reveals that, although the ancestral gene pattern is the same, no other splice patterns are shared within a phylogenetic lineage. Therefore, though most introns appear to have been inserted in the common ancestor gene for this lineage (an extra exon was inserted in <italic>P. patens </italic>gene EDQ878340), intron losses appear to have occurred independently in rice, <italic>Arabidopsis </italic>and <italic>P. patens</italic>. The same may be true for the C-terminal domain, which appears to have been lost only once in rice (in the putative ancestor of <italic>OsBgal2 </italic>and <italic>OsBgal7</italic>), though it appears to have been lost at least 4 times in <italic>Arabidopsis </italic>GH35 <italic>Bgals </italic>[##REF##17466346##12##].</p> <title>Analysis of expression by RT-PCR</title> <p><italic>OsBgal </italic>transcript expression in various tissues and stages of growth was analyzed by RT-PCR with primers specific to the 3'UTR of each gene to gain further insight into possible functions. The relative expression levels derived from normalization with the <italic>UBQ6 </italic>polyubiquitin gene gave the same patterns as those shown for normalization with the <italic>β-Actin </italic>gene in Figure ##FIG##3##4## and Additional file ##SUPPL##0##1##, indicating the patterns likely reflect the expression of the <italic>OsBgal </italic>genes, rather than that of the control.</p> <p>The rice tissues/organs and developmental stages chosen for RT-PCR analysis are shown in Figure ##FIG##3##4##. All 15 <italic>OsBgal </italic>genes were found to be expressed with different but overlapping expression patterns, confirming that all genes are active. These relative expression levels were used to define genes as having expression specific to a given tissue, if the relative transcript levels of the gene at that stage are significantly higher (color scale above 7, where 10 is the level in the tissue with the highest expression of each gene) over the levels at all other stages. The expression profiles revealed that at least 9 genes were highly expressed in one of the stages of vegetative tissues and 10 were highly expressed in reproductive tissues. Only <italic>OsBgal13 </italic>appears to be expressed in nearly every organ and developmental stage, though at different levels (Figure ##FIG##3##4##).</p> <p>Nine <italic>OsBgals </italic>(<italic>OsBgal1, OsBgal2, OsBgal3, OsBgal4, OsBgal6, OsBgal7, OsBgal8 and OsBgal9</italic>) displayed their highest transcript levels in young and early-mature stages of shoot, root, and leaf sheath. Among these genes, <italic>OsBgal6, OsBgal7, OsBgal8 and OsBgal9 </italic>showed highest intensity in leaf sheath at different times, especially, 15-day and 1-month old, which are both fast cell elongation stages. The high accumulation of OsBgals in this stage may act in sugar remodeling, since the carbohydrate remobilized from the leaf sheath and culm to grain can contribute as much as 38% to rice yield [##UREF##10##41##]. In contrast, low levels of most gene transcripts were detected in leaf blade, except for <italic>OsBgal13 </italic>transcripts, which were detected throughout the growth of the blade.</p> <p>Six genes, <italic>OsBgal5, OsBgal10, OsBgal11, OsBgal12, OsBgal14</italic>, and <italic>OsBgal15 </italic>appear to be expressed primarily in reproductive tissues. The relative transcript levels for these genes were at least five times higher in reproductive tissues than in vegetative tissues (Figure ##FIG##3##4##). In addition, no expression was detected at the initiating panicle stage, but it was high in later panicle development and flowering stages. Only <italic>OsBgal10, OsBgal11 </italic>and <italic>OsBgal15 </italic>transcripts were found to accumulate in anther. These genes appear to have diverged from two ancestor genes, since <italic>OsBgal5</italic>, <italic>OsBgal12</italic>, <italic>OsBgal14</italic>, and <italic>OsBgal15 </italic>have sequences and gene structures that are similar to each other, as do <italic>OsBgal10 </italic>and <italic>OsBgal11 </italic>(Figures ##FIG##0##1## and ##FIG##2##3##). Interestingly, <italic>AtBGAL7 </italic>and <italic>AtBGAL15</italic>, which fall into the group b with <italic>OsBgal5</italic>, <italic>OsBgal12</italic>, <italic>OsBgal14</italic>, and <italic>OsBgal15 </italic>(Figure ##FIG##0##1##), were reported to function in the early stages of microspore and pollen development [##REF##15517348##31##]. Likewise, <italic>AtBGAL11 </italic>and <italic>AtBGAL13</italic>, which fall in cluster c1 with <italic>OsBgal10 </italic>and <italic>OsBgal11 </italic>(Figures ##FIG##0##1##), show maximum expression in mature pollen [##REF##15517348##31##]. So, related rice and <italic>Arabidopsis </italic>Bgals function in reproductive tissues, suggesting their roles may be somewhat conserved in the two plants.</p> <p>Only the <italic>OsBgal4 </italic>transcript was expressed abundantly in seed at the imbibition stage, while <italic>OsBgal2 </italic>and <italic>OsBgal13 </italic>were also present at low levels. Five <italic>OsBgals </italic>were expressed in milk grain and six in grain during dry down. A high abundance of <italic>OsBgal2 </italic>and <italic>OsBgal13 </italic>transcripts was observed in grain during dry down. This evidence implies that they may function in grain development and senescence or are stored for roles in germination, such as cell wall remodeling. <italic>OsBgal1</italic>, <italic>OsBgal2, OsBgal6, OsBgal8</italic>, and <italic>OsBgal13 </italic>transcripts were also found in shoots and roots of seedlings after germination. Chantarangsee et al. [##UREF##9##38##] reported β-galactosidase activity in rice seeds, roots and shoots at 0–7 days after seed-soaking, and detected OsBgal1 and OsBgal2 proteins in embryos and seedling roots and shoots. In barley (<italic>Hordeum vulgare</italic>), Giannakouros and colleagues [##UREF##2##17##] were able to separate 4 isozymes of β-galactosidase from germinated grain. Thus, Bgals are present and may play roles in the growth and emergence of root and shoot primordia from the hull.</p> <p>To complement the RT-PCR expression profiles of the <italic>OsBals</italic>, ESTs retrieved from the UniGene database (UniGene December, 2007, NCBI) were compared with our findings. Full-length or partial cDNA and ESTs for most <italic>OsBgals </italic>are available in the database, except <italic>OsBgal6</italic>, for which only one EST from unspecified tissues is available, and <italic>OsBgal15</italic>, for which no corresponding ESTs or cDNA are available (Table ##TAB##0##1##). However, the other <italic>OsBgal </italic>genes had cDNA and ESTs present in the database, with numbers of ESTs that varied from 1 to 253.</p> <p>The <italic>OsBgal8 </italic>(Unigene cluster Os.22360, locus <italic>Os03g0255100</italic>) has the highest number of ESTs in the database: 253 ESTs from various tissues, including flower, stem, leaf, panicle, root, callus, seed, whole plant and unspecified tissues. The Unigene profiles indicate relatively abundant accumulation in flower, leaf, panicle, root, seed and stem with 761, 133, 120, 190, 216 and 228 transcripts per million (TPM), respectively (Table ##TAB##0##1##). The Unigene profile is somewhat similar to the RT-PCR pattern of <italic>OsBgal8</italic>, in which the transcripts accumulated significantly throughout rice plant, including in seed during dry-down, but not in seedling and milk grain.</p> <p><italic>OsBgal10 </italic>(Os.18310), <italic>OsBgal11</italic>, <italic>OsBgal12 </italic>(Os.46702) and <italic>OsBgal14 </italic>(Os.22528) expression data included mainly flower and panicle clones, or had relatively high expression in these tissues, which correlates to our RT-PCR analysis. For instance,<italic>OsBgal11 </italic>shows 249 TPM in flower and 158 TPM in panicle, while <italic>OsBgal14 </italic>shows 358 TPM in flower (Table ##TAB##0##1##).</p> <p>The EST and cDNA clones for the remaining <italic>OsBgals </italic>had diverse origins, which is consistent with the evaluation of <italic>OsBgals </italic>by RT-PCR, though the frequencies of expression derived from these two forms of expression data were not exactly the same. Therefore, Unigene database analysis could be used for an initial guideline of gene expression, though the expression patterns based on ESTs are fragmentary, since the tissues are not equally or completely represented and the source tissue annotation of many EST libraries is unclear or incomplete. Microarray experiments, such as that by Ma et al. [##REF##16140994##42##], who assessed expression of 37,132 non cross-hybridizing gene models in rice seedling shoots, tillering-stage shoots and roots, heading and filling stages of panicles and suspension cells, offer a more comprehensive picture of gene expression within a tissue. Although the tissue coverage in public microarray databases is as yet limited, the GEO expression data in Table ##TAB##0##1## shows that OsBgal expression levels may be affected by many treatments.</p> <title>Rice β-galactosidase cDNAs sequencing and identification of alternative splicing</title> <p>Although the expression of all <italic>OsBgals </italic>can be observed by RT-PCR and most have corresponding cDNAs and ESTs in the database, complete full-length cDNA clones were not available for <italic>OsBgal6</italic>, <italic>OsBgal11 </italic>and <italic>OsBgal15</italic>. Therefore, to confirm the mRNA sequence and gene structure, these genes were amplified by PCR and sequenced. The accession numbers for these sequences are listed in Table ##TAB##0##1##.</p> <p>The full-length cDNA of <italic>OsBgal</italic>6 was amplified by a set of primers corresponding to the transcribed locus <ext-link ext-link-type="gen" xlink:href="XM_475258">XM_475258</ext-link>. The cDNA was isolated from 15-day leaf sheaths, a tissue with abundant expression, and exerting panicles, a tissue with poor expression, based on RT-PCR expression analysis. Sequencing of the cDNA from the highly expressing tissue revealed an open reading frame (ORF) of 2,436 nucleotides encoding a protein of 812 amino acids (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU602310">EU602310</ext-link>), which confirmed the our predicted β-<italic>galactosidase </italic>sequence from chromosome 5 (Genbank accession no. <ext-link ext-link-type="gen" xlink:href="AC135419">AC135419</ext-link> and <ext-link ext-link-type="gen" xlink:href="AC135429">AC135429</ext-link>). A search through the dbEST database of <italic>O. sativa </italic>found 99% identity with a 520 bp EST from fertile panicle (<ext-link ext-link-type="gen" xlink:href="CK048087">CK048087</ext-link>). On the contrary, the Os <italic>Bgal</italic>6 cDNA amplified from exerting panicle, where expression is low, contains only 2,338 bp (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU603285">EU603285</ext-link>) which is 101 bp shorter than the cDNA from high expression tissues. The open reading frame beginning at the initiation codon of this transcript covers only exons 1 to 7 and encodes only 244 amino acids. Alternative splicing of this transcript occurred to utilize an alternative acceptor site (CAG) 80 nt downstream of functional acceptor site for intron 7, which introduced a deletion of 27 amino acid and a frameshift that leads to an in-frame stop codon (Figure ##FIG##4##5##). A second missplicing of this gene occurred at the intron 10/exon 11 junction, which again used an alternate acceptor site (CAG) located 20 nt downstream of the regular acceptor site. So, alternative splicing is likely to be one of the gene regulation mechanisms, since the transcript would be subject to nonsense mediated decay (nmd) [##REF##17194304##43##] and the protein produced would be useless as a β-galactosidase.</p> <p><italic>OsBgal11 </italic>(locus Os09g0539200) is annotated as a pseudogene in GenBank, because of its apparent inability to produce full-length protein due to a frameshift and subsequent stop codon in the AK119414 cDNA. We hypothesized that such phenomenon probably resulted from missplicing, as seen in <italic>OsBgal6</italic>. A 1343-nucleotide fragment from the 3' half of the <italic>OsBgal11 </italic>cDNA was amplified from exerting panicle (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU603286">EU603286</ext-link>), in which its expression is high based on RT-PCR. This cDNA shares 98% identity with the AK119414 cDNA, but had a substitution of 75 nucleotides in place of 95 nucleotides at positions 1282 to 1377 of AK119414. These substituted sequences allow the full-length translation of the predicted 839 amino acids of OsBgal11, whereas only 426 amino acids are translated from AK119414. Comparison with the genomic sequence and our amplified cDNA sequence indicates that the AK119414 transcript results from retention of intron 10, which extends exon 10 by 101 nt, thereby introducing an in-frame stop codon. In addition, an alternative donor site located at the functional acceptor site of intron 10 and the first nucleotide of the functional exon 11 produces an alternate intron 10 with the acceptor site located 79 nt downstream, within the functional exon 11, thereby reducing the size of exon 11. The shortened (456-residue) protein, includes all of the catalytic domain, based on the <italic>Penicillium </italic>sp. β-galactosidase structure [##REF##15491613##44##], so it is possible that it could still function as a Bgal.</p> <p>No experimentally determined cDNA or ESTs corresponding to <italic>OsBgal15 </italic>were available in public databases, despite the fact that our RT-PCR demonstrated its expression. Therefore, the full-length <italic>OsBgal15 </italic>cDNA was amplified, cloned and sequenced. The <italic>OsBgal15 </italic>cDNA sequence (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU051629">EU051629</ext-link>) contains a 2487-nucleotide ORF, which encodes 829 amino acids. The encoded protein shares over 99% identity with a predicted β-galactosidase of <italic>Oryza sativa </italic>(japonica cultivar-group) (accession no. <ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link>), except that 7 extra amino acid residues were inserted at residue 293, compared to the predicted sequence from <ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link>. Therefore, the correct splice donor site of exon 8 is located 21 nucleotides downstream of that predicted for the mRNA of <ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link>.</p> <p>The alternative splicing observed in this study introduces in-frame stop codons into the mRNA of <italic>OsBgal6 </italic>in the low expression tissue, but not for the high expression tissues, while it can also introduce a premature stop codon in <italic>OsBgal11</italic>, based on the AK119414 cDNA sequence. In vertebrates, there are a number of well-known alternatively splicing genes, for instance the fibronectin gene, which can generate over 20 different proteins, some of which have different patterns of localization and slightly different functions in human cells [##REF##7226224##45##]. Until relatively recently, there were very few examples of alternatively spliced genes in higher plants and still fewer examples known to generate functionally distinct proteins. Generally, alternative splicing yields two polypeptides of different size that are identical, except for the presence of a number of additional amino acids at the C-terminus of the larger isoform, eg., ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) activase [##REF##10430961##46##]; diacylglycerol kinase in tomato [##REF##11069705##47##]; and chloroplast ascorbate peroxidase in spinach [##REF##9931296##48##]. Missplicing events were observed in maize (<italic>Zea mays</italic>) <italic>Viviparous 1 </italic>(<italic>Vp1</italic>) transcripts, the majority of which are spliced incorrectly and contain insertions of intron sequences or deletions of coding region, so that they do not encode full-length proteins [##REF##12119408##49##]. Although alternative splicing is a widespread process used in higher eukaryotes to regulate gene expression and functional diversification of proteins [##REF##11173120##50##,##REF##12600935##51##], the putative control of tissue-specific expression by missplicing we observed has not been well characterized in rice.</p> <title>Preliminary Characterization of OsBgal13</title> <p>OsBgal13 is in cluster a2 of group a of the phyloenetic tree in Figure ##FIG##0##1##, along with AtBGAL9 and two <italic>P. patens </italic>Bgals. This and its general expression pattern suggest it may have a conserved maintenance function in plants. To assess what its function might be, OsBgal13 was expressed as a thioredoxin fusion protein in <italic>Escherichia coli </italic>strain OrigamiB (DE3), which we previously found to have little background β-galactosidase acitivity [##UREF##9##38##]. The protein was purified by immobilized metal affinity chromatography on nickel resin, and most of the OsBgal13 activity eluted in 5–20 mM imidazole. Further purification was achieved by gel filtration chromatography and the activity of the resulting protein was assessed. As shown in Table ##TAB##3##4##, among <italic>p</italic>-nitrophenol (<italic>p</italic>NP) glycosides, the OsBgal13 enzyme had highest activity toward <italic>p</italic>NP-β-D-galactoside (<italic>p</italic>NPGal), but also had 30% of this activity toward <italic>p</italic>NP-α-L-arabinoside (<italic>p</italic>NPAra), and low activity toward <italic>p</italic>NP-β-D-fucoside and <italic>p</italic>NP-β-D-mannoside, but no detectable activity toward <italic>p</italic>NP-β-D-glucoside and <italic>p</italic>NP-β-D-xyloside and negligible activity toward <italic>p</italic>NP-α-D-galactoside. OsBgal13 also hydrolyzed β-(1→3)-, β-(1→4)- and β-(1→6)-linked galactobiose and galactotriose, but no release of galactose or arabinose could be detected when the protein was incubated for 24 h with rice seedling alcohol insoluble residue (AIR), larchwood arabinogalactan, galactan, apple pectin, citrus pectin or oat xylan. So, despite the intriguing observation that OsBgal13 can hydrolyze both β-D-galactosyl residues and α-L-arabinosyl residues, components of type I and II arabinogalactans and xyloglucan and glucuronoarabinoxylan side chains [##REF##15012297##52##], only short galacto-oligosaccharides can be identified as possible natural substrates to date.</p> <title>Possible Functions of Other Isozymes</title> <p>The activity of mammal β-galactosidase specifically releases terminal β-galactosyl residues from glycosaminoglycans and the glycolipid GM1 ganglioside. Similarly, the digestion of the galactolipid monogalactosyldiacylglyerol by vacuolar and chloroplast β-galactosidases was reported in wheat [##REF##16663831##53##]. Likewise, β-galacosidase of tomato was reported to act in galactolipid turnover and degradation, which occurs in chloroplasts and chromoplasts during tomato fruit development [##UREF##11##54##,##UREF##12##55##]. However, relatively little is known about metabolism of other glycolipids, glycopoteins and proteoglycans involving plant β-galactosidases, though it is well characterized in mammals. The fact that OsBgal9 is similar to animal β-galactosidases and is predicted to localize to a lysosome-like vacuole, suggests it may play a similar role in rice.</p> <p>Among the plant-like OsBgals, OsBgal10 and OsBgal11 are also predicted to localize to lysosome-like vacuoles, so they could also play a role similar to animal-type β-galactosidases. However, these closely related isozymes have reproductive tissue specific expression, as are the isozymes of cluster b, OsBgal5, OsBgal12, OsBgal14 and OsBgal15. Similarly, AtBGAL11 and AtBGAL13, which clustered in group c1 with OsBgal10 and OsBgal11, and AtBGAL7 and AtBGAL15, which clustered in group b, were found to be expressed in flower [##REF##17466346##12##] and pollen [##REF##15517348##31##]. It may be that the ancestral genes for these enzymes developed reproductive-tissue specific roles before the ancestors of monocots and dicots diverged. If these functions are conserved, these isozymes may have similar roles in the two plants, though these functions remain to be determined.</p> <p>Among the other plant-like genes, it is likely many of them act in cell wall metabolism, but these roles are yet to be discerned. OsBgal1 was shown to hydrolyze β-(1→3)-, β-(1→4)- and β-(1→6)-linked oligosaccharides and, at a low level, arabinogalactan and be expressed in a range of vegetative and reproductive tissues [##UREF##9##38##], so it may have a general role on similar substrates in the cell wall. OsBgal2, OsBgal3 and OsBgal7 are closely related to it, as are 4 tomato Bgals involved in fruit ripening [##REF##10889266##13##] and an asparagus Bgal thought to act in senescence [##REF##12228457##56##], so a cell wall function seems likely for these other cluster c1 isozymes as well, though such a role must be proven. Kaneko and Kobayashi [##REF##12723614##37##] showed that OsBgal8, which is secreted from rice tissue culture cells, could release galactose from galactoxyloglucans and pectic galactans. So, a role in cell wall metabolism is quite likely for this group a isozyme.</p> <p>The role of OsBgal6 remains obscure, as no closely related enzyme has been characterized, however it seems to be specific to vegetative tissues of the young plant. The fact that the transcript found in panicle was misspliced to prevent production of active enzyme is intriguing, but the functional implications of the alternative splicing of both OsBgal6 and OsBgal11 remain to be clarified.</p>
<title>Results and Discussion</title> <title>Identification of rice GH35 genes and their protein products</title> <p>A total of 15 GH35 genes (defined <italic>OsBgal</italic>) were identified in rice genome databases (Table ##TAB##0##1##), each of which encodes a protein that contains the GH35 β-galactosidase active site consensus sequence G-G-P- [LIVM]-x-Q-x-E-N-E- [FY] [##REF##9649738##39##]. Most of the encoded proteins contain a GH35 domain at the N-terminus and a Gal-lectin-like domain at the C-terminus. Although the speculated carbohydrate-binding function of these Gal-lectin domains is not yet proven, their existence and conserved nature has led to the suggestion that they may increase the catalytic efficiency on polysaccharide substrates [##REF##17466346##12##]. The proteins from <italic>OsBgal</italic>2, <italic>OsBgal</italic>7 and <italic>OsBgal9 </italic>lack the Gal-lectin-like domain, suggesting that it is not necessary for the functions of these rice GH35 members.</p> <p>The deduced amino acid sequences of the OsBgal proteins were used to predict their putative signal peptides, protein lengths, molecular masses, pI values, possible N-glycosylation sites, and cellular destinations (Table ##TAB##1##2##). All OsBgals contain putative signal peptides, which range in length from 20 to 36 amino acids, except for that of OsBgal8 which is predicted to have 62 amino acids. The mature OsBgal proteins were predicted to contain from 653 (OsBgal9) to 894 (OsBgal8) amino acids, corresponding to molecular masses of 73.5 to 97.9 kDa. Most of the proteins were predicted to have pI in the acidic range (5.56–6.76), except those of OsBgal2, OsBgal7 and OsBgal10 are in the basic range (7.65–9.1) (Table ##TAB##1##2##). Only OsBgal2 and OsBgal7 lacked putative N-glycosylation sites, while the other enzymes were predicted to contain from 2 to 10 sites. Most OsBgal members are predicted to localize to the organelles of the secretory pathway, for instance, the Golgi apparatus, endoplasmic reticulum, or vacuole, or to be secreted. Only 3 OsBgals, including OsBgal9, OsBgal10 and OsBgal11, were predicted to be localized in the lysosome-like vacuole, which correlates to the location of mammalian β-galactosidases in the lysosome (Table ##TAB##1##2##). Although OsBgal10 and OsBgal11 are very similar in predicted mature protein length and possible cellular destinations, their predicted pI values are quite different. OsBgal5, OsBgal12, OsBgal14 and OsBgal15 all have high numbers of glycosylation sites, similar MW, pI and predicted possible destinations, which suggests that they may have redundant or overlapping functions.</p> <title>Phylogenetic analysis</title> <p>The multiple alignment of full-length protein sequences was used to construct an unrooted phylogenetic tree. The phylogenetic tree, which includes all the GH35 genes identified in the rice (monocot), <italic>A. thaliana </italic>(dicot) and <italic>Physcomitrella patens </italic>(a bryophyte moss) genomes, as well as representatives of animals, fungi, protists, archaea, and eubacteria, has three major branches, one of which is nearly plant specific (except for one <italic>Dictyostelium discoideum </italic>gene product), one of which contains representatives of animals, eubacteria and plants, and one of which contains archaeal, eubacterial and fungal enzymes (Figure ##FIG##0##1##). Almost all the rice, <italic>A. thaliana</italic>, and <italic>P. patens </italic>enzymes fall within the plant-specific cluster, except <italic>OsBgal9</italic>, <italic>AtBGAL</italic>17 and <italic>P. patens </italic>EDQ62875, which fall in the cluster with animal β-galactosidases. Previously, it was suggested that the plant enzymes that cluster with animals might have been transferred to plants by horizontal gene transfer [##REF##17466346##12##], but the fact that the bryophyte, which is thought to have diverged from the vascular plants early in plant evolution, has this type of β-galactosidase and the broad distribution of organisms with this type of β-galactosidase suggest that plants maintained a copy of this gene when plants and animals diverged. Perhaps the more relevant question is how the plant specific Bgal ancestor came to be, and how the C-terminal Gal-lectin-like domain was acquired by this lineage. It appears that the fungal-type and plant-type lineages may have diverged from the animal-type β-galactosidases early in GH35 evolution, but plants retained the animal-type β-galactosidase, as well.</p> <p>Within the plant-type β-galactosidases, OsBgals can be divided into three distinct groups (a, b and c in Figure ##FIG##0##1##). The major group (group a) contains 5 clusters with total 7 OsBgal members. The cluster a1 contains 4 rice genes, including OsBgal1, OsBgal2, OsBgal3, and OsBgal7, along with 6 AtBGALs, the <italic>Asparagus officinalis </italic>Bgal, and 4 tomato and 5 chickpea Bgals that are not shown, some of which contain C-terminal galactose-binding lectin-like domains and some of which do not. Another cluster, a2, contains OsBgal13, AtBGAL9 and two putative <italic>P. patens </italic>Bgals (EDQ59880 and EDQ81397), which suggests this cluster is of ancient origin. Two other putative <italic>P. patens </italic>Bgals fall in cluster a3, which appears to be closely related to cluster a2. OsBgal8 groups with AtBGAL8 and RsBgal1 of radish, which specifically hydrolyzes the β-(1→3)- and β-(1→6)-linked oligogalactans of AGPs [##REF##15980190##18##], in subsgroup a4. OsBgal4 groups with AtBGAL10 in cluster a5. The second biggest group, b, contains 4 rice Bgals with high amino acid sequence similarity, namely OsBgal5, OsBgal12, OsBgal14, and OsBgal15. These isozymes have nearly the same protein lengths and gene structures, as described below. Within group c, OsBgal10 and OsBgal11 (cluster c1) are closely related, while OsBgal6 (cluster c2) is less clearly associated, though both distance-based and maximum parsimony trees back this model.</p> <title>Chromosomal locations of OsBgals</title> <p>As shown in the map of their chromosomal locations and directions of transcription in Figure ##FIG##1##2##, the <italic>OsBgals </italic>are distributed over all chromosomes, except chromosomes 4, 7 and 11. Three <italic>OsBgals </italic>are located on chromosome 1, while 2 <italic>OsBgal </italic>each are found on chromosomes 3, 5, 6, and 10, and one <italic>OsBgal </italic>each is present on chromosomes 2, 8, 9, and 12. Similar to other gene families, <italic>OsBgals </italic>appear to have undergone gene duplication, as the rice genome has undergone genome-wide duplication events, including polyploidy, which promote the amplification of gene family members. Segmental duplication analysis identified four <italic>OsBgals </italic>located on the duplicated segmental regions in rice chromosomes. <italic>OsBgal7 </italic>(Rice genome project locus Os02g0219200) and <italic>OsBgal2 </italic>(Os06g0573600) are located on a chromosomal segment duplicated between chromosomes 2 and 6, while <italic>OsBgal10 </italic>(Os08g0549200) and <italic>OsBgal11 </italic>(Os09g0539200) are located on segments duplicated between chromosomes 8 and 9. These duplications are supported by the close phylogenetic relationships of these genes. Therefore, chromosomal segment duplication appears to have played a role in multiplication of <italic>OsBgals </italic>within the rice genome. In contrast, <italic>OsBgal12 </italic>and <italic>OsBgal14</italic>, which have high sequence similarity, are located close together on the same chromosome. Thus, these two genes may have diverged by tandem duplication.</p> <title>Exon-intron organization and primary gene structure analysis</title> <p>Automated annotations for all <italic>OsBgal </italic>genes were available in the public databases, except for <italic>OsBgal6</italic>, <italic>OsBgal14 </italic>and <italic>OsBgal15</italic>. However, we confirmed and corrected exon-intron organization manually by comparing the corresponding full length cDNA and EST sequences, which were either already available in publicly accessible databases or were obtained during this study, with genomic sequences.</p> <p>Gene structure comparisons showed eight different splicing patterns. The sizes of the majority of the coding exons are conserved, but some appear to have had intron-loss events (Figure ##FIG##2##3##). The pattern with the highest number of exons, found in <italic>OsBgal1, OsBgal3, OsBgal4</italic>, <italic>OsBgal6 and OsBgal13</italic>, contains 19 exons with 18 introns, with exons 1 through 10 encoding the GH35 domain, and exons 18 and 19 encoding the Gal-lectin-like domain. Ahn et al [##REF##17466346##12##] observed this same pattern in <italic>A. thaliana Bgals </italic>and surmised from parsimony that it is likely to be the ancestral gene pattern and other patterns were derived from it, mainly by intron loss. The presence of this pattern in diverse rice genes and in a <italic>P. patens </italic>gene (EDQ59880) supports this conclusion. The second most similar pattern is 18 exons separated by 17 introns with the loss of intron 1, which occurred in <italic>OsBgal11</italic>. A pattern of two-intron loss was found in <italic>OsBgal5, OsBgal12</italic>, <italic>OsBgal14</italic>, and <italic>OsBgal15</italic>, which contain 17 exons with the loss of introns 17 and 18.</p> <p>The <italic>OsBgal2</italic>, <italic>OsBgal7 </italic>and <italic>OsBgal9 </italic>genes all lack the galactose-binding, lectin-like domain at the C-terminus, however, their splicing patterns were different. <italic>OsBgal2 </italic>and <italic>OsBgal7 </italic>have patterns similar to other rice GH35 genes and retain the 5' end of exon 18, but <italic>OsBgal2 </italic>contains 16 exons with the loss of introns 11 and 16, while intron 15 has been lost in <italic>OsBgal7</italic>, leaving 17 exons (Figure ##FIG##2##3##). In contrast, <italic>OsBgal9 </italic>appeared to be distantly diverged from other rice GH35 members since its exon/intron organization does not correspond to the other <italic>OsBgals</italic>, except that an intron is found in the same position as intron 1 in the other genes.</p> <p>Although the pattern and lengths of exons in rice GH35 genes are similar, <italic>OsBgal13 </italic>spans approximately 20.8 kb, because it contains over three-fold longer introns than the other genes in the family. Long terminal repeat (LTR)-retrotransposons elements and other repeated sequences were found within <italic>OsBgal13 </italic>gene introns 1, 2, 8, 9, 10, 13, 14 and 16, which contain one GC rice region, three AT rich regions and eight transposable elements (Table ##TAB##2##3##). Likewise, <italic>OsBgal6</italic>, the second largest <italic>OsBgal </italic>gene due to its long (4 kb) intron 1, had one AT rich region and seven transposable elements, all of which were located within intron 1. In <italic>OsBgal </italic>genes of normal length, such as <italic>OsBgal1</italic>, no transposable elements were observed. These results may reflect the observation that rice genome expansion is forced by transposable element amplification [##REF##16819716##40##].</p> <p>Comparison of gene splice patterns and phylogenetic relationships between plant-type Bgals from rice, <italic>A. thaliana</italic>, and <italic>P. patens </italic>reveals that, although the ancestral gene pattern is the same, no other splice patterns are shared within a phylogenetic lineage. Therefore, though most introns appear to have been inserted in the common ancestor gene for this lineage (an extra exon was inserted in <italic>P. patens </italic>gene EDQ878340), intron losses appear to have occurred independently in rice, <italic>Arabidopsis </italic>and <italic>P. patens</italic>. The same may be true for the C-terminal domain, which appears to have been lost only once in rice (in the putative ancestor of <italic>OsBgal2 </italic>and <italic>OsBgal7</italic>), though it appears to have been lost at least 4 times in <italic>Arabidopsis </italic>GH35 <italic>Bgals </italic>[##REF##17466346##12##].</p> <title>Analysis of expression by RT-PCR</title> <p><italic>OsBgal </italic>transcript expression in various tissues and stages of growth was analyzed by RT-PCR with primers specific to the 3'UTR of each gene to gain further insight into possible functions. The relative expression levels derived from normalization with the <italic>UBQ6 </italic>polyubiquitin gene gave the same patterns as those shown for normalization with the <italic>β-Actin </italic>gene in Figure ##FIG##3##4## and Additional file ##SUPPL##0##1##, indicating the patterns likely reflect the expression of the <italic>OsBgal </italic>genes, rather than that of the control.</p> <p>The rice tissues/organs and developmental stages chosen for RT-PCR analysis are shown in Figure ##FIG##3##4##. All 15 <italic>OsBgal </italic>genes were found to be expressed with different but overlapping expression patterns, confirming that all genes are active. These relative expression levels were used to define genes as having expression specific to a given tissue, if the relative transcript levels of the gene at that stage are significantly higher (color scale above 7, where 10 is the level in the tissue with the highest expression of each gene) over the levels at all other stages. The expression profiles revealed that at least 9 genes were highly expressed in one of the stages of vegetative tissues and 10 were highly expressed in reproductive tissues. Only <italic>OsBgal13 </italic>appears to be expressed in nearly every organ and developmental stage, though at different levels (Figure ##FIG##3##4##).</p> <p>Nine <italic>OsBgals </italic>(<italic>OsBgal1, OsBgal2, OsBgal3, OsBgal4, OsBgal6, OsBgal7, OsBgal8 and OsBgal9</italic>) displayed their highest transcript levels in young and early-mature stages of shoot, root, and leaf sheath. Among these genes, <italic>OsBgal6, OsBgal7, OsBgal8 and OsBgal9 </italic>showed highest intensity in leaf sheath at different times, especially, 15-day and 1-month old, which are both fast cell elongation stages. The high accumulation of OsBgals in this stage may act in sugar remodeling, since the carbohydrate remobilized from the leaf sheath and culm to grain can contribute as much as 38% to rice yield [##UREF##10##41##]. In contrast, low levels of most gene transcripts were detected in leaf blade, except for <italic>OsBgal13 </italic>transcripts, which were detected throughout the growth of the blade.</p> <p>Six genes, <italic>OsBgal5, OsBgal10, OsBgal11, OsBgal12, OsBgal14</italic>, and <italic>OsBgal15 </italic>appear to be expressed primarily in reproductive tissues. The relative transcript levels for these genes were at least five times higher in reproductive tissues than in vegetative tissues (Figure ##FIG##3##4##). In addition, no expression was detected at the initiating panicle stage, but it was high in later panicle development and flowering stages. Only <italic>OsBgal10, OsBgal11 </italic>and <italic>OsBgal15 </italic>transcripts were found to accumulate in anther. These genes appear to have diverged from two ancestor genes, since <italic>OsBgal5</italic>, <italic>OsBgal12</italic>, <italic>OsBgal14</italic>, and <italic>OsBgal15 </italic>have sequences and gene structures that are similar to each other, as do <italic>OsBgal10 </italic>and <italic>OsBgal11 </italic>(Figures ##FIG##0##1## and ##FIG##2##3##). Interestingly, <italic>AtBGAL7 </italic>and <italic>AtBGAL15</italic>, which fall into the group b with <italic>OsBgal5</italic>, <italic>OsBgal12</italic>, <italic>OsBgal14</italic>, and <italic>OsBgal15 </italic>(Figure ##FIG##0##1##), were reported to function in the early stages of microspore and pollen development [##REF##15517348##31##]. Likewise, <italic>AtBGAL11 </italic>and <italic>AtBGAL13</italic>, which fall in cluster c1 with <italic>OsBgal10 </italic>and <italic>OsBgal11 </italic>(Figures ##FIG##0##1##), show maximum expression in mature pollen [##REF##15517348##31##]. So, related rice and <italic>Arabidopsis </italic>Bgals function in reproductive tissues, suggesting their roles may be somewhat conserved in the two plants.</p> <p>Only the <italic>OsBgal4 </italic>transcript was expressed abundantly in seed at the imbibition stage, while <italic>OsBgal2 </italic>and <italic>OsBgal13 </italic>were also present at low levels. Five <italic>OsBgals </italic>were expressed in milk grain and six in grain during dry down. A high abundance of <italic>OsBgal2 </italic>and <italic>OsBgal13 </italic>transcripts was observed in grain during dry down. This evidence implies that they may function in grain development and senescence or are stored for roles in germination, such as cell wall remodeling. <italic>OsBgal1</italic>, <italic>OsBgal2, OsBgal6, OsBgal8</italic>, and <italic>OsBgal13 </italic>transcripts were also found in shoots and roots of seedlings after germination. Chantarangsee et al. [##UREF##9##38##] reported β-galactosidase activity in rice seeds, roots and shoots at 0–7 days after seed-soaking, and detected OsBgal1 and OsBgal2 proteins in embryos and seedling roots and shoots. In barley (<italic>Hordeum vulgare</italic>), Giannakouros and colleagues [##UREF##2##17##] were able to separate 4 isozymes of β-galactosidase from germinated grain. Thus, Bgals are present and may play roles in the growth and emergence of root and shoot primordia from the hull.</p> <p>To complement the RT-PCR expression profiles of the <italic>OsBals</italic>, ESTs retrieved from the UniGene database (UniGene December, 2007, NCBI) were compared with our findings. Full-length or partial cDNA and ESTs for most <italic>OsBgals </italic>are available in the database, except <italic>OsBgal6</italic>, for which only one EST from unspecified tissues is available, and <italic>OsBgal15</italic>, for which no corresponding ESTs or cDNA are available (Table ##TAB##0##1##). However, the other <italic>OsBgal </italic>genes had cDNA and ESTs present in the database, with numbers of ESTs that varied from 1 to 253.</p> <p>The <italic>OsBgal8 </italic>(Unigene cluster Os.22360, locus <italic>Os03g0255100</italic>) has the highest number of ESTs in the database: 253 ESTs from various tissues, including flower, stem, leaf, panicle, root, callus, seed, whole plant and unspecified tissues. The Unigene profiles indicate relatively abundant accumulation in flower, leaf, panicle, root, seed and stem with 761, 133, 120, 190, 216 and 228 transcripts per million (TPM), respectively (Table ##TAB##0##1##). The Unigene profile is somewhat similar to the RT-PCR pattern of <italic>OsBgal8</italic>, in which the transcripts accumulated significantly throughout rice plant, including in seed during dry-down, but not in seedling and milk grain.</p> <p><italic>OsBgal10 </italic>(Os.18310), <italic>OsBgal11</italic>, <italic>OsBgal12 </italic>(Os.46702) and <italic>OsBgal14 </italic>(Os.22528) expression data included mainly flower and panicle clones, or had relatively high expression in these tissues, which correlates to our RT-PCR analysis. For instance,<italic>OsBgal11 </italic>shows 249 TPM in flower and 158 TPM in panicle, while <italic>OsBgal14 </italic>shows 358 TPM in flower (Table ##TAB##0##1##).</p> <p>The EST and cDNA clones for the remaining <italic>OsBgals </italic>had diverse origins, which is consistent with the evaluation of <italic>OsBgals </italic>by RT-PCR, though the frequencies of expression derived from these two forms of expression data were not exactly the same. Therefore, Unigene database analysis could be used for an initial guideline of gene expression, though the expression patterns based on ESTs are fragmentary, since the tissues are not equally or completely represented and the source tissue annotation of many EST libraries is unclear or incomplete. Microarray experiments, such as that by Ma et al. [##REF##16140994##42##], who assessed expression of 37,132 non cross-hybridizing gene models in rice seedling shoots, tillering-stage shoots and roots, heading and filling stages of panicles and suspension cells, offer a more comprehensive picture of gene expression within a tissue. Although the tissue coverage in public microarray databases is as yet limited, the GEO expression data in Table ##TAB##0##1## shows that OsBgal expression levels may be affected by many treatments.</p> <title>Rice β-galactosidase cDNAs sequencing and identification of alternative splicing</title> <p>Although the expression of all <italic>OsBgals </italic>can be observed by RT-PCR and most have corresponding cDNAs and ESTs in the database, complete full-length cDNA clones were not available for <italic>OsBgal6</italic>, <italic>OsBgal11 </italic>and <italic>OsBgal15</italic>. Therefore, to confirm the mRNA sequence and gene structure, these genes were amplified by PCR and sequenced. The accession numbers for these sequences are listed in Table ##TAB##0##1##.</p> <p>The full-length cDNA of <italic>OsBgal</italic>6 was amplified by a set of primers corresponding to the transcribed locus <ext-link ext-link-type="gen" xlink:href="XM_475258">XM_475258</ext-link>. The cDNA was isolated from 15-day leaf sheaths, a tissue with abundant expression, and exerting panicles, a tissue with poor expression, based on RT-PCR expression analysis. Sequencing of the cDNA from the highly expressing tissue revealed an open reading frame (ORF) of 2,436 nucleotides encoding a protein of 812 amino acids (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU602310">EU602310</ext-link>), which confirmed the our predicted β-<italic>galactosidase </italic>sequence from chromosome 5 (Genbank accession no. <ext-link ext-link-type="gen" xlink:href="AC135419">AC135419</ext-link> and <ext-link ext-link-type="gen" xlink:href="AC135429">AC135429</ext-link>). A search through the dbEST database of <italic>O. sativa </italic>found 99% identity with a 520 bp EST from fertile panicle (<ext-link ext-link-type="gen" xlink:href="CK048087">CK048087</ext-link>). On the contrary, the Os <italic>Bgal</italic>6 cDNA amplified from exerting panicle, where expression is low, contains only 2,338 bp (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU603285">EU603285</ext-link>) which is 101 bp shorter than the cDNA from high expression tissues. The open reading frame beginning at the initiation codon of this transcript covers only exons 1 to 7 and encodes only 244 amino acids. Alternative splicing of this transcript occurred to utilize an alternative acceptor site (CAG) 80 nt downstream of functional acceptor site for intron 7, which introduced a deletion of 27 amino acid and a frameshift that leads to an in-frame stop codon (Figure ##FIG##4##5##). A second missplicing of this gene occurred at the intron 10/exon 11 junction, which again used an alternate acceptor site (CAG) located 20 nt downstream of the regular acceptor site. So, alternative splicing is likely to be one of the gene regulation mechanisms, since the transcript would be subject to nonsense mediated decay (nmd) [##REF##17194304##43##] and the protein produced would be useless as a β-galactosidase.</p> <p><italic>OsBgal11 </italic>(locus Os09g0539200) is annotated as a pseudogene in GenBank, because of its apparent inability to produce full-length protein due to a frameshift and subsequent stop codon in the AK119414 cDNA. We hypothesized that such phenomenon probably resulted from missplicing, as seen in <italic>OsBgal6</italic>. A 1343-nucleotide fragment from the 3' half of the <italic>OsBgal11 </italic>cDNA was amplified from exerting panicle (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU603286">EU603286</ext-link>), in which its expression is high based on RT-PCR. This cDNA shares 98% identity with the AK119414 cDNA, but had a substitution of 75 nucleotides in place of 95 nucleotides at positions 1282 to 1377 of AK119414. These substituted sequences allow the full-length translation of the predicted 839 amino acids of OsBgal11, whereas only 426 amino acids are translated from AK119414. Comparison with the genomic sequence and our amplified cDNA sequence indicates that the AK119414 transcript results from retention of intron 10, which extends exon 10 by 101 nt, thereby introducing an in-frame stop codon. In addition, an alternative donor site located at the functional acceptor site of intron 10 and the first nucleotide of the functional exon 11 produces an alternate intron 10 with the acceptor site located 79 nt downstream, within the functional exon 11, thereby reducing the size of exon 11. The shortened (456-residue) protein, includes all of the catalytic domain, based on the <italic>Penicillium </italic>sp. β-galactosidase structure [##REF##15491613##44##], so it is possible that it could still function as a Bgal.</p> <p>No experimentally determined cDNA or ESTs corresponding to <italic>OsBgal15 </italic>were available in public databases, despite the fact that our RT-PCR demonstrated its expression. Therefore, the full-length <italic>OsBgal15 </italic>cDNA was amplified, cloned and sequenced. The <italic>OsBgal15 </italic>cDNA sequence (Genbank Accession <ext-link ext-link-type="gen" xlink:href="EU051629">EU051629</ext-link>) contains a 2487-nucleotide ORF, which encodes 829 amino acids. The encoded protein shares over 99% identity with a predicted β-galactosidase of <italic>Oryza sativa </italic>(japonica cultivar-group) (accession no. <ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link>), except that 7 extra amino acid residues were inserted at residue 293, compared to the predicted sequence from <ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link>. Therefore, the correct splice donor site of exon 8 is located 21 nucleotides downstream of that predicted for the mRNA of <ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link>.</p> <p>The alternative splicing observed in this study introduces in-frame stop codons into the mRNA of <italic>OsBgal6 </italic>in the low expression tissue, but not for the high expression tissues, while it can also introduce a premature stop codon in <italic>OsBgal11</italic>, based on the AK119414 cDNA sequence. In vertebrates, there are a number of well-known alternatively splicing genes, for instance the fibronectin gene, which can generate over 20 different proteins, some of which have different patterns of localization and slightly different functions in human cells [##REF##7226224##45##]. Until relatively recently, there were very few examples of alternatively spliced genes in higher plants and still fewer examples known to generate functionally distinct proteins. Generally, alternative splicing yields two polypeptides of different size that are identical, except for the presence of a number of additional amino acids at the C-terminus of the larger isoform, eg., ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) activase [##REF##10430961##46##]; diacylglycerol kinase in tomato [##REF##11069705##47##]; and chloroplast ascorbate peroxidase in spinach [##REF##9931296##48##]. Missplicing events were observed in maize (<italic>Zea mays</italic>) <italic>Viviparous 1 </italic>(<italic>Vp1</italic>) transcripts, the majority of which are spliced incorrectly and contain insertions of intron sequences or deletions of coding region, so that they do not encode full-length proteins [##REF##12119408##49##]. Although alternative splicing is a widespread process used in higher eukaryotes to regulate gene expression and functional diversification of proteins [##REF##11173120##50##,##REF##12600935##51##], the putative control of tissue-specific expression by missplicing we observed has not been well characterized in rice.</p> <title>Preliminary Characterization of OsBgal13</title> <p>OsBgal13 is in cluster a2 of group a of the phyloenetic tree in Figure ##FIG##0##1##, along with AtBGAL9 and two <italic>P. patens </italic>Bgals. This and its general expression pattern suggest it may have a conserved maintenance function in plants. To assess what its function might be, OsBgal13 was expressed as a thioredoxin fusion protein in <italic>Escherichia coli </italic>strain OrigamiB (DE3), which we previously found to have little background β-galactosidase acitivity [##UREF##9##38##]. The protein was purified by immobilized metal affinity chromatography on nickel resin, and most of the OsBgal13 activity eluted in 5–20 mM imidazole. Further purification was achieved by gel filtration chromatography and the activity of the resulting protein was assessed. As shown in Table ##TAB##3##4##, among <italic>p</italic>-nitrophenol (<italic>p</italic>NP) glycosides, the OsBgal13 enzyme had highest activity toward <italic>p</italic>NP-β-D-galactoside (<italic>p</italic>NPGal), but also had 30% of this activity toward <italic>p</italic>NP-α-L-arabinoside (<italic>p</italic>NPAra), and low activity toward <italic>p</italic>NP-β-D-fucoside and <italic>p</italic>NP-β-D-mannoside, but no detectable activity toward <italic>p</italic>NP-β-D-glucoside and <italic>p</italic>NP-β-D-xyloside and negligible activity toward <italic>p</italic>NP-α-D-galactoside. OsBgal13 also hydrolyzed β-(1→3)-, β-(1→4)- and β-(1→6)-linked galactobiose and galactotriose, but no release of galactose or arabinose could be detected when the protein was incubated for 24 h with rice seedling alcohol insoluble residue (AIR), larchwood arabinogalactan, galactan, apple pectin, citrus pectin or oat xylan. So, despite the intriguing observation that OsBgal13 can hydrolyze both β-D-galactosyl residues and α-L-arabinosyl residues, components of type I and II arabinogalactans and xyloglucan and glucuronoarabinoxylan side chains [##REF##15012297##52##], only short galacto-oligosaccharides can be identified as possible natural substrates to date.</p> <title>Possible Functions of Other Isozymes</title> <p>The activity of mammal β-galactosidase specifically releases terminal β-galactosyl residues from glycosaminoglycans and the glycolipid GM1 ganglioside. Similarly, the digestion of the galactolipid monogalactosyldiacylglyerol by vacuolar and chloroplast β-galactosidases was reported in wheat [##REF##16663831##53##]. Likewise, β-galacosidase of tomato was reported to act in galactolipid turnover and degradation, which occurs in chloroplasts and chromoplasts during tomato fruit development [##UREF##11##54##,##UREF##12##55##]. However, relatively little is known about metabolism of other glycolipids, glycopoteins and proteoglycans involving plant β-galactosidases, though it is well characterized in mammals. The fact that OsBgal9 is similar to animal β-galactosidases and is predicted to localize to a lysosome-like vacuole, suggests it may play a similar role in rice.</p> <p>Among the plant-like OsBgals, OsBgal10 and OsBgal11 are also predicted to localize to lysosome-like vacuoles, so they could also play a role similar to animal-type β-galactosidases. However, these closely related isozymes have reproductive tissue specific expression, as are the isozymes of cluster b, OsBgal5, OsBgal12, OsBgal14 and OsBgal15. Similarly, AtBGAL11 and AtBGAL13, which clustered in group c1 with OsBgal10 and OsBgal11, and AtBGAL7 and AtBGAL15, which clustered in group b, were found to be expressed in flower [##REF##17466346##12##] and pollen [##REF##15517348##31##]. It may be that the ancestral genes for these enzymes developed reproductive-tissue specific roles before the ancestors of monocots and dicots diverged. If these functions are conserved, these isozymes may have similar roles in the two plants, though these functions remain to be determined.</p> <p>Among the other plant-like genes, it is likely many of them act in cell wall metabolism, but these roles are yet to be discerned. OsBgal1 was shown to hydrolyze β-(1→3)-, β-(1→4)- and β-(1→6)-linked oligosaccharides and, at a low level, arabinogalactan and be expressed in a range of vegetative and reproductive tissues [##UREF##9##38##], so it may have a general role on similar substrates in the cell wall. OsBgal2, OsBgal3 and OsBgal7 are closely related to it, as are 4 tomato Bgals involved in fruit ripening [##REF##10889266##13##] and an asparagus Bgal thought to act in senescence [##REF##12228457##56##], so a cell wall function seems likely for these other cluster c1 isozymes as well, though such a role must be proven. Kaneko and Kobayashi [##REF##12723614##37##] showed that OsBgal8, which is secreted from rice tissue culture cells, could release galactose from galactoxyloglucans and pectic galactans. So, a role in cell wall metabolism is quite likely for this group a isozyme.</p> <p>The role of OsBgal6 remains obscure, as no closely related enzyme has been characterized, however it seems to be specific to vegetative tissues of the young plant. The fact that the transcript found in panicle was misspliced to prevent production of active enzyme is intriguing, but the functional implications of the alternative splicing of both OsBgal6 and OsBgal11 remain to be clarified.</p>
<title>Conclusion</title> <p>Fifteen rice glycosyl hydrolase family 35 genes encoding putative β-galactosidases were identified in this study. This number is similar to the 17 seen in <italic>Arabidopsis </italic>and these plants appear to contain 9 common gene lineages present in their ancestors before they diverged [##REF##17466346##12##]. OsBgal9 was clustered into the same group as Bgals from animal species, such as <italic>H. sapiens </italic>and <italic>D. melanogaster </italic>in the phylogenetic tree, while the other rice BGals fall in a nearly plant-specific subfamily of Bgals, most of which contain a C-terminal lectin-like domain. The presence of both types of β-galactosidase in the moss <italic>P. patens</italic>, a nonvascular plant, suggests that both types of genes were present early in plant evolution. Within the plant-type Bgals, group a (Figure ##FIG##0##1##), contains 7 OsBgals and many Bgals that have been characterized from other plants (for example, 9 of <italic>Arabidopsis</italic>, 4 of <italic>C. arietinum </italic>and 3 of <italic>L. esculentum</italic>), including some with and some without C-terminal galactose-binding-lectin-like domains. Many of these proteins, including OsBgal8 [##REF##12723614##37##], appear to act on the cell-wall-derived substrates. The isozymes of phylogenetic groups b and c1, OsBgal5, OsBgal10, OsBgal11, OsBgal14, and OsBgal15, appear to have reproductive-tissue specific functions, while OsBgal6 (group c2) is predominantly in young leaves and roots. <italic>OsBgal6 </italic>appeared to undergo alternative splicing to prevent its production in panicle, and alternative splicing was also found for <italic>OsBgal11</italic>. Understanding the functions of this gene family and significance of alternative splicing in this will require further functional investigation.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Many plant β-galactosidases (Bgals) have been well characterized and their deduced biological functions mainly involve degradation of structural pectins, xyloglucans or arabinogalactoproteins in plant cell walls. However, gene multiplicity in glycosyl hydrolase family 35 (GH35), to which these proteins belong, implies diverse functions. In this study, the gene multiplicity, apparent evolutionary relationships and transcript expression of rice Bgal genes were examined, in order to predict their biological functions.</p> <title>Results</title> <p>Fifteen rice Bgal genes were identified in the plant genome, one of which encodes a protein similar to animal Bgals (<italic>OsBgal9</italic>), and the remaining 14 fall in a nearly plant-specific subfamily of Bgals. The presence of both classes of Bgals in bryophytes, as well as vascular plants, suggests both gene lineages were present early in plant evolution. All 15 proteins were predicted to contain secretory signal sequences, suggesting they have secretory pathway or external roles. RT-PCR and database analysis found two distinct lineages to be expressed nearly exclusively in reproductive tissues and to be closely related to <italic>Arabidopsis </italic>Bgals expressed most highly in flower and pollen. On the other hand, <italic>OsBgal6 </italic>is expressed primarily in young vegetative tissues, and alternative splicing in panicle prevents its production of full-length protein in this reproductive tissue. <italic>OsBgal11 </italic>also showed alternative splicing to produce different length proteins. OsBgal13 produced by recombinant expression in <italic>Escherichia coli </italic>hydrolyzed α-L-arabinoside in addition to β-D-galactoside and β-(1→3)-, β-(1→4)- and β-(1→6)- linked galacto-oligosaccharides.</p> <title>Conclusion</title> <p>Rice <italic>GH35 </italic>contains fifteen genes with a diversity of protein sequences, predicted locations and expression and splicing patterns that suggest that OsBgals enzymes may play a variety of roles in metabolism of cell wall polysaccharides, glycoproteins and glycolipids.</p>
<title>Availability &amp; requirements</title> <p>NCBI databases: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sites/entrez"/></p> <p>CAZY website: <ext-link ext-link-type="uri" xlink:href="http://www.cazy.org/"/></p> <p>SignalP program: <ext-link ext-link-type="uri" xlink:href="http://www.cbs.dtu.dk/services/SignalP/"/></p> <p>NetNGlyc: <ext-link ext-link-type="uri" xlink:href="http://www.cbs.dtu.dk/services/NetNGlyc/"/></p> <p>NCBI map viewer: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/projects/mapview"/></p> <p>PHYLIP: <ext-link ext-link-type="uri" xlink:href="http://evolution.genetics.washington.edu/phylip.html"/></p> <p>TIGR rice segmental duplication database: <ext-link ext-link-type="uri" xlink:href="http://www.tigr.org/tdb/e2k1/osa1/segmental_dup/500kb/segdup_500kb.shtml"/></p> <p> RepeatMasker: <ext-link ext-link-type="uri" xlink:href="http://www.repeatmasker.org/cgi-bin/WEBRepeatMasker"/></p> <p>Genedoc: <ext-link ext-link-type="uri" xlink:href="http://www.nrbsc.org/gfx/genedoc/index.html"/></p> <title>Authors' contributions</title> <p>WT performed rice plant cultivation, RT-PCR analysis, cDNA cloning and sequencing, and drafted the manuscript. WT, JK–C and MC carried out gene structural analysis. WT and JK–C participated in data analysis, and gene and protein analysis. JM expressed, purified and characterized the recombinant OsBgal13 protein. JK–C organized and directed the project and helped to draft the manuscript. All authors read and approved the final manuscript.</p> <title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>Tassanee Onkoksong is gratefully acknowledged for OsBgal13 expression vector construction and Rodjana Opassiri for helpful comments and discussion. Prof. Yoichi Tsumuraya and Prof. Toshihisa Kotake are thanked for providing galacto-oligosaccharides. Funding for this work was provided by the National Center for Genetic Engineering and Biotechnology of the National Science and Technology Development Agency, Thailand, grant number BT-B-06-RG-19-4608. WT was supported by postdoctoral grant of the same Agency.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Phylogenetic relationship among the rice and other glycosyl hydrolase family 35 proteins.</bold> Deduced amino acid sequences were aligned using ClustalX program [##REF##9810230##63##] and edited with Gendoc [##UREF##14##62##], then the tree were created using neighbor joining method and analyzed with 1000 bootstrap replicates, for which the percent reproducibility is shown on the branches that gave higher than 50% reproducibility. A tree was also made by maximum parsimony with the PHYLIP Protpars program, and branches supported by this analysis are shown with thicker lines. The rice sequences are named Os and a number according to the OsBgal numbers in Tables 1 and 2. <italic>Arabidopsis thaliana </italic>sequences are indicated as At, followed by the AtBGAL numbers of Ahn et al. [##REF##17466346##12##]. Sequences derived from the genome of the bryophyte <italic>Physcomitrella patens subsp. patens are given as Pp and the Genbank accession number without the initial "EDQ." Other β-galactosidase or related sequences (and their accession numbers) shown are: Ao, Asparagus officinalis </italic>(CAA54525); Aso, <italic>Aspergillus oryzae </italic>(BAE58662); Bac, <italic>Bacillus circulans </italic>(O31341); Baf, <italic>Bacteroides fragilis </italic>(CAH09349 and CAH06328); Cac, <italic>Caulobacter crescentus </italic>(Q9A4M9); Cam, <italic>Caldivirga_maquiingensis </italic>(ABW01734); Ce, <italic>Caenorhabditis elegans </italic>(O76632); Dd (<italic>Dictyostelium discoideum</italic>) (EAL64656 and XP_635837); Dm, <italic>Drosophila melanogaster </italic>(Q9VMJ5); Enf, <italic>Enterococcus faecalis </italic>(AAO81613); Hs, human (<italic>Homo sapiens) b-galactosidase </italic>(P16278); HsBglk1, human b-galactosidase-1-like protein, (Q6UWU2); Hs Bglk 3, human b-galactosidase-1-like protein 3, (Q8NCI6); Hs_4, human locus 89944 (Q8IW92); Nv, sea anemone (<italic>Nematostella vectensis</italic>, XP_001631933); Pensp, <italic>Penicillium </italic>sp. (CAF32457); Rs, radish (<italic>Raphanus sativus</italic>) (BAD20774); Xac, <italic>Xanthomonas campestris </italic>(AAM41682); Xao, <italic>Xanthomonas oryzae </italic>(BAE68425); Xl African three-toed frog (<italic>Xenopus laevis</italic>): Xl 1 (AAI24928), Xl 2 (AAH74351), Xl 3 (AAH46858).</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Distribution and direction of <italic>OsBgal </italic>genes on rice chromosomes.</bold> The chromosome numbers are indicated at the bottom of each bar. Genes lying on duplicated segments of the genome have been joined by dashed lines. Os is short for OsBgals. The gene map was generated with the NCBI map viewer.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>The exon-intron organization of <italic>OsBgal </italic>genes.</bold> Exons are shown as boxes with corresponding exons having the same pattern. Introns are represented by lines. Groupings of the genes with the same splice pattern are enclosed by red boxes.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Relative expression levels of 15 <italic>OsBgal </italic>genes in different tissues determined by semiquantitative RT-PCR.</bold> Signals were quantified and normalized to the expression of <italic>β-Actin</italic>, and the highest observed expression level of each gene was designated as 10 (black), with other tissues expression set relative to this maximal level. See Additional file ##SUPPL##0##1## for actual values and standard deviations. IP, initiating panicle; DP, differentiating panicle; EP, exerting panicle; F, flower; A, Anther; MG, milk grain; GD, grain during dry down; d, day(s); m, month(s).</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Alternative splicing of <italic>OsBgal6 </italic>(A) and <italic>OsBgal11 </italic>(B).</bold> The high expression tissue for <italic>OsBgal6 </italic>is 15-day-old rice leaf sheath, while the low expression tissue is exerting panicle. For OsBgl11, the splicing pattern of the cDNA amplified from exerting panicle of KDML105 indica rice (Genbank accession <ext-link ext-link-type="gen" xlink:href="EU603286">EU603286</ext-link>) is shown above, while the database sequence from Nipponbare japonica rice (<ext-link ext-link-type="gen" xlink:href="AK119414">AK119414</ext-link>) is shown below. Abbreviations include AAS: alternative acceptor site, E: Exon. AAS1: alternative acceptor site (CAG) located 80 nt downstream of the functional acceptor site; AAS2: alternative acceptor site (CAG) located 20 nt downstream of the functional acceptor site; AAS3: intron retention and alternative splice site within exon 11. The data for <italic>OsBgal1 </italic>and <italic>OsBgal2 </italic>were previously reported by Chantarangsee et al. [##UREF##9##38##].</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Summary of identified genes homologous to glycosyl hydrolase family 35 galactosidase</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Gene name</td><td align="center">BGI entries (AAAA...)<sup>a</sup></td><td align="center">Genomic Contig.</td><td align="center">Chr.<sup>b</sup></td><td align="center">UniGene</td><td align="center">Gene Locus (NCBI, RGP)</td><td align="center">Gene Locus (TIGR)</td><td align="center">Corresponding cDNA<sup>c</sup></td><td align="center" colspan="2">Number<sup>d</sup></td><td align="center">Tissue librarys<sup>e</sup></td><td align="center" colspan="8">Expression profile (TPM) Transcript per millon<sup>f</sup></td><td align="center">GEO profile<sup>g</sup></td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td colspan="2"><hr/></td><td/><td colspan="8"><hr/></td><td/></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">mRNA</td><td align="center">ESTs</td><td/><td align="center">Cl.</td><td align="center">Fw.</td><td align="center">Lf.</td><td align="center">P</td><td align="center">Rt.</td><td align="center">Sd.</td><td align="center">St.</td><td align="center">Vm</td><td/></tr></thead><tbody><tr><td align="center"><italic>OsBgal1</italic></td><td align="center">1002951</td><td align="center">AP008209.1</td><td align="center">3</td><td align="center">Os.13559</td><td align="center">Os03g0165400</td><td align="center">LOC_Os03g06940</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AF499737">AF499737</ext-link></td><td align="center">5</td><td align="center">99</td><td align="center">Fw, Cl, P, Mx, Wp, UnT, Vm, Sd, St, Rt</td><td align="center">121</td><td align="center">241</td><td align="center">46</td><td align="center">120</td><td align="center">14</td><td align="center">92</td><td align="center">23</td><td align="center">217</td><td align="center">Up: CKl, AB&amp;GBd, AnX, Dw: AB&amp;GB, NaCl, Na<sub>3</sub>AsO<sub>4</sub>, CKrl</td></tr><tr><td align="center"><italic>OsBgal2</italic></td><td align="center">1003471</td><td align="center">AP004729.3, AP008212.1</td><td align="center">6</td><td align="center">Os.87715</td><td align="center">Os06g0573600</td><td align="center">LOC_Os06g37560</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK102756">AK102756</ext-link></td><td align="center">4</td><td align="center">201</td><td align="center">P, Lf, Mx, St, Rt, Cl, Fw, Wp, UnT</td><td align="center">66</td><td align="center">80</td><td align="center">291</td><td align="center">391</td><td align="center">205</td><td align="center">0</td><td align="center">118</td><td align="center">0</td><td align="center">Up: CKl, AB&amp;GBd, NaCl Dw: AnX, AB&amp;GB, Na<sub>3</sub>AsO<sub>4</sub>, CKrl</td></tr><tr><td align="center"><italic>OsBgal3</italic></td><td align="center">01012445, 01022298</td><td align="center">AP008207.1, AP003546.4</td><td align="center">1</td><td align="center">Os.18562</td><td align="center">Os01g0580200</td><td align="center">LOC_Os01g39830</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK103045">AK103045</ext-link></td><td align="center">2</td><td align="center">44</td><td align="center">Cl, Wp, P, St, Fw, Mx, UnT, Vm</td><td align="center">115</td><td align="center">43</td><td align="center">0</td><td align="center">37</td><td align="center">0</td><td align="center">30</td><td align="center">47</td><td align="center">217</td><td align="center">Up: CKl, AB&amp;GBd, AB&amp;GB, Na<sub>3</sub>AsO<sub>4</sub>, NaCl Dw: AnX, CKrl</td></tr><tr><td align="center"><italic>OsBgal4</italic></td><td align="center">ND</td><td align="center">AP003297.3, AP008207.1</td><td align="center">1</td><td align="center">Os.9892</td><td align="center">Os01g0875500</td><td align="center">LOC_Os01g65460</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK101399">AK101399</ext-link></td><td align="center">3</td><td align="center">107</td><td align="center">Cl, Fw, Mx, P, Rt, UnT, St, Wp</td><td align="center">279</td><td align="center">263</td><td align="center">0</td><td align="center">45</td><td align="center">29</td><td align="center">30</td><td align="center">0</td><td align="center">0</td><td align="center">Up: CKl, AnX, AB&amp;GBd, Na3AsO4, NaCl Dw: AB&amp;GB</td></tr><tr><td align="center"><italic>OsBgal5</italic></td><td align="center">1000473.1</td><td align="center">AP003447.2, AP008207.1, AP003445.2</td><td align="center">1</td><td align="center">Os.41043</td><td align="center">Os01g0533400</td><td align="center">LOC_Os01g34920</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK119447">AK119447</ext-link></td><td align="center">2</td><td align="center">3</td><td align="center">Mx, St</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">Up: Na<sub>3</sub>AsO<sub>4</sub>, CKrl Dw: CKl, AnX, NaCl</td></tr><tr><td align="center"><italic>OsBgal6</italic></td><td align="center">101013714.1</td><td align="center">AC135419.2, AC135429.2</td><td align="center">5</td><td align="center">Os.82841</td><td align="center">P0636F09.15</td><td align="center">LOC_Os05g35360</td><td align="center"><ext-link ext-link-type="gen" xlink:href="EU602310">EU602310</ext-link></td><td align="center">ND</td><td align="center">1</td><td align="center">UnT</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">-</td></tr><tr><td align="center"><italic>OsBgal7</italic></td><td align="center">01011567, 01017443</td><td align="center">AP008208.1</td><td align="center">2</td><td align="center">Os.14358</td><td align="center">Os02g0219200</td><td align="center">LOC_Os02g12730</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK059059">AK059059</ext-link></td><td align="center">7</td><td align="center">57</td><td align="center">P, Mx, St, Cl, UnT, Fw, Wp, Rt</td><td align="center">24</td><td align="center">21</td><td align="center">0</td><td align="center">135</td><td align="center">14</td><td align="center">92</td><td align="center">78</td><td align="center">0</td><td align="center">Up: AnX, CKl, AB&amp;GB, AB&amp;GBd, NaCl, Na3AsO4 Dw: CKrl</td></tr><tr><td align="center"><italic>OsBgal8</italic></td><td align="center">1001024</td><td align="center">AP008209.1, DP000009.2</td><td align="center">3</td><td align="center">Os.22360</td><td align="center">Os03g0255100</td><td align="center">LOC_Os03g15020</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK067479">AK067479</ext-link></td><td align="center">3</td><td align="center">253</td><td align="center">Fw, Mx, St, Lf, P, Rt, Cl, UnT, Sd, Wp</td><td align="center">42</td><td align="center">761</td><td align="center">133</td><td align="center">120</td><td align="center">190</td><td align="center">216</td><td align="center">228</td><td align="center">0</td><td align="center">Up: AnX, CKl, AB&amp;GBd, NaCl Dw: Na<sub>3</sub>AsO<sub>4</sub></td></tr><tr><td align="center"><italic>OsBgal9</italic></td><td align="center">ND</td><td align="center">AP008211.1</td><td align="center">5</td><td align="center">Os.14570</td><td align="center">Os05g0539400</td><td align="center">LOC_Os05g46200</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK068572">AK068572</ext-link></td><td align="center">2</td><td align="center">66</td><td align="center">Lf, St, P, Mx, Cl, UnT, Wp, Sd, Rt</td><td align="center">12</td><td align="center">0</td><td align="center">122</td><td align="center">52</td><td align="center">14</td><td align="center">247</td><td align="center">110</td><td align="center">0</td><td align="center">Up: CKl, AB&amp;GBd, AB&amp;GB, NaCl Dw: AnX, Na<sub>3</sub>AsO<sub>4</sub></td></tr><tr><td align="center"><italic>OsBgal10</italic></td><td align="center">1005991</td><td align="center">AP008214.1, AP003912.3</td><td align="center">8</td><td align="center">Os.18310</td><td align="center">Os08g0549200</td><td align="center">LOC_Os08g43570</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK069066">AK069066</ext-link></td><td align="center">3</td><td align="center">137</td><td align="center">Mx, Fw, P, St, UnT</td><td align="center">0</td><td align="center">249</td><td align="center">0</td><td align="center">158</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">Up: AB&amp;GBd, AB&amp;GB, AnX Dw: Na<sub>3</sub>AsO<sub>4</sub>, CKrl</td></tr><tr><td align="center"><italic>OsBgal11</italic></td><td align="center">1013958</td><td align="center">AP008215.1</td><td align="center">9</td><td align="center">Os.49945</td><td align="center">Os09g0539200</td><td align="center">LOC_Os09g36810</td><td align="center"><ext-link ext-link-type="gen" xlink:href="EU603286">EU603286</ext-link></td><td align="center">2</td><td align="center">7</td><td align="center">Mx, Sd</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">61</td><td align="center">0</td><td align="center">0</td><td align="center">Up: AnX, Na<sub>3</sub>AsO<sub>4</sub>, NaCl Dw: CKl</td></tr><tr><td align="center"><italic>OsBgal12</italic></td><td align="center">10005757.1</td><td align="center">AP008216.1, DP000086.1</td><td align="center">10</td><td align="center">Os.46702</td><td align="center">Os10g0330600</td><td align="center">LOC_Os10g18400</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK119350">AK119350</ext-link></td><td align="center">3</td><td align="center">3</td><td align="center">Mx, Fw</td><td align="center">0</td><td align="center">7</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">Up: AnX</td></tr><tr><td align="center"><italic>OsBgal13</italic></td><td align="center">ND</td><td align="center">AP008218.1</td><td align="center">12</td><td align="center">Os.52193</td><td align="center">Os12g0429200</td><td align="center">LOC_Os12g24170</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK065546">AK065546</ext-link></td><td align="center">2</td><td align="center">36</td><td align="center">Mx, St, Lf, Fw, UnT, P, Wp, Cl</td><td align="center">6</td><td align="center">21</td><td align="center">34</td><td align="center">15</td><td align="center">0</td><td align="center">0</td><td align="center">55</td><td align="center">0</td><td align="center">Up: CKl, AB&amp;GBd, AB&amp;GB Dw: AnX, Na<sub>3</sub>AsO<sub>4</sub>, NaCl</td></tr><tr><td align="center"><italic>OsBgal14</italic></td><td align="center">10006090</td><td align="center">ND</td><td align="center">10</td><td align="center">Os.22528</td><td align="center">J090043H02</td><td align="center">LOC_Os10g19960</td><td align="center"><ext-link ext-link-type="gen" xlink:href="AK242960">AK242960</ext-link></td><td align="center">2</td><td align="center">58</td><td align="center">Fw, P, UnT, Sd, Wp</td><td align="center">0</td><td align="center">358</td><td align="center">0</td><td align="center">7</td><td align="center">0</td><td align="center">61</td><td align="center">0</td><td align="center">0</td><td align="center">-</td></tr><tr><td align="center"><italic>OsBgal15</italic></td><td align="center">1003798.1</td><td align="center">AP004733</td><td align="center">6</td><td align="center">NF</td><td align="center">AP004733</td><td align="center">LOC_Os06g42310</td><td align="center"><ext-link ext-link-type="gen" xlink:href="EU051629">EU051629</ext-link></td><td align="center">ND</td><td/><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">-</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Predicted rice GH family 35 β-galactosidase protein properties and locations.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene Name</td><td align="left">Gene Locus <break/>(NCBI, RGP)</td><td align="center" colspan="3">Pre-protein</td><td align="center" colspan="5">Mature Protein</td></tr><tr><td/><td/><td colspan="3"><hr/></td><td colspan="5"><hr/></td></tr><tr><td/><td/><td align="left">MW<sup>a </sup>(kDa)</td><td align="left">AA<sup>b</sup></td><td align="left">Cleavage site<sup>c</sup></td><td align="left">MW<sup>a</sup></td><td align="left">AA<sup>b</sup></td><td align="left">pI<sup>a</sup></td><td align="center">N-gly site<sup>d</sup></td><td align="left">Possible destination<sup>e</sup></td></tr></thead><tbody><tr><td align="left"><italic>OsBgal1</italic></td><td align="left">Os03g0165400</td><td align="left">92.2</td><td align="left">841</td><td align="left">25–26</td><td align="left">89.8</td><td align="left">816</td><td align="left">5.78</td><td align="center">2</td><td align="left">Gol, Vac, ER ext</td></tr><tr><td align="left"><italic>OsBgal2</italic></td><td align="left">Os06g0573600</td><td align="left">78.3</td><td align="left">715</td><td align="left">20–21</td><td align="left">76.5</td><td align="left">695</td><td align="left">7.65</td><td align="center">0</td><td align="left">Gol, PM, IMP, MII, ER ext</td></tr><tr><td align="left"><italic>OsBgal3</italic></td><td align="left">Os01g0580200</td><td align="left">91.5</td><td align="left">827</td><td align="left">24–25</td><td align="left">89.5</td><td align="left">803</td><td align="left">5.59</td><td align="center">2</td><td align="left">Gol, PM, IMP, MII, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal4</italic></td><td align="left">Os01g0875500</td><td align="left">94.7</td><td align="left">851</td><td align="left">29–30</td><td align="left">91.7</td><td align="left">822</td><td align="left">6.76</td><td align="center">6</td><td align="left">Gol, Cis-Gol, MII, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal5</italic></td><td align="left">Os01g0533400</td><td align="left">91.7</td><td align="left">827</td><td align="left">25–26</td><td align="left">88.9</td><td align="left">802</td><td align="left">5.66</td><td align="center">8</td><td align="left">Gol, MII, ER</td></tr><tr><td align="left"><italic>OsBgal6</italic></td><td align="left">P0636F09.15</td><td align="left">90.6</td><td align="left">811</td><td align="left">17–18</td><td align="left">88.9</td><td align="left">794</td><td align="left">6.61</td><td align="center">5</td><td align="left">IMP, MC in, Sec, Cy, N</td></tr><tr><td align="left"><italic>OsBgal7</italic></td><td align="left">Os02g0219200</td><td align="left">78</td><td align="left">729</td><td align="left">36–37</td><td align="left">76.8</td><td align="left">693</td><td align="left">8.71</td><td align="center">0</td><td align="left">Gol, MII, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal8</italic></td><td align="left">Os03g0255100</td><td align="left">103.7</td><td align="left">956</td><td align="left">62–63</td><td align="left">96.5</td><td align="left">894</td><td align="left">5.62</td><td align="center">7</td><td align="left">PM, IMP, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal9</italic></td><td align="left">Os05g0539400</td><td align="left">75.7</td><td align="left">673</td><td align="left">20–21</td><td align="left">73.5</td><td align="left">653</td><td align="left">5.57</td><td align="center">6</td><td align="left">Vac-Lys, IMP, MC in, Per, ERl, N, Sec, Cy</td></tr><tr><td align="left"><italic>OsBgal10</italic></td><td align="left">Os08g0549200</td><td align="left">94.7</td><td align="left">848</td><td align="left">23–24</td><td align="left">92.8</td><td align="left">825</td><td align="left">9.1</td><td align="center">5</td><td align="left">Vac-Lys, PM, IMP, MC in, ER, Cy, Sec</td></tr><tr><td align="left"><italic>OsBgal11</italic></td><td align="left">Os09g0539200</td><td align="left">94.11</td><td align="left">838</td><td align="left">25–26</td><td align="left">91.9</td><td align="left">813</td><td align="left">6.52</td><td align="center">3</td><td align="left">Vac-Lys, Gol, PM, IMP, MII, Mic, ER, MC in, Sec</td></tr><tr><td align="left"><italic>OsBgal12</italic></td><td align="left">Os10g0330600</td><td align="left">92</td><td align="left">828</td><td align="left">23–24</td><td align="left">89.7</td><td align="left">805</td><td align="left">6.07</td><td align="center">10</td><td align="left">Gol, MII, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal13</italic></td><td align="left">Os12g0429200</td><td align="left">101.0</td><td align="left">919</td><td align="left">31–32</td><td align="left">97.9</td><td align="left">888</td><td align="left">5.56</td><td align="center">5</td><td align="left">IMP, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal14</italic></td><td align="left">J090043H02</td><td align="left">92.2</td><td align="left">828</td><td align="left">28–29</td><td align="left">89.6</td><td align="left">800</td><td align="left">5.99</td><td align="center">9</td><td align="left">PM, IMP, MII, ER, Sec</td></tr><tr><td align="left"><italic>OsBgal15</italic></td><td align="left">AP004733</td><td align="left">89.6</td><td align="left">809</td><td align="left">23–24</td><td align="left">87.3</td><td align="left">789</td><td align="left">6.17</td><td align="center">8</td><td align="left">Gol, MII, ER, Sec</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Identification of LTR-retrotransposon elements and DNA sequences for interspersed repeats and low complexity DNA sequences within <italic>OsBgal6 </italic>and <italic>OsBgal13</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Genes</td><td align="center" colspan="7">LTR-retrotransposon elements and DNA sequences for interspersed repeats</td></tr><tr><td/><td colspan="7"><hr/></td></tr><tr><td/><td align="left">intron 1</td><td align="left">intron 3</td><td align="left">intron 8</td><td align="left">intron 10</td><td align="left">intron 13</td><td align="left">intron 14</td><td align="left">intron 16</td></tr></thead><tbody><tr><td align="left"><italic>OsBgal6</italic></td><td align="left">Stowaway41_OS, OLO24C, OLO24B, RPO_OS, AT rich Low Complexity, AMYLTP, TA rich Low Complexity, OSTE15</td><td/><td/><td/><td/><td/><td/></tr><tr><td align="left"><italic>OsBgal13</italic></td><td align="left">GC rich Low complexity</td><td align="left">AT rich Low complexity</td><td align="left">MUDRN2</td><td align="left">OSTE28, OSTE33, CRM-I_OS LTR/Gypsy,</td><td align="left">OSLINE1-5, AT rich Low complexity</td><td align="left">AT rich Low complexity</td><td align="left">RdSpm875A_3.1, ID-2 DNA/Tourist</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Hydrolysis of <italic>p</italic>-nitrophenyl glycosides, oligosaccharides and polysaccharides by OsBgal13.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Substrate</td><td align="center">Activity</td></tr></thead><tbody><tr><td align="left"><italic>p</italic>NP-glycosides</td><td align="center">(Percent Relative Activity)</td></tr><tr><td align="left"><italic>p</italic>NP-β-D-galactopyranoside</td><td align="center">100</td></tr><tr><td align="left"><italic>p</italic>NP-β-D-glucopyranoside</td><td align="center">N.D.</td></tr><tr><td align="left"><italic>p</italic>NP-β-D-fucopyranoside</td><td align="center">7.4</td></tr><tr><td align="left"><italic>p</italic>NP-β-D-mannopyranoside</td><td align="center">4.5</td></tr><tr><td align="left"><italic>p</italic>NP-β-D-xylopyranoside</td><td align="center">N.D.</td></tr><tr><td align="left"><italic>p</italic>NP-β-L-arabinopyranoside</td><td align="center">30</td></tr><tr><td align="left"><italic>p</italic>NP-β-D-galactopyranoside</td><td align="center">1.5</td></tr><tr><td colspan="2"><hr/></td></tr><tr><td align="left">Oligosaccharides</td><td align="center">(estimated relative activity)</td></tr><tr><td align="left">β-1,3-galactobiose</td><td align="center">+++</td></tr><tr><td align="left">β-1,3-galactotriose</td><td align="center">++</td></tr><tr><td align="left">β-1,4-galactobiose</td><td align="center">++ *</td></tr><tr><td align="left">β-1,4-galactotriose</td><td align="center">++ *</td></tr><tr><td align="left">β-1,6-galactobiose</td><td align="center">+++</td></tr><tr><td align="left">β-1,6-galactotriose</td><td align="center">++</td></tr><tr><td align="left">Arabinogalactan, xylan from oat spelts, birchwood xylan, polygalacturonic acid, apple pectin, citric pectin, galactan</td><td align="center">N.D.</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>Sequences of primers used for semi-quantitative RT-PCR analysis of β-galactosidase gene expression.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene name</td><td align="left">GenBank ID</td><td align="left">Sense primer</td><td align="left">Antisense primer</td><td align="left">Annealing Tm (°C)</td></tr></thead><tbody><tr><td align="left">Actin</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK100267">AK100267</ext-link></td><td align="left">ACTCTGGTGATGGTGTCAGCC</td><td align="left">GTCAGCAATGCCAGGGAACATA</td><td align="left">54</td></tr><tr><td align="left"><italic>UBQ6</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM464194">XM464194</ext-link></td><td align="left">TCCTCCGTCTCAGGGGAG</td><td align="left">CTTGCCAGCGAAGATCAGAC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal1</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK102192">AK102192</ext-link></td><td align="left">CAAAGCACACAGAAAGCGAT</td><td align="left">TGCTCACCGCACAATCAACGA</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal2</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK102756">AK102756</ext-link></td><td align="left">AGGAAAGTGGGGGCGTATAG</td><td align="left">TCCCATTTACAACTCAACGT</td><td align="left">56</td></tr><tr><td align="left"><italic>OsBgal3</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK103045">AK103045</ext-link></td><td align="left">CCAAGGGGCTGTATGTATGGTC</td><td align="left">CCTGAGAGAATTCATTCACATACGG</td><td align="left">56</td></tr><tr><td align="left"><italic>OsBgal4</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK102715">AK102715</ext-link></td><td align="left">TCCATCGCTACAGATTCGCTC</td><td align="left">TCCAGAAATATCATGACGCGAC</td><td align="left">56</td></tr><tr><td align="left"><italic>OsBgal5</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK119447">AK119447</ext-link></td><td align="left">CTCATCTGCTTGCTTCATC</td><td align="left">CTAAAGTTGCCCTTCTCATC</td><td align="left">50</td></tr><tr><td align="left"><italic>OsBgal6</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="CR292731">CR292731</ext-link></td><td align="left">AGTCTTGCATAGGCAGGAG</td><td align="left">TCTGAACGAAGGTATCGCAC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal7</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK059059">AK059059</ext-link></td><td align="left">CCACCATTTGATACAGTCGTCG</td><td align="left">TTCCCGAGCAACGCAAAC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal8</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK067479">AK067479</ext-link></td><td align="left">CCTGACAGGTTTGATAGTGCTCG</td><td align="left">TGCTTTTCTTCACACAAGTGCATC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal9</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK068572">AK068572</ext-link></td><td align="left">GAAGGATCCAGATTTCACATGC</td><td align="left">TGCTGTTCATGTCATCATGTGC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal10</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK069066">AK069066</ext-link></td><td align="left">TCCAAGAGGCCTCCTCGTTC</td><td align="left">CTACATATAAAACCATGGACGATGGTG</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal11</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK119414">AK119414</ext-link></td><td align="left">GTGCAGGTGAGATGCAAGGTATC</td><td align="left">CAAACTGTCTGTCAACCTGTGATGG</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal12</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK119350">AK119350</ext-link></td><td align="left">CAACACAGCAAAACCATCT</td><td align="left">TTCCACGAAACAAAGTAAAGACA</td><td align="left">56</td></tr><tr><td align="left"><italic>OsBgal13</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK065546">AK065546</ext-link></td><td align="left">CCGAGGAGTCCTCAAAGATTTAGC</td><td align="left">GGAAATCTCCTTTGCATTTTTATTCAC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal14</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AK242960">AK242960</ext-link></td><td align="left">GTCGAAGGTGGCTTATGACG</td><td align="left">AATCGACAGTGCGGTATCTC</td><td align="left">54</td></tr><tr><td align="left"><italic>OsBgal15</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AP004733">AP004733</ext-link></td><td align="left">CGTACAAGGCTTTCACAGAAG</td><td align="left">CTACGATTACTTGATCACACTC</td><td align="left">54</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Relative expression levels of 15 <italic>OsBgal </italic>genes in different tissues determined by semiquantitative RT-PCR</bold>. Signals were quantified and normalized to the expression of β-actin. The values given are means with standard deviations from triplicate experiments.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>a</sup>Contig number in Beijing Genome Institute (the numbers given are the Genbank Accession without the initial 'AAAA').</p><p><sup>b</sup>Chromosome locations from mapping of the corresponding genes on the 12 rice chromosomes in GenBank.</p><p><sup>c</sup>Accessions of full-length cDNA clones of japonica rice databases [##UREF##15##64##] or indica clones from this study (<italic>OsBgal6</italic>, <italic>OsBgal11 </italic>and <italic>OsBgal15</italic>).</p><p><sup>d</sup>The number of mRNA and EST identified by UniGene from the NCBI Genbank database.</p><p><sup>e</sup>Tissue libraries means the source tissue from which the corresponding ESTs were isolated. Tissue abbreviations are: Cl: callus, Fw: Flower, Lf: leaf, P: Panicle, Rt: root, Sd: Seed, St: Stem, Mx: Mixed Tissues, Vm: Vegetative meristem, Wp: whole plant, UnT: unspecified tissue.</p><p><sup>f</sup>Expression profile means transcript frequency per million calculated from the number of ESTs in a certain library. ND: not determined.</p><p><sup>g</sup>GEO (Gene Expression Omnibus public repository of microarray and other forms of high-throughput data submitted by the scientific community). Abbreviations are: Up: up-regulated, Dw:down-regulate, CKrl: Cytokinin effect on rice roots and leaves: time course, CKl: Cytokinin-inducible type-A response regulator OsRR6 overexpression effect on rice leaves, AB&amp;GB: Abscisic acid and gibberellin effect on calluses, gene expression array-based, log ratio, AB&amp;GBd: Abscisic acid and gibberellin effect on calluses (dye-swap), AnX: Anoxia effect on rice coleoptiles, NaCl: Salinity stress response, Na<sub>3</sub>AsO<sub>4</sub>: Sodium arsenate effect on Azucena and Bala varieties.</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>determined by ProtParam. <sup>b </sup>AA means number of amino acids. <sup>c</sup>predicted by SignalP [##UREF##7##35##]. <sup>d</sup>predicted by NetNGlyc at the Expasy proteomics server <ext-link ext-link-type="uri" xlink:href="http://www.expasy.org"/> and manually curated to remove NPS/T sites, <sup>e</sup>cellular locations predicted by PSORT [##REF##2440339##65##], abbreviated as PM: plasma membrane, MII: Type II membrane protein, Gol: Golgi, IMP: Integral Membrane Protein, Sec: Secreted, ER: Endoplasmic reticulum; ext: extracellular; l: lumen, MC: Mitochondrial; in: inner membrane, Cy: Cytoplasmic, N: nuclear, Vac: vacuoles, Vac-Lys: Lysosome-like vacuoles, Per: Peroxisomal, and Mic: microsomal.</p></table-wrap-foot> <table-wrap-foot><p>RepeatMasker <ext-link ext-link-type="uri" xlink:href="http://www.repeatmasker.org/cgi-bin/WEBRepeatMasker"/> was used to identify these sequences.</p></table-wrap-foot> <table-wrap-foot><p>N.D. means not detectable. * For β-1,4-galactotriose and β-1,4-galactobiose transglycosylation was observed as much or more than hydrolysis.</p></table-wrap-foot>
<graphic xlink:href="1471-2229-8-84-1"/> <graphic xlink:href="1471-2229-8-84-2"/> <graphic xlink:href="1471-2229-8-84-3"/> <graphic xlink:href="1471-2229-8-84-4"/> <graphic xlink:href="1471-2229-8-84-5"/>
<media xlink:href="1471-2229-8-84-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["O'Brien", "Scriver CR, Beaudet AL, Sly WS, Valle D"], "given-names": ["JS"], "article-title": ["\u03b2-Galactosidase deficiency (GM1 gangliosidosis, galactosialidosis, and Morquio syndrome type B); ganglioside sialidase deficiency (mucolipidosis IV)"], "source": ["The Metabolic Basis of Inherited Disease"], "year": ["1989"], "publisher-name": ["New York: McGraw-Hill Book Co"], "fpage": ["1797"], "lpage": ["1806"]}, {"surname": ["En\u00e9as-Filho", "Sud\u00e9rio", "Gomes-Filho", "Prisco"], "given-names": ["J", "FB", "E", "JT"], "article-title": ["Multiple forms of cotyledonary \u03b2-galactosidases from "], "italic": ["Vigna unguiculata "], "source": ["Rev Bras Bot"], "year": ["2000"], "volume": ["23"], "fpage": ["69"], "lpage": ["76"]}, {"surname": ["Giannakouros", "Karagiorgos", "Simos"], "given-names": ["T", "A", "G"], "article-title": ["Expression of \u03b2-galactosidase multiple forms during barley ("], "italic": ["Hordeum vulgare"], "source": ["Physiol Plant"], "year": ["1991"], "volume": ["82"], "fpage": ["413"], "lpage": ["418"]}, {"surname": ["Hirano", "Tsumuraya", "Hashimoto"], "given-names": ["Y", "Y", "Y"], "article-title": ["Characterization of spinach leaf \u03b1-L-arabinofuranosidases and \u03b2-galactosidases and their synergistic action on an endogenous arabinogalactan-protein"], "source": ["Physiol Plant"], "year": ["1994"], "volume": ["92"], "fpage": ["286"], "lpage": ["296"]}, {"surname": ["Ross", "Redgwell", "Macrae"], "given-names": ["GS", "RJ", "EA"], "article-title": ["beta-galactosidase: Isolation and activity against specific fruit cell-wall polysaccharides"], "source": ["Planta"], "year": ["1993"], "volume": ["189"], "fpage": ["499"], "lpage": ["506"]}, {"surname": ["Andrews", "Li"], "given-names": ["PK", "S"], "article-title": ["Partial purification and characterization of beta-D-galactosidase from sweet cherry, a nonclimacteric fruit"], "source": ["J Agric Food Chem"], "year": ["1994"], "volume": ["2"], "fpage": ["2177"], "lpage": ["2182"]}, {"surname": ["Alcantara", "Dietrich", "Buckeridge"], "given-names": ["PHN", "SMC", "MS"], "article-title": ["Xyloglucan mobilisation and purification of a (XLLG/XLXG) specific-galactosidase from cotyledons of "], "italic": ["Copaifera langsdorffii"], "source": ["Plant Physiol Biochem"], "year": ["1999"], "volume": ["37"], "fpage": ["653"], "lpage": ["663"]}, {"surname": ["Tin\u00e9", "Silva", "Lima", "Carpita", "Buckeridge"], "given-names": ["MA", "CO", "DU", "NC", "MS"], "article-title": ["Fine structure of a mixed-oligomer storage xyloglucan from seeds of "], "italic": ["Hymenaea courbaril"], "source": ["Carbohydr Polym"], "year": ["2006"], "volume": ["66"], "fpage": ["444"], "lpage": ["454"]}, {"surname": ["Konno", "Tsumuki"], "given-names": ["H", "H"], "article-title": ["Purification of a \u03b2-galactosidase from rice shoots and its involvement in hydrolysis of the natural substrate in cell walls"], "source": ["Physiol Plant"], "year": ["1993"], "volume": ["89"], "fpage": ["40"], "lpage": ["47"]}, {"surname": ["Chantarangsee", "Tanthanuch", "Fujimura", "Fry", "Cairns"], "given-names": ["M", "W", "T", "SC", "JK"], "article-title": ["Molecular characterization of \u03b2-galactosidases from germinating rice ("], "italic": ["Oryza sativa"], "source": ["Plant Sci"], "year": ["2007"], "volume": ["173"], "fpage": ["118"], "lpage": ["134"]}, {"surname": ["Yoshida", "Aha"], "given-names": ["S", "SB"], "article-title": ["The accumulation process of carbohydrate in rice varieties in relation to their response to nitrogen in the tropics"], "source": ["Soil Sci Plant Nutr"], "year": ["1968"], "volume": ["14"], "fpage": ["153"], "lpage": ["162"]}, {"surname": ["G\u00fc\u00e7l\u00fc", "Paulin", "Soudain"], "given-names": ["J", "A", "P"], "article-title": ["Changes in polar lipids during ripening and senescence of cherry tomato ("], "italic": ["Lycopersicon esculentum"], "source": ["Physiol Plant"], "year": ["1989"], "volume": ["77"], "fpage": ["413"], "lpage": ["419"]}, {"surname": ["Whitaker"], "given-names": ["BD"], "article-title": ["Changes in galactolipid and phospholipid levels of tomato fruits stored at chilling and nonchilling temperatures"], "source": ["Phytochem"], "year": ["1992"], "volume": ["31"], "fpage": ["2627"], "lpage": ["2630"]}, {"surname": ["Gupta", "Jung", "Brunak"], "given-names": ["R", "E", "S"], "article-title": ["Prediction of N-glycosylation sites in human proteins"], "source": ["preparation"], "year": ["2004"]}, {"surname": ["Nicholas", "Nicholas", "Deerfield"], "given-names": ["KB", "HB", "DW"], "suffix": ["Jr", "II"], "article-title": ["GeneDoc: Analysis and visualization of genetic variation"], "source": ["EMBNEW News"], "year": ["1997"], "volume": ["4"], "fpage": ["14"]}, {"surname": ["Felsenstein"], "given-names": ["J"], "source": ["PHYLIP (Phylogeny Inference Package) version 36"], "year": ["2005"], "publisher-name": ["Distributed by the author, Department of Genome Sciences, University of Washington, Seattle"]}, {"surname": ["Sambrook", "Fitsch", "Maniatis"], "given-names": ["J", "EF", "T"], "source": ["Molecular Cloning: A Laboratory Manual"], "year": ["1989"], "publisher-name": ["Cold Spring Harbor, Cold Spring Harbor Press"]}]
BMC Plant Biol. 2008 Jul 30; 8:84