An RNA-Seq transcriptome analysis of orthophosphate-deficient white lupin reveals novel insights into phosphorus acclimation in plants

Potential functions of microRNAs in starch metabolism and development revealed by miRNA transcriptome profiling of cassava cultivars and their wild progenitor

BACKGROUND

MicroRNAs (miRNAs) are small (approximately 21 nucleotide) non-coding RNAs that are key post-transcriptional gene regulators in eukaryotic organisms. More than 100 cassava miRNAs have been identified in a conservation analysis and a repertoire of cassava miRNAs have also been characterised by next-generation sequencing (NGS) in recent studies. Here, using NGS, we profiled small non-coding RNAs and mRNA genes in two cassava cultivars and their wild progenitor to identify and characterise miRNAs that are potentially involved in plant growth and starch biosynthesis.

RESULTS

Six small RNA and six mRNA libraries from leaves and roots of the two cultivars, KU50 and Arg7, and their wild progenitor, W14, were subjected to NGS. Analysis of the sequencing data revealed 29 conserved miRNA families and 33 new miRNA families. Together, these miRNAs potentially targeted a total of 360 putative target genes. Whereas 16 miRNA families were highly expressed in cultivar leaves, another 13 miRNA families were highly expressed in storage roots of cultivars. Co-expression analysis revealed that the expression level of some targets had negative relationship with their corresponding miRNAs in storage roots and leaves; these targets included MYB33, ARF10, GRF1, RD19, APL2, NF-YA3 and SPL2, which are known to be involved in plant development, starch biosynthesis and response to environmental stimuli.

CONCLUSION

The identified miRNAs, target mRNAs and target gene ontology annotation all shed light on the possible functions of miRNAs in Manihot species. The differential expression of miRNAs between cultivars and their wild progenitor, together with our analysis of GO annotation and confirmation of miRNA: target pairs, might provide insight into know the differences between wild progenitor and cultivated cassava.

Structure, expression profile and phylogenetic inference of chalcone isomerase-like genes from the narrow-leafed lupin (Lupinus angustifolius L.) genome

Lupins, like other legumes, have a unique biosynthesis scheme of 5-deoxy-type flavonoids and isoflavonoids. A key enzyme in this pathway is chalcone isomerase (CHI), a member of CHI-fold protein family, encompassing subfamilies of CHI1, CHI2, CHI-like (CHIL), and fatty acid-binding (FAP) proteins. Here, two Lupinus angustifolius (narrow-leafed lupin) CHILs, LangCHIL1 and LangCHIL2, were identified and characterized using DNA fingerprinting, cytogenetic and linkage mapping, sequencing and expression profiling. Clones carrying CHIL sequences were assembled into two contigs. Full gene sequences were obtained from these contigs, and mapped in two L. angustifolius linkage groups by gene-specific markers. Bacterial artificial chromosome fluorescence in situ hybridization approach confirmed the localization of two LangCHIL genes in distinct chromosomes. The expression profiles of both LangCHIL isoforms were very similar. The highest level of transcription was in the roots of the third week of plant growth; thereafter, expression declined. The expression of both LangCHIL genes in leaves and stems was similar and low. Comparative mapping to reference legume genome sequences revealed strong syntenic links; however, LangCHIL2 contig had a much more conserved structure than LangCHIL1. LangCHIL2 is assumed to be an ancestor gene, whereas LangCHIL1 probably appeared as a result of duplication. As both copies are transcriptionally active, questions arise concerning their hypothetical functional divergence. Screening of the narrow-leafed lupin genome and transcriptome with CHI-fold protein sequences, followed by Bayesian inference of phylogeny and cross-genera synteny survey, identified representatives of all but one (CHI1) main subfamilies. They are as follows: two copies of CHI2, FAPa2 and CHIL, and single copies of FAPb and FAPa1. Duplicated genes are remnants of whole genome duplication which is assumed to have occurred after the divergence of Lupinus, Arachis, and Glycine.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

The Role of Ethylene in Plant Adaptations for Phosphate Acquisition in Soils - A Review

Although a role of ethylene in the regulation of senescence and plant stress responses in general has a long history, a possible involvement in the regulation of adaptive responses to nutrient deficiencies has been mainly investigated since the last two decades. In the case of plant responses to phosphate (Pi) starvation, ethylene was identified as a modulator of adaptive responses in root growth and morphology. The molecular base of these adaptations has been elucidated in supplementation studies with ethylene precursors and antagonists, as well as analysis of mutants and transgenic plants with modified ethylene biosynthesis and responsiveness, using mainly Arabidopsis thaliana as a model plant. However, increasing evidence suggests that apart from root growth responses, ethylene may be involved in various additional plant adaptations to Pi limitation including Pi mobilization in the rhizosphere, Pi uptake and internal Pi recycling. The ethylene-mediated responses are frequently characterized by high genotypic variability and may partially share common pathways in different nutrient limitations.

Remnants of the Legume Ancestral Genome Preserved in Gene-Rich Regions: Insights from Lupinus angustifolius Physical, Genetic, and Comparative Mapping

