RefFinder: a web-based tool for comprehensively analyzing and identifying reference genes.

Although many genes may serve as reference genes, they may cause different expression patterns by selecting different reference genes because no single gene is expressed consistently in all tested tissues of an organism under all environmental and developmental conditions. Thus, it is becoming increasingly important and necessary to identify suitable reference genes before performing gene expression analysis. Currently, there are several computational tools available for evaluating the stability of candidate reference genes. These tools are based on different statistical algorithms and may produce different rankings in stability within the same reference gene study. To date, the RefFinder is the only web-based tool available for comparing and evaluating housekeeping genes as candidates to be reference genes. In this tool, we integrated the four currently available computational programs (geNorm, NormFinder, BestKeeper, and the comparative ΔCt method) into a web-based tool for evaluating the stability and reliability of reference genes. According to the gene stability rankings derived from the four programs, we assigned an appropriate weight to each gene and calculated the geometric mean of weights for the final rankings. Aside from the overall ranking, a single program or combination of the four programs can be selected for evaluating the ranking of candidate reference genes. This tool has been widely used and validated by many research laboratories around the world. You may use this tool at http://www.heartcure.com.au/reffinder/  or  https://blooge.cn/RefFinder/ . You can also download this algorithm program from https://github.com/fulxie/RefFinder and setup on your own computer. RefFinder is developed by PHP. Users can deploy it to a Php-based server (Apache + PHP) and run it.

Evaluation of reference genes for miRNA and mRNA normalization in tobacco infected with PVYNTN, PVYN-Wi and PVYZ-NTN strains.

This is the first report on identification of the most suitable reference genes for RT-qPCR quantification of miRNA and mRNA in tobacco response to the prevalent recombinant potato virus Y (PVY) strains PVYNTN, PVYN-Wi and the newly identified PVYZ-NTN. Of 10 tested genes, the expression levels of neIF5C, nU2af and nPP2A were the most stable in samples taken from non-inoculated, mock-inoculated, and infected plants at three days post-inoculation (dpi) and 14 dpi. While the homologues of eIF5 were most stably expressed in tobacco in this study and in potato in our previous study (Yin et al., 2021) following inoculation with the same three PVY strains, the homologues of other two genes PP2A and U2af were stably expressed only in tobacco but unstable in potato. The tobacco homologue of PP2A, which was the most stably expressed one in tobacco interaction with PVYNTN, PVYN-Wi and PVYZ-NTN strains in this study, was the least stable one in tobacco interaction with the non-recombinant PVYO strain in a previous study (Baek et al., 2017). This study provides evidence on the influence of host species on expression of housekeeping genes and points out virus strain as a new factor influencing expression stability of reference gene. Caution should be taken when choosing reference genes in gene expression study in Solanaceae hosts response to different PVY strains.

Reference gene selection for miRNA and mRNA normalization in potato in response to potato virus Y.

This was the first report on evaluating candidate reference genes for quantifying the expression profiles of both coding (e.g., mRNA) and non-coding (e.g., miRNA) genes in potato response to potato virus Y (PVY) inoculation. The reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) method was employed to quantify the expression profiles of eight selected candidate reference genes; their expression stability was analyzed by four statistical algorithms, i.e., geNorm, BestKeeper, NormFinder and RefFinder. The most stable reference genes were sEF1a, sTUBb and seIF5 with a high stability. The least stable ones were sPP2A, sSUI1 and sGAPDH. The same reference gene allows for normalization of both miRNA and mRNA levels from a single RNA sample using cDNAs synthesized in a single RT reaction, in which a stem-loop primer was used for miRNAs and the oligo (dT) for mRNAs.

High-throughput transcriptome sequencing reveals extremely high doses of ionizing radiation-response genes in Caenorhabditis elegans.

