An RNA isolation system for plant tissues rich in secondary metabolites

Gene expression analysis for selection and validation of suitable housekeeping gene(s) in cadmium exposed pigeonpea plants inoculated with arbuscular mycorrhizae

The expression stability of six commonly used housekeeping genes (18S rRNA-18S ribosomal RNA, EF1α-elongation factor 1α, ACT1-Actin 1, GAPDH-Glyceraldehyde-3-phosphate dehydrogenase, TUB6-Tubulin/FtsZ family and UBC-Ubiquitin-conjugating enzyme) were scrutinized in leaves and roots of Cd stressed pigeonpea plants inoculated with arbuscular mycorrhizal (AM) species- Rhizoglomus intraradices (Ri), Funneliformis mosseae (Fm), Claroideoglomus etunicatum (Ce), C. claroideum (Cc). The stability profile of each gene was assessed using ΔC_t, BestKeeper, NormFinder, RefFinder and geNorm algorithmic programs, which ranked different genes as most and least stable according to the tissues analysed. All the statistical algorithms ranked TUB6 as most stable and EF1α least stable housekeeping (HK) genes in both the plant tissues. The selected HK genes were verified using metallothionein (CcMT1) i.e. a stress responsive gene, whose expression altered under conditions of metal stress and AM inoculation. The expression pattern of CcMT1 varied highly when least stable reference gene was used for normalization as compared to most stable gene, under different treatments. Thus, there is a need of selecting suitable reference gene to achieve reliable results in gene expression studies using quantitative real time PCR (qRT-PCR). The study conducted will help future gene expression analysis in pigeonpea under specific stress.

Underpinning the molecular programming attributing heat stress associated thermotolerance in tea (Camellia sinensis (L.) O. Kuntze)

The most daunting issue of global climate change is the deleterious impact of extreme temperatures on tea productivity and quality, which has resulted in a quest among researchers and growers. The current study aims to unravel molecular programming underpinning thermotolerance by characterizing heat tolerance and sensitivity response in 20 tea cultivars. The significantly higher negative influence of heat stress was recorded in a sensitive cultivar with reduced water retention (47%), chlorophyll content (33.79%), oxidation potential (32.48%), and increase in membrane damage (76.4%). Transcriptional profiling of most tolerant and sensitive cultivars identified 78 differentially expressed unigenes with chaperon domains, including low and high molecular weight heat shock protein (HSP) and heat shock transcription factors (HSFs) involved in heat shock response (HSR). Further, predicted transcriptional interactome network revealed their key role in thermotolerance via well-co-ordinated transcriptional regulation of aquaporins, starch metabolism, chlorophyll biosynthesis, calcium, and ethylene mediated plant signaling system. The study identified the key role of HSPs (CsHSP90) in regulating HSR in tea, wherein, structure-based molecular docking revealed the inhibitory role of geldanamycin (GDA) on CsHSP90 by blocking ATP binding site at N-terminal domain of predicted structure. Subsequently, GDA mediated leaf disc inhibitor assay further affirmed enhanced HSR with higher expression of CsHSP17.6, CsHSP70, HSP101, and CsHSFA2 genes in tea. Through the current study, efforts were made to extrapolate a deeper understanding of chaperons mediated regulation of HSR attributing thermotolerance in tea.

Integrated Metabolomics and Transcriptomics Using an Optimised Dual Extraction Process to Study Human Brain Cancer Cells and Tissues

The integration of untargeted metabolomics and transcriptomics from the same population of cells or tissue enhances the confidence in the identified metabolic pathways and understanding of the enzyme–metabolite relationship. Here, we optimised a simultaneous extraction method of metabolites/lipids and RNA from ependymoma cells (BXD-1425). Relative to established RNA (mirVana kit) or metabolite (sequential solvent addition and shaking) single extraction methods, four dual-extraction techniques were evaluated and compared (methanol:water:chloroform ratios): cryomill/mirVana (1:1:2); cryomill-wash/Econospin (5:1:2); rotation/phenol-chloroform (9:10:1); Sequential/mirVana (1:1:3). All methods extracted the same metabolites, yet rotation/phenol-chloroform did not extract lipids. Cryomill/mirVana and sequential/mirVana recovered the highest amounts of RNA, at 70 and 68% of that recovered with mirVana kit alone. sequential/mirVana, involving RNA extraction from the interphase of our established sequential solvent addition and shaking metabolomics-lipidomics extraction method, was the most efficient approach overall. Sequential/mirVana was applied to study a) the biological effect caused by acute serum starvation in BXD-1425 cells and b) primary ependymoma tumour tissue. We found (a) 64 differentially abundant metabolites and 28 differentially expressed metabolic genes, discovering four gene-metabolite interactions, and (b) all metabolites and 62% lipids were above the limit of detection, and RNA yield was sufficient for transcriptomics, in just 10 mg of tissue.

Comparative transcriptome analysis of Rheum australe, an endangered medicinal herb, growing in its natural habitat and those grown in controlled growth chambers

Rheum australe is an endangered medicinal herb of high altitude alpine region of Himalayas and is known to possess anti-cancerous properties. Unlike many herbs of the region, R. australe has broad leaves. The species thrives well under the environmental extremes in its niche habitat, therefore an understanding of transcriptome of R. australe to environmental cues was of significance. Since, temperature is one of the major environmental variables in the niche of R. australe, transcriptome was studied in the species growing in natural habitat and those grown in growth chambers maintained at 4 °C and 25 °C to understand genes associated with different temperatures. A total of 39,136 primarily assembled transcripts were obtained from 10,17,74,336 clean read, and 21,303 unigenes could match to public databases. An analysis of transcriptome by fragments per kilobase of transcript per million, followed by validation through qRT-PCR showed 22.4% up- and 22.5% down-regulated common differentially expressed genes in the species growing under natural habitat and at 4 °C as compared to those at 25 °C. These genes largely belonged to signaling pathway, transporters, secondary metabolites, phytohormones, and those associated with cellular protection, suggesting their importance in imparting adaptive advantage to R. australe in its niche.

Identification and characterization of SlbHLH, SlDof and SlWRKY transcription factors interacting with SlDPD gene involved in costunolide biosynthesis in Saussurea lappa

The genes involved in costunolide biosynthesis in Saussurea lappa have been identified recently by our lab. However, the study of transcriptional regulators of these genes was lacking for better opportunities for engineering the pharmacologically important biosynthetic pathway. Therefore, we cloned the promoter region of diphosphomevalonate decarboxylase gene (DPD) and analyzed its cis-acting regulatory elements to reveal the potential transcription factor (TF) binding sites for Dof, bHLH and WRKY family proteins in the gene promoter. The transcriptome study approach followed by the hidden Markov model based search, digital gene expression, co-expression network analysis, conserved domain properties and evolutionary analyses were carried out to screen out seven putative TFs for the DPD-TF interaction studies. Yeast one-hybrid assays were performed and three TFs were reported, namely, SlDOF2, SlbHLH3 and SlWRKY2 from Dof, bHLH and WRKY families, respectively that interacted positively with the DPD gene of the costunolide biosynthetic pathway. The tissue specific relative gene expression studies also supported the linked co-expression of the gene and its interacting TFs The present report will improve the understanding of transcriptional regulation pattern of costunolide biosynthetic pathway.

Transcriptional analysis reveals key insights into seasonal induced anthocyanin degradation and leaf color transition in purple tea (Camellia sinensis (L.) O. Kuntze)

