De novo and comparative transcriptome analysis of cultivated and wild spinach

Effect of differences in light source environment on transcriptome of leaf lettuce (Lactuca sativa L.) to optimize cultivation conditions

When used in closed-type plant factories, light-emitting diode (LED) illumination systems have the particular advantages of low heat emission and high luminous efficiency. The effects of illumination quality and intensity on the growth and morphogenesis of many plant species have been examined, but improvements are needed to optimize the illumination systems for better plant products with lower resource investments. In particular, new strategies are needed to reduce the wastage of plant products related to leaf senescence, and to better control the ingredients and appearance of leafy vegetables. Although the quality of light is often altered to change the characteristics of plant products, the transcriptional status underlying the physiological responses of plants to light has not been established. Herein, we performed a comprehensive gene expression analysis using RNA-sequencing to determine how red, blue, and red/blue LEDs and fluorescent light sources affect transcriptome involved in the leaf aging of leaf lettuce. The RNA-sequencing profiling revealed clear differences in the transcriptome between young and old leaves. Red LED light caused large variation between the two age classes, while a pure or mixed blue LED light spectrum induced fewer transcriptome differences between young and old leaves. Collectively, the expression levels of genes that showed homology with those of other model organisms provide a detailed physiological overview, incorporating such characteristics as the senescence, nutrient deficiency, and anthocyanin synthesis of the leaf lettuce plants. Our findings suggest that transcriptome profiles of leaf lettuce grown under different light sources provide helpful information to achieve better growth conditions for marketable and efficient green-vegetable production, with improved wastage control and efficient nutrient inputs.

Comparative transcriptome mining for terpenoid biosynthetic pathway genes in wild and cultivated species of Plantago

Plantagos are important economical and medicinal plants that possess several bioactive secondary metabolites, such as phenolics, iridoids, triterpenes, and alkaloids. Triterpenoids are the ubiquitous and dynamic secondary metabolites that are deployed by plants for chemical interactions and protection under biotic/abiotic stress. Plantago ovata, a cultivated species, is the source of psyllium, while Plantago major, a wild species, has significant therapeutic potential. Wild species are considered more tolerant to stressful conditions in comparison to their cultivated allies. In view of this, the present study aimed to decipher the terpenoid biosynthetic pathway operative in P. ovata and P. major using a comparative transcriptomics approach. Majority of terpenoid biosynthetic genes were observed as upregulated in P. major including rate limiting genes of MVA (HMGR) and MEP (DXR) pathways and genes (α-AS, BAS, SM, and CYP716) involved in ursolic acid biosynthesis, an important triterpenoid prevalent in Plantago species. The HPLC output further confirmed the higher concentration of ursolic acid in P. major as compared to P. ovata leaf samples, respectively. In addition to terpenoid biosynthesis, KEGG annotation revealed the involvement of differentially expressed unigenes in several metabolic pathways, aminoacyl-tRNA biosynthesis, biosynthesis of antibiotics, and biosynthesis of secondary metabolites. MYB was found as the most abundant transcription factor family in Plantago transcriptome. We have been able to generate valuable information which can help in improving terpenoid production in Plantago. Additionally, the present study has laid a strong foundation for deciphering other important metabolic pathways in Plantago.

Comparative chloroplast genome analyses of cultivated spinach and two wild progenitors shed light on the phylogenetic relationships and variation

Spinacia is a genus of important leafy vegetable crops worldwide and includes cultivated Spinacia oleracea and two wild progenitors, Spinacia turkestanica and Spinacia tetrandra. However, the chloroplast genomes of the two wild progenitors remain unpublished, limiting our knowledge of chloroplast genome evolution among these three Spinacia species. Here, we reported the complete chloroplast genomes of S. oleracea, S. turkestanica, and S. tetrandra obtained via Illumina sequencing. The three chloroplast genomes exhibited a typical quadripartite structure and were 150,739, 150,747, and 150,680 bp in size, respectively. Only three variants were identified between S. oleracea and S. turkestanica, whereas 690 variants were obtained between S. oleracea and S. tetrandra, strongly demonstrating the close relationship between S. turkestanica and S. oleracea. This was further supported by phylogenetic analysis. We reported a comprehensive variant dataset including 503 SNPs and 83 Indels using 85 Spinacia accessions containing 61 S. oleracea, 16 S. turkestanica, and eight S. tetrandra accessions. Thirteen S. oleracea accessions were derived through introgression from S. turkestanica that acts as the maternal parent. Together, these results provide a valuable resource for spinach breeding programs and improve our understanding of the phylogenetic relationships within Amaranthaceae.

Genetic Diversity, Structure, and Selective Sweeps in Spinacia turkestanica Associated With the Domestication of Cultivated Spinach