This study sought novel ionizing radiation-response (IR-response) genes in Caenorhabditis elegans (C. elegans). C. elegans was divided into three groups and exposed to different high doses of IR: 0 gray (Gy), 200 Gy, and 400 Gy. Total RNA was extracted from each group and sequenced. When the transcriptomes were compared among these groups, many genes were shown to be differentially expressed, and these genes were significantly enriched in IR-related biological processes and pathways, including gene ontology (GO) terms related to cellular behaviours, cellular growth and purine metabolism and kyoto encyclopedia of genes and genomes (KEGG) pathways related to ATP binding, GTPase regulator activity, and RNA degradation. Quantitative reverse-transcription PCR (qRT-PCR) confirmed that these genes displayed differential expression across the treatments. Further gene network analysis showed a cluster of novel gene families, such as the guanylate cyclase (GCY), Sm-like protein (LSM), diacylglycerol kinase (DGK), skp1-related protein (SKR), and glutathione S-transferase (GST) gene families which were upregulated. Thus, these genes likely play important roles in IR response. Meanwhile, some important genes that are well known to be involved in key signalling pathways, such as phosphoinositide-specific phospholipase C-3 (PLC-3), phosphatidylinositol 3-kinase age-1 (AGE-1), Raf homolog serine/threonine-protein kinase (LIN-45) and protein cbp-1 (CBP-1), also showed differential expression during IR response, suggesting that IR response might perturb these key signalling pathways. Our study revealed a series of novel IR-response genes in Caenorhabditis elegans that might act as regulators of IR response and represent promising markers of IR exposure.

microRNA response in potato virus Y infected tobacco shows strain-specificity depending on host and symptom severity.

The present study demonstrates how different potato virus Y (PVY) strains affect the miRNA balance in tobacco cv. Samsun. The two prevalent strains PVYNTN and PVYN-Wi caused severe and mild veinal necrosis (VN) respectively, and the unique PVYZ-NTN strain induced milder vein clearing (VCl) in the upper non-inoculated leaves. A single amino acid polymorphisms (SAPs) I252V and a Q412 to R412 substitution in the HC-Pro cistron of the PVYZ-NTN strain might relate to the loss of VN in tobacco. The abundance of 18 out of the 26 tested miRNAs was increased upon infection by the severe strains PVYNTN and PVYN-Wi. Expression of a group of defense related transcripts were increased accordingly. Two miRNAs, nta-miR6020a-5p and nta-miR6164a/b, which target the TIR-NBS-LRR type resistant TMV N genes involving in signal transduction, might correlate with the PVYNTN and PVYN-Wi induced VN. The down-regulated mRNAs, e.g., RAP2-7 and TOE3, PXC3, LRR-RLK, ATHB-14 and TCP4 targeted by nta-miR172, nta-miR390, nta-miR482, nta-miR166 and nta-miR319/159 respectively, were related to regulation of transcription, protein phosphorylation and cell differentiation. The observed strain-specific alteration of miRNAs and their targets are host dependent and corresponds to the symptom severity and the viral HC-Pro RNA levels.

Genome of wild olive and the evolution of oil biosynthesis.

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Identification, characterization, and gene expression analysis of nucleotide binding site (NB)-type resistance gene homologues in switchgrass.

BACKGROUND: Switchgrass (Panicum virgatum L.) is a warm-season perennial grass that can be used as a second generation bioenergy crop. However, foliar fungal pathogens, like switchgrass rust, have the potential to significantly reduce switchgrass biomass yield. Despite its importance as a prominent bioenergy crop, a genome-wide comprehensive analysis of NB-LRR disease resistance genes has yet to be performed in switchgrass. RESULTS: In this study, we used a homology-based computational approach to identify 1011 potential NB-LRR resistance gene homologs (RGHs) in the switchgrass genome (v 1.1). In addition, we identified 40 RGHs that potentially contain unique domains including major sperm protein domain, jacalin-like binding domain, calmodulin-like binding, and thioredoxin. RNA-sequencing analysis of leaf tissue from 'Alamo', a rust-resistant switchgrass cultivar, and 'Dacotah', a rust-susceptible switchgrass cultivar, identified 2634 high quality variants in the RGHs between the two cultivars. RNA-sequencing data from field-grown cultivar 'Summer' plants indicated that the expression of some of these RGHs was developmentally regulated. CONCLUSIONS: Our results provide useful insight into the molecular structure, distribution, and expression patterns of members of the NB-LRR gene family in switchgrass. These results also provide a foundation for future work aimed at elucidating the molecular mechanisms underlying disease resistance in this important bioenergy crop.

