RNA-Seq reveals genotype-specific molecular responses to water deficit in eucalyptus

The relative impact of parental and current environment on plant transcriptomes depends on type of stress and genotype

Through developmental plasticity, an individual organism integrates influences from its immediate environment with those due to the environment of its parents. While both effects on phenotypes are well documented, their relative impact has been little studied in natural systems, especially at the level of gene expression. We examined this issue in four genotypes of the annual plant Persicaria maculosa by varying two key resources-light and soil moisture-in both generations. Transcriptomic analyses showed that the relative effects of parent and offspring environment on gene expression (i.e. the number of differentially expressed transcripts, DETs) varied both for the two types of resource stress and among genotypes. For light, immediate environment induced more DETs than parental environment for all genotypes, although the precise proportion of parental versus immediate DETs varied among genotypes. By contrast, the relative effect of soil moisture varied dramatically among genotypes, from 8-fold more DETs due to parental than immediate conditions to 10-fold fewer. These findings provide evidence at the transcriptomic level that the relative impacts of parental and immediate environment on the developing organism may depend on the environmental factor and vary strongly among genotypes, providing potential for the interplay of these developmental influences to evolve.

Transcriptional reprogramming during recovery from drought stress in Eucalyptus grandis

The importance of drought as a constraint to agriculture and forestry is increasing with climate change. Genetic improvement of plants' resilience is one of the mitigation strategies to curb this threat. Although recovery from drought stress is important to long-term drought adaptation and has been considered as an indicator of dehydration tolerance in annual crops, this has not been well explored in forest trees. Thus, we aimed to investigate the physiological and transcriptional changes during drought stress and rewatering in Eucalyptus grandis W. Hill ex Maiden. We set up a greenhouse experiment where we imposed drought stress on 2-year-old seedlings and rewatered the recovery group after 17 days of drought. Our measurement of leaf stomatal conductance (gs) showed that, while gs was reduced by drought stress, it fully recovered after 5 days of rewatering. The RNA-seq analysis from stem samples revealed that genes related to known stress responses such as phytohormone and reactive oxygen species signaling were upregulated, while genes involved in metabolism and growth were downregulated due to drought stress. We observed reprogramming of signal transduction pathways and metabolic processes at 1 day of rewatering, indicating a quick response to rewatering. Our results suggest that recovery from drought stress may entail alterations in the jasmonic acid, salicylic acid, ethylene and brassinosteroid signaling pathways. Using co-expression network analysis, we identified hub genes, including the putative orthologs of ABI1, ABF2, ABF3, HAI2, BAM1, GolS2 and SIP1 during drought and CAT2, G6PD1, ADG1 and FD-1 during recovery. Taken together, by highlighting the molecular processes and identifying key genes, this study gives an overview of the mechanisms underlying the response of E. grandis to drought stress and recovery that trees may face repeatedly throughout their long life cycle. This provides a useful reference to the identification and further investigation of signaling pathways and target genes for future tree improvement.

The Eucalyptus grandis chloroplast proteome: Seasonal variations in leaf development

Chloroplast metabolism is very sensitive to environmental fluctuations and is intimately related to plant leaf development. Characterization of the chloroplast proteome dynamics can contribute to a better understanding on plant adaptation to different climate scenarios and leaf development processes. Herein, we carried out a discovery-driven analysis of the Eucalyptus grandis chloroplast proteome during leaf maturation and throughout different seasons of the year. The chloroplast proteome from young leaves differed the most from all assessed samples. Most upregulated proteins identified in mature and young leaves were those related to catabolic-redox signaling and biogenesis processes, respectively. Seasonal dynamics revealed unique proteome features in the fall and spring periods. The most abundant chloroplast protein in humid (wet) seasons (spring and summer) was a small subunit of RuBisCO, while in the dry periods (fall and winter) the proteins that showed the most pronounced accumulation were associated with photo-oxidative damage, Calvin cycle, shikimate pathway, and detoxification. Our investigation of the chloroplast proteome dynamics during leaf development revealed significant alterations in relation to the maturation event. Our findings also suggest that transition seasons induced the most pronounced chloroplast proteome changes over the year. This study contributes to a more comprehensive understanding on the subcellular mechanisms that lead to plant leaf adaptation and ultimately gives more insights into Eucalyptus grandis phenology.

RNA-Seq reveals different responses to drought in Neotropical trees from savannas and seasonally dry forests

BACKGROUND

Water is one of the main limiting factors for plant growth and crop productivity. Plants constantly monitor water availability and can rapidly adjust their metabolism by altering gene expression. This leads to phenotypic plasticity, which aids rapid adaptation to climate changes. Here, we address phenotypic plasticity under drought stress by analyzing differentially expressed genes (DEG) in four phylogenetically related neotropical Bignoniaceae tree species: two from savanna, Handroanthus ochraceus and Tabebuia aurea, and two from seasonally dry tropical forests (SDTF), Handroanthus impetiginosus and Handroanthus serratifolius. To the best of our knowledge, this is the first report of an RNA-Seq study comparing tree species from seasonally dry tropical forest and savanna ecosystems.

