De novo characterization of the Chinese fir (Cunninghamia lanceolata) transcriptome and analysis of candidate genes involved in cellulose and lignin biosynthesis

Harnessing de novo transcriptome sequencing to identify and characterize genes regulating carbohydrate biosynthesis pathways in Salvia guaranitica L

Introduction

Carbohydrate compounds serve multifaceted roles, from energy sources to stress protectants, found across diverse organisms including bacteria, fungi, and plants. Despite this broad importance, the molecular genetic framework underlying carbohydrate biosynthesis pathways, such as starch, sucrose, and glycolysis/gluconeogenesis in Salvia guaranitica, remains largely unexplored.

Methods

In this study, the Illumina-HiSeq 2500 platform was used to sequence the transcripts of S. guaranitica leaves, generating approximately 8.2 Gb of raw data. After filtering and removing adapter sequences, 38 million reads comprising 210 million high-quality nucleotide bases were obtained. De novo assembly resulted in 75,100 unigenes, which were annotated to establish a comprehensive database for investigating starch, sucrose, and glycolysis biosynthesis. Functional analyses of glucose-6-phosphate isomerase (SgGPI), trehalose-6-phosphate synthase/phosphatase (SgT6PS), and sucrose synthase (SgSUS) were performed using transgenic Arabidopsis thaliana.

Results

Among the unigenes, 410 were identified as putatively involved in these metabolic pathways, including 175 related to glycolysis/gluconeogenesis and 235 to starch and sucrose biosynthesis. Overexpression of SgGPI, SgT6PS, and SgSUS in transgenic A. thaliana enhanced leaf area, accelerated flower formation, and promoted overall growth compared to wild-type plants.

Discussion

These findings lay a foundation for understanding the roles of starch, sucrose, and glycolysis biosynthesis genes in S. guaranitica, offering insights into future metabolic engineering strategies for enhancing the production of valuable carbohydrate compounds in S. guaranitica or other plants.

An Optimized and Cost-Effective RNA Extraction Method for Secondary Metabolite-Enriched Tissues of Norway Spruce (Picea abies)

Since the development of next-generation sequencing techniques and with the growing interest in transcriptomic studies, there is a demand for high-throughput RNA extraction techniques. General RNA extraction protocols are unreliable when it comes to the quality and quantity of isolated RNA obtained from different tissue types of different plant species. Despite Norway spruce (Picea abies) being one of the most significant and commercially valuable tree species in European forests, only limited genetic research is available. In this study, we developed a cetyltrimethylammonium bromide (CTAB) protocol by modifying the original method. We compared this CTAB protocol with other widely used methods for extracting RNA from different tissues (needle, phloem, and root) of Norway spruce, known for its richness in polyphenols, polysaccharides, and secondary metabolites. The modified CTAB method proves to be superior to the kit-based and TRIzol-based methods for extracting RNA from the metabolite-rich tissues of Norway spruce, resulting in high RNA quality and integrity values (RIN~7-9). The modified CTAB RNA extraction method is rapid, cost-effective, and relatively simple in yielding the desired RNA quality from Norway spruce tissues. It is optimal for RNA sequencing and other downstream molecular applications.

Pyrolysis Characteristics and Determination of Kinetic and Thermodynamic Parameters of Raw and Torrefied Chinese Fir

Torrefaction influences the structural and physicochemical properties of biomass, thus further altering its thermal degradation behavior. In this study, the pyrolysis characteristics, reaction kinetics, and thermodynamic parameters of raw and torrefied Chinese fir (CF) were investigated. The torrefaction was conducted at 220 °C (mild) and 280 °C (severe), the pyrolysis was performed from ambient temperature to 600 °C, and four different heating rates (i.e., 5, 15, 25, and 35 °C/min) were adopted. The activation energy for pyrolysis was estimated by adopting three isoconversional methods. The master-plot method was employed to analyze the reaction mechanism. Furthermore, thermodynamic parameters, i.e., the enthalpy change (ΔH), Gibbs free energy change (ΔG), and entropy change (ΔS), were calculated. The average activation energy increased with the torrefaction temperature, whose values estimated by using different methods ranged from 88.57 to 97.70, from 121.04 to 126.35, and from 167.51 to 179.74 kJ/mol for raw, mildly, and severely torrefied CF samples, respectively. A compensation effect between the activation energy and pre-exponential factor was observed for all samples. The degradation process was characterized as endothermic, involving the formation of activated complexes and requiring extra energy for torrefied samples.

DNA-free high-quality RNA extraction from 39 difficult-to-extract plant species (representing seasonal tissues and tissue types) of 32 families, and its validation for downstream molecular applications

BACKGROUND

High-purity RNA serves as the basic requirement for downstream molecular analysis of plant species, especially the differential expression of genes to various biotic and abiotic stimuli. But, the extraction of high-quality RNA is usually difficult from plants rich in polysaccharides and polyphenols, and their presence usually interferes with the downstream applications. The aim of the study is to optimize the extraction of high-quality RNA from diverse plant species/tissues useful for downstream molecular applications.

RESULTS

Extraction of RNA using commercially available RNA extraction kits and routine hexadecyltrimethylammonium bromide (CTAB) methods did not yield good quality DNA-free RNA from Prosopis cineraria, Conocarpus erectus, and Phoenix dactylifera. A reliable protocol for the extraction of high-quality RNA from mature leaves of these difficult-to-extract trees was optimized after screening nine different methods. The DNase I-, and proteinase K treatment-free modified method, consisting of extraction with CTAB method followed by TRIzol, yielded high-quality DNA-free RNA with an A260/A280 and A260/A230 ratios > 2.0. Extraction of RNA from Conocarpus, the most difficult one, was successful by avoiding the heat incubation of ground tissue in a buffer at 65 oC. Pre-warming of the buffer for 5-10 min was sufficient to extract good-quality RNA. RNA integrity number of the extracted RNA samples ranged between 7 and 9.1, and the gel electrophoresis displayed intact bands of 28S and 18S RNA. A cDNA library constructed from the RNA of P. cineraria was used for the downstream applications. Real-time qPCR analysis using the cDNA from P. cineraria RNA confirmed the quality. The extraction of good quality RNA from samples of the desert-growing P. cineraria (> 20-years-old) collected in alternate months of the year 2021 (January to December covering winter, spring, autumn, and the very dry and hot summer) proved the efficacy of the protocol. The protocol's broad applicability was further validated by extracting good-quality RNA from 36 difficult-to-extract plant species, including tissues such as roots, flowers, floral organs, fruits, and seeds.

CONCLUSIONS

The modified DNase I and Proteinase K treatment-free protocol enables to extract DNA-free, high-quality, intact RNA from a total of 39 difficult-to-extract plant species belonging to 32 angiosperm families is useful to extract good-quality RNA from dicots and monocots irrespective of tissue types and growing seasons.

Transcriptome analysis identifies differentially expressed genes involved in lignin biosynthesis in barley

Lignin is an essential metabolite for plant growth but negatively affects the quality of forage barley. Genetic modification of quality traits to improve the forage digestibility requires an understanding of the molecular mechanism of lignin biosynthesis. RNA-Seq was used to quantify transcripts differentially expressed among leaf, stem and spike tissues from two barley genotypes. A total of 13,172 differentially expressed genes (DEGs) were identified, of which much more up-regulated DEGs were detected from the contrasting groups of leaf vs spike (L-S) and stem vs spike (S-S), and down-regulated DEGs were dominant in the group of stem vs leaf (S-L). 47 DEGs were successfully annotated to the monolignol pathway and six of them were candidate genes regulating the lignin biosynthesis. The qRT-PCR assay verified the expression profiles of the six candidate genes. Among them, four genes might positively regulate the lignin biosynthesis during forage barley development in terms of the consistency of their expression levels and changes of lignin content among the tissues, while the other two genes may have the reverse effects. These findings provide target genes for further investigations on molecular regulatory mechanisms of lignin biosynthesis and genetic resources for improvement of forage quality in barley molecular breeding programme.

Comparative transcriptome analysis reveals an insight into the candidate genes involved in anthocyanin and scent volatiles biosynthesis in colour changing flowers of Combretum indicum

Combretum indicum is a widely cultivated ornamental species displaying the distinct phenomenon of floral colour change. Flowers display a gradual colour change from white to red, attributed to increased cyanidin 3-O glucoside in petal tissues. The differently coloured flowers also emanate a complex blend of VOCs with trans-linalool oxide (furanoid) as the major compound in the emission profile. To understand molecular mechanisms regulating floral colour shifts and scent biosynthesis, we performed Illumina transcriptome sequencing, including de novo assembly and functional annotation, for the two stages of floral maturation (white and red). Homology analysis with functional classification identified 84 and 42 candidate genes associated with pigment and scent biosynthesis, respectively. Genes encoding transcription factors, such as MYB, ERF, WD40, WRKY, NAC, bHLH and bZIP, that play critical roles in regulating specialized metabolism were also identified in the transcriptome data. Differences in expression of genes were consistent with accumulation patterns of anthocyanins in the two different flower colours. A clear upregulation of flavonoid biosynthesis genes in red flower tissue is associated with increased pigment content. RT-qPCR-based expression analyses gave results consistent with the RNA-Seq data, suggesting the sequencing data are consistent and reliable. This study presents the first report of genetic information for C. indicum. Gene sequences generated from RNA-Seq, along with candidate genes identified by pathway mapping and their expression profiles, provide a valuable resource for subsequent studies towards molecular understanding of specialized metabolism in C. indicum flowers.