Genome-Wide Identification of R2R3-MYB Genes and Expression Analyses During Abiotic Stress in Gossypium raimondii.

The R2R3-MYB is one of the largest families of transcription factors, which have been implicated in multiple biological processes. There is great diversity in the number of R2R3-MYB genes in different plants. However, there is no report on genome-wide characterization of this gene family in cotton. In the present study, a total of 205 putative R2R3-MYB genes were identified in cotton D genome (Gossypium raimondii), that are much larger than that found in other cash crops with fully sequenced genomes. These GrMYBs were classified into 13 groups with the R2R3-MYB genes from Arabidopsis and rice. The amino acid motifs and phylogenetic tree were predicted and analyzed. The sequences of GrMYBs were distributed across 13 chromosomes at various densities. The results showed that the expansion of the G. Raimondii R2R3-MYB family was mainly attributable to whole genome duplication and segmental duplication. Moreover, the expression pattern of 52 selected GrMYBs and 46 GaMYBs were tested in roots and leaves under different abiotic stress conditions. The results revealed that the MYB genes in cotton were differentially expressed under salt and drought stress treatment. Our results will be useful for determining the precise role of the MYB genes during stress responses with crop improvement.

A cotton miRNA is involved in regulation of plant response to salt stress.

The present study functionally identified a new microRNA (microRNA ovual line 5, miRNVL5) with its target gene GhCHR from cotton (Gossypium hirsutum). The sequence of miRNVL5 precursor is 104 nt long, with a well developed secondary structure. GhCHR contains two DC1 and three PHD Cys/His-rich domains, suggesting that GhCHR encodes a zinc-finger domain-containing transcription factor. miRNVL5 and GhCHR express at various developmental stages of cotton. Under salt stress (50-400 mM NaCl), miRNVL5 expression was repressed, with concomitant high expression of GhCHR in cotton seedlings. Ectopic expression of GhCHR in Arabidopsis conferred salt stress tolerance by reducing Na(+) accumulation in plants and improving primary root growth and biomass. Interestingly, Arabidopsis constitutively expressing miRNVL5 showed hypersensitivity to salt stress. A GhCHR orthorlous gene At2g44380 from Arabidopsis that can be cleaved by miRNVL5 was identified by degradome sequencing, but no confidential miRNVL5 homologs in Arabidopsis have been identified. Microarray analysis of miRNVL5 transgenic Arabidopsis showed six downstream genes (CBF1, CBF2, CBF3, ERF4, AT3G22920, and AT3G49200), which were induced by salt stress in wild-type but repressed in miRNVL5-expressing Arabidopsis. These results indicate that miRNVL5 is involved in regulation of plant response to salt stress.

Identification and characterization of microRNAs in the plant parasitic root-knot nematode Meloidogyne incognita using deep sequencing.

The root-knot nematode Meloidogyne incognita is among the most damaging plant-parasitic pests of several crops including cotton (Gossypium hirsutum) and tomato (Lycopersicon escultentum). Recently, a genome has become available for M. incognita, which greatly facilitates investigation of the interactions between M. incognita and its plant hosts at the molecular level and enables formation of hypotheses concerning development at the cellular level. MicroRNAs (miRNAs) are a class of small RNA molecules that serve as endogenous gene regulators. They regulate many biological processes including reproduction, the sequencing of morphological development, and potentially of parasitism as well. Certain miRNAs regulate fundamental metabolism pathways and stress responses in M. incognita. Since a list of miRNAs has not been generated for M. incognita, we employed a bioinformatics tool called mirDeepFinder to identify miRNAs from the small RNA database of M. incognita (GSM611102) that was generated from deep sequencing. A total of 254 conserved miRNAs belonging to 161 miRNA families were identified, as were 35 novel miRNAs belonging to 31 families. The 16 most commonly found miRNAs in order of abundance were min-miR-100a, min-miR-124, min-miR-71a, min-miR-1, min-miR-228, min-miR-92, min-miR-72, min-miR-49b, min-miR-58, min-miR-252, min-miR-lin-4, min-miR-87, min-miR-2a, min-miR-34a, min-miR-50a, and min-miR-279a. The length of the pre-miRNAs varied greatly from 50 to 197 nt, with an average of 88 ± 39 nt. The average minimal folding free energy (MFE) and MFE index (MFEI) of the identified miRNAs were -30.3 Kcal/mol and 0.92, respectively, indicating that these miRNAs can readily fold into a typical hairpin secondary structure.

