Lotus Cell Walls and the Genes Involved in its Synthesis and Modification

The lotus genome (Nelumbo nucifera (Gaertn.)) lacks the paleo-triplication found in other eudicots and has evolved remarkably slowly with fewer nucleotide mutations. It is thought to have greater retention of duplicated genes than other angiosperms. We evaluated the potential genes involved in cell wall synthesis and its modification, and ethylene synthesis and response. In many cell wall transferases and hydrolases families, lotus had fewer members in most families when compared to Arabidopsis. Lotus had similar or fewer members in each family as found in poplar, grape and papaya. The exceptions were in the sialyl and beta-glucuronsyl transferases where similar number were found as in the core eudicots. Lotus had similar numbers of polygalacturonase and pectin methyl esterases as found in Arabidopsis but fewer in all other hydrolases families. For starch degradation, lotus had only two alpha amylases predicted genes versus eight to ten in other eudicots, with similar numbers of beta amylase genes predicted. Lotus also had less than half the number of genes predicted for the enzymes involved in lignin and tannin synthesis compared to Arabidopsis. The stress plant growth regulator ethylene’s synthesis, reception and response predicted genes were fewer in lotus than other eudicots. Only two ethylene receptor genes were predicted in lotus with five reported for Arabidopsis and six for tomato. Our analysis does not supports the conclusion that this species has greater retention of duplicated genes though our data does support the conclusion that lotus split occurred at the base of the eudicots.     

The glycosyltransferase repertoire of the spikemoss Selaginella moellendorffii and a comparative study of its cell wall

Spike mosses are among the most basal vascular plants, and one species, Selaginella moellendorffii, was recently selected for full genome sequencing by the Joint Genome Institute (JGI). Glycosyltransferases (GTs) are involved in many aspects of a plant life, including cell wall biosynthesis, protein glycosylation, primary and secondary metabolism. Here, we present a comparative study of the S. moellendorffii genome across 92 GT families and an additional family (DUF266) likely to include GTs. The study encompasses the moss Physcomitrella patens, a non-vascular land plant, while rice and Arabidopsis represent commelinid and non-commelinid seed plants. Analysis of the subset of GT-families particularly relevant to cell wall polysaccharide biosynthesis was complemented by a detailed analysis of S. moellendorffii cell walls. The S. moellendorffii cell wall contains many of the same components as seed plant cell walls, but appears to differ somewhat in its detailed architecture. The S. moellendorffii genome encodes fewer GTs (287 GTs including DUF266s) than the reference genomes. In a few families, notably GT51 and GT78, S. moellendorffii GTs have no higher plant orthologs, but in most families S. moellendorffii GTs have clear orthologies with Arabidopsis and rice. A gene naming convention of GTs is proposed which takes orthologies and GT-family membership into account. The evolutionary significance of apparently modern and ancient traits in S. moellendorffii is discussed, as is its use as a reference organism for functional annotation of GTs.     

Integrating genes and phenotype: a wheat–Arabidopsis–rice glycosyltransferase database for candidate gene analyses

Glycosyltransferases (GTs) constitute a very large multi-gene superfamily, containing several thousand members identified in sequenced organisms especially in plants. GTs are key enzymes involved in various biological processes such as cell wall formation, storage polysaccharides biosynthesis, and glycosylation of various metabolites. GTs have been identified in rice (Oryza sativa) and Arabidopsis thaliana, but their precise function has been demonstrated biochemically for only a few. In this work we have established a repertoire of virtually all the wheat (Triticum aestivum) GT sequences, using the large publicly available banks of expressed sequences. Based on sequence similarity with Arabidopsis and rice GTs compiled in the carbohydrate active enzyme database (CAZY), we have identified and classified these wheat sequences. The results were used to feed a searchable database available on the web () that can be used for initiating an exhaustive candidate gene survey in wheat applied to a particular biological process. This is illustrated through the identification of GT families which are expressed during cell wall formation in wheat grain maturation. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. This work was funded by a grant of the French ministry of research.     

