Transcriptome analysis by Illumina high-throughout paired-end sequencing reveals the complexity of differential gene expression during in vitro plantlet growth and flowering in Amaranthus tricolor L

Comparative Transcriptomic Analysis of the Metabolism of Betalains and Flavonoids in Red Amaranth Hypocotyl under Blue Light and Dark Conditions

Amaranth plants contain abundant betalains and flavonoids. Anthocyanins are important flavonoids; however, they cannot coexist in the same plant with betalains. Blue light influences metabolite synthesis and hypocotyl elongation; accordingly, analyses of its effects on betalain and flavonoid biosynthesis in Amaranthus tricolor may provide insight into the distribution of these plant pigments. We analyzed the betalain and flavonoid content and transcriptome profiles in amaranth hypocotyls under blue light and dark conditions. Furthermore, we analyzed the expression patterns of key genes related to betalains and flavonoids. Amaranth hypocotyls were shorter and redder and showed higher betalain and flavonoid content under blue light than in dark conditions. Key genes involved in the synthesis of betalains and flavonoids were upregulated under blue light. The gene encoding DELLA was also upregulated. These results suggest that blue light favors the synthesis of both betalains and flavonoids via the suppression of bioactive gibberellin and the promotion of DELLA protein accumulation, which also suppresses hypocotyl elongation. The metabolite profiles differed between plants under blue light and dark conditions. These findings improve our understanding of the environmental cues and molecular mechanisms underlying pigment variation in Amaranthus.

Transcriptome Analysis to Understand Salt Stress Regulation Mechanism of Chromohalobacter salexigens ANJ207

Soil salinity is one of the major global issues affecting soil quality and agricultural productivity. The plant growth-promoting halophilic bacteria that can thrive in regions of high salt (NaCl) concentration have the ability to promote the growth of plants in salty environments. In this study, attempts have been made to understand the salinity adaptation of plant growth-promoting moderately halophilic bacteria Chromohalobacter salexigens ANJ207 at the genetic level through transcriptome analysis. In order to identify the stress-responsive genes, the transcriptome sequencing of C. salexigens ANJ207 under different salt concentrations was carried out. Among the 8,936 transcripts obtained, 93 were upregulated while 1,149 were downregulated when the NaCl concentration was increased from 5 to 10%. At 10% NaCl concentration, genes coding for lactate dehydrogenase, catalase, and OsmC-like protein were upregulated. On the other hand, when salinity was increased from 10 to 25%, 1,954 genes were upregulated, while 1,287 were downregulated. At 25% NaCl, genes coding for PNPase, potassium transporter, aconitase, excinuclease subunit ABC, and transposase were found to be upregulated. The quantitative real-time PCR analysis showed an increase in the transcript of genes related to the biosynthesis of glycine betaine coline genes (gbcA, gbcB, and L-pro) and in the transcript of genes related to the uptake of glycine betaine (OpuAC, OpuAA, and OpuAB). The transcription of the genes involved in the biosynthesis of L-hydroxyproline (proD and proS) and one stress response proteolysis gene for periplasmic membrane stress sensing (serP) were also found to be increased. The presence of genes for various compatible solutes and their increase in expression at the high salt concentration indicated that a coordinated contribution by various compatible solutes might be responsible for salinity adaptation in ANJ207. The investigation provides new insights into the functional roles of various genes involved in salt stress tolerance and oxidative stress tolerance produced by high salt concentration in ANJ207 and further support the notion regarding the utilization of bacterium and their gene(s) in ameliorating salinity problem in agriculture.

Analysis of Biosynthetic Gene Clusters, Secretory, and Antimicrobial Peptides Reveals Environmental Suitability of Exiguobacterium profundum PHM11

Halotolerant bacteria produce a wide range of bioactive compounds with important applications in agriculture for abiotic stress amelioration and plant growth promotion. In the present study, 17 biosynthetic gene clusters (BGCs) were identified in Exiguobacterium profundum PHM11 belonging to saccharides, desmotamide, pseudaminic acid, dipeptide aldehydes, and terpene biosynthetic pathways representing approximately one-sixth of genomes. The terpene biosynthetic pathway was conserved in Exiguobacterium spp. while the E. profundum PHM11 genome confirms the presence of the 1-deoxy-d-xylulose 5-phosphate (DXP) pathway for the isopentenyl diphosphate (IPP) synthesis. Further, 2,877 signal peptides (SPs) were identified using the PrediSi server, out of which 592 proteins were prophesied for the secretion having a transmembrane helix (TMH). In addition, antimicrobial peptides (AMPs) were also identified using BAGEL4. The transcriptome analysis of PHM11 under salt stress reveals the differential expression of putative secretion and transporter genes having SPs and TMH. Priming of the rice, wheat and maize seeds with PHM11 under salt stress led to improvement in the root length, root diameters, surface area, number of links and forks, and shoot length. The study shows that the presence of BGCs, SPs, and secretion proteins constituting TMH and AMPs provides superior competitiveness in the environment and make E. profundum PHM11 a suitable candidate for plant growth promotion under salt stress.