microRNA evolution and expression analysis in polyploidized cotton genome.

Cotton (Gossypium hirsutum L.), the most important fibre plant in the world, is a tetraploid species, originating from the reunion of two ancestral cotton species ~1-2 million years ago. It has been reported that a great number of genes were quickly erased or preferentially remained after whole-genome duplication, ultimately leading to morphogenesis evolution. However, microRNAs (miRNAs), a new class of gene regulators, have not been well studied in polyploidization. Here, we systematically investigated miRNA evolution amongst cultivated upland cotton G. hirsutum (AADD) and its two ancestors, G. arboreum (AA) and G. raimondii (DD). Our results show that certain highly conserved miRNAs were likely to be lost, whereas certain were remained after genome polyploidization. Cotton-specific miRNAs might undergo remarkably expansion, resulting in overall miRNA increase in upland cotton. Based on the sequenced genomes of G. arboreum and G. raimondii, we are capable for the first time to categorize the origin of miRNAs and coding genes in upland cotton. Different genome-derived miRNAs and miRNA*s displayed asymmetric expression pattern, implicating their diverse functions in upland cotton. No miRNA targeting preference was observed between different genome-derived miRNAs. The origin of miRNAs and coding genes has no impact on becoming miRNAs and their targets, despite some miRNAs and their targets are extremely conserved in the three cotton species. GO- and KEGG-based analysis of conserved miRNAs show that conserved miRNAs and their targets participate in a series of important biological processes and metabolism pathways. Additionally, A-derived miRNAs might be more responsible for ovule and fibre development.

Small RNA sequencing identifies miRNA roles in ovule and fibre development.

MicroRNAs (miRNAs) have been found to be differentially expressed during cotton fibre development. However, which specific miRNAs and how they are involved in fibre development is unclear. Here, using deep sequencing, 65 conserved miRNA families were identified and 32 families were differentially expressed between leaf and ovule. At least 40 miRNAs were either leaf or ovule specific, whereas 62 miRNAs were shared in both leaf and ovule. qRT-PCR confirmed these miRNAs were differentially expressed during fibre early development. A total of 820 genes were potentially targeted by the identified miRNAs, whose functions are involved in a series of biological processes including fibre development, metabolism and signal transduction. Many predicted miRNA-target pairs were subsequently validated by degradome sequencing analysis. GO and KEGG analyses showed that the identified miRNAs and their targets were classified to 1027 GO terms including 568 biological processes, 324 molecular functions and 135 cellular components and were enriched to 78 KEGG pathways. At least seven unique miRNAs participate in trichome regulatory interaction network. Eleven trans-acting siRNA (tasiRNA) candidate genes were also identified in cotton. One has never been found in other plant species and two of them were derived from MYB and ARF, both of which play important roles in cotton fibre development. Sixteen genes were predicted to be tasiRNA targets, including sucrose synthase and MYB2. Together, this study discovered new miRNAs in cotton and offered evidences that miRNAs play important roles in cotton ovule/fibre development. The identification of tasiRNA genes and their targets broadens our understanding of the complicated regulatory mechanism of miRNAs in cotton.

Deep sequencing reveals important roles of microRNAs in response to drought and salinity stress in cotton.

