Transcriptome analysis of rice mature root tissue and root tips in early development by massive parallel sequencing

De novo transcriptome sequencing of drought tolerance-associated genes in little millet (Panicum sumatrense L.)

The genome size of the little millet Panicum sumatrense is unknown, although its genome is fairly diploid (2n = 4x = 36). Despite tremendous nutritional value and adaptability to adverse climatic conditions, P. sumatrense use was limited by their low palatability, coarse grain, and lack of variety of culinary preparations. Hence, understanding how to vary their usage to offer food and nutritional security in the continuously changing modern world, the proposed study was aimed to determine potential genes and metabolites implicated in drought resistance. The drought-resistant genotype of tiny millet OLM-203/Tarini was offered in pots under both relaxed and demanding circumstances. The experimental seedlings were 32 days old and had been under water stress for 23 days. A total of 7606 genes were compared between 23 and 32 days for roots and 7264 total genes were compared between 23 and 32 days for leaves, according to a research on differential expression genes (DEGs). Twenty essential genes for drought tolerance were up-or down-regulated in the control and treated roots of the OLM-203 genotype. For instance, the genes RS193 and XB34 were up-regulated in leaves while, WLIM1 was found to be down-regulated. Gene SKI35 was up-regulated in roots, whereas MPK6 and TCMOp1 were down-regulated in root samples. The roots and leaves of the tiny millet OLM-203 genotype expressed 36 up-regulated and 21 down-regulated serine transcripts, respectively. Gene annotations for leaf samples were classified as having "molecular function" (46%), "cellular component" (19%), and "biological process" (35%), while root sample gene annotations were categorized as having "biological process" (573 contigs), "molecular function" (401 contigs), and "cellular components" (166 contigs). Noteworthy, polyamines play a crucial role in drought stress tolerance in the genotype, and it was found that top ten DEGs encoding for polyamines were common in two tissues (leaf and root). Collectively, transcriptomics profiling (RNA-seq) unveiled transcriptional stability drought stress provide a new insight in underlying modus of operandi in little millet genotype "OLM-203/Tarini" in response to heat stress.

Analysis of oxidase activity and transcriptomic changes related to cutting propagation of hybrid larch

Hybrid larch is the main timber and afforestation tree species in Northeast China. To solve the problem of rooting difficulties in larch cutting propagation, enzyme activity determination and transcriptome sequencing were carried out on the rooting tissues at five timepoints after cutting. peroxidase (POD), indole acetic acid oxidase (IAAO) and polyphenol oxidase (PPO) play important roles in the larch rooting process after cutting. A total of 101.20 Gb of clean data was obtained by transcriptome sequencing, and 43,246 unigenes were obtained after further screening and assembly. According to GO analysis and KEGG enrichment analysis, we think that plant hormones play an important role in the rooting process of larch stem cuttings. in the plant hormone signal transduction pathway, a larch gene c141104.graph_c0 that is homologous to the Arabidopsis AUX1 was found to be significantly up-regulated. We suggest that AUX1 may promote IAA transport in larch, thus affecting adventitious root development. According to the results of POD, PPO IAAO indexes and GO analysis, we think s1 and s2 periods may be important periods in the rooting process of larch stem cuttings, so we built a gene regulatory network, a total of 14genes, including LBD, NAC, AP2/ERF, bHLH and etc., may be important in different stages of cutting propagation. As the rooting rate after cutting inhibits the development of larch clone propagation, identifying the genes that regulate rooting could help us to preliminarily understand the molecular mechanism of adventitious root formation and select a better treatment method for cutting propagation.

Genome-Wide Expression Analysis of Root Tips in Contrasting Rice Genotypes Revealed Novel Candidate Genes for Water Stress Adaptation

Root system architecture (RSA) is an important agronomic trait with vital roles in plant productivity under water stress conditions. A deep and branched root system may help plants to avoid water stress by enabling them to acquire more water and nutrient resources. Nevertheless, our knowledge of the genetics and molecular control mechanisms of RSA is still relatively limited. In this study, we analyzed the transcriptome response of root tips to water stress in two well-known genotypes of rice: IR64, a high-yielding lowland genotype, which represents a drought-susceptible and shallow-rooting genotype; and Azucena, a traditional, upland, drought-tolerant and deep-rooting genotype. We collected samples from three zones (Z) of root tip: two consecutive 5 mm sections (Z1 and Z2) and the following next 10 mm section (Z3), which mainly includes meristematic and maturation regions. Our results showed that Z1 of Azucena was enriched for genes involved in cell cycle and division and root growth and development whereas in IR64 root, responses to oxidative stress were strongly enriched. While the expansion of the lateral root system was used as a strategy by both genotypes when facing water shortage, it was more pronounced in Azucena. Our results also suggested that by enhancing meristematic cell wall thickening for insulation purposes as a means of confronting stress, the sensitive IR64 genotype may have reduced its capacity for root elongation to extract water from deeper layers of the soil. Furthermore, several members of gene families such as NAC, AP2/ERF, AUX/IAA, EXPANSIN, WRKY, and MYB emerged as main players in RSA and drought adaptation. We also found that HSP and HSF gene families participated in oxidative stress inhibition in IR64 root tip. Meta-quantitative trait loci (QTL) analysis revealed that 288 differentially expressed genes were colocalized with RSA QTLs previously reported under drought and normal conditions. This finding warrants further research into their possible roles in drought adaptation. Overall, our analyses presented several major molecular differences between Azucena and IR64, which may partly explain their differential root growth responses to water stress. It appears that Azucena avoided water stress through enhancing growth and root exploration to access water, whereas IR64 might mainly rely on cell insulation to maintain water and antioxidant system to withstand stress. We identified a large number of novel RSA and drought associated candidate genes, which should encourage further exploration of their potential to enhance drought adaptation in rice.

