The NAC transcription factor family in maritime pine (Pinus Pinaster): molecular regulation of two genes involved in stress responses

Transcriptome profiles reveal response mechanisms and key role of PsNAC1 in Pinus sylvestris var. mongolica to drought stress

BACKGROUND

Drought stress severely impedes plant growth, and only a limited number of species exhibit long-term resistance to such conditions. Pinus sylvestris var. mongolica, a dominant tree species in arid and semi-arid regions of China, exhibits strong drought resistance and plays a crucial role in the local ecosystem. However, the molecular mechanisms underlying this resistance remain poorly understood.

RESULTS

Here, we conducted transcriptome sequence and physiological indicators analysis of needle samples during drought treatment and rehydration stages. De-novo assembly yielded approximately 114,152 unigenes with an N50 length of 1,363 bp. We identified 6,506 differentially expressed genes (DEGs), with the majority being concentrated in the heavy drought stage (4,529 DEGs). Functional annotation revealed enrichment of drought-related GO terms such as response to water (GO:0009415: enriched 108 genes) and response to water deprivation (GO:0009414: enriched 106 genes), as well as KEGG categories including MAPK signaling pathway (K04733: enriched 35 genes) and monoterpenoid biosynthesis (K21374: enriched 27 genes). Multiple transcription factor families and functional protein families were differentially expressed during drought treatment. Co-expression network analysis identified a potential drought regulatory network between cytochrome P450 genes (Unigene4122_c1_g1) and a core regulatory transcription factor Unigene9098_c3_g1 (PsNAC1) with highly significant expression differences. We validated PsNAC1 overexpression in Arabidopsis and demonstrated enhanced drought resistance.

CONCLUSIONS

These findings provide insight into the molecular basis of drought resistance in P. sylvestris var. mongolica and lay the foundation for further exploration of its regulatory network.

Co-Expression of JcNAC1- and JcZFP8-Improved Agronomic Traits of Tobacco and Enhanced Drought Resistance through NbbHLH1 and NbbHLH2

Previous studies have identified numerous transcription factors involved in drought response, each of which play different roles in plants. The objective of the present study was to evaluate the effectiveness of two transcription factors on drought response in Jatropha curcas L., JcNAC1 and JcZFP8. The overexpression of these transcription factors in tobacco (Nicotiana benthamiana L.) improved drought resistance, but JcZFP8 delayed germination and JcNAC1 reduced biomass and yield. By constitutively co-expressing these two genes in tobacco, drought resistance was improved, and the negative effects of each of them were overcome. The transgenic plants with double-gene co-expression showed stronger drought tolerance with 1.76-fold greater accumulation of proline and lower H2O2 and malondialdehyde (MDA) content to 43 and 65% of wildtype (WT) levels, respectively. The expression levels of NbbHLH1 and NbbHLH2 genes upregulated linearly with the increased drought tolerance of double genes co-expression plants. In drought conditions, the leaf water contents of bhlh1, bhlh2, and bhlh1bhlh2 deletion mutants obtained by CRISPR-CAS9 knockout technique were maintained at 99%, 97%, and 97% of WT. The bhlh1bhlh2 was found with lower germination rate but with higher reactive oxygen levels (1.64-fold H2O2 and 1.41-fold MDA levels). Thus, the co-expression of two transcription factors with different functions overcame the adverse traits brought by a single gene and enhanced the shared drought-tolerant traits, which can provide guidance on theory and selection of gene combinations for the application of multi-gene co-expression in agriculture in the future.

An ethylene-induced NAC transcription factor acts as a multiple abiotic stress responsor in conifer

The proper response to various abiotic stresses is essential for plants' survival to overcome their sessile nature, especially for perennial trees with very long-life cycles. However, in conifers, the molecular mechanisms that coordinate multiple abiotic stress responses remain elusive. Here, the transcriptome response to various abiotic stresses like salt, cold, drought, heat shock and osmotic were systematically detected in Pinus tabuliformis (P. tabuliformis) seedlings. We found that four transcription factors were commonly induced by all tested stress treatments, while PtNAC3 and PtZFP30 were highly up-regulated and co-expressed. Unexpectedly, the exogenous hormone treatment assays and the content of the endogenous hormone indicates that the upregulation of PtNAC3 and PtZFP30 are mediated by ethylene. Time-course assay showed that the treatment by ethylene immediate precursor, 1-aminocyclopropane-1-carboxylic acid (ACC), activated the expression of PtNAC3 and PtZFP30 within 8 hours. We further confirm that the PtNAC3 can directly bind to the PtZFP30 promoter region and form a cascade. Overexpression of PtNAC3 enhanced unified abiotic stress tolerance without growth penalty in transgenic Arabidopsis and promoted reproductive success under abiotic stress by shortening the lifespan, suggesting it has great potential as a biological tool applied to plant breeding for abiotic stress tolerance. This study provides novel insights into the hub nodes of the abiotic stresses response network as well as the environmental adaptation mechanism in conifers, and provides a potential biofortification tool to enhance plant unified abiotic stress tolerance.