The narrow-leafed lupin (Lupinus angustifolius) was recently considered as a legume reference species. Genetic resources have been developed, including a draft genome sequence, linkage maps, nuclear DNA libraries, and cytogenetic chromosome-specific landmarks. Here, we used a complex approach, involving DNA fingerprinting, sequencing, genetic mapping, and molecular cytogenetics, to localize and analyze L. angustifolius gene-rich regions (GRRs). A L. angustifolius genomic bacterial artificial chromosome (BAC) library was screened with short sequence repeat (SSR)-based probes. Selected BACs were fingerprinted and assembled into contigs. BAC-end sequence (BES) annotation allowed us to choose clones for sequencing, targeting GRRs. Additionally, BESs were aligned to the scaffolds of the genome sequence. The genetic map was supplemented with 35 BES-derived markers, distributed in 14 linkage groups and tagging 37 scaffolds. The identified GRRs had an average gene density of 19.6 genes/100 kb and physical-to-genetic distance ratios of 11 to 109 kb/cM. Physical and genetic mapping was supported by multi-BAC-fluorescence in situ hybridization (FISH), and five new linkage groups were assigned to the chromosomes. Syntenic links to the genome sequences of five legume species (Medicago truncatula, Glycine max, Lotus japonicus, Phaseolus vulgaris, and Cajanus cajan) were identified. The comparative mapping of the two largest lupin GRRs provides novel evidence for ancient duplications in all of the studied species. These regions are conserved among representatives of the main clades of Papilionoideae. Furthermore, despite the complex evolution of legumes, some segments of the nuclear genome were not substantially modified and retained their quasi-ancestral structures. Cytogenetic markers anchored in these regions constitute a platform for heterologous mapping of legume genomes.

Root architecture responses: in search of phosphate

Soil phosphate represents the only source of phosphorus for plants and, consequently, is its entry into the trophic chain. This major component of nucleic acids, phospholipids, and energy currency of the cell (ATP) can limit plant growth because of its low mobility in soil. As a result, root responses to low phosphate favor the exploration of the shallower part of the soil, where phosphate tends to be more abundant, a strategy described as topsoil foraging. We will review the diverse developmental strategies that can be observed among plants by detailing the effect of phosphate deficiency on primary and lateral roots. We also discuss the formation of cluster roots: an advanced adaptive strategy to cope with low phosphate availability observed in a limited number of species. Finally, we will put this work into perspective for future research directions.

Comparative transcriptomic analysis of the response to cold acclimation in Eucalyptus dunnii

Eucalyptus dunnii is an important macrophanerophyte with high economic value. However, low temperature stress limits its productivity and distribution. To study the cold response mechanisms of E. dunnii, 5 cDNA libraries were constructed from mRNA extracted from leaves exposed to cold stress for varying lengths of time and were evaluated by RNA-Seq analysis. The assembly of the Illumina datasets was optimized using various assembly programs and parameters. The final optimized assembly generated 205,325 transcripts with an average length of 1,701 bp and N50 of 2,627 bp, representing 349.38 Mb of the E. dunnii transcriptome. Among these transcripts, 134,358 transcripts (65.4%) were annotated in the Nr database. According to the differential analysis results, most transcripts were up-regulated as the cold stress prolonging, suggesting that these transcripts may be involved in the response to cold stress. In addition, the cold-relevant GO categories, such as ‘response to stress’ and ‘translational initiation’, were the markedly enriched GO terms. The assembly of the E. dunnii gene index and the GO classification performed in this study will serve as useful genomic resources for the genetic improvement of E. dunnii and also provide insights into the molecular mechanisms of cold acclimation in E. dunnii.

Expression of novel cytosolic malate dehydrogenases (cMDH) in Lupinus angustifolius nodules during phosphorus starvation

During P deficiency, the increased activity of malate dehydrogenase (MDH, EC 1.1.1.37) can lead to malate accumulation. Cytosolic- and nodule-enhanced MDH (cMDH and neMDH, respectively) are known isoforms, which contribute to MDH activity in root nodules. The aim of this study was to investigate the role of the cMDH isoforms in nodule malate supply under P deficiency. Nodulated lupins (Lupinus angustifolius var. Tanjil) were hydroponically grown at adequate P (+P) or low P (-P). Total P concentration in nodules decreased under P deficiency, which coincided with an increase in total MDH activity. A consequence of higher MDH activity was the enhanced accumulation of malate derived from dark CO2 fixation via PEPC and not from pyruvate. Although no measurable neMDH presence could be detected via PCR, gene-specific primers detected two 1kb amplicons of cMDH, designated LangMDH1 (corresponding to +P, HQ690186) and LangMDH2 (corresponding to -P, HQ690187), respectively. Sequencing analyses of these cMDH amplicons showed them to be 96% identical on an amino acid level. There was a high degree of diversification between proteins detected in this study and other known MDH proteins, particularly those from other leguminous plants. Enhanced malate synthesis in P-deficient nodules was achieved via increased anaplerotic CO2 fixation and subsequent higher MDH activities. Novel isoforms of cytosolic MDH may be involved, as shown by gene expression of specific genes under P deficiency.

An RNA-Seq based gene expression atlas of the common bean

BACKGROUND

Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species.

RESULTS

Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form.

CONCLUSIONS

These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

Transcriptomic profiling revealed an important role of cell wall remodeling and ethylene signaling pathway during salt acclimation in Arabidopsis