RESULTS

Using a completely randomized block design with 4 species × 2 treatments (drought and wet) × 3 blocks (24 plants) and an RNA-seq approach, we detected a higher number of DEGs between treatments for the SDTF species H. serratifolius (3153 up-regulated and 2821 down-regulated under drought) and H. impetiginosus (332 and 207), than for the savanna species. H. ochraceus showed the lowest number of DEGs, with only five up and nine down-regulated genes, while T. aurea exhibited 242 up- and 96 down-regulated genes. The number of shared DEGs among species was not related to habitat of origin or phylogenetic relationship, since both T. aurea and H impetiginosus shared a similar number of DEGs with H. serratifolius. All four species shared a low number of enriched gene ontology (GO) terms and, in general, exhibited different mechanisms of response to water deficit. We also found 175 down-regulated and 255 up-regulated transcription factors from several families, indicating the importance of these master regulators in drought response.

CONCLUSION

Our findings show that phylogenetically related species may respond differently at gene expression level to drought stress. Savanna species seem to be less responsive to drought at the transcriptional level, likely due to morphological and anatomical adaptations to seasonal drought. The species with the largest geographic range and widest edaphic-climatic niche, H. serratifolius, was the most responsive, exhibiting the highest number of DEG and up- and down-regulated transcription factors (TF).

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Molecular bases of responses to abiotic stress in trees

Trees are constantly exposed to climate fluctuations, which vary with both time and geographic location. Environmental changes that are outside of the physiological favorable range usually negatively affect plant performance and trigger responses to abiotic stress. Long-living trees in particular have evolved a wide spectrum of molecular mechanisms to coordinate growth and development under stressful conditions, thus minimizing fitness costs. The ongoing development of techniques directed at quantifying abiotic stress has significantly increased our knowledge of physiological responses in woody plants. However, it is only within recent years that advances in next-generation sequencing and biochemical approaches have enabled us to begin to understand the complexity of the molecular systems that underlie these responses. Here, we review recent progress in our understanding of the molecular bases of drought and temperature stresses in trees, with a focus on functional, transcriptomic, epigenetic, and population genomic studies. In addition, we highlight topics that will contribute to progress in our understanding of the plastic and adaptive responses of woody plants to drought and temperature in a context of global climate change.

Proteomic analyses unraveling water stress response in two Eucalyptus species originating from contrasting environments for aridity

Eucalyptus are widely cultivated in several regions of the world due to their adaptability to different climatic conditions and amenable to tree breeding programs. With changes in environmental conditions pointing to an increase in aridity in many areas of the globe, the demand for genetic materials that adapt to this situation is required. Therefore, the aim of this work was to identify contrasting differences between two Eucalyptus species under water stress through the identification of differentially abundant proteins. For this, total protein extraction was proceeded from leaves of both species maintained at 40 and 80% of field capacity (FC). The 80% FC water regime was considered as the control and the 40% FC, severe water stress. The proteins were separated by 2-DE with subsequent identification of those differentially abundant by liquid nanocromatography coupled to high resolution MS (Q-Exactive). Comparative proteomics allowed to identify four proteins (ATP synthase gamma and alpha, glutamine synthetase and a vacuolar protein) that were more abundant in drought-tolerant species and simultaneously less abundant or unchanged in the drought- sensitive species, an uncharacterized protein found exclusively in plants under drought stress and also 10 proteins (plastid-lipid, ruBisCO activase, ruBisCO, protease ClpA, transketolase, isoflavone reductase, ferredoxin-NADP reductase, malate dehydrogenase, aminobutyrate transaminase and sedoheptulose-1-bisphosphatase) induced exclusively in the drought-tolerant species in response to water stress. These results suggest that such proteins may play a crucial role as potential markers of water stress tolerance through the identification of species-specific proteins, and future targets for genetic engineering.

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

Transcriptome analysis reveals metabolic alteration due to consecutive monoculture and abiotic stress stimuli in Rehamannia glutinosa Libosch

We deeply investigated the mechanism underlying metabolic regulation in response to consecutive monoculture (replanting disease) and different abiotic stresses that unfolded the response mechanism to consecutive monoculture problem through RNA-seq analysis. The consecutive monoculture problem (CMP) resulted of complex environmental stresses mediated by multiple factors. Previous studies have noted that multiple stress factors in consecutive monoculture soils or plants severely limited the interpretation of the critical molecular mechanism, and made a predict that the specifically responding factor was autotoxic allelochemicals. To identify the specifically responding genes, we compared transcriptome changes in roots of Rehamannia glutinosa Libosch using consecutive monoculture, salt, drought, and ferulic acid as stress factors. Comparing with normal growth, 2502, 2672, 2485, and 1956 genes were differentially expressed in R. glutinosa under consecutive monoculture practice, salt, drought, and ferulic acid stress, respectively. In addition, 510 genes were specifically expressed under consecutive monoculture, which were not present under the other stress conditions. Integrating the biological and enrichment analyses of the differentially expressed genes, the result demonstrated that the plants could alter enzyme genes expression to reconstruct the complicated metabolic pathways, which used to tolerate the CMP and abiotic stresses. Furthermore, most of the affected pathway genes were closely related to secondary metabolic processes, and the influence of consecutive monoculture practice on the transcriptome genes expression profile was very similar to the profile under salt stress and then to the profile under drought stress. The outlined schematic diagram unfolded the putative signal regulation mechanism in response to the CMP. Genes that differentially up- or down-regulated under consecutive monoculture practice may play important roles in the CMP or replanting disease in R. glutinosa.