Evolutionary history of the GH3 family of acyl adenylases in rosids

GH3 amino acid conjugases have been identified in many plant and bacterial species. The evolution of GH3 genes in plant species is explored using the sequenced rosids Arabidopsis, papaya, poplar, and grape. Analysis of the sequenced non-rosid eudicots monkey flower and columbine, the monocots maize and rice, as well as spikemoss and moss is included to provide further insight into the origin of GH3 clades. Comparison of co-linear genes in regions surrounding GH3 genes between species helps reconstruct the evolutionary history of the family. Combining analysis of synteny with phylogenetics, gene expression and functional data redefines the Group III GH3 genes, of which AtGH3.12/PBS3, a regulator of stress-induced salicylic acid metabolism and plant defense, is a member. Contrary to previous reports that restrict PBS3 to Arabidopsis and its close relatives, PBS3 syntelogs are identified in poplar, grape, columbine, maize and rice suggesting descent from a common ancestral chromosome dating to before the eudicot/monocot split. In addition, the clade containing PBS3 has undergone a unique expansion in Arabidopsis, with expression patterns for these genes consistent with specialized and evolving stress-responsive functions.     

Carbohydrate-active enzymes involved in the secondary cell wall biogenesis in hybrid aspen

Wood formation is a fundamental biological process with significant economic interest. While lignin biosynthesis is currently relatively well understood, the pathways leading to the synthesis of the key structural carbohydrates in wood fibers remain obscure. We have used a functional genomics approach to identify enzymes involved in carbohydrate biosynthesis and remodeling during xylem development in the hybrid aspen Populus tremula x tremuloides. Microarrays containing cDNA clones from different tissue-specific libraries were hybridized with probes obtained from narrow tissue sections prepared by cryosectioning of the developing xylem. Bioinformatic analyses using the sensitive tools developed for carbohydrate-active enzymes allowed the identification of 25 xylem-specific glycosyltransferases belonging to the Carbohydrate-Active EnZYme families GT2, GT8, GT14, GT31, GT43, GT47, and GT61 and nine glycosidases (or transglycosidases) belonging to the Carbohydrate-Active EnZYme families GH9, GH10, GH16, GH17, GH19, GH28, GH35, and GH51. While no genes encoding either polysaccharide lyases or carbohydrate esterases were found among the secondary wall-specific genes, one putative O-acetyltransferase was identified. These wood-specific enzyme genes constitute a valuable resource for future development of engineered fibers with improved performance in different applications.     

Comparative analysis of GT14/GT14-like gene family in Arabidopsis, Oryza, Populus, Sorghum and Vitis

Glycosyltransferase family14 (GT14) belongs to the glycosyltransferase (GT) superfamily that plays important roles in the biosynthesis of cell walls, the most abundant source of cellulosic biomass for bioethanol production. It has been hypothesized that DUF266 proteins are a new class of GTs related to GT14. In this study, we identified 62 GT14 and 106 DUF266 genes (named GT14-like herein) in Arabidopsis, Oryza, Populus, Sorghum and Vitis. Our phylogenetic analysis separated GT14 and GT14-like genes into two distinct clades, which were further divided into eight and five groups, respectively. Similarities in protein domain, 3D structure and gene expression were uncovered between the two phylogenetic clades, supporting the hypothesis that GT14 and GT14-like genes belong to one family. Therefore, we proposed a new family name, GT14/GT14-like family that combines both subfamilies. Variation in gene expression and protein subcellular localization within the GT14-like subfamily were greater than those within the GT14 subfamily. One-half of the Arabidopsis and Populus GT14/GT14-like genes were found to be preferentially expressed in stem/xylem, indicating that they are likely involved in cell wall biosynthesis. This study provided new insights into the evolution and functional diversification of the GT14/GT14-like family genes.     

A collection of glycosyltransferases from rice (Oryza sativa) exposed to atrazine

