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

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

De Novo Transcriptome Assembly, Gene Annotations, and Characterization of Functional Profiling Reveal Key Genes for Lead Alleviation in the Pb Hyperaccumulator Greek Mustard (Hirschfeldia incana L.)

Lead (Pb) contamination is a widespread environmental problem due to its toxicity to living organisms. Hirschfeldia incana L., a member of the Brassicaceae family, commonly found in the Mediterranean regions, is characterized by its ability to tolerate and accumulate Pb in soils and hydroponic cultures. This plant has been reported as an excellent model to assess the response of plants to Pb. However, the lack of genomic data for H. incana hinders research at the molecular level. In the present study, we carried out RNA deep transcriptome sequencing (RNA-seq) of H. incana under two conditions, control without Pb(NO₃)₂ and treatment with 100 µM of Pb(NO₃)₂ for 15 days. A total of 797.83 million reads were generated using Illumina sequencing technology. We assembled 77,491 transcript sequences with an average length of 959 bp and N50 of 1330 bp. Sequence similarity analyses and annotation of these transcripts were performed against the Arabidopsis thaliana nr protein database, Gene Ontology (GO), and KEGG databases. As a result, 13,046 GO terms and 138 KEGG maps were created. Under Pb stress, 577 and 270 genes were differentially expressed in roots and aboveground parts, respectively. Detailed elucidation of regulation of metal transporters, transcription factors (TFs), and plant hormone genes described the role of actors that allow the plant to fine-tune Pb stress responses. Our study revealed that several genes related to jasmonic acid biosynthesis and alpha-linoleic acid were upregulated, suggesting these components’ implication in Hirschfeldia incana L responses to Pb stress. This study provides data for further genomic analyses of the biological and molecular mechanisms leading to Pb tolerance and accumulation in Hirschfeldia incana L.

Global leaf and root transcriptome in response to cadmium reveals tolerance mechanisms in Arundo donax L

The expected increase of sustainable energy demand has shifted the attention towards bioenergy crops. Due to their know tolerance against abiotic stress and relatively low nutritional requirements, they have been proposed as election crops to be cultivated in marginal lands without disturbing the part of lands employed for agricultural purposes. Arundo donax L. is a promising bioenergy crop whose behaviour under water and salt stress has been recently studied at transcriptomic levels. As the anthropogenic activities produced in the last years a worrying increase of cadmium contamination worldwide, the aim of our work was to decipher the global transcriptomic response of A. donax leaf and root in the perspective of its cultivation in contaminated soil. In our study, RNA-seq libraries yielded a total of 416 million clean reads and 10.4 Gb per sample. De novo assembly of clean reads resulted in 378,521 transcripts and 126,668 unigenes with N50 length of 1812 bp and 1555 bp, respectively. Differential gene expression analysis revealed 5,303 deregulated transcripts (3,206 up- and 2,097 down regulated) specifically observed in the Cd-treated roots compared to Cd-treated leaves. Among them, we identified genes related to “Protein biosynthesis”, “Phytohormone action”, “Nutrient uptake”, “Cell wall organisation”, “Polyamine metabolism”, “Reactive oxygen species metabolism” and “Ion membrane transport”. Globally, our results indicate that ethylene biosynthesis and the downstream signal cascade are strongly induced by cadmium stress. In accordance to ethylene role in the interaction with the ROS generation and scavenging machinery, the transcription of several genes (NADPH oxidase 1, superoxide dismutase, ascorbate peroxidase, different glutathione S-transferases and catalase) devoted to cope the oxidative stress is strongly activated. Several small signal peptides belonging to ROTUNDIFOLIA, CLAVATA3, and C-TERMINALLY ENCODED PEPTIDE 1 (CEP) are also among the up-regulated genes in Cd-treated roots functioning as messenger molecules from root to shoot in order to communicate the stressful status to the upper part of the plants. Finally, the main finding of our work is that genes involved in cell wall remodelling and lignification are decisively up-regulated in giant reed roots. This probably represents a mechanism to avoid cadmium uptake which strongly supports the possibility to cultivate giant cane in contaminated soils in the perspective to reserve agricultural soil for food and feed crops.

Identification of Known and Novel Arundo donax L. MicroRNAs and Their Targets Using High-Throughput Sequencing and Degradome Analysis

MicroRNAs (miRNAs) are a class of non-coding molecules involved in the regulation of a variety of biological processes. They have been identified and characterized in several plant species, but only limited data are available for Arundo donax L., one of the most promising bioenergy crops. Here we identified, for the first time, A. donax conserved and novel miRNAs together with their targets, through a combined analysis of high-throughput sequencing of small RNAs, transcriptome and degradome data. A total of 134 conserved miRNAs, belonging to 45 families, and 27 novel miRNA candidates were identified, along with the corresponding primary and precursor miRNA sequences. A total of 96 targets, 69 for known miRNAs and 27 for novel miRNA candidates, were also identified by degradome analysis and selected slice sites were validated by 5′-RACE. The identified set of conserved and novel candidate miRNAs, together with their targets, extends our knowledge about miRNAs in monocots and pave the way to further investigations on miRNAs-mediated regulatory processes in A. donax, Poaceae and other bioenergy crops.