Toward Integrated Multi-Omics Intervention: Rice Trait Improvement and Stress Management

Rice (Oryza sativa) is an imperative staple crop for nearly half of the world's population. Challenging environmental conditions encompassing abiotic and biotic stresses negatively impact the quality and yield of rice. To assure food supply for the unprecedented ever-growing world population, the improvement of rice as a crop is of utmost importance. In this era, "omics" techniques have been comprehensively utilized to decipher the regulatory mechanisms and cellular intricacies in rice. Advancements in omics technologies have provided a strong platform for the reliable exploration of genetic resources involved in rice trait development. Omics disciplines like genomics, transcriptomics, proteomics, and metabolomics have significantly contributed toward the achievement of desired improvements in rice under optimal and stressful environments. The present review recapitulates the basic and applied multi-omics technologies in providing new orchestration toward the improvement of rice desirable traits. The article also provides a catalog of current scenario of omics applications in comprehending this imperative crop in relation to yield enhancement and various environmental stresses. Further, the appropriate databases in the field of data science to analyze big data, and retrieve relevant information vis-à-vis rice trait improvement and stress management are described.

Integrated Analysis of mRNA and Non-coding RNA Transcriptome in Pepper (Capsicum chinense) Hybrid at Seedling and Flowering Stages

Pepper is an important vegetable in the world. In this work, mRNA and ncRNA transcriptome profiles were applied to understand the heterosis effect on the alteration in the gene expression at the seedling and flowering stages between the hybrid and its parents in Capsicum chinense. Our phenotypic data indicated that the hybrid has dominance in leaf area, plant scope, plant height, and fruit-related traits. Kyoto Encyclopedia of Genes and Genomes analysis showed that nine members of the plant hormone signal transduction pathway were upregulated in the seedling and flowering stages of the hybrid, which was supported by weighted gene coexpression network analysis and that BC332_23046 (auxin response factor 8), BC332_18317 (auxin-responsive protein IAA20), BC332_13398 (ethylene-responsive transcription factor), and BC332_27606 (ethylene-responsive transcription factor WIN1) were candidate hub genes, suggesting the important potential role of the plant hormone signal transduction in pepper heterosis. Furthermore, some transcription factor families, including bHLH, MYB, and HSF were greatly over-dominant. We also identified 2,525 long ncRNAs (lncRNAs), 47 micro RNAs (miRNAs), and 71 circle RNAs (circRNAs) in the hybrid. In particular, downregulation of miR156, miR169, and miR369 in the hybrid suggested their relationship with pepper growth vigor. Moreover, we constructed some lncRNA–miRNA–mRNA regulatory networks that showed a multi-dimension to understand the ncRNA relationship with heterosis. These results will provide guidance for a better understanding of the molecular mechanism involved in pepper heterosis.

Dynamic formation and transcriptional regulation mediated by phytohormones during chalkiness formation in rice

BACKGROUND

Rice (Oryza sativa L.) Chalkiness, the opaque part in the kernel endosperm formed by loosely piled starch and protein bodies. Chalkiness is a complex quantitative trait regulated by multiple genes and various environmental factors. Phytohormones play important roles in the regulation of chalkiness formation but the underlying molecular mechanism is still unclear at present.

RESULTS

In this research, Xiangzaoxian24 (X24, pure line of indica rice with high-chalkiness) and its origin parents Xiangzaoxian11 (X11, female parent, pure line of indica rice with high-chalkiness) and Xiangzaoxian7 (X7, male parent, pure line of indica rice with low-chalkiness) were used as materials. The phenotype, physiological and biochemical traits combined with transcriptome analysis were conducted to illustrate the dynamic process and transcriptional regulation of rice chalkiness formation. Impressively, phytohormonal contents and multiple phytohormonal signals were significantly different in chalky caryopsis, suggesting the involvement of phytohormones, particularly ABA and auxin, in the regulation of rice chalkiness formation, through the interaction of multiple transcription factors and their downstream regulators.

CONCLUSION

These results indicated that chalkiness formation is a dynamic process associated with multiple genes, forming a complex regulatory network in which phytohormones play important roles. These results provided informative clues for illustrating the regulatory mechanisms of chalkiness formation in rice.

Plasmodesmata play pivotal role in sucrose supply to Meloidogyne graminicola-caused giant cells in rice