Overexpression of SmMYC2 enhances salt resistance in Arabidopsis thaliana and Salvia miltiorrhiza hairy roots

Soil salinity significantly affects both Salvia miltiorrhiza growth and development as well as seed germination throughout field cultivation and production. The basic helix-loop-helix (bHLH) transcription factor (TF) MYC2 contributes significantly to plant stress resistance as a key regulator of the jasmonic acid signaling pathway. In transgenic S. miltiorrhiza hairy roots, SmMYC2 has been shown to promote the accumulation of tanshinone and salvianolic acid, but its role in S. miltiorrhiza of resistance to abiotic stress is unclear. Herein, we found methyl jasmonate (MeJA), NaCl, and PEG treatment all significantly increased SmMYC2 expression. In response to salt stress, SmMYC2 overexpression in yeast increased its rate of growth. Additionally, overexpression of SmMYC2 transgenic Arabidopsis thaliana and S. miltiorrhiza hairy root showed that it might improve salt resistance in transgenic plant. In particular, compared to WT, overexpression of SmMYC2 transgenic Arabidopsis had higher levels of three antioxidant enzymes (superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT)), proline (Pro) content, and ABA-dependent and ABA-independent genes expression. They also had lower levels of malondialdehyde (MDA) and reactive oxygen species (ROS) accumulation. What's more, overexpression of SmMYC2 increases the expression of flavonoid synthesis genes and the accumulation of related components in Arabidopsis. These findings imply that SmMYC2 functions as a positive regulator that regulates plant tolerance to salt through ABA-dependent and independent signaling pathways.

MicroRNA-mediated post-transcriptional regulation of Pinus pinaster response and resistance to pinewood nematode

Pine wilt disease (PWD), caused by the parasitic nematode Bursaphelenchus xylophilus, or pinewood nematode (PWN), is a serious threat to pine forests in Europe. Pinus pinaster is highly susceptible to the disease and it is currently the most affected European pine species. In this work, we investigated the role of small RNAs (sRNAs) in regulating P. pinaster–PWN interaction in an early stage of infection. After performing an artificial PWN inoculation assay, we have identified 105 plant microRNAs (miRNAs) responsive to PWN. Based on their predicted targets, part of these miRNAs was associated with roles in jasmonate-response pathway, ROS detoxification, and terpenoid biosynthesis. Furthermore, by comparing resistant and susceptible plants, eight miRNAs with putative functions in plant defence and resistance to PWN have been identified. Finally, we explored the possibility of bidirectional trans-kingdom RNA silencing, identifying several P. pinaster genes putatively targeted by PWN miRNAs, which was supported by degradome analysis. Targets for P. pinaster miRNAs were also predicted in PWN, suggesting a role for trans-kingdom miRNA transfer and gene silencing both in PWN parasitism as in P. pinaster resistance to PWD. Our results provide new insights into previously unexplored roles of sRNA post-transcriptional regulation in P. pinaster response and resistance to PWN.

Transcriptomic and Metabolomic Analyses Reveal Key Metabolites, Pathways and Candidate Genes in Sophora davidii (Franch.) Skeels Seedlings Under Drought Stress

Soil aridification and desertification are particularly prominent in China's karst areas, severely limiting crop yields and vegetation restoration. Therefore, it is very important to identify naturally drought-tolerant plant species. Sophora davidii (Franch.) Skeels is resistant to drought and soil infertility, is deeply rooted and is an excellent plant material for soil and water conservation. We studied the transcriptomic and metabolomic changes in S. davidii in response to drought stress (CK, control; LD, mild drought stress; MD, moderate drought stress; and SD, severe drought stress). Sophora davidii grew normally under LD and MD stress but was inhibited under SD stress; the malondialdehyde (MDA), hydrogen peroxide (H₂O₂), soluble sugar, proline, chlorophyll a, chlorophyll b and carotenoid contents and ascorbate peroxidase (APX) activity significantly increased, while the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities and soluble protein content significantly decreased. In the LD/CK, MD/CK and SD/CK comparison groups, there were 318, 734 and 1779 DEGs, respectively, and 100, 168 and 281 differentially accumulated metabolites, respectively. Combined analysis of the transcriptomic and metabolomic data revealed the metabolic regulation of S. davidii in response to drought stress. First, key candidate genes such as PRR7, PRR5, GI, ELF3, PsbQ, PsaK, INV, AMY, E2.4.1.13, E3.2.1.2, NCED, PP2C, PYL, ABF, WRKY33, P5CS, PRODH, AOC3, HPD, GPX, GST, CAT and SOD1 may govern the drought resistance of S. davidii. Second, three metabolites (oxidised glutathione, abscisic acid and phenylalanine) were found to be related to drought tolerance. Third, several key candidate genes and metabolites involved in 10 metabolic pathways were identified, indicating that these metabolic pathways play an important role in the response to drought in S. davidii and possibly other plant species.