Purple-tea, an anthocyanin rich cultivar has recently gained popularity due to its health benefits and captivating leaf appearance. However, the sustainability of purple pigmentation and anthocyanin content during production period is hampered by seasonal variation. To understand seasonal dependent anthocyanin pigmentation in purple tea, global transcriptional and anthocyanin profiling was carried out in tea shoots with two leaves and a bud harvested during in early (reddish purple: S1_RP), main (dark gray purple: S2_GP) and backend flush (moderately olive green: S3_G) seasons. Of the three seasons, maximum accumulation of total anthocyanin content was recorded in S2_GP, while least amount was recorded during S3_G. Reference based transcriptome assembly of 412 million quality reads resulted into 71,349 non-redundant transcripts with 6081 significant differentially expressed genes. Interestingly, key DEGs involved in anthocyanin biosynthesis [PAL, 4CL, F3H, DFR and UGT/UFGT], vacuolar trafficking [ABC, MATE and GST] transcriptional regulation [MYB, NAC, bHLH, WRKY and HMG] and Abscisic acid signaling pathway [PYL and PP2C] were significantly upregulated in S2_GP. Conversely, DEGs associated with anthocyanin degradation [Prx and lac], repressor TFs and key components of auxin and ethylene signaling pathways [ARF, AUX/IAA/SAUR, ETR, ERF, EBF1/2] exhibited significant upregulation in S3_G, correlating positively with reduced anthocyanin content and purple coloration. The present study for the first-time elucidated genome-wide transcriptional insights and hypothesized the involvement of anthocyanin biosynthesis activators/repressor and anthocyanin degrading genes via peroxidases and laccases during seasonal induced leaf color transition in purple tea. Futuristically, key candidate gene(s) identified here can be used for genetic engineering and molecular breeding of seasonal independent anthocyanin-rich tea cultivars.

RNA Isolation from Articular Cartilage Tissue

Isolation of high-quality RNA directly from tissues is desirable to obtain precise information of in vivo gene expression profiles in cells embedded within their extracellular matrix (ECM). It is well known that purification of RNA from cartilage tissues is particularly challenging due to low cell (chondrocyte) content and its dense ECM rich in large negatively charged proteoglycans that can copurify with RNA. Older methodologies to purify RNA from cartilage involved the use of concentrated denaturing solutions containing guanidinium isothiocyanate followed by ultracentrifugation in cesium trifluoroacetate. Such ultracentrifugation approaches are rarely used now since the emergence of more user-friendly mini spin column chromatography kits. For this chapter, we tested and compared three methods to isolate RNA from immature murine articular (femoral head) cartilage and found that the combination of TRIzol® reagent and spin column chromatography (Norgen Total RNA Purification Kit) was the best approach to generate higher quality RNA. Here, the average RNA Integrity Number (RIN), as determined by Bioanalyzer technology, was 7.1. We then applied this method to attempt to isolate RNA directly from human articular cartilage harvested from three osteoarthritic (OA) knee joint specimens. As expected, the concentration and quality of RNA obtained differed between samples. However, from one specimen, we were able to isolate approximately 3 μg of total RNA (including small noncoding RNAs) from 100 mg of human OA cartilage with a RIN = 7.9. Despite the patient-to-patient variabilities that are known to exist between cartilage specimens from OA joints, we have demonstrated that it is possible to obtain reasonably high levels of RNA from human OA articular cartilage at a quality suitable for downstream analyses including microarray and RNA-Seq. A detailed description of our preferred RNA purification methodology, which can be used to isolate RNA from human, bovine, or rodent cartilage tissue, is provided in this chapter.

Multiple-genotypes transcriptional analysis revealed candidates genes and nucleotide variants for improvement of quality characteristics in tea (Camellia sinensis (L.) O. Kuntze)

Tea quality is a polygenic trait that exhibits tremendous genetic variability due to accumulation of array of secondary metabolites. To elucidate global molecular insights controlling quality attributes, metabolite profiling and transcriptome sequencing of twelve diverse tea cultivars was performed in tea shoots harvested during quality season. RP-HPLC-DAD analysis of quality parameters revealed significant difference in catechins, theanine and caffeine contents. Transcriptome sequencing resulted into 50,107 non-redundant transcripts with functional annotations of 81.6% (40,847) of the transcripts. Interestingly, 2872 differentially expressed transcripts exhibited significant enrichment in 38 pathways (FDR ≤ 0.05) including secondary metabolism, amino acid and carbon metabolism. Thirty-eight key candidates reportedly involved in biosynthesis of fatty acid derived volatiles, volatile terpenes, glycoside hydrolysis and key quality related pathways (flavonoid, caffeine and theanine-biosynthesis) were highly expressed in catechins-rich tea cultivars. Furthermore, enrichment of candidates involved in flavonoid biosynthesis, transcriptional regulation, volatile terpene and biosynthesis of fatty acid derived volatile in Protein-Protein Interactome network revealed well-coordinated regulation of quality characteristics in tea. Additionally, ascertainment of 23,649 non-synonymous SNPs and validation of candidate SNPs present in quality related genes suggests their potential utility in genome-wide mapping and marker development for expediting breeding of elite compound-rich tea cultivars.

De novo transcriptome analysis of the critically endangered alpine Himalayan herb Nardostachys jatamansi reveals the biosynthesis pathway genes of tissue-specific secondary metabolites

The study is the first report on de novo transcriptome analysis of Nardostachys jatamansi, a critically endangered medicinal plant of alpine Himalayas. Illumina GAIIx sequencing of plants collected during end of vegetative growth (August) yielded 48,411 unigenes. 74.45% of these were annotated using UNIPROT. GO enrichment analysis, KEGG pathways and PPI network indicated simultaneous utilization of leaf photosynthates for flowering, rhizome fortification, stress response and tissue-specific secondary metabolites biosynthesis. Among the secondary metabolite biosynthesis genes, terpenoids were predominant. UPLC-PDA analysis of in vitro plants revealed temperature-dependent, tissue-specific differential distribution of various phenolics. Thus, as compared to 25 °C, the phenolic contents of both leaves (gallic acid and rutin) and roots (p-coumaric acid and cinnamic acid) were higher at 15 °C. These phenolics accounted for the therapeutic properties reported in the plant. In qRT-PCR of in vitro plants, secondary metabolite biosynthesis pathway genes showed higher expression at 15 °C and 14 h/10 h photoperiod (conditions representing end of vegetative growth period). This provided cues for in vitro modulation of identified secondary metabolites. Such modulation of secondary metabolites in in vitro systems can eliminate the need for uprooting N. jatamansi from wild. Hence, the study is a step towards effective conservation of the plant.

Full-length transcriptome sequencing analysis and development of EST-SSR markers for the endangered species Populus wulianensis

Populus wulianensis is an endangered species endemic to Shandong Province, China. Despite the economic and ornamental value of this species, few genomics and genetic studies have been performed. In this study, we performed a relevant analysis of the full-length transcriptome sequencing data of P. wulianensis and obtained expressed sequence tag (EST)-simple sequence repeat (SSR) markers with polymorphisms that can be used for further genetic research. In total, 8.18 Gb (3,521,665) clean reads with an average GC content of 42.12% were obtained. From the corrected 64,737 high-quality isoforms, 42,323 transcript sequences were obtained after redundancy analysis with CD-HIT. Among these transcript sequences, 41,876 sequences were annotated successfully. A total of 23,539 potential EST-SSRs were identified from 16,057 sequences. Excluding mononucleotides, the most abundant motifs were trinucleotide SSRs (47.80%), followed by di- (46.80%), tetra- (2.98%), hexa- (1.58%) and pentanucleotide SSRs (0.84%). Among the 100 designed EST-SSRs, 18 were polymorphic with high PIC values (0.721 and 0.683) and could be used for analyses of the genetic diversity and population structure of P. wulianensis. These full-length transcriptome sequencing data will facilitate gene discovery and functional genomics research in P. wulianensis, and the novel EST-SSRs developed in our study will promote molecular-assisted breeding, genetic diversity and conservation biology research in this species.

PkGPPS.SSU interacts with two PkGGPPS to form heteromeric GPPS in Picrorhiza kurrooa: Molecular insights into the picroside biosynthetic pathway