On infection, plant-parasitic nematodes establish feeding sites in roots from which they take up carbohydrates among other nutrients. Knowledge on how carbohydrates are supplied to the nematodes' feeding sites is limited. Here, gene expression analyses showed that RNA levels of OsSWEET11 to OsSWEET15 were extremely low in both Meloidogyne graminicola (Mg)-caused galls and noninoculated roots. All the rice sucrose transporter genes, OsSUT1 to OsSUT5, were either down-regulated in Mg-caused galls compared with noninoculated rice roots or had very low transcript abundance. OsSUT1 was the only gene up-regulated in galls, at 14 days postinoculation (dpi), after being highly down-regulated at 3 and 7 dpi. OsSUT4 was down-regulated at 3 dpi. No noticeable OsSUTs promoter activities were detected in Mg-caused galls of pOsSUT1 to -5::GUS rice lines. Loading experiments with carboxyfluorescein diacetate (CFDA) demonstrated that symplastic connections exist between phloem and Mg-caused giant cells (GCs). According to data from OsGNS5- and OsGSL2-overexpressing rice plants that had decreased and increased callose deposition, respectively, callose negatively affected Mg parasitism and sucrose supply to Mg-caused GCs. Our results suggest that plasmodesmata-mediated sucrose transport plays a pivotal role in sucrose supply from rice root phloem to Mg-caused GCs, and OsSWEET11 to -15 and OsSUTs are not major players in it, although further functional analysis is needed for OsSUT1 and OsSUT4.

Transcriptome analysis of biological pathways associated with heterosis in Chinese cabbage

Chinese cabbage is an important vegetable in Asia, and high-yielding hybrids are needed to cope with the growing demand. A comparative transcriptome profiling was conducted to reveal the differentially expressed genes (DEGs) associated with heterosis in two hybrids relative to their parents. Our data suggests that heterosis is underlined by a significant upregulation of gene expression. High expression of DEGs in glycolysis and photosynthesis pathways in hybrids depicted their relation with growth and hybrid vigor. Besides, DEGs related to auxin, abscisic acid, ethylene and gibberellin were identified, implying that these hormones may boost the mechanisms of growth and developmental processes in the hybrids. Furthermore, transcription factors, including bHLH, ERF, MYB and WRKY were predicted to regulate downstream genes linked to hybrid vigor. Collectively, the present study will be helpful for a better understanding of the regulation mechanisms of heterosis to aid cabbage yield improvement.

Transcriptome and GWAS analyses reveal candidate gene for seminal root length of maize seedlings under drought stress

Water deficits are a major constraint on maize growth and yield, and deep roots are one of the major mechanisms of drought tolerance. In this study, four root and shoot traits were evaluated within an association panel consisting of 209 diverse maize accessions under well-watered (WW) and water-stressed (WS) conditions. A significant positive correlation was observed between seminal root length (SRL) under WS treatment and the drought tolerance index (DI) of maize seedlings. The transcriptome profiles of maize seminal roots were compared between four drought-tolerant lines and four drought-sensitive lines under both water conditions to identify genes associated with the drought stress response. After drought stress, 343 and 177 common differentially expressed genes (DEGs) were identified in the drought-tolerant group and drought-sensitive group, respectively. In parallel, a coexpression network underlying SRL was constructed on the basis of transcriptome data, and 10 hub genes involved in two significant associated modules were identified. Additionally, a genome-wide association study (GWAS) of the SRL revealed 62 loci for the two water treatments. By integrating the results of the GWAS, the common DEGs and the coexpression network analysis, 7 promising candidate genes were prioritized for further research. Together, our results provide a foundation for the enhanced understanding of seminal root changes in response to drought stress in maize.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Comparing RNA-Seq and microarray gene expression data in two zones of the Arabidopsis root apex relevant to spaceflight

PREMISE OF THE STUDY

The root apex is an important region involved in environmental sensing, but comprises a very small part of the root. Obtaining root apex transcriptomes is therefore challenging when the samples are limited. The feasibility of using tiny root sections for transcriptome analysis was examined, comparing RNA sequencing (RNA-Seq) to microarrays in characterizing genes that are relevant to spaceflight.

METHODS

Arabidopsis thaliana Columbia ecotype (Col-0) roots were sectioned into Zone 1 (0.5 mm; root cap and meristematic zone) and Zone 2 (1.5 mm; transition, elongation, and growth-terminating zone). Differential gene expression in each was compared.

RESULTS

Both microarrays and RNA-Seq proved applicable to the small samples. A total of 4180 genes were differentially expressed (with fold changes of 2 or greater) between Zone 1 and Zone 2. In addition, 771 unique genes and 19 novel transcriptionally active regions were identified by RNA-Seq that were not detected in microarrays. However, microarrays detected spaceflight-relevant genes that were missed in RNA-Seq.

DISCUSSION

Single root tip subsections can be used for transcriptome analysis using either RNA-Seq or microarrays. Both RNA-Seq and microarrays provided novel information. These data suggest that techniques for dealing with small, rare samples from spaceflight can be further enhanced, and that RNA-Seq may miss some spaceflight-relevant changes in gene expression.

Holm oak proteomic response to water limitation at seedling establishment stage reveals specific changes in different plant parts as well as interaction between roots and cotyledons

Quercus ilex is a dominant tree species in the Mediterranean region with double economic and ecological importance and increasing use in reforestation. Seedling establishment is extremely vulnerable to environmental stresses, particularly drought. A time course study on physiological and proteomic response of holm oak to water limitation stress and recovery during early heterotrophic growth is reported. Applied stress led to diminution in plant water content and root growth, oxidative stress in roots and some alterations in the anti-oxidative protection. Plant parts differed substantially in soluble sugar and free phenolic content, and in their changes during stress and recovery. Proteomic response in holm oak roots and cotyledons was estimated using combined 1-DE/2-DE approach and protein identification by MALDI TOF-TOF PMF and MS/MS. A total of 127 differentially abundant protein species (DAPs) were identified. DAPs related to starch metabolism, lipid to sugar conversion, reserve proteins and their mobilization were typical for cotyledons. DAPs in roots were involved in sugar utilization, secondary metabolism and defense, including pathogenesis related proteins from PR-5 and PR-10 families. Results emphasize specific proteome signatures of separate plant parts as well as importance of sink-source interaction between root and cotyledon in the time course of stress and in recovery.

