An RNA isolation system for plant tissues rich in secondary metabolites

Host Response of Arabidopsis thaliana Interaction with Fungal Endophytes Involves microRNAs

Plant and fungus interaction is a complex process involving many molecular factors determining the nature of relationship. The enigmatic methodology by which fungal endophytes are able to colonise a plant harmoniously is still inexplicable. Small RNAs have been identified as major regulatory elements under various biotic interactions. However, their role in endophytic plant-fungal interactions remain to be elucidated. Therefore, transcript expression data available on Gene Expression Omnibus for Arabidopsis thaliana was utilised for miRNAs identification under endophytism. The analysis predicted 15 miRNAs with differential expression of which the ath-miRNA398b modulation was significant. Application of psRNAtarget, C-mii, pmiREN, and TarDB provided a pool of 357 target genes for these miRNAs. Protein-protein interaction analysis identified major hub proteins, including BTB/POZ domain-containing protein, beta-Xylosidase-2 (AtBXL2), and Copper/Zinc Superoxide Dismutase-2 (AtSOD2). The quantitative real-time PCR validated the computational prediction and expression for selected target genes AtSOD2, AtBXL2, and AtRCA along with ath-miRNA398b under endophytism. Overall, results indicate that miRNAs have a significant role in regulating Arabidopsis thaliana-endophytic fungal interaction.

High-quality haplotype-resolved chromosome assembly provides evolutionary insights and targeted steviol glycosides (SGs) biosynthesis in Stevia rebaudiana Bertoni

Stevia rebaudiana Bertoni is popular source of plant-derived low/no-calorie natural sweeteners (LNCSs), collectively known as steviol glycosides (SGs). Nevertheless, genetic predisposition for targeted biosynthesis of SGs is complex due to multi-substrate functionality of key uridine diphosphate glycosyltransferases (UGTs). Here, we created a high-quality monoploid assembly of 1.34 Gb with N50 value of 110 Mb, 55 551 predicted protein-coding genes, and ~80% repetitive regions in Rebaudioside-A (Reb-A) enriched cultivar of S. rebaudiana. Additionally, a haplotype-based chromosome assembly consisting of haplotype A and haplotype B with an overall genome size of 2.33Gb was resolved, harbouring 639 634 variants including single nucleotide polymorphisms (SNPs), indels and structural variants (SVs). Furthermore, a lineage-specific whole genome duplication analysis revealed that gene families encoding UGTs and Cytochrome-P450 (CYPs) were tandemly duplicated. Additionally, expression analysis revealed five tandemly duplicated gene copies of UGT76G1 having significant correlations with Reb-A content, and identified key residue (leu200val) in the glycosylation of Reb-A. Furthermore, missense variations identified in the acceptor region of UGT76G1 in haplotype resolve genome, transcriptional and molecular docking analysis were confirmed with resequencing of 10 diverse stevia genotypes (~25X). Gene regulatory network analysis identified key transcription factors (MYB, bHLH, bZIP and AP2-ERF) as potential regulators of SG biosynthesis. Overall, this study provides haplotype-resolved chromosome-level genome assembly for genome editing and enhancing breeding efforts for targeted biosynthesis of SGs in S. rebaudiana.

MYB transcription factors, their regulation and interactions with non-coding RNAs during drought stress in Brassica juncea

BACKGROUND

Brassica juncea (L.) Czern is an important oilseed crop affected by various abiotic stresses like drought, heat, and salt. These stresses have detrimental effects on the crop's overall growth, development and yield. Various Transcription factors (TFs) are involved in regulation of plant stress response by modulating expression of stress-responsive genes. The myeloblastosis (MYB) TFs is one of the largest families of TFs associated with various developmental and biological processes such as plant growth, secondary metabolism, stress response etc. However, MYB TFs and their regulation by non-coding RNAs (ncRNAs) in response to stress have not been studied in B. juncea. Thus, we performed a detailed study on the MYB TF family and their interactions with miRNAs and Long non coding RNAs.

RESULTS

Computational investigation of genome and proteome data presented a comprehensive picture of the MYB genes and their protein architecture, including intron-exon organisation, conserved motif analysis, R2R3 MYB DNA-binding domains analysis, sub-cellular localization, protein-protein interaction and chromosomal locations. Phylogenetically, BjuMYBs were further classified into different subclades on the basis of topology and classification in Arabidopsis. A total of 751 MYBs were identified in B. juncea corresponding to 297 1R-BjuMYBs, 440 R2R3-BjuMYBs, 12 3R-BjuMYBs, and 2 4R-BjuMYBs types. We validated the transcriptional profiles of nine selected BjuMYBs under drought stress through RT-qPCR. Promoter analysis indicated the presence of drought-responsive cis-regulatory components. Additionally, the miRNA-MYB TF interactions was also studied, and most of the microRNAs (miRNAs) that target BjuMYBs were involved in abiotic stress response and developmental processes. Regulatory network analysis and expression patterns of lncRNA-miRNA-MYB indicated that selected long non-coding RNAs (lncRNAs) acted as strong endogenous target mimics (eTMs) of the miRNAs regulated expression of BjuMYBs under drought stress.

CONCLUSIONS

The present study has established preliminary groundwork of MYB TFs and their interaction with ncRNAs in B. juncea and it will help in developing drought- tolerant Brassica crops.

Quantitative and spatially resolved detection of multiplexed microRNA from plant tissue via hybridization to hydrogel-bound DNA probes in nanoliter well arrays

Understanding complex regulatory networks in plant systems requires elucidating the roles of various gene regulators under a spatial landscape. MicroRNA are key regulators that impart high information value through their tissue specificity and stability when using expression patterns for evaluating network outcomes. However, current techniques that utilize spatial multiplexing and quantitation of microRNA are limited to primarily mammalian systems. Here, we present a method to spatially resolve and quantify multiple endogenous microRNA in situ using ethanol fixed, paraffin embedded model plant species. This method utilizes target-specific microRNA capture along with universal ligating and labelling, all within functionalized hydrogel posts containing DNA probes in nanoliter well arrays. We demonstrate the platform's multiplexing capabilities through analyzing three endogenous microRNA in Arabidopsis thaliana rosettes which provide useful answers to fundamental plant growth and development from the unique expression patterns. The spatial tissue technique is also validated using non-spatial small RNA assays to demonstrate the versatility of the well array platform. Our new platform expands the toolkit of spatial omics technologies for plants.

Optimizing total RNA extraction method for human and mice samples

Background

Extracting high-quality total RNA is pivotal for advanced RNA molecular studies, such as Next-generation sequencing and expression microarrays where RNA is hybridized. Despite the development of numerous extraction methods in recent decades, like the cetyl-trimethyl ammonium bromide (CTAB) and the traditional TRIzol reagent methods, their complexity and high costs often impede their application in small-scale laboratories. Therefore, a practical and economical method for RNA extraction that maintains high standards of efficiency and quality needs to be provided to optimize RNA extraction from human and mice tissues.

Method

This study proposes enhancements to the TRIzol method by incorporating guanidine isothiocyanate (GITC-T method) and sodium dodecyl sulfate (SDS-T method). We evaluated the effectiveness of these modified methods compared to the TRIzol method using a micro-volume UV-visible spectrophotometer, electrophoresis, q-PCR, RNA-Seq, and whole transcriptome sequencing.

Result

The micro-volume UV-visible spectrophotometer, electrophoresis, and RNA-Seq demonstrated that the GITC-T method yielded RNA with higher yields, integrity, and purity, while the consistency in RNA quality between the two methods was confirmed. Taking mouse cerebral cortex tissue as a sample, the yield of total RNA extracted by the GITC-T method was 1,959.06 ± 49.68 ng/mg, while the yield of total RNA extracted by the TRIzol method was 1,673.08 ± 86.39 ng/mg. At the same time, the OD260/280 of the total RNA samples extracted by the GITC-T method was 2.03 ± 0.012, and the OD260/230 was 2.17 ± 0.031, while the OD260/280 of the total RNA samples extracted by the TRIzol method was 2.013 ± 0.041 and the OD260/230 was 2.11 ± 0.062. Furthermore, q-PCR indicated that the GITC-T method achieved higher yields, purity, and greater transcript abundance of total RNA from the same types of animal samples than the TRIzol method.

Conclusion

The GITC-T method not only yields higher purity and quantity of RNA but also reduces reagent consumption and overall costs, thereby presenting a more feasible option for small-scale laboratory settings.

CTAB Protocol for Obtaining High-Quality Total RNA from Hibiscus rosa-sinensis and Other Related Plants of the Family Malvaceae

Malvaceae family, also known as the Mallow family, is a family of flowering plants containing Hibiscus rosa-sinensis and other plants of high medicinal value. This study focuses on the challenges associated with high-quality RNA extraction from Hibiscus rosa-sinensis and its related plants characterized by high levels of mucilage and phenolic compounds in their tissues. High mucilage and secondary metabolite content pose obstacles in obtaining high-quality RNA, negatively impacting downstream applications, such as gene expression analysis. Our research aimed to develop an efficient RNA extraction method tailored to the unique characteristics of Malvaceae family plants especially Hibiscus rosa-sinensis. Through the substitution of NaCl with KCl, a crucial component of the CTAB buffer, our methodology successfully addressed the challenges posed by high mucilage and phenolic compound levels. This modification led to a significant reduction in sample viscosity, which is because of the high mucilage in these plants. Our modified CTAB extraction method yielded significantly more RNA with higher purity than the conventional CTAB methods alone. The extracted RNA was largely intact, as indicated by 28S/18S ratios and RIN values, yielding high-quality RNA with improved purity suggested by the 260/280 and 260/230 ratios. The proposed approach not only serves as a solution to the specific challenges encountered in Hibiscus rosa-sinensis but also holds promise for broader applications across different plants within the family.