Geranyl geranyl pyrophosphate synthase (GGPPS) is known to form an integral subunit of the heteromeric GPPS (geranyl pyrophosphate synthase) complex and catalyzes the biosynthesis of monoterpene in plants. Picrorhiza kurrooa Royle ex Benth., a medicinally important high altitude plant is known for picroside biomolecules, the monoterpenoids. However, the significance of heteromeric GPPS in P. kurrooa still remains obscure. Here, transient silencing of PkGGPPS was observed to reduce picroside-I (P-I) content by more than 60% as well as picroside-II (P-II) by more than 75%. Thus, PkGGPPS was found to be involved in the biosynthesis of P-I and P-II besides other terpenoids. To unravel the mechanism, small subunit of GPPS (PkGPPS.SSU) was identified from P. kurrooa. Protein-protein interaction studies in yeast as well as bimolecular fluorescence complementation (BiFC) in planta have indicated that large subunit of GPPS PkGPPS.LSUs (PkGGPPS1 and PkGGPPS2) and PkGPPS.SSU form a heteromeric GPPS. Presence of similar conserved domains such as light responsive motifs, low temperature responsive elements (LTRE), dehydration responsive elements (DREs), W Box and MeJA responsive elements in the promoters of PkGPPS.LSU and PkGPPS.SSU documented their involvement in picroside biosynthesis. Further, the tissue specific transcript expression analysis vis-à-vis epigenetic regulation (DNA methylation) of promoters as well as coding regions of PkGPPS.LSU and PkGPPS.SSU has strongly suggested their role in picroside biosynthesis. Taken together, the newly identified PkGPPS.SSU formed the heteromeric GPPS by interacting with PkGPPS.LSUs to synthesize P-I and P-II in P. kurrooa.

Double lysis: an integrative time-saving method yielding high-quality RNA from strawberry

Background

The isolation of high-quality RNA from strawberry leaves and fruits is notoriously cumbersome. Both tissues are extremely rich in active secondary metabolites, as phenolics and pigments that inevitably perturb the isolation of RNA. Many protocols have been developed to address this problem. However, they are either costly, or time-consuming, in particular for high number of many plant samples. We describe here a new method with an easy-to-handle approach to obtain high-quality RNA from strawberry leaves and fruits. The method, referred to as double lysis, uses a novel combination of CTAB and Trizol protocols.

Results

Compared to conventional Trizol-dependent protocols, either used individually, or in a commercial spin-column kit, the new method yields RNA at lower costs, but of significantly improved quality. The RNA obtained by double lysis was very pure as indicated by 260/280 ratios of 2.06 (leaves) and 2.07 (fruits), while 260/230 ratio was 2.07 and 1.75, respectively. Also visually, RNA sediments from double lysis showed a white color, indicative of efficient elimination of polyphenolics and pigments. In contrast, traditional Trizol method produced reddish precipitates. The purity of RNA isolated by double lysis enabled successful removal of genomic DNA and, thus, allowed for more efficient cDNA synthesis and RT-PCR. In addition, the suggested method is able to yield RNA with fully preserved integrity from strawberry leaves, fruits in addition to petals and roots.

Conclusion

Double lysis is a new RNA isolation protocol developed through the integrative coupling of CTAB and Trizol reagents. The method is cost-effective, robust, time-saving, and can cope even with recalcitrant plant tissues with high contents of phenolics and pigments, such as strawberry leaves and fruits.

Two light responsive WRKY genes exhibit positive and negative correlation with picroside content in Picrorhiza kurrooa Royle ex Benth, an endangered medicinal herb

Picrorhiza kurrooa is an endangered herb known to produce the medicinally important picrosides through isoprenoid pathway. The present work showed the functionality of WRKY motifs (TGAC cis-acting elements) present in the promoters of regulatory genes 3-hydroxy-3-methylglutaryl coenzyme A reductase (Pkhmgr) and 1-deoxy-d-xylulose-5-phosphate synthase (Pkdxs) of the picrosides biosynthetic pathway by electrophoretic mobility shift assay. Also, the two WRKY genes, PkdWRKY and PksWRKY, were characterized and found to contain double and single characteristic WRKY domains, respectively along with a zinc-finger motif in each domain. Expression analysis revealed that PkdWRKY and PksWRKY exhibited a positive and negative correlation, respectively, with picrosides content under the environment of light and in different tissues. Functional evaluation in yeast showed DNA binding ability of both PksWRKY and PkdWRKY; however, only PkdWRKY exhibited transcriptional activation ability. Transient overexpression of PkdWRKY and PksWRKY in tobacco modulated the expression of selected native genes of tobacco involved in MVA and MEP pathway suggesting functionality of PkdWRKY and PksWRKY in planta. Collectively, data suggested that PkdWRKY and PksWRKY might be positive and negative regulators, respectively in the picrosides biosynthetic pathway.

Hybrid de novo genome assembly of red gromwell (Lithospermum erythrorhizon) reveals evolutionary insight into shikonin biosynthesis

Lithospermum erythrorhizon (red gromwell; zicao) is a medicinal and economically valuable plant belonging to the Boraginaceae family. Roots from L. erythrorhizon have been used for centuries based on the antiviral and wound-healing properties produced from the bioactive compound shikonin and its derivatives. More recently, shikonin, its enantiomer alkannin, and several other shikonin/alkannin derivatives have collectively emerged as valuable natural colorants and as novel drug scaffolds. Despite several transcriptomes and proteomes having been generated from L. erythrorhizon, a reference genome is still unavailable. This has limited investigations into elucidating the shikonin/alkannin pathway and understanding its evolutionary and ecological significance. In this study, we obtained a de novo genome assembly for L. erythrorhizon using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies. The resulting genome is ∼367.41 Mb long, with a contig N50 size of 314.31 kb and 27,720 predicted protein-coding genes. Using the L. erythrorhizon genome, we identified several additional p-hydroxybenzoate:geranyltransferase (PGT) homologs and provide insight into their evolutionary history. Phylogenetic analysis of prenyltransferases suggests that PGTs originated in a common ancestor of modern shikonin/alkannin-producing Boraginaceous species, likely from a retrotransposition-derived duplication event of an ancestral prenyltransferase gene. Furthermore, knocking down expression of LePGT1 in L. erythrorhizon hairy root lines revealed that LePGT1 is predominantly responsible for shikonin production early in culture establishment. Taken together, the reference genome reported in this study and the provided analysis on the evolutionary origin of shikonin/alkannin biosynthesis will guide elucidation of the remainder of the pathway.

Comparative transcriptome analysis revealed gamma-irradiation mediated disruption of floral integrator gene(s) leading to prolonged vegetative phase in Stevia rebaudiana Bertoni

Stevia rebaudiana Bert. is getting global attention because of its ability to synthesize commercially important low/no calorie natural sweeteners (LNCSs) steviol glycosides (SGs). Considering, higher accumulation of SGs in vegetative phase followed by decrement during reproductive phase necessitate the understanding of different molecular components of floral transition to develop superior varieties/cultivars with prolonged vegetative phase in Stevia. Current comparative transcriptional analysis of low dose (5 kR) gamma-irradiated mutant genotype (SMG) with prolonged vegetative phase vis-à-vis background genotype (SBG) identified DGEs of major floral transition pathways, and expressed according to their physiological fate irrespective to SMG & SBG. Contrarily, reduced expression of floral integrator genes (FT and LEAFY) in mutant genotype suggests their involvement in prolonged vegetative phase phenotype. Likewise, GO and KEGG enrichment of photosynthesis and carbon assimilation efficiency might be associated with prolonged vegetative phase and higher accumulation of Stevioside content in mutant genotype. Furthermore, deviation of flowering related transcription factors (higher expressions except MIKS-type MADS-box SMG_PV compared to SBG_F) may possibly be correlated with low expression of floral integrator genes. Findings of current studies will facilitate the genetic manipulations and crop improvement efforts in Stevia through conventional breeding and genome editing approaches for increased SGs biosynthesis.

Identification of novel microRNAs and their targets in Chlorophytum borivilianum by small RNA and degradome sequencing