Differential transcriptome patterns associated with early seedling development in a wild and a domesticated common bean (Phaseolus vulgaris L.) accession

Genes that control "Domestication Syndrome" traits were direct targets of selection, like those controlling increased seed size in the common bean. However, selection for this trait brought about unintentional selection on genes controlling seedling growth. We hypothesized that wild and domesticated plants have different early seedling growth patterns as an indirect consequence of selection for a larger seed size during domestication, and those differences resulted from changes in gene expression patterns of the wild ancestor. Large seeds pose a challenge to reserve remobilization during early heterotrophic growth, particularly during a transition towards more fertile alluvial soils. To address our hypothesis, we characterized the patterns of gene expression of cotyledon, root, and leaf tissues of 7-day old seedlings of a wild and a landrace accession of the common bean. Differential expression analyses detected genes with contrasting patterns of expression between the two genotypes in all three tissues. Some of the differentially expressed genes with contrasting genotypic patterns are known to have domestication-related signatures of selection. Among these genes were some transcription factors associated with key roles in development. These genes may represent targets of indirect selection and ultimately explain the growth phenotypic differences between wild and domesticated seedlings.

Analysis of Genes with Alternatively Spliced Transcripts in the Leaf, Root, Panicle and Seed of Rice Using a Long Oligomer Microarray and RNA-Seq

Pre-mRNA splicing further increases protein diversity acquired through evolution. The underlying driving forces for this phenomenon are unknown, especially in terms of gene expression. A rice alternatively spliced transcript detection microarray (ASDM) and RNA sequencing (RNA-Seq) were applied to differentiate the transcriptome of 4 representative organs of Oryza sativa L. cv. Ilmi: leaves, roots, 1-cm-stage panicles and young seeds at 21 days after pollination. Comparison of data obtained by microarray and RNA-Seq showed a bell-shaped distribution and a co-lineation for highly expressed genes. Transcripts were classified according to the degree of organ enrichment using a coefficient value (CV, the ratio of the standard deviation to the mean values): highly variable (CVI), variable (CVII), and constitutive (CVIII) groups. A higher index of the portion of loci with alternatively spliced transcripts in a group (IAST) value was observed for the constitutive group. Genes of the highly variable group showed the characteristics of the examined organs, and alternatively spliced transcripts tended to exhibit the same organ specificity or less organ preferences, with avoidance of 'organ distinctness'. In addition, within a locus, a tendency of higher expression was found for transcripts with a longer coding sequence (CDS), and a spliced intron was the most commonly found type of alternative splicing for an extended CDS. Thus, pre-mRNA splicing might have evolved to retain maximum functionality in terms of organ preference and multiplicity.

Transcriptome-wide identification and characterization of CAD isoforms specific for podophyllotoxin biosynthesis from Podophyllum hexandrum

Podophyllotoxin (ptox) is a therapeutically important lignan derived from Podophyllum hexandrum and is used as a precursor for the synthesis of anticancer drugs etoposide, teniposide and etopophose. In spite of its enormous economic significance, genomic information on this endangered medicinal herb is scarce. We have performed de novo transcriptome analysis of methyl jasmonate (MeJA)-treated P. hexandrum cell cultures exhibiting enhanced ptox accumulation. The results revealed the maximum up-regulation of several isoforms of cinnamyl alcohol dehydrogenase (CAD). CAD catalyzes the synthesis of coniferyl alcohol and sinapyl alcohol from coniferaldehyde (CAld) and sinapaldehyde respectively. Coniferyl alcohol can produce both lignin and lignan while sinapyl alcohol produces only lignin. To isolate the CAD isoforms favoring ptox, we deduced full length cDNA sequences of four CAD isoforms: PhCAD1, PhCAD2, PhCAD3 and PhCAD4 from the contigs of the transcriptome data. In vitro enzyme assays indicated a higher affinity for CAld over sinapaldehyde for each isoform. In silico molecular docking analyses also suggested that PhCAD3 has a higher binding preference with CAld over sinapaldehyde, followed by PhCAD4, PhCAD2, and PhCAD1, respectively. The transgenic cell cultures overexpressing these isoforms independently revealed that PhCAD3 favored the maximum accumulation of ptox as compared to lignin followed by PhCAD4 and PhCAD2, whereas, PhCAD1 favored both equally. Together, our study reveals transcriptome-wide identification and characterization of ptox specific CAD isoforms from P. hexandrum. It provides a useful resource for future research not only on the ptox biosynthetic pathway but on overall P. hexandrum, an endangered medicinal herb with immense therapeutic importance.

High-Throughput Sequencing Reveals Single Nucleotide Variants in Longer-Kernel Bread Wheat

The transcriptomes of bread wheat Yunong 201 and its ethyl methanesulfonate derivative Yunong 3114 were obtained by next-sequencing technology. Single nucleotide variants (SNVs) in the wheat strains were explored and compared. A total of 5907 and 6287 non-synonymous SNVs were acquired for Yunong 201 and 3114, respectively. A total of 4021 genes with SNVs were obtained. The genes that underwent non-synonymous SNVs were significantly involved in ATP binding, protein phosphorylation, and cellular protein metabolic process. The heat map analysis also indicated that most of these mutant genes were significantly differentially expressed at different developmental stages. The SNVs in these genes possibly contribute to the longer kernel length of Yunong 3114. Our data provide useful information on wheat transcriptome for future studies on wheat functional genomics. This study could also help in illustrating the gene functions of the non-synonymous SNVs of Yunong 201 and 3114.