A modified CTAB method for the extraction of high-quality RNA from mono-and dicotyledonous plants rich in secondary metabolites

BACKGROUND

High-quality RNA extraction from woody plants is difficult because of the presence of polysaccharides and polyphenolics that bind or co-precipitate with the RNA. The CTAB (cetyl trimethylammonium bromide) based method is widely used for the isolation of nucleic acids from polysaccharide-rich plants. Despite the widespread use of the CTAB method, it is necessary to adapt it to particular plant species, tissues and organs. Here we described a simple and generalized method for RNA isolation from mature leaf tissues of several economically important woody (17) and herbaceous plants (2) rich in secondary metabolites. High yields were achieved from small amount (up to 50 mg) of plant material. Two main modifications were applied to the basic protocol: an increase in β-mercaptoethanol concentration (to 10%v/v) and the use of an effective DNase treatment. As opposed to similar studies, we tried to describe a more detailed protocol for isolating RNA, including the exact quantity and concentration of the reagents were used.

RESULTS

Our modified CTAB method is proved to be efficient in extracting the total RNA from a broad range of woody and herbaceous species. The RNA yield was ranged from 2.37 to 91.33 µg/µl. The A260:A280 and A260:A230 absorbance ratios were measured from 1.77 to 2.13 and from 1.81 to 2.22. The RIN value (RNA Integrity Number) of the samples fell between 7.1 and 8.1, which indicated that a small degree of RNA degradation occurred during extraction. The presence of a single peak in the melt curve analyses and low standard errors of the Ct values of replicated measurements indicated the specificity of the primers to bind to the cDNA.

CONCLUSIONS

Our RNA isolation method, with fine-tuned and detailed instructions, can produce high quality RNA from a small amount of starting plant material that is suitable for use in downstream transcriptional analyses. The use of an increased concentration of the reducing agent β-mercaptoethanol in the extraction buffer, as well as the application of DNaseI-treatment resulted in a method suitable for a wide range of plants without the need of further optimalization, especially in Rhus typhina (Staghorn sumac), for which molecular-genetic studies have not yet been sufficiently explored.

Quality control on oil palm RNA samples for efficient genomic downstream applications

Quality control (QC) is primordial for determining the efficiency in any downstream genomic applications. There are several steps in the verification of the quality of RNA samples destined for genomic studies. The aim of this research was to determine whether RNA should be discarded at the level of the field lab if it fails preliminary quality control using Optical Density (OD) measurements. In this study, all samples were submitted to rigorous quality control in every stage of work. RNA samples showing poor OD values still gave excellent results in downstream QC and genomic applications. At the end of the quality control exercise, it was observed that the original samples were the same and had not undergone any deterioration along the different stages of handling and manipulation. This paper shows the different and most important stages of quality control on RNA samples (RIN) for an effective down stream application in genomic studies. RNA samples should not be discarded based on preliminary QC from our field labs.

An Optimized and Cost-Effective RNA Extraction Method for Secondary Metabolite-Enriched Tissues of Norway Spruce (Picea abies)

Since the development of next-generation sequencing techniques and with the growing interest in transcriptomic studies, there is a demand for high-throughput RNA extraction techniques. General RNA extraction protocols are unreliable when it comes to the quality and quantity of isolated RNA obtained from different tissue types of different plant species. Despite Norway spruce (Picea abies) being one of the most significant and commercially valuable tree species in European forests, only limited genetic research is available. In this study, we developed a cetyltrimethylammonium bromide (CTAB) protocol by modifying the original method. We compared this CTAB protocol with other widely used methods for extracting RNA from different tissues (needle, phloem, and root) of Norway spruce, known for its richness in polyphenols, polysaccharides, and secondary metabolites. The modified CTAB method proves to be superior to the kit-based and TRIzol-based methods for extracting RNA from the metabolite-rich tissues of Norway spruce, resulting in high RNA quality and integrity values (RIN~7-9). The modified CTAB RNA extraction method is rapid, cost-effective, and relatively simple in yielding the desired RNA quality from Norway spruce tissues. It is optimal for RNA sequencing and other downstream molecular applications.

Methyl jasmonate reduces cadmium toxicity by enhancing phenol and flavonoid metabolism and activating the antioxidant defense system in pigeon pea (Cajanus cajan)

Cadmium (Cd) is absorbed by plant roots from soil along with essential nutrients and affects plant growth and productivity. Methyl jasmonate (Me-JA) play important roles to mitigate Cd toxicity in plants. We have investigated the role of Me-JA to ameliorate Cd toxicity in Pigeon pea (Cajanus cajan). Plant root growth, biomass, cellular antioxidant defense system and expression of key regulatory genes in molecular and signaling process have been analyzed. Two Cajanus cajan varieties AL-882 and PAU-881 were grown at 25 °C, 16/8h light/dark conditions in three biological replicates at 5 mM Cd concentration, three concentration of Me-JA (0, 10 nM, 100 nM) and two concentrations in combination of Me-JA + Cd (10 nM Me-JA +5 mM Cd, 100 nM Me-JA +5 mM Cd). The seedlings were exposed to Cd stress consequently plants showed decrease in primary root growth (60.71%, in AL-882 and 8.33%, in PAU-881), shoot and root biomass and antioxidant enzymes activities. Me-JA treatment resulted in increased primary root growth (63.64%, in AL-882) and overall plant biomass. Oxidative stress generated due to Cd stress was counter balanced by Me-JA treatment. Me-JA reduced H2O2 free radicals formation and enhanced antioxidant enzyme activities and phenolic content in stressed seedlings. Me-JA treatment increased expression of CALM, IP3, CDPK2, MPKs (involved in calcium and kinase signaling pathways) and reduced expression of metal transporters (IRT1 and HMA3) genes. This reduction in metal transporters gene expression is a probable reason for low toxicity effect of Cd in root after Me-JA treatment which has potential implications in reducing the risk of Cd in the food chain.

A universal protocol for high-quality DNA and RNA isolation from diverse plant species

Next-generation sequencing demands high-quality nucleic acid, yet isolating DNA and RNA is often challenging, particularly from plant tissues. Despite advances in developing various kits and reagents, these products are tailored to isolation of nucleic acid from model plant tissues. Here we introduce a universal lysis buffer to separate nucleic acid from various plant species, including recalcitrant plants, to facilitate molecular analyses, such as quantitative PCR (qPCR), transcriptomics, and whole-genome sequencing (WGS). The protocol is a modification of the original CTAB methods, which leads to nucleic acid isolation from many plant species, including monocots and eudicots. The lysis buffer consists of hexadecyltrimethylammonium bromide (CTAB), sodium chloride (NaCl), Tris base, ethylenediaminetetraacetic acid (EDTA) and β-mercaptoethanol (βME). The modified CTAB method enables the isolation of nucleic acid from small amounts of plant tissues (e.g., 15-100 mg) in a timely manner, which is well-suited for a large number of samples and also when adequate sample collection is a limiting factor. The protocol isolates not only DNA from various plant species but also RNA. This makes it highly effective for molecular analyses compared to previously described CTAB methods optimised for DNA isolation. The appropriate concentration of the components enables high-quality DNA and RNA isolation from plant tissues simultaneously. Additionally, this protocol is compatible with commercially available columns. For DNA and RNA to be qualified for next-generation sequencing platforms, the protocol is supplemented with columns to purify either DNA or RNA from the same tissue to meet high standards for sequencing analyses. This protocol provides an ideal approach to overcome potential obstacles in isolating high-quality DNA or RNA from a wide range of plant species for downstream molecular analysis.