Plant specific miRNAs (Novel miRNAs) are well known to perform distinctive functions in biological processes. Identification of new miRNAs is necessary to understand their gene regulation. Degradome provides an opportunity to explore the miRNA functions by comparing the miRNA population and their degraded products. In the present study, Small RNA sequencing data was used to identify novel miRNAs. Further, degradome sequencing was carried out to identify miRNAs targets in the plant, Chlorophytum borivilianum. The present study supplemented 40 more novel miRNAs correlating degradome data with smallRNAome. Novel miRNAs, complementary to mRNA partial sequences obtained from degradome sequencing were actually targeting the later. A big pool of miRNA was established by using Oryza sativa, Arabidopsis thaliana, Populus trichocarpa, Ricinus communis, and Vitis vinifera genomic data. Targets were identified for novel miRNAs and total 109 targets were predicted. BLAST2GO analysis elaborate about localization of novel miRNAs’ targets and their corresponding KEGG (Kyoto Encyclopedia for Genes and Genomes) pathways. Identified targets were annotated and were found to be involved in significant biological processes like Nitrogen metabolism, Pyruvate metabolism, Citrate cycle (TCA cycle), photosynthesis, and Glycolysis/Gluconeogenesis. The present study provides an overall view of the miRNA regulation in multiple metabolic pathways that are involved in plant growth, pathogen resistance and secondary metabolism of C. borivilianum.

Recurrent gene duplication in the angiosperm tribe Delphinieae (Ranunculaceae) inferred from intracellular gene transfer events and heteroplasmic mutations in the plastid matK gene

The study of intracellular gene transfer may allow for the detection of interesting evolutionary processes such as ancient polyploidization. We compared 24 plastid genomes (plastomes) from tribe Delphinieae, one from tribe Nigelleae and one from tribe Ranunculeae, including five newly sequenced genomes. The functional transfers of the plastids rpl32 and rps16 to the nucleus in tribe Delphinieae were identified. Unexpectedly, we discovered multiple divergent copies of the nuclear-encoded plastid rpl32 in the genus Aconitum. Phylogenetic and synonymous substitution rate analyses revealed that the nuclear-encoded plastid rpl32 underwent two major duplication events. These ancient gene duplication events probably occurred via multiple polyploidization events in Aconitum between 11.9 and 24.7 Mya. Furthermore, our sequence rate analysis indicated that the eight plastid-encoded rpl subunits in Aconitum had a significantly accelerated evolutionary rate compared to those in other genera, suggesting that highly divergent paralogs targeted to the plastid may contribute to an elevated rate of evolution in plastid rpl genes. In addition, heteroplasmy of the plastid matK from two Aconitum species suggested the existence of potentially functional plastid maturases in its plastome. Our results provide insight into the evolutionary history of the tribe Delphinieae.

Regulation of color transition in purple tea (Camellia sinensis)

Comparative proteomics and metabolomics study of juvenile green, light purple and dark purple leaf to identify key proteins and metabolites that putatively govern color transition in Camellia sinensis. Color transition from juvenile green to dark purple leaf in Camellia sinensis is a complex process and thought to be regulated by an intricate balance of genes, proteins and metabolites expression. A molecular-level understanding of proteins and metabolites expression is needed to define metabolic process underpinning color transition in C. sinensis. Here, purple leaf growth of C. sinensis cultivar was divided into three developmental stages viz. juvenile green (JG), light purple (LP) and dark purple (DP) leaf. Scanning electron microscope (SEM) analysis revealed a clear morphological variation such as cell size, shape and texture as tea leaf undergoing color transition. Proteomic and metabolomic analyses displayed the temporal changes in proteins and metabolites that occur in color transition process. In total, 211 differentially expressed proteins (DEPs) were identified presumably involved in secondary metabolic processes particularly, flavonoids/anthocyanin biosynthesis, phytohormone regulation, carbon and nitrogen assimilation and photosynthesis, among others. Subcellular localization of three candidate proteins was further evaluated by their transient expression in planta. Interactome study revealed that proteins involved in primary metabolism, precursor metabolite, photosynthesis, phytohormones, transcription factor and anthocyanin biosynthesis were found to be interact directly or indirectly and thus, regulate color transition from JG to DP leaf. The present study not only corroborated earlier findings but also identified novel proteins and metabolites that putatively govern color transition in C. sinensis. These findings provide a platform for future studies that may be utilized for metabolic engineering/molecular breeding in an effort to develop more desirable traits.

Proteome Analysis of the Gametophytes of a Western Himalayan Fern Diplazium maximum Reveals Their Adaptive Responses to Changes in Their Micro-Environment

Ferns have survived changing habitats and environmental extremes of different eras, wherein, the exploratory haploid gametophytes are believed to have played a major role. Therefore, the proteome of in vitro grown gametophytes of a temperate Himalayan fern, Diplazium maximum in response to 0 (G0), 1 (G1), and 3% (G3) sucrose was studied. A total of 110 differentially abundant protein spots (DAPs) were obtained. Among these, only 67 could be functionally categorized as unique proteins involved in various metabolic processes. Calcium dependent proteins, receptor like kinases, G proteins, proteins related to hormonal signaling and their interaction with other pathways, and regulatory proteins were recorded indicating the involvement of five different signaling pathways. DAPs involved in the activation of genes and transcription factors of signaling and transduction pathways, transport and ion channels, cell-wall and structural proteins, defense, chaperons, energy metabolism, protein synthesis, modification, and turnover were identified. The gametophytes responded to changes in their micro-environment. There was also significant increase in prothallus biomass and conversion of two-dimensional prothalli into three-dimensional prothallus clumps at 3% sucrose. The three-D clumps had higher photosynthetic surface area and also closer proximity for sexual reproduction and sporophyte formation. Highest accumulation of proline, enhanced scavenging of reactive oxygen species (ROS) and DAPs of mostly, abiotic stress tolerance, secondary metabolite synthesis, and detoxification at 3% sucrose indicated an adaptive response of gametophytes. Protein Protein Interaction network and Principal Component analyses, and qRT-PCR validation of genes encoding 12 proteins of various metabolic processes indicated differential adjustment of gametophytes to different levels of sucrose in the culture medium. Therefore, a hypothetical mechanism was proposed to show that even slight changes in the micro-environment of D. maximum gametophytes triggered multiple mechanisms of adaptation. Many DAPs identified in the study have potential use in crop improvement and metabolic engineering programs, phytoremediation and environmental protection.

Spatially resolved transcriptomics reveals plant host responses to pathogens

BACKGROUND

Thorough understanding of complex model systems requires the characterisation of processes in different cell types of an organism. This can be achieved with high-throughput spatial transcriptomics at a large scale. However, for plant model systems this is still challenging as suitable transcriptomics methods are sparsely available. Here we present GaST-seq (Grid-assisted, Spatial Transcriptome sequencing), an easy to adopt, micro-scale spatial-transcriptomics workflow that allows to study expression profiles across small areas of plant tissue at a fraction of the cost of existing sequencing-based methods.

RESULTS

We compare the GaST-seq method with widely used library preparation methods (Illumina TruSeq). In spatial experiments we show that the GaST-seq method is sensitive enough to identify expression differences across a plant organ. We further assess the spatial transcriptome response of Arabidopsis thaliana leaves exposed to the bacterial molecule flagellin-22, and show that with eukaryotic (Albugo laibachii) infection both host and pathogen spatial transcriptomes are obtained.

CONCLUSION

We show that our method can be used to identify known, rapidly flagellin-22 elicited genes, plant immune response pathways to bacterial attack and spatial expression patterns of genes associated with these pathways.

De novo transcriptome provides insights into the growth behaviour and resveratrol and trans-stilbenes biosynthesis in Dactylorhiza hatagirea - An endangered alpine terrestrial orchid of western Himalaya