Plants can successfully improve their resistance to previously lethal salinity stress by a short exposure to low levels of salt stress, a process known as salt acclimation (SA). In spite of its fundamental significance in theoretical study and agricultural practice, the molecular mechanisms underlying plant SA remain elusive. In this study, we found that salt acclimated Arabidopsis young seedlings can survive subsequent 200 mM NaCl stress. RNA-seq was performed to analyze the genome-wide transcriptional response under SA conditions. Among 518 differentially expressed genes (DEGs) under SA, 366 up-regulated genes were enriched for cell wall biosynthesis, osmoregulation, oxidative stress, or transcription factors. Seven DEGs participate in the synthesis of lignin and 24 DEGs encode plant cell wall proteins, suggesting the importance of cell wall remodeling under SA. Furthermore, in comparison to non-acclimated salt stress, 228 of 245 DEGs were repressed by acclimated salt stress, including many genes related to ethylene biosynthesis and signaling pathway. In addition, MAPK6, a major component of the ethylene signaling pathway, was found to play a crucial role in SA. Our transcriptomic analysis has provided important insight on the roles of transcription factors, cell wall remodeling, and the ethylene biosynthesis and signaling pathways during SA in Arabidopsis.

De novo assembly and characterization of stress transcriptome and regulatory networks under temperature, salt and hormone stresses in Lilium lancifolium

Plants have continually confrontation with different abiotic stresses, including salt, low temperature, drought or hormone stress. The plants acclimate to the environmental stresses relating with the falls of the molecular mesh including the stress signal receiver, signal transcriptional regulation and the expression of functional and structure genes. Using the RNA-seq, we carried out a transcriptional analysis under cold treatment for investigating a profound comprehension of the signal network and molecular metabolisms reaction included in abiotic stress reaction for Lilium lancifolium. Our study identified 18,722 unigenes had demonstrated the resemblance to the known exact proteins in the Swiss-Prot protein database and classified them by Gene ontology into three primary kinds: cellular component, biological process, and molecular function, and then 15,898 unigenes aligned to existing sequences in the KEGG databases. Based on the transcriptome results of cold stress, more stress-related genes were identified and analyzed of their expressions in other abiotic stress treatments as 37 °C, ABA, JA and Na. Meanwhile, bioinformatics qRT-PCR analyses of stress genes as LlDREB1, LlAP2, LlNAC1, LlHOT, LlR2R3-MYB and LlCDPK revealed that novel candidate genes encoding ethylene responsive transporters and serine/threonine receptor-like kinases, which contributed to speculate the signal regulation pathway during the abiotic stresses; engineering genes could also boost the tolerance to stress, as protected and maintained the function and structure of cellular components. Our research conjectured the abiotic stress signal transduction pathway and identified the expected key ingredients regulating the stress tolerance in Lilium lancifolium, which would enable the in-depth molecular exploration of stress-tolerance mechanisms in lily.

Genomics and molecular breeding in lesser explored pulse crops: current trends and future opportunities

Pulses are multipurpose crops for providing income, employment and food security in the underprivileged regions, notably the FAO-defined low-income food-deficit countries. Owing to their intrinsic ability to endure environmental adversities and the least input/management requirements, these crops remain central to subsistence farming. Given their pivotal role in rain-fed agriculture, substantial research has been invested to boost the productivity of these pulse crops. To this end, genomic tools and technologies have appeared as the compelling supplement to the conventional breeding. However, the progress in minor pulse crops including dry beans (Vigna spp.), lupins, lablab, lathyrus and vetches has remained unsatisfactory, hence these crops are often labeled as low profile or lesser researched. Nevertheless, recent scientific and technological breakthroughs particularly the next generation sequencing (NGS) are radically transforming the scenario of genomics and molecular breeding in these minor crops. NGS techniques have allowed de novo assembly of whole genomes in these orphan crops. Moreover, the availability of a reference genome sequence would promote re-sequencing of diverse genotypes to unlock allelic diversity at a genome-wide scale. In parallel, NGS has offered high-resolution genetic maps or more precisely, a robust genetic framework to implement whole-genome strategies for crop improvement. As has already been demonstrated in lupin, sequencing-based genotyping of the representative sample provided access to a number of functionally-relevant markers that could be deployed straight away in crop breeding programs. This article attempts to outline the recent progress made in genomics of these lesser explored pulse crops, and examines the prospects of genomics assisted integrated breeding to enhance and stabilize crop yields.

The temporal transcriptomic response of Pinus massoniana seedlings to phosphorus deficiency

BACKGROUND

Phosphorus (P) is an essential macronutrient for plant growth and development. Several genes involved in phosphorus deficiency stress have been identified in various plant species. However, a whole genome understanding of the molecular mechanisms involved in plant adaptations to low P remains elusive, and there is particularly little information on the genetic basis of these acclimations in coniferous trees. Masson pine (Pinus massoniana) is grown mainly in the tropical and subtropical regions in China, many of which are severely lacking in inorganic phosphate (Pi). In previous work, we described an elite P. massoniana genotype demonstrating a high tolerance to Pi-deficiency.

METHODOLOGY/PRINCIPAL FINDINGS

