Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding

Unraveling the genetic and molecular basis of heat stress in cotton

Human activities and climate change have resulted in frequent and intense weather fluctuations, leading to diverse abiotic stresses on crops which hampers greatly their metabolic activities. Heat stress, a prevalent abiotic factor, significantly influences cotton plant biological activities resulting in reducing yield and production. We must deepen our understanding of how plants respond to heat stress across various dimensions, encompassing genes, RNAs, proteins, metabolites for effective cotton breeding. Multi-omics methods, primarily genomics, transcriptomics, proteomics, metabolomics, and phenomics, proves instrumental in studying cotton's responses to abiotic stresses. Integrating genomics, transcriptomics, proteomics, and metabolomic is imperative for our better understanding regarding genetics and molecular basis of heat tolerance in cotton. The current review explores fundamental omics techniques, covering genomics, transcriptomics, proteomics, and metabolomics, to highlight the progress made in cotton omics research.

Progress in Rice Breeding Based on Genomic Research

The role of rice genomics in breeding progress is becoming increasingly important. Deeper research into the rice genome will contribute to the identification and utilization of outstanding functional genes, enriching the diversity and genetic basis of breeding materials and meeting the diverse demands for various improvements. Here, we review the significant contributions of rice genomics research to breeding progress over the last 25 years, discussing the profound impact of genomics on rice genome sequencing, functional gene exploration, and novel breeding methods, and we provide valuable insights for future research and breeding practices.

Optimizing an efficient ensemble approach for high-quality de novo transcriptome assembly of Thymus daenensis

Non-erroneous and well-optimized transcriptome assembly is a crucial prerequisite for authentic downstream analyses. Each de novo assembler has its own algorithm-dependent pros and cons to handle the assembly issues and should be specifically tested for each dataset. Here, we examined efficiency of seven state-of-art assemblers on ~ 30 Gb data obtained from mRNA-sequencing of Thymus daenensis. In an ensemble workflow, combining the outputs of different assemblers associated with an additional redundancy-reducing step could generate an optimized outcome in terms of completeness, annotatability, and ORF richness. Based on the normalized scores of 16 benchmarking metrics, EvidentialGene, BinPacker, Trinity, rnaSPAdes, CAP3, IDBA-trans, and Velvet-Oases performed better, respectively. EvidentialGene, as the best assembler, totally produced 316,786 transcripts, of which 235,730 (74%) were predicted to have a unique protein hit (on uniref100), and also half of its transcripts contained an ORF. The total number of unique BLAST hits for EvidentialGene was approximately three times greater than that of the worst assembler (Velvet-Oases). EvidentialGene could even capture 17% and 7% more average BLAST hits than BinPacker and Trinity. Although BinPacker and CAP3 produced longer transcripts, the EvidentialGene showed a higher collinearity between transcript size and ORF length. Compared with the other programs, EvidentialGene yielded a higher number of optimal transcript sets, further full-length transcripts, and lower possible misassemblies. Our finding corroborates that in non-model species, relying on a single assembler may not give an entirely satisfactory result. Therefore, this study proposes an ensemble approach of accompanying EvidentialGene pipelines to acquire a superior assembly for T. daenensis.

Comparative transcriptomic analysis provides key genetic resources in clove basil (Ocimum gratissimum) under cadmium stress