Improving the quality of barley transcriptome de novo assembling by using a hybrid approach for lines with varying spike and stem coloration

De novo transcriptome assembly is an important stage of RNA-seq data computational analysis. It allows the researchers to obtain the sequences of transcripts presented in the biological sample of interest. The availability of accurate and complete transcriptome sequence of the organism of interest is, in turn, an indispensable condition for further analysis of RNA-seq data. Through years of transcriptomic research, the bioinformatics community has developed a number of assembler programs for transcriptome reconstruction from short reads of RNA-seq libraries. Different assemblers makes it possible to conduct a de novo transcriptome reconstruction and a genome-guided reconstruction. The majority of the assemblers working with RNA-seq data are based on the De Bruijn graph method of sequence reconstruction. However, specif ics of their procedures can vary drastically, as do their results. A number of authors recommend a hybrid approach to transcriptome reconstruction based on combining the results of several assemblers in order to achieve a better transcriptome assembly. The advantage of this approach has been demonstrated in a number of studies, with RNA-seq experiments conducted on the Illumina platform. In this paper, we propose a hybrid approach for creating a transcriptome assembly of the barley Hordeum vulgare isogenic line Bowman and two nearly isogenic lines contrasting in spike pigmentation, based on the results of sequencing on the IonTorrent platform. This approach implements several de novo assemblers: Trinity, Trans-ABySS and rnaSPAdes. Several assembly metrics were examined: the percentage of reference transcripts observed in the assemblies, the percentage of RNA-seq reads involved, and BUSCO scores. It was shown that, based on the summation of these metrics, transcriptome meta-assembly surpasses individual transcriptome assemblies it consists of.

The First De Novo Transcriptome Assembly and Transcriptomic Dynamics of the Mangrove Tree Rhizophora stylosa Griff. (Rhizophoraceae)

Mangroves are salt-tolerant plant species that grow in coastal saline water and are adapted to harsh environmental conditions. In this study, we de novo assembled and functionally annotated the transcriptome of Rhizophora stylosa, the widely distributed mangrove from the largest mangrove family (Rhizophoraceae). The final transcriptome consists of 200,491 unigenes with an average length, and N50 of 912.7 and 1334 base pair, respectively. We then compared the genome-wide expression profiles between the two morphologically distinct natural populations of this species growing under different levels of salinity depending on their distance from the ocean. Among the 200,491 unigenes, 40,253 were identified as differentially expressed between the two populations, while 15,741 and 24,512 were up- and down-regulated, respectively. Functional annotation assigned thousands of upregulated genes in saline environment to the categories related to abiotic stresses such as response to salt-, osmotic-, and oxidative-stress. Validation of those genes may contribute to a better understanding of adaptation in mangroves species. This study reported, for the first time, the transcriptome of R. stylosa, and the dynamic of it in response to salt stress and provided a valuable resource for elucidation of the molecular mechanism underlying the salt stress response in mangroves and other plants that live under stress.

Unveiling the transcriptomic complexity of Miscanthus sinensis using a combination of PacBio long read- and Illumina short read sequencing platforms

Background

Miscanthus sinensis Andersson is a perennial grass that exhibits remarkable lignocellulose characteristics suitable for sustainable bioenergy production. However, knowledge of the genetic resources of this species is relatively limited, which considerably hampers further work on its biology and genetic improvement.

Results

In this study, through analyzing the transcriptome of mixed samples of leaves and stems using the latest PacBio Iso-Seq sequencing technology combined with Illumina HiSeq, we report the first full-length transcriptome dataset of M. sinensis with a total of 58.21 Gb clean data. An average of 15.75 Gb clean reads of each sample were obtained from the PacBio Iso-Seq system, which doubled the data size (6.68 Gb) obtained from the Illumina HiSeq platform. The integrated analyses of PacBio- and Illumina-based transcriptomic data uncovered 408,801 non-redundant transcripts with an average length of 1,685 bp. Of those, 189,406 transcripts were commonly identified by both methods, 169,149 transcripts with an average length of 619 bp were uniquely identified by Illumina HiSeq, and 51,246 transcripts with an average length of 2,535 bp were uniquely identified by PacBio Iso-Seq. Approximately 96 % of the final combined transcripts were mapped back to the Miscanthus genome, reflecting the high quality and coverage of our sequencing results. When comparing our data with genomes of four species of Andropogoneae, M. sinensis showed the closest relationship with sugarcane with up to 93 % mapping ratios, followed by sorghum with up to 80 % mapping ratios, indicating a high conservation of orthologs in these three genomes. Furthermore, 306,228 transcripts were successfully annotated against public databases including cell wall related genes and transcript factor families, thus providing many new insights into gene functions. The PacBio Iso-Seq data also helped identify 3,898 alternative splicing events and 2,963 annotated AS isoforms within 10 function categories.

