Shoot transcriptome of the giant reed, Arundo donax

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.

Functional Diversity of HemN-like Proteins

HemN is a radical S-adenosylmethionine (SAM) enzyme that catalyzes the anaerobic oxidative decarboxylation of coproporphyrinogen III to produce protoporphyrinogen IX, a key intermediate in heme biosynthesis. Proteins homologous to HemN (HemN-like proteins) are widespread in both prokaryotes and eukaryotes. Although these proteins are in most cases annotated as anaerobic coproporphyrinogen III oxidases (CPOs) in the public database, many of them are actually not CPOs but have diverse functions such as methyltransferases, cyclopropanases, heme chaperones, to name a few. This Perspective discusses the recent advances in the understanding of HemN-like proteins, and particular focus is placed on the diverse chemistries and functions of this growing protein family.

Novel genome characteristics contribute to the invasiveness of Phragmites australis (common reed)

The rapid invasion of the non‐native Phragmites australis (Poaceae, subfamily Arundinoideae) is a major threat to native wetland ecosystems in North America and elsewhere. We describe the first reference genome for P. australis and compare invasive (ssp. australis) and native (ssp. americanus) genotypes collected from replicated populations across the Laurentian Great Lakes to deduce genomic bases driving its invasive success. Here, we report novel genomic features including a Phragmites lineage‐specific whole genome duplication, followed by gene loss and preferential retention of genes associated with transcription factors and regulatory functions in the remaining duplicates. Comparative transcriptomic analyses revealed that genes associated with biotic stress and defence responses were expressed at a higher basal level in invasive genotypes, but native genotypes showed a stronger induction of defence responses when challenged by a fungal endophyte. The reference genome and transcriptomes, combined with previous ecological and environmental data, add to our understanding of mechanisms leading to invasiveness and support the development of novel, genomics‐assisted management approaches for invasive Phragmites.

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.

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.

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.

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.

In silico identification and characterization of a diverse subset of conserved microRNAs in bioenergy crop Arundo donax L

MicroRNAs (miRNAs) are small non-coding RNA molecules involved in the post-transcriptional regulation of gene expression in plants. Arundo donax L. is a perennial C₃ grass considered one of the most promising bioenergy crops. Despite its relevance, many fundamental aspects of its biology still remain to be elucidated. In the present study we carried out the first in silico mining and tissue-specific characterization of microRNAs and their putative targets in A. donax. We identified a total of 141 miRNAs belonging to 14 families along with the corresponding primary miRNAs, precursor miRNAs and a total of 462 high-confidence predicted targets and novel target sites were validated by 5′-race. Gene Ontology functional annotation showed that miRNA targets are constituted mainly by transcription factors, but three of the newly validated targets are enzymes involved in novel functions like RNA editing, acyl lipid metabolism and post-Golgi trafficking. Folding variability of pre-miRNA loops and phylogenetic analyses indicate variable selective pressure acting on the different miRNA families. The set of miRNAs identified in this study will pave the road to further miRNA research in Arundo donax and contribute towards a better understanding of miRNA-mediated gene regulatory processes in other bioenergy crops.

Selection of reference genes suitable for normalization of qPCR data under abiotic stresses in bioenergy crop Arundo donax L

Suitable reference gene selection in qRT-PCR is a key pre-requisite to produce reliable data in gene expression analyses. In this study, novel primers for six commonly used reference genes (AC1, TLF, Act2, TUB α, EF-1α and GAPDH) plus two new candidates (pDUF221 and RPN6) were designed and comparatively tested for expression stability under abiotic stresses (osmotic, heavy metal and heat shock) in shoot, root and their combination of Arundo donax L., a raising non-food energy crop. Expression stability rankings from the most to the least stable gene in each condition and in two tissues (young shoots and roots) were generated with geNorm, NormFinder and BestKeeper programs. All programs provided similar rankings and, strikingly, in most cases identified one of the new candidates, RPN6, as the most suitable reference gene. This novel set of reliable references allows to choose either the best combination of reference genes across multiple stress/organ conditions or to select condition-specific genes that can improve the quality of qRT-PCR analysis. This work provides a solid basis for the functional characterization of A. donax, by enabling accurate quantification of the transcriptional responsiveness under a series of common stress conditions of any gene of interest in this promising biomass/bioenergy species.

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.

Methods for suspension culture, protoplast extraction, and transformation of high-biomass yielding perennial grass Arundo donax

Arundo donax L. is a promising biofuel feedstock in the Mediterranean region. Despite considerable interest in its genetic improvement, Arundo tissue culture and transformation remains arduous. The authors developed methodologies for cell- and tissue culture and genetic engineering in Arundo. A media screen was conducted, and a suspension culture was established using callus induced from stem axillary bud explants. DBAP medium, containing 9 µM 2,4-D and 4.4 µM BAP, was found to be the most effective medium among those tested for inducing cell suspension cultures, which resulted in a five-fold increase in tissue mass over 14 days. In contrast, CIM medium containing 13 µM 2,4-D, resulted in just a 1.4-fold increase in mass over the same period. Optimized suspension cultures were superior to previously-described solidified medium-based callus culture methods for tissue mass increase. Suspension cultures proved to be very effective for subsequent protoplast isolation. Protoplast electroporation resulted in a 3.3 ± 1.5% transformation efficiency. A dual fluorescent reporter gene vector enabled the direct comparison of the CAMV 35S promoter with the switchgrass ubi2 promoter in single cells of Arundo. The switchgrass ubi2 promoter resulted in noticeably higher reporter gene expression compared with that conferred by the 35S promoter in Arundo.