To further investigate the mechanism of tolerance to low P, RNA-seq was performed to give an idea of extent of expression from the two mixed libraries, and microarray whose probes were designed based on the unigenes obtained from RNA-seq was used to elucidate the global gene expression profiles for the long-term phosphorus starvation. A total of 70,896 unigenes with lengths ranging from 201 to 20,490 bp were assembled from 112,108,862 high quality reads derived from RNA-Seq libraries. We identified 1,396 and 943 transcripts that were differentially regulated (P<0.05) under P1 (0.01 mM P) and P2 (0.06 mM P) Pi-deficiency conditions, respectively. Numerous transcripts were consistently differentially regulated under Pi deficiency stress, many of which were also up- or down-regulated in other species under the corresponding conditions, and are therefore ideal candidates for monitoring the P status of plants. The results also demonstrated the impact of different Pi starvation levels on global gene expression in Masson pine.

CONCLUSIONS/SIGNIFICANCE

To our knowledge, this work provides the first insight into the molecular mechanisms involved in acclimation to long-term Pi starvation and different Pi availability levels in coniferous trees.

The regulatory network of cluster-root function and development in phosphate-deficient white lupin (Lupinus albus) identified by transcriptome sequencing

Lupinus albus serves as model plant for root-induced mobilization of sparingly soluble soil phosphates via the formation of cluster-roots (CRs) that mediate secretion of protons, citrate, phenolics and acid phosphatases (APases). This study employed next-generation sequencing to investigate the molecular mechanisms behind these complex adaptive responses at the transcriptome level. We compared different stages of CR development, including pre-emergent (PE), juvenile (JU) and the mature (MA) stages. The results confirmed that the primary metabolism underwent significant modifications during CR maturation, promoting the biosynthesis of organic acids, as had been deduced from physiological studies. Citrate catabolism was downregulated, associated with citrate accumulation in MA clusters. Upregulation of the phenylpropanoid pathway reflected the accumulation of phenolics. Specific transcript expression of ALMT and MATE transporter genes correlated with the exudation of citrate and flavonoids. The expression of transcripts related to nucleotide degradation and APases in MA clusters coincided with the re-mobilization and hydrolysis of organic phosphate resources. Most interestingly, hormone-related gene expression suggested a central role of ethylene during CR maturation. This was associated with the upregulation of the iron (Fe)-deficiency regulated network that mediates ethylene-induced expression of Fe-deficiency responses in other species. Finally, transcripts related to abscisic acid and jasmonic acid were upregulated in MA clusters, while auxin- and brassinosteroid-related genes and cytokinin receptors were most strongly expressed during CR initiation. Key regulations proposed by the RNA-seq data were confirmed by quantitative real-time polymerase chain reaction (RT-qPCR) and some physiological analyses. A model for the gene network regulating CR development and function is presented.

Suppression of Photosynthetic Gene Expression in Roots Is Required for Sustained Root Growth under Phosphate Deficiency

Plants cope with inorganic phosphate (Pi) deficiencies in their environment by adjusting their developmental programs and metabolic activities. For Arabidopsis (Arabidopsis thaliana), the developmental responses include the inhibition of primary root growth and the enhanced formation of lateral roots and root hairs. Pi deficiency also inhibits photosynthesis by suppressing the expression of photosynthetic genes. Early studies showed that photosynthetic gene expression was also suppressed in Pi-deficient roots, a nonphotosynthetic organ; however, the biological relevance of this phenomenon remains unknown. In this work, we characterized an Arabidopsis mutant, hypersensitive to Pi starvation7 (hps7), that is hypersensitive to Pi deficiency; the hypersensitivity includes an increased inhibition of root growth. HPS7 encodes a tyrosylprotein sulfotransferase. Accumulation of HPS7 proteins in root tips is enhanced by Pi deficiency. Comparative RNA sequencing analyses indicated that the expression of many photosynthetic genes is activated in roots of hps7. Under Pi deficiency, the expression of photosynthetic genes in hps7 is further increased, which leads to enhanced accumulation of chlorophyll, starch, and sucrose. Pi-deficient hps7 roots also produce a high level of reactive oxygen species. Previous research showed that the overexpression of GOLDEN-like (GLK) transcription factors in transgenic Arabidopsis activates photosynthesis in roots. The GLK overexpressing (GLK OX) lines also exhibit increased inhibition of root growth under Pi deficiency. The increased inhibition of root growth in hps7 and GLK OX lines by Pi deficiency was completely reversed by growing the plants in the dark. Based on these results, we propose that suppression of photosynthetic gene expression is required for sustained root growth under Pi deficiency.

Comparative transcriptome profiling of two Tibetan wild barley genotypes in responses to low potassium

Potassium (K) deficiency is one of the major factors affecting crop growth and productivity. Development of low-K tolerant crops is an effective approach to solve the nutritional deficiency in agricultural production. Tibetan annual wild barley is rich in genetic diversity and can grow normally under poor soils, including low-K supply. However, the molecular mechanism about low K tolerance is still poorly understood. In this study, Illumina RNA-Sequencing was performed using two Tibetan wild barley genotypes differing in low K tolerance (XZ153, tolerant and XZ141, sensitive), to determine the genotypic difference in transcriptome profiling. We identified a total of 692 differentially expressed genes (DEGs) in two genotypes at 6 h and 48 h after low-K treatment, including transcription factors, transporters and kinases, oxidative stress and hormone signaling related genes. Meanwhile, 294 low-K tolerant associated DEGs were assigned to transporter and antioxidant activities, stimulus response, and other gene ontology (GO), which were mainly involved in starch and sucrose metabolism, lipid metabolism and ethylene biosynthesis. Finally, a hypothetical model of low-K tolerance mechanism in XZ153 was presented. It may be concluded that wild barley accession XZ153 has a higher capability of K absorption and use efficiency than XZ141 under low K stress. A rapid response to low K stress in XZ153 is attributed to its more K uptake and accumulation in plants, resulting in higher low K tolerance. The ethylene response pathway may account for the genotypic difference in low-K tolerance.