The rice (Oryza sativa) GTs belong to a super family possibly with hundreds of members. However, which GTs are involved in plant response to toxic chemicals is unknown. Here, we demonstrated 59 novel GT genes screened from our recent genome-wide sequencing datasets of rice crops exposed to atrazine (a herbicide persistent in ecosystems). Analysis of GT genes showed that most of the GTs contain functional domains typically found in proteins transferring glycosyl moieties to their target compounds. A phylogenetic analysis revealed that many GT genes from different families have diverse cis-elements necessary for response to biotic and environmental stresses. Experimental validation for the GTs was undertaken through a microarray, and 36 GT genes were significantly detected with an expression pattern similar to that from deep-sequencing datasets. Furthermore, 12 GT genes were randomly selected and confirmed by quantitative real-time RT-PCR. Finally, the special activity of total GTs was determined in rice roots and shoots, with an increased activity under the atrazine exposure. This response was closely associated with atrazine absorption in the rice tissues. These results indicate that exposure to atrazine can trigger specific GT genes and enzyme activities in rice.     

Fruit Development, Ripening and Quality Related Genes in the Papaya Genome

Papaya (Carica papaya L.) is the first fleshy fruit with a climacteric ripening pattern to be sequenced. As a member of the Rosids superorder in the order Brassicales, papaya apparently lacks the genome duplication that occurred twice in Arabidopsis. The predicted papaya genes that are homologous to those potentially involved in fruit growth, development, and ripening were investigated. Genes homologous to those involved in tomato fruit size and shape were found. Fewer predicted papaya expansin genes were found and no Expansin Like-B genes were predicted. Compared to Arabidopsis and tomato, fewer genes that may impact sugar accumulation in papaya, ethylene synthesis and response, respiration, chlorophyll degradation and carotenoid synthesis were predicted. Similar or fewer genes were found in papaya for the enzymes leading to volatile production than so far determined for tomato. The presence of fewer papaya genes in most fruit development and ripening categories suggests less subfunctionalization of gene action. The lack of whole genome duplication and reductions in most gene families and biosynthetic pathways make papaya a valuable and unique tool to study the evolution of fruit ripening and the complex regulatory networks active in fruit ripening.     

The plant glycosyltransferase clone collection for functional genomics

The glycosyltransferases (GTs) are an important and functionally diverse family of enzymes involved in glycan and glycoside biosynthesis. Plants have evolved large families of GTs which undertake the array of glycosylation reactions that occur during plant development and growth. Based on the Carbohydrate‐Active enZymes (CAZy) database, the genome of the reference plant Arabidopsis thaliana codes for over 450 GTs, while the rice genome (Oryza sativa) contains over 600 members. Collectively, GTs from these reference plants can be classified into over 40 distinct GT families. Although these enzymes are involved in many important plant specific processes such as cell‐wall and secondary metabolite biosynthesis, few have been functionally characterized. We have sought to develop a plant GTs clone resource that will enable functional genomic approaches to be undertaken by the plant research community. In total, 403 (88%) of CAZy defined Arabidopsis GTs have been cloned, while 96 (15%) of the GTs coded by rice have been cloned. The collection resulted in the update of a number of Arabidopsis GT gene models. The clones represent full‐length coding sequences without termination codons and are Gateway^® compatible. To demonstrate the utility of this JBEI GT Collection, a set of efficient particle bombardment plasmids (pBullet) was also constructed with markers for the endomembrane. The utility of the pBullet collection was demonstrated by localizing all members of the Arabidopsis GT14 family to the Golgi apparatus or the endoplasmic reticulum (ER). Updates to these resources are available at the JBEI GT Collection website http://www.addgene.org/ .     

Composition and Expression of Genes Encoding Carbohydrate-Active Enzymes in the Straw-Degrading Mushroom Volvariella volvacea