Comparative full-length transcriptome analysis by Oxford Nanopore Technologies reveals genes involved in anthocyanin accumulation in storage roots of sweet potatoes (Ipomoea batatas L.)

Background

Storage roots of sweet potatoes (Ipomoea batatas L.) with different colors vary in anthocyanin content, indicating different economically agronomic trait. As the newest DNA/RNA sequencing technology, Oxford Nanopore Technologies (ONT) have been applied in rapid transcriptome sequencing for investigation of genes related to nutrient metabolism. At present, few reports concern full-length transcriptome analysis based on ONT for study on the molecular mechanism of anthocyanin accumulation leading to color change of tuberous roots of sweet potato cultivars.

Results

The storage roots of purple-fleshed sweet potato (PFSP) and white-fleshed sweet potato (WFSP) at different developmental stages were subjected to anthocyanin content comparison by UV-visible spectroscopy as well as transcriptome analysis at ONT MinION platform. UV-visible spectrophotometric measurements demonstrated the anthocyanin content of PFSP was much higher than that of WFSP. ONT RNA-Seq results showed each sample generated average 2.75 GB clean data with Full-Length Percentage (FL%) over 70% and the length of N50 ranged from 1,192 to 1,395 bp, indicating reliable data for transcriptome analysis. Subsequent analysis illustrated intron retention was the most prominent splicing event present in the resulting transcripts. As compared PFSP with WFSP at the relative developmental stages with the highest (PH vs. WH) and the lowest (PL vs. WL) anthocyanin content, 282 and 216 genes were up-regulated and two and 11 genes were down-regulated respectively. The differential expression genes involved in flavonoid biosynthesis pathway include CCoAOMT, PpLDOX, DFR, Cytochrome P450, CHI, and CHS. The genes encoding oxygenase superfamily were significantly up-regulated when compared PFSP with WFSP at the relative developmental stages.

Conclusions

Comparative full-length transcriptome analysis based on ONT serves as an effective approach to detect the differences in anthocyanin accumulation in the storage roots of different sweet potato cultivars at transcript level, with noting that some key genes can now be closely related to flavonoids biosynthesis. This study helps to improve understanding of molecular mechanism for anthocyanin accumulation in sweet potatoes and also provides a theoretical basis for high-quality sweet potato breeding.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analyses Provide Comprehensive Insight Into Molecular Mechanisms of Cellulose and Lignin Biosynthesis in Cunninghamia lanceolata

Cunninghamia lanceolata is an essential timber species that provide 20%-30% raw materials for China's timber industry. Although a few transcriptomes have been published in C. lanceolata, full-length mRNA transcripts and regulatory mechanisms behind the cellulose and lignin biosynthesis have not been thoroughly investigated. Here, PacBio Iso-seq and RNA-seq analyses were adapted to identify the full-length and differentially expressed transcripts along a developmental gradient from apex to base of C. lanceolata shoots. A total of 48,846 high-quality full-length transcripts were obtained, of which 88.0% are completed transcriptome based on benchmarking universal single-copy orthologs (BUSCO) assessment. Along stem developmental gradient, 18,714 differentially expressed genes (DEGs) were detected. Further, 28 and 125 DEGs were identified as enzyme-coding genes of cellulose and lignin biosynthesis, respectively. Moreover, 57 transcription factors (TFs), including MYB and NAC, were identified to be involved in the regulatory network of cellulose and lignin biosynthesis through weighted gene co-expression network analysis (WGCNA). These TFs are composed of a comparable regulatory network of secondary cell wall formation in angiosperms, revealing a similar mechanism may exist in gymnosperms. Further, through qRT-PCR, we also investigated eight specific TFs involved in compression wood formation. Our findings provide a comprehensive and valuable source for molecular genetics breeding of C. lanceolata and will be beneficial for molecular-assisted selection.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

BACKGROUND

R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line.

RESULTS

The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms.

CONCLUSIONS

This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Genome-wide transcriptional analysis unveils the molecular basis of organ-specific expression of isosteroidal alkaloids biosynthesis in critically endangered Fritillaria roylei Hook

Fritillaria roylei Hook. is a critically endangered high altitude Himalayan medicinal plant species with rich source of pharmaceutically active structurally diverse steroidal alkaloids. Nevertheless, except few marker compounds, the chemistry of the plant remains unexplored. Therefore, in the current study, transcriptome sequencing efforts were made to elucidate isosteroidal alkaloids biosynthesis by creating first organ-specific genomic resource using bulb, stem, and leaf tissues derived from natural populations of Indian Himalayan region. Overall, 349.9 million high quality paired-end reads obtained using NovaSeq 6000 platform were assembled (de novo) into 82,848 unigenes and 31,061 isoforms. Functional annotation and organ specific differential expression (DE) analysis identified 2488 significant DE transcripts, grouped into three potential sub-clusters (sub-cluster I: 728 transcripts; sub-cluster II: 446 transcripts and Sub-cluster III: 1314 transcripts). Subsequently, pathway enrichment (GO, KEGG) and protein-protein network analysis revealed significantly higher enrichment of phenyl-propanoid and steroid backbone including terpenoid, sesquiterpenoid and triterpenoid biosynthesis in bulb. Additionally, upregulated expression of cytochrome P450, UDP-dependent Glucuronosyltransferase families and key transcription factor families (FAR1, bHLH, GRAS, C2H2, TCP and MYB) suggests 'bulb' as a primary site of MVA mediated isosteroidal alkaloids biosynthesis. The comprehensive elucidation of molecular insights in this study is a first step towards the understanding of isosteroidal alkaloid biosynthesis pathway in F. roylei. Furthermore, key genes and regulators identified here can facilitate metabolic engineering of potential bioactive compounds at industrial scale.

The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development

BACKGROUND

Cunninghamia lanceolata (Chinese fir), a member of the conifer family Cupressaceae, is one of the most popular cultivated trees for wood production in China. Continuous research is being performed to improve C. lanceolata breeding values. Given the high rate of seed abortion (one of the reasons being the failure of ovule and pollen development) in C. lanceolata, the proper formation of female/male cones could theoretically increase the number of offspring in future generations. MIKC MADS-box genes are well-known for their roles in the flower/cone development and comprise the typical/atypical floral development model for both angiosperms and gymnosperms.

RESULTS

We performed a transcriptomic analysis to find genes differentially expressed between female and male cones at a single, carefully determined developmental stage, focusing on the MIKC MADS-box genes. We finally obtained 47 unique MIKC MADS-box genes from C. lanceolata and divided these genes into separate branches. 27 out of the 47 MIKC MADS-box genes showed differential expression between female and male cones, and most of them were not expressed in leaves. Out of these 27 genes, most B-class genes (AP3/PI) were up-regulated in the male cone, while TM8 genes were up-regulated in the female cone. Then, with no obvious overall preference for AG (class C + D) genes in female/male cones, it seems likely that these genes are involved in the development of both cones. Finally, a small number of genes such as GGM7, SVP, AGL15, that were specifically expressed in female/male cones, making them candidate genes for sex-specific cone development.

CONCLUSIONS

Our study identified a number of MIKC MADS-box genes showing differential expression between female and male cones in C. lanceolata, illustrating a potential link of these genes with C. lanceolata cone development. On the basis of this, we postulated a possible cone development model for C. lanceolata. The gene expression library showing differential expression between female and male cones shown here, can be used to discover unknown regulatory networks related to sex-specific cone development in the future.

De novo transcriptome analysis of Rhizophora mucronata Lam. furnishes evidence for the existence of glyoxalase system correlated to glutathione metabolic enzymes and glutathione regulated transporter in salt tolerant mangroves