Genotype-by-sequencing (GBS) was used to explore the genetic diversity and structure of Spinacia turkestanica, and the selective sweeps involved in domestication of cultivated spinach, S. oleracea, from S. turkestanica. A total 7,065 single nucleotide polymorphisms (SNPs) generated for 16 Spinacia oleracea and 76 S. turkestanica accessions placed the S. oleracea accessions in one group, Q1, and the 76 S. turkestanica accessions, which originated from Central Asia, in two distinct groups, Q2 and Q3. The Q2 group shared greater genetic identity with the S. oleracea accessions, Q1, than the Q3 S. turkestanica group. Likewise, the S. oleracea Q1 group had a smaller Fst (0.008) with the Q2 group than with the Q3 group (Fst = 0.012), and a greater gene flow (Nm = 30.13) with the Q2 group than with the Q3 group (Nm = 21.83). The Q2 accessions originated primarily from Uzbekistan while the Q3 accessions originated mostly from Tajikistan. The Zarafshan Mountain Range appears to have served as a physical barrier that largely separated members of the Q2 and Q3 groups of S. turkestanica. Accessions with admixtures of Q2 and Q3 were collected primarily from lower elevations at the southern end of the Zarafshan Mountain Range in Uzbekistan. Selective sweep regions identified at 32, 49, and 52 Mb on chromosomes 1, 2, and 3, respectively, appear to have played a vital role in the domestication of S. oleracea as they are correlated with important domestication traits, including day length sensitivity for bolting (flowering). High XP-CLR scores at the 52 Mb genomic region of chromosome three suggest that a selective sweep at this region was responsible for early differentiation of S. turkestanica into two groups in Central Asia.

Digital RNA-seq transcriptome plus tissue anatomy analyses reveal the developmental mechanism of the calabash-shaped root in Tetrastigma hemsleyanum