Volvariella volvacea is one of a few commercial cultivated mushrooms mainly using straw as carbon source. In this study, the genome of V. volcacea was sequenced and assembled. A total of 285 genes encoding carbohydrate-active enzymes (CAZymes) in V. volvacea were identified and annotated. Among 15 fungi with sequenced genomes, V. volvacea ranks seventh in the number of genes encoding CAZymes. In addition, the composition of glycoside hydrolases in V. volcacea is dramatically different from other basidiomycetes: it is particularly rich in members of the glycoside hydrolase families GH10 (hemicellulose degradation) and GH43 (hemicellulose and pectin degradation), and the lyase families PL1, PL3 and PL4 (pectin degradation) but lacks families GH5b, GH11, GH26, GH62, GH93, GH115, GH105, GH9, GH53, GH32, GH74 and CE12. Analysis of genome-wide gene expression profiles of 3 strains using 3′-tag digital gene expression (DGE) reveals that 239 CAZyme genes were expressed even in potato destrose broth medium. Our data also showed that the formation of a heterokaryotic strain could dramatically increase the expression of a number of genes which were poorly expressed in its parental homokaryotic strains.     

保护性耕作对东北黑土微生物碳循环功能基因的影响

土壤微生物既参与土壤有机碳的分解也是土壤有机碳转化和固定的驱动者,是影响土壤碳循环和有机碳稳定性的关键因素。然而,保护性耕作（秸秆还田）如何通过调节土壤微生物碳循环功能基因组成来影响土壤CO₂释放的机理尚不明确。因此,依托中国科学院东北地理与农业生态研究所长春保护性耕作观测站,借助鸟枪法宏基因组测序技术,分析了玉米连作系统不同耕作方式（免耕（NT）、秋翻（MP）以及常规耕作（CT））对土壤CO₂释放速率、碳水化合物活性酶（CAZy）、碳循环功能基因（碳固定、甲烷代谢以及碳水化合物代谢）组成的影响。研究表明：基于生长季节土壤CO₂释放速率6年平均值分析发现,生长季前期免耕土壤的平均CO₂释放速率显著低于秋翻和常规耕作,分别比秋翻低28%（5月份）、11%（6月份）和23%（7月份）;比常规耕作低31%（5月份）、19%（6月份）和7%（7月份）。基于CAZy数据库注释结果,发现耕作处理显著影响一些糖苷水解酶（如GH102、GH5_38和GH13_17）、糖基转移酶（如GT39）和多糖裂解酶（如PL17和PL5_1）的基因丰度,与常规耕作相比,秸秆还田的免耕和秋翻处理的这些差异基因的相对丰度较高。基于京都基因与基因组百科全书（KEGG）数据库注释结果,发现耕作方式显著影响土壤碳循环功能基因组成（Adonis,多元方差分析,R²=0.45;P=0.006）,且免耕处理土壤的碳固定、甲烷代谢以及碳水化合物代谢功能基因组成不同于常规耕作和秋翻处理,单独聚为一类。免耕土壤上调的碳固定功能基因的相对丰度（所有上调功能基因相对丰度的平均值）分别比常规耕作和秋翻高17%和11%,而下调的2个功能基因（K01007和K00170）的丰度分别低19%（CT）、21%（MP）和14%（CT）、17%（MP）。免耕土壤上调的甲烷代谢基因相对丰度分别较常规耕作和秋翻高15%和10%;下调基因的丰度分别低13%（CT）和11%（MP）。免耕土壤上调的碳水化合物代谢功能基因丰度较常规耕作和秋翻高23%和14%;下调的基因丰度分别低25%（CT）和18%（MP）。冗余分析（db-RDA）表明土壤容重及土壤水溶性有机碳（DOC）是驱动土壤碳循环功能基因组成差异的主要因子（P<0.05）,且免耕土壤上调的碳固定功能基因（K00625、K01676、K09709、K00925和K14470等）、甲烷代谢基因（K03520、K00830、K10713、K15633和K00625等）和碳水化合物代谢功能基因（K00886、K00830、K01676、K00117和K00114等）与土壤DOC、容重或含水量呈显著正相关。此外,研究发现土壤CO₂释放速率与土壤碳循环功能基因组成显著相关（R²=0.80;P<0.01）,尤其是与一些碳水化合物代谢功能基因显著相关。这些结果说明免耕处理通过影响土壤理化性质改变土壤碳循环过程,且推断免耕秸秆还田和减少干扰的叠加效应通过调节碳循环功能基因组成来提高土壤固碳潜力。     