RNA-sequencing analysis reveals betalains metabolism in the leaf of Amaranthus tricolor L

Amaranth plants contain large amounts of betalains, including betaxanthins and betacyanins. Amaranthin is a betacyanin, and its molecular structure and associated metabolic pathway differ from those of betanin in beet plants. The chlorophyll, carotenoid, betalain, and flavonoid contents in amaranth leaves were analyzed. The abundance of betalain, betacyanin, and betaxanthin was 2-5-fold higher in the red leaf sectors than in the green leaf sectors. Moreover, a transcriptome database was constructed for the red and green sectors of amaranth leaves harvested from 30-day-old seedlings. 22 unigenes were selected to analyze the expression profiles in the two leaf sectors. The RNA-sequencing data indicated that many unigenes are involved in betalain metabolic pathways. The potential relationships between diverse metabolic pathways and betalain metabolism were analyzed. The validation of the expression of 22 selected unigenes in a qRT-PCR assay revealed the genes that were differentially expressed in the two leaf sectors. Betalains were biosynthesized in specific tissues of the red sectors of amaranth leaves. Almost all of the genes related to betalain metabolism were identified in the transcriptome database, and the expression profiles were different between the red sectors and green sectors in the leaf. Amaranth plants consist of diverse metabolic pathways, and the betalain metabolic pathway is linked to a group of other metabolic pathways.

Reference genes for RT-qPCR normalisation in different tissues, developmental stages and stress conditions of amaranth

Studies of gene expression are very important for the identification of genes that participate in different biological processes. Currently, reverse transcription quantitative real-time PCR (RT-qPCR) is a high-throughput, sensitive and widely used method for gene expression analysis. Nevertheless, RT-qPCR requires precise normalisation of data to avoid the misinterpretation of experimental data. In this sense, the selection of reference genes is critical for gene expression analysis. At this time, several studies focus on the selection of reference genes in several species. However, the identification and validation of reference genes for the normalisation of RT-qPCR have not been described in amaranth. A set of seven housekeeping genes were analysed using RT-qPCR, to determine the most stable reference genes in amaranth for normalisation of gene expression analysis. Transcript stability and gene expression level of candidate reference genes were analysed in different tissues, at different developmental stages and under different types of stress. The data were compared using the geNorm, NormFinder and Bestkeeper statistical methods. The reference genes optimum for normalisation of data varied with respect to treatment. The results indicate that AhyMDH, AhyGAPDH, AhyEF-1α and AhyACT would be optimum for accurate normalisation of experimental data, when all treatment are analysed in the same experiment. This study presents the most stable reference genes for normalisation of gene expression analysis in amaranth, which will contribute significantly to future gene studies of this species.

Comprehensive transcriptome profiling of Salvia miltiorrhiza for discovery of genes associated with the biosynthesis of tanshinones and phenolic acids

Tanshinones and phenolic acids are crucial bioactive compounds biosynthesized in Salvia miltiorrhiza. Methyl jasmonate (MeJA) is an effective elicitor to enhance the production of phenolic acids and tanshinones simultaneously, while yeast extract (YE) is used as a biotic elicitor that only induce tanshinones accumulation. However, little was known about the different molecular mechanism. To identify the downstream and regulatory genes involved in tanshinone and phenolic acid biosynthesis, we conducted comprehensive transcriptome profiling of S. miltiorrhiza hairy roots treated with either MeJA or YE. Total 55588 unigenes were assembled from about 1.72 billion clean reads, of which 42458 unigenes (76.4%) were successfully annotated. The expression patterns of 19 selected genes in the significantly upregulated unigenes were verified by quantitative real-time PCR. The candidate downstream genes and other cytochrome P450s involved in the late steps of tanshinone and phenolic acid biosynthesis pathways were screened from the RNA-seq dataset based on co-expression pattern analysis with specific biosynthetic genes. Additionally, 375 transcription factors were identified to exhibit a significant up-regulated expression pattern in response to induction. This study can provide us a valuable gene resource for elucidating the molecular mechanism of tanshinones and phenolic acids biosynthesis in hairy roots of S. miltiorrhiza.

