Cloning and characterization of homeologous cellulose synthase catalytic subunit 2 genes from allotetraploid cotton (Gossypium hirsutum L.)

Rhamnolipid-Enriched PA3 Fraction from Pseudomonas aeruginosa SWUC02 Primes Chili Plant Defense Against Anthracnose

Chili anthracnose, caused by Colletotrichum truncatum, causes significant yield loss in chili production. In this study, we investigated the elicitor properties of a rhamnolipid (RL)-enriched PA3 fraction derived from Pseudomonas aeruginosa SWUC02 in inducing systemic resistance in yellow chili seedlings and antifungal activity against C. truncatum CFPL01 (Col). Fractionation of the ethyl acetate extract yielded 12 fractions, with PA3 demonstrating the most effective disease suppression, reducing the disease severity index to 4 ± 7.35% at 7 days post-inoculation compared with Col inoculation alone (83 ± 23.57%). PA3 also exhibited direct antifungal activity, inhibiting Col mycelial growth by 41 ± 0.96% at 200 µg/mL. Subfractionation revealed PA3 as a mixture of mono- and di-RLs, confirmed by 1H nuclear magnetic resonance and electrospray ionization mass spectrometry data. Additionally, PA3 enhanced seed germination and promoted plant growth without causing phytotoxicity. Transcriptomics revealed that PA3 pre-treatment prior to Col infection primed the defense response, upregulating defense-related genes involved in the phenylpropanoid, flavonoid, and jasmonic acid biosynthesis pathways, as well as those associated with cell wall reinforcement. Our findings highlight the potential of RL-enriched PA3 as both an antifungal agent and a plant defense elicitor, with transcriptome data providing new insights into defense priming and resistance pathways in chili, offering an eco-friendly solution for sustainable anthracnose management.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Changes in secondary metabolites and fiber quality of cotton (Gossypium hirsutum) seed under consecutive water stress and in silico analysis of cellulose synthase and xyloglucan endotransglucosylase

Global warming has led to severe drought conditions. The selection of plant varieties that can withstand drought and produce increased yields are of utmost importance. In the current study, secondary metabolites, seed trait and fiber characteristic of cottonseeds (Gossypium hirsutum) exposed to double and third water stress exposure was investigated. Total phenol and tannin content in W₁S₃₃ increased significantly after third water stress exposure. Accumulation of wax was enhanced in seeds of W₃S₃₃ and W₃S₃₄ that were subjected to third water stress. Fiber quality parameters decreased when cottonseeds were rainfed. High irrigation resulted in fragile and delicate fiber. Seeds grown under 66% FC irrigation saved water and produced seeds that had the potential of producing high quality fibers. In silico analysis was performed on cellulose synthase A (CesA) and xyloglucan endotransglycosylase (XET) enzymes present in Gossypium hirsutum. The intracellular locations of the CesA and XET1 enzymes are the plasma membrane and cell wall, respectively. Proline is conserved in the C-terminal of the CesA enzyme and plays an important role in enzyme functionality. This study provides a better understanding as to the mechanisms by which the plant can tolerate and combat water stress conditions as well as reduce water consumption. In order to grow cotton seeds with desirable morphometric characteristics and optimal fibers under water stress exposure and in dry areas, it is better to use seeds that are irrigated under optimal irrigation conditions, ie 66% FC.

Advances about the Roles of Membranes in Cotton Fiber Development

Cotton fiber is an extremely elongated single cell derived from the ovule epidermis and is an ideal model for studying cell development. The plasma membrane is tremendously expanded and accompanied by the coordination of various physiological and biochemical activities on the membrane, one of the three major systems of a eukaryotic cell. This review compiles the recent progress and advances for the roles of the membrane in cotton fiber development: the functions of membrane lipids, especially the fatty acids, sphingolipids, and phytosterols; membrane channels, including aquaporins, the ATP-binding cassette (ABC) transporters, vacuolar invertase, and plasmodesmata; and the regulation mechanism of membrane proteins, such as membrane binding enzymes, annexins, and receptor-like kinases.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