Auxin-related gene families in abiotic stress response in Sorghum bicolor

Sorghum, a C4 model plant, has been studied to develop an understanding of the molecular mechanism of resistance to stress. The auxin-response genes, auxin/indole-3-acetic acid (Aux/IAA), auxin-response factor (ARF), Gretchen Hagen3 (GH3), small auxin-up RNAs, and lateral organ boundaries (LBD), are involved in growth/development and stress/defense responses in Arabidopsis and rice, but they have not been studied in sorghum. In the present paper, the chromosome distribution, gene duplication, promoters, intron/exon, and phylogenic relationships of Aux/IAA, ARF, GH3, and LBD genes in sorghum are presented. Furthermore, real-time PCR analysis demonstrated these genes are differently expressed in leaf/root of sorghum and indicated the expression profile of these gene families under IAA, brassinosteroid (BR), salt, and drought treatments. The SbGH3 and SbLBD genes, expressed in low level under natural condition, were highly induced by salt and drought stress consistent with their products being involved in both abiotic stresses. Three genes, SbIAA1, SbGH3-13, and SbLBD32, were highly induced under all the four treatments, IAA, BR, salt, and drought. The analysis provided new evidence for role of auxin in stress response, implied there are cross talk between auxin, BR and abiotic stress signaling pathways.     

Genome of papaya, a fast growing tropical fruit tree

Papaya is a major fruit crop in tropical and subtropical regions worldwide. It has long been recognized as a nutritious and healthy fruit rich in vitamins A and C. Its small genome, unique aspects of nascent sex chromosomes, and agricultural importance are justifications for sequencing the genome. A female plant of the transgenic variety SunUp was selected for sequencing its genome to avoid the complication of assembling the XY chromosomes in a male or hermaphrodite plant. The draft genome revealed fewer genes than sequenced genomes of flowering plants, partly due to its lack of genome wide duplication since the ancient triplication event shared by eudicots. Most gene families have fewer members in papaya, including significantly fewer disease resistance genes. However, striking amplifications in gene number were found in some functional groups, including MADS-box genes, starch synthases, and volatiles that might affect the speciation and adaptation of papaya. The draft genome was used to clone a gene controlling fruit flesh color and to accelerate the construction of physical maps of sex chromosomes in papaya. Finishing the papaya genome and re-sequencing selected genomes in the family will further facilitate papaya improvement and the study of genome and sex chromosome evolution in angiosperms, particularly in Caricaceae.     

A novel bioinformatics approach identifies candidate genes for the synthesis and feruloylation of arabinoxylan

Arabinoxylans (AXs) are major components of graminaceous plant cell walls, including those in the grain and straw of economically important cereals. Despite some recent advances in identifying the genes encoding biosynthetic enzymes for a number of other plant cell wall polysaccharides, the genes encoding enzymes of the final stages of AX synthesis have not been identified. We have therefore adopted a novel bioinformatics approach based on estimation of differential expression of orthologous genes between taxonomic divisions of species. Over 3 million public domain cereal and dicot expressed sequence tags were mapped onto the complete sets of rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana) genes, respectively. It was assumed that genes in cereals involved in AX biosynthesis would be expressed at high levels and that their orthologs in dicotyledonous plants would be expressed at much lower levels. Considering all rice genes encoding putative glycosyl transferases (GTs) predicted to be integral membrane proteins, genes in the GT43, GT47, and GT61 families emerged as much the strongest candidates. When the search was widened to all other rice or Arabidopsis genes predicted to encode integral membrane proteins, cereal genes in Pfam family PF02458 emerged as candidates for the feruloylation of AX. Our analysis, known activities, and recent findings elsewhere are most consistent with genes in the GT43 families encoding beta-1,4-xylan synthases, genes in the GT47 family encoding xylan alpha-1,2- or alpha-1,3-arabinosyl transferases, and genes in the GT61 family encoding feruloyl-AX beta-1,2-xylosyl transferases.     