Development of modified CTAB and Trizol protocols to isolate high molecular weight (HMW) RNA from polyphenol and polysaccharides rich pigeonpea (Cajanuscajan (L.) Millsp

Pigeonpea (Cajanuscajan L.) is a legume crop that contains high levels of polyphenolic compounds and polysaccharides that become a hindrance in extracting good-quality and enough amount of RNA from its tissues. With the existing methods of RNA isolation, the phenolic compounds may co-precipitate or bind to the RNA giving false results. Therefore, in the present study, we have modified conventional CTAB and Trizol-based methods which resulted in good quality with the absorbance A260/A280 ratios in the range of 1.83 to 1.98 and A260/230 ratios in the range of 2.0-2.23, revealed RNA to be of high purity and free of contaminants. Both of the proposed protocols yielded a good quantity of RNA ranging from 289 to 422μg per gram of tissue. Distinctly visible bands of 28S and 18S rRNA were observed without degradation or smear, which indicated the presence of intact RNA. RT-PCR analysis showed that isolated RNA was quantitatively sufficient and compliant for the subsequent gene expression analysis.

Comparison of five commercial kits for isolation of total RNA in samples of WSSV-infected shrimp

Viral diseases are the most serious threat to the expansion and development of shrimp aquaculture. Rapid diagnosis of the white spot syndrome virus (WSSV), a lethal shrimp pathogen, is essential to restrict its spread and reduce the mortality of infected shrimp. This virus has globally affected the shrimp farming industry, with a devastating economic impact. Several studies have focused on the expression of WSSV transcripts to understand the molecular mechanisms governing the pathological development of the disease. Since gene expression studies and molecular diagnostics at the early stages of infection depend on the efficient isolation of high-quality RNA, the extraction methods should be carefully selected. However, previous comparisons of the performance of RNA isolation kits have yet to be systematically investigated. In this study, 5 commercial RNA extraction methods were compared in WSSV-infected shrimp. The highest total RNA yield (ng mg-1 tissue) was obtained using TRIzol. Even though the 260/280 nm absorption ratios showed significant differences, the methods showed good purity values (>2.0). RNA integrity was evaluated in a denaturing agarose gel electrophoresis, and degradation was observed after the total RNA samples were treated with DNase I. Finally, the method that allowed the earlier detection of WSSV transcripts by qRT-PCR was the Zymo Direct-zol RNA MiniPrep kit. This study shows that the amount of observed (or estimated) WSSV transcripts might be affected because of the RNA isolation method. In addition, these results may contribute to improve the accuracy of the results obtained in gene expression studies, for more sensitive and robust detection of WSSV.

Genome-wide identification of 2-oxoglutarate and Fe (II)-dependent dioxygenase family genes and their expression profiling under drought and salt stress in potato

The 2-Oxoglutatrate-dependent dioxygenases (2OGDs) comprise the 2-Oxoglutatrate and Fe(II)-dependent dioxygenases (2ODD) enzyme families that facilitate the biosynthesis of various compounds like gibberellin, ethylene, etc. The 2OGDs are also involved in various catabolism pathways, such as auxin and salicylic acid catabolism. Despite their important roles, 2ODDs have not been studied in potato, which is the third most important crop globally. In this study, a comprehensive genome wide analysis was done to identify all 2ODDs in potatoes, and the putative genes were analysed for the presence of the signature 2OG-FeII_Oxy (PF03171) domain and the conserved DIOX_N (PF14226) domain. A total of 205 St2ODDs were identified and classified into eight groups based on their function. The physiochemical properties, gene structures, and motifs were analysed, and gene duplication events were also searched for St2ODDs. The active amino acid residues responsible for binding with 2-oxoglutarate and Fe (II) were conserved throughout the St2ODDs. The three-dimensional (3D) structures of the representative members of flavanol synthase (FNS), 1-aminocyclopropane-1-carboxylic acid oxidases (ACOs), and gibberellin oxidases (GAOXs) were made and docked with their respective substrates, and the potential interactions were visualised. The expression patterns of the St2ODDs under abiotic stressors such as heat, salt, and drought were also analysed. We found altered expression levels of St2ODDs under abiotic stress conditions, which was further confirmed for drought and salt stress using qRT-PCR. The expression levels of St2ODD115, St2ODD34, and St2ODD99 were found to be upregulated in drought stress with 2.2, 1.8, and 2.6 fold changes, respectively. After rewatering, the expression levels were normal. In salt stress, the expression levels of St2ODD151, St2ODD76, St2ODD91, and St2ODD34 were found to be upregulated after 24 hours (h), 48 hours (h), 72 hours (h), and 96 hours (h). Altogether, the elevated expression levels suggest the importance of St2ODDs under abiotic stresses, i.e., drought and salt. Overall, our study provided a knowledge base for the 2ODD gene family in potato, which can be used further to study the important roles of 2ODDs in potato plants.

Genome-wide transcriptional profiling and physiological investigation elucidating the molecular mechanism of multiple abiotic stress response in Stevia rebaudiana Bertoni

Considering the major source of plant-derived low/non-calorie steviol glycosides (SGs), comprehensive physiological, biochemical, and deep transcriptional investigations were conducted to explicit deeper insight into multiple abiotic stress responses in Stevia rebaudiana. The physiological indicators including photosynthesis, chlorophyll, relative water content, shoot growth, electrolyte leakage, and SG biosynthesis were negatively impacted under drought (DS), followed by salinity (SS) and waterlogging (WS). Global transcriptional analysis revealed significant upregulated expression of the genes encoding for ROS detoxification (GST, SOD, APX, glutathione peroxidase), osmotic adjustment (alpha-trehalose-phosphate and S-adenosylmethionine decarboxylase), ion transporters (CAX, NHX, CNGS, VPPase, VATPase), water channel (PIP1, TIP) and abiotic stress-responsive candidate genes (LEA, HSPs, and Dehydrins) regulating abiotic stress response in S. rebaudiana. These inferences were complemented with predicted interactome network that revealed regulation of energy metabolism by key stress-responsive genes (GST, HKT1, MAPKs, P5CSs, PIP), transcription factors (HSFA2, DREB1A, DREB2A), and abiotic stress responsive pathways (ABA, ethylene, ion stress). This is the first detailed study to comprehend the molecular regulation of stress response and their interplay under DS, SS, and WS. The key genes and regulators can be functionally validated, and will facilitate targeted gene editing for genetic improvement of crop sustainability under changing environmental conditions in S. rebaudiana.

Integration of miRNA dynamics and drought tolerant QTLs in rice reveals the role of miR2919 in drought stress response

To combat drought stress in rice, a major threat to global food security, three major quantitative trait loci for 'yield under drought stress' (qDTYs) were successfully exploited in the last decade. However, their molecular basis still remains unknown. To understand the role of secondary regulation by miRNA in drought stress response and their relation, if any, with the three qDTYs, the miRNA dynamics under drought stress was studied at booting stage in two drought tolerant (Sahbaghi Dhan and Vandana) and one drought sensitive (IR 20) cultivars. In total, 53 known and 40 novel differentially expressed (DE) miRNAs were identified. The primary drought responsive miRNAs were Osa-MIR2919, Osa-MIR3979, Osa-MIR159f, Osa-MIR156k, Osa-MIR528, Osa-MIR530, Osa-MIR2091, Osa-MIR531a, Osa-MIR531b as well as three novel ones. Sixty-one target genes that corresponded to 11 known and 4 novel DE miRNAs were found to be co-localized with the three qDTYs, out of the 1746 target genes identified. We could validate miRNA-mRNA expression under drought for nine known and three novel miRNAs in eight different rice genotypes showing varying degree of tolerance. From our study, Osa-MIR2919, Osa-MIR3979, Osa-MIR528, Osa-MIR2091-5p and Chr01_11911S14Astr and their target genes LOC_Os01g72000, LOC_Os01g66890, LOC_Os01g57990, LOC_Os01g56780, LOC_Os01g72834, LOC_Os01g61880 and LOC_Os01g72780 were identified as the most promising candidates for drought tolerance at booting stage. Of these, Osa-MIR2919 with 19 target genes in the qDTYs is being reported for the first time. It acts as a negative regulator of drought stress tolerance by modulating the cytokinin and brassinosteroid signalling pathway.

Physio-Biochemical Integrators and Transcriptome Analysis Reveal Nano-Elicitation Associated Response during Dendrocalamus asper (Schult. and Schult. F.) Backer ex K. Heyne Micropropagation

Bamboos are perennial, arborescent, monocarpic and industrially important non-timber plants. They are important for various purposes, such as carbon sequestration, biodiversity support, construction, and food and fiber production. However, traditional vegetative propagation is insufficient for bamboo multiplication. Moreover, little is known about the mechanism of gold nanoparticles (AuNPs) in vitro proliferation and regulation of physiological and biochemical properties. In this study, we investigated the impacts of citrate and cetyltrimethylammonium bromide (CTAB) coated AuNPs on in vitro proliferation, photosynthetic pigment content and antioxidant potential of Dendrocalamus asper (Schult. and Schult. F.) Backer ex K. Heyne. Various morpho-physiological and biochemical parameters were differentially affected along the citrate- and CTAB-coated AuNPs concentration gradients (200-600 µM). In vitro shoot proliferation, photosynthetic pigment content and antioxidant activities were higher in D. asper grown on Murashige and Skoog medium supplemented with 2 mg·L-1 benzyladenine and 400 µM citrate-coated AuNPs than in those grown on Murashige and Skoog medium supplemented with 600 µM CTAB- coated AuNPs. Identification of genes regulating in vitro D. asper proliferation will help understand the molecular regulation of AuNPs-mediated elicitation for modulating various physiological and biochemical activities during micropropagation. Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway analyses identified differentially expressed genes associated with in vitro modulation of AuNPs-regulated biological processes and molecular functions. The findings of this study provide new insight into AuNPs-mediated elicitation of in vitro mass scale bamboo propagation.

Systematic analyses with genomic and metabolomic insights reveal a new species, Ophiocordyceps indica sp. nov. from treeline area of Indian Western Himalayan region

Ophiocordyceps is a species-rich genus in the order Hypocreales (Sordariomycetes, Ascomycota) depicting a fascinating relationship between microbes and insects. In the present study, a new species, Ophiocordyceps indica sp. nov., is discovered infecting lepidopteran larvae from tree line locations (2,202-2,653 m AMSL) of the Kullu District, Himachal Pradesh, Indian Western Himalayan region, using combinations of morphological and molecular phylogenetic analyses. A phylogeny for Ophiocordyceps based on a combined multigene (nrSSU, nrLSU, tef-1α, and RPB1) dataset is provided, and its taxonomic status within Ophiocordycipitaceae is briefly discussed. Its genome size (~59 Mb) revealed 94% genetic similarity with O. sinensis; however, it differs from other extant Ophiocordyceps species based on morphological characteristics, molecular phylogenetic relationships, and genetic distance. O. indica is identified as the second homothallic species in the family Ophiocordycipitaceae, after O. sinensis. The presence of targeted marker components, viz. nucleosides (2,303.25 μg/g), amino acids (6.15%), mannitol (10.13%), and biological activity data, suggests it to be a new potential source of nutraceutical importance. Data generated around this economically important species will expand our understanding regarding the diversity of Ophiocordyceps-like taxa from new locations, thus providing new research avenues.

Integrated analysis of smRNAome, transcriptome, and degradome data to decipher microRNAs regulating costunolide biosynthesis in Saussurea lappa

Saussurea lappa (S. lappa) has been known to synthesize medicinally important, costunolide. Due to its immense therapeutic importance, understanding of regulatory mechanism associated with its biosynthesis is crucial. The identification of genes and transcription factors (TFs) in S. lappa, created a clear picture of costunolide biosynthesis pathways. Further to understand the regulation of costunolide biosynthesis by miRNAs, an integrated study of transcriptome, miRNAs, and degradome was performed. Identified candidate miRNAs and associated feed-forward loops (FFLs) illustrates their regulatory role in secondary metabolite biosynthesis. Small RNA and degradome sequencing were performed for leaf and root tissues to determine miRNAs-targets pairs. A total of 711 and 525 such targets were obtained for novel and known miRNAs respectively. This data was used to generate costunolide-specific miRNA-TF-gene interactome to perform systematic analyses through graph theoretical approach. Interestingly, miR171c.1 and sla-miR121 were identified as key regulators to connect and co-regulate both mevalonate and sesquiterpenoid pathways to bio-synthesize costunolide. Tissue-specific FFLs were identified to be involved in costunolide biosynthesis which further suggests the evolutionary co-relation of root-specific networks in synthesis of secondary metabolites in addition to leaf-specific networks. This integrative approach allowed us to determine candidate miRNAs and associated tissue-specific motifs involved in the diversification of secondary metabolites. MiRNAs identified in present study can provide alternatives for bioengineering tool to enhance the synthesis of costunolide and other secondary metabolites in S. lappa.

Heterologous expression of Chlorophytum borivilianum Squalene epoxidase in tobacco modulates stigmasterol production and alters vegetative and reproductive growth

KEYMESSAGE

CbSE overexpression increased stigmasterol levels and altered plant morphology. The genes upstream and downstream of CbSE were found to be upregulated, which confirms its regulatory role in the saponin biosynthetic pathway. Chlorophytum borivilianum is a high-value medicinal plant with many promising preclinical applications that include saponins as a major active ingredient. Squalene epoxidase (SE) is one of the major rate-limiting enzymes of the saponin biosynthetic pathway. Here, we functionally characterized C. borivilianum SE (CbSE) by over-expressing heterologously in Nicotiana tabacum. The heterologous expression of CbSE resulted in stunted pant growth with altered leaf and flower morphology. Next, RT-qPCR analysis of transgenic plants overexpressing CbSE revealed increased expression levels of Cycloartenol synthase (CAS), Beta amyrin synthase (βAS), and cytochrome P450 monooxygenase 51 (CYP51) (Cytochrome P450), which encode key enzymes for triterpenoid and phytosterol biosynthesis in C. borivilianum. Further, Methyl Jasmonate (MeJa) treatment upregulated Squalene synthase (SQS), SE, and Oxidosqualene cyclases (OSCs) to a significant level. GC-MS analysis of the leaf and hairy roots of the transformants showed an increased stigmasterol content (0.5-1.0 fold) compared to wild type (WT) plants. These results indicate that CbSE is a rate-limiting gene, which encodes an efficient enzyme responsible for phytosterol and triterpenoid production in C. borivilianum.

Genome-wide analysis of Catalase gene family reveal insights into abiotic stress response mechanism in Brassica juncea and B. rapa

Environmental stresses affect the yield and productivity of Brassica crops. Catalases are important antioxidant enzymes involved in reducing excess hydrogen peroxide produced by environmental stresses. In the present study, nine and seven CAT family members in two oilseed Brassica species (B. juncea and B. rapa) were identified with complete characterization based on gene and protein structure. Phylogenetic classification categorized CAT proteins into three classes and differentiated the monocot and dicot-specific CAT proteins. Further, the gene and protein characterizations revealed a high degree of conservation across the CAT family members. Differences were observed in the CAT-HEME binding affinity in CAT1, CAT2, and CAT3 isozymes, which could suggest their differential enzyme activities in different conditions. Furthermore, protein-protein interaction with other antioxidant proteins suggested their coordinated role in ROS scavenging mechanisms. Notably, the differential gene expression of BjuCATs and BraCATs and CAT enzyme activities suggested their crucial roles in major abiotic stresses faced by Brassica species. Promoter analysis in BjuCATs and BraCATs suggested the presence of abiotic-stress responsive cis-regulatory elements. Gene regulatory network analysis suggested miRNA and TF mediated stress response in BjuCATs and BraCATs. CAT family screening and characterization in Brassica sp. has established a basic ground for further functional validation in abiotic and heavy-metal stresses which can help in developing stress tolerant crops.

An optimized RNA extraction method for diverse leaves of Hawaiian Metrosideros, a hypervariable tree species complex

Premise

The isolation of RNA from trees is challenging due to the interference of polyphenols and polysaccharides with downstream processes. Furthermore, many RNA extraction protocols are time consuming and involve hazardous chemicals. To address these issues, we aimed to develop a safe protocol for high-quality RNA extraction from diverse Metrosideros taxa representing a broad range of leaf toughness, pubescence, and secondary metabolites.

Methods and Results

We tested popular RNA isolation kits and protocols that were effective on other recalcitrant trees, including a broad range of optimization and purification steps. We optimized a protocol involving two silica-membrane column-based kits that yielded high-quantity RNA with an RNA integrity number >7 and without DNA contamination. All RNA samples were used successfully in a follow-on RNA-Seq experiment.

Conclusions

We present an optimized high-throughput RNA extraction protocol that yielded high-quality and high-quantity RNA from three contrasting leaf phenotypes within a hyperdiverse woody species complex.

Transcriptome analysis of gall oak (Quercus infectoria): De novo assembly, functional annotation and metabolic pathways analysis

Gall oak (Quercus infectoria) is a native tree of Iran, whose gall extract is used to treat many diseases. The presence of abundant secondary metabolites with various bioactivities in this plant has made it medically important. Despite its medicinal value, due to the lack of genomic information, the biosynthetic pathways of these compounds in this species are still unknown. The current research was aimed at observing, characterizing, and investigating the biosynthetic pathways of these compounds in Q.infectoria. De novo transcriptome assembly was conducted using the RNA sequencing technique. A total of 89,335 unigenes were generated, of which 6928 unigenes showed differential expression in leaves compared to root tissue. Gene ontology examination of DEGs revealed GO-term enrichment was related to cellular processes and enzyme activity. KEGG enrichment analysis for DEGs showed that most unigenes were related to metabolic pathways and biosynthesis of secondary metabolites. Moreover, 39 families of transcription factors were identified, of which the C2H2, bZIP, bHLH, and ERF TFs had the highest frequency. In the absence of a reference genome, the overall study of transcriptome will provide a reference for future functional and comparative studies. Moreover, the data obtained from sequencing and de novo assembly can be a valuable scientific resource for Q.infectoria.

Lemon grass essential oil improves Gladiolus grandiflorus postharvest life by modulating water relations, microbial growth, biochemical activity, and gene expression

Gladiolus (Gladiolus grandiflorus Andrews) is a high-valued bulbous cut flower. However, the shorter postharvest life of the gladiolus, limits its marketing and commercial value. In the present investigation, the effect of lemon grass (LG) essential oil as an antimicrobial agent was studied towards increasing the vase life of gladiolus. The results revealed that as compared to control (distilled water), treatment with a lower concentration of 5 µL L^-1 LG essential oil prolonged the vase life of gladiolus up to 11 days (d). Scanning Electron Microscope (SEM) observation indicated that the sample treated with 5 µL L^-1 LG essential oil showed intact vasculature, suggesting reduced microbial blockage at the stem end which was further corroborated by microbial count. Biochemical analysis suggested an increased level of total soluble sugars, carotenoid content, lower MDA accumulation, and higher activity of antioxidant enzymes in LG treated flowers. Moreover, transcripts levels of genes associated with senescence viz., GgCyP1 and GgERS1a were downregulated, while expression of GDAD1 and antioxidant genes such as GgP5C5, GgPOD 1, GgMnSOD, and GgCAT1 were upregulated in LG treated cut spikes as compared to control. Among various treatments we have concluded that, the vase life of the gladiolus cut spike was improved along with the relative fresh flower weight and diameter of flower at the lower dose of 5 µL L^-1 LG oil in the vase solution. Thus, LG oil as an eco-friendly agent has the potential to extend the postharvest life of cut flowers.

Molecular Detection and Phylogenetic Analysis of Deformed Wing Virus and Sacbrood Virus Isolated from Pollen

Among many pathogens and pests, honey bee viruses are known as one of the most common cause of diseases in honey bee colonies. In this study, we demonstrate that pollen grains and bee bread are potential sources of viral DNA. We extracted DNA from 3 types of pollen samples: directly provided by beekeepers (n = 12), purchased from trade markets (n = 5), and obtained from honeycombs (bee bread, n = 10). The extracted DNA was used for molecular detection (RT-PCR analysis) of six of the most widely distributed honey bee viruses: deformed wing virus, sacbrood virus, acute bee paralysis virus, black queen cell virus, Kashmir bee virus, Israeli acute paralysis virus, and chronic bee paralysis virus. We successfully managed to establish only the deformed wing virus (DWV) and the sacbrood virus (SBV), with different distribution frequencies depending on the territory of the country. The phylogenetic analyses of Bulgarian isolates were performed with the most similar sequences available in molecular databases from other countries. Phylogenies of Bulgarian viral strains demonstrated genetically heterogeneous populations of DWV and relatively homogenous populations of SBV. In conclusion, the results obtained from the current study have shown that pollen is a valuable source for molecular detection of honey bee pathogens. This allows epidemiological monitoring of honey bee diseases at a regional and a national level.

Strategies to study the metabolic origins of specialized plant metabolites: The specialized 1,4-naphthoquinones

One of the hallmarks of specialized plant metabolites is that they are produced using precursors from central metabolism. Therefore, in addition to identifying and characterizing the pathway genes and enzymes involved in synthesizing a specialized compound, it is critical to study its metabolic origins. Identifying what primary metabolic pathways supply precursors to specialized metabolism and how primary metabolism has diversified to sustain fluxes to specialized metabolite pathways is imperative to optimizing synthetic biology strategies for producing high-value plant natural products in crops and microbial systems. Improved understanding of the metabolic origins of specialized plant metabolites has also revealed instances of recurrent evolution of the same compound, or nearly identical compounds, with similar ecological functions, thereby expanding knowledge about the factors driving the chemical diversity in the plant kingdom. In this chapter, we describe detailed methods for performing tracer studies, chemical inhibitor experiments, and reverse genetics. We use examples from investigations of the metabolic origins of specialized plant 1,4-naphthoquinones (1,4-NQs). The plant 1,4-NQs provide an excellent case study for illustrating the importance of investigating the metabolic origins of specialized metabolites. Over half a century of research by many groups has revealed that the pathways to synthesize plant 1,4-NQs are the result of multiple events of convergent evolution across several disparate plant lineages and that plant 1,4-NQ pathways are supported by extraordinary events of metabolic innovation and by various primary metabolic sources.

Polyamine metabolizing rhizobacteria Pseudomonas sp. GBPI_506 modulates hormone signaling to enhance lateral roots and nicotine biosynthesis in Nicotiana benthamiana

Beneficial rhizobacteria in the soil are important drivers of plant health and growth. In this study, we provide the draft genome of a root colonizing and auxin-producing Pseudomonas sp. strain GBPI_506. The bacterium was investigated for its contribution in the growth of Nicotiana benthamiana (Nb) and biosynthesis of nicotine. The bacterium showed chemotaxis towards root exudates potentially mediated by putrescine, a polyamine compound, to colonize the roots of Nb. Application of the bacterium with the roots of Nb, increased plant biomass and total soluble sugars in the leaves, and promoted lateral root (LR) development as compared to the un-inoculated plants. Confocal analysis using transgenic (DR5:GFP) Arabidopsis showed increased auxin trafficking in the LR of inoculated plants. Upregulation of nicotine biosynthesis genes and genes involved in salicylic acid (SA) and jasmonic acid (JA) signaling in the roots of inoculated plants suggested increased nicotine biosynthesis as a result of bacterial application. An increased JA content in roots and nicotine accumulation in leaves provided evidence on JA-mediated upregulation of nicotine biosynthesis in the bacterized plants. The findings suggested that the bacterial root colonization triggered networking between auxin, SA, and JA to facilitate LR development leading to enhanced plant growth and nicotine biosynthesis in Nb.

Seed priming with brassinolides improves growth and reinforces antioxidative defenses under normal and heat stress conditions in seedlings of Brassica juncea

Environmental stresses pose a major challenge for plant researchers to fulfill increasing food demand. Researchers are trying to generate high-yielding and stress-tolerant or resistant varieties using classical genetics and modern gene-editing tools; however, both approaches have limitations. Chemical treatments emerged as an alternative to improve yield and impart stress resilience. Brassinosteroids (BRs) are a group of phytohormones that regulate various biological processes, including stress management. With foliar spray methods, BR treatments showed promising results but are not economically feasible. We hypothesize that priming of seeds, which requires lesser amounts of BRs, could be equally effective in promoting growth and stress tolerance. Owing to this notion, we analyzed the impact of priming seeds with selected BRs, namely, 24-epibrassinolide (EBL) and 28-homobrassinolide (HBL), in Brassica juncea under normal and heat shock stress conditions. Seeds primed with BRs and grown until seedlings stage at normal conditions (20°C) were subjected to a heat shock (35°C) for a few hours, relating to what plants experience in natural conditions. Heat shock reduced the growth and biomass with an increased accumulation of reactive oxygen species. As anticipated, BRs treatments significantly improved the growth and physiological parameters with an enhanced antioxidant defense under both conditions. Transcriptional analyses revealed that BRs concomitantly induce growth and oxidative stress-responsive gene expression via the canonical BR-signaling pathway. Transfer of unstressed and heat-shock-treated seedlings to field conditions demonstrated the long-term effectivity of BR-priming. Our results showed seed priming with BRs could improve growth and resilience against heat shock; hence, it appears to be a viable strategy to enhance crop yields and stress tolerance.

Charged Gold Nanoparticles Promote In Vitro Proliferation in Nardostachys jatamansi by Differentially Regulating Chlorophyll Content, Hormone Concentration, and Antioxidant Activity

Nardostachys jatamansi is a critically endangered medicinal plant and endemic to the Himalayas, having high commercial demand globally. The accumulation of various secondary metabolites in its shoots and roots with antioxidant potential are well-documented in traditional as well as modern medicine systems. In the present study, we first attempted to investigate the impact of citrate (−ve charge, 11.1 ± 1.9 nm) and CTAB (+ve charge, 19.5 ± 3.2 nm) coated gold nanoparticles (AuNPs) on the in vitro proliferation and antioxidant activities of N. jatamansi. Both the nanoparticles differentially affected the morphological and biochemical parameters, chlorophyll content, internal hormone concentration, and antioxidant activities in a concentration-dependent (10–100 µM) manner. Vigorous shooting was observed in half strength MS medium supplemented with IAA (1 mg/L) with 60 µM citrate-AuNPs (46.4 ± 3.7 mm) and 40 µM CTAB-AuNPs (42.2 ± 3.2 mm). Similarly, the maximum number of roots (5.00 ± 0.67 and 5.33 ± 0.58) and root length (29.9 ± 1.5 mm and 27.3 ± 4.8 mm) was reported in half-strength MS medium with IAA (1 mg/L) supplemented with 60 µM citrate-AuNPs and 40 µM CTAB-AuNPs, respectively. In addition, plants growing on MS medium supplemented with 60 µM citrate-AuNPs and 40 µM CTAB-AuNPs showed significantly enhanced photosynthetic pigments (chlorophyll a and b, carotenoids, and total chlorophyll), internal hormone concentration (GA3, IAA, and ABA), and antioxidant activities (total phenolics, flavonoids, DPPH, and SOD enzyme activity). Moreover, the transcript analysis of ANR1, ARF18, PLY9, SAUR28, GID1A, GRF1, SOD, and CAT further confirmed the role of 60 µM citrate-AuNPs and 40 µM CTAB-AuNPs in the improvement in the growth and antioxidant activities of N. jatamansi. Bearing in mind the urgent requirements of the effective conservation measures of this endangered species, the present findings suggest the elicitation of citrate-AuNPs and CTAB-AuNPs would significantly improve the potential applications of N. jatamansi in the medicinal plant-based industry.

An Integrative Transcriptional Network Revealed Spatial Molecular Interplay Underlying Alantolactone and Inulin Biosynthesis in Inula racemosa Hook f

Inula racemosa Hook. f. (Pushkarmula), a perennial Himalayan herb known for its aromatic and phytopharmaceutical attributes, is not yet explored at genomic/transcriptomic scale. In this study, efforts were made to unveil the global transcriptional atlas underlying organ-specific specialized metabolite biosynthesis by integrating RNA-Seq analysis of 433 million sequenced reads with the phytochemical analysis of leaf, stem, and root tissues. Overall, 7242 of 83,772 assembled nonredundant unigenes were identified exhibiting spatial expression in leaf (3761), root (2748), and stem (733). Subsequently, integration of the predicted transcriptional interactome network of 2541 unigenes (71,841 edges) with gene ontology and KEGG pathway enrichment analysis revealed isoprenoid, terpenoid, diterpenoid, and gibberellin biosynthesis with antimicrobial activities in root tissue. Interestingly, the root-specific expression of germacrene-mediated alantolactone biosynthesis (GAS, GAO, G8H, IPP, DMAP, and KAO) and antimicrobial activities (BZR1, DEFL, LTP) well-supported with both quantitative expression profiling and phytochemical accumulation of alantolactones (726.08 μg/10 mg) and isoalantolactones (988.59 μg/10 mg), which suggests “roots” as the site of alantolactone biosynthesis. A significant interaction of leaf-specific carbohydrate metabolism with root-specific inulin biosynthesis indicates source (leaf) to sink (root) regulation of inulin. Our findings comprehensively demonstrate the source-sink transcriptional regulation of alantolactone and inulin biosynthesis, which can be further extended for upscaling the targeted specialized metabolites. Nevertheless, the genomic resource created in this study can also be utilized for development of genome-wide functionally relevant molecular markers to expedite the breeding strategies for genetic improvement of I. racemosa.

Spatial Genomic Resource Reveals Molecular Insights into Key Bioactive-Metabolite Biosynthesis in Endangered Angelica glauca Edgew

Angelica glauca Edgew, which is an endangered medicinal and aromatic herb, is a rich source of numerous industrially important bioactive metabolites, including terpenoids, phenolics, and phthalides. Nevertheless, genomic interventions for the sustainable utilization and restoration of its genetic resources are greatly offset due to the scarcity of the genomic resources and key regulators of the underlying specialized metabolism. To unravel the global atlas of the specialized metabolism, the first spatial transcriptome sequencing of the leaf, stem, and root generated 109 million high-quality paired-end reads, assembled de novo into 81,162 unigenes, which exhibit a 61.53% significant homology with the six public protein databases. The organ-specific clustering grouped 1136 differentially expressed unigenes into four subclusters differentially enriched in the leaf, stem, and root tissues. The prediction of the transcriptional-interactome network by integrating enriched gene ontology (GO) and the KEGG metabolic pathways identified the key regulatory unigenes that correspond to terpenoid, flavonoid, and carotenoid biosynthesis in the leaf tissue, followed by the stem and root tissues. Furthermore, the stem and root-specific significant enrichments of phenylalanine ammonia lyase (PAL), cinnamate-4-hydroxylase (C4H), and caffeic acid 3-O-methyltransferase (COMT) indicate that phenylalanine mediated the ferulic acid biosynthesis in the stem and root. However, the root-specific expressions of NADPH-dependent alkenal/one oxidoreductase (NADPH-AOR), S-adenosyl-L-methionine-dependent methyltransferases (SDMs), polyketide cyclase (PKC), and CYP72A15 suggest the “root” as the primary site of phthalide biosynthesis. Additionally, the GC-MS and UPLC analyses corresponded to the organ-specific gene expressions, with higher contents of limonene and phthalide compounds in the roots, while there was a higher accumulation of ferulic acid in the stem, followed by in the root and leaf tissues. The first comprehensive genomic resource with an array of candidate genes of the key metabolic pathways can be potentially utilized for the targeted upscaling of aromatic and pharmaceutically important bioactive metabolites. This will also expedite genomic-assisted conservation and breeding strategies for the revival of the endangered A. glauca.

A Novel Bacillus safensis-Based Formulation along with Mycorrhiza Inoculation for Controlling Alternaria alternata and Simultaneously Improving Growth, Nutrient Uptake, and Steviol Glycosides in Stevia rebaudiana under Field Conditions

The excess use of chemicals by farmers in the agroecosystems degrades soil quality, disturbs soil ecology, and increases soil salinity and health hazards in humans. Stevia rebaudiana is an important medicinal and aromatic crop whose leaves contain steviol glycosides (SGs). The Bacillus safensis NAIMCC-B-02323 strain STJP from the rhizosphere of S. rebaudiana producing salicylic acid (16.80 µg/mL), chitinase (75.58 U/mL), β-1,3-glucanase (220.36 U/mL), and cellulase (170 U/mL) was taken as a plant growth-promoting rhizobacteria (PGPR). The cell-free supernatant (CFS) from strain STJP showed significant biocontrol activity against Alternaria alternata (80%), suggesting the protective role of extracellular metabolite(s) against phytopathogens. Paneer whey-based bioformulation (P-WBF) was developed to exploit B. safensis STJP to enhance the growth, nutrient uptake, soil properties, stevioside content, and SGs biosynthesis in S. rebaudiana under an A. alternata-infested field. The combined treatment of P-WBF and mycorrhiza (Glomus fasciculatum ABTEC) significantly enhanced plant growth parameters after 90 days, in comparison with control. The symbiotic action (P-WBF and mycorrhiza) displayed much better results in terms of chlorophyll a and b (improved by 132.85% and 39.80%, respectively), protein (by 278.75%), flavonoid (by 86.99%), carbohydrate (by 103.84%), antioxidant (by 75.11%), and stevioside (by 120.62%) contents in plants as compared to the untreated set. Further, the augmentation of potassium (by 132.39%), phosphorous (by 94.22%), and zinc (by 111.11%) uptake in plant tissues and soil was also observed by the application of P-WBF and mycorrhiza. The expression of UGT74G1 and UGT85C2 genes related to SG biosynthesis was upregulated (2.7- and 3.2-fold, respectively) in plants treated with P-WBF and mycorrhiza as further confirmed by the accumulation of SGs. The results suggest that the application of P-WBF and mycorrhiza not only provides an ecofriendly and sustainable solution to improve stevioside content in S. rebaudiana by a nutrient-linked mechanism but also paves the way to enhanced production of stevioside.

Isolation of high-quality RNA for high throughput applications from secondary metabolite-rich Crocus sativus L

Objective

Isolating high-quality RNA is a basic requirement while performing high throughput sequencing, microarray, and various other molecular investigations. However, it has been quite challenging to isolate RNA with absolute purity from plants like Crocus sativus that are rich in secondary metabolites, polysaccharides, and other interfering compounds which often irreversibly co-precipitate with the RNA. While many methods have been proposed for RNA extraction including CTAB, TriZol, and SDS-based methods, which invariably yield less and poor quality RNA and hence it necessitated the isolation of high-quality RNA suitable for high throughput applications.

Results

In the present study we made certain adjustments to the available protocols including modifications in the extraction buffer itself and the procedure employed. Our method led to the isolation of clear and non-dispersive total RNA with an RNA Integrity Number (RIN) value greater than 7.5. The quality of the RNA was further assessed by qPCR-based amplification of mRNA and mature miRNAs such as Cs-MIR166c and Cs-MIR396a.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-06095-z.

Characterization of APX and APX-R gene family in Brassica juncea and B. rapa for tolerance against abiotic stresses

APX and APX-R gene families were identified and characterized in two important oilseed species of Brassica. Gene expression under abiotic stress conditions, recombinant protein expression, and analysis further divulged their drought, heat, and salt-responsive behavior. Ascorbate peroxidases (APX) are heme-dependent enzymes that rid the cells of H₂O₂ and regulate diverse biological processes. In the present study, we performed APX gene family characterization in two Brassica sp. (B. juncea and B. rapa) as these are commercially important oilseed crops and affected severely by abiotic stresses. We identified 16 BjuAPX and 9 BraAPX genes and 2 APX-R genes each in B. juncea and B. rapa genomes, respectively. Phylogenetic analysis divided the APX genes into five distinct clades, which exhibited conservation in the gene structure, motif organization, and sub-cellular location within the clade. Structural analysis of APX and APX-R proteins revealed the amino acid substitutions in conserved domains of APX-R proteins. The expression profiling of BjuAPX and BraAPX genes showed that 3 BjuAPX, 7BraAPX, and 2 BraAPX-R genes were drought and heat responsive. Notably, BjuAAPX1a, BjuAPX1d, BjuAAPX6, BraAAPX1a, BraAAPX2, and BraAAPX3b showed high expression levels in RT-qPCR. Cis-regulatory elements in APX and APX-R gene promoters supported the differential behavior of these genes. Further, two stress-responsive genes BjuAPX1d and BraAAPX2 were cloned, characterized, and their roles were validated under heat, drought, salt, and cold stress in bacterial expression system. This study for the first time reports the presence of APX activity in dimeric and LMW form of purified BraAAPX2 protein. The study may help pave way for developing abiotic stress-tolerant Brassica crops.

Genome-Wide Analysis of TIR-NBS-LRR Gene Family in Potato Identified StTNLC7G2 Inducing Reactive Oxygen Species in Presence of Alternaria solani

Resistance gene analogs (RGAs) comprising NBS-LRR gene family members are considered prominent candidates in the development of disease-resistant genotypes. NBS-LRR gene family comprised a very large number of genes; therefore, members of one subfamily TIR-NBS-LRR (TNL) are identified in the present study from Solanum tuberosum genome, followed by their bioinformatics characterization. The study identified a total of 44 genes encoding 60 TNL transcripts with two prominent clusters at chromosome 1 and chromosome 11. Expression analysis of 14 TNL genes after Alternaria solani infection at 1, 2, 3, 5, and 7 days post inoculation in two disease-tolerant varieties, Kufri Jyoti and Kufri Pukhraj, and one relatively susceptible variety, Kufri Chandramukhi, showed differential expression of many genes including a high expression (>15-fold) of StTNLC6G2T1 and StTNLC11G9T1. Functional characterization of one such gene, StTNLC7G2, reveals involvement in the generation of reactive oxygen species under A. solani attack, implicating its putative role in plant defense via hypersensitive response.

Genome-wide identification and characterization of functionally relevant microsatellite markers from transcription factor genes of Tea (Camellia sinensis (L.) O. Kuntze)

Tea, being one of the most popular beverages requires large set of molecular markers for genetic improvement of quality, yield and stress tolerance. Identification of functionally relevant microsatellite or simple sequence repeat (SSR) marker resources from regulatory “Transcription factor (TF) genes” can be potential targets to expedite molecular breeding efforts. In current study, 2776 transcripts encoding TFs harbouring 3687 SSR loci yielding 1843 flanking markers were identified from traits specific transcriptome resource of 20 popular tea cultivars. Of these, 689 functionally relevant SSR markers were successfully validated and assigned to 15 chromosomes (Chr) of CSS genome. Interestingly, 589 polymorphic markers including 403 core-set of TF-SSR markers amplified 2864 alleles in key TF families (bHLH, WRKY, MYB-related, C2H2, ERF, C3H, NAC, FAR1, MYB and G2-like). Their significant network interactions with key genes corresponding to aroma, quality and stress tolerance suggests their potential implications in traits dissection. Furthermore, single amino acid repeat reiteration in CDS revealed presence of favoured and hydrophobic amino acids. Successful deployment of markers for genetic diversity characterization of 135 popular tea cultivars and segregation in bi-parental population suggests their wider utility in high-throughput genotyping studies in tea.

Stress-Responsive cis-Regulatory Elements Underline Podophyllotoxin Biosynthesis and Better Performance of Sinopodophyllum hexandrum Under Water Deficit Conditions

Sinopodophyllum hexandrum is an endangered medicinal herb known for its bioactive lignan podophyllotoxin (PTOX), which is used for the preparation of anticancer drugs. In its natural habitat, S. hexandrum is exposed to a multitude of adversities, such as fluctuating temperatures, water deficit, and high UV radiations. Transcriptional regulation of genes, which is regulated by the condition-specific binding of transcriptional factors to precise motifs in the promoter region, underlines responses to an environmental cue. Therefore, analysis of promoter sequences could ascertain the spatio-temporal expression of genes and overall stress responses. Unavailability of genomic information does not permit such analysis in S. hexandrum, especially on regulation of PTOX pathway. Accordingly, this study describes isolation and in silico analysis of 5'-upstream regions of ShPLR (PINORESINOL-LARICIRESINOL REDUCTASE) and ShSLD (SECOISOLARICIRESINOL DEHYDROGENASE), the two key genes of the PTOX biosynthetic pathway. Data showed a range of motifs related to basal transcription, stress-responsive elements, such as those for drought, low temperature, and light, suggesting that the expression of these genes and resulting PTOX accumulation would be affected by, at least, these environmental cues. While the impact of temperature and light on PTOX accumulation is well studied, the effect of water deficit on the physiology of S. hexandrum and PTOX accumulation remains obscure. Given the presence of drought-responsive elements in the promoters of the key genes, the impact of water deficit on growth and development and PTOX accumulation was studied. The results showed decline in relative water content and net photosynthetic rate, and increase in relative electrolyte leakage with stress progression. Plants under stress exhibited a reduction in transpiration rate and chlorophyll content, with a gradual increase in osmoprotectant content. Besides, stressed plants showed an increase in the expression of genes involved in the phenylpropanoid pathway and PTOX biosynthesis, and an increase in PTOX accumulation. Upon re-watering, non-irrigated plants showed a significant improvement in biochemical and physiological parameters. Summarily, our results demonstrated the importance of osmoprotectants during water deficit and the revival capacity of the species from water deficit, wherein PTOX synthesis was also modulated. Moreover, isolated promoter sequences could be employed in genetic transformation to mediate the expression of stress-induced genes in other plant systems.

AtUSP17 negatively regulates salt stress tolerance through modulation of multiple signaling pathways in Arabidopsis

AtUSP17 is a multiple stress-inducible gene that encodes a universal stress protein (USP) in Arabidopsis thaliana. In the present study, we functionally characterized AtUSP17 using its knock-down mutant, Atusp17, and AtUSP17-overexpression lines (WTOE). The overexpression of AtUSP17 in wild-type and Atusp17 mutant Arabidopsis plants resulted in higher sensitivity to salt stress during seed germination than WT and Atusp17 mutant lines. In addition, the WTOE and FC lines exhibited higher abscisic acid (ABA) sensitivity than Atusp17 mutant during germination. The exogenous application of ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) was able to rescue the salt hypersensitive phenotype of WTOE lines. In contrast, AgNO₃ , an ethylene action inhibitor, further blocked the effect of ACC during germination. The addition of ACC under salt stress resulted in reduced reactive oxygen species (ROS) accumulation, expression of ABA-responsive genes, improved proline synthesis, increased expression of positive regulators of ethylene signaling and antioxidant defense genes with enhanced antioxidant enzyme activities. The WTOE lines exhibited salt sensitivity even at the adult plant stage, while Atusp17 mutant exhibited higher salt tolerance with higher chlorophyll, relative water content and lower electrolyte leakage as compared with WT. The BAR interaction viewer database and available literature mining identified AtUSP17-interacting proteins, which include RGS1, RACK1C and PRN1 involved in G-protein signaling, which play a crucial role in salt stress responses. Based on the present study and available literature, we proposed a model in which AtUSP17 negatively mediates salt tolerance in Arabidopsis through modulation of ethylene, ABA, ROS, and G-protein signaling and responses.