De Novo Transcriptome Assembly of Agave H11648 by Illumina Sequencing and Identification of Cellulose Synthase Genes in Agave Species

Agave plants are important crassulacean acid metabolism (CAM) plants with multiple agricultural uses, such as being used in tequila and fiber production. Agave hybrid H11648 ((A. amaniensis Trel. and Nowell × A. angustifolia Haw.) × A. amaniensis) is the main cultivated Agave species for fiber production in large tropical areas around the world. In this study, we conducted a transcriptome analysis of A. H11648. About 49.25 million clean reads were obtained by Illumina paired-end sequencing. De novo assembly produced 148,046 unigenes with more than 40% annotated in public databases, or matched homologs in model plants. More homologous gene pairs were found in Asparagus genome than in Arabidopsis or rice, which indicated a close evolutionary relationship between Asparagus and A. H11648. CAM-related gene families were also characterized as previously reported in A. americana. We further identified 12 cellulose synthase genes (CesA) in Asparagus genome and 38 CesA sequences from A. H11648, A. americana, A. deserti and A. tequilana. The full-length CesA genes were used as references for the cloning and assembly of their homologs in other Agave species. As a result, we obtained CesA1/3/4/5/7 genes with full-length coding region in the four Agave species. Phylogenetic and expression analysis revealed a conserved evolutionary pattern, which could not explain the distinct fiber traits in different Agave species. We inferred that transcriptional regulation might be responsible for Agave fiber development. This study represents the transcriptome of A. H11648, which would expand the number of Agave genes and benefit relevant studies of Agave fiber development.

The cotton XLIM protein (GhXLIM6) is required for fiber development via maintaining dynamic F-actin cytoskeleton and modulating cellulose biosynthesis

LIM domain proteins are cysteine-rich proteins, and are often considered as actin bundlers and transcription factors in plants. However, the roles of XLIM proteins in plants (especially in cotton) remain unexplored in detail so far. In this study, we identified a cotton XLIM protein (GhXLIM6) that is preferentially expressed in cotton fiber during whole elongation stage and early secondary cell wall (SCW) synthesis stage. The GhXLIM6-silenced transgenic cotton produces shorter fibers with thinner cell walls, compared with wild-type (WT). GhXLIM6 protein could directly bind F-actin and promote actin polymerization both in vitro and in vivo. It also acts as a transcription factor to suppress GhKNL1 expression through binding the PAL-box element of GhKNL1 promoter, and subsequently regulate the expression of CesA genes related to cellulose biosynthesis and deposition in SCWs of cotton fibers. The cellulose content in fibers of GhXLIM6RNAi cotton is lower than that in WT. Taken together, these data reveal the dual roles of GhXLIM6 in fiber development. On one hand, GhXLIM6 functions in fiber elongation through binding to F-actin to maintain the dynamic F-actin cytoskeleton. On the other hand, GhXLIM6 fine-tunes fiber SCW formation, probably through directly suppressing transcription of GhKNL1 to promote cellulose biosynthesis.

Draft genome analysis provides insights into the fiber yield, crude protein biosynthesis, and vegetative growth of domesticated ramie (Boehmeria nivea L. Gaud)

Plentiful bast fiber, a high crude protein content, and vigorous vegetative growth make ramie a popular fiber and forage crop. Here, we report the draft genome of ramie, along with a genomic comparison and evolutionary analysis. The draft genome contained a sequence of approximately 335.6 Mb with 42,463 predicted genes. A high-density genetic map with 4,338 single nucleotide polymorphisms (SNPs) was developed and used to anchor the genome sequence, thus, creating an integrated genetic and physical map containing a 58.2-Mb genome sequence and 4,304 molecular markers. A genomic comparison identified 1,075 unique gene families in ramie, containing 4,082 genes. Among these unique genes, five were cellulose synthase genes that were specifically expressed in stem bark, and 3 encoded a WAT1-related protein, suggesting that they are probably related to high bast fiber yield. An evolutionary analysis detected 106 positively selected genes, 22 of which were related to nitrogen metabolism, indicating that they are probably responsible for the crude protein content and vegetative growth of domesticated varieties. This study is the first to characterize the genome and develop a high-density genetic map of ramie and provides a basis for the genetic and molecular study of this crop.