Planting aromatic plant might be a promising strategy for safely utilizing heavy metal (HM)-contaminated soils, as HMs in essential oil could be completely excluded using some special technologies with ease. Clove basil (Ocimum gratissimum L.) is an important aromatic plant used in essential oil production. Improving cadmium (Cd) tolerance in clove basil can increase its production and improve the utilization efficiency of Cd-contaminated soils. However, the lack of genomic information on clove basil greatly restricts molecular studies and applications in phytoremediation. In this study, we demonstrated that high levels of Cd treatments (0.8, 1.6 and 6.5 mg/L) significantly impacted the growth and physiological attributes of clove basil. Cd contents in clove basil tissues increased with treatment concentrations. To identify Cd stress-responsive genes, we conducted a comparative transcriptomic analysis using seedlings cultured in the Hoagland's solution without Cd ion (control) or containing 1.6 mg/L CdCl2 (a moderate concentration of Cd stress for clove basil seedlings). A total of 104.38 Gb clean data with high-quality were generated in clove basil under Cd stress through Illumina sequencing. More than 1,800 differential expressed genes (DEGs) were identified after Cd treatment. The reliability and reproducibility of the transcriptomic data were validated through qRT-PCR analysis and Sanger sequencing. KEGG classification analysis identified the "MAPK signaling pathway," "plant hormone signal transduction" and "plant-pathogen interaction" as the top three pathways. DEGs were divided into five clusters based on their expression patterns during Cd stress. The functional annotation of DEGs indicated that downregulated DEGs were mainly involved in the "photosynthesis system," whereas upregulated DEGs were significantly assigned to the "MAPK signaling pathway" and "plant-pathogen interaction pathway." Furthermore, we identified a total of 78 transcription factors (TFs), including members of bHLH, WRKY, AP2/ERF, and MYB family. The expression of six bHLH genes, one WRKY and one ERF genes were significantly induced by Cd stress, suggesting that these TFs might play essential roles in regulating Cd stress responses. Overall, our study provides key genetic resources and new insights into Cd adaption mechanisms in clove basil.

Transcriptional reprogramming during recovery from drought stress in Eucalyptus grandis

The importance of drought as a constraint to agriculture and forestry is increasing with climate change. Genetic improvement of plants' resilience is one of the mitigation strategies to curb this threat. Although recovery from drought stress is important to long-term drought adaptation and has been considered as an indicator of dehydration tolerance in annual crops, this has not been well explored in forest trees. Thus, we aimed to investigate the physiological and transcriptional changes during drought stress and rewatering in Eucalyptus grandis W. Hill ex Maiden. We set up a greenhouse experiment where we imposed drought stress on 2-year-old seedlings and rewatered the recovery group after 17 days of drought. Our measurement of leaf stomatal conductance (gs) showed that, while gs was reduced by drought stress, it fully recovered after 5 days of rewatering. The RNA-seq analysis from stem samples revealed that genes related to known stress responses such as phytohormone and reactive oxygen species signaling were upregulated, while genes involved in metabolism and growth were downregulated due to drought stress. We observed reprogramming of signal transduction pathways and metabolic processes at 1 day of rewatering, indicating a quick response to rewatering. Our results suggest that recovery from drought stress may entail alterations in the jasmonic acid, salicylic acid, ethylene and brassinosteroid signaling pathways. Using co-expression network analysis, we identified hub genes, including the putative orthologs of ABI1, ABF2, ABF3, HAI2, BAM1, GolS2 and SIP1 during drought and CAT2, G6PD1, ADG1 and FD-1 during recovery. Taken together, by highlighting the molecular processes and identifying key genes, this study gives an overview of the mechanisms underlying the response of E. grandis to drought stress and recovery that trees may face repeatedly throughout their long life cycle. This provides a useful reference to the identification and further investigation of signaling pathways and target genes for future tree improvement.

Untargeted Metabolomic Approach to Determine the Regulatory Pathways on Salicylic Acid-Mediated Stress Response in Aphanamixis polystachya Seedlings

