Genomic and expression analysis of glycosyl hydrolase family 35 genes from rice (Oryza sativa L.)
© Tanthanuch et al; licensee BioMed Central Ltd. 2008
Received: 26 April 2008
Accepted: 30 July 2008
Published: 30 July 2008
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.
Fifteen rice Bgal genes were identified in the plant genome, one of which encodes a protein similar to animal Bgals (OsBgal9), 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 Arabidopsis Bgals expressed most highly in flower and pollen. On the other hand, OsBgal6 is expressed primarily in young vegetative tissues, and alternative splicing in panicle prevents its production of full-length protein in this reproductive tissue. OsBgal11 also showed alternative splicing to produce different length proteins. OsBgal13 produced by recombinant expression in Escherichia coli hydrolyzed α-L-arabinoside in addition to β-D-galactoside and β-(1→3)-, β-(1→4)- and β-(1→6)- linked galacto-oligosaccharides.
Rice GH35 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.
β-Galactosidases (Bgals) (EC 188.8.131.52) 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 [1, 2], thereby acting to control fruit softening during ripening , hydrolysis of dietary lactose [4, 5] and degradation of glycolipids and proteoglycans in mammals [6, 7], and metabolism of lactose and other galactosides in microorganisms [8, 9]. The extensive diversity of Bgals raises the question of what is the basis of their different substrate specificities, which help determine their functions.
Glycosyl hydrolases (GH) have been grouped and classified as families on the basis of structural similarity  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 , 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 Arabidopsis , tomato , papaya , apple , Vigna unguiculata  and barley , suggesting that GH35 gene multiplicity is ubiquitous in plants.
Analysis of plant Bgal activities led to the proposal of two classes, I and II . 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 [19, 20], as reported in spinach leaves  and hypocotyls and young leaves of radish . 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 [1, 2]. The apparent involvement of Bgals in pectin disassembly during fruit ripening has been reported in various plant species, including tomato [3, 21, 22], muskmelon , kiwifruit , persimmon , sweet cherry , mango , and peach . Bgals expression also accompanies many stages of plant development in other tissues, for example, spinach leaves , mungbean seedlings , radish hypocotyls and young leaves , and meristem zones of roots, cotyledons, vascular tissues, trichomes, and pollen of tobacco [30, 31]. Bgals that appear to catalyze galactose removal from xyloglucans during their disassembly were observed in cotyledons of nasturtium (Tropueolum mujus L.) seed , Copaifera langsdorffii seed , and Hymenaea courbaril seed, where they were reported to act in cooperation with α-xylosidase, β-glucosidase, and other enzymes to achieve xyloglucan degradation [34, 35]. Therefore, Bgals appear to play a role in cell wall remodeling in many plant processes.
In rice, GH35 has not yet been well characterized, and the biological function of rice Bgals remain mysterious. Konno and Tsumuki  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  isolated a Bgal with 40 and 47 kDa subunits from the medium of rice suspension cells. Recently, Chantarangsee et al.  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 Bgals have been determined and OsBgal13 was expressed in Escherichia coli to gain a clue to OsBgal physiological functions.
Results and Discussion
Identification of rice GH35 genes and their protein products
Summary of identified genes homologous to glycosyl hydrolase family 35 galactosidase
BGI entries (AAAA...)a
Gene Locus (NCBI, RGP)
Gene Locus (TIGR)
Expression profile (TPM) Transcript per millonf
Fw, Cl, P, Mx, Wp, UnT, Vm, Sd, St, Rt
Up: CKl, AB&GBd, AnX, Dw: AB&GB, NaCl, Na3AsO4, CKrl
P, Lf, Mx, St, Rt, Cl, Fw, Wp, UnT
Up: CKl, AB&GBd, NaCl Dw: AnX, AB&GB, Na3AsO4, CKrl
Cl, Wp, P, St, Fw, Mx, UnT, Vm
Up: CKl, AB&GBd, AB&GB, Na3AsO4, NaCl Dw: AnX, CKrl
Cl, Fw, Mx, P, Rt, UnT, St, Wp
Up: CKl, AnX, AB&GBd, Na3AsO4, NaCl Dw: AB&GB
AP003447.2, AP008207.1, AP003445.2
Up: Na3AsO4, CKrl Dw: CKl, AnX, NaCl
P, Mx, St, Cl, UnT, Fw, Wp, Rt
Up: AnX, CKl, AB&GB, AB&GBd, NaCl, Na3AsO4 Dw: CKrl
Fw, Mx, St, Lf, P, Rt, Cl, UnT, Sd, Wp
Up: AnX, CKl, AB&GBd, NaCl Dw: Na3AsO4
Lf, St, P, Mx, Cl, UnT, Wp, Sd, Rt
Up: CKl, AB&GBd, AB&GB, NaCl Dw: AnX, Na3AsO4
Mx, Fw, P, St, UnT
Up: AB&GBd, AB&GB, AnX Dw: Na3AsO4, CKrl
Up: AnX, Na3AsO4, NaCl Dw: CKl
Mx, St, Lf, Fw, UnT, P, Wp, Cl
Up: CKl, AB&GBd, AB&GB Dw: AnX, Na3AsO4, NaCl
Fw, P, UnT, Sd, Wp
Predicted rice GH family 35 β-galactosidase protein properties and locations.
Gol, Vac, ER ext
Gol, PM, IMP, MII, ER ext
Gol, PM, IMP, MII, ER, Sec
Gol, Cis-Gol, MII, ER, Sec
Gol, MII, ER
IMP, MC in, Sec, Cy, N
Gol, MII, ER, Sec
PM, IMP, ER, Sec
Vac-Lys, IMP, MC in, Per, ERl, N, Sec, Cy
Vac-Lys, PM, IMP, MC in, ER, Cy, Sec
Vac-Lys, Gol, PM, IMP, MII, Mic, ER, MC in, Sec
Gol, MII, ER, Sec
IMP, ER, Sec
PM, IMP, MII, ER, Sec
Gol, MII, ER, Sec
Within the plant-type β-galactosidases, OsBgals can be divided into three distinct groups (a, b and c in Figure 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 Asparagus officinalis 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 P. patens Bgals (EDQ59880 and EDQ81397), which suggests this cluster is of ancient origin. Two other putative P. patens 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 , 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.
Chromosomal locations of OsBgals
Exon-intron organization and primary gene structure analysis
Automated annotations for all OsBgal genes were available in the public databases, except for OsBgal6, OsBgal14 and OsBgal15. 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.
The OsBgal2, OsBgal7 and OsBgal9 genes all lack the galactose-binding, lectin-like domain at the C-terminus, however, their splicing patterns were different. OsBgal2 and OsBgal7 have patterns similar to other rice GH35 genes and retain the 5' end of exon 18, but OsBgal2 contains 16 exons with the loss of introns 11 and 16, while intron 15 has been lost in OsBgal7, leaving 17 exons (Figure 3). In contrast, OsBgal9 appeared to be distantly diverged from other rice GH35 members since its exon/intron organization does not correspond to the other OsBgals, except that an intron is found in the same position as intron 1 in the other genes.
Identification of LTR-retrotransposon elements and DNA sequences for interspersed repeats and low complexity DNA sequences within OsBgal6 and OsBgal13.
LTR-retrotransposon elements and DNA sequences for interspersed repeats
Stowaway41_OS, OLO24C, OLO24B, RPO_OS, AT rich Low Complexity, AMYLTP, TA rich Low Complexity, OSTE15
GC rich Low complexity
AT rich Low complexity
OSTE28, OSTE33, CRM-I_OS LTR/Gypsy,
OSLINE1-5, AT rich Low complexity
AT rich Low complexity
RdSpm875A_3.1, ID-2 DNA/Tourist
Comparison of gene splice patterns and phylogenetic relationships between plant-type Bgals from rice, A. thaliana, and P. patens 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 P. patens gene EDQ878340), intron losses appear to have occurred independently in rice, Arabidopsis and P. patens. 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 OsBgal2 and OsBgal7), though it appears to have been lost at least 4 times in Arabidopsis GH35 Bgals .
Analysis of expression by RT-PCR
The rice tissues/organs and developmental stages chosen for RT-PCR analysis are shown in Figure 4. All 15 OsBgal 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 OsBgal13 appears to be expressed in nearly every organ and developmental stage, though at different levels (Figure 4).
Nine OsBgals (OsBgal1, OsBgal2, OsBgal3, OsBgal4, OsBgal6, OsBgal7, OsBgal8 and OsBgal9) displayed their highest transcript levels in young and early-mature stages of shoot, root, and leaf sheath. Among these genes, OsBgal6, OsBgal7, OsBgal8 and OsBgal9 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 . In contrast, low levels of most gene transcripts were detected in leaf blade, except for OsBgal13 transcripts, which were detected throughout the growth of the blade.
Six genes, OsBgal5, OsBgal10, OsBgal11, OsBgal12, OsBgal14, and OsBgal15 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 4). In addition, no expression was detected at the initiating panicle stage, but it was high in later panicle development and flowering stages. Only OsBgal10, OsBgal11 and OsBgal15 transcripts were found to accumulate in anther. These genes appear to have diverged from two ancestor genes, since OsBgal5, OsBgal12, OsBgal14, and OsBgal15 have sequences and gene structures that are similar to each other, as do OsBgal10 and OsBgal11 (Figures 1 and 3). Interestingly, AtBGAL7 and AtBGAL15, which fall into the group b with OsBgal5, OsBgal12, OsBgal14, and OsBgal15 (Figure 1), were reported to function in the early stages of microspore and pollen development . Likewise, AtBGAL11 and AtBGAL13, which fall in cluster c1 with OsBgal10 and OsBgal11 (Figures 1), show maximum expression in mature pollen . So, related rice and Arabidopsis Bgals function in reproductive tissues, suggesting their roles may be somewhat conserved in the two plants.
Only the OsBgal4 transcript was expressed abundantly in seed at the imbibition stage, while OsBgal2 and OsBgal13 were also present at low levels. Five OsBgals were expressed in milk grain and six in grain during dry down. A high abundance of OsBgal2 and OsBgal13 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. OsBgal1, OsBgal2, OsBgal6, OsBgal8, and OsBgal13 transcripts were also found in shoots and roots of seedlings after germination. Chantarangsee et al.  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 (Hordeum vulgare), Giannakouros and colleagues  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.
To complement the RT-PCR expression profiles of the OsBals, ESTs retrieved from the UniGene database (UniGene December, 2007, NCBI) were compared with our findings. Full-length or partial cDNA and ESTs for most OsBgals are available in the database, except OsBgal6, for which only one EST from unspecified tissues is available, and OsBgal15, for which no corresponding ESTs or cDNA are available (Table 1). However, the other OsBgal genes had cDNA and ESTs present in the database, with numbers of ESTs that varied from 1 to 253.
The OsBgal8 (Unigene cluster Os.22360, locus Os03g0255100) 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 1). The Unigene profile is somewhat similar to the RT-PCR pattern of OsBgal8, in which the transcripts accumulated significantly throughout rice plant, including in seed during dry-down, but not in seedling and milk grain.
OsBgal10 (Os.18310), OsBgal11, OsBgal12 (Os.46702) and OsBgal14 (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,OsBgal11 shows 249 TPM in flower and 158 TPM in panicle, while OsBgal14 shows 358 TPM in flower (Table 1).
The EST and cDNA clones for the remaining OsBgals had diverse origins, which is consistent with the evaluation of OsBgals 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. , 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 1 shows that OsBgal expression levels may be affected by many treatments.
Rice β-galactosidase cDNAs sequencing and identification of alternative splicing
Although the expression of all OsBgals 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 OsBgal6, OsBgal11 and OsBgal15. 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 1.
OsBgal11 (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 OsBgal6. A 1343-nucleotide fragment from the 3' half of the OsBgal11 cDNA was amplified from exerting panicle (Genbank Accession EU603286), 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 Penicillium sp. β-galactosidase structure , so it is possible that it could still function as a Bgal.
No experimentally determined cDNA or ESTs corresponding to OsBgal15 were available in public databases, despite the fact that our RT-PCR demonstrated its expression. Therefore, the full-length OsBgal15 cDNA was amplified, cloned and sequenced. The OsBgal15 cDNA sequence (Genbank Accession EU051629) contains a 2487-nucleotide ORF, which encodes 829 amino acids. The encoded protein shares over 99% identity with a predicted β-galactosidase of Oryza sativa (japonica cultivar-group) (accession no. AP004733), except that 7 extra amino acid residues were inserted at residue 293, compared to the predicted sequence from AP004733. Therefore, the correct splice donor site of exon 8 is located 21 nucleotides downstream of that predicted for the mRNA of AP004733.
The alternative splicing observed in this study introduces in-frame stop codons into the mRNA of OsBgal6 in the low expression tissue, but not for the high expression tissues, while it can also introduce a premature stop codon in OsBgal11, 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 . 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 ; diacylglycerol kinase in tomato ; and chloroplast ascorbate peroxidase in spinach . Missplicing events were observed in maize (Zea mays) Viviparous 1 (Vp1) 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 . Although alternative splicing is a widespread process used in higher eukaryotes to regulate gene expression and functional diversification of proteins [50, 51], the putative control of tissue-specific expression by missplicing we observed has not been well characterized in rice.
Preliminary Characterization of OsBgal13
Hydrolysis of p-nitrophenyl glycosides, oligosaccharides and polysaccharides by OsBgal13.
(Percent Relative Activity)
(estimated relative activity)
Arabinogalactan, xylan from oat spelts, birchwood xylan, polygalacturonic acid, apple pectin, citric pectin, galactan
Possible Functions of Other Isozymes
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 . Likewise, β-galacosidase of tomato was reported to act in galactolipid turnover and degradation, which occurs in chloroplasts and chromoplasts during tomato fruit development [54, 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.
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  and pollen . 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.
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 , 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  and an asparagus Bgal thought to act in senescence , 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  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.
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.
Fifteen rice glycosyl hydrolase family 35 genes encoding putative β-galactosidases were identified in this study. This number is similar to the 17 seen in Arabidopsis and these plants appear to contain 9 common gene lineages present in their ancestors before they diverged . OsBgal9 was clustered into the same group as Bgals from animal species, such as H. sapiens and D. melanogaster 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 P. patens, a nonvascular plant, suggests that both types of genes were present early in plant evolution. Within the plant-type Bgals, group a (Figure 1), contains 7 OsBgals and many Bgals that have been characterized from other plants (for example, 9 of Arabidopsis, 4 of C. arietinum and 3 of L. esculentum), including some with and some without C-terminal galactose-binding-lectin-like domains. Many of these proteins, including OsBgal8 , 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. OsBgal6 appeared to undergo alternative splicing to prevent its production in panicle, and alternative splicing was also found for OsBgal11. Understanding the functions of this gene family and significance of alternative splicing in this will require further functional investigation.
Sequence data and database search
To find all GH35 genes in rice, tBLASTn searches  were performed in the National Center for Biotechnology Information (NCBI) Genbank nr, indica rice genome, and expressed sequence tags (dbEST) databases with OsBgal1 (AK102192) as the query. The retrieved OsBgal genes were used to identify their Unigene cluster, Gene Locus, conserved domains, gene position on the 12 Oryza sativa (rice) chromosomes, homolog genes, EST-based expression profiles and GEO profile in the NCBI databases (please see Availability & 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 & requirements for more information) or by BLASTp searches at NCBI. Putative signal peptides were predicted with the SignalP program (please see Availability & requirements for more information) ; and putative N-glycosylation sites were identified with the NetNGlyc ; (please see Availability & requirements for more information) and manually inspected to remove NPS/T sites.
Chromosomal locations of genes were identified on and the map drawn with the NCBI map viewer (please see Availability & requirements for more information). Segmental duplication analysis was done with DAGchainer  and the TIGR rice segmental duplication database (please see Availability & 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 & requirements for more information).
Construction of phylogenetic trees
The multiple sequence alignment of β-galactosidase protein sequences from rice and other organisms was made with ClustalW  and manually adjusted and edited to remove unconserved N- and C-terminal regions with Genedoc (please see Availability & requirements for more information) . Phylogenetic trees were constructed by the neighbor-joining facility in ClustalX 1.83  and the branching order verified for maximum parsimony with the Protpars program in the PHYLIP software suite .
Indica rice (Oryza sativa 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.
Semi-quantitative reverse transcription-polymerase chain reaction (qRT-PCR)
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 ), 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.
Sequences of primers used for semi-quantitative RT-PCR analysis of β-galactosidase gene expression.
Annealing Tm (°C)
Amplification of cDNA clones and sequencing
The first strand cDNA reverse-transcribed from RNA of developing and emerging panicle was used as a template to amplify OsBgal6, OsBgal11, and OsBgal15 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 Pfu 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).
Two pairs of primers were used to amplify OsBgal6 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 OsBgal6 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).
The oligonucleotide primers used for amplification of a partial OsBgal11 cDNA were: bgal11_Seq_r1 (TCATCGCATGTGCAGTG) and bgal11_STOP_R (CCTTCTTCCTAAGCCGCCTG).
The Osbgal15 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 EcoRIsite of pBluescript KS(+) by standard methods . The recombinant plasmid was sequenced with the T7, M13, bgal15QIENE_F, bgal15QIENE_REV, bgal15_1138_F (ACTCATCTTCTGCCTGCTTG), and Bgal15_1797_R (TACGGTGCCATTGTTG TTG) primers.
Expression of OsBgal13
A cDNA encoding the predicted mature OsBgal13 protein was amplified from the full-length cDNA clone of Genbank Accession AK065546, 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 Pfu 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  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. . The protein was extracted and purified on Ni-NTA IMAC resin (QIAGEN), as with OsBgal1 , 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 p-nitrophenyl (pNP) glycosides for 30 min at 30°C, followed by stopping the reaction with 2 volumes 0.4 M Na2CO3 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 .
Availability & requirements
NCBI databases: http://www.ncbi.nlm.nih.gov/sites/entrez
CAZY website: http://www.cazy.org/
SignalP program: http://www.cbs.dtu.dk/services/SignalP/
NCBI map viewer: http://www.ncbi.nlm.nih.gov/projects/mapview
TIGR rice segmental duplication database: http://www.tigr.org/tdb/e2k1/osa1/segmental_dup/500kb/segdup_500kb.shtml
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.
- McCartney L, Ormerod AP, Gidley MJ, Knox JP: Temporal and spatial regulation of pectic (1→4)-β-D-galactan in cell walls of developing pea cotyledons: implications for mechanical properties. Plant J. 2000, 22: 105-113.PubMedView ArticleGoogle Scholar
- Sørensen SO, Pauly M, Bush M, Skjøt M, McCann MC, Borkhardt B, Ulvskov P: Pectin engineering: modification of potato pectin by in vivo expression of an endo-1,4-β-D-galactanase. Proc Natl Acad Sci USA. 2000, 97: 7639-7644.View ArticleGoogle Scholar
- Moctezuma E, Smith DL, Gross KC: Antisense suppression of a β-galactosidase gene (TBG6) in tomato increases fruit cracking. J Exp Bot. 2003, 54 (390): 2025-2033.PubMedView ArticleGoogle Scholar
- Bayless TM, Rosensweig NS: A racial difference in incidence of lactase deficiency. JAMA. 1966, 197: 968-972.PubMedView ArticleGoogle Scholar
- Cuatrecasas P, Lockwood DH, Caldwell JR: Lactase deficiency in the adult: A common occurrence. Lancet. 1965, 1: 14-18.PubMedView ArticleGoogle Scholar
- O'Brien JS: β-Galactosidase deficiency (GM1 gangliosidosis, galactosialidosis, and Morquio syndrome type B); ganglioside sialidase deficiency (mucolipidosis IV). The Metabolic Basis of Inherited Disease. Edited by: Scriver CR, Beaudet AL, Sly WS, Valle D. 1989, New York: McGraw-Hill Book Co, 1797-1806.Google Scholar
- Conzelmann E, Sandhoff K: Glycolipid and glycoprotein degradation. Adv Enzymol Relat Areas Mol Biol. 1987, 60: 89-216.PubMedGoogle Scholar
- Campbell JH, Lengyel JA, Langridge J: Evolution of a second gene for β-galactosidase in Escherichia coli. Proc Natl Acad Sci USA. 1973, 70 (6): 1841-1845.PubMedPubMed CentralView ArticleGoogle Scholar
- Akasaki M, Suzuki M, Funakoshi I, Yamashina I: Characterization of beta-galactosidase from a special strain of Aspergillus oryzae. J Biochem (Tokyo). 1976, 80 (6): 1195-200.Google Scholar
- Henrissat B: A classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem J. 1991, 280: 309-316.PubMedPubMed CentralView ArticleGoogle Scholar
- Henrissat B, Bairoch A: New families in the classification of glycosyl hydrolases based on amino-acid sequence similarities. Biochem J. 1993, 293: 781-788.PubMedPubMed CentralView ArticleGoogle Scholar
- Ahn YO, Zheng M, Bevan DR, Esen A, Shiu SH, Benson J, Peng HP, Miller JT, Cheng LI, Poulton JE, Shih MC: Functional genomic analysis of Arabidopsis thaliana glycoside hydrolase family 35. Phytochemistry. 2007, 68 (11): 1510-1520.PubMedView ArticleGoogle Scholar
- Smith DL, Gross KC: A family of at least seven β-galactosidase genes is expressed during tomato fruit development. Plant Physiol. 2000, 23: 1173-1183.View ArticleGoogle Scholar
- Lazan H, Ng SY, Goh LY, Ali ZM: Papaya beta-galactosidase/galactanase isoforms in differential cell wall hydrolysis and fruit softening during ripening. Plant Physiol Biochem. 2004, 42 (11): 847-853.PubMedView ArticleGoogle Scholar
- Ross CS, Cavin S, Wegrzyn T, MacRae EA, Redgwell RJ: Apple β-galactosidase activity against cell wall polysaccharides and characterization of a related cDNA clone. Plant Physiol. 1994, 106: 521-528.PubMedPubMed CentralView ArticleGoogle Scholar
- Enéas-Filho J, Sudério FB, Gomes-Filho E, Prisco JT: Multiple forms of cotyledonary β-galactosidases from Vigna unguiculata quiescent seeds. Rev Bras Bot. 2000, 23: 69-76.View ArticleGoogle Scholar
- Giannakouros T, Karagiorgos A, Simos G: Expression of β-galactosidase multiple forms during barley (Hordeum vulgare) seed germination. Separation and characterization of enzyme isoforms. Physiol Plant. 1991, 82: 413-418.View ArticleGoogle Scholar
- Kotake T, Dina S, Konishi T, Kaneko S, Igarashi K, Samejima M, Watanabe Y, Kimura K, Tsumuraya Y: Molecular cloning of a β-galactosidase from radish that specifically hydrolyzes β-(1→3)- and β-(1→6)-galactosyl residues of arabinogalactan protein. Plant Physiol. 2005, 138: 1563-1576.PubMedPubMed CentralView ArticleGoogle Scholar
- Sekimata M, Ogura K, Tsumuraya Y, Hashimoto Y, Yamamoto S: A β-galactosidase from radish (Raphanus sativus L.) Seeds. Plant Physiol. 1989, 90: 567-574.PubMedPubMed CentralView ArticleGoogle Scholar
- Hirano Y, Tsumuraya Y, Hashimoto Y: Characterization of spinach leaf α-L-arabinofuranosidases and β-galactosidases and their synergistic action on an endogenous arabinogalactan-protein. Physiol Plant. 1994, 92 (2): 286-296.View ArticleGoogle Scholar
- Pressey R: β-Galactosidases in ripening tomatoes. Plant Physiol. 1983, 71: 132-135.PubMedPubMed CentralView ArticleGoogle Scholar
- Carey AT, Holt K, Picard S, Wilde R, Tucker GA, Bird CR, Schuch W, Seymour GR: Tomato exo-(1→4)-β-D-galactanase: Isolation, changes during ripening in normal and mutant tomato fruit, and characterization of a related cDNA clone. Plant Physiol. 1995, 108: 1099-1107.PubMedPubMed CentralView ArticleGoogle Scholar
- Ranwala AP, Suematsu C, Masuda H: The role of β-galactosidases in the modification of cell wall components during muskmelon fruit ripening. Plant Physiol. 1992, 100: 1318-1325.PubMedPubMed CentralView ArticleGoogle Scholar
- Ross GS, Redgwell RJ, Macrae EA: beta-galactosidase: Isolation and activity against specific fruit cell-wall polysaccharides. Planta. 1993, 189: 499-506.View ArticleGoogle Scholar
- Kang IK, Suh SG, Gross KC, Byun JK: N-terminal amino acid sequence of persimmon fruit β-galactosidase. Plant Physiol. 1994, 105: 975-979.PubMedPubMed CentralView ArticleGoogle Scholar
- Andrews PK, Li S: Partial purification and characterization of beta-D-galactosidase from sweet cherry, a nonclimacteric fruit. J Agric Food Chem. 1994, 2 (10): 2177-2182.View ArticleGoogle Scholar
- Ali ZM, Armugam S, Lazan H: β-Galactosidase and its significance in ripening mango fruit. Phytochemistry. 1995, 38 (5): 1109-1114.PubMedView ArticleGoogle Scholar
- Lee DH, Kang S-G, Suh S-G, Byun JK: Purification and characterization of a β-galactosidase from peach (Prunus persica). Mol Cells. 2003, 15 (1): 68-74.PubMedGoogle Scholar
- Li S-C, Han J-W, Chen K-C, Chen C-S: Purification and characterization of isoforms of β-galactosidases in mung bean seedlings. Phytochemistry. 2001, 57 (3): 349-359.PubMedView ArticleGoogle Scholar
- Wu A, Liu J: Isolation of the promoter of a cotton beta-galactosidase gene (GhGal1) and its expression in transgenic tobacco plants. Sci China C Life Sci. 2006, 49 (2): 105-114.PubMedView ArticleGoogle Scholar
- Hrubá P, Honys D, Twell D, Capková V, Tupy J: Expression of β-galactosidase and β-xylosidase genes during microspore and pollen development. Planta. 2005, 220: 931-940.PubMedView ArticleGoogle Scholar
- Edwards M, Bowman YL, Deagll ICM, Reid JSG: A β-D-Galactosidase from Nasturtium (Tropueolum mujus L.) cotyledons: Purification, properties, and demonstration that xyloglucan is the natural Substrate. J Biol Chem. 1988, 263 (9): 4333-4337.PubMedGoogle Scholar
- Alcantara PHN, Dietrich SMC, Buckeridge MS: Xyloglucan mobilisation and purification of a (XLLG/XLXG) specific-galactosidase from cotyledons of Copaifera langsdorffii. Plant Physiol Biochem. 1999, 37 (9): 653-663.View ArticleGoogle Scholar
- Tiné MA, Cortelazzo AL, Buckeridge MS: Xyloglucan mobilisation in cotyledons of developing plantlets of Hymenaea courbaril L. (Leguminosae caesalpinoideae). Plant Sci. 2000, 154 (2): 117-126.PubMedView ArticleGoogle Scholar
- Tiné MA, Silva CO, Lima DU, Carpita NC, Buckeridge MS: Fine structure of a mixed-oligomer storage xyloglucan from seeds of Hymenaea courbaril. Carbohydr Polym. 2006, 66: 444-454.View ArticleGoogle Scholar
- Konno H, Tsumuki H: Purification of a β-galactosidase from rice shoots and its involvement in hydrolysis of the natural substrate in cell walls. Physiol Plant. 1993, 89: 40-47.View ArticleGoogle Scholar
- Kaneko S, Kobayashi H: Purification and characterization of extracellular β-galactosidase secreated by suspension cultured rice (Oryza sativa L.) cells. Biosci Biotechnol Biochem. 2003, 67 (3): 627-630.PubMedView ArticleGoogle Scholar
- Chantarangsee M, Tanthanuch W, Fujimura T, Fry SC, Cairns JK: Molecular characterization of β-galactosidases from germinating rice (Oryza sativa). Plant Sci. 2007, 173 (2): 118-134.View ArticleGoogle Scholar
- Henrissat B: Glycosidase families. Biochem Soc Trans. 1998, 26: 153-156.PubMedView ArticleGoogle Scholar
- Guo X, Xu G, Zhang Y, Wen X, Hu W, Fan L: Incongruent evolution of chromosomal size in rice. Genet Mol Res. 2006, 5 (2): 373-389.PubMedGoogle Scholar
- Yoshida S, Aha SB: The accumulation process of carbohydrate in rice varieties in relation to their response to nitrogen in the tropics. Soil Sci Plant Nutr. 1968, 14: 153-162.View ArticleGoogle Scholar
- Ma L, Chen C, Liu X, Jiao Y, Su N, Li L, Wang X, Cao M, Sun N, Zhang X, Bao J, Li J, Pedersen S, Bolund L, Zhao H, Yuan L, Wong GK, Wang J, Deng XW, Wang J: A microarray analysis of the rice transcriptome and its comparison to Arabidopsis. Genome Res. 2005, 15: 1274-1283.PubMedPubMed CentralView ArticleGoogle Scholar
- Campbell MA, Haas BJ, Hamilton JP, Mount SM, Buell CR: Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis. BMC Genomics. 2006, 7: 327-PubMedPubMed CentralView ArticleGoogle Scholar
- Rojas AL, Nagem RA, Neustroev KN, Arand M, Adamska M, Eneyskaya EV, Kulminskaya AA, Garratt RC, Golubev AM, Polikarpov I: Crystal structures of beta-galactosidase from Penicillium sp. and its complex with galactose. J Mol Biol. 2004, 343 (5): 1281-1292.PubMedView ArticleGoogle Scholar
- Sharp PA: Speculations on RNA splicing. Cell. 1981, 23: 643-646.PubMedView ArticleGoogle Scholar
- Zhang N, Portis AR: Mechanism of light regulation of Rubisco: A specific role for the larger Rubisco activase isoform involving reductive activation by thioredoxin-f. Proc Natl Acad Sci USA. 1999, 96: 9438-9443.PubMedPubMed CentralView ArticleGoogle Scholar
- Snedden WA, Blumwald E: Alternative splicing of a novel diacylglycerol kinase in tomato leads to a calmodulin-binding isoform. Plant J. 2000, 24: 317-326.PubMedView ArticleGoogle Scholar
- Yoshimura K, Yabuta Y, Tamoi M, Ishikawa T, Shigeoka S: Alternatively spliced mRNA variants of chloroplast ascorbate peroxidase isoenzymes in spinach leaves. Biochem J. 1999, 338: 41-48.PubMedPubMed CentralView ArticleGoogle Scholar
- McKibbin RS, Wilkinson MD, Bailey PC, Flintham JE, Andrew LM, Lazzeri PA, Gale MD, Lenton JR, Holdsworth MJ: Transcripts of Vp-1 homeologues are misspliced in modern wheat and ancestral species. Proc Natl Acad Sci USA. 2002, 99 (15): 10203-10208.PubMedPubMed CentralView ArticleGoogle Scholar
- Graveley BR: Alternative splicing: increasing diversity in the proteomic World. Trends Genet. 2001, 17: 100-107.PubMedView ArticleGoogle Scholar
- Faustino NA, Cooper TA: Pre-mRNA splicing and human disease. Genes Dev. 2003, 17: 419-437.PubMedView ArticleGoogle Scholar
- Carpita NC: Structure and biogenesis of the cell walls of grasses. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 445-476.PubMedView ArticleGoogle Scholar
- Bhalla PL, Dalling MJ: Characteristics of a β-galactosidase associated with the stroma of chloroplasts prepared from mesophyll protoplasts of the primary leaf of wheat. Plant Physiol. 1984, 76: 92-95.PubMedPubMed CentralView ArticleGoogle Scholar
- Güçlü J, Paulin A, Soudain P: Changes in polar lipids during ripening and senescence of cherry tomato (Lycopersicon esculentum): relation to climacteric and ethylene increases. Physiol Plant. 1989, 77: 413-419.View ArticleGoogle Scholar
- Whitaker BD: Changes in galactolipid and phospholipid levels of tomato fruits stored at chilling and nonchilling temperatures. Phytochem. 1992, 31: 2627-2630.View ArticleGoogle Scholar
- King GA, Davies KM, Stewart RJ, Borst WM: Similarities in gene expression during the postharvest-induced senescence of spears and natural foliar senescence of Asparagus. Plant Physiol. 1995, 108: 125-128.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402.PubMedPubMed CentralView ArticleGoogle Scholar
- Bendtsen JD, Nielsen H, Heijne Gv, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795.PubMedView ArticleGoogle Scholar
- Gupta R, Jung E, Brunak S: Prediction of N-glycosylation sites in human proteins. preparation. 2004Google Scholar
- Haas BJ, Delcher AL, Wortman JR, Salzberg SL: DAGchainer: a tool for mining segmental genome duplications and synteny. Bioinformatics. 2004, 20: 3643-3646.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicholas KB, Nicholas HB, Deerfield DW: GeneDoc: Analysis and visualization of genetic variation. EMBNEW News. 1997, 4: 14-Google Scholar
- Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ: Multiple sequence alignment with Clustal X. Trends Biochem Sci. 1998, 23: 403-405.PubMedView ArticleGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.6. 2005, Distributed by the author, Department of Genome Sciences, University of Washington, SeattleGoogle Scholar
- Chomczynski P, Sacchi N: Signal-step method of RNA isolation by acid guanidinium thiocyanate -phenol -chloroform extraction. Anal Biochem. 1987, 162: 156-159.PubMedView ArticleGoogle Scholar
- Sambrook J, Fitsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, Cold Spring Harbor, Cold Spring Harbor PressGoogle Scholar
- Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H, Hotta I, Kojima K, Namiki T, Ohneda E, Yahagi W, Suzuki K, Li CJ, Ohtsuki K, Shishiki T, Otomo Y, Murakami K, Iida Y, Sugano S, Fujimura T, Suzuki Y, Tsunoda Y, Kurosaki T, Kodama T, Masuda H, Kobayashi M, Xie Q, Lu M, Narikawa R, Sugiyama A, Mizuno K, Yokomizo S, Niikura J, Ikeda R, Ishibiki J, Kawamata M, Yoshimura A, Miura J, Kusumegi T, Oka M, Ryu R, Ueda M, Matsubara K, Kawai J, Carninci P, Adachi J, Aizawa K, Arakawa T, Fukuda S, Hara A, Hashizume W, Hayatsu N, Imotani K, Ishii Y, Itoh M, Kagawa I, Kondo S, Konno H, Miyazaki A, Osato N, Ota Y, Saito R, Sasaki D, Sato K, Shibata K, Shinagawa A, Shiraki T, Yoshino M, Hayashizaki Y, Yasunishi A: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301 (5631): 376-379.PubMedView ArticleGoogle Scholar
- Opassiri R, Pomthong B, Onkoksoong T, Akiyama T, Esen A, Ketudat Cairns JR: Analysis of rice glycosyl hydrolase family 1 and expression of Os4bglu12 β-glucosidase. BMC Plant Biol. 2006, 6: 33-PubMedPubMed CentralView ArticleGoogle Scholar