Single-molecule sequencing and Hi-C-based proximity-guided assembly of amaranth (Amaranthus hypochondriacus) chromosomes provide insights into genome evolution

Background

Amaranth (Amaranthus hypochondriacus) was a food staple among the ancient civilizations of Central and South America that has recently received increased attention due to the high nutritional value of the seeds, with the potential to help alleviate malnutrition and food security concerns, particularly in arid and semiarid regions of the developing world. Here, we present a reference-quality assembly of the amaranth genome which will assist the agronomic development of the species.

Results

Utilizing single-molecule, real-time sequencing (Pacific Biosciences) and chromatin interaction mapping (Hi-C) to close assembly gaps and scaffold contigs, respectively, we improved our previously reported Illumina-based assembly to produce a chromosome-scale assembly with a scaffold N50 of 24.4 Mb. The 16 largest scaffolds contain 98% of the assembly and likely represent the haploid chromosomes (n = 16). To demonstrate the accuracy and utility of this approach, we produced physical and genetic maps and identified candidate genes for the betalain pigmentation pathway. The chromosome-scale assembly facilitated a genome-wide syntenic comparison of amaranth with other Amaranthaceae species, revealing chromosome loss and fusion events in amaranth that explain the reduction from the ancestral haploid chromosome number (n = 18) for a tetraploid member of the Amaranthaceae.

Conclusions

The assembly method reported here minimizes cost by relying primarily on short-read technology and is one of the first reported uses of in vivo Hi-C for assembly of a plant genome. Our analyses implicate chromosome loss and fusion as major evolutionary events in the 2n = 32 amaranths and clearly establish the homoeologous relationship among most of the subgenome chromosomes, which will facilitate future investigations of intragenomic changes that occurred post polyploidization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0412-4) contains supplementary material, which is available to authorized users.

Utility of next-generation RNA-sequencing in identifying chimeric transcription involving human endogenous retroviruses

Several studies have shown that human endogenous retroviruses and endogenous retrovirus-like repeats (here collectively HERVs) impose direct regulation on human genes through enhancer and promoter motifs present in their long terminal repeats (LTRs). Although chimeric transcription in which novel gene isoforms containing retroviral and human sequence are transcribed from viral promoters are commonly associated with disease, regulation by HERVs is beneficial in other settings; for example, in human testis chimeric isoforms of TP63 induced by an ERV9 LTR protect the male germ line upon DNA damage by inducing apoptosis, whereas in the human globin locus the γ- and β-globin switch during normal hematopoiesis is mediated by complex interactions of an ERV9 LTR and surrounding human sequence. The advent of deep sequencing or next-generation sequencing (NGS) has revolutionized the way researchers solve important scientific questions and develop novel hypotheses in relation to human genome regulation. We recently applied next-generation paired-end RNA-sequencing (RNA-seq) together with chromatin immunoprecipitation with sequencing (ChIP-seq) to examine ERV9 chimeric transcription in human reference cell lines from Encyclopedia of DNA Elements (ENCODE). This led to the discovery of advanced regulation mechanisms by ERV9s and other HERVs across numerous human loci including transcription of large gene-unannotated genomic regions, as well as cooperative regulation by multiple HERVs and non-LTR repeats such as Alu elements. In this article, well-established examples of human gene regulation by HERVs are reviewed followed by a description of paired-end RNA-seq, and its application in identifying chimeric transcription genome-widely. Based on integrative analyses of RNA-seq and ChIP-seq, data we then present novel examples of regulation by ERV9s of tumor suppressor genes CADM2 and SEMA3A, as well as transcription of an unannotated region. Taken together, this article highlights the high suitability of contemporary sequencing methods in future analyses of human biology in relation to evolutionary acquired retroviruses in the human genome.

De novo transcriptomic analysis of the female and male adults of the blood fluke Schistosoma turkestanicum

Background

Schistosoma turkestanicum is a parasite of considerable veterinary importance as an agent of animal schistosomiasis in many countries, including China. The S. turkestanicum cercariae can also infect humans, causing cercarial dermatitis in many countries and regions of the world. In spite of its significance as a pathogen of animals and humans, there is little transcriptomic and genomic data in the public databases.

Methods

Herein, we performed the transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females of S. turkestanicum and de novo transcriptome assembly.

Results