De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification

BACKGROUND

Red clover (Trifolium pratense L.) is a versatile forage crop legume, which can tolerate a variety of soils and is suitable for silage production for winter feed and for grazing. It is one of the most important forage legumes in temperate livestock agriculture. Its beneficial attributes include ability to fix nitrogen, improve soil and provide protein rich animal feed. It is however, a short-lived perennial providing good biomass yield for two or three years. Improved persistency is thus a major breeding target. Better water-stress tolerance is one of the key factors influencing persistency, but little is known about how red clover tolerates water stress.

RESULTS

Plants from a full sib mapping family were used in a drought experiment, in which the growth rate and relative water content (RWC) identified two pools of ten plants contrasting in their tolerance to drought. Key metabolites were measured and RNA-Seq analysis was carried out on four bulked samples: the two pools sampled before and after drought. Massively parallel sequencing was used to analyse the bulked RNA samples. A de novo transcriptome reconstruction based on the RNA-Seq data was made, resulting in 45181 contigs, representing 'transcript tags'. These transcript tags were annotated with gene ontology (GO) terms. One of the most striking results from the expression analysis was that the drought sensitive plants were characterised by having approximately twice the number of differentially expressed transcript tags than the tolerant plants after drought. This difference was evident in most of the major GO terms. Before onset of drought the sensitive plants overexpressed a number of genes annotated as senescence-related. Furthermore, the concentration of three metabolites, particularly pinitol, but also proline and malate increased in leaves after drought stress.

CONCLUSIONS

This de novo assembly of a red clover transcriptome from leaf material of droughted and non-droughted plants provides a rich source for gene identification, single nucleotide polymorphisms (SNP) and short sequence repeats (SSR). Comparison of gene expression levels between pools and treatments identified candidate genes for further analysis of the genetic basis of drought tolerance in red clover.

Legume genomics: understanding biology through DNA and RNA sequencing

BACKGROUND

The legume family (Leguminosae) consists of approx. 17 000 species. A few of these species, including, but not limited to, Phaseolus vulgaris, Cicer arietinum and Cajanus cajan, are important dietary components, providing protein for approx. 300 million people worldwide. Additional species, including soybean (Glycine max) and alfalfa (Medicago sativa), are important crops utilized mainly in animal feed. In addition, legumes are important contributors to biological nitrogen, forming symbiotic relationships with rhizobia to fix atmospheric N2 and providing up to 30 % of available nitrogen for the next season of crops. The application of high-throughput genomic technologies including genome sequencing projects, genome re-sequencing (DNA-seq) and transcriptome sequencing (RNA-seq) by the legume research community has provided major insights into genome evolution, genomic architecture and domestication.

SCOPE AND CONCLUSIONS

This review presents an overview of the current state of legume genomics and explores the role that next-generation sequencing technologies play in advancing legume genomics. The adoption of next-generation sequencing and implementation of associated bioinformatic tools has allowed researchers to turn each species of interest into their own model organism. To illustrate the power of next-generation sequencing, an in-depth overview of the transcriptomes of both soybean and white lupin (Lupinus albus) is provided. The soybean transcriptome focuses on analysing seed development in two near-isogenic lines, examining the role of transporters, oil biosynthesis and nitrogen utilization. The white lupin transcriptome analysis examines how phosphate deficiency alters gene expression patterns, inducing the formation of cluster roots. Such studies illustrate the power of next-generation sequencing and bioinformatic analyses in elucidating the gene networks underlying biological processes.

RNA-Seq transcriptome analysis to identify genes involved in metabolism-based diclofop resistance in Lolium rigidum

Weed control failures due to herbicide resistance are an increasing and worldwide problem that significantly affect crop yields. Metabolism-based herbicide resistance (referred to as metabolic resistance) in weeds is not well characterized at the genetic level. An RNA-Seq transcriptome analysis was used to find candidate genes that conferred metabolic resistance to the herbicide diclofop in a diclofop-resistant population (R) of the major global weed Lolium rigidum. A reference cDNA transcriptome (19 623 contigs) was assembled and assigned putative annotations. Global gene expression was measured using Illumina reads from untreated control, adjuvant-only control, and diclofop treatment of R and susceptible (S). Contigs that showed constitutive expression differences between untreated R and untreated S were selected for further validation analysis, including 11 contigs putatively annotated as cytochrome P450 (CytP450), glutathione transferase (GST), or glucosyltransferase (GT), and 17 additional contigs with annotations related to metabolism or signal transduction. In a forward genetics validation experiment, nine contigs had constitutive up-regulation in R individuals from a segregating F2 population, including three CytP450, one nitronate monooxygenase (NMO), three GST, and one GT. Principal component analysis using these nine contigs differentiated F2 -R from F2 -S individuals. In a physiological validation experiment in which 2,4-D pre-treatment induced diclofop protection in S individuals due to increased metabolism, seven of the nine genetically validated contigs were induced significantly. Four contigs (two CytP450, NMO, and GT) were consistently highly expressed in nine field-evolved metabolic resistant L. rigidum populations. These four contigs were strongly associated with the resistance phenotype and are major candidates for contributing to metabolic diclofop resistance.