This is the first report on de novo transcriptome of Dactylorhiza hatagirea, a critically-endangered, terrestrial orchid of alpine Himalayas. The plant is acclaimed for medicinal properties but little is known about its secondary-metabolites profile or cues regulating their biosynthesis. De novo transcriptome analysis was therefore, undertaken to gain basic understanding on these aspects, while circumventing the acute limitation of plant material availability. 65,384 transcripts and finally, 37,371 unigenes were assembled de novo from a total of 236 million reads obtained from shoot, tuber and leaves of the plant. Dominance of differentially-expressing-genes (DEGs) related to cold-stress-response and plant-hormone-signal-transduction; and those involved in photosynthesis, sugar-metabolism and secondary-metabolite-synthesis provided insights into carbohydrate-partitioning in the plant during its preparation for freezing winter at natural habitat. DEGs of glucomannan, ascorbic acid, carotenoids, phylloquinone/naphthoquinones, indole alkaloids, resveratrol and stilbene biosynthesis revealed the secondary-metabolite profile of D. hatagirea. UHPLC results confirmed appreciable amounts of resveratrol and trans-stilbene in D. hatagirea tubers, for the first time. Expression analysis of 15 selected genes including those of phenylpropanoid pathway confirmed the validity of RNA-seq data. Opportunistic growth, temperature- and tissue-specific-differential-expression of secondary metabolite biosynthesis and stress tolerant genes were confirmed using clonal plants growing at 8, 15 and 25 °C.

Transcriptional profiling of contrasting genotypes revealed key candidates and nucleotide variations for drought dissection in Camellia sinensis (L.) O. Kuntze

Tea is popular health beverage consumed by millions of people worldwide. Drought is among the acute abiotic stress severely affecting tea cultivation, globally. In current study, transcriptome sequencing of four diverse tea genotypes with inherent contrasting genetic response to drought (tolerant & sensitive) generated more than 140 million reads. De novo and reference-based assembly and functional annotation of 67,093 transcripts with multifarious public protein databases yielded 54,484 (78.2%) transcripts with significant enrichment of GO and KEGG drought responsive pathways in tolerant genotypes. Comparative DGE and qRT analysis revealed key role of ABA dependent & independent pathways, potassium & ABC membrane transporters (AtABCG22, AtABCG11, AtABCC5 & AtABCC4) and antioxidant defence system against oxidative stress in tolerant genotypes, while seems to be failed in sensitive genotypes. Additionally, highly expressed UPL3HECT E3 ligases and RING E3 ligases possibly enhance drought tolerance by actively regulating functional modification of stress related genes. Further, ascertainment of, 80803 high quality putative SNPs with functional validation of key non-synonymous SNPs suggested their implications for developing high-throughput genotyping platform in tea. Futuristically, functionally relevant genomic resources can be potentially utilized for gene discovery, genetic engineering and marker-assisted genetic improvement for better yield and quality in tea under drought conditions.

Genome-wide identification and characterization of abiotic-stress responsive SOD (superoxide dismutase) gene family in Brassica juncea and B. rapa

Background

Abiotic stresses like drought, heat, cold and salinity cause major productivity loss in the rapeseed-mustard crops (Brassica). Major efforts have been made in the past to identify genes that provide resistance against such stresses. Superoxide dismutase (SOD) proteins, member of the metallo-enzyme family play vital role in protecting plants against abiotic stresses. In the present study, genome-wide analysis of abiotic stress responsive SOD gene family has been done in B. juncea and B. rapa.

Results

A total of 29 and 18 SOD genes were identified in B. juncea and B. rapa respectively and chromosome location mapping indicated their wide distribution across genome. On the basis of domain composition, the SODs were phylogenetically classified into sub-groups which was also substantiated by the gene structure and sub-cellular locations of SOD proteins. Functional annotation of SODs was also done by Gene Ontology (GO) mapping and the result was corroborated by the identified cis-regulatory elements in the promoter region of SOD genes. Based on FPKM analysis of SRA data available for drought, heat and salt stress, we identified 14 and 10 abiotic stress responsive SOD genes in B. rapa and B. juncea respectively. The differential expression analysis under drought and heat stress of identified abiotic-stress responsive SOD genes was done through quantitative Real Time PCR.

Conclusion

We identified abiotic-stress responsive genes that could help in improving the plant tolerance against abiotic stresses. This was the first study to describe the genome-wide analysis of SOD gene family in B. rapa and B. juncea, and the results will help in laying basic ground for future work of cloning and functional validation of SOD genes during abiotic stresses leading to Brassica crop improvement.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5593-5) contains supplementary material, which is available to authorized users.

Transcriptome analysis of the curry tree (Bergera koenigii L., Rutaceae) during leaf development

The curry tree (Bergera koenigii L.) is a widely cultivated plant used in South Asian cooking. Next-generation sequencing was used to generate the transcriptome of the curry leaf to detect changes in gene expression during leaf development, such as those genes involved in the production of oils which lend the leaf its characteristic taste, aroma, and medicinal properties. Using abundance estimation (RSEM) and differential expression analysis, genes that were significantly differentially expressed were identified. The transcriptome was annotated with BLASTx using the non-redundant (nr) protein database, and Gene Ontology (GO) terms were assigned based on the top BLAST hit using Blast2GO. Lastly, functional enrichment of the assigned GO terms was analyzed for genes that were significantly differentially expressed. Of the most enriched GO categories, pathways involved in cell wall, membrane, and lignin synthesis were found to be most upregulated in immature leaf tissue, possibly due to the growth and expansion of the leaf tissue. Terpene synthases, which synthesize monoterpenes and sesquiterpenes, which comprise much of the curry essential oil, were found to be significantly upregulated in mature leaf tissue, suggesting that oil production increases later in leaf development. Enzymes involved in pigment production were also significantly upregulated in mature leaves. The findings were based on computational estimates of gene expression from RNA-seq data, and further study is warranted to validate these results using targeted techniques, such as quantitative PCR.

Efficient and High-Quality RNA Isolation from Metabolite-Rich Tissues of Stevia rebaudiana, an Important Commercial Crop

Stevia rebaudiana, a perennial herb native to northeastern Paraguay, has gained immense attention globally over the recent decades due to the natural sweetness of its leaves. Like in most plants, this particular species contains high amount of secondary metabolites, thus rendering the isolation of high quality and quantity RNA extract for molecular applications rather challenging. An effective, high-yield and high-quality RNA isolation protocol for this economically important plant species was devised here based on the cetyltrimethylammonium bromide (CTAB) extraction method, with an additional genomic DNA (gDNA) removal step. DNA and other contaminants that may affect downstream applications were effectively removed. Our results exhibited that RNA samples isolated from the leaves and stems of Stevia rebaudiana using this improvised method are high in integrity and quality with RNA integrity number (RIN) of more than 8 and low in contaminants.

Long Non-coding RNAs Coordinate Developmental Transitions and Other Key Biological Processes in Grapevine

Long non-coding RNAs (lncRNAs) are transcripts >200 nucleotides that have prominently surfaced as dynamic regulatory molecules. Using computational approaches, we identified and characterized 56,441 lncRNAs in grapevine (Vitis vinifera) by harnessing RNA-seq data from 10 developmental stages of leaf, inflorescence, and berry tissues. We conducted differential expression analysis and determined tissue- and developmental stage-specificity of lncRNAs in grapevine, which indicated their spatiotemporal regulation. Functional annotation using co-expression analysis revealed their involvement in regulation of developmental transitions in sync with transcription factors (TFs). Further, pathway enrichment analysis revealed lncRNAs associated with biosynthetic and secondary metabolic pathways. Additionally, we identified 115, 560, and 133 lncRNAs as putative miRNA precursors, targets, and endogenous target mimics, respectively, which provided an insight into the interplay of regulatory RNAs. We also explored lncRNA-mediated regulation of extra-chromosomal genes–i.e., mitochondrial and chloroplast coding sequences and observed their involvement in key biological processes like ‘photosynthesis’ and ‘oxidative phosphorylation’. In brief, these transcripts coordinate important biological functions via interactions with both coding and non-coding RNAs as well as TFs in grapevine. Our study would facilitate future experiments in unraveling regulatory mechanisms of development in this fruit crop of economic importance.

An optimized isolation protocol yields high-quality RNA from cassava tissues (Manihot esculenta Crantz)