The accumulation of a metabolic by product - methylglyoxal above a minimal range can be highly toxic in all organisms. Stress induced elevation in methylglyoxal inactivates proteins and nucleic acids. Glutathione dependent glyoxalase enzymes like glyoxalase I and glyoxalase II together with glutathione independent glyoxalase III play inevitable role in methylglyoxal detoxification. Glyoxalase genes are generally conserved but with obvious exceptions. Mangroves being potent harsh land inhabitants, their internal organelles are constantly been exposed to elevated levels of methylglyoxal. First and foremost it is important to detect the presence of glyoxalases in mangroves. De novo transcriptome analysis of mangrove species Rhizophora mucronata Lam., identified eleven putative glyoxalase proteins (RmGLYI-1 to 5, RmGLYII-1 to 5 and RmGLYIII). Molecular characterization proposed PLN02300 or PLN02367 as the key domains of RmGLYI proteins. They possess molecular weight ranging from 26.45 to 32.53 kDa and may localize in cytosol or chloroplast. RmGLYII proteins of molecular weight 28.64-36 kDa, carrying PLN02398 or PLN02469 domains are expected to be localized in diverse cellular compartments. Cytosolic RmGLYIII with DJ-1/PfpI domain carries a molecular weight 26.4 kDa. Detailed structural analysis revealed monomeric nature of RmGLYI-1 and RmGLYII-1 whereas RmGLYIII is found to be homodimer. Molecular phylogenetic analysis and multiple sequence alignment specified conserved metal ion/substrate binding residues of RmGLY proteins. Estimation of relative expression of glyoxalases under salt stress indicated the prominence of RmGLYI and RmGLYII over RmGLYIII. The aforementioned prominence is supported by salt induced expression difference of glutathione metabolic enzymes and glutathione regulated transporter protein.

Transcriptomic Profiling of Cryptomeria fortunei Hooibrenk Vascular Cambium Identifies Candidate Genes Involved in Phenylpropanoid Metabolism

Cryptomeria fortunei Hooibrenk (Chinese cedar) is a coniferous tree from southern China that has an important function in landscaping and timber production. Lignin is one of the key components of secondary cell walls, which have a crucial role in conducting water and providing mechanical support for the upward growth of plants. It is mainly biosynthesized via the phenylpropanoid metabolic pathway, of which the molecular mechanism remains so far unresolved in C. fortunei. In order to obtain further insight into this pathway, we performed transcriptome sequencing of the C. fortunei cambial zone at 5 successive growth stages. We generated 78,673 unigenes from transcriptome data, of which 45,214 (57.47%) were successfully annotated in the non-redundant protein database (NR). A total of 8975 unigenes were identified to be significantly differentially expressed between Sample_B and Sample_A after analyzing their expression profiles. Of the differentially expressed genes (DEGs), 6817 (75.96%) and 2158 (24.04%) were up- and down-regulated, respectively. 83 DEGs were involved in phenylpropanoid metabolism, 37 DEGs that encoded v-Myb avian myeloblastosis viral oncogene homolog (MYB) transcription factor (TF), and many candidates that encoded lignin synthesizing enzymes. These findings contribute to understanding the expression pattern of C. fortunei cambial zone transcriptome. Furthermore, our results provide additional insight towards understanding the molecular mechanisms of wood formation in C. fortunei.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

RNA Sequencing-Associated Study Identifies GmDRR1 as Positively Regulating the Establishment of Symbiosis in Soybean

In soybean (Glycine max)-rhizobium interactions, the type III secretion system (T3SS) of rhizobium plays a key role in regulating host specificity. However, the lack of information on the role of T3SS in signaling networks limits our understanding of symbiosis. Here, we conducted an RNA sequencing analysis of three soybean chromosome segment substituted lines, one female parent and two derived lines with different chromosome-substituted segments of wild soybean and opposite nodulation patterns. By analyzing chromosome-linked differentially expressed genes in the substituted segments and quantitative trait loci (QTL)-assisted selection in the substituted-segment region, genes that may respond to type III effectors to mediate plant immunity-related signaling were identified. To narrow down the number of candidate genes, QTL assistant was used to identify the candidate region consistent with the substituted segments. Furthermore, one candidate gene, GmDRR1, was identified in the substituted segment. To investigate the role of GmDRR1 in symbiosis establishment, GmDRR1-overexpression and RNA interference soybean lines were constructed. The nodule number increased in the former compared with wild-type soybean. Additionally, the T3SS-regulated effectors appeared to interact with the GmDDR1 signaling pathway. This finding will allow the detection of T3SS-regulated effectors involved in legume-rhizobium interactions.

Genome survey of Chinese fir (Cunninghamia lanceolata): Identification of genomic SSRs and demonstration of their utility in genetic diversity analysis

Chinese fir (Cunninghamia lanceolata) is an important coniferous species that accounts for 20–30% of the total commercial timber production in China. Though traditional breeding of Chinese fir has achieved remarkable success, molecular-assisted breeding has made little progress due to limited availability of genomic information. In this study, a survey of Chinese fir genome was performed using the Illumina HiSeq Xten sequencing platform. K-mer analysis indicated that Chinese fir has a large genome of approximately 11.6 Gb with 74.89% repetitive elements and is highly heterozygous. Meanwhile, its genome size was estimated to be 13.2 Gb using flow cytometry. A total of 778.02 Gb clean reads were assembled into 10,982,272 scaffolds with an N50 of 1.57 kb. In total, 362,193 SSR loci were detected with a frequency of 13.18 kb. Dinucleotide repeats were the most abundant (up to 73.6% of the total SSRs), followed by trinucleotide and tetranucleotide repeats. Forty-six polymorphic pairs were developed, and 298 alleles were successfully amplified from 199 Chinese fir clones. The average PIC value was 0.53, indicating that the identified genomic SSR (gSSR) markers have a high degree of polymorphism. In addition, these breeding resources were divided into three groups, and a limited gene flow existed among these inferred groups.

Physiological, biochemical and gene-expressional responses to water deficit in apple subjected to partial root-zone drying (PRD)

Water scarcity is one of the major factors limiting apple production. Partial root-zone drying (PRD) is a water-saving irrigation technique necessary to improve the efficiency of irrigation techniques to optimize the amount of fruit produced with the volume of water used. The apple trees cv. Red Delicious were exposed to four treatments, including (1) control with 100% of the crop evapotranspiration (ETc) needs; (2) alternate partial root-zone drying with 75% of the ETc needs (APRD75); (3) fixed partial root-zone drying with 75% of the ETc needs (FPRD75); (4) fixed partial root-zone irrigation with 50% of the ETc needs (FPRD50) in a semiarid region of Iran. Results showed that leaf water potential (Ψ _leaf), and chlorophyll were significantly decreased in FPRD50 compared to control and other PRD treatments. APRD75 and FPRD75 treatments significantly enhanced (+) -catechin (+C), epicatechin (EC), chlorogenic acid (CGA), caffeic acid (CA) as well as increased water use efficiency (WUE) (by 30-40% compared to control) without significant reduction of yield. PRD reduced gibberellic acid (GA3) and kinetin, while, increased the abscisic acid (ABA) and salicylic acid (SA) levels. The abiotic stress-responsive transcription factors (TFs) MdoMYB121, MdoMYB155, MdbZIP2, and MdbZIP48 were highly expressed in all PRD treatments. Our results demonstrated that APRD75 and FPRD75 have the potential to stimulate antioxidant defense mechanisms, hormonal signaling pathways, and expression of drought-tolerance TFs to improve WUE while maintaining crop yield. Therefore, APRD75andFPRD75 with water savings as compared to full irrigation might be a suitable strategy for irrigation apple trees under water scarcity.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

De novo transcriptome sequencing of radish (Raphanus sativus L.) fleshy roots: analysis of major genes involved in the anthocyanin synthesis pathway

Background

The HongXin radish (Raphanus sativus L.), which contains the natural red pigment (red radish pigment), is grown in the Fuling district of Chongqing City. However, the molecular mechanisms underlying anthocyanin synthesis for the formation of natural red pigment in the fleshy roots of HongXin radish are not well studied.

Results

De novo transcriptome of HX-1 radish, as well as that of the advanced inbred lines HX-2 and HX-3 were characterized using next generation sequencing (NGS) technology. In total, approximately 66.22 million paired-end reads comprising 34, 927 unigenes (N50 = 1, 621 bp) were obtained. Based on sequence similarity search with known proteins, total of 30, 127 (about 86.26%) unigenes were identified. Additionally, functional annotation and classification of these unigenes indicated that most of the unigenes were predominantly enriched in the metabolic process-related terms, especially for the biosynthetic pathways of secondary metabolites. Moreover, majority of the anthocyanin biosynthesis-related genes (ABRGs) involved in the regulation of anthocyanin biosynthesis were identified by targeted search for their annotation. Subsequently, the expression of 15 putative ABRGs involved in the anthocyanin synthesis-related pathways were validated using quantitative real-time polymerase chain reaction (qRT-PCR). Of those, RsPAL2, RsCHS-B2, RsDFR1, RsDFR2, RsFLS, RsMT3 and RsUFGT73B2-like were identified significantly associated with anthocyanin biosynthesis. Especially for RsDFR1, RsDFR2 and RsFLS, of those, RsDFR1 and RsDFR2 were highest enriched in the HX-3 and WG-3, but RsFLS were down-regulated in HX-3 and WG-3. We proposed that the transcripts of RsDFR1, RsDFR2 and RsFLS might be act as key regulators in anthocyanin biosynthesis pathway.

Conclusions