Conclusions

Taken together, the present study provides a rich data set of full-length transcripts that greatly enriches our understanding of M. sinensis transcriptomic resources, thus facilitating further genetic improvement and molecular studies of the Miscanthus species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07971-x.

De Novo Transcriptome Assembly, Functional Annotation, and Transcriptome Dynamics Analyses Reveal Stress Tolerance Genes in Mangrove Tree (Bruguiera gymnorhiza)

Their high adaptability to difficult coastal conditions makes mangrove trees a valuable resource and an interesting model system for understanding the molecular mechanisms underlying stress tolerance and adaptation of plants to the stressful environmental conditions. In this study, we used RNA sequencing (RNA-Seq) for de novo assembling and characterizing the Bruguiera gymnorhiza (L.) Lamk leaf transcriptome. B. gymnorhiza is one of the most widely distributed mangrove species from the biggest family of mangroves; Rhizophoraceae. The de novo assembly was followed by functional annotations and identification of individual transcripts and gene families that are involved in abiotic stress response. We then compared the genome-wide expression profiles between two populations of B. gymnorhiza, growing under different levels of stress, in their natural habitats. One population living in high salinity environment, in the shore of the Pacific Ocean- Japan, and the other population living about one kilometre farther from the ocean, and next to the estuary of a river; in less saline and more brackish condition. Many genes involved in response to salt and osmotic stress, showed elevated expression levels in trees growing next to the ocean in high salinity condition. Validation of these genes may contribute to future salt-resistance research in mangroves and other woody plants. Furthermore, the sequences and transcriptome data provided in this study are valuable scientific resources for future comparative transcriptome research in plants growing under stressful conditions.

Transcriptome Analysis of Tetraploid and Octoploid Common Reed (Phragmites australis)

Polyploidization in plants is thought to have occurred as coping mechanism with environmental stresses. Polyploidization-driven adaptation is often achieved through interplay of gene networks involved in differentially expressed genes, which triggers the plant to evolve special phenotypic traits for survival. Phragmites australis is a cosmopolitan species with highly variable phenotypic traits and high adaptation capacity to various habitats. The species' ploidy level varies from 3x to 12x, thus it is an ideal organism to investigate the molecular evolution of polyploidy and gene regulation mediated by different numbers of chromosome copies. In this study, we used high-throughput RNAseq data as a tool, to analyze the gene expression profiles in tetraploid and octoploid P. australis. The estimated divergence time between tetraploid and octoploid P. australis was dated to the border between Pliocene and Pleistocene. This study identified 439 up- and 956 down-regulated transcripts in tetraploids compared to octoploids. Gene ontology and pathway analysis revealed that tetraploids tended to express genes responsible for reproduction and seed germination to complete the reproduction cycle early, and expressed genes related to defense against UV-B light and fungi, whereas octoploids expressed mainly genes related to thermotolerance. Most differentially expressed genes were enriched in chaperones, folding catalysts and protein processing in endoplasmic reticulum pathways. Multiple biased isoform usage of the same gene was detected in differentially expressed genes, and the ones upregulated in octoploids were related to reduced DNA methylation. Our study provides new insights into the role of polyploidization on environmental responses and potential stress tolerance in grass species.

Interspecific and intraspecific phenotypic diversity for drought adaptation in bioenergy Arundo species

Biomass crops are commonly grown in low-grade land and selection of drought-tolerant accessions is of major importance to sustain productivity. In this work, we assess phenotypic variation under different environmental scenarios in a series of accessions of Arundo donax, and contrast it with two closely related species, Arundo donaciformis and Arundo plinii. Gas-exchange and stomatal anatomy analysis showed an elevated photosynthetic capacity in A. plinii compared to A. donax and A. donaciformis with a significant intraspecific variation in A. donax. The three species showed significantly contrasting behaviour of transpiration under developing water stress and increasing vapour pressure deficit (VPD), with A. donax being the most conservative while A. plinii showed an elevated degree of insensitivity to environmental cues. Under optimal conditions, A. donax had the highest estimated leaf area (projected leaf area) and plant dry weight although a significant reduction under water stress was observed for A. donax and A. donaciformis accessions while no differences were recorded for A. plinii between optimal growing conditions (well-watered [WW]) and reduced soil water availability (water-stressed [WS]). A. donax displayed a markedly conservative water use behaviour but elevated sensitivity of biomass accumulation under stress conditions. By contrast, in A. plinii, biomass and transpiration were largely insensitive to WS and increasing VPD, though biomass dry weight under optimal conditions was significantly lower than A. donax. We provide evidence of interspecific phenotypic variation within the Arundo genus while the intraspecific phenotypic plasticity may be exploited for further selection of superior clones under disadvantageous environmental conditions. The extensive trade-off between water use and biomass accumulation present in the three species under stress conditions provides a series of novel traits to be exploited in the selection of superior clones adapted to different environmental scenarios. Non-destructive approaches are provided to screen large populations for water-stress-tolerant A. donax clones.

