Qingming Tang

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

BackgroundRamie 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.ResultsUsing 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.ConclusionTo 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.

Identification of drought stress-responsive transcription factors in ramie (Boehmeria nivea L. Gaud)

BackgroundRamie fiber extracted from stem bark is one of the most important natural fibers. Drought is a main environment stress which severely inhibits the stem growth of ramie and leads to a decrease of the fiber yield. The drought stress-regulatory mechanism of ramie is poorly understood.ResultUsing Illumina sequencing, approximately 4.8 and 4.7 million (M) 21-nt cDNA tags were respectively sequenced in the cDNA libraries derived from the drought-stressed ramie (DS) and the control ramie under well water condition (CO). The tags generated from the two libraries were aligned with ramie transcriptome to annotate their function and a total of 23,912 and 22,826 ramie genes were matched by these tags of DS and CO library, respectively. Comparison of gene expression level between CO and DS ramie based on the differences of tag frequencies appearing in the two libraries revealed that there were 1516 potential drought stress-responsive genes, in which 24 genes function as transcription factor (TF). Among these 24 TFs, the unigene19721 encoding the DELLA protein which is a key negative regulator in gibberellins (GAs) signal pathway was probably markedly up-regulated under water stress for a increase of tag abundance in DS library, which is possibly responsible for the inhibition of the growth of drought-stressed ramie. In order to validate the change of expression of these potential stress-responsive TFs under water deficit condition, the unigene19721 and another eleven potential stress-responsive TFs were chosen for further expression analysis in well-watered and drought-stressed ramie by real-time quantitative PCR (qRT-PCR) and the result showed that all 12 TFs were authentically involved in the response of drought stress.ConclusionIn this study, twelve TFs involving in the response of drought stress were first found by Illumina tag-sequencing and qRT-PCR in ramie. The discovery of these drought stress-responsive TFs will be helpful for further understanding the drought stress-regulatory mechanism of ramie and improving the drought tolerance ability of ramie.

Development and characterization of 1,827 expressed sequence tag-derived simple sequence repeat markers for ramie (Boehmeria nivea L. Gaud).

Ramie (Boehmeria nivea L. Gaud) is one of the most important natural fiber crops, and improvement of fiber yield and quality is the main goal in efforts to breed superior cultivars. However, efforts aimed at enhancing the understanding of ramie genetics and developing more effective breeding strategies have been hampered by the shortage of simple sequence repeat (SSR) markers. In our previous study, we had assembled de novo 43,990 expressed sequence tags (ESTs). In the present study, we searched these previously assembled ESTs for SSRs and identified 1,685 ESTs (3.83%) containing 1,878 SSRs. Next, we designed 1,827 primer pairs complementary to regions flanking these SSRs, and these regions were designated as SSR markers. Among these markers, dinucleotide and trinucleotide repeat motifs were the most abundant types (36.4% and 36.3%, respectively), whereas tetranucleotide, pentanucleotide, and hexanucleotide motifs represented <10% of the markers. The motif AG/CT was the most abundant, accounting for 28.74% of the markers. One hundred EST-SSR markers (97 SSRs located in genes encoding transcription factors and 3 SSRs in genes encoding cellulose synthases) were amplified using polymerase chain reaction for detecting 24 ramie varieties. Of these 100 markers, 98 markers were successfully amplified and 81 markers were polymorphic, with 2–6 alleles among the 24 varieties. Analysis of the genetic diversity of all 24 varieties revealed similarity coefficients that ranged from 0.51 to 0.80. The EST-SSRs developed in this study represent the first large-scale development of SSR markers for ramie. These SSR markers could be used for development of genetic and physical maps, quantitative trait loci mapping, genetic diversity studies, association mapping, and cultivar fingerprinting.

Genome-wide transcriptomic profiling of ramie (Boehmeria nivea L. Gaud) in response to cadmium stress.