The assembled radish transcript sequences were analysed to identify the key ABRGs involved in the regulation of anthocyanin biosynthesis. Additionally, the expression patterns of candidate ABRGs involved in the anthocyanin biosynthetic pathway were validated by qRT-PCR. We proposed that the transcripts of RsDFR1, RsDFR2 and RsFLS might be acted as key regulators in anthocyanin biosynthesis pathway. This study will enhance our understanding of the biosynthesis and metabolism of anthocyanin in radish.

Effects of Sucrose, Boric Acid, pH, and Incubation Time on in Vitro Germination of Pollen and Tube Growth of Chinese fir (Cunnighamial lanceolata L.)

In vitro pollen germination provides a novel approach and strategy to accelerate genetic improvement of tree breeding. Studies about pollen germination and tube growth of Chinese fir are limited. Therefore, this study aimed to investigate the effects of sucrose, boric acid, pH, and time of incubation on pollen germination and tube growth. Pollen from 9 clones were selected. In vitro germination was performed in basic media as control, and in different concentrations of sucrose (0, 10 and 15%), boric acid (0.01, 0.1 and 0.2%), and pH levels (4.5, 5 and 7). Pollen germination rates and tube growth were recorded periodically at 1, 12, 24, and 48 h. The results showed that sucrose imposes significant effects on pollen germination and tube growth. The effects are most obvious at concentration of 15%. Boric acid significantly promoted germination and tube growth. The promotion was most notable in lower concentration of 0.01%. The media adjusted to pH 7.0 boosted the germination and pollen tube growth. The optimum time of incubation was 24 and 48 h for pollen germination and tube growth, respectively. Sucrose, pH, and time of incubation were positively correlated, whereas boric acid negatively correlated with pollen germination and tube growth. This study provided experimental evidences for selecting viable pollens for Chinese fir breeding.

Molecular Mechanism of Xylogenesis in Moso Bamboo (Phyllostachys edulis) Shoots during Cold Storage

A bamboo shoot is the immature stem of the woody grass and a nutritious and popular vegetable in East Asia. However, it undergoes a rapid xylogenesis process right after harvest, even being stored in a cold chamber. To investigate the molecular regulation mechanisms of xylogenesis in Moso bamboo (Phyllostachys edulis) shoots (MBSes) during cold storage, the measurement of cell wall polymers (cellulose, hemicellulose, and lignin) and related enzyme activities (phenylalanine ammonia lyase (PAL), cinnamyl alcohol dehydrogenase (CAD), peroxidase (POD), and xylan xylosyltransferase (XylT)) and transcriptomic analysis were performed during cold storage. It was noticed that cellulose and lignin contents increased, while hemicellulose content exhibited a downward trend. PAL, CAD, and POD activity presented an upward trend generally in MBS when stored at 4 °C for 16 days. XylT activity showed a descending trend during the stages of storage, but slightly increased during the 8th to 12th days after harvest at 4 °C. Transcriptomic analysis identified 72, 28, 44, and 31 functional unigenes encoding lignin, cellulose, xylan biosynthesis enzymes, and transcription factors (TFs), respectively. Many of these secondary cell wall (SCW)-related genes showed higher expression levels in the later period of cold storage. Quantitative RT-PCR analysis of the selected genes conformed to the expression pattern. Our study provides a comprehensive analysis of MBS secondary wall biosynthesis at the molecular level during the cold storage process. The results give insight into the xylogenesis process of this economically important vegetable and shed light on solving this problem of the post-harvest industry.

De novo transcriptome sequencing and metabolite profiling analyses reveal the complex metabolic genes involved in the terpenoid biosynthesis in Blue Anise Sage (Salvia guaranitica L.)

Many terpenoid compounds have been extracted from different tissues of Salvia guaranitica. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 4 Gb of raw data were generated from the transcriptome of S. guaranitica leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 32 million, comprising 186 million of high-quality nucleotide bases. A total of 61,400 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 267 unigenes that are putatively involved in terpenoid metabolism (including, 198 mevalonate and methyl-erythritol phosphate (MEP) pathways, terpenoid backbone biosynthesis genes and 69 terpene synthases genes). Moreover, three terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic Arabidopsis; most transgenic Arabidopsis plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich blue anise sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. guaranitica.

Identification and expression profiling of genes governing lignin biosynthesis in Casuarina equisetifolia L

Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.

RNA-Seq analysis of differential gene expression in Betula luminifera xylem during the early stages of tension wood formation

Background

Betula luminifera H. Winkler, which is widely distributed in southern China, is an economically important broadleaf tree species. However, little genomic information of B. luminifera is available, and little is known about the molecular mechanisms of wood formation in this species. Meanwhile, few efforts have focused on investigating the early transcriptional changes during tension wood formation in woody plants.

Results

A reference transcriptome dataset was first generated containing 45,700 Unigenes, and 35,135 (76.9%) Unigenes were annotated by a BLAST similarity search against four public databases. Then, based on an anatomical investigation, the global gene expression changes during the early stages of tension wood formation were analyzed. Gene expression profiling showed that a total of 13,273 Unigenes were differentially regulated during the early stages of tension wood formation. Most genes involved in cellulose and lignin biosynthesis were highlighted to reveal their biological importance in tension wood formation. In addition, the transcription levels of many genes involved in the auxin response pathway were significantly changed during the early stages of tension wood formation. Furthermore, 18 TFs co-expressed with key enzymes of cellulose synthesis were identified.

Conclusions

Our results revealed the transcriptional changes associated with TW formation and identified potential key genes in the regulation of this process. These results will help to dissect the molecular mechanism of wood formation and provide key candidate genes for marker-assisted selection in B. luminifera.

Transcriptome analysis provides insights into xylogenesis formation in Moso bamboo (Phyllostachys edulis) shoot

Maturation-related changes in cell wall composition and the molecular mechanisms underlying cell wall changes were investigated from the apical, middle and basal segments in moso bamboo shoot (MBS). With maturation extent from apical to basal regions in MBS, lignin and cellulose content increased, whereas heteroxylan exhibited a decreasing trend. Activities of phenylalanine amonnialyase (PAL), cinnamyl alcohol dehydrogenase (CAD) and cinnamate-4-hydroxylase (C4H), which are involved in lignin biosynthesis, increased rapidly from the apex to the base sections. The comparative transcriptomic analysis was carried out to identify some key genes involved in secondary cell walls (SCW) formation underlying the cell wall compositions changes including 63, 8, 18, and 31 functional unigenes encoding biosynthesis of lignin, cellulose, xylan and NAC-MYB-based transcription factors, respectively. Genes related to secondary cell wall formation and lignin biosynthesis had higher expression levels in the middle and basal segments compared to those in the apical segments. Furthermore, the expression profile of PePAL gene showed positive relationships with cellulose-related gene PeCESA4, xylan-related genes PeIRX9 and PeIRX10. Our results indicated that lignification occurred in the more mature middle and basal segments in MBS at harvest while lignification of MBS were correlated with higher expression levels of PeCESA4, PeIRX9 and PeIRX10 genes.

Desiccation Treatment and Endogenous IAA Levels Are Key Factors Influencing High Frequency Somatic Embryogenesis in Cunninghamia lanceolata (Lamb.) Hook

Cunninghamia lanceolata (Lamb.) Hook (Chinese fir) is an important tree, commercially and ecologically, in southern China. The traditional regenerating methods are based on organogenesis and cutting propagation. Here, we report the development of a high-frequency somatic embryogenesis (SE) regeneration system synchronized via a liquid culture from immature zygotic embryos. Following synchronization, PEM II cell aggregates were developmentally equivalent in appearance to cleaved zygotic embryos. Embryo and suspensor growth and subsequent occurrence of the apical and then the cotyledonary meristems were similar for zygotic and SE embryo development. However, SE proembryos exhibited a more reddish coloration than zygotic proembryos, and SE embryos were smaller than zygotic embryos. Mature somatic embryos gave rise to plantlets on hormone-free medium. For juvenile explants, low concentrations of endogenous indole-3-acetic acid in initial explants correlated with improved proembryogenic mass formation, and high SE competency. Analysis of karyotypes and microsatellites detected no major genetic variation in the plants regenerated via SE, and suggest a potential in the further development of this system as a reliable methodology for true-to-type seedling production. Treatment with polyethylene glycol (PEG) and abscisic acid (ABA) were of great importance to proembryo formation and complemented each other. ABA assisted the growth of embryonal masses, whereas PEG facilitated the organization of the proembryo-like structures. SOMATIC EMBRYOGENESIS RECEPTOR KINASE SERK) and the WUSCHEL homeobox (WOX) transcription factor served as molecular markers during early embryogenesis. Our results show that ClSERKs are conserved and redundantly expressed during SE. SERK and WOX transcript levels were highest during development of the proembryos and lowest in developed embryos. ClWOX13 expression correlates with the critical transition from proembryogenic masses to proembryos. Both SERK and WOX expression reveal their applicability in Chinese fir as markers of early embryogenesis. Overall, the findings provided evidence for the potential of this system in high fidelity Chinese fir seedlings production. Also, SE modification strategies were demonstrated and could be applied in other conifer species on the basis of our hormonal, morphological and molecular analyses.