Role of GhHDA5 in H3K9 deacetylation and fiber initiation in Gossypium hirsutum

Cotton fibers are single-celled trichomes that initiate from the epidermal cells of the ovules at or before anthesis. Here, we identified that the histone deacetylase (HDAC) activity is essential for proper cotton fiber initiation. We further identified 15 HDACs homoeologs in each of the A- and D-subgenomes of Gossypium hirsutum. Few of these HDAC homoeologs expressed preferentially during the early stages of fiber development [-1, 0 and 6 days post-anthesis (DPA)]. Among them, GhHDA5 expressed significantly at the time of fiber initiation (-1 and 0 DPA). The in vitro assay for HDAC activity indicated that GhHDA5 primarily deacetylates H3K9 acetylation marks. Moreover, the reduced expression of GhHDA5 also suppresses fiber initiation and lint yield in the RNA interference (RNAi) lines. The 0 DPA ovules of GhHDA5RNAi lines also showed alterations in reactive oxygen species homeostasis and elevated autophagic cell death in the developing fibers. The differentially expressed genes (DEGs) identified through RNA-seq of RNAi line (DEP12) and their pathway analysis showed that GhHDA5 modulates expression of many stress and development-related genes involved in fiber development. The reduced expression of GhHDA5 in the RNAi lines also resulted in H3K9 hyper-acetylation on the promoter region of few DEGs assessed by chromatin immunoprecipitation assay. The positively co-expressed genes with GhHDA5 showed cumulative higher expression during fiber initiation, and gene ontology annotation suggests their involvement in fiber development. Furthermore, the predicted protein interaction network in the positively co-expressed genes indicates HDA5 modulates fiber initiation-specific gene expression through a complex involving reported repressors.

De novo sequencing and assembly analysis of transcriptome in the Sodom apple (Calotropis gigantea)

Background

The Sodom apple (Calotropis gigantea), a member of the Asclepiadaceae family, is a large evergreen shrub native to continental Asia and northern Africa. As an important medicinal shrub and a fiber resource plant, there is an urgent need for developing molecular markers to facilitate breeding and genetic improvement of varieties.

Results

In this study, using the Illumina high throughput sequencing technique we obtained about 45 million paired end sequencing reads, De novo assembled and generated a total of 133,634 transcripts with a mean of 1837.47 bp in length. Based on protein homology searches against available databases, a total of 21,851 unigenes were functionally annotated. In particular, many transcripts that encode for putative proteins involved in fiber and secondary metabolite biosynthesis were identified and analyzed. Key fiber genes identified were validated experimentally through Real-Time PCR technique. Various transcription factors involved in regulating plant response to abiotic stress were also identified. In addition, based on the unigene sequences assembled, 11,623 microsatellites loci were detected, which provide very useful resources for developing microsatellite molecular markers.

Conclusion

This study is the first report on transcriptome information in the Calotropis species and provides rich gene transcript resources for conducting further studies on understanding the molecular basis of fiber and secondary metabolite biosynthesis, serving the genetic improvement and resource utilization in Calotropis plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1908-3) contains supplementary material, which is available to authorized users.

Two cotton fiber-associated glycosyltransferases, GhGT43A1 and GhGT43C1, function in hemicellulose glucuronoxylan biosynthesis during plant development

