Transcriptome analysis for molecular landscaping of genes controlling diterpene andrographolide biosynthesis in Andrographis paniculata (Burm . f.) Nees.

Kalmegh [Andrographis paniculata (Burm. f.) Nees.] is one of the essential medicinal plants due to an important terpenoid, i.e. andrographolide possesses immense therapeutic and pharmacological uses. The experiment was performed to elucidate the expression of candidate genes associated with andrographolide biosynthesis. Based on results obtained in chromatography for andrographolide content analysis of six genotypes, two contrast genotypes, i.e. IC-520361 (maximum andrographolide content-2.33%) and Anand Local (lowest andrographolide content-1.01%) were selected for the transcriptome analysis. A total of 1.04 Gb of raw data were produced using MiSeq Illumina platform, in which IC 520361 generated 645 million base pairs sequence along with 4,524,251 raw reads and Anand Local produced 419 million base pairs sequence along with 3,021,316 raw reads. The combined assembly of high quality reads generated for both the samples had 33,247,454 bp of total assembled bases and 38,292 of transcripts. The GC percent of assembled transcripts was 44.79%, an average read length was 800 bp and N50 value was 1186 bp. Species-specific distribution using BLAST X (Nr), showed the highest Blast hits with Sesamum indicum. Out of 23,346 transcripts, 87% of transcripts annotated in UniProt KB (Universal Protein Resource KnowledgeBase) database and only 0.21% of transcripts were annotated in TAIR (The Arabidopsis Information Resources). Biological processes gene ontology classified based on Blast2GO showed, out of 6853 transcripts, 1370 of transcripts were represented by terpenoid biosynthetic pathway, which involved in secondary metabolite andrographolide biosynthesis. The heat map showed 1016 transcripts were differentially expressed between two kalmegh genotypes, in which nine important differentially expressed transcripts related to MEP (2C methyl-d-erythritol 4-phosphate) and MVA (Mevalonic acid) andrographolide biosynthesis pathways such as, geranyl diphosphate synthase small subunit, Isopentenyl-diphosphate delta-isomerase i-like, 4, 13-hydroxy-3-methylglutaryl-coenzyme a reductase etc. were upregulated in IC 520361 as compared to Anand Local, which were validated through RT-qPCR. The highest expression of gene 13-hydroxy-3-methylglutaryl-coenzyme a reductase (HMGR) was reported, which is responsible for accumulation of andrographolide in leaf. This comparative transcriptome analysis confirmed the expression level of genes were higher in accession IC 520361 as compare to Anand Local related to andrographolide biosynthesis pathways i.e. MEP and MVA. These up-regulated genes could be over-expressed to enhance the andrographolide content using genetic engineering of these metabolic pathways. It will also give an idea to the breeder for development of molecular markers for direct screening of the genotypes.

Transcriptomic profiling of developing fiber in levant cotton (Gossypium herbaceum L.).

Cotton (Gossypium spp.) is an imperative economic crop of the globe due to its natural textile fiber. Molecular mechanisms of fiber development have been greatly revealed in allotetraploid cotton but remained unexplored in Gossypium herbaceum. G. herbaceum can withstand the rigors of nature like drought and pests but produce coarse lint. This undesirable characteristic strongly needs the knowledge of fiber development at molecular basis. The present study reported the transcriptome sequence of the developing fiber of G. herbaceum on pyrosequencing and its analysis. About 1.38 million raw and 1.12 million quality trimmed reads were obtained followed by de novo assembly-generated 20,125 unigenes containing 14,882 coding sequences (CDs). BLASTx-based test of homology indicated that A1-derived transcripts shared a high similarity with Gossypium arboreum (A2). Functional annotation of the CDs using the UniProt categorized them into biological processes, cellular components, and molecular function, COG classification showed that a large number of CDs have significant homology in COG database (6215 CDs), and mapping of CDs with Kyoto Encyclopedia of Genes and Genomes (KEGG) database generated 200 pathways ultimately showing predominant engagement in the fiber development process. Transcription factors were predicted by comparison with Plant Transcription Factor Database, and their differential expression between stages exposed their important regulatory role in fiber development. Differential expression analysis based on reads per kilobase of transcript per million mapped reads (RPKM) value revealed activities of specific gene related to carbohydrate and lipid synthesis, carbon metabolism, energy metabolism, signal transduction, etc., at four stages of fiber development, and was validated by qPCR. Overall, this study will help as a valuable foundation for diploid cotton fiber improvement.

De novo Transcriptome Sequencing to Dissect Candidate Genes Associated with Pearl Millet-Downy Mildew (Sclerospora graminicola Sacc.) Interaction.

Understanding the plant-pathogen interactions is of utmost importance to design strategies for minimizing the economic deficits caused by pathogens in crops. With an aim to identify genes underlying resistance to downy mildew, a major disease responsible for productivity loss in pearl millet, transcriptome analysis was performed in downy mildew resistant and susceptible genotypes upon infection and control on 454 Roche NGS platform. A total of ~685 Mb data was obtained with 1 575 290 raw reads. The raw reads were pre-processed into high-quality (HQ) reads making to ~82% with an average of 427 bases. The assembly was optimized using four assemblers viz. Newbler, MIRA, CLC and Trinity, out of which MIRA with a total of 14.10 Mb and 90118 transcripts proved to be the best for assembling reads. Differential expression analysis depicted 1396 and 936 and 1000 and 1591 transcripts up and down regulated in resistant inoculated/resistant control and susceptible inoculated/susceptible control respectively with a common of 3644 transcripts. The pathways for secondary metabolism, specifically the phenylpropanoid pathway was up-regulated in resistant genotype. Transcripts up-regulated as a part of defense response included classes of R genes, PR proteins, HR induced proteins and plant hormonal signaling transduction proteins. The transcripts for skp1 protein, purothionin, V type proton ATPase were found to have the highest expression in resistant genotype. Ten transcripts, selected on the basis of their involvement in defense mechanism were validated with qRT-PCR and showed positive co-relation with transcriptome data. Transcriptome analysis evoked potentials of hypersensitive response and systemic acquired resistance as possible mechanism operating in defense mechanism in pearl millet against downy mildew infection.

Validation of dbEST-SSRs and transferability of some other solanaceous species SSR in ashwagandha [Withania Somnifera (L.) Dunal].

Cross-species transferability and expressed sequence tags (ESTs) in public databases are cost-effective means for developing simple sequence repeats (SSRs) for less-studied species like medicinal plants. In this study, 11 EST-SSR markers developed from 742 available ESTs of Withania Somnifera EST sequences and 95 SSR primer pairs derived from other solanaceous crops (tomato, eggplant, chili, and tobacco) were utilized for their amplification and validation. Out of 11, 10 EST-SSRs showed good amplification quality and produced 13 loci with a product size ranging between 167 and 291 bp. Similarly, of the 95 cross-genera SSR loci assayed, 20 (21 %) markers showed the transferability of 5, 27, 32, and 14.2 % from eggplant, chili, tomato, and tobacco, respectively, to ashwagandha. In toto, these 30 SSR markers reported here will be valuable resources and may be applicable for the analysis of intra- and inter-specific genetic diversity in ashwagandha for which till date no information about SSR is available.