Differential transcriptome analysis between Populus and its synthesized allotriploids driven by second-division restitution

In this report, we compared transcriptomic differences between a synthetic Populus section Tacamahaca triploid driven by second-division restitution and its parents using a high-throughput RNA-seq method. A total of 4,080 genes were differentially expressed between the high-growth vigor allotriploids (SDR-H) and their parents, and 719 genes were non-additively expressed in SDR-H. Differences in gene expression between the allotriploid and male parent were more significant than those between the allotriploid and female parent, which may be caused by maternal effects. We observed 3,559 differentially expressed genes (DEGs) between the SDR-H and male parent. Notably, the genes were mainly involved in metabolic process, cell proliferation, DNA methylation, cell division, and meristem and developmental growth. Among the 1,056 DEGs between SDR-H and female parent, many genes were associated with metabolic process and carbon utilization. In addition, 1,789 DEGs between high- and low-growth vigor allotriploid were mainly associated with metabolic process, auxin poplar transport, and regulation of meristem growth. Our results indicated that the higher poplar ploidy level can generate extensive transcriptomic diversity compared with its parents. Overall, these results increased our understanding of the driving force for phenotypic variation and adaptation in allopolyploids driven by second-division restitution.

Wild soybean roots depend on specific transcription factors and oxidation reduction related genesin response to alkaline stress

Soil alkalinity is an important environmental problem limiting agricultural productivity. Wild soybean (Glycine soja) shows strong alkaline stress tolerance, so it is an ideal plant candidate for studying the molecular mechanisms of alkaline tolerance and identifying alkaline stress-responsive genes. However, limited information is available about G. soja responses to alkaline stress on a genomic scale. Therefore, in the present study, we used RNA sequencing to compare transcript profiles of G. soja root responses to sodium bicarbonate (NaHCO3) at six time points, and a total of 68,138,478 pairs of clean reads were obtained using the Illumina GAIIX. Expression patterns of 46,404 G. soja genes were profiled in all six samples based on RNA-seq data using Cufflinks software. Then, t12 transcription factors from MYB, WRKY, NAC, bZIP, C2H2, HB, and TIFY families and 12 oxidation reduction related genes were chosen and verified to be induced in response to alkaline stress by using quantitative real-time polymerase chain reaction (qRT-PCR). The GO functional annotation analysis showed that besides "transcriptional regulation" and "oxidation reduction," these genes were involved in a variety of processes, such as "binding" and "response to stress." This is the first comprehensive transcriptome profiling analysis of wild soybean root under alkaline stress by RNA sequencing. Our results highlight changes in the gene expression patterns and identify a set of genes induced by NaHCO3 stress. These findings provide a base for the global analyses of G. soja alkaline stress tolerance mechanisms.

The role of thionins in rice defence against root pathogens

Thionins are antimicrobial peptides that are involved in plant defence. Here, we present an in-depth analysis of the role of rice thionin genes in defence responses against two root pathogens: the root-knot nematode Meloidogyne graminicola and the oomycete Pythium graminicola. The expression of rice thionin genes was observed to be differentially regulated by defence-related hormones, whereas all analysed genes were consistently down-regulated in M. graminicola-induced galls, at least until 7 days post-inoculation (dpi). Transgenic lines of Oryza sativa cv. Nipponbare overproducing OsTHI7 revealed decreased susceptibility to M. graminicola infection and P. graminicola colonization. Taken together, these results demonstrate the role of rice thionin genes in defence against two of the most damaging root pathogens attacking rice.

Dynamic Analysis of Gene Expression in Rice Superior and Inferior Grains by RNA-Seq

Poor grain filling of inferior grains located on lower secondary panicle branch causes great drop in rice yield and quality. Dynamic gene expression patterns between superior and inferior grains were examined from the view of the whole transcriptome by using RNA-Seq method. In total, 19,442 genes were detected during rice grain development. Genes involved in starch synthesis, grain storage and grain development were interrogated in particular in superior and inferior grains. Of the genes involved in sucrose to starch transformation process, most were expressed at lower level in inferior grains at early filling stage compared to that of superior grains. But at late filling stage, the expression of those genes was higher in inferior grains and lower in superior grains. The same trends were observed in the expression of grain storage protein genes. While, evidence that genes involved in cell cycle showed higher expression in inferior grains during whole period of grain filling indicated that cell proliferation was active till the late filling stage. In conclusion, delayed expression of most starch synthesis genes in inferior grains and low capacity of sink organ might be two important factors causing low filling rate of inferior grain at early filling stage, and shortage of carbohydrate supply was a limiting factor at late filling stage.

Cytological and transcriptional dynamics analysis of host plant revealed stage-specific biological processes related to compatible rice-Ustilaginoidea virens interaction