Transcriptome and metabolite analyses reveal the complex metabolic genes involved in volatile terpenoid biosynthesis in garden sage (Salvia officinalis)

A large number of terpenoid compounds have been extracted from different tissues of S. officinalis. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 6.6 Gb of raw data were generated from the transcriptome of S. officinalis leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 21 million, comprising 98 million of high-quality nucleotide bases. 48,671 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 135 unigenes that are putatively involved in terpenoid metabolism, including 70 mevalonate and methyl-erythritol phosphate pathways, terpenoid backbone biosynthesis genes, and 65 terpene synthase genes. Moreover, five terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic tobacco; most transgenic tobacco plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. officinalis.

Identification of critical genes associated with lignin biosynthesis in radish (Raphanus sativus L.) by de novo transcriptome sequencing

Radish is an important root vegetable crop with high nutritional, economic, and medicinal value. Lignin is an important secondary metabolite possessing a great effect on plant growth and product quality. To date, lignin biosynthesis-related genes have been identified in some important plant species. However, little information on characterization of critical genes involved in plant lignin biosynthesis is available in radish. In this study, a total of 71,148 transcripts sequences were obtained from radish root, of which 66 assembled unigenes and ten candidate genes were identified to be involved in lignin monolignol biosynthesis. Full-length cDNA sequences of seven randomly selected genes were isolated and sequenced from radish root, and the assembled unigenes covered more than 80% of their corresponding cDNA sequences. Moreover, the lignin content gradually accumulated in leaf during the developmental stages, and it increased from pre-cortex to cortex splitting stage, followed by a decrease at thickening stage and then increased at mature stage in root. RT-qPCR analysis revealed that all these genes except RsF5H exhibited relatively low expression level in root at thickening stage. The expression profiles of Rs4CL5, RsCCoAOMT1, and RsCOMT genes were consistent with the changes of root lignin content, implying that these candidate genes may play important roles in lignin formation in radish root. These findings would provide valuable information for identification of lignin biosynthesis-related genes and facilitate dissection of molecular mechanism underlying lignin biosynthesis in radish and other root vegetable crops.

Investigating right ovary degeneration in chick embryos by transcriptome sequencing

In asymmetric chick gonads, the left and right female gonads undergo distinct programs during development, generating a functional ovary on the left side only. Despite some progress being made in recent years, the mechanisms of molecular regulation remain incompletely understood, and little genomic information is available regarding the degeneration of the right ovary in the chick embryo testis. In this study, we performed transcriptome sequencing to investigate differentially expressed genes in the left and right ovaries and gene functions at two critical time points; embryonic days 6 (E6) and 10 (E10). Using high-throughput RNA-sequencing technologies, 539 and 1046 genes were identified as being significantly differentially expressed between 6R-VS-6L and 10R-VS-10L. Gene ontology analysis of the differentially expressed genes revealed enrichment in functional pathways. Among these, candidate genes associated with degeneration of the right ovary in the chick embryo were identified. Identification of a pathway involved in ovarian degeneration provides an important resource for the further study of its molecular mechanisms and functions.

Comparative transcriptome analysis reveals candidate genes for the biosynthesis of natural insecticide in Tanacetum cinerariifolium

Background

Pyrethrins are monoterpenoids and consist of either a chrysanthemic acid or pyrethric acid with a rethrolone moiety. Natural pyrethrins are safe and eco-friendly while possessing strong insecticidal properties. Despite such advantages of commercial value coming with the eco-friendly tag, most enzymes/genes involved in the pyrethrin biosynthesis pathway remain unidentified and uncharacterized. Since the flowers of Tanacetum cinerariifolium are rich in major pyrethrins, next generation transcriptome sequencing was undertaken to compare the flowers and the leaves of the plant de novo to identify differentially expressed transcripts and ascertain which among them might be involved in and responsible for the differential accumulation of pyrethrins in T. cinerariifolium flowers.

Results

In this first tissue specific transcriptome analysis of the non-model plant T. cinerariifolium, a total of 23,200,000 and 28,500,110 high quality Illumina next generation sequence reads, with a length of 101 bp, were generated for the flower and leaf tissue respectively. After functional enrichment analysis and GO based annotation using public protein databases such as UniRef, PFAM, SMART, KEGG and NR, 4443 and 8901 unigenes were identified in the flower and leaf tissue respectively. These could be assigned to 13344 KEGG pathways and the pyrethrin biosynthesis contextualized. The 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was involved in the biosynthesis of acid moiety of pyrethrin and this pathway predominated in the flowers as compared to the leaves. However, enzymes related to oxylipin biosynthesis were found predominantly in the leaf tissue, which suggested that major steps of pyrethrin biosynthesis occurred in the flowers.

Conclusions

Transcriptome comparison between the flower and leaf tissue of T. cinerariifolium provided an elaborate list of tissue specific transcripts that was useful in elucidating the differences in the expression of the biosynthetic pathways leading to differential presence of pyrethrin in the flowers. The information generated on genes, pathways and markers related to pyrethrin biosynthesis in this study will be helpful in enhancing the production of these useful compounds for value added breeding programs. Related proteome comparison to overlay our transcriptome comparison can generate more relevant information to better understand flower specific accumulation of secondary metabolites in general and pyrethrin accumulation in particular.

Electronic supplementary material

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

De novo Transcriptome Assembly of Chinese Kale and Global Expression Analysis of Genes Involved in Glucosinolate Metabolism in Multiple Tissues

Chinese kale, a vegetable of the cruciferous family, is a popular crop in southern China and Southeast Asia due to its high glucosinolate content and nutritional qualities. However, there is little research on the molecular genetics and genes involved in glucosinolate metabolism and its regulation in Chinese kale. In this study, we sequenced and characterized the transcriptomes and expression profiles of genes expressed in 11 tissues of Chinese kale. A total of 216 million 150-bp clean reads were generated using RNA-sequencing technology. From the sequences, 98,180 unigenes were assembled for the whole plant, and 49,582~98,423 unigenes were assembled for each tissue. Blast analysis indicated that a total of 80,688 (82.18%) unigenes exhibited similarity to known proteins. The functional annotation and classification tools used in this study suggested that genes principally expressed in Chinese kale, were mostly involved in fundamental processes, such as cellular and molecular functions, the signal transduction, and biosynthesis of secondary metabolites. The expression levels of all unigenes were analyzed in various tissues of Chinese kale. A large number of candidate genes involved in glucosinolate metabolism and its regulation were identified, and the expression patterns of these genes were analyzed. We found that most of the genes involved in glucosinolate biosynthesis were highly expressed in the root, petiole, and in senescent leaves. The expression patterns of ten glucosinolate biosynthetic genes from RNA-seq were validated by quantitative RT-PCR in different tissues. These results provided an initial and global overview of Chinese kale gene functions and expression activities in different tissues.

De novo transcriptome sequencing of black pepper (Piper nigrum L.) and an analysis of genes involved in phenylpropanoid metabolism in response to Phytophthora capsici

BACKGROUND

Piper nigrum L., or "black pepper", is an economically important spice crop in tropical regions. Black pepper production is markedly affected by foot rot disease caused by Phytophthora capsici, and genetic improvement of black pepper is essential for combating foot rot diseases. However, little is known about the mechanism of anti- P. capsici in black pepper. The molecular mechanisms underlying foot rot susceptibility were studied by comparing transcriptome analysis between resistant (Piper flaviflorum) and susceptible (Piper nigrum cv. Reyin-1) black pepper species.

RESULTS

116,432 unigenes were acquired from six libraries (three replicates of resistant and susceptible black pepper samples), which were integrated by applying BLAST similarity searches and noted by adopting Kyoto Encyclopaedia of Genes and Gene Ontology (GO) genome orthology identifiers. The reference transcriptome was mapped using two sets of digital gene expression data. Using GO enrichment analysis for the differentially expressed genes, the majority of the genes associated with the phenylpropanoid biosynthesis pathway were identified in P. flaviflorum. In addition, the expression of genes revealed that after susceptible and resistant species were inoculated with P. capsici, the majority of genes incorporated in the phenylpropanoid metabolism pathway were up-regulated in both species. Among various treatments and organs, all the genes were up-regulated to a relatively high degree in resistant species. Phenylalanine ammonia lyase and peroxidase enzyme activity increased in susceptible and resistant species after inoculation with P. capsici, and the resistant species increased faster. The resistant plants retain their vascular structure in lignin revealed by histochemical analysis.

CONCLUSIONS

Our data provide critical information regarding target genes and a technological basis for future studies of black pepper genetic improvements, including transgenic breeding.

De Novo Transcriptome Assembly and Annotation of the Leaves and Callus of Cyclocarya Paliurus (Bata1) Iljinskaja