Silencing of a LIM gene in cotton exhibits enhanced resistance against Apolygus lucorum

Plant bugs (Miridae species) have become major agricultural pests that cause increasing and severe economic damage. Plant-mediated RNA interference (RNAi) is emerging as an eco-friendly, efficient, and reliable strategy for pest management. In this study, we isolated and characterized a lethal gene of Apolygus lucorum and named it Apolygus lucorum LIM (AlLIM), which produced A. lucorum mortality rates ranging from 38% to 81%. Downregulation of the AlLIM gene expression in A. lucorum by injection of a double-stranded RNA (dsRNA) led to muscle structural disorganization that resulted in metamorphosis deficiency and increased mortality. Then we constructed a plant expression vector that enabled transgenic cotton to highly and stably express dsRNA of AlLIM (dsAlLIM) by Agrobacterium-mediated genetic transformation. In the field bioassay, dsAlLIM transgenic cotton was protected from A. lucorum damage with high efficiency, with almost no detectable yield loss. Therefore, our study successfully provides a promising genetically modified strategy to overpower A. lucorum attack.

The development of SSR markers based on RNA-sequencing and its validation between and within Carex L. species

BACKGROUND

Carex L. is one of the largest genera in the Cyperaceae family and an important vascular plant in the ecosystem. However, the genetic background of Carex is complex and the classification is not clear. In order to investigate the gene function annotation of Carex, RNA-sequencing analysis was performed. Simple sequence repeats (SSRs) were generated based on the Illumina data and then were utilized to investigate the genetic characteristics of the 79 Carex germplasms.

RESULTS

In this study, 36,403 unigenes with a total length of 41,724,615 bp were obtained and annotated based on GO, KOG, KEGG, NR databases. The results provide a theoretical basis for gene function exploration. Out of 8776 SSRs, 96 pairs of primers were randomly selected. One hundred eighty polymorphic bands were amplified with a polymorphism rate of 100% based on 42 pairs of primers with higher polymorphism levels. The average band number was 4.3 per primer, the average distance value was 0.548, and the polymorphic information content was ranged from 0.133 to 0.494. The number of observed alleles (Na), effective alleles (Ne), Nei's (1973) gene diversity (H), and the Shannon information index (I) were 2.000, 1.376, 0.243, and 0.391, respectively. NJ clustering divided into three groups and the accessions from New Zealand showed a similar genetic attribute and clustered into one group. UPGMA and PCoA analysis also revealed the same result. The analysis of molecular variance (AMOVA) revealed a superior genetic diversity within accessions than between accessions based on geographic origin cluster and NJ cluster. What's more, the fingerprints of 79 Carex species are established in this study. Different combinations of primer pairs can be used to identify multiple Carex at one time, which overcomes the difficulties of traditional identification methods.

CONCLUSIONS

The transcriptomic analysis shed new light on the function categories from the annotated genes and will facilitate future gene functional studies. The genetic characteristics analysis indicated that gene flow was extensive among 79 Carex species. These markers can be used to investigate the evolutionary history of Carex and related species, as well as to serve as a guide in future breeding projects.

Genetic Improvement of Arundo donax L.: Opportunities and Challenges

Arundo donax L., the giant reed-being a long-duration, low-cost, non-food energy crop able to grow in marginal lands-has emerged as a potential alternative to produce biomass for both energy production, with low carbon emissions, and industrial bioproducts. In recent years, pioneering efforts have been made to genetically improve this very promising energy crop. This review analyses the recent advances and challenges encountered in using clonal selection, mutagenesis/somaclonal variation and transgenesis/genome editing. Attempts to improve crop yield, in vitro propagation efficiency, salt and heavy metal tolerance by clonal selection were carried out, although limited by the species' low genetic diversity and availability of mutants. Mutagenesis and somaclonal variation have also been attempted on this species; however, since Arundo donax is polyploid, it is very difficult to induce and select promising mutations. In more recent years, genomics and transcriptomics data are becoming available in Arundo, closing the gap to make possible the genetic manipulation of this energy crop in the near future. The challenge will regard the functional characterization of the genes/sequences generated by genomic sequencing and transcriptomic analysis in a complex polyploid genome.