Rice false smut, a fungal disease caused by Ustilaginoidea virens is becoming a severe detriment to rice production worldwide. However, little is known about the molecular response of rice to attacks by the smut pathogen. In this article, we define the initial infection process as having three stages: initial colonization on the pistil (stage 1, S1), amplification on the anther (stage 2, S2) and sporulation in the anther chambers (stage 3, S3). Based on the transcriptome of rice hosts in response to U. virens in two separate years, we identified 126, 204, and 580 specific regulated genes in their respective stages S1, S2, and S3, respectively, by excluding common expression patterns in other openly biotic/abiotic databases using bioinformatics. As the disease progresses, several stage-specific biological processes (BP) terms were distinctively enriched: "Phosphorylation" in stage S1, "PCD" in S2, and "Cell wall biogenesis" in S3, implying a concise signal cascade indicative of the tactics that smut pathogens use to control host rice cells during infection. 113 regulated genes were coexpressed among the three stages. They shared highly conserved promoter cis-element in the promoters in response to the regulation of WRKY and Myb for up-regulation, and ABA and Ca2+ for down regulation, indicating their potentially critical roles in signal transduction during rice-U. virens interaction. We further analyzed seven highly regulated unique genes; four were specific to pollen development, implying that pollen-related genes play critical roles in the establishment of rice susceptibility to U. virens. To my knowledge, this is the first report about probing of molecular response of rice to smut pathogen infection, which will greatly expand our understanding of the molecular events surrounding infection by rice false smut.

Transcriptome characterization and sequencing-based identification of drought-responsive genes in potato

Potato (Solanum tubersosum L.) is relatively vulnerable to abiotic stress conditions such as drought, but the tolerance mechanisms to such stress in potato are largely unknown. To gain insight into the transcriptome dynamics that are associated with drought stress, genome-wide gene expression profile was conducted by Solexa sequencing to generate a large dataset and a comprehensive transcriptome profile for potato. Here, we report a reference for the potato transcriptome using leaf tissues under drought-stressed condition from a local potato cultivar 'Longshu 3'. Analysis of 86,965,482 RNA-Seq reads permitted the detection and quantification of expression levels of 7,284 genes at transcriptional levels, among them, 6,754 genes were enriched in draught-treated leaves while 6,419 in control. We identified 842 drought-responsive up-regulated and 494 down-regulated candidate genes with significantly differentially expression under continued drought stress treatments. Those differently expressed genes were mostly enriched in 89 gene categories and 21 KEGG pathways. Drought-stressed leaves had increased expression of genes involved in stress response compared with control leaves. A subset of differentially expressed genes associated with drought response was examined using quantitative real-time PCR. These results provide a broad spectrum of candidate genes that are essential for understanding the molecular regulation of potato in response to abiotic stresses.

ABC transporters in Dictyostelium discoideum development

ATP-binding cassette (ABC) transporters can translocate a broad spectrum of molecules across the cell membrane including physiological cargo and toxins. ABC transporters are known for the role they play in resistance towards anticancer agents in chemotherapy of cancer patients. There are 68 ABC transporters annotated in the genome of the social amoeba Dictyostelium discoideum. We have characterized more than half of these ABC transporters through a systematic study of mutations in their genes. We have analyzed morphological and transcriptional phenotypes for these mutants during growth and development and found that most of the mutants exhibited rather subtle phenotypes. A few of the genes may share physiological functions, as reflected in their transcriptional phenotypes. Since most of the abc-transporter mutants showed subtle morphological phenotypes, we utilized these transcriptional phenotypes to identify genes that are important for development by looking for transcripts whose abundance was unperturbed in most of the mutants. We found a set of 668 genes that includes many validated D. discoideum developmental genes. We have also found that abcG6 and abcG18 may have potential roles in intercellular signaling during terminal differentiation of spores and stalks.

Transcriptional profiling of sweetpotato (Ipomoea batatas) roots indicates down-regulation of lignin biosynthesis and up-regulation of starch biosynthesis at an early stage of storage root formation

Background

The number of fibrous roots that develop into storage roots determines sweetpotato yield. The aim of the present study was to identify the molecular mechanisms involved in the initiation of storage root formation, by performing a detailed transcriptomic analysis of initiating storage roots using next-generation sequencing platforms. A two-step approach was undertaken: (1) generating a database for the sweetpotato root transcriptome using 454-Roche sequencing of a cDNA library created from pooled samples of two root types: fibrous and initiating storage roots; (2) comparing the expression profiles of initiating storage roots and fibrous roots, using the Illumina Genome Analyzer to sequence cDNA libraries of the two root types and map the data onto the root transcriptome database.

Results

Use of the 454-Roche platform generated a total of 524,607 reads, 85.6% of which were clustered into 55,296 contigs that matched 40,278 known genes. The reads, generated by the Illumina Genome Analyzer, were found to map to 31,284 contigs out of the 55,296 contigs serving as the database. A total of 8,353 contigs were found to exhibit differential expression between the two root types (at least 2.5-fold change). The Illumina-based differential expression results were validated for nine putative genes using quantitative real-time PCR. The differential expression profiles indicated down-regulation of classical root functions, such as transport, as well as down-regulation of lignin biosynthesis in initiating storage roots, and up-regulation of carbohydrate metabolism and starch biosynthesis. In addition, data indicated delicate control of regulators of meristematic tissue identity and maintenance, associated with the initiation of storage root formation.

Conclusions