Cyclocarya Paliurus (Bata1) Iljinskaja contains various bioactive secondary metabolites especially in leaves, such as triterpenes, flavonoids, polysaccharides and alkaloids, and its leaves are widely used as an hyperglycemic tea in China. In the present paper, we sequenced the transcriptome of the leaves and callus of Cyclocarya Paliurus using Illumina Hiseq 4000 platform. After sequencing and de novo assembly, a total of 65,654 unigenes were generated with an N50 length of 1,244bp. Among them, 35,041 (53.37%) unigenes were annotated in NCBI Non-Redundant database, 19,453 (29.63%) unigenes were classified into Gene Ontology (GO) database, and 7,259 (11.06%) unigenes were assigned to Clusters of Orthologous Group (COG) categories. Furthermore, 11,697 (17.81%) unigenes were mapped onto 335 pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG), among which 1,312 unigenes were identified to be involved in biosynthesis of secondary metabolites. In addition, a total of 11,247 putative simple sequence repeats (SSRs) were detected. This transcriptome dataset provides a comprehensive sequence resource for gene expression profiling, genetic diversity, evolution and further molecular genetics research on Cyclocarya Paliurus.

Screening and Validation of Housekeeping Genes of the Root and Cotyledon of Cunninghamia lanceolata under Abiotic Stresses by Using Quantitative Real-Time PCR

Cunninghamia lanceolata (Chinese fir) is a fast-growing and commercially important conifer of the Cupressaceae family. Due to the unavailability of complete genome sequences and relatively poor genetic background information of the Chinese fir, it is necessary to identify and analyze the expression levels of suitable housekeeping genes (HKGs) as internal reference for precise analysis. Based on the results of database analysis and transcriptome sequencing, we have chosen five candidate HKGs (Actin, GAPDH, EF1a, 18S rRNA, and UBQ) with conservative sequences in the Chinese fir and related species for quantitative analysis. The expression levels of these HKGs in roots and cotyledons under five different abiotic stresses in different time intervals were measured by qRT-PCR. The data were statistically analyzed using the following algorithms: NormFinder, BestKeeper, and geNorm. Finally, RankAggreg was applied to merge the sequences generated from three programs and rank these according to consensus sequences. The expression levels of these HKGs showed variable stabilities under different abiotic stresses. Among these, Actin was the most stable internal control in root, and GAPDH was the most stable housekeeping gene in cotyledon. We have also described an experimental procedure for selecting HKGs based on the de novo sequencing database of other non-model plants.

Transcriptome analysis of immature xylem in the Chinese fir at different developmental phases

Background.Chinese fir [Cunninghamia lanceolata (Lamb.) Hook.] is one of the most important native tree species for timber production in southern China. An understanding of overall fast growing stage, stem growth stage and senescence stage cambium transcriptome variation is lacking. We used transcriptome sequencing to identify the repertoire of genes expressed during development of xylem tissue in Chinese fir, aiming to delineate the molecular mechanisms of wood formation. Results. We carried out transcriptome sequencing at three different cultivation ages (7Y, 15Y and 21Y) generating 68.71 million reads (13.88 Gbp). A total of 140,486 unigenes with a mean size of 568.64 base pairs (bp) were obtained via de novo assembly. Of these, 27,427 unigenes (19.52%) were further annotated by comparison to public protein databases. A total of 5,331 (3.79%) unigenes were mapped into 118 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Differentially expressed genes (DEG) analysis identified 3, 16 and 5,899 DEGs from the comparison of 7Y vs. 15Y, 7Y vs. 21Y and 15Y vs. 21Y, respectively, in the immature xylem tissues, including 2,638 significantly up-regulated and 3,280 significantly down-regulated genes. Besides, five NAC transcription factors, 190 MYB transcription factors, and 34 WRKY transcription factors were identified respectively from Chinese fir transcriptome. Conclusion. Our results revealed the active transcriptional pathways and identified the DEGs at different cultivation phases of Chinese fir wood formation. This transcriptome dataset will aid in understanding and carrying out future studies on the molecular basis of Chinese fir wood formation and contribute to future artificial production and applications.

Transcriptomic analysis reveals the gene expression profile that specifically responds to IBA during adventitious rooting in mung bean seedlings

BACKGROUND

Auxin plays a critical role in inducing adventitious rooting in many plants. Indole-3-butyric acid (IBA) is the most widely employed auxin for adventitious rooting. However, the molecular mechanisms by which auxin regulate the process of adventitious rooting are less well known.

RESULTS

The RNA-Seq data analysis indicated that IBA treatment greatly increased the amount of clean reads and the amount of expressed unigenes by 24.29 % and 27.42 % and by 4.3 % and 5.04 % at two time points, respectively, and significantly increased the numbers of unigenes numbered with RPKM = 10-100 and RPKM = 500-1000 by 13.04 % and 3.12 % and by 24.66 % and 108.2 % at two time points, respectively. Gene Ontology (GO) enrichment analysis indicated that the enrichment of down-regulated GOs was 2.87-fold higher than that of up-regulated GOs at stage 1, suggesting that IBA significantly down-regulated gene expression at 6 h. The GO functional category indicated that IBA significantly up- or down-regulated processes associated with auxin signaling, ribosome assembly and protein synthesis, photosynthesis, oxidoreductase activity and extracellular region, secondary cell wall biogenesis, and the cell wall during the development process. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment indicated that ribosome biogenesis, plant hormone signal transduction, pentose and glucuronate interconversions, photosynthesis, phenylpropanoid biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, ribosome, cutin, flavonoid biosynthesis, and phenylalanine metabolism were the pathways most highly regulated by IBA. A total of 6369 differentially expressed (2-fold change > 2) unigenes (DEGs) with 3693 (58 %) that were up-regulated and 2676 (42 %) down-regulated, 5433 unigenes with 2208 (40.6 %) that were up-regulated and 3225 (59.4 %) down-regulated, and 7664 unigenes with 3187 (41.6 %) that were up-regulated and 4477 (58.4 %) down-regulated were detected at stage 1, stage 2, and between stage 1 and stage 2, respectively, suggesting that IBA treatment increased the number of DEGs. A total of 143 DEGs specifically involved in plant hormone signaling and 345 transcription factor (TF) genes were also regulated by IBA. qRT-PCR validation of the 36 genes with known functions indicated a strong correlation with the RNA-Seq data.

CONCLUSIONS

The changes in GO functional categories, KEGG pathways, and global DEG profiling during adventitious rooting induced by IBA were analyzed. These results provide valuable information about the molecular traits of IBA regulation of adventitious rooting.

Selecting Superior De Novo Transcriptome Assemblies: Lessons Learned by Leveraging the Best Plant Genome

Whereas de novo assemblies of RNA-Seq data are being published for a growing number of species across the tree of life, there are currently no broadly accepted methods for evaluating such assemblies. Here we present a detailed comparison of 99 transcriptome assemblies, generated with 6 de novo assemblers including CLC, Trinity, SOAP, Oases, ABySS and NextGENe. Controlled analyses of de novo assemblies for Arabidopsis thaliana and Oryza sativa transcriptomes provide new insights into the strengths and limitations of transcriptome assembly strategies. We find that the leading assemblers generate reassuringly accurate assemblies for the majority of transcripts. At the same time, we find a propensity for assemblers to fail to fully assemble highly expressed genes. Surprisingly, the instance of true chimeric assemblies is very low for all assemblers. Normalized libraries are reduced in highly abundant transcripts, but they also lack 1000s of low abundance transcripts. We conclude that the quality of de novo transcriptome assemblies is best assessed through consideration of a combination of metrics: 1) proportion of reads mapping to an assembly 2) recovery of conserved, widely expressed genes, 3) N50 length statistics, and 4) the total number of unigenes. We provide benchmark Illumina transcriptome data and introduce SCERNA, a broadly applicable modular protocol for de novo assembly improvement. Finally, our de novo assembly of the Arabidopsis leaf transcriptome revealed ~20 putative Arabidopsis genes lacking in the current annotation.

Comparative Transcriptome Analyses Reveal a Special Glucosinolate Metabolism Mechanism in Brassica alboglabra Sprouts

Brassica sprouts contain abundant phytochemicals, especially glucosinolates (GSs). Various methods have been used to enhance GS content in sprouts. However, the molecular basis of GS metabolism in sprouts remains an open question. Here we employed RNA-seq analysis to compare the transcriptomes of high-GS (JL-08) and low-GS (JL-09) Brassica alboglabra sprouts. Paired-end Illumina RNA-seq reads were generated and mapped to the Brassica oleracea reference genome. The differentially expressed genes were analyzed between JL-08 and JL-09. Among these, 1477 genes were up-regulated and 1239 down-regulated in JL-09 compared with JL-08. Enrichment analysis of these differentially expressed genes showed that the GS biosynthesis had the smallest enrichment factor and the highest Q-value of all metabolic pathways in Kyoto Encyclopedia of Genes and Genomes database, indicating the main metabolic difference between JL-08 and JL-09 is the GS biosynthetic pathway. Thirty-seven genes of the sequenced data were annotated as putatively involved in GS biosynthesis, degradation, and regulation, of which 11 were differentially expressed in JL-08 and JL-09. The expression level of GS degradation enzyme myrosinase in high-GS JL-08 was lower compared with low-GS JL-09. Surprisingly, in high-GS JL-08, the expression levels of GS biosynthesis genes were also lower than those in low-GS JL-09. As the GS contents in sprouts are determined by dynamic equilibrium of seed stored GS mobilization, de novo synthesis, degradation, and extra transport, the result of this study leads us to suggest that efforts to increase GS content should focus on either raising GS content in seeds or decreasing myrosinase activity, rather than improving the expression level of GS biosynthesis genes in sprouts.

