Gene transcript profiles in the desert plant Nitraria tangutorum during fruit development and ripening

Nitraria tangutorum Bobr., a valuable wild shrub distributed in Northwest China, produces edible and medicinal berries. However, little is known about the molecular mechanisms of its fruit development and ripening. We performed de novo transcriptome sequencing of N. tangutorum fruit using the Illumina HiSeq™ 2000 sequencing platform. More than 62.94 million reads were obtained and assembled into 69,306 unigenes (average length, 587 bp). These unigenes were annotated by querying against five databases (Nr, Swiss-Prot, GO, COG, and KEGG); 42,929 and 26,809 unigenes were found in the Nr and Swiss-Prot databases, respectively. In ortholog analyses, 33,363 unigenes were assigned with one or more GO terms, 15,537 hits were aligned to 25 COG classes, and 24,592 unigenes were classified into 128 KEGG pathways. Digital gene expression analyses were conducted on N. tangutorum fruit at the green (S1), yellow (S2), and red (S3) developmental stages. In total, 8240, 5985, and 4994 differentially expressed genes (DEGs) were detected for S1 vs. S2, S1 vs. S3, and S2 vs. S3, respectively. Cluster analyses showed that a large proportion of DEGs related to plant hormones and transcription factors (TFs) showed high expression in S1, down-regulated expression in S2, and up-regulated expression in S3. We analyzed the expression patterns of 23 genes encoding 12 putative enzymes involved in flavonoid biosynthesis. The expression profiles of 10 DEGs involved in flavonoid biosynthesis were validated by Q-PCR analysis. The assembled and annotated transcriptome sequences and gene expression profile analyses provide valuable genetic resources for research on N. tangutorum.