Whole-transcriptome response to water stress in a California endemic oak, Quercus lobata

Reduced water availability during drought can create major stress for many plant species. Within a species, populations with a history of seasonal drought may have evolved the ability to tolerate drought more than those in areas of high precipitation and low seasonality. In this study, we assessed response to water stress in a California oak species, Quercus lobata Née, by measuring changes in gene expression profiles before and after a simulated drought stress treatment through water deprivation of seedlings in a greenhouse setting. Using whole-transcriptome sequencing from nine samples from three collection localities, we identified which genes are involved in response to drought stress and tested the hypothesis that seedlings sampled from climatically different regions of the species range respond to water stress differently. We observed a surprisingly massive transcriptional response to drought: 35,347 of 68,434 contigs (52%) were differentially expressed before versus after drought treatment, of which 18,111 were down-regulated and 17,236 were up-regulated. Genes functionally associated with abiotic stresses and death were enriched among the up-regulated genes, whereas metabolic and cell part-related genes were enriched among the down-regulated. We found 56 contigs that exhibited significantly different expression responses to the drought treatment among the three populations (treatment × population interaction), suggesting that those genes may be involved in local adaptation to drought stress. These genes have stress response (e.g., WRKY DNA-binding protein 51 and HSP20-like chaperones superfamily protein), metabolic (e.g., phosphoglycerate kinase and protein kinase superfamily protein), transport/transfer (e.g., cationic amino acid transporter 7 and K+ transporter) and regulatory functions (e.g., WRKY51 and Homeodomain-like transcriptional regulator). Baseline expression levels of 1310 unique contigs also differed among pairs of populations, and they were enriched for metabolic and cell part-related genes. Out of the large fraction of the transcriptome that was differentially expressed in response to our drought treatment, we identified several novel genes that are candidates for involvement in local adaptation to drought.

Medium term water deficit elicits distinct transcriptome responses in Eucalyptus species of contrasting environmental origin

Background

Climatic and edaphic conditions over geological timescales have generated enormous diversity of adaptive traits and high speciation within the genus Eucalyptus (L. Hér.). Eucalypt species occur from high rainfall to semi-arid zones and from the tropics to latitudes as high as 43°S. Despite several morphological and metabolomic characterizations, little is known regarding gene expression differences that underpin differences in tolerance to environmental change. Using species of contrasting taxonomy, morphology and physiology (E. globulus and E. cladocalyx), this study combines physiological characterizations with ‘second-generation’ sequencing to identify key genes involved in eucalypt responses to medium-term water limitation.

Results

One hundred twenty Million high-quality HiSeq reads were created from 14 tissue samples in plants that had been successfully subjected to a water deficit treatment or a well-watered control. Alignment to the E. grandis genome saw 23,623 genes of which 468 exhibited differential expression (FDR < 0.01) in one or both ecotypes in response to the treatment. Further analysis identified 80 genes that demonstrated a significant species-specific response of which 74 were linked to the ‘dry’ species E. cladocalyx where 23 of these genes were uncharacterised. The majority (approximately 80%) of these differentially expressed genes, were expressed in stem tissue. Key genes that differentiated species responses were linked to photoprotection/redox balance, phytohormone/signalling, primary photosynthesis/cellular metabolism and secondary metabolism based on plant metabolic pathway network analysis.

Conclusion

These results highlight a more definitive response to water deficit by a ‘dry’ climate eucalypt, particularly in stem tissue, identifying key pathways and associated genes that are responsible for the differences between ‘wet’ and ‘dry’ climate eucalypts. This knowledge provides the opportunity to further investigate and understand the mechanisms and genetic variation linked to this important environmental response that will assist with genomic efforts in managing native populations as well as in tree improvement programs under future climate scenarios.

Electronic supplementary material

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

Application of Genomic Technologies to the Breeding of Trees