Elucidation of the biochemical and molecular basis of the differential disease expression in two cultivars of chili (Capsicum annuum) in response to Colletotrichum capsici infection

Chili plants are affected by the hemibiotrophic ascomycota fungus Colletotrichum capsici causing Anthracnose. Infection results in yield and marketability loss due to a decrease in the quality of fruits. The study of morphological symptom development in two cultivars, Bullet, and Beldanga, showed very different disease expression pattern. To understand the reasons behind such differential response, we investigated, in a time-dependent manner, biochemical activities of important defense enzymes, PR proteins, like peroxidase, polyphenol-oxidase, phenylalanine ammonia lyase, β-glucanase, chitinase, catalase, as well as phenols, flavonoids, chlorophyll and the key signaling molecule nitric oxide in their leaves. We further performed real-time nitric oxide (NO) detection studies. The results showed striking differences in the activity profile of these defense molecules through the course of the study. We monitored the gene expression levels of 12 important defense-related genes under in vivo condition. The transcription levels were mostly increased in the tolerant cultivar till 7 days post-infection (DPI), while downregulation of some of the genes were observed in the susceptible one. These data indicated that disease manifestation is a simulated response of these defense molecules which can nullify the effect of the pathogen and its products, when resistance occurs. Alternatively, the pathogen suppresses the host defense when the disease develops.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11738-021-03334-x.

Wood transcriptome analysis of Pinus densiflora identifies genes critical for secondary cell wall formation and NAC transcription factors involved in tracheid formation

Although conifers have significant ecological and economic value, information on transcriptional regulation of wood formation in conifers is still limited. Here, to gain insight into secondary cell wall (SCW) biosynthesis and tracheid formation in conifers, we performed wood tissue-specific transcriptome analyses of Pinus densiflora (Korean red pine) using RNA sequencing. In addition, to obtain full-length transcriptome information, PacBio single molecule real-time iso-sequencing was carried out using RNAs from 28 tissues of P. densiflora. Subsequent comparative tissue-specific transcriptome analysis successfully pinpointed critical genes encoding key proteins involved in biosynthesis of the major secondary wall components (cellulose, galactoglucomannan, xylan and lignin). Furthermore, we predicted a total of 62 NAC (NAM, ATAF1/2 and CUC2) family transcription factor members and identified seven PdeNAC genes preferentially expressed in developing xylem tissues in P. densiflora. Protoplast-based transcriptional activation analysis found that four PdeNAC genes, homologous to VND, NST and SND/ANAC075, upregulated GUS activity driven by an SCW-specific cellulose synthase promoter. Consistently, transient overexpression of the four PdeNACs induced xylem vessel cell-like SCW deposition in both tobacco (Nicotiana benthamiana) and Arabidopsis leaves. Taken together, our data provide a foundation for further research to unravel transcriptional regulation of wood formation in conifers, especially SCW formation and tracheid differentiation.

Insights Into the Mechanisms Implicated in Pinus pinaster Resistance to Pinewood Nematode

Pine wilt disease (PWD), caused by the plant-parasitic nematode Bursaphelenchus xylophilus, has become a severe environmental problem in the Iberian Peninsula with devastating effects in Pinus pinaster forests. Despite the high levels of this species' susceptibility, previous studies reported heritable resistance in P. pinaster trees. Understanding the basis of this resistance can be of extreme relevance for future programs aiming at reducing the disease impact on P. pinaster forests. In this study, we highlighted the mechanisms possibly involved in P. pinaster resistance to PWD, by comparing the transcriptional changes between resistant and susceptible plants after infection. Our analysis revealed a higher number of differentially expressed genes (DEGs) in resistant plants (1,916) when compared with susceptible plants (1,226). Resistance to PWN is mediated by the induction of the jasmonic acid (JA) defense pathway, secondary metabolism pathways, lignin synthesis, oxidative stress response genes, and resistance genes. Quantification of the acetyl bromide-soluble lignin confirmed a significant increase of cell wall lignification of stem tissues around the inoculation zone in resistant plants. In addition to less lignified cell walls, susceptibility to the pine wood nematode seems associated with the activation of the salicylic acid (SA) defense pathway at 72 hpi, as revealed by the higher SA levels in the tissues of susceptible plants. Cell wall reinforcement and hormone signaling mechanisms seem therefore essential for a resistance response.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression patterns during somatic embryogenesis in Dimocarpus longan Lour