UAV-Based LiDAR for High-Throughput Determination of Plant Height and Above-Ground Biomass of the Bioenergy Grass Arundo donax

Replacing fossil fuels with cellulosic biofuels is a valuable component of reducing the drivers of climate change. This leads to a requirement to develop more productive bioenergy crops, such as Arundo donax with the aim of increasing above-ground biomass (AGB). However, direct measurement of AGB is time consuming, destructive, and labor-intensive. Phenotyping of plant height and biomass production is a bottleneck in genomics- and phenomics-assisted breeding. Here, an unmanned aerial vehicle (UAV) for remote sensing equipped with light detection and ranging (LiDAR) was tested for remote plant height and biomass determination in A. donax. Experiments were conducted on three A. donax ecotypes grown in well-watered and moderate drought stress conditions. A novel UAV-LiDAR data collection and processing workflow produced a dense three-dimensional (3D) point cloud for crop height estimation through a normalized digital surface model (DSM) that acts as a crop height model (CHM). Manual measurements of crop height and biomass were taken in parallel and compared to LiDAR CHM estimates. Stepwise multiple regression was used to estimate biomass. Analysis of variance (ANOVA) tests and pairwise comparisons were used to determine differences between ecotypes and drought stress treatments. We found a significant relationship between the sensor readings and manually measured crop height and biomass, with determination coefficients of 0.73 and 0.71 for height and biomass, respectively. Differences in crop heights were detected more precisely from LiDAR estimates than from manual measurement. Crop biomass differences were also more evident in LiDAR estimates, suggesting differences in ecotypes’ productivity and tolerance to drought. Based on these results, application of the presented UAV-LiDAR workflow will provide new opportunities in assessing bioenergy crop morpho-physiological traits and in delivering improved genotypes for biorefining.

Transcriptomics of the Rooibos (Aspalathus linearis) Species Complex

Rooibos (Aspalathus linearis), widely known as a herbal tea, is endemic to the Cape Floristic Region of South Africa (SA). It produces a wide range of phenolic compounds that have been associated with diverse health promoting properties of the plant. The species comprises several growth forms that differ in their morphology and biochemical composition, only one of which is cultivated and used commercially. Here, we established methodologies for non-invasive transcriptome research of wild-growing South African plant species, including (1) harvesting and transport of plant material suitable for RNA sequencing; (2) inexpensive, high-throughput biochemical sample screening; (3) extraction of high-quality RNA from recalcitrant, polysaccharide- and polyphenol rich plant material; and (4) biocomputational analysis of Illumina sequencing data, together with the evaluation of programs for transcriptome assembly (Trinity, IDBA-Trans, SOAPdenovo-Trans, CLC), protein prediction, as well as functional and taxonomic transcript annotation. In the process, we established a biochemically characterized sample pool from 44 distinct rooibos ecotypes (1–5 harvests) and generated four in-depth annotated transcriptomes (each comprising on average ≈86,000 transcripts) from rooibos plants that represent distinct growth forms and differ in their biochemical profiles. These resources will serve future rooibos research and plant breeding endeavours.

Transcriptional response of giant reed (Arundo donax L.) low ecotype to long-term salt stress by unigene-based RNAseq

The giant reed is a fast growing herbaceous non-food crop considered as eligible alternative energy source to reduce the usage of fossil fuels. Tolerance of this plant to abiotic stress has been demonstrated across a range of stressful conditions, thus allowing cultivation in marginal or poorly cultivated land in order not to compromise food security and to overcome land use controversies. In this work, we de novo sequenced, assembled and analyzed the A. donax low G34 ecotype leaf transcriptome (RNAseq analysis) subjected to severe long-term salt stress (256.67 mM NaCl corresponding to 32 dS m^-1 electric conductibility). In order to shed light upon the response to high salinity of this non model plant, we analyzed clusters related to salt sensory and signaling transduction, transcription factors, hormone regulation, Reactive Oxygen Species (ROS) scavenging and osmolyte biosynthesis, all of them showing different regulation compared to untreated plants. The analysis of clusters related to ethylene biosynthesis and signaling indicated that gene transcription is modulated towards the minimization of ethylene negative effects upon plant growth. Certainly, the photosynthesis is strongly affected since genes involved in Rubisco biosynthesis and assembly are down-regulated. However, a shift towards C4 photosynthesis is likely to occur as gene regulation is aimed to activate the primary CO₂ fixation to PEP (phosphoenolpyruvate). The analysis of "carbon metabolism" category revealed that G34 ecotype under salt stress induces the expression of glycolysis and Krebs cycle related genes, this being consistent with the hypothesis that some sort of salt avoidance might be occurred in A. donax G34 low ecotype. By comparing our results with findings obtained with other giant reed ecotype, we identified several differences in the response to salt that are in accordance with the possibility that heritable phenotypic differences among clones of A. donax might be accumulated especially in ecotypes originating from distant geographical areas, despite their asexual reproduction modality. Additionally, 26,838 simple sequence repeat (SSR) markers were identified and validated. This SSR dataset definitely expands the marker catalogue of A. donax facilitating the genotypic characterization of this species.