Genome-Wide Transcriptional Profiling and Metabolic Analysis Uncover Multiple Molecular Responses of the Grass Species Lolium perenne Under Low-Intensity Xenobiotic Stress

Lolium perenne, which is a major component of pastures, lawns, and grass strips, can be exposed to xenobiotic stresses due to diffuse and residual contaminations of soil. L. perenne was recently shown to undergo metabolic adjustments in response to sub-toxic levels of xenobiotics. To gain insight in such chemical stress responses, a de novo transcriptome analysis was carried out on leaves from plants subjected at the root level to low levels of xenobiotics, glyphosate, tebuconazole, and a combination of the two, leading to no adverse physiological effect. Chemical treatments influenced significantly the relative proportions of functional categories and of transcripts related to carbohydrate processes, to signaling, to protein-kinase cascades, such as Serine/Threonine-protein kinases, to transcriptional regulations, to responses to abiotic or biotic stimuli and to responses to phytohormones. Transcriptomics-based expressions of genes encoding different types of SNF1 (sucrose non-fermenting 1)-related kinases involved in sugar and stress signaling or encoding key metabolic enzymes were in line with specific qRT-PCR analysis or with the important metabolic and regulatory changes revealed by metabolomic analysis. The effects of pesticide treatments on metabolites and gene expression strongly suggest that pesticides at low levels, as single molecule or as mixture, affect cell signaling and functioning even in the absence of major physiological impact. This global analysis of L. perenne therefore highlighted the interactions between molecular regulation of responses to xenobiotics, and also carbohydrate dynamics, energy dysfunction, phytohormones and calcium signaling.

Revealing crosstalk of plant and fungi in the symbiotic roots of sewage-cleaning Eichhornia crassipes using direct de novo metatranscriptomic analysis

Cultivation and environmental changes can induce development of novel phenotypes in plants. For example, the root morphology of cultivated purple root Eichhornia crassipes differs remarkably from normal Eichhornia crassipes and also shows an enhanced ability to absorb heavy metal from groundwater. However, the changes in gene expression associated with these processes are unknown because of the lack of information on its large and unsequenced genome and its complex plant-rhizosphere symbiotic system. To investigate these gene expression changes, we applied a new strategy, direct de novo metatranscriptome analysis. Using this approach, we assembled the metatranscriptome of the entire rhizosphere and identified species-specific differentially expressed genes (DEGs) via hyper-accurate algorithms, showing a polarized plant/fungus distribution: the plant genes were responsible for morphological changes to the root system, offering a greater volume and surface area that hosts more fungi; while genes associated with heavy metal response in the fungus Fusarium were upregulated more than 3600-fold. These results suggested a distinct and synergistic functional response by the plant and fungal transcriptomes, indicating significant plant/fungal crosstalk during environmental changes. This study demonstrates that the metatranscriptomic approach adopted here offers a cost-efficient strategy to study symbiosis systems without the need for a priori genomic knowledge.

Transcriptome sequencing and analysis of major genes involved in calcium signaling pathways in pear plants (Pyrus calleryana Decne.)

Background

Pears (Pyrus spp. L.) are an important genus of trees that produce one of the world’s oldest fruit crops. Salinity stress is a common limiting factor for plant productivity that significantly affects the flavor and nutritional quality of pear fruits. Much research has shown that calcium signaling pathways, mediated by Calcineurin B-like proteins (CBLs) and their interacting kinases (CIPKs), are closely associated with responses to stresses, including salt. However, little is known about the molecular mechanisms that govern the relationship between salt stress and calcium signaling pathways in pear plants. The available genomic information for pears has promoted much functional genomic analysis and molecular breeding of the genus. This provided an ample foundation for characterizing the transcriptome of pear under salt stress.

Results

A high-throughput Illumina RNA-seq technology was used to identify a total of 78,695 unigenes that were successfully annotated by BLASTX analysis, using the publicly available protein database. Additionally, 2,855 novel transcripts, 218,167 SNPs, 23,248 indels and 18,322 alternative splicing events occurred. Assembled unique sequences were annotated and classified with Gene Ontology (GO), Clusters of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, which revealed that the main activated genes in pear are predominately involved in functions such as basic physiological processes, metabolic pathways, operation of cellular components, signal transduction mechanisms, and other molecular activities. Through targeted searches of the annotations, the majority of the genes involved in calcium signaling pathways were identified, among which, four genes were validated by molecular cloning, while 11 were validated by RT-qPCR expression profiles under salt stress treatment.

Conclusions

These results facilitate a better understanding of the molecular genetics and functional genomic mechanisms of salt stress in pear plants. Furthermore, they provide a valuable foundation for additional research on the molecular biology and functional genomics of pear and related species.

Electronic supplementary material

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

De Novo Characterization of the Mung Bean Transcriptome and Transcriptomic Analysis of Adventitious Rooting in Seedlings Using RNA-Seq

Adventitious rooting is the most important mechanism underlying vegetative propagation and an important strategy for plant propagation under environmental stress. The present study was conducted to obtain transcriptomic data and examine gene expression using RNA-Seq and bioinformatics analysis, thereby providing a foundation for understanding the molecular mechanisms controlling adventitious rooting. Three cDNA libraries constructed from mRNA samples from mung bean hypocotyls during adventitious rooting were sequenced. These three samples generated a total of 73 million, 60 million, and 59 million 100-bp reads, respectively. These reads were assembled into 78,697 unigenes with an average length of 832 bp, totaling 65 Mb. The unigenes were aligned against six public protein databases, and 29,029 unigenes (36.77%) were annotated using BLASTx. Among them, 28,225 (35.75%) and 28,119 (35.62%) unigenes had homologs in the TrEMBL and NCBI non-redundant (Nr) databases, respectively. Of these unigenes, 21,140 were assigned to gene ontology classes, and a total of 11,990 unigenes were classified into 25 KOG functional categories. A total of 7,357 unigenes were annotated to 4,524 KOs, and 4,651 unigenes were mapped onto 342 KEGG pathways using BLAST comparison against the KEGG database. A total of 11,717 unigenes were differentially expressed (fold change>2) during the root induction stage, with 8,772 unigenes down-regulated and 2,945 unigenes up-regulated. A total of 12,737 unigenes were differentially expressed during the root initiation stage, with 9,303 unigenes down-regulated and 3,434 unigenes up-regulated. A total of 5,334 unigenes were differentially expressed between the root induction and initiation stage, with 2,167 unigenes down-regulated and 3,167 unigenes up-regulated. qRT-PCR validation of the 39 genes with known functions indicated a strong correlation (92.3%) with the RNA-Seq data. The GO enrichment, pathway mapping, and gene expression profiles reveal molecular traits for root induction and initiation. This study provides a platform for functional genomic research with this species.

RNA-seq Reveals Complicated Transcriptomic Responses to Drought Stress in a Nonmodel Tropic Plant, Bombax ceiba L

High-throughput transcriptome provides an unbiased approach for understanding the genetic basis and gene functions in response to different conditions. Here we sequenced RNA-seq libraries derived from a Bombax ceiba L. system under a controlled experiment. As a known medicinal and ornamental plant, B. ceiba grows mainly in hot-dry monsoon rainforests in Southeast Asia and Australia. Due to the specific growth environment, it has evolved a unique system that enables a physiologic response to drought stress. To date, few studies have characterized the genome-wide features of drought endurance in B. ceiba. In this study, we first attempted to characterize and identify the most differentially expressed genes and associated functional pathways under drought treatment and normal condition. Using RNA-seq technology, we generated the first transcriptome of B. ceiba and identified 59 differentially expressed genes with greater than 1,000-fold changes under two conditions. The set of upregulated genes implicates interplay among various pathways: plants growth, ubiquitin-mediated proteolysis, polysaccharides hydrolyzation, oxidative phosphorylation and photosynthesis, etc. In contrast, genes associated with stem growth, cell division, fruit ripening senescence, disease resistance, and proline synthesis are repressed. Notably, key genes of high RPKM levels in drought are AUX1, JAZ, and psbS, which are known to regulate the growth of plants, the resistance against abiotic stress, and the photosynthesis process. Furthermore, 16,656 microsatellite markers and 3,071 single-nucleotide polymorphisms (SNPs) were predicted by in silico methods. The identification and functional annotation of differentially expressed genes, microsatellites, and SNPs represent a major step forward and would serve as a valuable resource for understanding the complexity underlying drought endurance and adaptation in B. ceiba.