Drought and salinity are two major environmental factors adversely affecting plant growth and productivity. However, the regulatory mechanism is unknown. In this study, the potential roles of small regulatory microRNAs (miRNAs) in cotton response to those stresses were investigated. Using next-generation deep sequencing, a total of 337 miRNAs with precursors were identified, comprising 289 known miRNAs and 48 novel miRNAs. Of these miRNAs, 155 miRNAs were expressed differentially. Target prediction, Gene Ontology (GO)-based functional classification, and Kyoto Encyclopedia of Genes and Genomes (KEGG)-based functional enrichment show that these miRNAs might play roles in response to salinity and drought stresses through targeting a series of stress-related genes. Degradome sequencing analysis showed that at least 55 predicted target genes were further validated to be regulated by 60 miRNAs. CitationRank-based literature mining was employed to determinhe the importance of genes related to drought and salinity stress. The NAC, MYB, and MAPK families were ranked top under the context of drought and salinity, indicating their important roles for the plant to combat drought and salinity stress. According to target prediction, a series of cotton miRNAs are associated with these top-ranked genes, including miR164, miR172, miR396, miR1520, miR6158, ghr-n24, ghr-n56, and ghr-n59. Interestingly, 163 cotton miRNAs were also identified to target 210 genes that are important in fibre development. These results will contribute to cotton stress-resistant breeding as well as understanding fibre development.

Global microRNA modification in cotton (Gossypium hirsutum L.).

MicroRNAs (miRNAs) are small noncoding RNAs participating in versatile biological processes via post-transcriptionally gene regulation. However, how miRNAs are modified or degraded remains unknown, despite years of studies have unravelled much details of miRNA biogenesis and function. Here, we systematically investigated miRNA modification using six small RNA sequencing libraries generated from cotton seedling as well as cotton fibre at five developmental stages. Our results show that 1-2-nt truncation and addition on both 5' and 3' ends of miRNAs are the major modification forms. The 5' and 3' end miRNA modification was almost equal in the six development stages. Truncation was more common than addition on both 5' and 3' end. Structure analysis of the 5' and 3' ends of miRNAs and isomiRs shows that uridine is the preferential nucleotide at the first position of both 5' and 3' ends. According to analysis of nucleotides truncated and tailed from miRNAs, both miRNAs and isomiRs share a similar positional structure distribution at their 5' and 3' ends, respectively. Furthermore, opposite to previous reports, cytodine is more frequently truncated and tailed from the two ends of isomiRs, implying existence of a complex cytodine balance in isomiRs. Comparison of isomiR expression shows differential miRNA modification amongst the six developmental stages in terms of selective modification form, development-dependent modification and differential expression abundance. Our results globally uncovered miRNA modification features in cotton, which could contribute us to understanding miRNA's postmature modification and its regulatory function.

Genome-wide identification and expression analysis of TCP transcription factors in Gossypium raimondii.

Plant-specific TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors play versatile functions in multiple aspects of plant growth and development. However, no systematical study has been performed in cotton. In this study, we performed for the first time the genome-wide identification and expression analysis of the TCP transcription factor family in Gossypium raimondii. A total of 38 non-redundant cotton TCP encoding genes were identified. The TCP transcription factors were divided into eleven subgroups based on phylogenetic analysis. Most TCP genes within the same subfamily demonstrated similar exon and intron organization and the motif structures were highly conserved among the subfamilies. Additionally, the chromosomal distribution pattern revealed that TCP genes were unevenly distributed across 11 out of the 13 chromosomes; segmental duplication is a predominant duplication event for TCP genes and the major contributor to the expansion of TCP gene family in G. raimondii. Moreover, the expression profiles of TCP genes shed light on their functional divergence.

Expression of microRNAs and their targets regulates floral development in tobacco (Nicotiana tabacum).