Plants thrive under abiotic and biotic stress conditions with the changes in phytohormones like salicylic acid (SA), resulting in the synthesis of secondary metabolites. The present study determines the response of plants in producing secondary metabolites towards different SA concentrations at varying time intervals. Liquid chromatography-mass spectrometry-based metabolomics studies in Aphanamixis polystachya (Wall.) Parker seedlings are grown at 10 mM, 50 mM, and 100 mM SA concentrations, showed the differential expression of metabolites towards the stress. Alkaloids like amaranthin showed a 15-fold increase on the second day, and analog of androvinblastin showed a 20-fold increase on the sixth day in 10 mM SA compared with other groups. Flavanoid cyanidin 3-3 glucosyl was found to be with a 22-fold increment along with terpenoids betavulgaroside (18-fold), asiaticoside (17-fold), mubenin B (20-fold), and deslanoside (22-fold) increment in 50 mM SA on the sixth day. The shock exerted by 100 mM was too harsh, and the lowered metabolite production level was insufficient for the seedlings to survive at this higher SA condition. Arrangement of stressed groups using Pearson correlation studies, principal component analysis, and partial least square analysis placed 10 mM SA and controlled group closer and 50 mM SA and 100 mM SA groups closer to each other. The study observed that SA regulates metabolites that mediate biotic stress responses at low concentrations, and higher concentrations regulate abiotic stress regulating metabolites.

Transcriptome profiling revealed candidate genes, pathways and transcription factors related to nitrogen utilization and excessive nitrogen stress in perennial ryegrass

Ryegrass (Lolium perenne L.), a high-quality forage grass, is a good nutrient source for herbivorous livestock. However, improving nitrogen use efficiency and avoiding nitrate toxicity caused by excessive nitrogen are continual challenges in ryegrass production. The molecular mechanism underlying the response of ryegrass to nitrogen, especially excessive nitrogen, remains unclear. In this study, the transcriptomic changes under different nitrogen levels were investigated in perennial ryegrass by high-throughput next-generation RNA sequencing. Phenotypic characterization showed that treatment with half of the standard N concentration (N0.5) led to a better growth state than the other three treatments. The treatments with the standard N concentration (N1) and treatments with ten times higher than the standard N concentration (N10) contained excessive nitrogen, which placed stress on plant growth. Analysis of differentially expressed genes indicated that 345 and 104 genes are involved in the regulation of nitrogen utilization and excessive nitrogen stress, respectively. KEGG enrichment analysis suggested that "photosynthesis-antenna proteins" may respond positively to appropriate nitrogen conditions, whereas "steroid biosynthesis", "carotenoid biosynthesis" and "C5-branched dibasic acid metabolism" were identified as the top significantly enriched pathways in response to excessive nitrogen. Additionally, 21 transcription factors (TFs) related to nitrogen utilization were classified into 10 families, especially the AP2-EREBP and MYB TF families. Four TFs related to excessive nitrogen stress were identified, including LOBs, NACs, AP2-EREBPs and HBs. The expression patterns of these selected genes were also analyzed. These results provide new insight into the regulatory mechanism of ryegrass in response to nitrogen utilization and excessive nitrogen stress.

How to Choose a Good Marker to Analyze the Olive Germplasm (Olea europaea L.) and Derived Products

The olive tree (Olea europaea L.) is one of the most cultivated crops in the Mediterranean basin. Its economic importance is mainly due to the intense production of table olives and oil. Cultivated varieties are characterized by high morphological and genetic variability and present a large number of synonyms and homonyms. This necessitates the introduction of a rapid and accurate system for varietal identification. In the past, the recognition of olive cultivars was based solely on analysis of the morphological traits, however, these are highly influenced by environmental conditions. Therefore, over the years, several methods based on DNA analysis were developed, allowing a more accurate and reliable varietal identification. This review aims to investigate the evolving history of olive tree characterization approaches, starting from the earlier morphological methods to the latest technologies based on molecular markers, focusing on the main applications of each approach. Furthermore, we discuss the impact of the advent of next generation sequencing and the recent sequencing of the olive genome on the strategies used for the development of new molecular markers.

Selection and Validation of Reference Genes for qRT-PCR Analysis in the Oil-Rich Tuber Crop Tiger Nut (Cyperus esculentus) Based on Transcriptome Data