Cadmium (Cd) contamination in agricultural soils has become a major environmental problem in China. Ramie, a fiber crop, has frequently been proposed for use as a phytoremediation crop for the restoration of Cd-contaminated farmlands. However, high levels of Cd can greatly inhibit stem growth in ramie, which reduces its economic value as a crop. To understand the potential mechanisms behind this phenomenon, the ramie genes involved in the Cd stress response were identified using Illumina pair-end sequencing on two Cd-stressed plants (CdS1 and CdS2) and two control plants (CO1 and CO2). Approximately 48.7, 51.6, 41.2, and 47.1 million clean sequence reads were generated from the libraries of CO1, CO2, CdS1, and CdS2, respectively, and de novo assembled to yield 56,932 non-redundant unigenes. A total of 26,686 (46.9%) genes were annotated for their function. Comparison of gene expression levels in CO and CdS ramie revealed 155 differentially expressed genes (DEGs) between treatment and control conditions. Sixteen DEGs were further analyzed for expression differences by using real-time quantitative PCR (qRT-PCR). Among these 16 DEGs, 2 genes encoding GA2-oxidase (a major enzyme for deactivating bioactive gibberellins [GAs]) showed markedly up-regulated expression in Cd stressed ramie. This might be responsible for the growth inhibition of Cd-stressed ramie. Pathway enrichment analysis revealed that the cutin, suberine and wax biosynthesis pathway was markedly enriched by DEGs. The discovery of these Cd stress-responsive genes and pathways will be helpful in further understanding the mechanism of Cd-stress response and improving Cd stress tolerance in ramie.

Genome-wide transcriptional changes of ramie (Boehmeria nivea L. Gaud) in response to root-lesion nematode infection.

Ramie fiber extracted from stem bark is one of the most important natural fibers. The root-lesion nematode (RLN) Pratylenchus coffeae is a major ramie pest and causes large fiber yield losses in China annually. The response mechanism of ramie to RLN infection is poorly understood. In this study, we identified genes that are potentially involved in the RLN-resistance in ramie using Illumina pair-end sequencing in two RLN-infected plants (Inf1 and Inf2) and two control plants (CO1 and CO2). Approximately 56.3, 51.7, 43.4, and 45.0 million sequencing reads were generated from the libraries of CO1, CO2, Inf1, and Inf2, respectively. De novo assembly for these 196 million reads yielded 50,486 unigenes with an average length of 853.3bp. A total of 24,820 (49.2%) genes were annotated for their function. Comparison of gene expression levels between CO and Inf ramie revealed 777 differentially expressed genes (DEGs). The expression levels of 12 DEGs were further confirmed by real-time quantitative PCR (qRT-PCR). Pathway enrichment analysis showed that three pathways (phenylalanine metabolism, carotenoid biosynthesis, and phenylpropanoid biosynthesis) were strongly influenced by RLN infection. A series of candidate genes and pathways that may contribute to the defense response against RLN in ramie will be helpful for further improving resistance to RLN infection.

Identification of a CONSTANS homologous gene with distinct diurnal expression patterns in varied photoperiods in ramie (Boehmeria nivea L. Gaud)

Ramie is an important natural fiber and forage crop in China. Breeding of late- or non-flowering varieties, with higher vegetative yields, is an important goal in ramie improvement. However, the ramie genes involved in flowering regulation have not previously been identified. In model plants, such as rice and Arabidopsis, the CONSTANS (CO) and CONSTANS-like (COL) genes play key roles in flowering regulation. In the present study, six ramie COL genes (BnCOL1-BnCOL6) with a full-length open reading frame (ORF) were identified. Sequence alignment revealed that all six BnCOL proteins contained conserved CCT (CO, COL, TOC1) and B-box I domains, but that only four of these proteins contained the B-box II domain. Expression pattern analysis showed that BnCOL1-BnCOL6 were mainly expressed in the stem xylem, flowers, and leaves. Phylogenetic analysis classified the six newly identified BnCOL proteins, and also 16 COL proteins with known functions in other species, into three groups. The BnCOL2-encoded protein was assigned to the same group as the CO- and Hd1-encoded proteins, suggesting that this BnCOL2-encoded protein showed the highest level of homology with the CO/Hd1-encoded proteins. Photoperiodic experiments showed that BnCOL2 exhibited a diurnal expression pattern under long- and short-day conditions. Subcellular localization examination revealed that the BnCOL2 protein fused with YFP was localized in the nucleus. Because the homologous sequence and similar expression pattern between BnCOL2 and CO/Hd1, the BnCOL2 possibly has a role in flowering modulation, and can be used as a candidate gene for research in the flowering regulation of ramie.