Approximately 81.1 (female) and 80.5 (male) million high-quality clean reads were obtained and then 29,526 (female) and 41,346 (male) unigenes were assembled. A total of 34,624 unigenes were produced from S. turkestanicum females and males, with an average length of 878 nucleotides (nt) and N50 of 1480 nt. Of these unigenes, 25,158 (72.7 %) were annotated by blast searches against the NCBI non-redundant protein database. Among these, 21,995 (63.5 %), 22,189 (64.1 %) and 13,754 (39.7 %) of the unigenes had significant similarity in the NCBI non-redundant protein (NR), non-redundant nucleotide (NT) and Swiss-Prot databases, respectively. In addition, 3150 unigenes were identified to be expressed specifically in females and 1014 unigenes were identified to be expressed specifically in males. Interestingly, several pathways associated with gonadal development and sex maintenance were found, including the Wnt signaling pathway (103; 2 %) and progesterone-mediated oocyte maturation (77; 1.5 %).

Conclusions

The present study characterized and compared the transcriptomes of adult female and male blood fluke, S. turkestanicum. These results will not only serve as valuable resources for future functional genomics studies to understand the molecular aspects of S. turkestanicum, but also will provide essential information for ongoing whole genome sequencing efforts on this pathogenic blood fluke.

Comprehensive transcriptome analysis of Crocus sativus for discovery and expression of genes involved in apocarotenoid biosynthesis

Background

Crocus sativus stigmas form rich source of apocarotenoids like crocin, picrocrocin and saffranal which besides imparting color, flavour and aroma to saffron spice also have tremendous pharmacological properties. Inspite of their importance, the biosynthetic pathway of Crocus apocarotenoids is not fully elucidated. Moreover, the mechanism of their stigma specific accumulation remains unknown. Therefore, deep transcriptome sequencing of Crocus stigma and rest of the flower tissue was done to identify the genes and transcriptional regulators involved in the biosynthesis of these compounds.

Results

Transcriptome of stigma and rest of the flower tissue was sequenced using Illumina Genome Analyzer IIx platform which generated 64,604,402 flower and 51,350,714 stigma reads. Sequences were assembled de novo using trinity resulting in 64,438 transcripts which were classified into 32,204 unigenes comprising of 9853 clusters and 22,351 singletons. A comprehensive functional annotation and gene ontology (GO) analysis was carried out. 58.5 % of the transcripts showed similarity to sequences present in public databases while rest could be specific to Crocus. 5789 transcripts showed similarity to transcription factors representing 76 families out of which Myb family was most abundant. Many genes involved in carotenoid/apocarotenoid pathway were identified for the first time in this study which includes zeta-carotene isomerase and desaturase, carotenoid isomerase and lycopene epsilon-cyclase. GO analysis showed that the predominant classes in biological process category include metabolic process followed by cellular process and primary metabolic process. KEGG mapping analysis indicated that pathways involved in ribosome, carbon and starch and sucrose metabolism were highly represented. Differential expression analysis indicated that key carotenoid/apocarotenoid pathway genes including phytoene synthase, phytoene desaturase and carotenoid cleavage dioxygenase 2 are enriched in stigma thereby providing molecular proof for stigma to be the site of apocarotenoid biosynthesis.

Conclusions

This data would provide a rich source for understanding the carotenoid/apocarotenoid metabolism in Crocus. The database would also help in investigating many questions related to saffron biology including flower development.

Electronic supplementary material

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

De novo assembly, transcriptome characterization, lignin accumulation, and anatomic characteristics: novel insights into lignin biosynthesis during celery leaf development

Celery of the family Apiaceae is a biennial herb that is cultivated and consumed worldwide. Lignin is essential for cell wall structural integrity, stem strength, water transport, mechanical support, and plant pathogen defense. This study discussed the mechanism of lignin formation at different stages of celery development. The transcriptome profile, lignin distribution, anatomical characteristics, and expression profile of leaves at three stages were analyzed. Regulating lignin synthesis in celery growth development has a significant economic value. Celery leaves at three stages were collected, and Illumina paired-end sequencing technology was used to analyze large-scale transcriptome sequences. From Stage 1 to 3, the collenchyma and vascular bundles in the petioles and leaf blades thickened and expanded, whereas the phloem and the xylem extensively developed. Spongy and palisade mesophyll tissues further developed and were tightly arranged. Lignin accumulation increased in the petioles and the mesophyll (palisade and spongy), and the xylem showed strong lignification. Lignin accumulation in different tissues and at different stages of celery development coincides with the anatomic characteristics and transcript levels of genes involved in lignin biosynthesis. Identifying the genes that encode lignin biosynthesis-related enzymes accompanied by lignin distribution may help elucidate the regulatory mechanisms of lignin biosynthesis in celery.