Tiger nut (Cyperus esculentus), a perennial C₄ plant of the Cyperaceae family, is an unconventional crop that is distinguished by its oil-rich tubers, which also possesses the advantages of strong resistance, wide adaptability, short life periods, and large biomass. To facilitate studies on gene expression in this species, we identified and validated a series of reference genes (RGs) based on transcriptome data, which can be employed as internal controls for qRT-PCR analysis in tiger nut. Fourteen putative candidate RGs were identified and evaluated across nine different tissues of two cultivars, and the RGs were analyzed using three different algorithms (geNorm, NormFinder, and BestKeeper). The stability rankings of the candidate RGs were merged into consensus lists with RankAggreg. For the below-ground storage organ of tiger nut, the optimal RGs were TUB4 and UCE2 in different developmental stages of tubers. UCE2 and UBL5 were the most stably expressed RGs among all tissues, while Rubisco and PGK exhibited the lowest expression stability. UCE2, UBL5 and Rubisco were compared to normalize the expression levels of the caleosin (CLO) and diacylglycerol acyltransferase 2-2 (DGAT2-2) genes across the same tissues. Our results showed that the RGs identified in this study, which exhibit more uniform expression patterns, may be utilized for the normalization of qRT-PCR results, promoting further research on gene expression in various tissues of tiger nut.

DiCoExpress: a tool to process multifactorial RNAseq experiments from quality controls to co-expression analysis through differential analysis based on contrasts inside GLM models

BACKGROUND

RNAseq is nowadays the method of choice for transcriptome analysis. In the last decades, a high number of statistical methods, and associated bioinformatics tools, for RNAseq analysis were developed. More recently, statistical studies realised neutral comparison studies using benchmark datasets, shedding light on the most appropriate approaches for RNAseq data analysis.

RESULTS

DiCoExpress is a script-based tool implemented in R that includes methods chosen based on their performance in neutral comparisons studies. DiCoExpress uses pre-existing R packages including FactoMineR, edgeR and coseq, to perform quality control, differential, and co-expression analysis of RNAseq data. Users can perform the full analysis, providing a mapped read expression data file and a file containing the information on the experimental design. Following the quality control step, the user can move on to the differential expression analysis performed using generalized linear models thanks to the automated contrast writing function. A co-expression analysis is implemented using the coseq package. Lists of differentially expressed genes and identified co-expression clusters are automatically analyzed for enrichment of annotations provided by the user. We used DiCoExpress to analyze a publicly available RNAseq dataset on the transcriptional response of Brassica napus L. to silicon treatment in plant roots and mature leaves. This dataset, including two biological factors and three replicates for each condition, allowed us to demonstrate in a tutorial all the features of DiCoExpress.

CONCLUSIONS

DiCoExpress is an R script-based tool allowing users to perform a full RNAseq analysis from quality controls to co-expression analysis through differential analysis based on contrasts inside generalized linear models. DiCoExpress focuses on the statistical modelling of gene expression according to the experimental design and facilitates the data analysis leading the biological interpretation of the results.

Analysis of Rice Proteins with DLN Repressor Motif/S

Transcriptional regulation includes both activation and repression of downstream genes. In plants, a well-established class of repressors are proteins with an ERF-associated amphiphilic repression/EAR domain. They contain either DLNxxP or LxLxL as the identifying hexapeptide motif. In rice (Oryza sativa), we have identified a total of 266 DLN repressor proteins, with the former motif and its modifications thereof comprising 227 transcription factors and 39 transcriptional regulators. Apart from DLNxxP motif conservation, DLNxP and DLNxxxP motifs with variable numbers/positions of proline and those without any proline conservation have been identified. Most of the DLN repressome proteins have a single DLN motif, with higher relative percentage in the C-terminal region. We have designed a simple yeast-based experiment wherein a DLN motif can successfully cause strong repression of downstream reporter genes, when fused to a transcriptional activator of rice or yeast. The DLN hexapeptide motif is essential for repression, and at least two “DLN” residues cause maximal repression. Comparatively, rice has more DLN repressor encoding genes than Arabidopsis, and DLNSPP motif from rice is 40% stronger than the known Arabidopsis SRDX motif. The study reports a straightforward assay to analyze repressor activity, along with the identification of a strong DLN repressor from rice.