Comparative transcriptome analyses in contrasting onion (Allium cepa L.) genotypes for drought stress

Onion (Allium cepa L.) is an important vegetable crop widely grown for diverse culinary and nutraceutical properties. Being a shallow-rooted plant, it is prone to drought. In the present study, transcriptome sequencing of drought-tolerant (1656) and drought-sensitive (1627) onion genotypes was performed to elucidate the molecular basis of differential response to drought stress. A total of 123206 and 139252 transcripts (average transcript length: 690 bases) were generated after assembly for 1656 and 1627, respectively. Differential gene expression analyses revealed upregulation and downregulation of 1189 and 1180 genes, respectively, in 1656, whereas in 1627, upregulation and downregulation of 872 and 1292 genes, respectively, was observed. Genes encoding transcription factors, cytochrome P450, membrane transporters, and flavonoids, and those related to carbohydrate metabolism were found to exhibit a differential expression behavior in the tolerant and susceptible genotypes. The information generated can facilitate a better understanding of molecular mechanisms underlying drought response in onion.

De novo RNA sequencing analysis of Aeluropus littoralis halophyte plant under salinity stress

The study of salt tolerance mechanisms in halophyte plants can provide valuable information for crop breeding and plant engineering programs. The aim of the present study was to investigate whole transcriptome analysis of Aeluropus littoralis in response to salinity stress (200 and 400 mM NaCl) by de novo RNA-sequencing. To assemble the transcriptome, Trinity v2.4.0 and Bridger tools, were comparatively used with two k-mer sizes (25 and 32 bp). The de novo assembled transcriptome by Bridger (k-mer 32) was chosen as final assembly for subsequent analysis. In general, 103290 transcripts were obtained. The differential expression analysis (log₂^FC > 1 and FDR < 0.01) showed that 1861 transcripts expressed differentially, including169 up and 316 down-regulated transcripts in 200 mM NaCl treatment and 1035 up and 430 down-regulated transcripts in 400 mM NaCl treatment compared to control. In addition, 89 transcripts were common in both treatments. The most important over-represented terms in the GO analysis of differentially expressed genes (FDR < 0.05) were chitin response, response to abscisic acid, and regulation of jasmonic acid mediated signaling pathway under 400 mM NaCl treatment and cell cycle, cell division, and mitotic cell cycle process under 200 mM treatment. In addition, the phosphatidylcholine biosynthetic process term was common in both salt treatments. Interestingly, under 400 mM salt treatment, the PRC1 complex that contributes to chromatin remodeling was also enriched along with vacuole as a general salinity stress responsive cell component. Among enriched pathways, the MAPK signaling pathway (ko04016) and phytohormone signal transduction (ko04075) were significantly enriched in 400 mM NaCl treatment, whereas DNA replication (ko03032) was the only pathway that significantly enriched in 200 mM NaCl treatment. Finally, our findings indicate the salt-concentration depended responses of A. littoralis, which well-known salinity stress-related pathways are induced in 400 mM NaCl, while less considered pathways, e.g. cell cycle and DNA replication, are highlighted under 200 mM NaCl treatment.

Marginal Lands to Grow Novel Bio-Based Crops: A Plant Breeding Perspective

The biomass demand to fuel a growing global bio-based economy is expected to tremendously increase over the next decades, and projections indicate that dedicated biomass crops will satisfy a large portion of it. The establishment of dedicated biomass crops raises huge concerns, as they can subtract land that is required for food production, undermining food security. In this context, perennial biomass crops suitable for cultivation on marginal lands (MALs) raise attraction, as they could supply biomass without competing for land with food supply. While these crops withstand marginal conditions well, their biomass yield and quality do not ensure acceptable economic returns to farmers and cost-effective biomass conversion into bio-based products, claiming genetic improvement. However, this is constrained by the lack of genetic resources for most of these crops. Here we first review the advantages of cultivating novel perennial biomass crops on MALs, highlighting management practices to enhance the environmental and economic sustainability of these agro-systems. Subsequently, we discuss the preeminent breeding targets to improve the yield and quality of the biomass obtainable from these crops, as well as the stability of biomass production under MALs conditions. These targets include crop architecture and phenology, efficiency in the use of resources, lignocellulose composition in relation to bio-based applications, and tolerance to abiotic stresses. For each target trait, we outline optimal ideotypes, discuss the available breeding resources in the context of (orphan) biomass crops, and provide meaningful examples of genetic improvement. Finally, we discuss the available tools to breed novel perennial biomass crops. These comprise conventional breeding methods (recurrent selection and hybridization), molecular techniques to dissect the genetics of complex traits, speed up selection, and perform transgenic modification (genetic mapping, QTL and GWAS analysis, marker-assisted selection, genomic selection, transformation protocols), and novel high-throughput phenotyping platforms. Furthermore, novel tools to transfer genetic knowledge from model to orphan crops (i.e., universal markers) are also conceptualized, with the belief that their development will enhance the efficiency of plant breeding in orphan biomass crops, enabling a sustainable use of MALs for biomass provision.