Identification and validation of reference genes for qRT-PCR based studies in horse gram (Macrotyloma uniflorum)

The quantitative real-time polymerase chain reaction (qRT-PCR) is the most sensitive and commonly used technique for gene expression studies in biological systems. However, the reliability of qRT-PCR results depends on the selection of reference gene(s) for data normalization. Horse gram (Macrotyloma uniflorum) is an important legume crop on which several molecular studies have been reported. However, the stability of reference genes has not been evaluated. In the present study, nine candidate reference genes were identified from horse gram RNA-seq data and evaluated in two horse gram genotypes, HPK4 and HPKM317 under six abiotic stresses viz. cold, drought, salinity, heat, abscisic acid and methyl viologen-induced oxidative stress. The results were evaluated using geNorm, Bestkeeper, Normfinder and delta-delta Ct methods and comprehensive ranking was assigned using RefFinder and RankAggreg software. The overall result showed that TCTP was one of the most stable genes in all samples and in genotype HPK4, while in HPKM317 profilin was most stably expressed. However, PSMA5 was identified as least stable in all the experimental conditions. Expression of target genes dehydrin and early response to dehydration 6 under drought stress was also validated using TCTP and profilin for data normalization, either alone or in combination, which confirmed their suitability for qRT-PCR data normalization. Thus, TCTP and profilin genes may be used for qRT-PCR data normalization for molecular and genomic studies in horse gram.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01104-0.