Post-transcriptional regulation in root development

Plants constantly adapt their root system to the changing environmental conditions. This developmental plasticity is underpinned by changes in the profile of the mRNA expressed. Here we review how post-transcriptional modulation of gene expression control root development and growth. In particular we focus on the role of small RNA-mediated post-transcriptional regulation processes. Small RNAs play an important role in fine tuning gene expression during root formation and patterning, development of lateral organs and symbiosis, nutrient homeostasis, and other stress-related responses. We also highlight the impact of alternative splicing on root development and the establishment of symbiotic structures as well as the emerging role of long noncoding RNAs in root physiology.

RNA sequencing analysis reveals transcriptomic variations in tobacco (Nicotiana tabacum) leaves affected by climate, soil, and tillage factors

The growth and development of plants are sensitive to their surroundings. Although numerous studies have analyzed plant transcriptomic variation, few have quantified the effect of combinations of factors or identified factor-specific effects. In this study, we performed RNA sequencing (RNA-seq) analysis on tobacco leaves derived from 10 treatment combinations of three groups of ecological factors, i.e., climate factors (CFs), soil factors (SFs), and tillage factors (TFs). We detected 4980, 2916, and 1605 differentially expressed genes (DEGs) that were affected by CFs, SFs, and TFs, which included 2703, 768, and 507 specific and 703 common DEGs (simultaneously regulated by CFs, SFs, and TFs), respectively. GO and KEGG enrichment analyses showed that genes involved in abiotic stress responses and secondary metabolic pathways were overrepresented in the common and CF-specific DEGs. In addition, we noted enrichment in CF-specific DEGs related to the circadian rhythm, SF-specific DEGs involved in mineral nutrient absorption and transport, and SF- and TF-specific DEGs associated with photosynthesis. Based on these results, we propose a model that explains how plants adapt to various ecological factors at the transcriptomic level. Additionally, the identified DEGs lay the foundation for future investigations of stress resistance, circadian rhythm and photosynthesis in tobacco.

Transcriptomic analysis of Petunia hybrida in response to salt stress using high throughput RNA sequencing

Salinity and drought stress are the primary cause of crop losses worldwide. In sodic saline soils sodium chloride (NaCl) disrupts normal plant growth and development. The complex interactions of plant systems with abiotic stress have made RNA sequencing a more holistic and appealing approach to study transcriptome level responses in a single cell and/or tissue. In this work, we determined the Petunia transcriptome response to NaCl stress by sequencing leaf samples and assembling 196 million Illumina reads with Trinity software. Using our reference transcriptome we identified more than 7,000 genes that were differentially expressed within 24 h of acute NaCl stress. The proposed transcriptome can also be used as an excellent tool for biological and bioinformatics in the absence of an available Petunia genome and it is available at the SOL Genomics Network (SGN) http://solgenomics.net. Genes related to regulation of reactive oxygen species, transport, and signal transductions as well as novel and undescribed transcripts were among those differentially expressed in response to salt stress. The candidate genes identified in this study can be applied as markers for breeding or to genetically engineer plants to enhance salt tolerance. Gene Ontology analyses indicated that most of the NaCl damage happened at 24 h inducing genotoxicity, affecting transport and organelles due to the high concentration of Na+ ions. Finally, we report a modification to the library preparation protocol whereby cDNA samples were bar-coded with non-HPLC purified primers, without affecting the quality and quantity of the RNA-seq data. The methodological improvement presented here could substantially reduce the cost of sample preparation for future high-throughput RNA sequencing experiments.

The angiosperm phloem sieve tube system: a role in mediating traits important to modern agriculture

The plant vascular system serves a vital function by distributing water, nutrients and hormones essential for growth and development to the various organs of the plant. In this review, attention is focused on the role played by the phloem as the conduit for delivery of both photosynthate and information macromolecules, especially from the context of its mediation in traits that are important to modern agriculture. Resource allocation of sugars and amino acids, by the phloem, to specific sink tissues is of importance to crop yield and global food security. Current findings are discussed in the context of a hierarchical control network that operates to integrate resource allocation to competing sinks. The role of plasmodesmata that connect companion cells to neighbouring sieve elements and phloem parenchyma cells is evaluated in terms of their function as valves, connecting the sieve tube pressure manifold system to the various plant tissues. Recent studies have also revealed that plasmodesmata and the phloem sieve tube system function cooperatively to mediate the long-distance delivery of proteins and a diverse array of RNA species. Delivery of these information macromolecules is discussed in terms of their roles in control over the vegetative-to-floral transition, tuberization in potato, stress-related signalling involving miRNAs, and genetic reprogramming through the delivery of 24-nucleotide small RNAs that function in transcriptional gene silencing in recipient sink organs. Finally, we discuss important future research areas that could contribute to developing agricultural crops with engineered performance characteristics for enhance yield potential.

RNA-seq analysis identifies an intricate regulatory network controlling cluster root development in white lupin

Background