Isolation of High-Quality Total RNA from Chinese Fir (Cunninghamia lanceolata (Lamb.) Hook)

RNA isolation with RNA in a high quantity is a basic analytical method in plant genetics, molecular biology and related physiological investigations. To understand the genetic and molecular biology of Chinese fir, sufficient high-quality total RNA must be obtained for cDNA library construction and other downstream molecular applications. However, extracting RNA from Chinese fir is difficult and often requires the modification of existing protocols. Chinese fir tissues containing large amounts of polysaccharides and polyphenol compounds and are one of the most difficult plant tissues for RNA isolation. Therefore, we developed a simple method for extracting high-quality RNA from Chinese fir tissues. RNA isolations were performed within two hours, RNA quality was measured for yield and purity. Total RNA obtained from this procedure was successfully used for cDNA library construction, RT-PCR and transcriptome sequencing. It was proven that extracted RNA was intact and suitable for downstream molecular applications, including RT-PCR and qPCR, and other downstream molecular applications. Thus, this protocol represents a simple, efficient, and low-cost method.

De Novo Transcriptome Sequencing of Low Temperature-Treated Phlox subulata and Analysis of the Genes Involved in Cold Stress

Phlox subulata, a perennial herbaceous flower, can survive during the winter of northeast China, where the temperature can drop to −30 °C, suggesting that P. subulata is an ideal model for studying the molecular mechanisms of cold acclimation in plants. However, little is known about the gene expression profile of P. subulata under cold stress. Here, we examined changes in cold stress-related genes in P. subulata. We sequenced three cold-treated (CT) and control (CK) samples of P. subulata. After denovo assembly and quantitative assessment of the obtained reads, 99,174 unigenes were generated. Based on similarity searches with known proteins in public protein databases, 59,994 unigenes were functionally annotated. Among all differentially expressed genes (DEGs), 8302, 10,638 and 11,021 up-regulated genes and 9898, 17,876, and 12,358 down-regulated genes were identified after treatment at 4, 0, and −10 °C, respectively. Furthermore, 3417 up-regulated unigenes were expressed only in CT samples. Twenty major cold-related genes, including transcription factors, antioxidant enzymes, osmoregulation proteins, and Ca²⁺ and ABA signaling components, were identified, and their expression levels were estimated. Overall, this is the first transcriptome sequencing of this plant species under cold stress. Studies of DEGs involved in cold-related metabolic pathways may facilitate the discovery of cold-resistance genes.

De novo characterization of fall dormant and nondormant alfalfa (Medicago sativa L.) leaf transcriptome and identification of candidate genes related to fall dormancy

Alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage legumes worldwide. Fall dormancy is an adaptive character related to the biomass production and winter survival in alfalfa. The physiological, biochemical and molecular mechanisms causing fall dormancy and the related genes have not been well studied. In this study, we sequenced two standard varieties of alfalfa (dormant and non-dormant) at two time points and generated approximately 160 million high quality paired-end sequence reads using sequencing by synthesis (SBS) technology. The de novo transcriptome assembly generated a set of 192,875 transcripts with an average length of 856 bp representing about 165.1 Mb of the alfalfa leaf transcriptome. After assembly, 111,062 (57.6%) transcripts were annotated against the NCBI non-redundant database. A total of 30,165 (15.6%) transcripts were mapped to 323 Kyoto Encyclopedia of Genes and Genomes pathways. We also identified 41,973 simple sequence repeats, which can be used to generate markers for alfalfa, and 1,541 transcription factors were identified across 1,350 transcripts. Gene expression between dormant and non-dormant alfalfa at different time points were performed, and we identified several differentially expressed genes potentially related to fall dormancy. The Gene Ontology and pathways information were also identified. We sequenced and assembled the leaf transcriptome of alfalfa related to fall dormancy, and also identified some genes of interest involved in the fall dormancy mechanism. Thus, our research focused on studying fall dormancy in alfalfa through transcriptome sequencing. The sequencing and gene expression data generated in this study may be used further to elucidate the complete mechanisms governing fall dormancy in alfalfa.

De novo characterization of the Lycium chinense Mill. leaf transcriptome and analysis of candidate genes involved in carotenoid biosynthesis

Lycium chinense Mill. (Chinese wolfberry), enriching in carotenoids, is an important Chinese herbal medicine. However, studies on the functional genomics research, especially the carotenoid biosynthesis and accumulation, are limited because of insufficiently available datasets. RNA-Seq was performed by the Illumina sequencing platform. Approximately 26 million clean reads were generated after filtering. Clean reads were assembled by SOAPdenovo and subsequently annotated. Among all 61,595 unigenes, 37,816 (61.39%), 25,266 (41.02%), and 17,598 (28.57%) unigenes were annotated in NCBI non-redundant protein, Swiss-Prot, and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, respectively. A total of 16,073 and 11,394 unigenes were assigned to Gene Ontology and Cluster of Orthologous Group, respectively. Furthermore, the majority of genes encoding the enzymes in the carotenoid biosynthesis pathway were identified in the unigene datasets. We first found several genes related to L. chinense carotenoid biosynthesis. The expression levels and the biological functions of these genes involved in carotenoid biosynthesis in the leaf and the green ripening fruit were further confirmed by qPCR and high performance liquid chromatography (HPLC). In the present study, we first characterized the transcriptome of L. chinense leaf, which may provide useful data for functional genomics investigations in L. chinense in the future. And essential genes involved in the carotenoid biosynthesis pathway may contribute to elucidate the expression patterns in different stages of development and fruit ripening and the specific mechanisms of carotenoid biosynthesis/accumulation in L. chinense.

De novo transcriptome analysis of the Siberian apricot (Prunus sibirica L.) and search for potential SSR markers by 454 pyrosequencing

The Siberian apricot, an economically and ecologically important plant in China, contains seeds high in oil and can grow on marginal land. Although this species has multiple purposes and may be a feedstock of biofuel in China, transcriptome information and molecular research on this species remain limited. RNA-Seq technology has been widely applied to transcriptomics, genomics and the development of molecular markers, and functional gene studies. In this study, we obtained 1,243,067 high-quality reads with a mean size of 425 bp in a single run, totaling 528.4 Mb of sequence data using 454 GS FLX Titanium sequencing. All reads were assembled de novo into 46,940 unigenes with a mean size of 651 bp (range: 45-5566 bp). Assembled unigenes were annotated in multiple public databases based on similarity alignments to genes and proteins. 191 unigenes involving in lipid biosynthesis and metabolism were found, among them, expression patterns of two desaturase enzymes were analyzed by quantitative real-time polymerase chain reaction (qRT-PCR), based on six tissues from Siberian apricot, the seeds had the highest expression. 7304 simple sequence repeats (SSR) were identified from 6509 unigenes, a total of 9930 primer pairs were designed, 50 primer pairs were randomly selected to validate of the usefulness, and 24 (48%) primer pairs produced bands of the expected size. These data provide a base of sequence information to improve agronomic characters and molecular marker-assisted breeding to alter the composition of fatty acids in seeds from this plant, and hence, facilitate its utilization as a future biodiesel feedstock.

De novo transcriptome sequencing of radish (Raphanus sativus L.) and analysis of major genes involved in glucosinolate metabolism

Background

Radish (Raphanus sativus L.), is an important root vegetable crop worldwide. Glucosinolates in the fleshy taproot significantly affect the flavor and nutritional quality of radish. However, little is known about the molecular mechanisms underlying glucosinolate metabolism in radish taproots. The limited availability of radish genomic information has greatly hindered functional genomic analysis and molecular breeding in radish.

Results

In this study, a high-throughput, large-scale RNA sequencing technology was employed to characterize the de novo transcriptome of radish roots at different stages of development. Approximately 66.11 million paired-end reads representing 73,084 unigenes with a N50 length of 1,095 bp, and a total length of 55.73 Mb were obtained. Comparison with the publicly available protein database indicates that a total of 67,305 (about 92.09% of the assembled unigenes) unigenes exhibit similarity (e –value ≤ 1.0e^-5) to known proteins. The functional annotation and classification including Gene Ontology (GO), Clusters of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the main activated genes in radish taproots are predominately involved in basic physiological and metabolic processes, biosynthesis of secondary metabolite pathways, signal transduction mechanisms and other cellular components and molecular function related terms. The majority of the genes encoding enzymes involved in glucosinolate (GS) metabolism and regulation pathways were identified in the unigene dataset by targeted searches of their annotations. A number of candidate radish genes in the glucosinolate metabolism related pathways were also discovered, from which, eight genes were validated by T-A cloning and sequencing while four were validated by quantitative RT-PCR expression profiling.

Conclusions

The ensuing transcriptome dataset provides a comprehensive sequence resource for molecular genetics research in radish. It will serve as an important public information platform to further understanding of the molecular mechanisms involved in biosynthesis and metabolism of the related nutritional and flavor components during taproot formation in radish.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-836) contains supplementary material, which is available to authorized users.