An optimized method for RNA extraction from the polyurethane oligomer degrading strain Pseudomonas capeferrum TDA1 growing on aromatic substrates such as phenol and 2,4-diaminotoluene

Bacterial degradation of xenobiotic compounds is an intense field of research already for decades. Lately, this research is complemented by downstream applications including Next Generation Sequencing (NGS), RT-PCR, qPCR, and RNA-seq. For most of these molecular applications, high-quality RNA is a fundamental necessity. However, during the degradation of aromatic substrates, phenolic or polyphenolic compounds such as polycatechols are formed and interact irreversibly with nucleic acids, making RNA extraction from these sources a major challenge. Therefore, we established a method for total RNA extraction from the aromatic degrading Pseudomonas capeferrum TDA1 based on RNAzol® RT, glycogen and a final cleaning step. It yields a high-quality RNA from cells grown on TDA1 and on phenol compared to standard assays conducted in the study. To our knowledge, this is the first report tackling the problem of polyphenolic compound interference with total RNA isolation in bacteria. It might be considered as a guideline to improve total RNA extraction from other bacterial species.

Identification of defense related gene families and their response against powdery and downy mildew infections in Vitis vinifera

Background

Grapevine (Vitis vinifera) productivity has been severely affected by various bacterial, viral and fungal diseases worldwide. When a plant is infected with the pathogen, various defense mechanisms are subsequently activated in plants at various molecular levels. Thus, for substantiating the disease control in an eco-friendly way, it is essential to understand the molecular mechanisms governing pathogen resistance in grapes.