MicroRNAs (miRNAs) are an extensive class of endogenous posttranscriptional gene regulators that function to mediate gene expression by cleaving target mRNAs or by preventing protein translation. Because of their importance in mediating gene regulation, identifying and elucidating the function of miRNAs have been the primary focus of many researchers. Now that many miRNAs have been identified and assessed for their functionality, the next step is to create expression profiles for miRNAs, so that gene expression studies can be further enhanced with knowledge of the basal expression levels of miRNAs and their targets. In a previous study, 259 putative miRNAs were identified in tobacco, in which 11 of them were confirmed. The primary goal of this study was to further expand on that study and create an expression profile for nine miRNAs and their targets in a tissue-specific manner in tobacco. We chose to study miRNAs that largely play a role in floral development and nutrient stress response. The results of our study show that all tested miRNAs and their targets were expressed in a differential manner. The results of our study also show that out of the tested miRNAs and their targets, miR159, miR157, miR167, miR172, and superoxide dismutase were expressed the highest, suggesting that these genes may play a vital role in the growth and development of tobacco. Corrected expression of miRNAs and their targets regulates floral development.

Genome sequencing and analysis of the paclitaxel-producing endophytic fungus Penicillium aurantiogriseum NRRL 62431.

BACKGROUND: Paclitaxel (Taxol™) is an important anticancer drug with a unique mode of action. The biosynthesis of paclitaxel had been considered restricted to the Taxus species until it was discovered in Taxomyces andreanae, an endophytic fungus of T. brevifolia. Subsequently, paclitaxel was found in hazel (Corylus avellana L.) and in several other endophytic fungi. The distribution of paclitaxel in plants and endophytic fungi and the reported sequence homology of key genes in paclitaxel biosynthesis between plant and fungi species raises the question about whether the origin of this pathway in these two physically associated groups could have been facilitated by horizontal gene transfer. RESULTS: The ability of the endophytic fungus of hazel Penicillium aurantiogriseum NRRL 62431 to independently synthesize paclitaxel was established by liquid chromatography-mass spectrometry and proton nuclear magnetic resonance. The genome of Penicillium aurantiogriseum NRRL 62431 was sequenced and gene candidates that may be involved in paclitaxel biosynthesis were identified by comparison with the 13 known paclitaxel biosynthetic genes in Taxus. We found that paclitaxel biosynthetic gene candidates in P. aurantiogriseum NRRL 62431 have evolved independently and that horizontal gene transfer between this endophytic fungus and its plant host is unlikely. CONCLUSIONS: Our findings shed new light on how paclitaxel-producing endophytic fungi synthesize paclitaxel, and will facilitate metabolic engineering for the industrial production of paclitaxel from fungi.

High-throughput deep sequencing shows that microRNAs play important roles in switchgrass responses to drought and salinity stress.

MicroRNAs (miRNAs) are an important class of small regulatory RNAs. The goal of this study was to analyse stress-responsive miRNAs in switchgrass (Panicum virgatum), the emerging biofuel crop, to facilitate choosing gene targets for improving biomass and biofuel yield. After sequencing three small RNA libraries constructed from control, salt- and drought-treated switchgrass using Illumina sequencing technology, we identified 670 known miRNA families from a total of more than 50 million short reads. A total of 273 miRNAs were identified with precursors: 126 conserved miRNAs and 147 novel miRNAs. Of them, 265 miRNAs were found to have their opposite sequences (miRNA*) with 2-nt overhang on the 3' end. Of them, 194 were detected in switchgrass transcriptome sequences generated from 31 high-throughput RNA sequencing (RNA-Seq) data sets in NCBI. Many miRNAs were differentially or uniquely expressed during salinity or drought stress treatment. We also discovered 11 miRNA clusters containing 29 miRNAs. These identified miRNAs potentially targeted 28549 genes with a various function, including transcription factors, stress-response proteins and cellulose biosynthesis-related proteins. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the identified miRNAs and their targets were classified to 3779 GO terms including 1534 molecular functions, 1851 biological processes and 394 cellular components and were enriched to 147 KEGG pathways. Interestingly, 195 miRNA families and 450 targets were involved in the biosynthesis pathways of carbon, glucose, starch, fatty acid and lignin and in xylem formation, which could aid in designing next-generation switchgrass for biomass and biofuel.

Deep sequencing reveals transcriptome re-programming of Taxus × media cells to the elicitation with methyl jasmonate.