Highly adapted plant species are able to alter their root architecture to improve nutrient uptake and thrive in environments with limited nutrient supply. Cluster roots (CRs) are specialised structures of dense lateral roots formed by several plant species for the effective mining of nutrient rich soil patches through a combination of increased surface area and exudation of carboxylates. White lupin is becoming a model-species allowing for the discovery of gene networks involved in CR development. A greater understanding of the underlying molecular mechanisms driving these developmental processes is important for the generation of smarter plants for a world with diminishing resources to improve food security.

Results

RNA-seq analyses for three developmental stages of the CR formed under phosphorus-limited conditions and two of non-cluster roots have been performed for white lupin. In total 133,045,174 high-quality paired-end reads were used for a de novo assembly of the root transcriptome and merged with LAGI01 (Lupinus albus gene index) to generate an improved LAGI02 with 65,097 functionally annotated contigs. This was followed by comparative gene expression analysis. We show marked differences in the transcriptional response across the various cluster root stages to adjust to phosphate limitation by increasing uptake capacity and adjusting metabolic pathways. Several transcription factors such as PLT, SCR, PHB, PHV or AUX/IAA with a known role in the control of meristem activity and developmental processes show an increased expression in the tip of the CR. Genes involved in hormonal responses (PIN, LAX, YUC) and cell cycle control (CYCA/B, CDK) are also differentially expressed. In addition, we identify primary transcripts of miRNAs with established function in the root meristem.

Conclusions

Our gene expression analysis shows an intricate network of transcription factors and plant hormones controlling CR initiation and formation. In addition, functional differences between the different CR developmental stages in the acclimation to phosphorus starvation have been identified.

Ethylene Response Factor070 regulates root development and phosphate starvation-mediated responses

Inorganic phosphate (Pi) availability is a major factor determining growth and consequently the productivity of crops. However, it is one of the least available macronutrients due to its high fixation in the rhizospheres. To overcome this constraint, plants have developed adaptive responses to better acquire, utilize, and recycle Pi. Molecular determinants of these adaptive mechanisms include transcription factors (TFs) that play a major role in transcriptional control, thereby regulating genome-scale networks. In this study, we have characterized the biological role of Arabidopsis thaliana Ethylene Response Factor070 (AtERF070), a Pi starvation-induced TF belonging to the Apetala2/Ethylene Response Factor family of TFs in Arabidopsis (Arabidopsis thaliana). It is localized to the nucleus and induced specifically in Pi-deprived roots and shoots. RNA interference-mediated suppression of AtERF070 led to augmented lateral root development resulting in higher Pi accumulation, whereas there were reductions in both primary root length and lateral root number in 12-d-old transgenic seedlings overexpressing AtERF070. When the overexpressing lines were grown to maturity under greenhouse conditions, they revealed a stunted bushy appearance that could be rescued by gibberellic acid application. Furthermore, a number of Pi starvation-responsive genes were modulated in AtERF070-overexpressing and RNA interference lines, thereby suggesting a potential role for this TF in maintaining Pi homeostasis.

Molecular mechanisms underlying phosphate sensing, signaling, and adaptation in plants

As an essential plant macronutrient, the low availability of phosphorus (P) in most soils imposes serious limitation on crop production. Plants have evolved complex responsive and adaptive mechanisms for acquisition, remobilization and recycling of phosphate (Pi) to maintain P homeostasis. Spatio-temporal molecular, physiological, and biochemical Pi deficiency responses developed by plants are the consequence of local and systemic sensing and signaling pathways. Pi deficiency is sensed locally by the root system where hormones serve as important signaling components in terms of developmental reprogramming, leading to changes in root system architecture. Root-to-shoot and shoot-to-root signals, delivered through the xylem and phloem, respectively, involving Pi itself, hormones, miRNAs, mRNAs, and sucrose, serve to coordinate Pi deficiency responses at the whole-plant level. A combination of chromatin remodeling, transcriptional and posttranslational events contribute to globally regulating a wide range of Pi deficiency responses. In this review, recent advances are evaluated in terms of progress toward developing a comprehensive understanding of the molecular events underlying control over P homeostasis. Application of this knowledge, in terms of developing crop plants having enhanced attributes for P use efficiency, is discussed from the perspective of agricultural sustainability in the face of diminishing global P supplies.

RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis

BACKGROUND

One of the concerns of assembling de novo transcriptomes is determining the amount of read sequences required to ensure a comprehensive coverage of genes expressed in a particular sample. In this report, we describe the use of Illumina paired-end RNA-Seq (PE RNA-Seq) reads from Hevea brasiliensis (rubber tree) bark to devise a transcript mapping approach for the estimation of the read amount needed for deep transcriptome coverage.

FINDINGS

We optimized the assembly of a Hevea bark transcriptome based on 16 Gb Illumina PE RNA-Seq reads using the Oases assembler across a range of k-mer sizes. We then assessed assembly quality based on transcript N50 length and transcript mapping statistics in relation to (a) known Hevea cDNAs with complete open reading frames, (b) a set of core eukaryotic genes and (c) Hevea genome scaffolds. This was followed by a systematic transcript mapping process where sub-assemblies from a series of incremental amounts of bark transcripts were aligned to transcripts from the entire bark transcriptome assembly. The exercise served to relate read amounts to the degree of transcript mapping level, the latter being an indicator of the coverage of gene transcripts expressed in the sample. As read amounts or datasize increased toward 16 Gb, the number of transcripts mapped to the entire bark assembly approached saturation. A colour matrix was subsequently generated to illustrate sequencing depth requirement in relation to the degree of coverage of total sample transcripts.