Transcriptional, metabolic and DNA methylation changes underpinning the response of Arundo donax ecotypes to NaCl excess

Arundo donax ecotypes react differently to salinity, partly due to differences in constitutive defences and methylome plasticity. Arundo donax L. is a C3 fast-growing grass that yields high biomass under stress. To elucidate its ability to produce biomass under high salinity, we investigated short/long-term NaCl responses of three ecotypes through transcriptional, metabolic and DNA methylation profiling of leaves and roots. Prolonged salt treatment discriminated the sensitive ecotype 'Cercola' from the tolerant 'Domitiana' and 'Canneto' in terms of biomass. Transcriptional and metabolic responses to NaCl differed between the ecotypes. In roots, constitutive expression of ion transporter and stress-related transcription factors' genes was higher in 'Canneto' and 'Domitiana' than 'Cercola' and 21-day NaCl drove strong up-regulation in all ecotypes. In leaves, unstressed 'Domitiana' confirmed higher expression of the above genes, whose transcription was repressed in 'Domitiana' but induced in 'Cercola' following NaCl treatment. In all ecotypes, salinity increased proline, ABA and leaf antioxidants, paralleled by up-regulation of antioxidant genes in 'Canneto' and 'Cercola' but not in 'Domitiana', which tolerated a higher level of oxidative damage. Changes in DNA methylation patterns highlighted a marked capacity of the tolerant 'Domitiana' ecotype to adjust DNA methylation to salt stress. The reduced salt sensitivity of 'Domitiana' and, to a lesser extent, 'Canneto' appears to rely on a complex set of constitutively activated defences, possibly due to the environmental conditions of the site of origin, and on higher plasticity of the methylome. Our findings provide insights into the mechanisms of adaptability of A. donax ecotypes to salinity, offering new perspectives for the improvement of this species for cultivation in limiting environments.

Impact of high or low levels of phosphorus and high sodium in soils on productivity and stress tolerance of Arundo donax plants

The potential of Arundo donax to grow in degraded soils, characterized by excess of salinity (Na+), and phosphorus deficiency (-P) or excess (+P) also coupled with salinity (+NaP), was investigated by combining in vivo plant phenotyping, quantification of metabolites and ultrastructural imaging of leaves with a transcriptome-wide screening. Photosynthesis and growth were impaired by + Na, -P and + NaP. While + Na caused stomatal closure, enhanced biosynthesis of carotenoids, sucrose and isoprene and impaired anatomy of cell walls, +P negatively affected starch production and isoprene emission, and damaged chloroplasts. Finally, +NaP largely inhibited photosynthesis due to stomatal limitations, increased sugar content, induced/repressed a number of genes 10 time higher with respect to + P and + Na, and caused appearance of numerous and large plastoglobules and starch granules in chloroplasts. Our results show that A. donax is sensitive to unbalances of soil ion content, despite activation of defensive mechanisms that enhance plant resilience, growth and biomass production of A. donax under these conditions.

RNASeq analysis of giant cane reveals the leaf transcriptome dynamics under long-term salt stress

BACKGROUND

To compensate for the lack of information about the molecular mechanism involved in Arundo donax L. response to salt stress, we de novo sequenced, assembled and analyzed the A. donax leaf transcriptome subjected to two levels of long-term salt stress (namely, S3 severe and S4 extreme).

RESULTS

The picture that emerges from the identification of differentially expressed genes is consistent with a salt dose-dependent response. Hence, a deeper re-programming of the gene expression occurs in those plants grew at extreme salt level than in those subjected to severe salt stress, probably representing for them an "emergency" state. In particular, we analyzed clusters related to salt sensory and signaling, transcription factors, hormone regulation, Reactive Oxygen Species (ROS) scavenging, osmolyte biosynthesis and biomass production, all of them showing different regulation either versus untreated plants or between the two treatments. Importantly, the photosynthesis is strongly impaired in samples treated with both levels of salinity stress. However, in extreme salt conditions, a dramatic switch from C3 Calvin cycle to C4 photosynthesis is likely to occur, this probably being the more impressive finding of our work.