May the Fittest Protein Evolve: Favoring the Plant-Specific Origin and Expansion of NAC Transcription Factors

Plant-specific NAC transcription factors (TFs) evolve during the transition from aquatic to terrestrial plant life and are amplified to become one of the biggest TF families. This is because they regulate genes involved in water conductance and cell support. They also control flower and fruit formation. The review presented here focuses on various properties, regulatory intricacies, and developmental roles of NAC family members. Processes controlled by NACs depend majorly on their transcriptional properties. NACs can function as both activators and/or repressors. Additionally, their homo/hetero dimerization abilities can also affect DNA binding and activation properties. The active protein levels are dependent on the regulatory cascades. Because NACs regulate both development and stress responses in plants, in-depth knowledge about them has the potential to help guide future crop improvement studies.

Gene Discovery of Characteristic Metabolic Pathways in the Tea Plant (Camellia sinensis) Using 'Omics'-Based Network Approaches: A Future Perspective

Characteristic secondary metabolites, including flavonoids, theanine and caffeine, in the tea plant (Camellia sinensis) are the primary sources of the rich flavors, fresh taste, and health benefits of tea. The decoding of genes involved in these characteristic components is still significantly lagging, which lays an obstacle for applied genetic improvement and metabolic engineering. With the popularity of high-throughout transcriptomics and metabolomics, 'omics'-based network approaches, such as gene co-expression network and gene-to-metabolite network, have emerged as powerful tools for gene discovery of plant-specialized (secondary) metabolism. Thus, it is pivotal to summarize and introduce such system-based strategies in facilitating gene identification of characteristic metabolic pathways in the tea plant (or other plants). In this review, we describe recent advances in transcriptomics and metabolomics for transcript and metabolite profiling, and highlight 'omics'-based network strategies using successful examples in model and non-model plants. Further, we summarize recent progress in 'omics' analysis for gene identification of characteristic metabolites in the tea plant. Limitations of the current strategies are discussed by comparison with 'omics'-based network approaches. Finally, we demonstrate the potential of introducing such network strategies in the tea plant, with a prospects ending for a promising network discovery of characteristic metabolite genes in the tea plant.

Comparative Proteomic Analysis of Plant Acclimation to Six Different Long-Term Environmental Changes

Plants are constantly challenged in their natural environment by a range of changing conditions. We investigated the acclimation processes and adaptive plant responses to various long-term mild changes and compared them directly within one experimental set-up. Arabidopsis thaliana plants were grown in hydroponic culture for 10 d under controlled abiotic stress (15°C, 25°C, salt and osmotic) and in nutrient deficiency (nitrate and phosphate). Plant growth was monitored and proteomic experiments were performed. Resource allocation between tissues altered during the plants' response. The growth patterns and induced changes of the proteomes indicated that the underlying mechanisms of the adaptation processes are highly specific to the respective environmental condition. Our results indicated differential regulation of response to salt and osmotic treatment, while the proteins in the changed temperature regime showed an inverse, temperature-sensitive control. There was a high correlation of protein level between the nutrient-deficient treatments, but the enriched pathways varied greatly. The proteomic analysis also revealed new insights into the regulation of proteins specific to the shoot and the root. Our investigation revealed unique strategies of plant acclimation to the different applied treatments on a physiological and proteome level, and these strategies are quite distinct in tissues below and above ground.

rosettR: protocol and software for seedling area and growth analysis

BACKGROUND

Growth is an important parameter to consider when studying the impact of treatments or mutations on plant physiology. Leaf area and growth rates can be estimated efficiently from images of plants, but the experiment setup, image analysis, and statistical evaluation can be laborious, often requiring substantial manual effort and programming skills.