The recent introduction of next generation sequencing (NGS) technologies represents a major revolution in providing new tools for identifying the genes and/or genomic intervals controlling important traits for selection in breeding programs. In perennial fruit trees with long generation times and large sizes of adult plants, the impact of these techniques is even more important. High-throughput DNA sequencing technologies have provided complete annotated sequences in many important tree species. Most of the high-throughput genotyping platforms described are being used for studies of genetic diversity and population structure. Dissection of complex traits became possible through the availability of genome sequences along with phenotypic variation data, which allow to elucidate the causative genetic differences that give rise to observed phenotypic variation. Association mapping facilitates the association between genetic markers and phenotype in unstructured and complex populations, identifying molecular markers for assisted selection and breeding. Also, genomic data provide in silico identification and characterization of genes and gene families related to important traits, enabling new tools for molecular marker assisted selection in tree breeding. Deep sequencing of transcriptomes is also a powerful tool for the analysis of precise expression levels of each gene in a sample. It consists in quantifying short cDNA reads, obtained by NGS technologies, in order to compare the entire transcriptomes between genotypes and environmental conditions. The miRNAs are non-coding short RNAs involved in the regulation of different physiological processes, which can be identified by high-throughput sequencing of RNA libraries obtained by reverse transcription of purified short RNAs, and by in silico comparison with known miRNAs from other species. All together, NGS techniques and their applications have increased the resources for plant breeding in tree species, closing the former gap of genetic tools between trees and annual species.

Large-scale transcriptional profiling of lignified tissues in Tectona grandis

BACKGROUND

Currently, Tectona grandis is one of the most valuable trees in the world and no transcript dataset related to secondary xylem is available. Considering how important the secondary xylem and sapwood transition from young to mature trees is, little is known about the expression differences between those successional processes and which transcription factors could regulate lignin biosynthesis in this tropical tree. Although MYB transcription factors are one of the largest superfamilies in plants related to secondary metabolism, it has not yet been characterized in teak. These results will open new perspectives for studies of diversity, ecology, breeding and genomic programs aiming to understand deeply the biology of this species.

RESULTS

We present a widely expressed gene catalog for T. grandis using Illumina technology and the de novo assembly. A total of 462,260 transcripts were obtained, with 1,502 and 931 genes differentially expressed for stem and branch secondary xylem, respectively, during age transition. Analysis of stem and branch secondary xylem indicates substantial similarity in gene ontologies including carbohydrate enzymes, response to stress, protein binding, and allowed us to find transcription factors and heat-shock proteins differentially expressed. TgMYB1 displays a MYB domain and a predicted coiled-coil (CC) domain, while TgMYB2, TgMYB3 and TgMYB4 showed R2R3-MYB domain and grouped with MYBs from several gymnosperms and flowering plants. TgMYB1, TgMYB4 and TgCES presented higher expression in mature secondary xylem, in contrast with TgMYB2, TgHsp1, TgHsp2, TgHsp3, and TgBi whose expression is higher in young lignified tissues. TgMYB3 is expressed at lower level in secondary xylem.

CONCLUSIONS

Expression patterns of MYB transcription factors and heat-shock proteins in lignified tissues are dissimilar when tree development was evaluated, obtaining more expression of TgMYB1 and TgMYB4 in lignified tissues of 60-year-old trees, and more expression in TgHsp1, TgHsp2, TgHsp3 and TgBi in stem secondary xylem of 12-year-old trees. We are opening a door for further functional characterization by reverse genetics and marker-assisted selection with those genes. Investigation of some of the key regulators of lignin biosynthesis in teak, however, could be a valuable step towards understanding how rigidity of teak wood and extractives content are different from most other woods. The obtained transcriptome data represents new sequences of T. grandis deposited in public databases, representing an unprecedented opportunity to discover several related-genes associated with secondary xylem such as transcription factors and stress-related genes in a tropical tree.

Utilizing intraspecific variation in phenotypic plasticity to bolster agricultural and forest productivity under climate change

Climate change threatens the ability of agriculture and forestry to meet growing global demands for food, fibre and wood products. Information gathered from genotype-by-environment interactions (G × E), which demonstrate intraspecific variation in phenotypic plasticity (the ability of a genotype to alter its phenotype in response to environmental change), may prove important for bolstering agricultural and forest productivity under climate change. Nonetheless, very few studies have explicitly quantified genotype plasticity-productivity relationships in agriculture or forestry. Here, we conceptualize the importance of intraspecific variation in agricultural and forest species plasticity, and discuss the physiological and genetic factors contributing to intraspecific variation in phenotypic plasticity. Our discussion highlights the need for an integrated understanding of the mechanisms of G × E, more extensive assessments of genotypic responses to climate change under field conditions, and explicit testing of genotype plasticity-productivity relationships. Ultimately, further investigation of intraspecific variation in phenotypic plasticity in agriculture and forestry may prove important for identifying genotypes capable of increasing or sustaining productivity under more extreme climatic conditions.

Transcriptome Profiling of the Potato (Solanum tuberosum L.) Plant under Drought Stress and Water-Stimulus Conditions