The NAM, ATAF1/2, and CUC2 form a huge plant-specific gene family of NAC TFs that are involved in the growth, development, and regulation of biotic and abiotic stress responses. Although the draft genome of longan (Dimocarpus longan Lour.) has been published, however the comprehensive data regarding the functions, evolution, and expression patterns of the NAC family are still unavailable. In this study, a comprehensive analysis of the NAC transcription factor family in longan was performed, and a total of 114 NAC genes were found. We investigated the NAC gene family exploring the phylogeny, domain conservation, intron/exon, motifs, cis-regulatory elements, protein-protein interaction, and expression profiles of RNA-seq samples in different tissues and early somatic embryogenesis of longan. Phylogenetic analysis showed that the genes with similar gene structure and motif distribution were clustered in the same group. Cis-element identification indicates the possible role of NAC genes in biological and physiological processes. Protein-protein interaction identified the DlNACs homologous with Arabidopsis proteins. We further investigated the expression pattern of DlNAC genes in different tissues (pulp, stem, large fruit, young fruit, and flower) during somatic embryogenesis at embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC), and globular embryos (GE) stages. The qRT-PCR results showed that the DlNAC genes were expressed higher at EC and GE stage compared with ICpEC stage. In conclusion, our results provide insight into the evolution, diversity, and characterization of NAC genes in the longan and provide a base for understanding their biological roles and molecular mechanisms in plants.

Comprehensive genomic characterization of NAC transcription factor family and their response to salt and drought stress in peanut

BACKGROUND

Peanut is one of the most important oil crop species worldwide. NAC transcription factor (TF) genes play important roles in the salt and drought stress responses of plants by activating or repressing target gene expression. However, little is known about NAC genes in peanut.

RESULTS

We performed a genome-wide characterization of NAC genes from the diploid wild peanut species Arachis duranensis and Arachis ipaensis, which included analyses of chromosomal locations, gene structures, conserved motifs, expression patterns, and cis-acting elements within their promoter regions. In total, 81 and 79 NAC genes were identified from A. duranensis and A. ipaensis genomes. Phylogenetic analysis of peanut NACs along with their Arabidopsis and rice counterparts categorized these proteins into 18 distinct subgroups. Fifty-one orthologous gene pairs were identified, and 46 orthologues were found to be highly syntenic on the chromosomes of both A. duranensis and A. ipaensis. Comparative RNA sequencing (RNA-seq)-based analysis revealed that the expression of 43 NAC genes was up- or downregulated under salt stress and under drought stress. Among these genes, the expression of 17 genes in cultivated peanut (Arachis hypogaea) was up- or downregulated under both stresses. Moreover, quantitative reverse transcription PCR (RT-qPCR)-based analysis revealed that the expression of most of the randomly selected NAC genes tended to be consistent with the comparative RNA-seq results.

CONCLUSION

Our results facilitated the functional characterization of peanut NAC genes, and the genes involved in salt and drought stress responses identified in this study could be potential genes for peanut improvement.

AgNAC1, a celery transcription factor, related to regulation on lignin biosynthesis and salt tolerance

The NAC transcription factor participates in various biotic and abiotic stress responses and plays a critical role in plant development. Lignin is a water-insoluble dietary fiber, but it is second only to cellulose in abundance. Celery is the main source of dietary fiber, but its quality and production are limited by various abiotic stresses. Here, AgNAC1 containing the NAM domain was identified from celery. AgNAC1 was found to be a nuclear protein. Transgenic Arabidopsis thaliana plants hosting AgNAC1 have longer root lengths and stomatal axis lengths than the wide type (WT). The evidence from lignin determination and expression levels of lignin-related genes indicated that AgNAC1 plays a vital role in lignin biosynthesis. Furthermore, the results of the physiological characterization and the drought and salt treatments indicate that AgNAC1-overexpressing plants are significantly resistive to salt stress. Under drought and salt treatments, the AgNAC1 transgenic Arabidopsis thaliana plants presented increased superoxide dismutase (SOD) and peroxidase (POD) activities and decreased malondialdehyde (MDA) content and size of stomatal apertures relatively to the WT plants. The AgNAC1 served as a positive regulator in inducing the expression of stress-responsive genes. Overall, the overexpressing AgNAC1 enhanced the plants' resistance to salt stress and played a regulatory role in lignin accumulation.

Large-Scale Cloning and Comparative Analysis of TaNAC Genes in Response to Stripe Rust and Powdery Mildew in Wheat (Triticum aestivum L.)