Xylan is the major hemicellulosic constituent in dicot secondary cell walls. Cell wall composition of cotton fiber changes dynamically throughout development. Not only the amounts but also the molecular sizes of the hemicellulosic polysaccharides show substantial changes during cotton fiber development. However, none of the genes encoding glycosyltransferases (GTs) responsible for synthesizing xylan have been isolated and characterized in cotton fiber. In this study, we applied a bioinformatics approach and identified two putative GTs from cotton, designated GhGT43A1 and GhGT43C1, which belong to the CAZy GT43 family and are closely related to Arabidopsis IRX9 and IRX14, respectively. We show that GhGT43A1 is highly and preferentially expressed in 15 and 20 days post-anthesis (dpa) cotton fiber, whereas GhGT43C1 is ubiquitously expressed in most organs, with especially high expression in 15 dpa fiber and hypocotyl. Complementation analysis demonstrates that GhG43A1 and GhGT43C1 are orthologs of Arabidopsis IRX9 and IRX14, respectively. Furthermore, we show that overexpression of GhGT43A1 or GhGT43C1 in Arabidopsis results in increased xylan content. We also show that overexpression of GhGT43A1 or GhGT43C1 leads to more cellulose deposition. These findings suggest that GhGT43A1 and GhGT43C1 likely participate in xylan synthesis during fiber development.

Isolation of developing secondary xylem specific cellulose synthase genes and their expression profiles during hormone signalling in Eucalyptus tereticornis

Cellulose synthases (CesA) represent a group of β-1, 4 glycosyl transferases involved in cellulose biosynthesis. Recent reports in higher plants have revealed that two groups of CesA gene families exist, which are associated with either primary or secondary cell wall deposition. The present study aimed at identifying developing secondary xylem specific cellulose synthase genes from Eucalyptus tereticornis, a species predominantly used in paper and pulp industries in the tropics. The differential expression analysis of the three EtCesA genes using qRT-PCR revealed 49 to 87 fold relative expression in developing secondary xylem tissues. Three full length gene sequences of EtCesA1, EtCesA2 and EtCesA3 were isolated with the size of 2940, 3114 and 3123 bp, respectively. Phytohormone regulation of all three EtCesA genes were studied by exogenous application of gibberellic acid, naphthalene acetic acid, indole acetic acid and 2, 4-epibrassinolide in internode tissues derived from three-month-old rooted cuttings. All three EtCesA transcripts were upregulated by indole acetic acid and gibberellic acid. This study demonstrates that the increased cellulose deposition in the secondary wood induced by hormones can be attributed to the upregulation of xylem specific CesAs.

Cotton KNL1, encoding a class II KNOX transcription factor, is involved in regulation of fibre development

In this study, the GhKNL1 (KNOTTED1-LIKE) gene, encoding a classical class II KNOX protein was identified in cotton (Gossypium hirsutum). GhKNL1 was preferentially expressed in developing fibres at the stage of secondary cell wall (SCW) biosynthesis. GhKNL1 was localized in the cell nucleus, and could interact with GhOFP4, as well as AtOFP1, AtOFP4, and AtMYB75. However, GhKNL1 lacked transcriptional activation activity. Dominant repression of GhKNL1 affected fibre development of cotton. The expression levels of genes related to fibre elongation and SCW biosynthesis were altered in transgenic fibres of cotton. As a result, transgenic cotton plants produced aberrant, shrunken, and collapsed fibre cells. Length and cell-wall thickness of fibres of transgenic cotton plants were significantly reduced compared with the wild type. Furthermore, overexpression and dominant repression of GhKNL1 in Arabidopsis resulted in a reduction in interfascicular fibre cell-wall thickening of basal stems of transgenic plants. Complementation revealed that GhKNL1 rescued the defective phenotype of Arabidopsis knat7 mutant in some extent. These data suggest that GhKNL1, as a transcription factor, participates in regulating fibre development of cotton.

Transcriptome profiling of Gossypium barbadense inoculated with Verticillium dahliae provides a resource for cotton improvement