Drought stress can seriously affect tuberization, yield and quality of potato plant. However, the precise molecular mechanisms governing potato stolon’s response to drought stress and water supply are not very well understood. In this work, a potato (Solanum tuberosum L.) variant, Ningshu 4, was subjected to severe drought stress treatment (DT) and re-watering treatment (RWT) at tuber bulking stage. Strand-specific cDNA libraries of stolon materials were constructed for paired-end transcriptome sequencing analyses and differentially expressed gene (DEG) examination. In comparison to untreated-control (CT) plants, 3189 and 1797 DEGs were identified in DT and RWT plants and 4154 solely expressed DEGs were screened out from these two comparison groups. Interestingly, 263 genes showed opposite expression patterns in DT and RWT plants. Among them, genes homologous to Protein Phosphatase 2C (PP2C), Aspartic protease in guard cell 1 (ASPG1), auxin-responsive protein, Arabidopsis pseudo response regualtor 2 (APRR2), GA stimulated transcripts in Arabidopsis 6 (GASA6), Calmodulin-like protein 19 (CML19), abscisic acid 8'-hydroxylases and calcium-transporting ATPase, et al. were related with drought-stress and water stimulus response. Sixteen DEGs involved in starch synthesis, accumulation and tuber formation exhibited significantly different expression upon re-watering. In addition, 1630, 1527 and 1596 transcription factor encoding genes were detected in CT, DT and RWT. DEGs of ERF, bHLH, MYB, NAC, WRKY, C2H2, bZIP and HD-ZIP families accounted for 50% in three comparison groups, respectively. Furthermore, characteristics of 565 gene ontology (GO) and 108 Kyoto Encyclopedia of Genes and Genomes pathways (KEGG) were analyzed with the 4154 DEGs. All these results suggest that the drought- and water-stimulus response could be implemented by the regulated expression of metabolic pathway DEGs, and these genes were involved in the endogenous hormone biosynthesis and signal transduction pathways. Our data provide more direct information for future study on the interaction between key genes involved in various metabolic pathways under drought stress in potato.

Towards the identification of the loci of adaptive evolution

1. Establishing the genetic and molecular basis underlying adaptive traits is one of the major goals of evolutionary geneticists in order to understand the connection between genotype and phenotype and elucidate the mechanisms of evolutionary change. Despite considerable effort to address this question, there remain relatively few systems in which the genes shaping adaptations have been identified. 2. Here, we review the experimental tools that have been applied to document the molecular basis underlying evolution in several natural systems, in order to highlight their benefits, limitations and suitability. In most cases, a combination of DNA, RNA and functional methodologies with field experiments will be needed to uncover the genes and mechanisms shaping adaptation in nature.

Development of genetic markers in Eucalyptus species by target enrichment and exome sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.

Molecular and physiological responses to abiotic stress in forest trees and their relevance to tree improvement

Abiotic stresses, such as drought, salinity and cold, are the major environmental stresses that adversely affect tree growth and, thus, forest productivity, and play a major role in determining the geographic distribution of tree species. Tree responses and tolerance to abiotic stress are complex biological processes that are best analyzed at a systems level using genetic, genomic, metabolomic and phenomic approaches. This will expedite the dissection of stress-sensing and signaling networks to further support efficient genetic improvement programs. Enormous genetic diversity for stress tolerance exists within some forest-tree species, and due to advances in sequencing technologies the molecular genetic basis for this diversity has been rapidly unfolding in recent years. In addition, the use of emerging phenotyping technologies extends the suite of traits that can be measured and will provide us with a better understanding of stress tolerance. The elucidation of abiotic stress-tolerance mechanisms will allow for effective pyramiding of multiple tolerances in a single tree through genetic engineering. Here we review recent progress in the dissection of the molecular basis of abiotic stress tolerance in forest trees, with special emphasis on Populus, Pinus, Picea, Eucalyptus and Quercus spp. We also outline practices that will enable the deployment of trees engineered for abiotic stress tolerance to land owners. Finally, recommendations for future work are discussed.

RNA-Seq using two populations reveals genes and alleles controlling wood traits and growth in Eucalyptus nitens

Eucalyptus nitens is a perennial forest tree species grown mainly for kraft pulp production in many parts of the world. Kraft pulp yield (KPY) is a key determinant of plantation profitability and increasing the KPY of trees grown in plantations is a major breeding objective. To speed up the breeding process, molecular markers that can predict KPY are desirable. To achieve this goal, we carried out RNA-Seq studies on trees at extremes of KPY in two different trials to identify genes and alleles whose expression correlated with KPY. KPY is positively correlated with growth measured as diameter at breast height (DBH) in both trials. In total, six RNA bulks from two treatments were sequenced on an Illumina HiSeq platform. At 5% false discovery rate level, 3953 transcripts showed differential expression in the same direction in both trials; 2551 (65%) were down-regulated and 1402 (35%) were up-regulated in low KPY samples. The genes up-regulated in low KPY trees were largely involved in biotic and abiotic stress response reflecting the low growth among low KPY trees. Genes down-regulated in low KPY trees mainly belonged to gene categories involved in wood formation and growth. Differential allelic expression was observed in 2103 SNPs (in 1068 genes) and of these 640 SNPs (30%) occurred in 313 unique genes that were also differentially expressed. These SNPs may represent the cis-acting regulatory variants that influence total gene expression. In addition we also identified 196 genes which had Ka/Ks ratios greater than 1.5, suggesting that these genes are under positive selection. Candidate genes and alleles identified in this study will provide a valuable resource for future association studies aimed at identifying molecular markers for KPY and growth.

Water stress and recovery in the performance of two Eucalyptus globulus clones: physiological and biochemical profiles