RESULTS

Here we present rosettR, a non-destructive and high-throughput phenotyping protocol for the measurement of total rosette area of seedlings grown in plates in sterile conditions. We demonstrate that our protocol can be used to accurately detect growth differences among different genotypes and in response to light regimes and osmotic stress. rosettR is implemented as a package for the statistical computing software R and provides easy to use functions to design an experiment, analyze the images, and generate reports on quality control as well as a final comparison across genotypes and applied treatments. Experiment procedures are included as part of the package documentation.

CONCLUSIONS

Using rosettR it is straight-forward to perform accurate, reproducible measurements of rosette area and relative growth rate with high-throughput using inexpensive equipment. Suitable applications include screening mutant populations for growth phenotypes visible at early growth stages and profiling different genotypes in a wide variety of treatments.

Genome-Wide Identification, Localization, and Expression Analysis of Proanthocyanidin-Associated Genes in Brassica

Proanthocyanidins (PA) is a type of prominent flavonoid compound deposited in seed coats which controls the pigmentation in all Brassica species. Annotation of Brassica juncea genome survey sequences showed 72 PA genes; however, a functional description of these genes, especially how their interactions regulate seed pigmentation, remains elusive. In the present study, we designed 19 primer pairs to screen a bacterial artificial chromosome (BAC) library of B. juncea. A total of 284 BAC clones were identified and sequenced. Alignment of the sequences confirmed that 55 genes were cloned, with every Arabidopsis PA gene having 2-7 homologs in B. juncea. BLAST analysis using the recently released B. rapa or B. napus genome database identified 31 and 58 homologous genes, respectively. Mapping and phylogenetic analysis indicated that 30 B. juncea PA genes are located in the A-genome chromosomes except A04, whereas the remaining 25 genes are mapped to the B-genome chromosomes except B05 and B07. RNA-seq data and Fragments Per Kilobase of a transcript per Million mapped reads (FPKM) analysis showed that most of the PA genes were expressed in the seed coat of B. juncea and B. napus, and that BjuTT3, BjuTT18, BjuANR, BjuTT4-2, BjuTT4-3, BjuTT19-1, and BjuTT19-3 are transcriptionally regulated, and not expressed or downregulated in yellow-seeded testa. Importantly, our study facilitates in better understanding of the molecular mechanism underlying Brassica PA profiles and accumulation, as well as in further characterization of PA genes.

Rice Improvement Through Genome-Based Functional Analysis and Molecular Breeding in India

Rice is one of the main pillars of food security in India. Its improvement for higher yield in sustainable agriculture system is also vital to provide energy and nutritional needs of growing world population, expected to reach more than 9 billion by 2050. The high quality genome sequence of rice has provided a rich resource to mine information about diversity of genes and alleles which can contribute to improvement of useful agronomic traits. Defining the function of each gene and regulatory element of rice remains a challenge for the rice community in the coming years. Subsequent to participation in IRGSP, India has continued to contribute in the areas of diversity analysis, transcriptomics, functional genomics, marker development, QTL mapping and molecular breeding, through national and multi-national research programs. These efforts have helped generate resources for rice improvement, some of which have already been deployed to mitigate loss due to environmental stress and pathogens. With renewed efforts, Indian researchers are making new strides, along with the international scientific community, in both basic research and realization of its translational impact.

Stress-induced and epigenetic-mediated maize transcriptome regulation study by means of transcriptome reannotation and differential expression analysis