This study adds a valuable resource of sweetpotato root transcript sequences to available data, facilitating the identification of genes of interest. This resource enabled us to identify genes that are involved in the earliest stage of storage root formation, highlighting the reduction in carbon flow toward phenylpropanoid biosynthesis and its delivery into carbohydrate metabolism and starch biosynthesis, as major events involved in storage root initiation. The novel transcripts related to storage root initiation identified in this study provide a starting point for further investigation into the molecular mechanisms underlying this process.

microRNA profiling of root tissues and root forming explant cultures in Medicago truncatula

Plant root architecture is regulated by the initiation and modulation of cell division in regions containing pluripotent stem cells known as meristems. In roots, meristems are formed early in embryogenesis, in the case of the root apical meristem (RAM), and during organogenesis at the site of lateral root or, in legumes, nodule formation. Root meristems can also be generated in vitro from leaf explants cultures supplemented with auxin. microRNAs (miRNAs) have emerged as regulators of many key biological functions in plants including root development. To identify key miRNAs involved in root meristem formation in Medicago truncatula, we used deep sequencing to compare miRNA populations. Comparisons were made between: (1) the root tip (RT), containing the RAM and the elongation zone (EZ) tissue and (2) root forming callus (RFC) and non-root forming callus (NRFC). We identified 83 previously reported miRNAs, 24 new to M. truncatula, in 44 families. For the first time in M. truncatula, members of conserved miRNA families miR165, miR181 and miR397 were found. Bioinformatic analysis identified 38 potential novel miRNAs. Selected miRNAs and targets were validated using Taqman miRNA assays and 5' RACE. Many miRNAs were differentially expressed between tissues, particularly RFC and NRFC. Target prediction revealed a number of miRNAs to target genes previously shown to be differentially expressed between RT and EZ or RFC and NRFC and important in root development. Additionally, we predict the miRNA/target relationships for miR397 and miR160 to be conserved in M. truncatula. Amongst the predictions, were AUXIN RESPONSE FACTOR 10, targeted by miR160 and a LACCASE-like gene, targeted by miR397, both are miRNA/target pairings conserved in other species.

Global transcriptome profiling of Salicornia europaea L. shoots under NaCl treatment

BACKGROUND

Soil salinity is a major abiotic stress that limits agriculture productivity worldwide. Salicornia europaea is well adapted to extreme saline environments with more than 1,000 mM NaCl in the soil, so it could serve as an important model species for studying halophilic mechanisms in euhalophytes. To obtain insights into the molecular basis of salt tolerance, we present here the first extensive transcriptome analysis of this species using the Illumina HiSeq™ 2000.

PRINCIPAL FINDINGS

A total of 41 and 39 million clean reads from the salt-treated (Se200S) and salt-free (SeCKS) tissues of S. europaea shoots were obtained, and de novo assembly produced 97,865 and 101,751 unigenes, respectively. Upon further assembly with EST data from both Se200S and SeCKS, 109,712 high-quality non-redundant unigenes were generated with a mean unigene size of 639 bp. Additionally, a total of 3,979 differentially expressed genes (DEGs) were detected between the Se200S and SeCKS libraries, with 348 unigenes solely expressed in Se200S and 460 unigenes solely expressed in SeCKS. Furthermore, we identified a large number of genes that are involved in ion homeostasis and osmotic adjustment, including cation transporters and proteins for the synthesis of low-molecular compounds. All unigenes were functionally annotated within the COG, GO and KEGG pathways, and 10 genes were validated by qRT-PCR.

CONCLUSION

Our data contains the extensive sequencing and gene-annotation analysis of S. europaea. This genetic knowledge will be very useful for future studies on the molecular adaptation to abiotic stress in euhalophytes and will facilitate the genetic manipulation of other economically important crops.

Transcriptome analysis of rice root heterosis by RNA-Seq

Background

Heterosis is a phenomenon in which hybrids exhibit superior performance relative to parental phenotypes. In addition to the heterosis of above-ground agronomic traits on which most existing studies have focused, root heterosis is also an indispensable component of heterosis in the entire plant and of major importance to plant breeding. Consequently, systematic investigations of root heterosis, particularly in reproductive-stage rice, are needed. The recent advent of RNA sequencing technology (RNA-Seq) provides an opportunity to conduct in-depth transcript profiling for heterosis studies.

Results

Using the Illumina HiSeq 2000 platform, the root transcriptomes of the super-hybrid rice variety Xieyou 9308 and its parents were analyzed at tillering and heading stages. Approximately 391 million high-quality paired-end reads (100-bp in size) were generated and aligned against the Nipponbare reference genome. We found that 38,872 of 42,081 (92.4%) annotated transcripts were represented by at least one sequence read. A total of 829 and 4186 transcripts that were differentially expressed between the hybrid and its parents (DG_HP) were identified at tillering and heading stages, respectively. Out of the DG_HP, 66.59% were down-regulated at the tillering stage and 64.41% were up-regulated at the heading stage. At the heading stage, the DG_HP were significantly enriched in pathways related to processes such as carbohydrate metabolism and plant hormone signal transduction, with most of the key genes that are involved in the two pathways being up-regulated in the hybrid. Several significant DG_HP that could be mapped to quantitative trait loci (QTLs) for yield and root traits are also involved in carbohydrate metabolism and plant hormone signal transduction pathways.

Conclusions

An extensive transcriptome dataset was obtained by RNA-Seq, giving a comprehensive overview of the root transcriptomes at tillering and heading stages in a heterotic rice cross and providing a useful resource for the rice research community. Using comparative transcriptome analysis, we detected DG_HP and identified a group of potential candidate transcripts. The changes in the expression of the candidate transcripts may lay a foundation for future studies on molecular mechanisms underlying root heterosis.