Tetrastigma hemsleyanum Diels & Gilg ex Diels is a liana plant with promising medicinal and ornamental values. Its calabash-shaped roots (CRs) are served as a traditional Chinese herb. However, it takes a long growth period to form CRs. In this study, three types of architectural roots, including fibrous roots (FRs), bar-shaped roots (BRs) and CRs, were employed as materials, and the characteristics of histo-anatomy and digital RNA-seq transcriptome profiles were analyzed. Among the three types of roots, the vascular bundles in FRs were intact, while some of the vascular bundles degenerated in BRs, and only few traces of vascular bundles existed in CRs. Meanwhile, no obvious cell inclusions were found in the cytoplasm of FRs, while a few inclusions were found in BRs, and abundant inclusions were detected in CRs, which might be the main source of medicinal components in roots. The transcriptome profiles and qRT-PCR validation indicated that seven upregulated genes, encoding xyloglucan glycosyltransferase, ACC oxidase, CYP711A1, SHORT-ROOT transcript factor, galacturonosyltransferas, WAT1 and WRKY, and two downregulated genes, encoding LRR receptor-like serine/threonine-protein kinase and CYP83B1, were probably involved in the formation and development of CRs. In addition, Gene Ontology terms of intrinsic component of membrane, integral component of membrane, cell periphery, membrane part, plasma membrane, membrane, intrinsic component of plasma membrane, cellular chemical homeostasis and plasma membrane part were probably related to the formation of CRs. Kyoto Encyclopedia of Genes and Genomes pathways related to the development of CRs probably included MAPK signaling pathway-plant, plant hormone signal transduction and circadian rhythm-plant. Our finding suggested a probable mode for the formation of CRs.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Spinacia oleracea under Abiotic Stress

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is an accurate and convenient method for mRNA quantification. Selection of optimal reference gene(s) is an important step in RT-qPCR experiments. However, the stability of housekeeping genes in spinach (Spinacia oleracea) under various abiotic stresses is unclear. Evaluating the stability of candidate genes and determining the optimal gene(s) for normalization of gene expression in spinach are necessary to investigate the gene expression patterns during development and stress response. In this study, ten housekeeping genes, 18S ribosomal RNA (18S rRNA), actin, ADP ribosylation factor (ARF), cytochrome c oxidase subunit 5C (COX), cyclophilin (CYP), elongation factor 1-alpha (EF1α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone H3 (H3), 50S ribosomal protein L2 (RPL2), and tubulin alpha chain (TUBα) from spinach, were selected as candidates in roots, stems, leaves, flowers, and seedlings in response to high temperature, CdCl₂, NaCl, NaHCO₃, and Na₂CO₃ stresses. The expression of these genes was quantified by RT-qPCR and evaluated by NormFinder, BestKeeper, and geNorm. 18S rRNA, actin, ARF, COX, CYP, EF1α, GAPDH, H3, and RPL2 were detected as optimal reference genes for gene expression analysis of different organs and stress responses. The results were further confirmed by the expression pattern normalized with different reference genes of two heat-responsive genes. Here, we optimized the detection method of the gene expression pattern in spinach. Our results provide the optimal candidate reference genes which were crucial for RT-qPCR analysis.

Variability in somatic embryo-forming capacity of spinach

High variability in somatic embryo (SE)-forming capacity has previously been observed in several spinach cultivars. Such variability frequently accounted for more variation in embryogenic response of the explants than the factor being investigated. Hence, the variability in embryogenic capacity was examined in the present study at both the population and the single-seedling level, using seeds of spinach cultivar Matador obtained from nine European seed companies. Seed population obtained from Slovenia (Sl) was superior to others, with the highest regeneration frequency (100%) and the highest mean SE number (14.4). A total of 82% of these seedlings had 80–100% of regenerating explants, while in populations with intermediate embryogenic capacity approximately 40% of seedlings had 20–60% of regenerating explants. The explants from the majority of seedlings (52–100%) in the least responsive populations were irresponsive. Furthermore, the explants from Sl seedlings regenerated from 10–20 (43.5%) up to > 20 (27.6%) SEs on average, while the explants from the majority of seedlings belonging to other populations regenerated 1–10 SEs. The present study strongly indicates that the variability of plant material must not be overlooked, because choosing more responsive individuals for one treatment and less responsive ones for another may lead to misinterpretation of the data.

Transcriptomic evidence that insulin signalling pathway regulates the ageing of subterranean termite castes

Insulin is a protein hormone that controls the metabolism of sugar, fat and protein via signal transduction in cells, influencing growth and developmental processes such as reproduction and ageing. From nematodes to fruit flies, rodents and other animals, glucose signalling mechanisms are highly conserved. Reproductive termites (queens and kings) exhibit an extraordinarily long lifespan relative to non-reproductive individuals such as workers, despite being generated from the same genome, thus providing a unique model for the investigation of longevity. The key reason for this molecular mechanism, however, remains unclear. To clarify the molecular mechanism underlying this phenomenon, we sequenced the transcriptomes of the primary kings (PKs), primary queens (PQs), male (WMs) and female (WFs) workers of the lower subterranean termite Reticulitermes chinensis. We performed RNA sequencing and identified 33 insulin signalling pathway-related genes in R. chinensis. RT-qPCR analyses revealed that EIF4E and RPS6 genes were highly expressed in WMs and WFs workers, while mTOR expression was lower in PKs and PQs than in WMs and WFs. PQs and PKs exhibited lower expression of akt2-a than female workers. As the highly conserved insulin signalling pathway can significantly prolong the healthspan and lifespan, so we infer that the insulin signalling pathway regulates ageing in the subterranean termite R. chinensis. Further studies are recommended to reveal the biological function of insulin signalling pathway-related genes in the survival of termites to provide new insights into biomolecular homeostasis maintenance and its relationship to remarkable longevity.

Application of cucumber mosaic virus to efficient induction and long-term maintenance of virus-induced gene silencing in spinach

Virus-induced gene silencing (VIGS) is a useful tool for functional genomics in plants. In this study, we tried to apply cucumber mosaic virus (CMV) to efficient induction of VIGS in spinach. Although VIGS for spinach had been previously developed based on two viruses (beet curly top virus and tobacco rattle virus), they still have some problems with systemic movement and long-term maintenance of VIGS in spinach. Although ordinary CMV strains infect spinach inducing distinct mosaic symptoms, using a CMV pseudorecombinant, we can modify the viral pathogenicity to attenuate viral symptoms that may mask the silencing phenotype. We here successfully demonstrated the viral ability to silence the phytoene desaturase (PDS) and the dihydroflavonol 4-reductase (DFR) genes in spinach. Because CMV could quickly induce VIGS even at 7-10 days postinoculation and the virus did not disappear even at the flowering stage, this CMV-based VIGS system would contribute to functional genomics in spinach and especially to the elucidation of molecular mechanisms for some properties unique to spinach such as plasticity of sex expression; the CMV-induced VIGS can last until the flowering stage after the virus was inoculated onto the seedling.

SpinachBase: a central portal for spinach genomics

Spinach (Spinacia oleracea L.) is a nutritious vegetable enriched with many essential minerals and vitamins. A reference spinach genome has been recently released, and additional spinach genomic resources are being rapidly developed. Therefore, there is an urgent need of a central database to store, query, analyze and integrate various resources of spinach genomic data. To this end, we developed SpinachBase (http://spinachbase.org), which provides centralized public accesses to genomic data as well as analytical tools to assist research and breeding in spinach. The database currently stores the spinach reference genome sequence, and sequences and comprehensive functional annotations of protein-coding genes predicted from the genome. The database also contains gene expression profiles derived from RNA-Seq experiments as well as highly co-expressed genes and genetic variants called from transcriptome sequences of 120 cultivated and wild Spinacia accessions. Biochemical pathways have been predicted from spinach protein-coding genes and are available through a pathway database (SpinachCyc) within SpinachBase. SpinachBase provides a suite of analysis and visualization tools including a genome browser, sequence similarity searches with BLAST, functional enrichment and functional classification analyses and functions to query and retrieve gene sequences and annotations.

Transcriptome analysis reveals plasticity in gene regulation due to environmental cues in Primula sikkimensis, a high altitude plant species

BACKGROUND

Studying plasticity in gene expression in natural systems is crucial, for predicting and managing the effects of climate change on plant species. To understand the contribution of gene expression level variations to abiotic stress compensation in a Himalaya plant (Primula sikkimensis), we carried out a transplant experiment within (Ambient), and beyond (Below Ambient and Above Ambient) the altitudinal range limit of species. We sequenced nine transcriptomes (three each from each altitudinal range condition) using Illumina sequencing technology. We compared the fitness variation of transplants among three transplant conditions.

RESULTS

A large number of significantly differentially expressed genes (DEGs) between below ambient versus ambient (109) and above ambient versus ambient (85) were identified. Transcripts involved in plant growth and development were mostly up-regulated in below ambient conditions. Transcripts involved in signalling, defence, and membrane transport were mostly up-regulated in above ambient condition. Pathway analysis revealed that most of the genes involved in metabolic processes, secondary metabolism, and flavonoid biosynthesis were differentially expressed in below ambient conditions, whereas most of the genes involved in photosynthesis and plant hormone signalling were differentially expressed in above ambient conditions. In addition, we observed higher reproductive fitness in transplant individuals at below ambient condition compared to above ambient conditions; contrary to what we expect from the cold adaptive P. sikkimensis plants.

CONCLUSIONS

We reveal P. sikkimensis's capacity for rapid adaptation to climate change through transcriptome variation, which may facilitate the phenotypic plasticity observed in morphological and life history traits. The genes and pathways identified provide a genetic resource for understanding the temperature stress (both the hot and cold stress) tolerance mechanism of P. sikkimensis in their natural environment.

Comparative transcriptome analysis of pigeonpea, Cajanus cajan (L.) and one of its wild relatives Cajanus platycarpus (Benth.) Maesen

Pigeonpea is a major source of dietary protein to the vegetarian population of the Indian sub-continent. Crop improvement to mitigate biotic and abiotic stresses for realization of its potential yield and bridging yield gap is the need of the hour. Availability of limited genomic resources in the cultivated germplasm, however, is a serious bottleneck towards successful molecular breeding for the development of superior genotypes in pigeonpea. In view of this, improvement of pigeonpea can be attempted through transgenesis or by exploiting genetic resources from its wild relatives. Pigeonpea wild relatives are known to be bestowed with agronomic traits of importance; discovery and deployment of genes from them can provide a lucrative option for crop improvement. Understanding molecular signatures of wild relatives would not only provide information about the mechanism behind desired traits but also enable us to extrapolate the information to cultivated pigeonpea. The present study deals with the characterization of leaf transcriptomes of Cajanus cajan and one of its wild relatives, Cajanus platycarpus. Illumina sequencing revealed 0.11 million transcripts in both the species with an annotation of 0.09 million (82%) transcripts using BLASTX. Comparative transcriptome analyses on the whole, divulged cues about the wild relative being vigilant and agile. Gene ontology and Mapman analysis depicted higher number of transcripts in the wild relative pertaining to signaling, transcription factors and stress responsive genes. Further, networking between the differentially expressed MapMan bins demonstrated conspicuous interactions between different bins through 535 nodes (512 Genes and 23 Pathways) and 1857 edges. The authenticity of RNA-seq analysis was confirmed by qRT-PCR. The information emanating from this study can provide valuable information and resource for future translational research including genome editing to alleviate varied stresses. Further, this learning can be a platform for in-depth investigations to decipher molecular mechanisms for mitigation of various stresses in the wild relative.

Molecular evidence for recent divergence of X- and Y-linked gene pairs in Spinacia oleracea L

Dioecy has evolved recently and independently from cosexual populations in many angiosperm lineages, providing opportunities to understand the evolutionary process underlying this transition. Spinach (Spinacia oleracea) is a dioecious plant with homomorphic sex chromosomes (XY). Although most of the spinach Y chromosome recombines with the X chromosome, a region around the male-determining locus on Y does not recombine with its X counterpart, suggesting that this region might be related to the evolution of dioecy in the species. To identify genes located in the non-recombining region (MSY, male-specific region of Y), RNA-seq analysis of male and female progeny plants (eight each) from a sib-cross of a dioecious line was performed. We discovered only 354 sex-chromosomal SNPs in 219 transcript sequences (genes). We randomly selected 39 sex-chromosomal genes to examine the reproducibility of the RNA-seq results and observed tight linkage to the male-determining locus in a spinach segregating population (140 individuals). Further analysis using a large-scale population (>1400) and over 100 spinach germplasm accessions and cultivars showed that SNPs in at least 12 genes are fully linked to the male-determining locus, suggesting that the genes reside in the spinach MSY. Synonymous substitution rates of the MSY genes and X homologues predict a recent divergence (0.40 ± 0.08 Mya). Furthermore, synonymous divergence between spinach and its wild relative (S. tetrandra), whose sex chromosomes (XY) originated from a common ancestral chromosome, predicted that the species diverged around 5.7 Mya. Assuming that dioecy in Spinacia evolved before speciation within the genus and has a monophyletic origin, our data suggest that recombination around the spinach sex-determining locus might have stopped significantly later than the evolution of dioecy in Spinacia.

Genome-Wide Transcriptomic Analysis Reveals Insights into the Response to Citrus bark cracking viroid (CBCVd) in Hop (Humulus lupulus L.)

Viroids are smallest known pathogen that consist of non-capsidated, single-stranded non-coding RNA replicons and they exploits host factors for their replication and propagation. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) is a serious threat, which spreads rapidly within hop gardens. In this study, we employed comprehensive transcriptome analyses to dissect host-viroid interactions and identify gene expression changes that are associated with disease development in hop. Our analysis revealed that CBCVd-infection resulted in the massive modulation of activity of over 2000 genes. Expression of genes associated with plant immune responses (protein kinase and mitogen-activated protein kinase), hypersensitive responses, phytohormone signaling pathways, photosynthesis, pigment metabolism, protein metabolism, sugar metabolism, and modification, and others were altered, which could be attributed to systemic symptom development upon CBCVd-infection in hop. In addition, genes encoding RNA-dependent RNA polymerase, pathogenesis-related protein, chitinase, as well as those related to basal defense responses were up-regulated. The expression levels of several genes identified from RNA sequencing analysis were confirmed by qRT-PCR. Our systematic comprehensive CBCVd-responsive transcriptome analysis provides a better understanding and insights into complex viroid-hop plant interaction. This information will assist further in the development of future measures for the prevention of CBCVd spread in hop fields.

Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors

Background

The hop plant (Humulus lupulus L.) is a valuable source of several secondary metabolites, such as flavonoids, bitter acids, and essential oils. These compounds are widely implicated in the beer brewing industry and are having potential biomedical applications. Several independent breeding programs around the world have been initiated to develop new cultivars with enriched lupulin and secondary metabolite contents but met with limited success due to several constraints. In the present work, a pioneering attempt has been made to overexpress master regulator binary transcription factor complex formed by HlWRKY1 and HlWDR1 using a plant expression vector to enhance the level of prenylflavonoid and bitter acid content in the hop. Subsequently, we performed transcriptional profiling using high-throughput RNA-Seq technology in leaves of resultant transformants and wild-type hop to gain in-depth information about the genome-wide functional changes induced by HlWRKY1 and HlWDR1 overexpression.

Results

The transgenic WW-lines exhibited an elevated expression of structural and regulatory genes involved in prenylflavonoid and bitter acid biosynthesis pathways. In addition, the comparative transcriptome analysis revealed a total of 522 transcripts involved in 30 pathways, including lipids and amino acids biosynthesis, primary carbon metabolism, phytohormone signaling and stress responses were differentially expressed in WW-transformants. It was apparent from the whole transcriptome sequencing that modulation of primary carbon metabolism and other pathways by HlWRKY1 and HlWDR1 overexpression resulted in enhanced substrate flux towards secondary metabolites pathway. The detailed analyses suggested that none of the pathways or genes, which have a detrimental effect on physiology, growth and development processes, were induced on a genome-wide scale in WW-transgenic lines.

Conclusions

Taken together, our results suggest that HlWRKY1 and HlWDR1 simultaneous overexpression positively regulates the prenylflavonoid and bitter acid biosynthesis pathways in the hop and thus these transgenes are presented as prospective candidates for achieving enhanced secondary metabolite content in the hop.

Electronic supplementary material

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

Detection of betacyanin in red-tube spinach (Spinacia oleracea) and its biofortification by strategic hydroponics

Betacyanins have been reported as water-soluble, nitrogenous pigments found in the order Caryophyllales, and they are known for powerful natural antioxidant. The biofortification of secondary metabolites such as anthocyanins and betacyanins has recently been performed in food crops by metabolic engineering through genetic modification. However, the distribution of genetically modified foods is strictly regulated. Therefore, we aimed to develop a new method for biofortifying natural antioxidants, betacyanins, without genetic modification. We first detected the presence of betacyanins in red-tube spinach (Spinacia oleracea) through ultraviolet-visible spectroscopy and mass spectrometry. We then hydroponically cultivated plants in the presence of three candidate compounds for betacyanin biofortification: dopamine, Ca2+, and sucrose. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) and antioxidant activity analyses showed that sucrose was most successful in biofortifying betacyanins, and reverse transcription polymerase chain reaction (RT-PCR) indicated that several genes involved in betacyanin biosynthesis were induced by sucrose. Therefore, strategic hydroponics represents a new approach for cultivating betacyanin-enriched vegetables.