Identification of Putative Rhamnogalacturonan-II Specific Glycosyltransferases in Arabidopsis Using a Combination of Bioinformatics Approaches

Rhamnogalacturonan-II (RG-II) is a complex plant cell wall polysaccharide that is composed of an α(1,4)-linked homogalacturonan backbone substituted with four side chains. It exists in the cell wall in the form of a dimer that is cross-linked by a borate di-ester. Despite its highly complex structure, RG-II is evolutionarily conserved in the plant kingdom suggesting that this polymer has fundamental functions in the primary wall organisation. In this study, we have set up a bioinformatics strategy aimed at identifying putative glycosyltransferases (GTs) involved in RG-II biosynthesis. This strategy is based on the selection of candidate genes encoding type II membrane proteins that are tightly coexpressed in both rice and Arabidopsis with previously characterised genes encoding enzymes involved in the synthesis of RG-II and exhibiting an up-regulation upon isoxaben treatment. This study results in the final selection of 26 putative Arabidopsis GTs, including 10 sequences already classified in the CAZy database. Among these CAZy sequences, the screening protocol allowed the selection of α-galacturonosyltransferases involved in the synthesis of α4-GalA oligogalacturonides present in both homogalacturonans and RG-II, and two sialyltransferase-like sequences previously proposed to be involved in the transfer of Kdo and/or Dha on the pectic backbone of RG-II. In addition, 16 non-CAZy GT sequences were retrieved in the present study. Four of them exhibited a GT-A fold. The remaining sequences harbored a GT-B like fold and a fucosyltransferase signature. Based on homologies with glycosyltransferases of known functions, putative roles in the RG-II biosynthesis are proposed for some GT candidates.     

Family 6 Glycosyltransferases in Vertebrates and Bacteria: Inactivation and Horizontal Gene Transfer May Enhance Mutualism between Vertebrates and Bacteria

Glycosyltransferases (GTs) control the synthesis and structures of glycans. Inactivation and intense allelic variation in members of the GT6 family generate species-specific and individual variations in carbohydrate structures, including histo-blood group oligosaccharides, resulting in anti-glycan antibodies that target glycan-decorated pathogens. GT6 genes are ubiquitous in vertebrates but are otherwise rare, existing in a few bacteria, one protozoan, and cyanophages, suggesting lateral gene transfer. Prokaryotic GT6 genes correspond to one exon of vertebrate genes, yet their translated protein sequences are strikingly similar. Bacterial and phage GT6 genes influence the surface chemistry of bacteria, affecting their interactions, including those with vertebrate hosts.     

Comparison of broad-scope assays of nucleotide sugar-dependent glycosyltransferases

Glycosyltransferases (GTs) are abundant in nature and diverse in their range of substrates. Application of GTs is, however, often complicated by their narrow substrate specificity. GTs with tailored specificities are highly demanded for targeted glycosylation reactions. Engineering of such GTs is, however, restricted by lack of practical and broad-scope assays currently available. Here we present an improvement of an inexpensive and simple assay that relies on the enzymatic detection of inorganic phosphate cleaved from nucleoside phosphate products released in GT reactions. This phosphatase-coupled assay (PCA) is compared with other GT assays: a pH shift assay and a commercially available immunoassay in Escherichia coli cell-free extract (CE). Furthermore, we probe PCA with three GTs with different specificities. Our results demonstrate that PCA is a versatile and apparently general GT assay with a detection limit as low as 1 mU. The detection limit of the pH shift assay is roughly 4 times higher. The immunoassay, by contrast, detected only nucleoside diphosphates (NDPs) but had the lowest detection limit. Compared with these assays, PCA showed superior robustness and, therefore, appears to be a suitable general screening assay for nucleotide sugar-dependent GTs.