We developed and modified a precise, rapid, and reproducible protocol isolating high-quality RNA from tissues of multiple varieties of cassava plants (Manihot esculenta Crantz). The resulting method is suitable for use in mini, midi, and maxi preparations and rapidly achieves high total RNA yields (170-600 μg·g^-1) using low-cost chemicals and consumables and with minimal contamination from polysaccharides, polyphenols, proteins, and other secondary metabolites. In particular, A ₂₆₀ : A ₂₈₀ ratios were > 2.0 for RNA from various tissues, and all of the present RNA samples yielded ribosomal integrity number values of greater than six. The resulting high purity and quality of isolated RNA will facilitate downstream applications (quantitative reverse transcriptase-polymerase chain reaction or RNA sequencing) in cassava molecular breeding.

Ectopic expression of SOD and APX genes in Arabidopsis alters metabolic pools and genes related to secondary cell wall cellulose biosynthesis and improve salt tolerance

Hydrogen peroxide (H₂O₂) is known to accumulate in plants during abiotic stress conditions and also acts as a signalling molecule. In this study, Arabidopsis thaliana transgenics overexpressing cytosolic CuZn-superoxide dismutase (PaSOD) from poly-extremophile high-altitude Himalayan plant Potentilla atrosanguinea, cytosolic ascorbate peroxidase (RaAPX) from Rheum australe and dual transgenics overexpressing both the genes were developed and analyzed under salt stress. In comparison to wild-type (WT) or single transgenics, the performance of dual transgenics under salt stress was better with higher biomass accumulation and cellulose content. We identified genes involved in cell wall biosynthesis, including nine cellulose synthases (CesA), seven cellulose synthase-like proteins together with other wall-related genes. RNA-seq analysis and qPCR revealed differential regulation of genes (CesA 4, 7 and 8) and transcription factors (MYB46 and 83) involved in secondary cell wall cellulose biosynthesis, amongst which most of the cellulose biosynthesis gene showed upregulation in single (PaSOD line) and dual transgenics at 100 mM salt stress. A positive correlation between cellulose content and H₂O₂ accumulation was observed in these transgenic lines. Further, cellulose content was 1.6-2 folds significantly higher in PaSOD and dual transgenic lines, 1.4 fold higher in RaAPX lines as compared to WT plants under stress conditions. Additionally, transgenics overexpressing PaSOD and RaAPX also displayed higher amounts of phenolics as compared to WT. The novelty of present study is that H₂O₂ apart from its role in signalling, it also provides mechanical strength to plants and aid in plant biomass production during salt stress by transcriptional activation of cellulose biosynthesis pathway. This modulation of the cellulose biosynthetic machinery in plants has the potential to provide insight into plant growth, morphogenesis and to create plants with enhanced cellulose content for biofuel use.

Global Transcriptional Insights of Pollen-Pistil Interactions Commencing Self-Incompatibility and Fertilization in Tea [Camellia sinensis (L.) O. Kuntze]

This study explicates molecular insights commencing Self-Incompatibility (SI) and CC (cross-compatibility/fertilization) in self (SP) and cross (CP) pollinated pistils of tea. The fluorescence microscopy analysis revealed ceased/deviated pollen tubes in SP, while successful fertilization occurred in CP at 48 HAP. Global transcriptome sequencing of SP and CP pistils generated 109.7 million reads with overall 77.9% mapping rate to draft tea genome. Furthermore, concatenated de novo assembly resulted into 48,163 transcripts. Functional annotations and enrichment analysis (KEGG & GO) resulted into 3793 differentially expressed genes (DEGs). Among these, de novo and reference-based expression analysis identified 195 DEGs involved in pollen-pistil interaction. Interestingly, the presence of 182 genes [PT germination & elongation (67), S-locus (11), fertilization (43), disease resistance protein (30) and abscission (31)] in a major hub of the protein-protein interactome network suggests a complex signaling cascade commencing SI/CC. Furthermore, tissue-specific qRT-PCR analysis affirmed the localized expression of 42 DE putative key candidates in stigma-style and ovary, and suggested that LSI initiated in style and was sustained up to ovary with the active involvement of csRNS, SRKs & SKIPs during SP. Nonetheless, COBL10, RALF, FERONIA-rlk, LLG and MAPKs were possibly facilitating fertilization. The current study comprehensively unravels molecular insights of phase-specific pollen-pistil interaction during SI and fertilization, which can be utilized to enhance breeding efficiency and genetic improvement in tea.

Transcriptome-wide mining, characterization, and development of microsatellite markers in Lychnis kiusiana (Caryophyllaceae)

BACKGROUND

Lychnis kiusiana Makino is an endangered perennial herb native to wetland areas in Korea and Japan. Despite its conservational and evolutionary significance, population genetic resources are lacking for this species. Next-generation sequencing has been accepted as a rapid and cost-effective solution for the identification of microsatellite markers in nonmodel plants.

RESULTS

Using Illumina HiSeq 2000 sequencing technology, we assembled 67,498,600 reads into 91,900 contigs and identified 11,403 microsatellite repeat motifs in 9563 contigs. A total of 4510 microsatellite-containing transcripts had Gene Ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 124 pathways with significant scores. Many microsatellites in the L. kiusiana leaf transcriptome were linked to genes involved in the plant response to light intensity, salt stress, temperature stimulus, and nutrient and water deprivation. A total of 12,486 single-nucleotide polymorphisms (SNPs) were identified on transcripts harboring microsatellites. The analysis of nucleotide substitution rates for 2389 unigenes indicated that 39 genes were under strong positive selection. The primers of 6911 microsatellites were designed, and 40 of 50 selected primer pairs were consistently and successfully amplified from 51 individuals. Twenty-five of these were polymorphic, and the average number of alleles per SSR locus was 6.96, with a range from 2 to 15. The observed and expected heterozygosities ranged from 0.137 to 0.902 and 0.131 to 0.827, respectively, and locus-specific F_IS estimates ranged from - 0.116 to 0.290. Eleven of the 25 primer pairs were successfully amplified in three additional species of Lychnis: 56% in L. wilfordii, 64% in L. cognata and 80% in L. fulgens.

CONCLUSIONS

The transcriptomic SSR markers of Lychnis kiusiana provide a valuable resource for understanding the population genetics, evolutionary history, and effective conservation management of this species. Furthermore, the identified microsatellite loci linked to the annotated genes should be useful for developing functional markers of L. kiusiana. The developed markers represent a potentially valuable source of transcriptomic SSR markers for population genetic analyses with moderate levels of cross-taxon portability.

Validation of reference genes for qRT-PCR data normalisation in lentil (Lens culinaris) under leaf developmental stages and abiotic stresses

Lentil (Lens culinaris) is one of the most important staple food crops of developing countries. Transcriptome based global gene expression profiling followed by validation of expression of important genes through quantitative real time-PCR (qRT-PCR) has achieved significance in recent years. However, there is a severe scarcity of information regarding stable reference genes in lentil, which is mandatory for qRT-PCR data normalisation. Hence, the present study was under-taken to identify the most stable reference gene(s) in lentil. Expression stability of eight candidate genes viz. ribulose 1,5-bisphosphate carboxylase large subunit (Rbcl), ribosomal protein L2 (RPL2), 18S rRNA, tubulin (Tub), elongation factor 1α (EF1α), glyceraldehydes-3-phosphate dehydrogenase (GAPDH), heat shock protein (HSP70), and Maturase (mat K) was evaluated in five varieties of lentil at three different stages of leaf development and abiotic stress conditions using qRT-PCR. The results were analysed using four types of statistical software viz., geNorm, BestKeeper, NormFinder and RefFinder; all softwares identified RPL2 as most stable under abiotic stress conditions and developmental stages followed by Tub and Rbcl; while, HSP70 was identified as least stable. Relative expression of the target genes, defensin and PR4, was evaluated under abiotic stress conditions and data normalisation was done using two stable reference genes, RPL2 and Tub, either alone or in combination and with two least stable genes, HSP70 and 18S. The present work provides a list of potential reference genes in lentil, which will help in selection of appropriate reference gene for qRT-PCR data normalization depending upon the experiment.