Eucalyptus plantations are among the most productive forest stands in Portugal and Spain, being mostly used for pulp production and, more recently, as an energy crop. However, the region's Mediterranean climate, with characteristic severe summer drought, negatively affects eucalypt growth and increases mortality. Although the physiological response to water shortage is well characterized for this species, evidence about the plants' recovery ability remains scarce. In order to assess the physiological and biochemical response of Eucalyptus globulus during the recovery phase, two genotypes (AL-18 and AL-10) were submitted to a 3-week water stress period at two different intensities (18 and 25% of field capacity), followed by 1 week of rewatering. Recovery was assessed 1 day and 1 week after rehydration. Drought reduced height, biomass, water potential, NPQ and gas exchange in both genotypes. Contrarily, the levels of pigments, chlorophyll fluorescence parameters (F(v) /F(m) and (φPSII)), MDA and ABA increased. During recovery, the physiological and biochemical profile of stressed plants showed a similar trend: they experienced reversion of altered traits (MDA, ABA, E, g(s), pigments), while other parameters did not recover ((φPSII), NPQ). Furthermore, an overcompensation of CO(2) assimilation was achieved 1 week after rehydration, which was accompanied by greater growth and re-establishment of oxidative balance. Both genotypes were tolerant to the tested conditions, although clonal differences were found. AL-10 was more productive and showed a more rapid and dynamic response to rehydration (namely in carotenoid content, (φPSII) and NPQ) compared to clone AL-18.

Biochemical and molecular changes associated with heteroxylan biosynthesis in Neolamarckia cadamba (Rubiaceae) during xylogenesis

Wood, derived from plant secondary growth, is a commercially important material. Both cellulose and lignin assembly have been well studied during wood formation (xylogenesis), but heteroxylan biosynthesis is less well defined. Elucidation of the heteroxylan biosynthetic pathway is crucial to understand the mechanism of wood formation. Here, we use Neolamarckia cadamba, a fast-growing tropical tree, as a sample to analyze heteroxylan formation at the biochemical and molecular levels during wood formation. Analysis of the non-cellulosic polysaccharides isolated from N. cadamba stems shows that heteroxylans dominate non-cellulosic polysaccharides and increase with xylogenesis. Microsomes isolated from stems of 1-year-old N. cadamba exhibited UDP-Xyl synthase and xylosyltransferase activities with the highest activity present in the middle and basal stem regions. To further understand the genetic basis of heteroxylan synthesis, RNA sequencing (RNA-seq) was used to generate transcriptomes of N. cadamba during xylogenesis. The RNA-seq results showed that genes related to heteroxylan synthesis had higher expression levels in the middle and basal part of the stem compared to the apical part. Our results describe the heteroxylan distribution and heteroxylan synthesis trait in N. cadamba and give a new example for understanding the mechanism of heteroxylan synthesis in tropical tree species in future.

Assessment of the geographic origins of pinewood nematode isolates via single nucleotide polymorphism in effector genes

The pinewood nematode, Bursaphelenchus xylophilus, is native to North America but it only causes damaging pine wilt disease in those regions of the world where it has been introduced. The accurate detection of the species and its dispersal routes are thus essential to define effective control measures. The main goals of this study were to analyse the genetic diversity among B. xylophilus isolates from different geographic locations and identify single nucleotide polymorphism (SNPs) markers for geographic origin, through a comparative transcriptomic approach. The transcriptomes of seven B. xylophilus isolates, from Continental Portugal (4), China (1), Japan (1) and USA (1), were sequenced in the next generation platform Roche 454. Analysis of effector gene transcripts revealed inter-isolate nucleotide diversity that was validated by Sanger sequencing in the genomic DNA of the seven isolates and eight additional isolates from different geographic locations: Madeira Island (2), China (1), USA (1), Japan (2) and South Korea (2). The analysis identified 136 polymorphic positions in 10 effector transcripts. Pairwise comparison of the 136 SNPs through Neighbor-Joining and the Maximum Likelihood methods and 5-mer frequency analysis with the alignment-independent bilinear multivariate modelling approach correlated the SNPs with the isolates geographic origin. Furthermore, the SNP analysis indicated a closer proximity of the Portuguese isolates to the Korean and Chinese isolates than to the Japanese or American isolates. Each geographic cluster carried exclusive alleles that can be used as SNP markers for B. xylophilus isolate identification.

Analysis and annotation of the hexaploid oat seed transcriptome

BACKGROUND

Next generation sequencing provides new opportunities to explore transcriptomes. However, challenges remain for accurate differentiation of homoeoalleles and paralogs, particularly in polyploid organisms with no supporting genome sequence. In this study, RNA-Seq was employed to generate and characterize the first gene expression atlas for hexaploid oat.

RESULTS