The NAM, ATAF1/2, and CUC2 (NAC) transcription factors (TFs) constitute the largest plant-specific TF superfamily, and play important roles in various physiological processes, including stress responses. Stripe rust and powdery mildew are the most damaging of the fungal diseases that afflict wheat (Triticum aestivum L.). However, studies on Triticum aestivum NAC (TaNAC)s’ role in resistance to the two diseases are still limited, especially in an overall comparative analysis of TaNACs responding or not to fungal stress. In the present study, 186 TaNAC transcripts were obtained from the resistant hexaploid wheat line N9134 under fungal stress, and 180 new transcripts were submitted to GenBank. Statistical results show that 35.1% (54/154) of TaNAC genes responded to stripe rust and powdery mildew in the seedling stage. “Abnormal” coding transcripts of differentially expressed (DE)-TaNAC genes in wheat responding to fungal stress were found in a significantly higher proportion (24/117 vs. 8/69, p = 0.0098) than in non-DE-NACs. This hinted that the alternative splicing of TaNAC genes was active in transcriptional or post-transcriptional regulation during plant-pathogen interactions. Full-length NAC proteins were classified into nine groups via phylogenetic analysis. Multiple-sequence alignment revealed diversity in the C-terminal structural organization, but the differentially expressed gene (DEG)-encoding proteins enriched in Subgroups VI and VII were conserved, with WV[L/V]CR amino acid residues in Motif 7 following the NAM domain. Our data that showed TaNAC TFs responded to fungal disease, which was affected by expression levels and by the regulation of multifarious transcript variants. These data for TaNAC responses to stripe rust and/or powdery mildew and their numerous structural variants provide a good resource for NAC function–mechanism analysis in the context of biotic-stress tolerance in wheat.

Stress induced dynamic adjustment of conserved miR164:NAC module

Aims including the rationale

Salinity and drought are the two major stresses limiting the productivity of economically important crops such as Glycine max (soybean). The incidence of these stresses during the pod development stages affects the quality and quantity of seeds, which compromise the yield of soybean. The miR164:NAC module has been shown to play a critical role in regulating the response to salt and drought stress in several plant species. However, biological role of miR164:NAC module in salt stress in soybean is not fully understood.

Methods

In this study, we identified 215 salt responsive miRNAs, using miScript miRNA array with a sensitive and a tolerant soybean genotype, William82 and INCASoy36, respectively. The targets of these salt regulated miRNAs were searched in the degradome datasets.

Key results

It was found that four salt stress deregulated miRNAs targeted the NAC transcription factor and among these miR164k and miR408d showed antagonistic expression in the two soybean genotypes. The expression of miR164k was higher in salt tolerant INCASoy36 as compared to salt sensitive William82, under unstressed conditions. However under salt stress, miR164k was downregulated in INCASoy36 (-2.65 fold), whereas it was upregulated in William82 (4.68 fold). A transient co-expression assay validated that gma-miR164k directs the cleavage of GmNAC1 transcript. Bioinformatics analysis revealed that the regulation of NAC transcription factor family by members of miR164 family is conserved across many species. The dynamic expression profiles of miR164 and NAC-TFs were captured in different tissues of rice, tobacco, and two soybean genotypes under drought and salt stress conditions.

Main conclusion

Collectively, our results suggest that genetically determined dynamic modulation of the conserved miR164:NAC-TF module may play an important role in determining the adaptive response of plants to stress.

Measuring viability selection from prospective cohort mortality studies: A case study in maritime pine

By changing the genetic background available for selection at subsequent life stages, stage-specific selection can define adaptive potential across the life cycle. We propose and evaluate here a neutrality test and a Bayesian method to infer stage-specific viability selection coefficients using sequential random genotypic samples drawn from a longitudinal cohort mortality study, within a generation. The approach is suitable for investigating selective mortality in large natural or experimental cohorts of any organism in which individual tagging and tracking are unfeasible. Numerical simulation results indicate that the method can discriminate loci under strong viability selection, and provided samples are large, yield accurate estimates of the corresponding selection coefficients. Genotypic frequency changes are largely driven by sampling noise under weak selection, however, compromising inference in that case. We apply the proposed methods to analyze viability selection operating at early recruitment stages in a natural maritime pine (Pinus pinaster Ait.) population. We measured temporal genotypic frequency changes at 384 candidate-gene SNP loci among seedlings sampled from the time of emergence in autumn until the summer of the following year, a period with high elimination rates. We detected five loci undergoing allele frequency changes larger than expected from stochastic mortality and sampling, with putative functions that could influence survival at early seedling stages. Our results illustrate how new statistical and sampling schemes can be used to conduct genomic scans of contemporary selection on specific life stages.

Evaluation of the diversity and phylogenetic implications of NAC transcription factor members of four reference species from the different embryophytic plant groups