BACKGROUND: Plant cell culture represents an alternative source for producing high-value secondary metabolites including paclitaxel (Taxol®), which is mainly produced in Taxus and has been widely used in cancer chemotherapy. The phytohormone methyl jasmonate (MeJA) can significantly increase the production of paclitaxel, which is induced in plants as a secondary metabolite possibly in defense against herbivores and pathogens. In cell culture, MeJA also elicits the accumulation of paclitaxel; however, the mechanism is still largely unknown. METHODOLOGY/PRINCIPAL FINDINGS: To obtain insight into the global regulation mechanism of MeJA in the steady state of paclitaxel production (7 days after MeJA addition), especially on paclitaxel biosynthesis, we sequenced the transcriptomes of MeJA-treated and untreated Taxus × media cells and obtained ∼ 32.5 M high quality reads, from which 40,348 unique sequences were obtained by de novo assembly. Expression level analysis indicated that a large number of genes were associated with transcriptional regulation, DNA and histone modification, and MeJA signaling network. All the 29 known genes involved in the biosynthesis of terpenoid backbone and paclitaxel were found with 18 genes showing increased transcript abundance following elicitation of MeJA. The significantly up-regulated changes of 9 genes in paclitaxel biosynthesis were validated by qRT-PCR assays. According to the expression changes and the previously proposed enzyme functions, multiple candidates for the unknown steps in paclitaxel biosynthesis were identified. We also found some genes putatively involved in the transport and degradation of paclitaxel. Potential target prediction of miRNAs indicated that miRNAs may play an important role in the gene expression regulation following the elicitation of MeJA. CONCLUSIONS/SIGNIFICANCE: Our results shed new light on the global regulation mechanism by which MeJA regulates the physiology of Taxus cells and is helpful to understand how MeJA elicits other plant species besides Taxus.

miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs.

miRDeepFinder is a software package developed to identify and functionally analyze plant microRNAs (miRNAs) and their targets from small RNA datasets obtained from deep sequencing. The functions available in miRDeepFinder include pre-processing of raw data, identifying conserved miRNAs, mining and classifying novel miRNAs, miRNA expression profiling, predicting miRNA targets, and gene pathway and gene network analysis involving miRNAs. The fundamental design of miRDeepFinder is based on miRNA biogenesis, miRNA-mediated gene regulation and target recognition, such as perfect or near perfect hairpin structures, different read abundances of miRNA and miRNA*, and targeting patterns of plant miRNAs. To test the accuracy and robustness of miRDeepFinder, we analyzed a small RNA deep sequencing dataset of Arabidopsis thaliana published in the GEO database of NCBI. Our test retrieved 128 of 131 (97.7%) known miRNAs that have a more than 3 read count in Arabidopsis. Because many known miRNAs are not associated with miRNA*s in small RNA datasets, miRDeepFinder was also designed to recover miRNA candidates without the presence of miRNA*. To mine as many miRNAs as possible, miRDeepFinder allows users to compare mature miRNAs and their miRNA*s with other small RNA datasets from the same species. Cleaveland software package was also incorporated into miRDeepFinder for miRNA target identification using degradome sequencing analysis. Using this new computational tool, we identified 13 novel miRNA candidates with miRNA*s from Arabidopsis and validated 12 of them experimentally. Interestingly, of the 12 verified novel miRNAs, a miRNA named AC1 spans the exons of two genes (UTG71C4 and UGT71C3). Both the mature AC1 miRNA and its miRNA* were also found in four other small RNA datasets. We also developed a tool, "miRNA primer designer" to design primers for any type of miRNAs. miRDeepFinder provides a powerful tool for analyzing small RNA datasets from all species, with or without the availability of genome information. miRDeepFinder and miRNA primer designer are freely available at http://www.leonxie.com/DeepFinder.php and at http://www.leonxie.com/miRNAprimerDesigner.php , respectively. A program (called RefFinder: http://www.leonxie.com/referencegene.php ) was also developed for assessing the reliable reference genes for gene expression analysis, including miRNAs.