Expression Analysis of Oxalate Metabolic Pathway Genes Reveals Oxalate Regulation Patterns in Spinach

Spinach (Spinacia oleracea L.) is one of most important leafy vegetables because of its high nutritional value and high oxalate content, which can be toxic with negative effects on human nutrition. Ammonium and nitrate can effectively regulate oxalate accumulation, although the mechanisms underlying the oxalate biosynthesis and regulation are still undetermined in plants. In the present study, we identified 25 putative genes that are involved in the oxalate biosynthetic and degradation pathway, before analyzing the oxalate content and the expression levels of the corresponding proteins under normal growth conditions, with or without ammonium and nitrate treatments, using high and low oxalate-accumulated spinach genotypes. The two cultivars exhibited different profiles of total oxalate and soluble oxalate accumulation. The high oxalate concentrations in spinach were as a result of the high transcription levels of the genes that are involved in oxalate biosynthesis under normal growth conditions, such as SoGLO2, SoGLO3, three SoOXACs, SoMLS, SoMDH1, SoMDH2, and SoMDH4. The results revealed that the ammonium and nitrate were able to control the oxalate content in leaves, possibly because of the different transcription levels of the genes. The oxalate content is regulated by complex regulatory mechanisms and is varied in the different varieties of spinach. The results from this research may be used to assist the investigation of the mechanism of oxalate regulation and breeding for reduced oxalate content in spinach.