Plant's response and adaptation to abiotic stresses involve sophisticated genetic and epigenetic regulatory systems. To obtain a global view of molecular response to osmotic stresses, including the non-coding portion of genome, we conducted a total leaf transcriptome analysis on maize plants subjected to prolonged drought and salt stresses. Stress application to both B73 wild type and the epiregulator mutant rpd1-1/rmr6 allowed dissection of the epigenetic component of stress response. Coupling total RNA-Seq and transcriptome re-assembly we annotated thousands of new maize transcripts, together with 13,387 lncRNAs that may play critical roles in regulating gene expression. Differential expression analysis revealed hundreds of genes modulated by long-term stress application, including also many lncRNAs and transposons specifically induced by stresses. The amplitude and dynamic of the stress-modulated gene sets are very different between B73 and rpd1-1/rmr6 mutant plants, as result of stress-like effect on genome regulation caused by the mutation itself, which activates many stress-related genes even in control condition. The analyzed extensive set of total RNA-Seq data, together with the improvement of the transcriptome and the identification of the non-coding portion of the transcriptome give a revealing insight into the genetic and epigenetic mechanism responsible for maize molecular response to abiotic stresses.

Selection of reliable reference genes for RT-qPCR analysis during developmental stages and abiotic stress in Setaria viridis

Real-time PCR (RT-qPCR) expression analysis is a powerful analytical technique, but reliable results depend on the use of stable reference genes for proper normalization. This study proposed to test the expression stability of 13 candidate reference genes in Setaria viridis, a monocot species recently proposed as a new C4 model plant. Gene expression stability of these genes was assayed across different tissues and developmental stages of Setaria and under drought or aluminum stress. In general, our results showed Protein Kinase, RNA Binding Protein and SDH as the most stable genes. Moreover, pairwise analysis showed that two reference genes were sufficient to normalize the gene expression data under each condition. By contrast, GAPDH and ACT were the least stably expressed genes tested. Validation of suitable reference genes was carried out to profile the expression of P5CS and GolS during abiotic stress. In addition, normalization of gene expression of SuSy, involved in sugar metabolism, was assayed in the developmental dataset. This study provides a list of reliable reference genes for transcript normalization in S. viridis in different tissues and stages of development and under abiotic stresses, which will facilitate genetic studies in this monocot model plant.

Ascribing Functions to Genes: Journey Towards Genetic Improvement of Rice Via Functional Genomics

Rice, one of the most important cereal crops for mankind, feeds more than half the world population. Rice has been heralded as a model cereal owing to its small genome size, amenability to easy transformation, high synteny to other cereal crops and availability of complete genome sequence. Moreover, sequence wealth in rice is getting more refined and precise due to resequencing efforts. This humungous resource of sequence data has confronted research fraternity with a herculean challenge as well as an excellent opportunity to functionally validate expressed as well as regulatory portions of the genome. This will not only help us in understanding the genetic basis of plant architecture and physiology but would also steer us towards developing improved cultivars. No single technique can achieve such a mammoth task. Functional genomics through its diverse tools viz. loss and gain of function mutants, multifarious omics strategies like transcriptomics, proteomics, metabolomics and phenomics provide us with the necessary handle. A paradigm shift in technological advances in functional genomics strategies has been instrumental in generating considerable amount of information w.r.t functionality of rice genome. We now have several databases and online resources for functionally validated genes but despite that we are far from reaching the desired milestone of functionally characterizing each and every rice gene. There is an urgent need for a common platform, for information already available in rice, and collaborative efforts between researchers in a concerted manner as well as healthy public-private partnership, for genetic improvement of rice crop better able to handle the pressures of climate change and exponentially increasing population.

FRAMA: from RNA-seq data to annotated mRNA assemblies

Background

Advances in second-generation sequencing of RNA made a near-complete characterization of transcriptomes affordable. However, the reconstruction of full-length mRNAs via de novo RNA-seq assembly is still difficult due to the complexity of eukaryote transcriptomes with highly similar paralogs and multiple alternative splice variants. Here, we present FRAMA, a genome-independent annotation tool for de novo mRNA assemblies that addresses several post-assembly tasks, such as reduction of contig redundancy, ortholog assignment, correction of misassembled transcripts, scaffolding of fragmented transcripts and coding sequence identification.

Results