Results

In our study, we performed genome-wide identification of various defensive genes expressed during powdery mildew (PM) and downy mildew (DM) infections in grapevine. Consequently, we identified 6, 21, 2, 5, 3 and 48 genes of Enhanced Disease Susceptibility 1 (EDS1), Non-Race-specific Disease Resistance (NDR1), Phytoalexin deficient 4 (PAD4), Nonexpressor of PR Gene (NPR), Required for Mla-specified resistance (RAR) and Pathogenesis Related (PR), respectively, in the grapevine genome. The phylogenetic study revealed that V. vinifera defensive genes are evolutionarily related to Arabidopsis thaliana. Differential expression analysis resulted in identification of 2, 4, 7, 2, 4, 1 and 7 differentially expressed Nucleotide-binding leucine rich repeat receptor (NLR), EDS1, NDR1, PAD4, NPR, RAR1 and PR respectively against PM infections and 28, 2, 5, 4, 1 and 19 differentially expressed NLR, EDS1, NDR1, NPR, RAR1 and PR respectively against DM infections in V. vinifera. The co-expression study showed the occurrence of closely correlated defensive genes that were expressed during PM and DM stress conditions.

Conclusion

The PM and DM responsive defensive genes found in this study can be characterized in future for impelling studies relaying fungal and oomycete resistance in plants, and the functionally validated genes would then be available for conducting in-planta transgenic gene expression studies for grapes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08081-4.