The software packages Trinity and Oases were used to produce a transcript assembly from nearly 134 million 100-bp paired-end reads from developing oat seeds. Based on the quality-parameters employed, Oases assemblies were superior. The Oases 67-kmer assembly, denoted dnOST (de novo Oat Seed Transcriptome), is over 55 million nucleotides in length and the average transcript length is 1,043 nucleotides. The 74.8× sequencing depth was adequate to differentiate a large proportion of putative homoeoalleles and paralogs. To assess the robustness of dnOST, we successfully identified gene transcripts associated with the biosynthetic pathways of three compounds with health-promoting properties (avenanthramides, tocols, β-glucans), and quantified their expression.

CONCLUSIONS

To our knowledge, this study provides the first direct performance comparison between two major assemblers in a polyploid organism. The workflow we developed provides a useful guide for comparable analyses in other organisms. The transcript assembly developed here is a major advance. It expands the number of oat ESTs 3-fold, and constitutes the first comprehensive transcriptome study in oat. This resource will be a useful new tool both for analysis of genes relevant to nutritional enhancement of oat, and for improvement of this crop in general.

The plasticity of the grapevine berry transcriptome

BACKGROUND

Phenotypic plasticity refers to the range of phenotypes a single genotype can express as a function of its environment. These phenotypic variations are attributable to the effect of the environment on the expression and function of genes influencing plastic traits. We investigated phenotypic plasticity in grapevine by comparing the berry transcriptome in a single clone of the vegetatively-propagated common grapevine species Vitis vinifera cultivar Corvina through 3 consecutive growth years cultivated in 11 different vineyards in the Verona area of Italy.

RESULTS

Most of the berry transcriptome clustered by year of growth rather than common environmental conditions or viticulture practices, and transcripts related to secondary metabolism showed high sensitivity towards different climates, as confirmed also by metabolomic data obtained from the same samples. When analyzed in 11 vineyards during 1 growth year, the environmentally-sensitive berry transcriptome comprised 5% of protein-coding genes and 18% of the transcripts modulated during berry development. Plastic genes were particularly enriched in ontology categories such as transcription factors, translation, transport, and secondary metabolism. Specific plastic transcripts were associated with groups of vineyards sharing common viticulture practices or environmental conditions, and plastic transcriptome reprogramming was more intense in the year characterized by extreme weather conditions. We also identified a set of genes that lacked plasticity, showing either constitutive expression or similar modulation in all berries.

CONCLUSIONS

Our data reveal candidate genes potentially responsible for the phenotypic plasticity of grapevine and provide the first step towards the characterization of grapevine transcriptome plasticity under different agricultural systems.

Xylem transcription profiles indicate potential metabolic responses for economically relevant characteristics of Eucalyptus species

Background

Eucalyptus is one of the most important sources of industrial cellulose. Three species of this botanical group are intensively used in breeding programs: E. globulus, E. grandis and E. urophylla. E. globulus is adapted to subtropical/temperate areas and is considered a source of high-quality cellulose; E. grandis grows rapidly and is adapted to tropical/subtropical climates; and E. urophylla, though less productive, is considered a source of genes related to robustness. Wood, or secondary xylem, results from cambium vascular differentiation and is mostly composed of cellulose, lignin and hemicelluloses. In this study, the xylem transcriptomes of the three Eucalyptus species were investigated in order to provide insights on the particularities presented by each of these species.

Results

Data analysis showed that (1) most Eucalyptus genes are expressed in xylem; (2) most genes expressed in species-specific way constitutes genes with unknown functions and are interesting targets for future studies; (3) relevant differences were observed in the phenylpropanoid pathway: E. grandis xylem presents higher expression of genes involved in lignin formation whereas E. urophylla seems to deviates the pathway towards flavonoid formation; (4) stress-related genes are considerably more expressed in E. urophylla, suggesting that these genes may contribute to its robustness.

Conclusions

The comparison of these three transcriptomes indicates the molecular signatures underlying some of their distinct wood characteristics. This information may contribute to the understanding of xylogenesis, thus increasing the potential of genetic engineering approaches aiming at the improvement of Eucalyptus forest plantations productivity.

Reference genes for high-throughput quantitative reverse transcription-PCR analysis of gene expression in organs and tissues of Eucalyptus grown in various environmental conditions

Interest in the genomics of Eucalyptus has skyrocketed thanks to the recent sequencing of the genome of Eucalyptus grandis and to a growing number of large-scale transcriptomic studies. Quantitative reverse transcription-PCR (RT-PCR) is the method of choice for gene expression analysis and can now also be used as a high-throughput method. The selection of appropriate internal controls is becoming of utmost importance to ensure accurate expression results in Eucalyptus. To this end, we selected 21 candidate reference genes and used high-throughput microfluidic dynamic arrays to assess their expression among a large panel of developmental and environmental conditions with a special focus on wood-forming tissues. We analyzed the expression stability of these genes by using three distinct statistical algorithms (geNorm, NormFinder and ΔCt), and used principal component analysis to compare methods and rankings. We showed that the most stable genes identified depended not only on the panel of biological samples considered but also on the statistical method used. We then developed a comprehensive integration of the rankings generated by the three methods and identified the optimal reference genes for 17 distinct experimental sets covering 13 organs and tissues, as well as various developmental and environmental conditions. The expression patterns of Eucalyptus master genes EgMYB1 and EgMYB2 experimentally validated our selection. Our findings provide an important resource for the selection of appropriate reference genes for accurate and reliable normalization of gene expression data in the organs and tissues of Eucalyptus trees grown in a range of conditions including abiotic stresses.