De Novo Transcriptomic Analysis and Development of EST–SSRs for Styrax japonicus

Styrax japonicus sieb. et Zucc. is widely distributed in China with ornamental and medicinal values. However, the transcriptome of S. japonicus has not yet been reported. In this study, we carried out the first transcriptome analysis of S. japonicus and developed a set of expressed sequence tag–simple sequence repeats (EST–SSRs). We obtained 338,570,222 clean reads in total, of which the mean GC content was 41.58%. In total, 136,071 unigenes were obtained having an average length of 611 bp and 71,226 unigenes were favorably annotated in the database. In total, we identified 55,977 potential EST–SSRs from 38,611 unigenes, of which there was 1 SSR per 6.73 kb. The di-nucleotide repeats (40.40%) were the most identified SSRs. One set of 60 primer pairs was randomly selected, and the amplified products in S. japonicus were validated; 28 primer pairs successfully produced clear amplicons. A total of 21 (35%) polymorphic genic SSR markers were identified between two populations. In total, 15 alleles were detected and the average number was 6. The average of observed heterozygosity and expected heterozygosity was 0.614 and 0.552, respectively. The polymorphism information content (PIC) value fluctuated between 0.074 and 0.855, with a mean value of 0.504, which was also the middle level. This study provides useful information for diversity studies and resource assessments of S. japonicus.

Elucidating genes involved in sesquiterpenoid and flavonoid biosynthetic pathways in Saussurea lappa by de novo leaf transcriptome analysis

Saussurea lappa (family Asteraceae) possesses immense pharmacological potential mainly due to the presence of sesquiterpene lactones. In spite of its medicinal importance, S. lappa has been poorly explored at the molecular level. We initiated leaf transcriptome sequencing of S. lappa using the illumina highseq 2000 platform and generated 62,039,614 raw reads. Trinity assembler generated 122,434 contigs with an N50 value of 1053 bp. The assembled transcripts were compared against the non-redundant protein database at NCBI. The Blast2GO analysis assigned gene ontology (GO) terms, categorized into molecular functions (3132), biological processes (4477) and cellular components (1.927). Using KEGG, around 476 contigs were assigned to 39 pathways. For secondary metabolic pathways, we identified transcripts encoding genes involved in sesquiterpenoid and flavonoid biosynthesis. Relatively low number of transcripts were also found encoding for genes involved in the alkaloid pathway. Our data will contribute to functional genomics and metabolic engineering studies in this plant.

Reconfiguration of the plastid genome in Lamprocapnos spectabilis: IR boundary shifting, inversion, and intraspecific variation

We generated a complete plastid genome (plastome) sequence for Lamprocapnos spectabilis, providing the first complete plastome from the subfamily Fumarioideae (Papaveraceae). The Lamprocapnos plastome shows large differences in size, structure, gene content, and substitution rates compared with two sequenced Papaveraceae plastomes. We propose a model that explains the major rearrangements observed, involving at least six inverted repeat (IR) boundary shifts and five inversions, generating a number of gene duplications and relocations, as well as a two-fold expansion of the IR and miniaturized small single-copy region. A reduction in the substitution rates for genes transferred from the single-copy regions to the IR was observed. Accelerated substitution rates of plastid accD and clpP were detected in the Lamprocapnos plastome. The accelerated substitution rate for the accD gene was correlated with a large insertion of amino acid repeat (AAR) motifs in the middle region, but the forces driving the higher substitution rate of the clpP gene are unclear. We found a variable number of AARs in Lamprocapnos accD and ycf1 genes within individuals, and the repeats were associated with coiled-coil regions. In addition, comparative analysis of three Papaveraceae plastomes revealed loss of rps15 in Papaver, and functional replacement to the nucleus was identified.

Heterologous Expression of Key C and N Metabolic Enzymes Improves Re-assimilation of Photorespired CO2 and NH3, and Growth

We investigated the effect of the heterologous expression of phosphoenolpyruvate carboxylase (ZmPepcase), aspartate aminotransferase (GmAspAT), and glutamine synthetase (NtGS) on carbon (C) and nitrogen (N) metabolism in Arabidopsis (Arabidopsis thaliana). These transgenes were expressed either separately or in different combinations. The highest gains in shoot dry weight were observed in transgenic lines coexpressing all three genes. Tracer experiments using NaH14CO3 suggested that the coexpression of ZmPepcase, GmAspAT, and NtGS resulted in a higher flux of assimilated CO2 toward sugars and amino acids. Upon feeding the leaf discs with glycine-1-14C, transgenic lines evolved significantly lower 14CO2 levels than the wild type, suggesting that a higher reassimilation of CO2 evolved during photorespiration. Leaves of transgenic plants accumulated significantly lower ammonium without any significant difference in the levels of photorespiratory ammonium relative to the wild type, suggesting a higher reassimilation of photorespired NH3 Transgenic lines also showed improved photosynthetic rates, higher shoot biomass accumulation, and improved seed yield in comparison with wild-type plants under both optimum and limiting N conditions. This work demonstrates that the heterologous coexpression of ZmPepcase, GmAspAT, and NtGS reduced the photorespiratory loss of C and N with concomitant enhancements in shoot biomass and seed yield.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

Extraction of RNA from Recalcitrant Tree Species Paulownia elongata

Isolation of pure RNA is the basic requisite for most molecular biology work. Plants contain polyphenols and polysaccharides, which can interfere with isolation of pure RNA from them. Especially hardwood tree species like Paulownia elongata have surplus amount of RNA-binding alkaloids, proteins and secondary metabolites that can further complicate the process of RNA extraction. Paulownia elongata is a fast-growing tree species which is known for its role in environmental adaptability and biofuel research. Here we describe an economical, efficient and time-saving method (2 h) to extract RNA from leaf tissues of the tree Paulownia elongata. Lack of DNA contamination and good RNA integrity were confirmed using RNA Gel electrophoresis. The purity of RNA was confirmed using Nanodrop spectrophotometer that revealed an A₂₆₀:A₂₈₀ ratio of about 2.0. The purified RNA was successfully used in the downstream applications such as RT-PCR (Reverse Transcription PCR) and qPCR (quantitative PCR). This method could be used for RNA extraction from several other recalcitrant tree species.

Biofortification of safflower: an oil seed crop engineered for ALA-targeting better sustainability and plant based omega-3 fatty acids