De novo transcriptome analysis provides insights into formation of in vitro adventitious root from leaf explants of Arnebia euchroma

BACKGROUND

Adventitious root formation is considered a major developmental step during the propagation of difficult to root plants, especially in horticultural crops. Recently, adventitious roots induced through plant tissue culture methods have also been used for production of phytochemicals such as flavonoids, anthocyanins and anthraquinones. It is rather well understood which horticultural species will easily form adventitious roots, but the factors affecting this process at molecular level or regulating the induction process in in vitro conditions are far less known. The present study was conducted to identify transcripts involved in in vitro induction and formation of adventitious roots using Arnebia euchroma leaves at different time points (intact leaf (control), 3 h, 12 h, 24 h, 3 d, 7 d, 10 d and 15 d). A. euchroma is an endangered medicinal Himalayan herb whose root contains red naphthoquinone pigments. These phytoconstituents are widely used as an herbal ingredient in Asian traditional medicine as well as natural colouring agent in food and cosmetics.

RESULTS

A total of 137.93 to 293.76 million raw reads were generated and assembled to 54,587 transcripts with average length of 1512.27 bps and N50 of 2193 bps, respectively. In addition, 50,107 differentially expressed genes were identified and found to be involved in plant hormone signal transduction, cell wall modification and wound induced mitogen activated protein kinase signalling. The data exhibited dominance of auxin responsive (AUXIN RESPONSE FACTOR8, IAA13, GRETCHEN HAGEN3.1) and sucrose translocation (BETA-31 FRUCTOFURANOSIDASE and MONOSACCHARIDE-SENSING protein1) genes during induction phase. In the initiation phase, the expression of LATERAL ORGAN BOUNDARIES DOMAIN16, EXPANSIN-B15, ENDOGLUCANASE25 and LEUCINE-rich repeat EXTENSION-like proteins was increased. During the expression phase, the same transcripts, with exception of LATERAL ORGAN BOUNDARIES DOMAIN16 were identified. Overall, the transcriptomic analysis revealed a similar patterns of genes, however, their expression level varied in subsequent phases of in vitro adventitious root formation in A. euchroma.

CONCLUSION

The results presented here will be helpful in understanding key regulators of in vitro adventitious root development in Arnebia species, which may be deployed in the future for phytochemical production at a commercial scale.

Efficient method for isolation of high-quality RNA from Psidium guajava L. tissues

Acquiring high-quality RNA in sufficient amounts is crucial in plant molecular biology and genetic studies. Several methods for RNA extraction from plants are available in the literature, mainly due to the great biochemical diversity present in each species and tissue, which can complicate or prevent the extraction. Psidium guajava (Myrtaceae family) is a perennial fruit tree of medicinal and economic value; nevertheless, only a few molecular studies are available for the species. One reason is the difficulty in obtaining RNA due to the content of the samples, which are rich in polyphenols, polysaccharides, and secondary metabolites. Furthermore, there are few studies available for the isolation of RNA from guava or Psidium samples, which hampers advances in the study of the genus. Here, quality and yields of RNA isolates were compared using six extraction protocols: two protocols based on the application of cetyltrimethylammonium bromide (CTAB) lysis buffer, one protocol which uses the TRIzol reagent, one which applies guanidine thiocyanate lysis buffer followed by organic phase extraction, and two commercial kits (PureLink RNA Mini Kit and RNeasy Plant Mini Kit). The CTAB-based method provided the highest RNA yields and quality for five different tissues (flower bud, immature leaf, young leaf, mature leaf, and root), genotypes, and stress conditions. For the most efficient protocol, the average yield of RNA from guava leaves was 203.06 μg/g of tissue, and the A₂₆₀/A₂₈₀ and A₂₆₀/A₂₃₀ ratios were 2.1 and 2.2, respectively. RT-qPCR analysis demonstrated that the purity of the samples was sufficient for molecular biology experiments. CTAB-based methods for RNA isolation were found to be the most efficient, providing the highest RNA yields and quality for tissues from P. guajava. Additionally, they were compatible for downstream RNA-based applications, besides being simple and cost-effective.