The Huperzia selago Shoot Tip Transcriptome Sheds New Light on the Evolution of Leaves

Lycopodiophyta—consisting of three orders, Lycopodiales, Isoetales and Selaginellales, with different types of shoot apical meristems (SAMs)—form the earliest branch among the extant vascular plants. They represent a sister group to all other vascular plants, from which they differ in that their leaves are microphylls—that is, leaves with a single, unbranched vein, emerging from the protostele without a leaf gap—not megaphylls. All leaves represent determinate organs originating on the flanks of indeterminate SAMs. Thus, leaf formation requires the suppression of indeterminacy, that is, of KNOX transcription factors. In seed plants, this is mediated by different groups of transcription factors including ARP and YABBY.

We generated a shoot tip transcriptome of Huperzia selago (Lycopodiales) to examine the genes involved in leaf formation. Our H. selago transcriptome does not contain any ARP homolog, although transcriptomes of Selaginella spp. do. Surprisingly, we discovered a YABBY homolog, although these transcription factors were assumed to have evolved only in seed plants.

The existence of a YABBY homolog in H. selago suggests that YABBY evolved already in the common ancestor of the vascular plants, and subsequently was lost in some lineages like Selaginellales, whereas ARP may have been lost in Lycopodiales. The presence of YABBY in the common ancestor of vascular plants would also support the hypothesis that this common ancestor had a simplex SAM. Furthermore, a comparison of the expression patterns of ARP in shoot tips of Selaginella kraussiana (Harrison CJ, etal. 2005. Independent recruitment of a conserved developmental mechanism during leaf evolution. Nature 434(7032):509–514.) and YABBY in shoot tips of H. selago implies that the development of microphylls, unlike megaphylls, does not seem to depend on the combined activities of ARP and YABBY. Altogether, our data show that Lycopodiophyta are a diverse group; so, in order to understand the role of Lycopodiophyta in evolution, representatives of Lycopodiales, Selaginellales, as well as of Isoetales, have to be examined.