Deep sequencing reveals differences in the transcriptional landscapes of fibers from two cultivated species of cotton

Cotton (Gossypium) fiber is the most prevalent natural product used in the textile industry. The two major cultivated species, G. hirsutum (Gh) and G. barbadense (Gb), are allotetraploids with contrasting fiber quality properties. To better understand the molecular basis for their fiber differences, EST pyrosequencing was used to document the fiber transcriptomes at two key development stages, 10 days post anthesis (dpa), representing the peak of fiber elongation, and 22 dpa, representing the transition to secondary cell wall synthesis. The 617,000 high quality reads (89% of the total 692,000 reads) from 4 libraries were assembled into 46,072 unigenes, comprising 38,297 contigs and 7,775 singletons. Functional annotation of the unigenes together with comparative digital gene expression (DGE) revealed a diverse set of functions and processes that were partly linked to specific fiber stages. Globally, 2,770 contigs (7%) showed differential expression (>2-fold) between 10 and 22 dpa (irrespective of genotype), with 70% more highly expressed at 10 dpa, while 2,248 (6%) were differentially expressed between the genotypes (irrespective of stage). The most significant genes with differential DGE at 10 dpa included expansins and lipid transfer proteins (higher in Gb), while at 22 dpa tubulins, cellulose, and sucrose synthases showed higher expression in Gb. DGE was compared with expression data of 10 dpa-old fibers from Affymetrix microarrays. Among 543 contigs showing differential expression on both platforms, 74% were consistent in being either over-expressed in Gh (242 genes) or in Gb (161 genes). Furthermore, the unigene set served to identify 339 new SSRs and close to 21,000 inter-genotypic SNPs. Subsets of 88 SSRs and 48 SNPs were validated through mapping and added 65 new loci to a RIL genetic map. The new set of fiber ESTs and the gene-based markers complement existing available resources useful in basic and applied research for crop improvement in cotton.

A gene-phenotype network based on genetic variability for drought responses reveals key physiological processes in controlled and natural environments

Identifying the connections between molecular and physiological processes underlying the diversity of drought stress responses in plants is key for basic and applied science. Drought stress response involves a large number of molecular pathways and subsequent physiological processes. Therefore, it constitutes an archetypical systems biology model. We first inferred a gene-phenotype network exploiting differences in drought responses of eight sunflower (Helianthus annuus) genotypes to two drought stress scenarios. Large transcriptomic data were obtained with the sunflower Affymetrix microarray, comprising 32423 probesets, and were associated to nine morpho-physiological traits (integrated transpired water, leaf transpiration rate, osmotic potential, relative water content, leaf mass per area, carbon isotope discrimination, plant height, number of leaves and collar diameter) using sPLS regression. Overall, we could associate the expression patterns of 1263 probesets to six phenotypic traits and identify if correlations were due to treatment, genotype and/or their interaction. We also identified genes whose expression is affected at moderate and/or intense drought stress together with genes whose expression variation could explain phenotypic and drought tolerance variability among our genetic material. We then used the network model to study phenotypic changes in less tractable agronomical conditions, i.e. sunflower hybrids subjected to different watering regimes in field trials. Mapping this new dataset in the gene-phenotype network allowed us to identify genes whose expression was robustly affected by water deprivation in both controlled and field conditions. The enrichment in genes correlated to relative water content and osmotic potential provides evidence of the importance of these traits in agronomical conditions.

Accelerating the domestication of forest trees in a changing world

In light of impending water and arable land shortages, population growth and climate change, it is more important than ever to examine how forest tree domestication can be accelerated to sustainably meet future demands for wood, biomass, paper, fuel and biomaterials. Because of long breeding cycles, tree domestication cannot be rapidly achieved through traditional genetic improvement methods alone. Integrating modern genetic and genomic techniques with conventional breeding will expedite tree domestication. Breeders will only embrace these technologies if they are cost-effective and readily accessible, and forest landowners will only adopt end-products that meet with regulatory approval and public acceptance. All parties involved must work together to achieve these objectives for the benefit of society.

RNA-Seq Assembly - Are We There Yet?

Transcriptomic sequence resources represent invaluable assets for research, in particular for non-model species without a sequenced genome. To date, the Next Generation Sequencing technologies 454/Roche and Illumina have been used to generate transcriptome sequence databases by mRNA-Seq for more than fifty different plant species. While some of the databases were successfully used for downstream applications, such as proteomics, the assembly parameters indicate that the assemblies do not yet accurately reflect the actual plant transcriptomes. Two different assembly strategies have been used, overlap consensus based assemblers for long reads and Eulerian path/de Bruijn graph assembler for short reads. In this review, we discuss the challenges and solutions to the transcriptome assembly problem. A list of quality control parameters and the necessary scripts to produce them are provided.