BACKGROUND

Verticillium wilt, caused by the fungal pathogen Verticillium dahliae, is the most severe disease in cotton (Gossypium spp.), causing great lint losses worldwide. Disease management could be achieved in the field if genetically improved, resistant plants were used. However, the interaction between V. dahliae and cotton is a complicated process, and its molecular mechanism remains obscure. To understand better the defense response to this pathogen as a means for obtaining more tolerant cultivars, we monitored the transcriptome profiles of roots from resistant plants of G. barbadense cv. Pima90-53 that were challenged with V. dahliae.

RESULTS

In all, 46,192 high-quality expressed sequence tags (ESTs) were generated from a full-length cDNA library of G. barbadense. They were clustered and assembled into 23126 unigenes that comprised 2661 contigs and 20465 singletons. Those unigenes were assigned Gene Ontology terms and mapped to 289 KEGG pathways. A total of 3027 unigenes were found to be homologous to known defense-related genes in other plants. They were assigned to the functional classification of plant-pathogen interactions, including disease defenses and signal transduction. The branch of "SA→NPR1→TGA→PR-1→Disease resistance" was first discovered in the interaction of cotton-V. dahliae, indicating that this wilt process includes both biotrophic and necrotrophic stages. In all, 4936 genes coding for putative transcription factors (TF) were identified in our library. The most abundant TF family was the NAC group (527), followed by G2-like (440), MYB (372), BHLH (331), bZIP (271) ERF, C3H, and WRKY. We also analyzed the expression of genes involved in pathogen-associated molecular pattern (PAMP) recognition, the activation of effector-triggered immunity, TFs, and hormone biosynthesis, as well as genes that are pathogenesis-related, or have roles in signaling/regulatory functions and cell wall modification. Their differential expression patterns were compared among mock-/inoculated- and resistant/susceptible cotton. Our results suggest that the cotton defense response has significant transcriptional complexity and that large accumulations of defense-related transcripts may contribute to V. dahliae resistance in cotton. Therefore, these data provide a resource for cotton improvement through molecular breeding approaches.

CONCLUSIONS

This study generated a substantial amount of cotton transcript sequences that are related to defense responses against V. dahliae. These genomics resources and knowledge of important related genes contribute to our understanding of host-pathogen interactions and the defense mechanisms utilized by G. barbadense, a non-model plant system. These tools can be applied in establishing a modern breeding program that uses marker-assisted selections and oligonucleotide arrays to identify candidate genes that can be linked to valuable agronomic traits in cotton, including disease resistance.

An integrative analysis of four CESA isoforms specific for fiber cellulose production between Gossypium hirsutum and Gossypium barbadense

Cotton fiber is an excellent model system of cellulose biosynthesis; however, it has not been widely studied due to the lack of information about the cellulose synthase (CESA) family of genes in cotton. In this study, we initially identified six full-length CESA genes designated as GhCESA5-GhCESA10. Phylogenetic analysis and gene co-expression profiling revealed that CESA1, CESA2, CESA7, and CESA8 were the major isoforms for secondary cell wall biosynthesis, whereas CESA3, CESA5, CESA6, CESA9, and CESA10 should involve in primary cell wall formation for cotton fiber initiation and elongation. Using integrative analysis of gene expression patterns, CESA protein levels, and cellulose biosynthesis in vivo, we detected that CESA8 could play an enhancing role for rapid and massive cellulose accumulation in Gossypium hirsutum and Gossypium barbadense. We found that CESA2 displayed a major expression in non-fiber tissues and that CESA1, a housekeeping gene like, was predominantly expressed in all tissues. Further, a dynamic alteration was observed in cell wall composition and a significant discrepancy was observed between the cotton species during fiber elongation, suggesting that pectin accumulation and xyloglucan reduction might contribute to cell wall transition. In addition, we discussed that callose synthesis might be regulated in vivo for massive cellulose production during active secondary cell wall biosynthesis in cotton fibers.