Higher anthocyanin accumulation associated with higher transcription levels of anthocyanin biosynthesis genes in spinach

Spinach (Spinacia oleracea L.) is widely cultivated as an economically important green leafy vegetable crop for fresh and processing consumption. The red-purple spinach shows abundant anthocyanin accumulation in the leaf and leaf petiole. However, the molecular mechanisms of anthocyanin synthesis in this species are still undetermined. In the present study, we investigated pigment formation and identified anthocyanin biosynthetic genes in spinach. We also analyzed the expression of these genes in purple and green cultivars by quantitative PCR. The accumulation of anthocyanin showed that it was the dominant pigment resulting in the red coloration in spinach. In total, 22 biosynthesis genes and 25 regulatory genes were identified in spinach, based on the spinach genomic and transcriptomic database. Furthermore, the expression patterns of genes encoding enzymes indicated that SoPAL, SoUFGT3, and SoUFGT4 were possible candidate genes for anthocyanin biosynthesis in red-purple spinach. The expression patterns of transcription factors indicated that two SoMYB genes, three SobHLH genes, and one SoWD40 gene were drastically up-regulated and co-expression in red-purple spinach, suggesting an essential role of regulatory genes in the anthocyanin biosynthesis of spinach. These results will enhance our understanding of the molecular mechanisms of anthocyanin biosynthesis in purple spinach.

Proteomics and Phosphoproteomics of Heat Stress-Responsive Mechanisms in Spinach

Elevated temperatures limit plant growth and reproduction and pose a growing threat to agriculture. Plant heat stress response is highly conserved and fine-tuned in multiple pathways. Spinach (Spinacia oleracea L.) is a cold tolerant but heat sensitive green leafy vegetable. In this study, heat adaptation mechanisms in a spinach sibling inbred heat-tolerant line Sp75 were investigated using physiological, proteomic, and phosphoproteomic approaches. The abundance patterns of 911 heat stress-responsive proteins, and phosphorylation level changes of 45 phosphoproteins indicated heat-induced calcium-mediated signaling, ROS homeostasis, endomembrane trafficking, and cross-membrane transport pathways, as well as more than 15 transcription regulation factors. Although photosynthesis was inhibited, diverse primary and secondary metabolic pathways were employed for defense against heat stress, such as glycolysis, pentose phosphate pathway, amino acid metabolism, fatty acid metabolism, nucleotide metabolism, vitamin metabolism, and isoprenoid biosynthesis. These data constitute a heat stress-responsive metabolic atlas in spinach, which will springboard further investigations into the sophisticated molecular mechanisms of plant heat adaptation and inform spinach molecular breeding initiatives.

Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions

Spinach is an important leafy vegetable enriched with multiple necessary nutrients. Here we report the draft genome sequence of spinach (Spinacia oleracea, 2n=12), which contains 25,495 protein-coding genes. The spinach genome is highly repetitive with 74.4% of its content in the form of transposable elements. No recent whole genome duplication events are observed in spinach. Genome syntenic analysis between spinach and sugar beet suggests substantial inter- and intra-chromosome rearrangements during the Caryophyllales genome evolution. Transcriptome sequencing of 120 cultivated and wild spinach accessions reveals more than 420 K variants. Our data suggests that S. turkestanica is likely the direct progenitor of cultivated spinach and spinach domestication has a weak bottleneck. We identify 93 domestication sweeps in the spinach genome, some of which are associated with important agronomic traits including bolting, flowering and leaf numbers. This study offers insights into spinach evolution and domestication and provides resources for spinach research and improvement.

Spinach is an economically important vegetable crop but previous genomic resources were of limited use for comparative and functional analyses. Here, Xu et al. present a high quality draft spinach genome and transcriptome data for multiple Spinacia accessions providing insight into Caryophyllales genome evolution.

Towards a whole-genome sequence for rye (Secale cereale L.)

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Transcriptome Analysis and Differential Gene Expression on the Testis of Orange Mud Crab, Scylla olivacea, during Sexual Maturation

Adequate genetic information is essential for sustainable crustacean fisheries and aquaculture management. The commercially important orange mud crab, Scylla olivacea, is prevalent in Southeast Asia region and is highly sought after. Although it is a suitable aquaculture candidate, full domestication of this species is hampered by the lack of knowledge about the sexual maturation process and the molecular mechanisms behind it, especially in males. To date, data on its whole genome is yet to be reported for S. olivacea. The available transcriptome data published previously on this species focus primarily on females and the role of central nervous system in reproductive development. De novo transcriptome sequencing for the testes of S. olivacea from immature, maturing and mature stages were performed. A total of approximately 144 million high-quality reads were generated and de novo assembled into 160,569 transcripts with a total length of 142.2 Mb. Approximately 15–23% of the total assembled transcripts were annotated when compared to public protein sequence databases (i.e. UniProt database, Interpro database, Pfam database and Drosophila melanogaster protein database), and GO-categorised with GO Ontology terms. A total of 156,181 high-quality Single-Nucleotide Polymorphisms (SNPs) were mined from the transcriptome data of present study. Transcriptome comparison among the testes of different maturation stages revealed one gene (beta crystallin like gene) with the most significant differential expression—up-regulated in immature stage and down-regulated in maturing and mature stages. This was further validated by qRT-PCR. In conclusion, a comprehensive transcriptome of the testis of orange mud crabs from different maturation stages were obtained. This report provides an invaluable resource for enhancing our understanding of this species’ genome structure and biology, as expressed and controlled by their gonads.

De novo assembly, functional annotation, and analysis of the giant reed (Arundo donax L.) leaf transcriptome provide tools for the development of a biofuel feedstock

BACKGROUND

Arundo donax has attracted renewed interest as a potential candidate energy crop for use in biomass-to-liquid fuel conversion processes and biorefineries. This is due to its high productivity, adaptability to marginal land conditions, and suitability for biofuel and biomaterial production. Despite its importance, the genomic resources currently available for supporting the improvement of this species are still limited.

RESULTS

We used RNA sequencing (RNA-Seq) to de novo assemble and characterize the A. donax leaf transcriptome. The sequencing generated 1249 million clean reads that were assembled using single-k-mer and multi-k-mer approaches into 62,596 unique sequences (unitranscripts) with an N50 of 1134 bp. TransDecoder and Trinotate software suites were used to obtain putative coding sequences and annotate them by mapping to UniProtKB/Swiss-Prot and UniRef90 databases, searching for known transcripts, proteins, protein domains, and signal peptides. Furthermore, the unitranscripts were annotated by mapping them to the NCBI non-redundant, GO and KEGG pathway databases using Blast2GO. The transcriptome was also characterized by BLAST searches to investigate homologous transcripts of key genes involved in important metabolic pathways, such as lignin, cellulose, purine, and thiamine biosynthesis and carbon fixation. Moreover, a set of homologous transcripts of key genes involved in stomatal development and of genes coding for stress-associated proteins (SAPs) were identified. Additionally, 8364 simple sequence repeat (SSR) markers were identified and surveyed. SSRs appeared more abundant in non-coding regions (63.18%) than in coding regions (36.82%). This SSR dataset represents the first marker catalogue of A. donax. 53 SSRs (PolySSRs) were then predicted to be polymorphic between ecotype-specific assemblies, suggesting genetic variability in the studied ecotypes.

CONCLUSIONS

This study provides the first publicly available leaf transcriptome for the A. donax bioenergy crop. The functional annotation and characterization of the transcriptome will be highly useful for providing insight into the molecular mechanisms underlying its extreme adaptability. The identification of homologous transcripts involved in key metabolic pathways offers a platform for directing future efforts in genetic improvement of this species. Finally, the identified SSRs will facilitate the harnessing of untapped genetic diversity. This transcriptome should be of value to ongoing functional genomics and genetic studies in this crop of paramount economic importance.