We applied FRAMA to assemble and annotate the transcriptome of the naked mole-rat and assess the quality of the obtained compilation of transcripts with the aid of publicy available naked mole-rat gene annotations.

Based on a de novo transcriptome assembly (Trinity), FRAMA annotated 21,984 naked mole-rat mRNAs (12,100 full-length CDSs), corresponding to 16,887 genes. The scaffolding of 3488 genes increased the median sequence information 1.27-fold. In total, FRAMA detected and corrected 4774 misassembled genes, which were predominantly caused by fusion of genes. A comparison with three different sources of naked mole-rat transcripts reveals that FRAMA’s gene models are better supported by RNA-seq data than any other transcript set. Further, our results demonstrate the competitiveness of FRAMA to state of the art genome-based transcript reconstruction approaches.

Conclusion

FRAMA realizes the de novo construction of a low-redundant transcript catalog for eukaryotes, including the extension and refinement of transcripts. Thereby, results delivered by FRAMA provide the basis for comprehensive downstream analyses like gene expression studies or comparative transcriptomics. FRAMA is available at https://github.com/gengit/FRAMA.

Electronic supplementary material

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

Three Rice NAC Transcription Factors Heteromerize and Are Associated with Seed Size

NACs are plant-specific transcription factors (TFs) involved in multiple aspects of development and stress. In rice, three NAC TF encoding genes, namely ONAC020, ONAC026, and ONAC023 express specifically during seed development, at extremely high levels. They exhibit significantly strong association with seed size/weight with the sequence variations located in the upstream regulatory region. Concomitantly, their expression pattern/levels during seed development vary amongst different accessions with variation in seed size. The alterations in the promoter sequences of the three genes, amongst the five rice accessions, correlate with the expression levels to a certain extent only. In terms of transcriptional properties, the three NAC TFs can activate and/or suppress downstream genes, though to different extents. Only ONAC026 is localized to the nucleus while ONAC020 and ONAC023 are targeted to the ER and cytoplasm, respectively. Interestingly, these two proteins interact with ONAC026 and the dimers localize in the nucleus. Trans-splicing between ONAC020 and ONAC026 results in three additional forms of ONAC020. The transcriptional properties including activation, repression, subcellular localization and heterodimerization of trans-spliced forms of ONAC020 and ONAC026 are different, indicating toward their role as competitors. The analysis presented in this paper helps to conclude that the three NAC genes, which are associated with seed size, have independent as well as overlapping roles during the process and can be exploited as potential targets for crop improvement.

Role of Protein Tyrosine Phosphatases in Plants

Reversible protein phosphorylation is a crucial regulatory mechanism that controls many biological processes in eukaryotes. In plants, phosphorylation events primarily occur on serine (Ser) and threonine (Thr) residues, while in certain cases, it was also discovered on tyrosine (Tyr) residues. In contrary to plants, extensive reports on Tyr phosphorylation regulating a large numbers of biological processes exist in animals. Despite of such prodigious function in animals, Tyr phosphorylation is a least studied mechanism of protein regulation in plants. Recently, various chemical analytical procedures have strengthened the view that Tyr phosphorylation is equally prevalent in plants as in animals. However, regardless of Tyr phosphorylation events occuring in plants, no evidence could be found for the existence of gene encoding for Tyr phosphorylation i.e. the typical Tyr kinases. Various methodologies have suggested that plant responses to stress signals and developmental processes involved modifications in protein Tyr phosphorylation. Correspondingly, various reports have established the role of PTPs (Protein Tyrosine Phosphatases) in the dephosphorylation and inactivation of mitogen activated protein kinases (MAPKs) hence, in the regulation of MAPK signaling cascade. Besides this, many dual specificity protein phosphatases (DSPs) are also known to bind starch and regulate starch metabolism through reversible phosphorylation. Here, we are emphasizing the significant progress on protein Tyr phosphatases to understand the role of these enzymes in the regulation of post-translational modification in plant physiology and development.