De novo assembly and characterization of transcriptome using Illumina paired-end sequencing and identification of CesA gene in ramie (Boehmeria nivea L. Gaud)

BACKGROUND

Ramie fiber, extracted from vegetative organ stem bast, is one of the most important natural fibers. Understanding the molecular mechanisms of the vegetative growth of the ramie and the formation and development of bast fiber is essential for improving the yield and quality of the ramie fiber. However, only 418 expressed tag sequences (ESTs) of ramie deposited in public databases are far from sufficient to understand the molecular mechanisms. Thus, high-throughput transcriptome sequencing is essential to generate enormous ramie transcript sequences for the purpose of gene discovery, especially genes such as the cellulose synthase (CesA) gene.

RESULTS

Using Illumina paired-end sequencing, about 53 million sequencing reads were generated. De novo assembly yielded 43,990 unigenes with an average length of 824 bp. By sequence similarity searching for known proteins, a total of 34,192 (77.7%) genes were annotated for their function. Out of these annotated unigenes, 16,050 and 13,042 unigenes were assigned to gene ontology and clusters of orthologous group, respectively. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) indicated that 19,846 unigenes were mapped to 126 KEGG pathways, and 565 genes were assigned to http://starch and sucrose metabolic pathway which was related with cellulose biosynthesis. Additionally, 51 CesA genes involved in cellulose biosynthesis were identified. Analysis of tissue-specific expression pattern of the 51 CesA genes revealed that there were 36 genes with a relatively high expression levels in the stem bark, which suggests that they are most likely responsible for the biosynthesis of bast fiber.

CONCLUSION

To the best of our knowledge, this study is the first to characterize the ramie transcriptome and the substantial amount of transcripts obtained will accelerate the understanding of the ramie vegetative growth and development mechanism. Moreover, discovery of the 36 CesA genes with relatively high expression levels in the stem bark will present an opportunity to understand the ramie bast fiber formation and development mechanisms.

De novo assembly and characterization of transcriptome using Illumina paired-end sequencing and identification of CesA gene in ramie (Boehmeria nivea L. Gaud)

BACKGROUND

Ramie fiber, extracted from vegetative organ stem bast, is one of the most important natural fibers. Understanding the molecular mechanisms of the vegetative growth of the ramie and the formation and development of bast fiber is essential for improving the yield and quality of the ramie fiber. However, only 418 expressed tag sequences (ESTs) of ramie deposited in public databases are far from sufficient to understand the molecular mechanisms. Thus, high-throughput transcriptome sequencing is essential to generate enormous ramie transcript sequences for the purpose of gene discovery, especially genes such as the cellulose synthase (CesA) gene.

RESULTS

Using Illumina paired-end sequencing, about 53 million sequencing reads were generated. De novo assembly yielded 43,990 unigenes with an average length of 824 bp. By sequence similarity searching for known proteins, a total of 34,192 (77.7%) genes were annotated for their function. Out of these annotated unigenes, 16,050 and 13,042 unigenes were assigned to gene ontology and clusters of orthologous group, respectively. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) indicated that 19,846 unigenes were mapped to 126 KEGG pathways, and 565 genes were assigned to http://starch and sucrose metabolic pathway which was related with cellulose biosynthesis. Additionally, 51 CesA genes involved in cellulose biosynthesis were identified. Analysis of tissue-specific expression pattern of the 51 CesA genes revealed that there were 36 genes with a relatively high expression levels in the stem bark, which suggests that they are most likely responsible for the biosynthesis of bast fiber.

CONCLUSION

To the best of our knowledge, this study is the first to characterize the ramie transcriptome and the substantial amount of transcripts obtained will accelerate the understanding of the ramie vegetative growth and development mechanism. Moreover, discovery of the 36 CesA genes with relatively high expression levels in the stem bark will present an opportunity to understand the ramie bast fiber formation and development mechanisms.