Characterization of the transcriptomes and cuticular protein gene expression of alate adult, brachypterous neotenic and adultoid reproductives of Reticulitermes labralis

The separation of primary reproductive and secondary reproductive roles based on the differentiation of alate adults and neotenic reproductives is the most prominent characteristic of termites. To clarify the mechanism underlying this differentiation, we sequenced the transcriptomes of alate adults (ARs), brachypterous neotenics (BNs) and adultoid reproductives (ANs) from the last instar nymphs of Reticulitermes labralis. A total of 404,152,188 clean sequencing reads was obtained and 61,953 unigenes were assembled. Of the 54 identified cuticular protein (CP) genes of the reproductives, 22 were classified into the CPR family and 7 were classified into the CPG family. qRT-PCR analyses of the 6 CP genes revealed that the CP genes involved in exocuticle sclerotization were highly expressed in the ARs and RR-1 involved in soft endocuticle was highly expressed in the ARs and ANs. These results suggest that the alate adults might increase cuticular component deposition to adapt to new or changing environments and that the development of reproductive individuals into primary or secondary reproductives is controlled by the expression of cuticular protein genes involved in the hardening of the exocuticle. In addition, the AN caste is a transitional type between the BN and AR castes in the process of evolution.

Unravelling molecular mechanisms from floral initiation to lipid biosynthesis in a promising biofuel tree species, Pongamia pinnata using transcriptome analysis

Pongamia pinnata (L.) (Fabaceae) is a promising biofuel tree species which is underexploited in the areas of both fundamental and applied research, due to the lack of information either on transcriptome or genomic data. To investigate the possible metabolic pathways, we performed whole transcriptome analysis of Pongamia through Illumina NextSeq platform and generated 2.8 GB of paired end sequence reads. The de novo assembly of raw reads generated 40,000 contigs and 35,000 transcripts, representing leaf, flower and seed unigenes. Spatial and temporal expression profiles of photoperiod and floral homeotic genes in Pongamia, identified GIGANTEA (GI) - CONSTANS (CO) - FLOWERING LOCUS T (FT) as active signal cascade for floral initiation. Four prominent stages of seed development were selected in a high yielding Pongamia accession (TOIL 1) to follow the temporal expression patterns of key fatty acid biosynthetic genes involved in lipid biosynthesis and accumulation. Our results provide insights into an array of molecular events from flowering to seed maturity in Pongamia which will provide substantial basis for modulation of fatty acid composition and enhancing oil yields which should serve as a potential feedstock for biofuel production.

Genetic diversity and association analysis of leafminer (Liriomyza langei) resistance in spinach (Spinacia oleracea)

Leafminer (Liriomyza langei) is a major insect pest of many important agricultural crops, including spinach (Spinacia oleracea). Use of genetic resistance is an efficient, economic, and environment-friendly method to control this pest. The objective of this research was to conduct association analysis and identify single nucleotide polymorphism (SNP) markers associated with leafminer resistance in spinach germplasm. A total of 300 USDA spinach germplasm accessions were used for the association analysis of leafminer resistance. Genotyping by sequencing (GBS) was used for genotyping and 783 SNPs from GBS were used for association analysis. The leafminer resistance showed a near normal distribution with a wide range from 1.1 to 11.7 stings per square centimeter leaf area, suggesting that the leafminer resistance in spinach is a complex trait controlled by multiple genes with minor effect in this spinach panel. Association analysis indicated that five SNP markers, AYZV02040968_7171, AYZV02076752_412, AYZV02098618_4615, AYZV02147304_383, and AYZV02271373_398, were associated with the leafminer resistance with LOD 2.5 or higher. The SNP markers may be useful for breeders to select plants and lines for leafminer resistance in spinach breeding programs through marker-assisted selection.

SpinachDB: A Well-Characterized Genomic Database for Gene Family Classification and SNP Information of Spinach

Spinach (Spinacia oleracea L.), which originated in central and western Asia, belongs to the family Amaranthaceae. Spinach is one of most important leafy vegetables with a high nutritional value as well as being a perfect research material for plant sex chromosome models. As the completion of genome assembly and gene prediction of spinach, we developed SpinachDB (http://222.73.98.124/spinachdb) to store, annotate, mine and analyze genomics and genetics datasets efficiently. In this study, all of 21702 spinach genes were annotated. A total of 15741 spinach genes were catalogued into 4351 families, including identification of a substantial number of transcription factors. To construct a high-density genetic map, a total of 131592 SSRs and 1125743 potential SNPs located in 548801 loci of spinach genome were identified in 11 cultivated and wild spinach cultivars. The expression profiles were also performed with RNA-seq data using the FPKM method, which could be used to compare the genes. Paralogs in spinach and the orthologous genes in Arabidopsis, grape, sugar beet and rice were identified for comparative genome analysis. Finally, the SpinachDB website contains seven main sections, including the homepage; the GBrowse map that integrates genome, genes, SSR and SNP marker information; the Blast alignment service; the gene family classification search tool; the orthologous and paralogous gene pairs search tool; and the download and useful contact information. SpinachDB will be continually expanded to include newly generated robust genomics and genetics data sets along with the associated data mining and analysis tools.