The first draft genome of Picrorhiza kurrooa, an endangered medicinal herb from Himalayas

Picrorhiza kurrooa is an endangered medicinal herb which is distributed across the Himalayan region at an altitude between 3000–5000 m above mean sea level. The medicinal properties of P. kurrooa are attributed to monoterpenoid picrosides present in leaf, rhizome and root of the plant. However, no genomic information is currently available for P. kurrooa, which limits our understanding about its molecular systems and associated responses. The present study brings the first assembled draft genome of P. kurrooa by using 227 Gb of raw data generated by Illumina and PacBio RS II sequencing platforms. The assembled genome has a size of n = ~ 1.7 Gb with 12,924 scaffolds. Four pronged assembly quality validations studies, including experimentally reported ESTs mapping and directed sequencing of the assembled contigs, confirmed high reliability of the assembly. About 76% of the genome is covered by complex repeats alone. Annotation revealed 24,798 protein coding and 9789 non-coding genes. Using the assembled genome, a total of 710 miRNAs were discovered, many of which were found responsible for molecular response against temperature changes. The miRNAs and targets were validated experimentally. The availability of draft genome sequence will aid in genetic improvement and conservation of P. kurrooa. Also, this study provided an efficient approach for assembling complex genomes while dealing with repeats when regular assemblers failed to progress due to repeats.

Genome-wide transcriptional analysis unveils the molecular basis of organ-specific expression of isosteroidal alkaloids biosynthesis in critically endangered Fritillaria roylei Hook

Fritillaria roylei Hook. is a critically endangered high altitude Himalayan medicinal plant species with rich source of pharmaceutically active structurally diverse steroidal alkaloids. Nevertheless, except few marker compounds, the chemistry of the plant remains unexplored. Therefore, in the current study, transcriptome sequencing efforts were made to elucidate isosteroidal alkaloids biosynthesis by creating first organ-specific genomic resource using bulb, stem, and leaf tissues derived from natural populations of Indian Himalayan region. Overall, 349.9 million high quality paired-end reads obtained using NovaSeq 6000 platform were assembled (de novo) into 82,848 unigenes and 31,061 isoforms. Functional annotation and organ specific differential expression (DE) analysis identified 2488 significant DE transcripts, grouped into three potential sub-clusters (sub-cluster I: 728 transcripts; sub-cluster II: 446 transcripts and Sub-cluster III: 1314 transcripts). Subsequently, pathway enrichment (GO, KEGG) and protein-protein network analysis revealed significantly higher enrichment of phenyl-propanoid and steroid backbone including terpenoid, sesquiterpenoid and triterpenoid biosynthesis in bulb. Additionally, upregulated expression of cytochrome P450, UDP-dependent Glucuronosyltransferase families and key transcription factor families (FAR1, bHLH, GRAS, C2H2, TCP and MYB) suggests 'bulb' as a primary site of MVA mediated isosteroidal alkaloids biosynthesis. The comprehensive elucidation of molecular insights in this study is a first step towards the understanding of isosteroidal alkaloid biosynthesis pathway in F. roylei. Furthermore, key genes and regulators identified here can facilitate metabolic engineering of potential bioactive compounds at industrial scale.

Investigation of long non-coding RNAs as regulatory players of grapevine response to powdery and downy mildew infection

BACKGROUND

Long non-coding RNAs (lncRNAs) are regulatory transcripts of length > 200 nt. Owing to the rapidly progressing RNA-sequencing technologies, lncRNAs are emerging as considerable nodes in the plant antifungal defense networks. Therefore, we investigated their role in Vitis vinifera (grapevine) in response to obligate biotrophic fungal phytopathogens, Erysiphe necator (powdery mildew, PM) and Plasmopara viticola (downy mildew, DM), which impose huge agro-economic burden on grape-growers worldwide.

RESULTS

Using computational approach based on RNA-seq data, 71 PM- and 83 DM-responsive V. vinifera lncRNAs were identified and comprehensively examined for their putative functional roles in plant defense response. V. vinifera protein coding sequences (CDS) were also profiled based on expression levels, and 1037 PM-responsive and 670 DM-responsive CDS were identified. Next, co-expression analysis-based functional annotation revealed their association with gene ontology (GO) terms for 'response to stress', 'response to biotic stimulus', 'immune system process', etc. Further investigation based on analysis of domains, enzyme classification, pathways enrichment, transcription factors (TFs), interactions with microRNAs (miRNAs), and real-time quantitative PCR of lncRNAs and co-expressing CDS pairs suggested their involvement in modulation of basal and specific defense responses such as: Ca²⁺-dependent signaling, cell wall reinforcement, reactive oxygen species metabolism, pathogenesis related proteins accumulation, phytohormonal signal transduction, and secondary metabolism.

CONCLUSIONS

Overall, the identified lncRNAs provide insights into the underlying intricacy of grapevine transcriptional reprogramming/post-transcriptional regulation to delay or seize the living cell-dependent pathogen growth. Therefore, in addition to defense-responsive genes such as TFs, the identified lncRNAs can be further examined and leveraged to candidates for biotechnological improvement/breeding to enhance fungal stress resistance in this susceptible fruit crop of economic and nutritional importance.

Transcriptomic resources for prairie grass (Bromus catharticus): expressed transcripts, tissue-specific genes, and identification and validation of EST-SSR markers

BACKGROUND

Prairie grass (Bromus catharticus) is a typical cool-season forage crop with high biomass production and fast growth rate during winter and spring. However, its genetic research and breeding has remained stagnant due to limited available genomic resources. The aim of this study was to generate large-scale genomic data using high-throughput transcriptome sequencing, and perform a preliminary validation of EST-SSR markers of B. catharticus.

RESULTS

Eleven tissue samples including seeds, leaves, and stems were collected from a new high-yield strain of prairie grass BCS1103. A total of 257,773 unigenes were obtained, of which 193,082 (74.90%) were annotated. Comparison analysis between tissues identified 1803, 3030, and 1570 genes specifically and highly expressed in seed, leaf, and stem, respectively. A total of 37,288 EST-SSRs were identified from unigene sequences, and more than 80,000 primer pairs were designed. We synthesized 420 primer pairs and selected 52 ones with high polymorphisms to estimate genetic diversity and population structure in 24 B. catharticus accessions worldwide. Despite low diversity indicated by an average genetic distance of 0.364, the accessions from South America and Asia and wild accessions showed higher genetic diversity. Moreover, South American accessions showed a pure ancestry, while Asian accessions demonstrated mixed internal relationships, which indicated a different probability of gene flow. Phylogenetic analysis clustered the studied accessions into four clades, being consistent with phenotypic clustering results. Finally, Mantel analysis suggested the total phenotypic variation was mostly contributed by genetic component. Stem diameter, plant height, leaf width, and biomass yield were significantly correlated with genetic data (r > 0.6, P < 0.001), and might be used in the future selection and breeding.

CONCLUSION

A genomic resource was generated that could benefit genetic and taxonomic studies, as well as molecular breeding for B. catharticus and its relatives in the future.

Prickle morphogenesis in rose is coupled with secondary metabolite accumulation and governed by canonical MBW transcriptional complex

Rose is an economically important flowering plant that holds an essential place in cut flower, medicinal, and aromatic industries. The presence of prickles, epidermal outgrowths resembling trichomes, on rose is highly undesirable as these make harvesting and transportation difficult. Attempts were made for generating rose varieties lacking prickles via breeding and natural selections; however, these approaches obtained only chimeric and genetically unstable prickle-less mutants. The alternative way to get rid of prickles is via genetic manipulations, but the molecular mechanisms of prickle initiation and development in rose are almost unexplored. Therefore, the present study was carried out to understand the morphological, molecular, and correlated metabolic changes underlining prickle morphogenesis in a prickle-bearing Rosa hybrida L. cv. "First Red (FR)". The histological and metabolomic analyses at three distinct stages of the prickle morphogenesis, namely, emerging tiny initiating prickles, partially greenish soft prickles, and brownish hard prickles, demonstrated a gradually increasing deposition of phenolic compounds and lignification with development. Corresponding RNAseq analysis revealed an upregulation of the genes involved in secondary metabolism, especially in the phenylpropanoid biosynthetic pathway. A set of genes encoding a transcriptional network similar to the one regulating epidermal cell differentiation leading to phenylpropanoid accumulation and trichome development, was also upregulated. Differential expression of this transcriptional network in prickle-less R. hybrida L. cv. "Himalayan Wonder" compared to prickly FR plants substantiated its involvement in prickle morphogenesis. The results collectively supported the proposition that prickles are evolved from trichomes and provided molecular clues towards engineering prickle-less roses.

SIGNIFICANCE STATEMENT

Prickles, the vasculature less epidermal outgrowths resembling trichomes, are defense organs protecting plants against herbivory. Despite biological significance, the mechanism of prickle morphogenesis remains obscure. Here, we show that like trichomes, prickles accumulate secondary metabolites, especially lignin and flavonoids, during morphogenesis. Cognate transcriptome analysis demonstrated that upregulation of a hormone-regulated transcriptional activation-inhibition network, known to govern trichome morphogenesis, likely triggers the differentiation of epidermal cells to outgrow into prickle.