NAC transcription factors (TFs) are one of the largest and important TF family that are involved in the regulation of plant growth and development. They are characterized by a highly conserved N-terminal domain and a variable C-terminal domain. In the present study, the amino acid sequences of NAC TFs from four embryophytic plant species viz. Arabidopsis thaliana (Angiosperm), Picea abies (Gymnosperm), Selaginella moellendorffii (Pteridophyte) and Physcomitrella patens (Bryophyte) as reference of the different plant groups were collected from the Plant Transcription Factor Database (PTFD) and the phylogenetic relationships were evaluated. The phylogenetic tree revealed that the majority of the NAC members were interspersed in the major subgroups that indicated the expansion of the NAC members predates the speciation events. Thirty one (31), five (05), one (1) and ten (10) paralog pairs were determined respectively for Arabidopsis, Picea, Selaginella and Physcomitrella. The structure-function relationship of paralog pairs were inferred from the phylogenetic tree of combined set of paralogous gene pairs by studying the prevalence of flanking regions and motif analysis of the NAC proteins. The motif analysis revealed the presence of an N-terminal conserved domain, a characteristic of the majority of NAC family proteins. Conserved motifs in the C-terminal region were absent in the majority of the protein sequences except few members in Arabidopsis and Physcomitrella. Also the time of gene duplication of the paralog pairs were calculated that revealed the duplication events occurred between 4.48 and 45.94 MYA Arabidopsis, 167.57-532.86 MYA in Picea, and 29.12-53.53 MYA in Physcomitrella.

Cytological, Biochemical and Molecular Events of the Embryogenic State in Douglas-fir (Pseudotsuga menziesii [Mirb.])

Somatic embryogenesis techniques have been developed for most coniferous species, but only using very juvenile material. To extend the techniques' scope, better integrated understanding of the key biological, physiological and molecular characteristics of embryogenic state is required. Therefore, embryonal masses (EMs) and non-embryogenic calli (NECs) have been compared during proliferation at multiple levels. EMs and NECs originating from a single somatic embryo (isogenic lines) of each of three unrelated genotypes were used in the analyses, which included comparison of the lines' anatomy by transmission light microscopy, transcriptomes by RNAseq Illumina sequencing, proteomes by free-gel analysis, contents of endogenous phytohormones (indole-3-acetic acid, cytokinins and ABA) by LC-MS analysis, and soluble sugar contents by HPLC. EMs were characterized by upregulation (relative to levels in NECs) of transcripts, proteins, transcription factors and active cytokinins associated with cell differentiation accompanied by histological, carbohydrate content and genetic markers of cell division. In contrast, NECs were characterized by upregulation (relative to levels in EMs) of transcripts, proteins and products associated with responses to stimuli (ABA, degradation forms of cytokinins, phenols), oxidative stress (reactive oxygen species) and carbohydrate storage (starch). Sub-Network Enrichment Analyses that highlighted functions and interactions of transcripts and proteins that significantly differed between EMs and NECs corroborated these findings. The study shows the utility of a novel approach involving integrated multi-scale transcriptomic, proteomic, biochemical, histological and anatomical analyses to obtain insights into molecular events associated with embryogenesis and more specifically to the embryogenic state of cell in Douglas-fir.

Unique resistance traits against downy mildew from the center of origin of grapevine (Vitis vinifera)

The Eurasian grapevine (Vitis vinifera), an Old World species now cultivated worldwide for high-quality wine production, is extremely susceptible to the agent of downy mildew, Plasmopara viticola. The cultivation of resistant V. vinifera varieties would be a sustainable way to reduce the damage caused by the pathogen and the impact of disease management, which involves the economic, health and environmental costs of frequent fungicide application. We report the finding of unique downy mildew resistance traits in a winemaking cultivar from the domestication center of V. vinifera, and characterize the expression of a range of genes associated with the resistance mechanism. Based on comparative experimental inoculations, confocal microscopy and transcriptomics analyses, our study shows that V. vinifera cv. Mgaloblishvili, native to Georgia (South Caucasus), exhibits unique resistance traits against P. viticola. Its defense response, leading to a limitation of P. viticola growth and sporulation, is determined by the overexpression of genes related to pathogen recognition, the ethylene signaling pathway, synthesis of antimicrobial compounds and enzymes, and the development of structural barriers. The unique resistant traits found in Mgaloblishvili highlight the presence of a rare defense system in V. vinifera against P. viticola which promises fresh opportunities for grapevine genetic improvement.

PpNAC1, a main regulator of phenylalanine biosynthesis and utilization in maritime pine

The transcriptional regulation of phenylalanine metabolism is particularly important in conifers, long-lived species that use large amounts of carbon in wood. Here, we show that the Pinus pinaster transcription factor, PpNAC1, is a main regulator of phenylalanine biosynthesis and utilization. A phylogenetic analysis classified PpNAC1 in the NST proteins group and was selected for functional characterization. PpNAC1 is predominantly expressed in the secondary xylem and compression wood of adult trees. Silencing of PpNAC1 in P. pinaster results in the alteration of stem vascular radial patterning and the down-regulation of several genes associated with cell wall biogenesis and secondary metabolism. Furthermore, transactivation and EMSA analyses showed that PpNAC1 is able to activate its own expression and PpMyb4 promoter, while PpMyb4 is able to activate PpMyb8, a transcriptional regulator of phenylalanine and lignin biosynthesis in maritime pine. Together, these results suggest that PpNAC1 is a functional ortholog of the ArabidopsisSND1 and NST1 genes and support the idea that key regulators governing secondary cell wall formation could be conserved between gymnosperms and angiosperms. Understanding the molecular switches controlling wood formation is of paramount importance for fundamental tree biology and paves the way for applications in conifer biotechnology.