The Rice Genome Knowledgebase (RGKbase): an annotation database for rice comparative genomics and evolutionary biology

Over the past 10 years, genomes of cultivated rice cultivars and their wild counterparts have been sequenced although most efforts are focused on genome assembly and annotation of two major cultivated rice (Oryza sativa L.) subspecies, 93-11 (indica) and Nipponbare (japonica). To integrate information from genome assemblies and annotations for better analysis and application, we now introduce a comparative rice genome database, the Rice Genome Knowledgebase (RGKbase, http://rgkbase.big.ac.cn/RGKbase/). RGKbase is built to have three major components: (i) integrated data curation for rice genomics and molecular biology, which includes genome sequence assemblies, transcriptomic and epigenomic data, genetic variations, quantitative trait loci (QTLs) and the relevant literature; (ii) User-friendly viewers, such as Gbrowse, GeneBrowse and Circos, for genome annotations and evolutionary dynamics and (iii) Bioinformatic tools for compositional and synteny analyses, gene family classifications, gene ontology terms and pathways and gene co-expression networks. RGKbase current includes data from five rice cultivars and species: Nipponbare (japonica), 93-11 (indica), PA64s (indica), the African rice (Oryza glaberrima) and a wild rice species (Oryza brachyantha). We are also constantly introducing new datasets from variety of public efforts, such as two recent releases—sequence data from ∼1000 rice varieties, which are mapped into the reference genome, yielding ample high-quality single-nucleotide polymorphisms and insertions–deletions.

Transcriptional reprogramming by root knot and migratory nematode infection in rice

Rice is one of the most important staple crops worldwide, but its yield is compromised by different pathogens, including plant-parasitic nematodes. In this study we have characterized specific and general responses of rice (Oryza sativa) roots challenged with two endoparasitic nematodes with very different modes of action. Local transcriptional changes in rice roots upon root knot (Meloidogyne graminicola) and root rot nematode (RRN, Hirschmanniella oryzae) infection were studied at two time points (3 and 7 d after infection, dai), using mRNA-seq. Our results confirm that root knot nematodes (RKNs), which feed as sedentary endoparasites, stimulate metabolic pathways in the root, and enhance nutrient transport towards the induced root gall. The migratory RRNs, on the other hand, induce programmed cell death and oxidative stress, and obstruct the normal metabolic activity of the root. While RRN infection causes up-regulation of biotic stress-related genes early in the infection, the sedentary RKNs suppress the local defense pathways (e.g. salicylic acid and ethylene pathways). Interestingly, hormone pathways mainly involved in plant development were strongly induced (gibberellin) or repressed (cytokinin) at 3 dai. These results uncover previously unrecognized nematode-induced expression profiles related to their specific infection strategy.

AP2/ERF Transcription Factor in Rice: Genome-Wide Canvas and Syntenic Relationships between Monocots and Eudicots

The transcription factor family intimately regulates gene expression in response to hormones, biotic and abiotic factors, symbiotic interactions, cell differentiation, and stress signalling pathways in plants. In this study, 170 AP2/ERF family genes are identified by phylogenetic analysis of the rice genome (Oryza sativa l. japonica) and they are divided into a total of 11 groups, including four major groups (AP2, ERF, DREB, and RAV), 10 subgroups, and two soloists. Gene structure analysis revealed that, at position-6, the amino acid threonine (Thr-6) is conserved in the double domain AP2 proteins compared to the amino acid arginine (Arg-6), which is preserved in the single domain of ERF proteins. In addition, the histidine (His) amino acid is found in both domains of the double domain AP2 protein, which is missing in single domain ERF proteins. Motif analysis indicates that most of the conserved motifs, apart from the AP2/ERF domain, are exclusively distributed among the specific clades in the phylogenetic tree and regulate plausible functions. Expression analysis reveals a widespread distribution of the rice AP2/ERF family genes within plant tissues. In the vegetative organs, the transcripts of these genes are found most abundant in the roots followed by the leaf and stem; whereas, in reproductive tissues, the gene expression of this family is observed high in the embryo and lemma. From chromosomal localization, it appears that repetition and tandem-duplication may contribute to the evolution of new genes in the rice genome. In this study, interspecies comparisons between rice and wheat reveal 34 rice loci and unveil the extent of collinearity between the two genomes. It was subsequently ascertained that chromosome-9 has more orthologous loci for CRT/DRE genes whereas chromosome-2 exhibits orthologs for ERF subfamily members. Maximum conserved synteny is found in chromosome-3 for AP2 double domain subfamily genes. Macrosynteny between rice and Arabidopsis, a distant, related genome, uncovered 11 homologs/orthologs loci in both genomes. The distribution of AP2/ERF family gene paralogs in Arabidopsis was most frequent in chromosome-1 followed by chromosome-5. In Arabidopsis, ERF subfamily gene orthologs are found on chromosome-1, chromosome-3, and chromosome-5, whereas DRE subfamily genes are found on chromosome-2 and chromosome-5. Orthologs for RAV and AP2 with double domains in Arabidopsis are located on chromosome-1 and chromosome-3, respectively. In conclusion, the data generated in this survey will be useful for conducting genomic research to determine the precise role of the AP2/ERF gene during stress responses with the ultimate goal of improving crops.