CONCLUSIONS

We devised a procedure, the "transcript mapping saturation test", to estimate the amount of RNA-Seq reads needed for deep coverage of transcriptomes. For Hevea de novo assembly, we propose generating between 5-8 Gb reads, whereby around 90% transcript coverage could be achieved with optimized k-mers and transcript N50 length. The principle behind this methodology may also be applied to other non-model plants, or with reads from other second generation sequencing platforms.

Phosphate nutrition: improving low-phosphate tolerance in crops

Phosphorus is an essential nutrient that is required for all major developmental processes and reproduction in plants. It is also a major constituent of the fertilizers required to sustain high-yield agriculture. Levels of phosphate--the only form of phosphorus that can be assimilated by plants--are suboptimal in most natural and agricultural ecosystems, and when phosphate is applied as fertilizer in soils, it is rapidly immobilized owing to fixation and microbial activity. Thus, cultivated plants use only approximately 20-30% of the applied phosphate, and the rest is lost, eventually causing water eutrophication. Recent advances in the understanding of mechanisms by which wild and cultivated species adapt to low-phosphate stress and the implementation of alternative bacterial pathways for phosphorus metabolism have started to allow the design of more effective breeding and genetic engineering strategies to produce highly phosphate-efficient crops, optimize fertilizer use, and reach agricultural sustainability with a lower environmental cost. In this review, we outline the current advances in research on the complex network of plant responses to low-phosphorus stress and discuss some strategies used to manipulate genes involved in phosphate uptake, remobilization, and metabolism to develop low-phosphate-tolerant crops, which could help in designing more efficient crops.

Trait-directed de novo population transcriptome dissects genetic regulation of a balanced polymorphism in phosphorus nutrition/arsenate tolerance in a wild grass, Holcus lanatus

The aim of this study was to characterize the transcriptome of a balanced polymorphism, under the regulation of a single gene, for phosphate fertilizer responsiveness/arsenate tolerance in wild grass Holcus lanatus genotypes screened from the same habitat. De novo transcriptome sequencing, RNAseq (RNA sequencing) and single nucleotide polymorphism (SNP) calling were conducted on RNA extracted from H. lanatus. Roche 454 sequencing data were assembled into c. 22,000 isotigs, and paired-end Illumina reads for phosphorus-starved (P-) and phosphorus-treated (P+) genovars of tolerant (T) and nontolerant (N) phenotypes were mapped to this reference transcriptome. Heatmaps of the gene expression data showed strong clustering of each P+/P- treated genovar, as well as clustering by N/T phenotype. Statistical analysis identified 87 isotigs to be significantly differentially expressed between N and T phenotypes and 258 between P+ and P- treated plants. SNPs and transcript expression that systematically differed between N and T phenotypes had regulatory function, namely proteases, kinases and ribonuclear RNA-binding protein and transposable elements. A single gene for arsenate tolerance led to distinct phenotype transcriptomes and SNP profiles, with large differences in upstream post-translational and post-transcriptional regulatory genes rather than in genes directly involved in P nutrition transport and metabolism per se.

Functional genomics to study stress responses in crop legumes: progress and prospects

Legumes are important food crops worldwide, contributing to more than 33% of human dietary protein. The production of crop legumes is frequently impacted by abiotic and biotic stresses. It is therefore important to identify genes conferring resistance to biotic stresses and tolerance to abiotic stresses that can be used to both understand molecular mechanisms of plant response to the environment and to accelerate crop improvement. Recent advances in genomics offer a range of approaches such as the sequencing of genomes and transcriptomes, gene expression microarray as well as RNA-seq based gene expression profiling, and map-based cloning for the identification and isolation of biotic and abiotic stress-responsive genes in several crop legumes. These candidate stress associated genes should provide insights into the molecular mechanisms of stress tolerance and ultimately help to develop legume varieties with improved stress tolerance and productivity under adverse conditions. This review provides an overview on recent advances in the functional genomics of crop legumes that includes the discovery as well as validation of candidate genes.

Comparative Analysis of the Symbiotic Efficiency of Medicago truncatula and Medicago sativa under Phosphorus Deficiency

Phosphorus (P)-deficiency is a major abiotic stress that limits legume growth in many types of soils. The relationship between Medicago and Sinorhizobium, is known to be affected by different environmental conditions. Recent reports have shown that, in combination with S. meliloti 2011, Medicago truncatula had a lower symbiotic efficiency than Medicago sativa. However, little is known about how Medicago-Sinorhizobium is affected by P-deficiency at the whole-plant level. The objective of the present study was to compare and characterize the symbiotic efficiency of N2 fixation of M. truncatula and M. sativa grown in sand under P-limitation. Under this condition, M. truncatula exhibited a significantly higher rate of N2 fixation. The specific activity of the nodules was much higher in M. truncatula in comparison to M. sativa, partially as a result of an increase in electron allocation to N2 versus H+. Although the main organic acid, succinate, exhibited a strong tendency to decrease under P-deficiency, the more efficient symbiotic ability observed in M. truncatula coincided with an apparent increase in the content of malate in its nodules. Our results indicate that the higher efficiency of the M. truncatula symbiotic system is related to the ability to increase malate content under limited P-conditions.