Evolution of Hormone Signaling Networks in Plant Defense

Studies with model plants such as Arabidopsis thaliana have revealed that phytohormones are central regulators of plant defense. The intricate network of phytohormone signaling pathways enables plants to activate appropriate and effective defense responses against pathogens as well as to balance defense and growth. The timing of the evolution of most phytohormone signaling pathways seems to coincide with the colonization of land, a likely requirement for plant adaptations to the more variable terrestrial environments, which included the presence of pathogens. In this review, we explore the evolution of defense hormone signaling networks by combining the model plant-based knowledge about molecular components mediating phytohormone signaling and cross talk with available genome information of other plant species. We highlight conserved hubs in hormone cross talk and discuss evolutionary advantages of defense hormone cross talk. Finally, we examine possibilities of engineering hormone cross talk for improvement of plant fitness and crop production.

Sugarcane transcriptome analysis in response to infection caused by Acidovorax avenae subsp. avenae

Sugarcane is an important tropical crop mainly cultivated to produce ethanol and sugar. Crop productivity is negatively affected by Acidovorax avenae subsp avenae (Aaa), which causes the red stripe disease. Little is known about the molecular mechanisms triggered in response to the infection. We have investigated the molecular mechanism activated in sugarcane using a RNA-seq approach. We have produced a de novo transcriptome assembly (TR7) from sugarcane RNA-seq libraries submitted to drought and infection with Aaa. Together, these libraries present 247 million of raw reads and resulted in 168,767 reference transcripts. Mapping in TR7 of reads obtained from infected libraries, revealed 798 differentially expressed transcripts, of which 723 were annotated, corresponding to 467 genes. GO and KEGG enrichment analysis showed that several metabolic pathways, such as code for proteins response to stress, metabolism of carbohydrates, processes of transcription and translation of proteins, amino acid metabolism and biosynthesis of secondary metabolites were significantly regulated in sugarcane. Differential analysis revealed that genes in the biosynthetic pathways of ET and JA PRRs, oxidative burst genes, NBS-LRR genes, cell wall fortification genes, SAR induced genes and pathogenesis-related genes (PR) were upregulated. In addition, 20 genes were validated by RT-qPCR. Together, these data contribute to a better understanding of the molecular mechanisms triggered by the Aaa in sugarcane and opens the opportunity for the development of molecular markers associated with disease tolerance in breeding programs.

Transcript profiling for early stages during embryo development in Scots pine

BACKGROUND

Characterization of the expression and function of genes regulating embryo development in conifers is interesting from an evolutionary point of view. However, our knowledge about the regulation of embryo development in conifers is limited. During early embryo development in Pinus species the proembyo goes through a cleavage process, named cleavage polyembryony, giving rise to four embryos. One of these embryos develops to a dominant embryo, which will develop further into a mature, cotyledonary embryo, while the other embryos, the subordinate embryos, are degraded. The main goal of this study has been to identify processes that might be important for regulating the cleavage process and for the development of a dominant embryo.

RESULTS

RNA samples from embryos and megagametophytes at four early developmental stages during seed development in Pinus sylvestris were subjected to high-throughput sequencing. A total of 6.6 million raw reads was generated, resulting in 121,938 transcripts, out of which 36.106 contained ORFs. 18,638 transcripts were differentially expressed (DETs) in embryos and megagametophytes. GO enrichment analysis of transcripts up-regulated in embryos showed enrichment for different cellular processes, while those up-regulated in megagametophytes were enriched for accumulation of storage material and responses to stress. The highest number of DETs was detected during the initiation of the cleavage process. Transcripts related to embryogenic competence, cell wall modifications, cell division pattern, axis specification and response to hormones and stress were highly abundant and differentially expressed during early embryo development. The abundance of representative DETs was confirmed by qRT-PCR analyses.

CONCLUSION

Based on the processes identified in the GO enrichment analyses and the expression of the selected transcripts we suggest that (i) processes related to embryogenic competence and cell wall loosening are involved in activating the cleavage process; (ii) apical-basal polarization is strictly regulated in dominant embryos but not in the subordinate embryos; (iii) the transition from the morphogenic phase to the maturation phase is not completed in subordinate embryos. This is the first genome-wide transcript expression profiling of the earliest stages during embryo development in a Pinus species. Our results can serve as a framework for future studies to reveal the functions of identified genes.

Transcriptome Analysis of the Response to NaCl in Suaeda maritima Provides an Insight into Salt Tolerance Mechanisms in Halophytes