CONCLUSIONS

Considered the distinct response to salt doses, genes either involved in severe or in extreme salt response could constitute useful markers of the physiological status of A. donax to deepen our understanding of its biology and productivity in salinized soil. Finally, many of the unigenes identified in the present study have the potential to be used for the development of A. donax varieties with improved productivity and stress tolerance, in particular the knock out of the GTL1 gene acting as negative regulator of water use efficiency has been proposed as good target for genome editing.

Transcriptome analysis and codominant markers development in caper, a drought tolerant orphan crop with medicinal value

Caper (Capparis spinosa L.) is a xerophytic shrub cultivated for its flower buds and fruits, used as food and for their medicinal properties. Breeding programs and even proper taxonomic classification of the genus Capparis has been hampered so far by the lack of reliable genetic information and molecular markers. Here, we present the first genomic resource for C. spinosa, generated by transcriptomic approach and de novo assembly. The sequencing effort produced nearly 80 million clean reads assembled into 124,723 unitranscripts. Careful annotation and comparison with public databases revealed homologs to genes with a key role in important metabolic pathways linked to abiotic stress tolerance and bio-compounds production, such purine, thiamine and phenylpropanoid biosynthesis, α-linolenic acid and lipid metabolism. Additionally, a panel of genes involved in stomatal development/distribution and encoding for Stress Associated Proteins (SAPs) was also identified. We also used the transcriptomic data to uncover novel molecular markers for caper. Out of 50 SSRs tested, 14 proved polymorphic and represent the first set of SSR markers for the genus Capparis. This transcriptome will be an important contribution to future studies and breeding programs for this orphan crop, aiding to the development of improved varieties to sustain agriculture in arid conditions.

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.

Selection and evaluation of reference genes for qRT-PCR analysis in Euscaphis konishii Hayata based on transcriptome data

BACKGROUND

Quantitative real-time reverse transcription-polymerase chain reaction has been widely used in gene expression analysis, however, to have reliable and accurate results, reference genes are necessary to normalize gene expression under different experimental conditions. Several reliable reference genes have been reported in plants of Traditional Chinese Medicine, but none have been identified for Euscaphis konishii Hayata.

RESULTS

In this study, 12 candidate reference genes, including 3 common housekeeping genes and 9 novel genes based on E. konishii Hayata transcriptome data were selected and analyzed in different tissues (root, branch, leaf, capsule and seed), capsule and seed development stages. Expression stability was calculated using geNorm and NormFinder, the minimal number of reference genes required for accurate normalization was calculated by Vn/Vn + 1 using geNorm. EkEEF-5A-1 and EkADF2 were the two most stable reference genes for all samples, while EkGSTU1 and EkGAPDH were the most stable reference genes for tissue samples. For seed development stages, EkGAPDH and EkEEF-5A-1 were the most stable genes, whereas EkGSTU1 and EkGAPDH were identified as the two most stable genes in the capsule development stages. Two reference genes were sufficient to normalize gene expression across all sample sets.

CONCLUSION

Results of this study revealed that suitable reference genes should be selected for different experimental samples, and not all the common reference genes are suitable for different tissue samples and/or experimental conditions. In this study, we present the first data of reference genes selection for E. konishii Hayata based on transcriptome data, our data will facilitate further studies in molecular biology and gene function on E. konishii Hayata and other closely related species.

RNA-seq-based genome annotation and identification of long-noncoding RNAs in the grapevine cultivar 'Riesling'

BACKGROUND

The technological advances of RNA-seq and de novo transcriptome assembly have enabled genome annotation and transcriptome profiling in highly heterozygous species such as grapevine (Vitis vinifera L.). This work is an attempt to utilize a de novo-assembled transcriptome of the V. vinifera cultivar 'Riesling' to improve annotation of the grapevine reference genome sequence.

RESULTS

Here we show that the transcriptome assembly of a single V. vinifera cultivar is insufficient for a complete genome annotation of the grapevine reference genome constructed from V. vinifera PN40024. Further, we provide evidence that the gene models we identified cannot be completely anchored to the previously published V. vinifera PN40024 gene models. In addition to these findings, we present a computational pipeline for the de novo identification of lncRNAs. Our results demonstrate that, in grapevine, lncRNAs are significantly different from protein coding transcripts in such metrics as length, GC-content, minimum free energy, and length-corrected minimum free energy.

CONCLUSIONS

In grapevine, high-level heterozygosity necessitates that transcriptome characterization be based on cultivar-specific reference genome sequences. Our results strengthen the hypothesis that lncRNAs have thermodynamically different properties than protein-coding RNAs. The analyses of both coding and non-coding RNAs will be instrumental in uncovering inter-cultivar variation in wild and cultivated grapevine species.