Alpha-linolenic acid (ALA) deficiency and a skewed n6:n3 fatty acid ratio in the diet is a major explanation for the prevalence of cardiovascular diseases and inflammatory/autoimmune diseases. There is mounting evidence of the health benefits associated with omega-3 long chain polyunsaturated fatty acids (LC PUFA's). Although present in abundance in fish, a number of factors limit our consumption of fish based omega-3 PUFA's. To name a few, overexploitation of wild fish stocks has reduced their sustainability due to increased demand of aquaculture for fish oil and meal; the pollution of marine food webs has raised concerns over the ingestion of toxic substances such as heavy metals and dioxins; vegetarians do not consider fish-based sources for supplemental nutrition. Thus alternative sources are being sought and one approach to the sustainable supply of LC-PUFAs is the metabolic engineering of transgenic plants with the capacity to synthesize n3 LC-PUFAs. The present investigation was carried out with the goal of developing transgenic safflower capable of producing pharmaceutically important alpha-linolenic acid (ALA, C18:3, n3). This crop was selected as the seeds accumulate ~ 78% of the total fatty acids as linoleic acid (LA, C18:2, n6), the immediate precursor of ALA. In the present work, ALA production was achieved successfully in safflower seeds by transforming safflower hypocotyls with Arabidopsis specific delta 15 desaturase (FAD3) driven by truncated seed specific promoter. Transgenic safflower fortified with ALA is not only potentially valuable nutritional superior novel oil but also has reduced ratio of LA to ALA which is required for good health.

A method for extracting high-quality total RNA from plant rich in polysaccharides and polyphenols using Dendrobium huoshanense

Acquiring high quality RNA is the basis of plant molecular biology research, plant genetics and other physiological investigations. At present, a large number of nucleotide isolation methods have been exploited or modified, such as commercial kits, CTAB, SDS methods and so on. Due to the nature of different plants, extraction methods vary. Moreover, efficiency of certain approach cannot be guaranteed due to composition of different plants and extracting high quality RNA from plants rich in polysaccharides and polyphenols are often difficult. The physical and chemical properties of polysaccharides which are similar to nucleic acids and other secondary metabolites will be coprecipitated with RNA irreversibly. Therefore, how to remove polysaccharides and other secondary metabolites during RNA extraction is the primary challenge. Dendrobium huoshanense is an Orchidaceae perennial herb that is rich in polysaccharides and other secondary metabolites. By using D. huoshanense as the subject, we improved the method originated from CHAN and made it suitable for plants containing high amount of polysaccharides and polyphenols. The extracted total RNA was clear and non-dispersive, with good integrity and no obvious contamination with DNA and other impurities. And it was also evaluated by gel electrophoresis, nucleic acid quantitative detector and PCR assessment. Thus, as a simple approach, it is suitable and efficient in RNA isolation for plants rich in polysaccharides and polyphenols.

AV2 protein of tomato leaf curl Palampur virus promotes systemic necrosis in Nicotiana benthamiana and interacts with host Catalase2

Tomato leaf curl Palampur virus (ToLCPalV) is a whitefly-transmitted, bipartite begomovirus. Here, we demonstrated that ectopic expression of AV2 from a Potato virus X (PVX)-based vector accelerated systemic necrosis and reactive oxygen species (ROS) accumulation in Nicotiana benthamiana. Furthermore, 10 amino acids from N-terminal region of AV2 were found to be associated with the systemic necrosis symptom/phenotype. Mutational studies of ToLCPalV infectious clones lacking the AV2 revealed that AV2 is essential for the systemic movement of DNA-A, symptom severity and viral DNA accumulation. In a yeast two-hybrid assay, Catalase2 (Cat2) was found to associate with AV2 protein. Further, silencing of Cat2 resulted in appearance of necrotic lesions on N. benthamiana and these plants were highly susceptible to ToLCPalV infection in comparison to control plants. Infection ToLCPalV on Solanum lycopersicum resulted in downregulation of Cat2 transcripts, followed by accumulation of ROS and stress marker transcripts. The AV2 protein also suppressed virus-induced gene silencing (VIGS) of the Phytoene desaturase (PDS) gene. Our results show that AV2 is essential for the pathogenicity, systemic movement and suppression of gene silencing in the host. Altogether, our findings suggest that interactions between AV2 and Cat2 might play a crucial role in the establishment of ToLCPalV infection.

Small RNA profiling for identification of miRNAs involved in regulation of saponins biosynthesis in Chlorophytum borivilianum

BACKGROUND

MicroRNAs act as molecular regulator of cell signaling, plant growth and development, and regulate various primary and secondary plant metabolic processes. In the present study, deep sequencing of small RNAs was carried out to identify known and novel miRNAs from pharmaceutically important plant, Chlorophytum borivilianum.

RESULTS

Total 442 known miRNAs and 5 novel miRNAs were identified from young leaf small RNA library. Experimental validation with stem loop RT-PCR confirmed the in silico identification. Based on transcriptome data of root and leaf of C. borivilianum, Oryza sativa, and Arabidopsis thaliana target gene prediction was done using psRNAtarget and mirRanda. BLAST2GO helped in localization of predicted targets and KEGG (Kyoto Encyclopedia for Genes and Genomes) pathway analysis concluded that miR9662, miR894, miR172, and miR166 might be involved in regulating saponin biosynthetic pathway. The correlation between miRNA and its target gene was further validated by RT-qPCR analysis.

CONCLUSION

This study provides first elaborated glimpse of miRNA pool of C. borivilianum, which can help to understand the miRNA dependent regulation of saponin biosynthesis and to design further metabolic engineering experiment to enhance their contents in the plant.

An Improved Total RNA Extraction Method for White Jelly Mushroom Tremella fuciformis Rich in Polysaccharides

An improved method for extracting high quality and quantity RNA from a jelly mushroom and a dimorphic fungus-Tremella fuciformis which is especially rich in polysaccharides, is described. RNA was extracted from T. fuciformis mycelium M1332 and its parental monokaryotic yeast-like cells Y13 and Y32. The A260/280 and A260/230 ratios were both approximately 2, and the RNA integrity number was larger than 8.9. The yields of RNA were between 108 and 213 µg/g fresh wt. Downstream molecular applications including reverse transcriptional PCR and quantitative real-time PCR were also performed. This protocol is reliable and may be widely applicable for total RNA extraction from other jelly mushrooms or filamentous fungi rich in polysaccharides.

Molecular dissection of transcriptional reprogramming of steviol glycosides synthesis in leaf tissue during developmental phase transitions in Stevia rebaudiana Bert

Stevia is a natural source of commercially important steviol glycosides (SGs), which share biosynthesis route with gibberellic acids (GAs) through plastidal MEP and cytosolic MVA pathways. Ontogeny-dependent deviation in SGs biosynthesis is one of the key factor for global cultivation of Stevia, has not been studied at transcriptional level. To dissect underlying molecular mechanism, we followed a global transcriptome sequencing approach and generated more than 100 million reads. Annotation of 41,262 de novo assembled transcripts identified all the genes required for SGs and GAs biosynthesis. Differential gene expression and quantitative analysis of important pathway genes (DXS, HMGR, KA13H) and gene regulators (WRKY, MYB, NAC TFs) indicated developmental phase dependent utilization of metabolic flux between SGs and GAs synthesis. Further, identification of 124 CYPs and 45 UGTs enrich the genomic resources, and their PPI network analysis with SGs/GAs biosynthesis proteins identifies putative candidates involved in metabolic changes, as supported by their developmental phase-dependent expression. These putative targets can expedite molecular breeding and genetic engineering efforts to enhance SGs content, biomass and yield. Futuristically, the generated dataset will be a useful resource for development of functional molecular markers for diversity characterization, genome mapping and evolutionary studies in Stevia.

Molecular and functional characterization of two isoforms of chalcone synthase and their expression analysis in relation to flavonoid constituents in Grewia asiatica L

Chalcone synthase constitutes a functionally diverse gene family producing wide range of flavonoids by catalyzing the initial step of the phenylpropanoid pathway. There is a pivotal role of flavonoids in pollen function as they are imperative for pollen maturation and pollen tube growth during sexual reproduction in flowering plants. Here we focused on medicinally important fruit-bearing shrub Grewia asiatica. It is a rich repository of flavonoids. The fruits are highly acclaimed for various putative health benefits. Despite its importance, full commercial exploitation is hampered due to two drawbacks which include short shelf life of its fruits and larger seed volume. To circumvent these constraints, seed abortion is one of the viable options. Molecular interventions tested in a number of economic crops have been to impair male reproductive function by disrupting the chalcone synthase (CHS) gene activity. Against this backdrop the aim of the present study included cloning and characterization of two full-length cDNA clones of GaCHS isoforms from the CHS multigene family. These included GaCHS1 (NCBI acc. KX129910) and GaCHS2 (NCBI acc. KX129911) with an ORF of 1176 and 1170 bp, respectively. GaCHSs were heterologously expressed and purified in E. coli to validate their functionality. Functionality of CHS isoforms was also characterized via enzyme kinetic studies using five different substrates. We observed differential substrate specificities in terms of their Km and Vmax values. Accumulation of flavonoid constituents naringenin and quercetin were also quantified and their relative concentrations corroborated well with the expression levels of GaCHSs. Further, our results demonstrate that GaCHS isoforms show differential expression patterns at different reproductive phenological stages. Transcript levels of GaCHS2 were more than its isoform GaCHS1 at the anthesis stage of flower development pointing towards its probable role in male reproductive maturity.