Although salt tolerance is a feature representative of halophytes, most studies on this topic in plants have been conducted on glycophytes. Transcriptome profiles are also available for only a limited number of halophytes. Hence, the present study was conducted to understand the molecular basis of salt tolerance through the transcriptome profiling of the halophyte Suaeda maritima, which is an emerging plant model for research on salt tolerance. Illumina sequencing revealed 72,588 clustered transcripts, including 27,434 that were annotated using BLASTX. Salt application resulted in the 2-fold or greater upregulation of 647 genes and downregulation of 735 genes. Of these, 391 proteins were homologous to proteins in the COGs (cluster of orthologous groups) database, and the majorities were grouped into the poorly characterized category. Approximately 50% of the genes assigned to MapMan pathways showed homology to S. maritima. The majority of such genes represented transcription factors. Several genes also contributed to cell wall and carbohydrate metabolism, ion relation, redox responses and G protein, phosphoinositide and hormone signaling. Real-time PCR was used to validate the results of the deep sequencing for the most of the genes. This study demonstrates the expression of protein kinase C, the target of diacylglycerol in phosphoinositide signaling, for the first time in plants. This study further reveals that the biochemical and molecular responses occurring at several levels are associated with salt tolerance in S. maritima. At the structural level, adaptations to high salinity levels include the remodeling of cell walls and the modification of membrane lipids. At the cellular level, the accumulation of glycinebetaine and the sequestration and exclusion of Na+ appear to be important. Moreover, this study also shows that the processes related to salt tolerance might be highly complex, as reflected by the salt-induced enhancement of transcription factor expression, including hormone-responsive factors, and that this process might be initially triggered by G protein and phosphoinositide signaling.

Biosynthesis and Metabolic Fate of Phenylalanine in Conifers

The amino acid phenylalanine (Phe) is a critical metabolic node that plays an essential role in the interconnection between primary and secondary metabolism in plants. Phe is used as a protein building block but it is also as a precursor for numerous plant compounds that are crucial for plant reproduction, growth, development, and defense against different types of stresses. The metabolism of Phe plays a central role in the channeling of carbon from photosynthesis to the biosynthesis of phenylpropanoids. The study of this metabolic pathway is particularly relevant in trees, which divert large amounts of carbon into the biosynthesis of Phe-derived compounds, particularly lignin, an important constituent of wood. The trunks of trees are metabolic sinks that consume a considerable percentage of carbon and energy from photosynthesis, and carbon is finally immobilized in wood. This paper reviews recent advances in the biosynthesis and metabolic utilization of Phe in conifer trees. Two alternative routes have been identified: the ancient phenylpyruvate pathway that is present in microorganisms, and the arogenate pathway that possibly evolved later during plant evolution. Additionally, an efficient nitrogen recycling mechanism is required to maintain sustained growth during xylem formation. The relevance of phenylalanine metabolic pathways in wood formation, the biotic interactions, and ultraviolet protection is discussed. The genetic manipulation and transcriptional regulation of the pathways are also outlined.

Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize

NAC proteins are plant-specific transcription factors (TFs). Although they play a pivotal role in regulating distinct biological processes, TFs in maize are yet to be investigated comprehensively. Within the maize genome, we identified 152 putative NAC domain-encoding genes (ZmNACs), including eight membrane-bound members, by systematic sequence analysis and physically mapped them onto ten chromosomes of maize. In silico analysis of the ZmNACs and comparison with similar genes in other plants such as Arabidopsis, rice, and soybean, revealed a similar NAC sequence architecture. Phylogenetically, the ZmNACs were arranged into six distinct subgroups (I-VI) possessing conserved motifs. Phylogenetic analysis using stress-related NAC TFs from Arabidopsis, rice, and soybean as seeding sequences identified 24 of the 152 ZmNACs (all from Group II) as putative stress-responsive genes, including one dehydration-responsive ZmSNAC1 gene reported earlier. One drought-tolerant genotype (HKI577) and one susceptible genotype (PC13T-3) were used for studying the expression pattern of the NAC genes during drought stress. qRT-PCR based expression profiles of 11 genes predicted to be related to stress confirmed strong differential gene expression during drought stress. Phylogenetic analyses revealed that ZmNAC18, ZmNAC51, ZmNAC145, and ZmNAC72, which were up-regulated in the tolerant genotype and down-regulated in the susceptible genotype, belonged to the same group to which also belong other drought-responsive genes, namely SNAC1, OsNAC6, ANAC019, and ANAC055, which act as a transcriptional activator and are strongly induced under stress from various abiotic sources. Differentially expressed ZmNAC genes, alone or in combination with each other or with other type(s) of TFs, may control the general cellular machinery and regulate stress-responsive downstream genes. Alternatively, they may serve as a platform to regulate a broad set of genes, which are subsequently fine-tuned by specific regulators. This genome-wide identification and expression profiling opens new avenues for systematic functional analysis of new members of the NAC gene family, which may be exploited in developing lines that are better adapted to drought.