Transcriptional regulation of secondary growth in Arabidopsis thaliana

Spatiotemporal transcriptome atlas reveals gene regulatory patterns during the organogenesis of the rapid growing bamboo shoots

Bamboo with its remarkable growth rate and economic significance, offers an ideal system to investigate the molecular basis of organogenesis in rapidly growing plants, particular in monocots, where gene regulatory networks governing the maintenance and differentiation of shoot apical and intercalary meristems remain a subject of controversy. We employed both spatial and single-nucleus transcriptome sequencing on 10× platform to precisely dissect the gene functions in various tissues and early developmental stages of bamboo shoots. Our comprehensive analysis reveals distinct cell trajectories during shoot development, uncovering critical genes and pathways involved in procambium differentiation, intercalary meristem formation, and vascular tissue development. Spatial and temporal expression patterns of key regulatory genes, particularly those related to hormone signaling and lipid metabolism, strongly support the hypothesis that intercalary meristem origin from surrounded parenchyma cells. Specific gene expressions in intercalary meristem exhibit regular and dispersed distribution pattern, offering clues for understanding the intricate molecular mechanisms that drive the rapid growth of bamboo shoots. The single-nucleus and spatial transcriptome analysis reveal a comprehensive landscape of gene activity, enhancing the understanding of the molecular architecture of organogenesis and providing valuable resources for future genomic and genetic studies relying on identities of specific cell types.

Dissection of transcriptional events in graft incompatible reactions of "Bearss" lemon (Citrus limon) and "Valencia" sweet orange (C. sinensis) on a novel citrandarin (C. reticulata × Poncirus trifoliata) rootstock

Citrus is commercially propagated via grafting, which ensures trees have consistent fruit traits combined with favorable traits from the rootstock such as soil adaptability, vigor, and resistance to soil pathogens. Graft incompatibility can occur when the scion and rootstock are not able to form a permanent, healthy union. Understanding and preventing graft incompatibility is of great importance in the breeding of new fruit cultivars and in the choice of scion and rootstock by growers. The rootstock US-1283, a citrandarin generated from a cross of "Ninkat" mandarin (Citrus reticulata) and "Gotha Road" #6 trifoliate orange (Poncirus trifoliata), was released after years of field evaluation because of its superior productivity and good fruit quality on "Hamlin" sweet orange (C. sinensis) under Florida's growing conditions. Subsequently, it was observed that trees of "Bearss" lemon (C. limon) and "Valencia" sweet orange (C. sinensis) grafted onto US-1283 exhibited unhealthy growth near the graft union. The incompatibility manifested as stem grooving and necrosis underneath the bark on the rootstock side of the graft. Another citrandarin rootstock, US-812 (C. reticulata "Sunki" × P. trifoliata "Benecke"), is fully graft compatible with the same scions. Transcriptome analysis was performed on the vascular tissues above and below the graft union of US-812 and US-1283 graft combinations with "Bearss" and "Valencia" to identify expression networks associated with incompatibility and help understand the processes and potential causes of incompatibility. Transcriptional reprogramming was stronger in the incompatible rootstock than in the grafted scions. Differentially expressed genes (DEGs) in US-1283, but not the scions, were associated with oxidative stress and plant defense, among others, similar to a pathogen-induced immune response localized to the rootstock; however, no pathogen infection was detected. Therefore, it is hypothesized that this response could have been triggered by signaling miscommunications between rootstock and scion either through (1) unknown molecules from the scion that were perceived as danger signals by the rootstock, (2) missing signals from the scion or missing receptors in the rootstock necessary for the formation of a healthy graft union, (3) the overall perception of the scion by the rootstock as non-self, or (4) a combination of the above.

ACC SYNTHASE4 inhibits gibberellin biosynthesis and FLOWERING LOCUS T expression during citrus flowering

Flowering is an essential process in fruit trees. Flower number and timing have a substantial impact on the yield and maturity of fruit. Ethylene and gibberellin (GA) play vital roles in flowering, but the mechanism of coordinated regulation of flowering in woody plants by GA and ethylene is still unclear. In this study, a lemon (Citrus limon L. Burm) 1-aminocyclopropane-1-carboxylic acid synthase gene (CiACS4) was overexpressed in Nicotiana tabacum and resulted in late flowering and increased flower number. Further transformation of citrus revealed that ethylene and starch content increased, and soluble sugar content decreased in 35S:CiACS4 lemon. Inhibition of CiACS4 in lemon resulted in effects opposite to that of 35S:CiACS4 in transgenic plants. Overexpression of the CiACS4-interacting protein ETHYLENE RESPONSE FACTOR3 (CiERF3) in N. tabacum resulted in delayed flowering and more flowers. Further experiments revealed that the CiACS4-CiERF3 complex can bind the promoters of FLOWERING LOCUS T (CiFT) and GOLDEN2-LIKE (CiFE) and suppress their expression. Moreover, overexpression of CiFE in N. tabacum led to early flowering and decreased flowers, and ethylene, starch, and soluble sugar contents were opposite to those in 35S:CiACS4 transgenic plants. Interestingly, CiFE also bound the promoter of CiFT. Additionally, GA3 and 1-aminocyclopropanecarboxylic acid (ACC) treatments delayed flowering in adult citrus, and treatment with GA and ethylene inhibitors increased flower number. ACC treatment also inhibited the expression of CiFT and CiFE. This study provides a theoretical basis for the application of ethylene to regulate flower number and mitigate the impacts of extreme weather on citrus yield due to delayed flowering.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Comparing modeling methods of genomic prediction for growth traits of a tropical timber species, Shorea macrophylla

Introduction

Shorea macrophylla is a commercially important tropical tree species grown for timber and oil. It is amenable to plantation forestry due to its fast initial growth. Genomic selection (GS) has been used in tree breeding studies to shorten long breeding cycles but has not previously been applied to S. macrophylla.

Methods

To build genomic prediction models for GS, leaves and growth trait data were collected from a half-sib progeny population of S. macrophylla in Sari Bumi Kusuma forest concession, central Kalimantan, Indonesia. 18037 SNP markers were identified in two ddRAD-seq libraries. Genomic prediction models based on these SNPs were then generated for diameter at breast height and total height in the 7th year from planting (D7 and H7).

Results and discussion

These traits were chosen because of their relatively high narrow-sense genomic heritability and because seven years was considered long enough to assess initial growth. Genomic prediction models were built using 6 methods and their derivatives with the full set of identified SNPs and subsets of 48, 96, and 192 SNPs selected based on the results of a genome-wide association study (GWAS). The GBLUP and RKHS methods gave the highest predictive ability for D7 and H7 with the sets of selected SNPs and showed that D7 has an additive genetic architecture while H7 has an epistatic genetic architecture. LightGBM and CNN1D also achieved high predictive abilities for D7 with 48 and 96 selected SNPs, and for H7 with 96 and 192 selected SNPs, showing that gradient boosting decision trees and deep learning can be useful in genomic prediction. Predictive abilities were higher in H7 when smaller number of SNP subsets selected by GWAS p-value was used, However, D7 showed the contrary tendency, which might have originated from the difference in genetic architecture between primary and secondary growth of the species. This study suggests that GS with GWAS-based SNP selection can be used in breeding for non-cultivated tree species to improve initial growth and reduce genotyping costs for next-generation seedlings.

Association Mapping and Expression Analysis of the Genes Involved in the Wood Formation of Poplar

Xylogenesis is a complex and sequential biosynthetic process controlled by polygenes. Deciphering the genetic architecture of this complex quantitative trait could provide valuable information for increasing wood biomass and improving its properties. Here, we performed genomic resequencing of 64 24-year-old trees (64 hybrids of section Aigeiros and their parents) grown in the same field and conducted full-sib family-based association analyses of two growth and six woody traits using GEMMA as a choice of association model selection. We identified 1342 significantly associated single nucleotide polymorphisms (SNPs), 673 located in the region upstream and downstream of 565 protein-encoding genes. The transcriptional regulation network of secondary cell wall (SCW) biosynthesis was further constructed based on the published data of poplar miRNA, transcriptome, and degradome. These provided a certain scientific basis for the in-depth understanding of the mechanism of poplar timber formation and the molecular-assisted breeding in the future.

Integrated QTL Mapping, Meta-Analysis, and RNA-Sequencing Reveal Candidate Genes for Maize Deep-Sowing Tolerance

Synergetic elongation of mesocotyl and coleoptile are crucial in governing maize seedlings emergence, especially for the maize sown in deep soil. Studying the genomic regions controlling maize deep-sowing tolerance would aid the development of new varieties that are resistant to harsh conditions, such as drought and low temperature during seed germination. Using 346 F2:3 maize population families from W64A × K12 cross at three sowing depths, we identified 33 quantitative trait loci (QTLs) for the emergence rate, mesocotyl, coleoptile, and seedling lengths via composite interval mapping (CIM). These loci explained 2.89% to 14.17% of phenotypic variation in a single environment, while 12 of 13 major QTLs were identified at two or more sowing environments. Among those, four major QTLs in Bin 1.09, Bin 4.08, Bin 6.01, and Bin 7.02 supported pleiotropy for multiple deep-sowing tolerant traits. Meta-analysis identified 17 meta-QTLs (MQTLs) based on 130 original QTLs from present and previous studies. RNA-Sequencing of mesocotyl and coleoptile in both parents (W64A and K12) at 3 cm and 20 cm sowing environments identified 50 candidate genes expressed differentially in all major QTLs and MQTLs regions: six involved in the circadian clock, 27 associated with phytohormones biosynthesis and signal transduction, seven controlled lignin biosynthesis, five regulated cell wall organization formation and stabilization, three were responsible for sucrose and starch metabolism, and two in the antioxidant enzyme system. These genes with highly interconnected networks may form a complex molecular mechanism of maize deep-sowing tolerance. Findings of this study will facilitate the construction of molecular modules for deep-sowing tolerance in maize. The major QTLs and MQTLs identified could be used in marker-assisted breeding to develop elite maize varieties.

SbMYB3 transcription factor promotes root-specific flavone biosynthesis in Scutellaria baicalensis

Scutellaria baicalensis Georgi produces abundant root-specific flavones (RSFs), which provide various benefits to human health. We have elucidated the complete biosynthetic pathways of baicalein and wogonin. However, the transcriptional regulation of flavone biosynthesis in S. baicalensis remains unclear. We show that the SbMYB3 transcription factor functions as a transcriptional activator involved in the biosynthesis of RSFs in S. baicalensis. Yeast one-hybrid and transcriptional activation assays showed that SbMYB3 binds to the promoter of flavone synthase II-2 (SbFNSII-2) and enhances its transcription. In S. baicalensis hairy roots, RNAi of SbMYB3 reduced the accumulation of baicalin and wogonoside, and SbMYB3 knockout decreased the biosynthesis of baicalein, baicalin, wogonin, and wogonoside, whereas SbMYB3 overexpression enhanced the contents of baicalein, baicalin, wogonin, and wogonoside. Transcript profiling by qRT-PCR demonstrated that SbMYB3 activates SbFNSII-2 expression directly, thus leading to more abundant accumulation of RSFs. This study provides a potential target for metabolic engineering of RSFs.

Comparative transcriptomics of tropical woody plants supports fast and furious strategy along the leaf economics spectrum in lianas

Lianas, climbing woody plants, influence the structure and function of tropical forests. Climbing traits have evolved multiple times, including ancestral groups such as gymnosperms and pteridophytes, but the genetic basis of the liana strategy is largely unknown. Here, we use a comparative transcriptomic approach for 47 tropical plant species, including ten lianas of diverse taxonomic origins, to identify genes that are consistently expressed or downregulated only in lianas. Our comparative analysis of full-length transcripts enabled the identification of a core interactomic network common to lianas. Sets of transcripts identified from our analysis reveal features related to functional traits pertinent to leaf economics spectrum in lianas, include upregulation of genes controlling epidermal cuticular properties, cell wall remodeling, carbon concentrating mechanism, cell cycle progression, DNA repair and a large suit of downregulated transcription factors and enzymes involved in ABA-mediated stress response as well as lignin and suberin synthesis. All together, these genes are known to be significant in shaping plant morphologies through responses such as gravitropism, phyllotaxy and shade avoidance.

Summary: The full-length fraction of liana transcriptomes mapped on a protein–protein interactome revealed the nature of their convergence through distinct sets of expressed and downregulated genes not observed in free-standing plants.

Arabidopsis SnRK1 negatively regulates phenylpropanoid metabolism via Kelch domain-containing F-box proteins

Phenylpropanoid metabolism represents a substantial metabolic sink for photosynthetically fixed carbon. The evolutionarily conserved Sucrose Non-Fermenting Related Kinase 1 (SnRK1) is a major metabolic sensor that reprograms metabolism upon carbon deprivation. However, it is not clear if and how the SnRK1-mediated sugar signaling pathway controls phenylpropanoid metabolism. Here, we show that Arabidopsis SnRK1 negatively regulates phenylpropanoid biosynthesis via a group of Kelch domain-containing F-box (KFB) proteins that are responsible for the ubiquitination and degradation of phenylalanine ammonia lyase (PAL). Downregulation of AtSnRK1 significantly promoted the accumulation of soluble phenolics and lignin polymers and drastically increased PAL cellular accumulation but only slightly altered its transcription level. Co-expression of SnRK1α with PAL in Nicotiana benthamiana leaves resulted in the severe attenuation of the latter's protein level, but protein interaction assays suggested PAL is not a direct substrate of SnRK1. Furthermore, up or downregulation of AtSnRK1 positively affected KFB^PALs gene expression, and energy starvation upregulated KFB^PAL expression, which partially depends on AtSnRK1. Collectively, our study reveals that SnRK1 negatively regulates phenylpropanoid biosynthesis, and KFB^PALs act as regulatory components of the SnRK1 signaling network, transcriptionally regulated by SnRK1 and subsequently mediate proteasomal degradation of PAL in response to the cellular carbon availability.

Histone methyltransferase ATX1 dynamically regulates fiber secondary cell wall biosynthesis in Arabidopsis inflorescence stem

Secondary wall thickening in the sclerenchyma cells is strictly controlled by a complex network of transcription factors in vascular plants. However, little is known about the epigenetic mechanism regulating secondary wall biosynthesis. In this study, we identified that ARABIDOPSIS HOMOLOG of TRITHORAX1 (ATX1), a H3K4-histone methyltransferase, mediates the regulation of fiber cell wall development in inflorescence stems of Arabidopsis thaliana. Genome-wide analysis revealed that the up-regulation of genes involved in secondary wall formation during stem development is largely coordinated by increasing level of H3K4 tri-methylation. Among all histone methyltransferases for H3K4me3 in Arabidopsis, ATX1 is markedly increased during the inflorescence stem development and loss-of-function mutant atx1 was impaired in secondary wall thickening in interfascicular fibers. Genetic analysis showed that ATX1 positively regulates secondary wall deposition through activating the expression of secondary wall NAC master switch genes, SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN1 (SND1) and NAC SECONDARY WALL THICKENING PROMOTING FACTOR1 (NST1). We further identified that ATX1 directly binds the loci of SND1 and NST1, and activates their expression by increasing H3K4me3 levels at these loci. Taken together, our results reveal that ATX1 plays a key role in the regulation of secondary wall biosynthesis in interfascicular fibers during inflorescence stem development of Arabidopsis.

Microarray analysis of Arabidopsis thaliana exposed to single and mixed infections with Cucumber mosaic virus and turnip viruses

Cucumber mosaic virus (CMV), Turnip mosaic virus (TuMV) and Turnip crinkle virus (TCV) are important plant infecting viruses. In the present study, whole transcriptome alteration of Arabidopsis thaliana in response to CMV, TuMV and TCV, individual as well as mixed infections of CMV and TuMV/CMV and TCV were investigated using microarray data. In response to CMV, TuMV and TCV infections, a total of 2517, 3985 and 277 specific differentially expressed genes (DEGs) were up-regulated, while 2615, 3620 and 243 specific DEGs were down-regulated, respectively. The number of 1222 and 30 common DEGs were up-regulated during CMV and TuMV as well as CMV and TCV infections, while 914 and 24 common DEGs were respectively down-regulated. Genes encoding immune response mediators, signal transducer activity, signaling and stress response functions were among the most significantly upregulated genes during CMV and TuMV or CMV and TCV mixed infections. The NAC, C3H, C2H2, WRKY and bZIP were the most commonly presented transcription factor (TF) families in CMV and TuMV infection, while AP2-EREBP and C3H were the TF families involved in CMV and TCV infections. Moreover, analysis of miRNAs during CMV and TuMV and CMV and TCV infections have demonstrated the role of miRNAs in the down regulation of host genes in response to viral infections. These results identified the commonly expressed virus-responsive genes and pathways during plant-virus interaction which might develop novel antiviral strategies for improving plant resistance to mixed viral infections.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00925-3.

Gene Expression Dynamics in Rice Peduncles at the Heading Stage

Improving grain yield in the staple food crop rice has been long sought goal of plant biotechnology. One of the traits with significant impact on rice breeding programs is peduncle elongation at the time of heading failing which leads to significant reduction in grain yield due to incomplete panicle exsertion. To decipher transcriptional dynamics and molecular players underlying peduncle elongation, we performed RNA sequencing analysis of elongating and non-elongating peduncles in two Indian cultivars, Swarna and Pokkali, at the time of heading. Along with genes associated with cell division and cell wall biosynthesis, we observed significant enrichment of genes associated with auxins, gibberellins, and brassinosteroid biosynthesis/signaling in the elongating peduncles before heading in both the genotypes. Similarly, genes associated with carbohydrate metabolism and mobilization, abiotic stress response along with cytokinin, abscisic acid, jasmonic acid, and ethylene biosynthesis/signaling were enriched in non-elongating peduncles post heading. Significant enrichment of genes belonging to key transcription factor families highlights their specialized roles in peduncle elongation and grain filling before and after heading, respectively. A comparison with anther/pollen development-related genes provided 76 candidates with overlapping roles in anther/pollen development and peduncle elongation. Some of these are important for carbohydrate remobilization to the developing grains. These can be engineered to combat with incomplete panicle exsertion in male sterile lines and manipulate carbohydrate dynamics in grasses. Overall, this study provides baseline information about potential target genes for engineering peduncle elongation with implications on plant height, biomass composition and grain yields in rice.

MiR319a-targeted PtoTCP20 regulates secondary growth via interactions with PtoWOX4 and PtoWND6 in Populus tomentosa

Secondary growth is a key characteristic of trees, which requires the coordination of multiple regulatory mechanisms including transcriptional regulators and microRNAs (miRNAs). However, the roles of microRNAs in the regulation of secondary growth need to be explored in depth. Here, the role of miR319a and its target, PtoTCP20, in the secondary growth of Populus tomentosa stem was investigated using genetic and molecular analyses. The expression level of miR319a gradually decreased from primary to secondary growth in P. tomentosa, while that of PtoTCP20 gradually increased. MiR319a overexpression in seedlings resulted in delayed secondary growth and decreased xylem production, while miR319a knockdown and PtoTCP20 overexpression promoted secondary growth and increased xylem production. Further analysis showed that PtoTCP20 interacted with PtoWOX4a and activated PtoWND6 transcription in vitro and in vivo. Our data show that PtoTCP20 controls vascular cambium proliferation by binding to PtoWOX4a, and promotes secondary xylem differentiation by activating PtoWND6 transcription, thereby regulating secondary growth in P. tomentosa. Our findings provide insight into the molecular mechanisms underlying secondary growth in trees.

Weeping candidate genes screened using comparative transcriptomic analysis of weeping and upright progeny in an F1 population of Prunus mume

Weeping is a specific plant architecture with high ornamental value. Despite the considerable importance of the weeping habit to landscaping applications and knowledge of plant architecture biology, little is known regarding the underlying molecular mechanisms. In this study, growth and phytohormone content were analyzed among the progeny of different branch types in an F1 mapping population of Prunus mume with varying plant architecture. Bulked segregant RNA sequencing was conducted to compare differences among progeny at a transcriptional level. The weeping habit appears to be a complex process regulated by a series of metabolic pathways, with photosynthesis and flavonoid biosynthesis highly enriched in differentially expressed genes between weeping and upright progeny. Based on functional annotation and homologous analyses, we identified 30 candidate genes related to weeping that merit further analysis, including 10 genes related to IAA and GA₃ biosynthesis, together with 6 genes related to secondary branch growth. The results of this study will facilitate further studies of the molecular mechanisms underlying the weeping habit in P. mume.

Regulation of Lignin Biosynthesis by Post-translational Protein Modifications

Post-translational modification of proteins exerts essential roles in many biological processes in plants. The function of these chemical modifications has been extensively characterized in many physiological processes, but how these modifications regulate lignin biosynthesis for wood formation remained largely unknown. Over the past decade, post-translational modification of several proteins has been associated with lignification. Phosphorylation, ubiquitination, glycosylation, and S-nitrosylation of transcription factors, monolignol enzymes, and peroxidases were shown to have primordial roles in the regulation of lignin biosynthesis. The main discoveries of post-translational modifications in lignin biosynthesis are discussed in this review.

Purine and pyrimidine dinucleoside polyphosphates differentially affect the phenylpropanoid pathway in Vitis vinifera L. cv. Monastrell suspension cultured cells

It is known that the concentration of dinucleoside polyphosphates (Np_nN's) in cells increases under stress and that adverse environmental factors induce biosynthesis of phenylpropanoids, which protect the plant against stress. Previously, we showed that purine Np_nN's such as Ap₃A and Ap₄A induce both the activity of enzymes of the phenylpropanoid pathway and the expression of relevant genes in Arabidopsis seedlings. Moreover, we showed that Ap₃A induced stilbene biosynthesis in Vitis vinifera cv. Monastrell suspension cultured cells. Data presented in this paper show that pyrimidine-containing Np_nN's also modify the biosynthesis of stilbenes, affecting the transcript level of genes encoding key enzymes of the phenylpropanoid pathway and of these, Up₄U caused the most effective accumulation of trans-resveratrol in the culture media. Similar effect was caused by Ap₃A and Gp₃G. Other pyrimidine Np_nN's, such as Cp₃C, Cp₄C, and Ap₄C, strongly inhibited the biosynthesis of stilbenes, but markedly (6- to 8-fold) induced the expression of the cinnamoyl-CoA reductase gene that controls lignin biosynthesis. Purine counterparts also clearly induced biosynthesis of trans-resveratrol and trans-piceid, but only slightly induced the expression of genes involved in lignin biosynthesis. In cells, Up₃U caused a greater accumulation of trans-resveratrol and trans-piceid than did Up₄U. Each of the Np_nN's studied induced expression of the gene encoding the resveratrol transporter VvABCG44, which operates within the Vitis vinifera cell membrane. AMP, GMP, UMP, and CMP, potential products of Np_nN degradation, did not affect the accumulation of stilbenes. The results obtained strongly support that Np_nN's play a role as signaling molecules in plants.

Affinity-based high-resolution analysis of DNA binding by VASCULAR-RELATED NAC-DOMAIN7 via fluorescence correlation spectroscopy

VASCULAR-RELATED NAC-DOMAIN7 (VND7) is the master transcription factor for vessel element differentiation in Arabidopsis thaliana. To identify the cis-acting sequence(s) bound by VND7, we employed fluorescence correlation spectroscopy (FCS) to find VND7-DNA interactions quantitatively. This identified an 18-bp sequence from the promoter of XYLEM CYSTEINE PEPTIDASE1 (XCP1), a direct target of VND7. A quantitative assay for binding affinity between VND7 and the 18-bp sequence revealed the core nucleotides contributing to specific binding between VND7 and the 18-bp sequence. Moreover, by combining the systematic evolution of ligands by exponential enrichment (SELEX) technique with known consensus sequences, we defined a motif termed the Ideal Core Structure for binding by VND7 (ICSV). We also used FCS to search for VND7 binding sequences in the promoter regions of other direct targets. Taking these data together, we proposed that VND7 preferentially binds to the ICSV sequence. Additionally, we found that substitutions among the core nucleotides affected transcriptional regulation by VND7 in vivo, indicating that the core nucleotides contribute to vessel-element-specific gene expression. Furthermore, our results demonstrate that FCS is a powerful tool for unveiling the DNA-binding properties of transcription factors.

SOBIR1/EVR prevents precocious initiation of fiber differentiation during wood development through a mechanism involving BP and ERECTA

In plants, secondary growth results in radial expansion of stems and roots, generating large amounts of biomass in the form of wood. Using genome-wide association studies (GWAS)-guided reverse genetics in Arabidopsis thaliana, we discovered SOBIR1/EVR, previously known to control plant immunoresponses and abscission, as a regulator of secondary growth. We present anatomical, genetic, and molecular evidence indicating that SOBIR1/EVR prevents the precocious differentiation of xylem fiber, a key cell type for wood development. SOBIR1/EVR acts through a mechanism that involves BREVIPEDICELLUS (BP) and ERECTA (ER), 2 proteins previously known to regulate xylem fiber development. We demonstrate that BP binds SOBIR1/EVR promoter and that SOBIR1/EVR expression is enhanced in bp mutants, suggesting a direct, negative regulation of BP over SOBIR1/EVR expression. We show that SOBIR1/EVR physically interacts with ER and that defects caused by the sobir1/evr mutation are aggravated by mutating ER, indicating that SOBIR1/EVR and ERECTA act together in the control of the precocious formation of xylem fiber development.

The Arabidopsis thaliana transcription factor MYB59 regulates calcium signalling during plant growth and stress response

Transcription factor MYB59 is involved in plant growth and stress responses by acting as negative regulator of Ca signalling and homeostasis. The Arabidopsis thaliana transcription factor MYB59 is induced by cadmium (Cd) and plays a key role in the regulation of cell cycle progression and root elongation, but its mechanism of action is poorly understood. We investigated the expression of MYB59 and differences between wild-type plants, the myb59 mutant and MYB59-overexpressing lines (obtained by transformation in the mutant genotype) during plant growth and in response to various forms of stress. We also compared the transcriptomes of wild-type and myb59 mutant plants to determine putative MYB59 targets. The myb59 mutant has longer roots, smaller leaves and smaller cells than wild-type plants and responds differently to stress in germination assay. Transcriptomic analysis revealed the upregulation in the myb59 mutant of multiple genes involved in calcium (Ca) homeostasis and signalling, including those encoding calmodulin-like proteins and Ca transporters. Notably, MYB59 was strongly induced by Ca deficiency, and the myb59 mutant was characterized by higher levels of cytosolic Ca in root cells and showed a modest alteration of Ca transient frequency in guard cells, associated with the absence of Ca-induced stomatal closure. These results indicate that MYB59 negatively regulates Ca homeostasis and signalling during Ca deficiency, thus controlling plant growth and stress responses.

Genome-wide identification and expression analysis of the MYB transcription factor in moso bamboo (Phyllostachys edulis)

The MYB family, one of the largest transcription factor (TF) families in the plant kingdom, plays vital roles in cell formation, morphogenesis and signal transduction, as well as responses to biotic and abiotic stresses. However, the underlying function of bamboo MYB TFs remains unclear. To gain insight into the status of these proteins, a total of 85 PeMYBs, which were further divided into 11 subgroups, were identified in moso bamboo (Phyllostachys edulis) by using a genome-wide search strategy. Gene structure analysis showed that PeMYBs were significantly different, with exon numbers varying from 4 to 13. Phylogenetic analysis indicated that PeMYBs clustered into 27 clades, of which the function of 18 clades has been predicted. In addition, almost all of the PeMYBs were differently expressed in leaves, panicles, rhizomes and shoots based on RNA-seq data. Furthermore, qRT-PCR analysis showed that 12 PeMYBs related to the biosynthesis and deposition of the secondary cell wall (SCW) were constitutively expressed, and their transcript abundance levels have changed significantly with increasing height of the bamboo shoots, for which the degree of lignification continuously increased. This result indicated that these PeMYBs might play fundamental roles in SCW thickening and bamboo shoot lignification. The present comprehensive and systematic study on the members of the MYB family provided a reference and solid foundation for further functional analysis of MYB TFs in moso bamboo.

Meta-analysis of transcriptomic responses to biotic and abiotic stress in tomato

A wide range of biotic stresses (BS) and abiotic stresses (AS) adversely affect plant growth and productivity worldwide. The study of individual genes cannot be considered as an effective approach for the understanding of tolerance mechanisms, since these stresses are frequent and often in combination with each other, and a large number of genes are involved in these mechanisms. The availability of high-throughput genomic data has enabled the discovery of the role of transcription factors (TFs) in regulatory networks. A meta-analysis of BS and AS responses was performed by analyzing a total of 391 microarray samples from 23 different experiments and 2,336 differentially expressed genes (DEGs) involved in multiple stresses were identified. We identified 1,862 genes differentially regulated in response to BS was much greater than that regulated by AS, 835 genes, and found 15.4% or 361 DEGs with the conserved expression between AS and BS. The greatest percent of genes related to the cellular process (>76% genes), metabolic process (>76% genes) and response to stimulus (>50%). About 4.2% of genes involved in BS and AS responses belonged to the TF families. We identified several genes, which encode TFs that play an important role in AS and BS responses. These proteins included Jasmonate Ethylene Response Factor 1 (JERF1), SlGRAS6, MYB48, SlERF4, EIL2, protein LATE ELONGATED HYPOCOTYL (LHY), SlERF1, WRKY 26, basic leucine zipper TF, inducer of CBF expression 1-like, pti6, EIL3 and WRKY 11. Six of these proteins, JERF1, MYB48, protein LHY, EIL3, EIL2 and SlGRAS6, play central roles in these mechanisms. This research promoted a new approach to clarify the expression profiles of various genes under different conditions in plants, detected common genes from differentially regulated in response to these conditions and introduced them as candidate genes for improving plant tolerance through genetic engineering approach.

Transcription factors PRE3 and WOX11 are involved in the formation of new lateral roots from secondary growth taproot in A. thaliana

The spatial deployment of lateral roots determines the ability of a plant to interact with the surrounding environment for nutrition and anchorage. This paper shows that besides the pericycle, the vascular cambium becomes active in Arabidopsis thaliana taproot at a later stage of development and is also able to form new lateral roots. To demonstrate the above, we implemented a two-step approach in which the first step leads to development of a secondary structure in A. thaliana taproot, and the second applies a mechanical stress on the vascular cambium to initiate formation of a new lateral root primordium. GUS staining showed PRE3, DR5 and WOX11 signals in the cambial zone of the root during new lateral root formation. An advanced level of wood formation, characterized by the presence of medullar rays, was achieved. Preliminary investigations suggest the involvement of auxin and two transcription factors (PRE3/ATBS1/bHLH135/TMO7 and WOX11) in the transition of some vascular cambium initials from a role as producers of xylem/phloem mother cells to founder cells of a new lateral root primordium.

Transcriptome profiling of Brassica napus stem sections in relation to differences in lignin content

BACKGROUND

Brassica crops are cultivated widely for human consumption and animal feed purposes, and oilseed rape/canola (Brassica napus and rapa) is the second most important oilseed worldwide. Because of its natural diversity and genetic complexity, genomics studies on oilseed rape will be a useful resource base to modify the quantity and quality of biomass in various crops, and therefore, should have a positive impact on lignocellulosic biofuel production. The objective of this study was to perform microarray analysis on two variable lignin containing oilseed rape cultivars to target novel genes and transcription factors of importance in Brassica lignin regulation for applied research.

RESULTS

To gain insight into the molecular networks controlling cell wall biosynthetic and regulatory events, we conducted lignin and microarray analysis of top and basal stem sections of brown seeded Brassica napus DH12075 and yellow seeded YN01-429 cultivars. A total of 9500 genes were differentially expressed 2-fold or higher in the stem between the cultivars, with a higher number of expressed genes in the basal section. Of the upregulated genes, many were transcription factors and a considerable number of these were associated with secondary wall synthesis and lignification in B. napus and other plant species. The three largest groups of transcription factors with differential expression were C2H2 and C3HC4 zinc fingers and bHLH. A significant number of genes related to lignin and carbohydrate metabolism also showed differential expression patterns between the stem sections of the two cultivars. Within the same cultivar, the number of upregulated genes was higher in the top section relative to the basal one.

CONCLUSION

In this study, we identified and established expression patterns of many new genes likely involved in cell wall biosynthesis and regulation. Some genes with known roles in other biochemical pathways were also identified to have a potential role in cell wall biosynthesis. This stem transcriptome profiling will allow for selecting novel regulatory and structural genes for functional characterization, a strategy which may provide tools for modifying cell wall composition to facilitate fermentation for biofuel production.

The role of plant hormones during grafting

For millennia, people have cut and joined different plant tissues together through a process known as grafting. By creating a chimeric organism, desirable properties from two plants combine to enhance disease resistance, abiotic stress tolerance, vigour or facilitate the asexual propagation of plants. In addition, grafting has been extremely informative in science for studying and identifying the long-distance movement of molecules. Despite its increasing use in horticulture and science, how plants undertake the process of grafting remains elusive. Here, we discuss specifically the role of eight major plant hormones during the wound healing and vascular formation process, two phenomena involved in grafting. We furthermore present the roles of these hormones during graft formation and highlight knowledge gaps and future areas of interest for the field of grafting biology.

Native vascular related NAC transcription factors are efficient regulator of multiple classes of secondary wall associated genes in banana

Secondary-wall deposition in xylem vessel elements is regulated by vascular-related NAC transcription factors (VNDs). We show that three banana VNDs (MusaVND1, MusaVND2 and MusaVND3) directly regulate multiple secondary-wall associated genes by binding to their 5'-upstream regulatory region. Transgenic banana harboring either P_MusaVND1:GUS, P_MusaVND2:GUS or P_MusaVND3:GUS showed specific GUS staining in lignified tissues. MusaVND1, MusaVND2 and MusaVND3 encodes transcriptional-activators as its C-terminal region drive expression of reporter genes in vivo in yeast. Purified MusaVND1, MusaVND2 and MusaVND3 proteins in gel shift assay bind to 19-bp secondary-wall NAC binding element (SNBE) while it fails to bind mutated SNBE. Putative SNBE sites in the 5'-upstream regulatory region of important secondary-wall associated genes related to programmed cell death (XCP1), cell-wall modification (IRX1/CesA8, IRX3/CesA7,IRX5/CesA4, IRX8, IRX10 and IRX12) and transcriptional regulation (MYB52, MYB48/59, MYB85, MYB58/72, MYB46, and MYB83) in banana was identified and mobility of these regulatory regions got retarded by MusaVND1, MusaVND2 and MusaVND3. Transcript level of these important secondary wall associated genes were elevated in transgenic banana overexpressing either MusaVND1, MusaVND2 or MusaVND3. Present study suggested promoters with prospective utilization in wall modification in banana (a potential biofuel crop) and suggest a complex transcriptional regulation of secondary wall deposition in plants.

Drought Tolerance Conferred in Soybean (Glycine max. L) by GmMYB84, a Novel R2R3-MYB Transcription Factor

MYB-type transcription factors (MYB TFs) play diverse roles in plant development and stress responses. However, the mechanisms underlying the actions of MYB TFs during stress response remain unclear. In this study we identified a R2R3-MYB TF in soybean (Glycine max), denoted GmMYB84, which contributes to drought resistance. Expression of GmMYB84 was induced by drought, salt stress, H2O2 and ABA. Compared with the wild type (WT), GmMYB84-overexpressing soybean mutants (OE lines) exhibited enhanced drought resistance with a higher survival rate, longer primary root length, greater proline and reactive oxygen species (ROS) contents, higher antioxidant enzyme activities [peroxidase (POD), catalase (CAT) and superoxide dismutase (SOD)], a lower dehydration rate and reduced malondialdehyde (MDA) content. We also found that ROS could induce SOD/POD/CAT activity in OE lines. In particular, we found that the optimal level of ROS is required for GmMYB84 to modulate primary root elongation. Some ROS-related genes were up-regulated under abiotic stress in GmMYB84 transgenic plants compared with the WT. Furthermore, electrophoretic mobility shift assay and luciferase reporter analysis demonstrated that GmMYB84 binds directly to the promoter of GmRBOHB-1 and GmRBOHB-2 genes. Based on this evidence, we propose a model for how GmMYB84, H2O2 and antioxidant enzymes work together to control root growth under both optimal and drought stress conditions.

Functional network analysis of genes differentially expressed during xylogenesis in soc1ful woody Arabidopsis plants

Many plant genes are known to be involved in the development of cambium and wood, but how the expression and functional interaction of these genes determine the unique biology of wood remains largely unknown. We used the soc1ful loss of function mutant - the woodiest genotype known in the otherwise herbaceous model plant Arabidopsis - to investigate the expression and interactions of genes involved in secondary growth (wood formation). Detailed anatomical observations of the stem in combination with mRNA sequencing were used to assess transcriptome remodeling during xylogenesis in wild-type and woody soc1ful plants. To interpret the transcriptome changes, we constructed functional gene association networks of differentially expressed genes using the STRING database. This analysis revealed functionally enriched gene association hubs that are differentially expressed in herbaceous and woody tissues. In particular, we observed the differential expression of genes related to mechanical stress and jasmonate biosynthesis/signaling during wood formation in soc1ful plants that may be an effect of greater tension within woody tissues. Our results suggest that habit shifts from herbaceous to woody life forms observed in many angiosperm lineages could have evolved convergently by genetic changes that modulate the gene expression and interaction network, and thereby redeploy the conserved wood developmental program.

Short day transcriptomic programming during induction of dormancy in grapevine

Bud dormancy in grapevine is an adaptive strategy for the survival of drought, high and low temperatures and freeze dehydration stress that limit the range of cultivar adaptation. Therefore, development of a comprehensive understanding of the biological mechanisms involved in bud dormancy is needed to promote advances in selection and breeding, and to develop improved cultural practices for existing grape cultivars. The seasonally indeterminate grapevine, which continuously develops compound axillary buds during the growing season, provides an excellent system for dissecting dormancy, because the grapevine does not transition through terminal bud development prior to dormancy. This study used gene expression patterns and targeted metabolite analysis of two grapevine genotypes that are short photoperiod responsive (Vitis riparia) and non-responsive (V. hybrid, Seyval) for dormancy development to determine differences between bud maturation and dormancy commitment. Grapevine gene expression and metabolites were monitored at seven time points under long (LD, 15 h) and short (SD, 13 h) day treatments. The use of age-matched buds and a small (2 h) photoperiod difference minimized developmental differences and allowed us to separate general photoperiod from dormancy specific gene responses. Gene expression profiles indicated three distinct phases (perception, induction and dormancy) in SD-induced dormancy development in V. riparia. Different genes from the NAC DOMAIN CONTAINING PROTEIN 19 and WRKY families of transcription factors were differentially expressed in each phase of dormancy. Metabolite and transcriptome analyses indicated ABA, trehalose, raffinose and resveratrol compounds have a potential role in dormancy commitment. Finally, a comparison between V. riparia compound axillary bud dormancy and dormancy responses in other species emphasized the relationship between dormancy and the expression of RESVERATROL SYNTHASE and genes associated with C3HC4-TYPE RING FINGER and NAC DOMAIN CONTAINING PROTEIN 19 transcription factors.

Exogenous adenosine 5'-phosphoramidate behaves as a signal molecule in plants; it augments metabolism of phenylpropanoids and salicylic acid in Arabidopsis thaliana seedlings

Cells contain various congeners of the canonical nucleotides. Some of these accumulate in cells under stress and may function as signal molecules. Their cellular levels are enzymatically controlled. Previously, we demonstrated a signaling function for diadenosine polyphosphates and cyclic nucleotides in Arabidopsis thaliana and grape, Vitis vinifera. These compounds increased the expression of genes for and the specific activity of enzymes of phenylpropanoid pathways resulting in the accumulation of certain products of these pathways. Here, we show that adenosine 5'-phosphoramidate, whose level can be controlled by HIT-family proteins, induced similar effects. This natural nucleotide, when added to A. thaliana seedlings, activated the genes for phenylalanine:ammonia lyase, 4-coumarate:coenzyme A ligase, cinnamate-4-hydroxylase, chalcone synthase, cinnamoyl-coenzyme A:NADP oxidoreductase and isochorismate synthase, which encode proteins catalyzing key reactions of phenylpropanoid pathways, and caused accumulation of lignins, anthocyanins and salicylic acid. Adenosine 5'-phosphofluoridate, a synthetic congener of adenosine 5'-phosphoramidate, behaved similarly. The results allow us to postulate that adenosine 5'-phosphoramidate should be considered as a novel signaling molecule.

Small tandem target mimic-mediated blockage of microRNA858 induces anthocyanin accumulation in tomato

Our work strongly suggests that microRNA858 regulates anthocyanin biosynthesis in tomato by modulating the expression of two R2R3 MYB transcription factors, underscoring the importance of microRNAs in the gene regulatory network controlling specialized metabolism in plants. The biological functions of microRNA858 (miR858), a recently identified small RNA, are not well understood. Here, we identified miR858 as a negative regulator of anthocyanin biosynthesis in tomato (Solanum lycopersicum). RNA ligase-mediated 5'RACE cleavage assay showed that miR858 mediates the cleavage of SlMYB7-like and SlMYB48-like transcripts in tomato. Expression analysis revealed an inverse correlation between the accumulation of miR858 and its target SlMYB7-like mRNA, in different tissues of tomato. Subsequently, a small tandem target mimic construct for blocking miR858 (STTM858) was generated and transformed into tomato. The majority of endogenous miR858 was blocked in STTM858 over-expressing tomato plants, whereas SlMYB7-like transcripts increased significantly. Concomitantly, upregulated expression was detected for several anthocyanin biosynthetic genes, including PAL, CHS, DFR, ANS and 3GT. As a result, anthocyanins were highly accumulated in young seedlings, leaves, stems and leaf buds of the transgenic plants under normal growth conditions. In addition, over-expression of STTM858 in tomato activated another MYB transcription factor, SlMYB48, implicating the possible involvement of SlMYB48 in anthocyanin biosynthesis.

Multifaceted regulations of gateway enzyme phenylalanine ammonia-lyase in the biosynthesis of phenylpropanoids

Phenylpropanoid biosynthesis in plants engenders a vast variety of aromatic metabolites critically important for their growth, development, and environmental adaptation. Some of these aromatic compounds have high economic value. Phenylalanine ammonia-lyase (PAL) is the first committed enzyme in the pathway; it diverts the central flux of carbon from the primary metabolism to the synthesis of myriad phenolics. Over the decades, many studies have shown that exquisite regulatory mechanisms at multiple levels control the transcription and the enzymatic activity of PALs. In this review, a current overview of our understanding of the complicated regulatory mechanisms governing the activity of PAL is presented; recent progress in unraveling its post-translational modifications, its metabolite feedback regulation, and its enzyme organization is highlighted.

Temporal aspects of copper homeostasis and its crosstalk with hormones

To cope with the dual nature of copper as being essential and toxic for cells, plants temporarily adapt the expression of copper homeostasis components to assure its delivery to cuproproteins while avoiding the interference of potential oxidative damage derived from both copper uptake and photosynthetic reactions during light hours. The circadian clock participates in the temporal organization of coordination of plant nutrition adapting metabolic responses to the daily oscillations. This timely control improves plant fitness and reproduction and holds biotechnological potential to drive increased crop yields. Hormonal pathways, including those of abscisic acid, gibberellins, ethylene, auxins, and jasmonates are also under direct clock and light control, both in mono and dicotyledons. In this review, we focus on copper transport in Arabidopsis thaliana and Oryza sativa and the presumable role of hormones in metal homeostasis matching nutrient availability to growth requirements and preventing metal toxicity. The presence of putative hormone-dependent regulatory elements in the promoters of copper transporters genes suggests hormonal regulation to match special copper requirements during plant development. Spatial and temporal processes that can be affected by hormones include the regulation of copper uptake into roots, intracellular trafficking and compartmentalization, and long-distance transport to developing vegetative and reproductive tissues. In turn, hormone biosynthesis and signaling are also influenced by copper availability, which suggests reciprocal regulation subjected to temporal control by the central oscillator of the circadian clock. This transcriptional regulatory network, coordinates environmental and hormonal signaling with developmental pathways to allow enhanced micronutrient acquisition efficiency.

Cellular events during interfascicular cambium ontogenesis in inflorescence stems of Arabidopsis

Development of cambium and its activity is important for our knowledge of the mechanism of secondary growth. Arabidopsis thaliana emerges as a good model plant for such a kind of study. Thus, this paper reports on cellular events taking place in the interfascicular regions of inflorescence stems of A. thaliana, leading to the development of interfascicular cambium from differentiated interfascicular parenchyma cells (IPC). These events are as follows: appearance of auxin accumulation, PIN1 gene expression, polar PIN1 protein localization in the basal plasma membrane and periclinal divisions. Distribution of auxin was observed to be higher in differentiating into cambium parenchyma cells compared to cells within the pith and cortex. Expression of PIN1 in IPC was always preceded by auxin accumulation. Basal localization of PIN1 was already established in the cells prior to their periclinal division. These cellular events initiated within parenchyma cells adjacent to the vascular bundles and successively extended from that point towards the middle region of the interfascicular area, located between neighboring vascular bundles. The final consequence of which was the closure of the cambial ring within the stem. Changes in the chemical composition of IPC walls were also detected and included changes of pectic epitopes, xyloglucans (XG) and extensins rich in hydroxyproline (HRGPs). In summary, results presented in this paper describe interfascicular cambium ontogenesis in terms of successive cellular events in the interfascicular regions of inflorescence stems of Arabidopsis.

Arabidopsis thaliana peroxidases involved in lignin biosynthesis: in silico promoter analysis and hormonal regulation

Phytohormones such as auxins, cytokinins, and brassinosteroids, act by means of a signaling cascade of transcription factors of the families NAC, MYB, AP2 (APETALA2), MADS and class III HD (homeodomain) Zip, regulating secondary growth. When the hormonal regulation of Zinnia elegans peroxidase (ZePrx), an enzyme involved in lignin biosynthesis, was studied, it was found that this peroxidase is sensitive to a plethora of hormones which control xylem lignification. In a previous study we sought Arabidopsis thaliana homologues to ZePrx. Peroxidases 4, 52, 49 and 72 are the four peroxidases that fulfill the restrictive conditions that a peroxidase involved in lignification must have. In the present study, we focus our attention on hormonal regulation in order to establish the minimal structural and regulatory elements contained in the promoter region which an AtPrx involved in lignification must have. The results indicate that of the four peroxidases selected in our previous study, the one most likely to be homologous to ZePrx is AtPrx52. The results suggest that hormones such as auxins, cytokinins and BRs directly regulate AtPrx52, and that the AtPrx52 promoter may be the target of the set of transcription factors (NAC, MYB, AP2 and class I and III HD Zip) which are up-regulated by these hormones during secondary growth. In addition, the AtPrx52 promoter contains multiple copies of all the putative cis-elements (the ACGT box, the OCS box, the OPAQ box, the L1BX, the MYCL box and the W box) known to confer regulation by NO and H2O2.

Induced transcriptional profiling of phenylpropanoid pathway genes increased flavonoid and lignin content in Arabidopsis leaves in response to microbial products

BACKGROUND

The production and use of biologically derived soil additives is one of the fastest growing sectors of the fertilizer industry. These products have been shown to improve crop yields while at the same time reducing fertilizer inputs to and nutrient loss from cropland. The mechanisms driving the changes in primary productivity and soil processes are poorly understood and little is known about changes in secondary productivity associated with the use of microbial products. Here we investigate secondary metabolic responses to a biologically derived soil additive by monitoring changes in the phenlypropanoid (PP) pathway in Arabidopsis thaliana.

RESULTS

This study was designed to test the influence of one of these products (Soil Builder™-AF, SB) on secondary metabolism after being applied at different times. One time (TI) application of SB to Arabidopsis increased the accumulation of flavonoids compared to multiple (TII) applications of the same products. Fourteen phenolic compounds including flavonols and anothocyanins were identified by mass spectrometry. Kaempferol-3,7-O-bis-α-L-rhamnoside and quercetin 3,7-dirhamnoside, the major compounds, increased 3-fold and 4-fold, respectively compared to control in the TI treatment. The most abundant anthocyanin was cyanidin 3-rhamnoglucoside, which increased 3-fold and 2-fold in TI compared to the control and TII, respectively. Simultaneously, the expression of genes coding for key enzymes in the PP pathway (phenylalanine ammonia lyase, cinnamate 4-hydroxylase, chalcone synthase, flavonoid-3'-O-hydroxylase, flavonol synthase1 and dihydroflavonol-4-reductase) and regulatory genes (production of anthocyanin pigment2, MYB12, MYB113, MYB114, EGL3, and TT8) were up-regulated in both treatments (TI and TII). Furthermore, application of TI and TII induced expression of the lignin pathway genes (hydroxyl cinamyl transferase, caffeyl-CoA O-methyl transferase, cinnamyl alcohol dehydrogenase, cinnamyl-CoA reductase, secondary wall-associated NAC domain protein1, MYB58 and MYB63 resulting in higher accumulation of lignin content compared to the control.

CONCLUSIONS

These results indicate that the additions of microbially based soil additives have a perceptible influence on phenylpropanoid pathway gene regulation and its production of secondary metabolites. These findings open an avenue of research to investigate the mode of action of microbially-based soil additives which may assist in the sustainable production of food, feed, fuel and fiber.

Transcriptome sequencing and profiling of expressed genes in cambial zone and differentiating xylem of Japanese cedar (Cryptomeria japonica)

BACKGROUND

Forest trees have ecological and economic importance, and Japanese cedar has highly valued wood attributes. Thus, studies of molecular aspects of wood formation offer practical information that may be used for screening and forward genetics approaches to improving wood quality.

RESULTS

After identifying expressed sequence tags in Japanese cedar tissue undergoing xylogenesis, we designed a custom cDNA microarray to compare expression of highly regulated genes throughout a growing season. This led to identification of candidate genes involved both in wood formation and later cessation of growth and dormancy. Based on homology to orthologous protein groups, the genes were assigned to functional classes. A high proportion of sequences fell into functional classes related to posttranscriptional modification and signal transduction, while transcription factors and genes involved in the metabolism of sugars, cell-wall synthesis and lignification, and cold hardiness were among other classes of genes identified as having a potential role in xylem formation and seasonal wood formation.

CONCLUSIONS

We obtained 55,051 unique sequences by next-generation sequencing of a cDNA library prepared from cambial meristem and derivative cells. Previous studies on conifers have identified unique sequences expressed in developing xylem, but this is the first comprehensive study utilizing a collection of expressed sequence tags for expression studies related to xylem formation in Japanese cedar, which belongs to a different lineage than the Pinaceae. Our characterization of these sequences should allow comparative studies of genome evolution and functional genetics of wood species.

Both cyclic-AMP and cyclic-GMP can act as regulators of the phenylpropanoid pathway in Arabidopsis thaliana seedlings

Cyclic nucleotides (cAMP and cGMP) are important signaling molecules that control a range of cellular functions and modulate different reactions. It is known that under abiotic or biotic stress plant cells synthesize these nucleotides and that they also enhance the activity of the phenylpropanoid pathway. Wondering what is the relation between these two facts, we investigated how the exogenously applied membrane-permeable derivatives, 8-Br-cAMP or 8-Br-cGMP, which are believed to act as the original cyclic nucleotides, affect the expression of the genes for and the specific activity of three enzymes of the phenylpropanoid pathway in Arabidopsis thaliana seedlings. We found that the expression of the genes of phenylalanine ammonia-lyase (PAL2), 4-coumarate:coenzyme A ligase (4CL1) and chalcone synthase (CHS), and the specific activities of PAL (EC 4.3.1.5), 4CL (EC 6.2.1.12) and CHS (EC 2.3.1.74) were induced in the same way by either of these cyclic nucleotides used at 5 μM concentration. None of the possible cAMP and cGMP degradation products (AMP, GMP, adenosine or guanosine) evoked such effects. Expression of PAL1, 4CL2 and 4CL3 were practically not affected. Although the investigated nucleotides induced rapid expression of the aforementioned enzymes, they did not affect the level of anthocyanins within the same period. We discuss the effects exerted by the exogenously administered cyclic nucleotides, their relation with stress and the role which the phenylpropanoid pathways the cyclic nucleotides may play in plants.

Sequencing, de novo assembly, annotation and SSR and SNP detection of sabaigrass (Eulaliopsis binata) transcriptome

Eulaliopsis binata is one of the best fiber grass plants for its high fiber quality and production. Large scale trancriptome sequencing of E. binata was first performed using mixed leaf samples of 20 wild clusters. A total of 26,438,832 clean reads were generated and were assembled into 59,134 isogenes with an average length of 845bp. BLAST against the NCBI non-redundant protein, KEGG and GO databases has classified these isogenes into functional categories for understanding gene functions and regulation pathways. Only 15.0% of the assembled isogenes were similar to known proteins and 24.4% has no hits in the nr protein database. The total isogenes and 5306 highly expressed isogenes were performed by BLASTx with the MAIZEWALL, the cell wall navigator and the PlantTFDB databases. A total of 6681 simple sequence repeats (SSRs) and 147,177 single nucleotide polymorphisms (SNPs) were detected in the isogenes and 5723 pairs of SSR primers were designed.

Transcriptional programming during cell wall maturation in the expanding Arabidopsis stem

BACKGROUND

Plant cell walls are complex dynamic structures that play a vital role in coordinating the directional growth of plant tissues. The rapid elongation of the inflorescence stem in the model plant Arabidopsis thaliana is accompanied by radical changes in cell wall structure and chemistry, but analysis of the underlying mechanisms and identification of the genes that are involved has been hampered by difficulties in accurately sampling discrete developmental states along the developing stem.

RESULTS

By creating stem growth kinematic profiles for individual expanding Arabidopsis stems we have been able to harvest and pool developmentally-matched tissue samples, and to use these for comparative analysis of global transcript profiles at four distinct phases of stem growth: the period of elongation rate increase, the point of maximum growth rate, the point of stem growth cessation and the fully matured stem. The resulting profiles identify numerous genes whose expression is affected as the stem tissues pass through these defined growth transitions, including both novel loci and genes identified in earlier studies. Of particular note is the preponderance of highly active genes associated with secondary cell wall deposition in the region of stem growth cessation, and of genes associated with defence and stress responses in the fully mature stem.

CONCLUSIONS

The use of growth kinematic profiling to create tissue samples that are accurately positioned along the expansion growth continuum of Arabidopsis inflorescence stems establishes a new standard for transcript profiling analyses of such tissues. The resulting expression profiles identify a substantial number of genes whose expression is correlated for the first time with rapid cell wall extension and subsequent fortification, and thus provide an important new resource for plant biologists interested in gene discovery related to plant biomass accumulation.

Adventitious root induction in Arabidopsis thaliana as a model for in vitro root organogenesis

Adventitious root formation, the development of roots on non-root tissue (e.g. leaves, hypocotyls and stems) is a critical step during micropropagation. Although root induction treatments are routinely used for a large number of species micropropagated in vitro as well as for in vivo cuttings, the mechanisms controlling adventitious rooting are still poorly understood. Researchers attempt to gain better insight into the molecular aspects by studying adventitious rooting in Arabidopsis thaliana. The existing assay involves etiolation of seedlings and measurements of de novo formed roots on the elongated hypocotyl. The etiolated hypocotyls express a novel auxin-controlled signal transduction pathway in which auxin response factors (ARFs), microRNAs and environmental conditions that drive adventitious rooting are integrated. An alternative assay makes use of so-called thin cell layers (TCL), excised strips of cells from the inflorescence stem of Arabidopsis thaliana. However, both the etiolated seedling system and the TCL assay are only distantly related to industrial rooting processes in which roots are induced on adult stem tissue. Here, we describe an adventitious root induction system that uses segments of the inflorescence stems of Arabidopsis thaliana, which have a histological structure similar to cuttings or in vitro micropropagated shoots. The system allows multiple treatments with chemicals as well as the evaluation of different environmental conditions on a large number of explants. It is therefore suitable for high throughput chemical screenings and experiments that require numerous data points for statistical analysis. Using this assay, the adventitious root induction capacity of classical auxins was evaluated and a differential response to the different auxins could be demonstrated. NAA, IBA and IAA stimulated adventitious rooting on the stem segment, whereas 2,4-D and picloram did not. Light conditions profoundly influenced the root induction capacity of the auxins. Additionally to the environmental control of adventitious root formation, we also investigated the spatial and temporal aspects of stem-based adventitious root organogenesis. To determine the cells involved in de novo root initiation on the adult stems, we adopted scanning electron microscopy, which allows the visualization of the auxin responsive stem tissue. Using this technique, direct (without callus interface) and indirect (with intermediate callus phase) organogenesis was readily distinguished. The described micro-stem segment system is also suitable for other non-woody species and it is a valuable tool to perform fast evaluations of different treatments to study adventitious root induction.

Transcriptome analysis of intraspecific competition in Arabidopsis thaliana reveals organ-specific signatures related to nutrient acquisition and general stress response pathways

BACKGROUND

Plants are sessile and therefore have to perceive and adjust to changes in their environment. The presence of neighbours leads to a competitive situation where resources and space will be limited. Complex adaptive responses to such situation are poorly understood at the molecular level.

RESULTS

Using microarrays, we analysed whole-genome expression changes in Arabidopsis thaliana plants subjected to intraspecific competition. The leaf and root transcriptome was strongly altered by competition. Differentially expressed genes were enriched in genes involved in nutrient deficiency (mainly N, P, K), perception of light quality, and responses to abiotic and biotic stresses. Interestingly, performance of the generalist insect Spodoptera littoralis on densely grown plants was significantly reduced, suggesting that plants under competition display enhanced resistance to herbivory.

CONCLUSIONS

This study provides a comprehensive list of genes whose expression is affected by intraspecific competition in Arabidopsis. The outcome is a unique response that involves genes related to light, nutrient deficiency, abiotic stress, and defence responses.

Characterization of microRNAs expressed during secondary wall biosynthesis in Acacia mangium

MicroRNAs (miRNAs) play critical regulatory roles by acting as sequence specific guide during secondary wall formation in woody and non-woody species. Although thousands of plant miRNAs have been sequenced, there is no comprehensive view of miRNA mediated gene regulatory network to provide profound biological insights into the regulation of xylem development. Herein, we report the involvement of six highly conserved amg-miRNA families (amg-miR166, amg-miR172, amg-miR168, amg-miR159, amg-miR394, and amg-miR156) as the potential regulatory sequences of secondary cell wall biosynthesis. Within this highly conserved amg-miRNA family, only amg-miR166 exhibited strong differences in expression between phloem and xylem tissue. The functional characterization of amg-miR166 targets in various tissues revealed three groups of HD-ZIP III: ATHB8, ATHB15, and REVOLUTA which play pivotal roles in xylem development. Although these three groups vary in their functions, -psRNA target analysis indicated that miRNA target sequences of the nine different members of HD-ZIP III are always conserved. We found that precursor structures of amg-miR166 undergo exhaustive sequence variation even within members of the same family. Gene expression analysis showed three key lignin pathway genes: C4H, CAD, and CCoAOMT were upregulated in compression wood where a cascade of miRNAs was downregulated. This study offers a comprehensive analysis on the involvement of highly conserved miRNAs implicated in the secondary wall formation of woody plants.

Transcriptional analysis of differentially expressed genes in response to stem inclination in young seedlings of pine

The gravitropic response in trees is a widely studied phenomenon, however understanding of the molecular mechanism involved remains unclear. The purpose of this work was to identify differentially expressed genes in response to inclination using a comparative approach for two conifer species. Young seedlings were subjected to inclination and samples were collected at four different times points. First, suppression subtractive hybridisation (SSH) was used to identify differentially regulated genes in radiata pine (Pinus radiata D. Don). cDNA libraries were constructed from the upper and lower part of inclined stems in a time course experiment, ranging from 2.5 h to 1 month. From a total of 3092 sequences obtained, 2203 elements were assembled, displaying homology to a public database. A total of 942 unigene elements were identified using bioinformatic tools after redundancy analysis. Of these, 614 corresponded to known function genes and 328 to unknown function genes, including hypothetical proteins. Comparative analysis between radiata pine and maritime pine (Pinus pinaster Ait.) was performed to validate the differential expression of relevant candidate genes using qPCR. Selected genes were involved in several functional categories: hormone regulation, phenylpropanoid pathway and signal transduction. This comparative approach for the two conifer species helped determine the molecular gene pattern generated by inclination, providing a set of Pinus gene signatures that may be involved in the gravitropic stress response. These genes may also represent relevant candidate genes involved in the gravitropic response and potentially in wood formation.

Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression

Background

Conifers have very large genomes (13 to 30 Gigabases) that are mostly uncharacterized although extensive cDNA resources have recently become available. This report presents a global overview of transcriptome variation in a conifer tree and documents conservation and diversity of gene expression patterns among major vegetative tissues.

Results

An oligonucleotide microarray was developed from Picea glauca and P. sitchensis cDNA datasets. It represents 23,853 unique genes and was shown to be suitable for transcriptome profiling in several species. A comparison of secondary xylem and phelloderm tissues showed that preferential expression in these vascular tissues was highly conserved among Picea spp. RNA-Sequencing strongly confirmed tissue preferential expression and provided a robust validation of the microarray design. A small database of transcription profiles called PiceaGenExpress was developed from over 150 hybridizations spanning eight major tissue types. In total, transcripts were detected for 92% of the genes on the microarray, in at least one tissue. Non-annotated genes were predominantly expressed at low levels in fewer tissues than genes of known or predicted function. Diversity of expression within gene families may be rapidly assessed from PiceaGenExpress. In conifer trees, dehydrins and late embryogenesis abundant (LEA) osmotic regulation proteins occur in large gene families compared to angiosperms. Strong contrasts and low diversity was observed in the dehydrin family, while diverse patterns suggested a greater degree of diversification among LEAs.

Conclusion

Together, the oligonucleotide microarray and the PiceaGenExpress database represent the first resource of this kind for gymnosperm plants. The spruce transcriptome analysis reported here is expected to accelerate genetic studies in the large and important group comprised of conifer trees.

Down-regulation of SlIAA15 in tomato altered stem xylem development and production of volatile compounds in leaf exudates

The Aux/IAA family genes encode short-lived nuclear proteins that function as transcriptional regulators in auxin signal transduction. Aux/IAA genes have been reported to control many processes of plant development. Our recent study showed that down-regulation of SlIAA15 in tomato reduced apical dominance, altered pattern of axillary shoot development, increased lateral root formation and leaves thickness. The SlIAA15 suppressed lines display strong reduction of trichome density, suggesting that SlIAA15 is involved in trichome formation. Here, we reported that SlIAA15-suppressed transgenic lines display increased number of xylem cells compared to wild-type plants. Moreover, the monoterpene content in trichome exudates are significantly reduced in SlIAA15 down-regulated leaves. The results provide the roles of SlIAA15 in production of volatile compounds in leaf exudates and xylem development, clearly indicating that members of the Aux/IAA gene family can play distinct and specific functions. 

Analysis of 2,297 expressed sequence tags (ESTs) from a cDNA library of flax (Linum ustitatissimum L.) bark tissue

Bast fibre crops are the second most important natural fibre crops following cotton. Of these, flax (Linum ustitatissimum L.) is the most widely planted in the world, with its fibre used for high quality linen textile. A cDNA library of flax bark tissues was constructed with the purpose of identifying genes involved in the Bast fibre development. A total of 2,297 unigene sequences were obtained from 3,200 randomly selected clones of the cDNA library. These sequences were grouped into 155 clusters and 2,142 singletons, which have been submitted to the GenBank databases. By putative functional annotation, 23.3% of these sequences were similar to known proteins in GenBank, 44.0% of these sequences were similar to unknown proteins, and 32.7% of these sequences showed no significant similarity to any other protein sequences in existing databases. Classified by the Gene Ontology, 24.8, 23.1 and 14.3% were assigned to molecular function, biological process, and cellular component GO terms, respectively. By further bioinformatics approaches, about 110 ESTs matched cell wall related genes in the MAIZEWALL database, representing 16 functional categories of all 19 categories, of which, the most abundant category was protein synthesis. Based on the PlantTFDB database, 39 of the 64 transcription factor families in the Arabidopsis thaliana genome were identified as being involved in flax cell wall formation. The sequences and bioinformatics analysis data generated in this paper will be useful for gene expression, cloning and genetic engineering studies to characterize bast fibre development and improve the properties of the bast fibres.

Comprehensive analysis of differential genes and miRNA profiles for discovery of topping-responsive genes in flue-cured tobacco roots

Decapitation/topping is an important cultivating measure for flue-cured tobacco, and diverse biology processes are changed to respond to the topping, such as hormonal balance, root development, source-sink relationship, ability of nicotine synthesis and stress tolerance. The purpose of this study was to clarify the molecular mechanism involved in the response of flue-cured tobacco to topping. The differentially expressed genes and micro RNAs (miRNAs) before and after topping were screened with a combination of suppression subtractive hybridization (SSH) and miRNA deep sequencing. In all, 560 differently expressed clones were sequenced by SSH, and then 129 high quality expressed sequence tags were acquired. These expressed sequence tags were mainly involved in secondary metabolism (13.5%), hormone metabolism (4%), signaling/transcription (17.5%), stress/defense (20%), protein metabolism (13%), carbon metabolism (7%), other metabolism (12%) and unknown function (13%). The results contribute new data to the list of possible candidate genes involved in the response of flue-cured tobacco to topping. NAC transcription factor, a differential gene identified by SSH, had been proved to have a role in the regulation of nicotine biosynthesis. High-throughput sequencing of two small RNA libraries in combination with SSH screening revealed 15 differential miRNAs whose target genes were identical to some differential genes identified in SSH, suggesting that miRNAs play a critical role in post-transcriptional gene regulation in the response of flue-cured tobacco to decapitation. Based on the role of these miRNAs and differential genes identified from SSH in response to topping, an miRNA mediated model for flue-cured tobacco in response to topping is proposed.

The multiple fuzzy origins of woodiness within Balsaminaceae using an integrated approach. Where do we draw the line?

BACKGROUND AND AIMS

The family Balsaminaceae is essentially herbaceous, except for some woodier species that can be described as 'woody' herbs or small shrubs. The family is nested within the so-called balsaminoid clade of Ericales, including the exclusively woody families Tetrameristaceae and Marcgraviaceae, which is sister to the remaining families of the predominantly woody order. A molecular phylogeny of Balsaminaceae is compared with wood anatomical observations to find out whether the woodier species are derived from herbaceous taxa (i.e. secondary woodiness), or whether woodiness in the family represents the ancestral state for the order (i.e. primary woodiness).

METHODS

Wood anatomical observations of 68 Impatiens species and Hydrocera triflora, of which 47 are included in a multigene phylogeny, are carried out using light and scanning electron microscopy and compared with the molecular phylogenetic insights.

KEY RESULTS

There is much continuous variation in wood development between the Impatiens species studied, making the distinction between herbaceousness and woodiness difficult. However, the most woody species, unambiguously considered as truly woody shrubs, all display paedomorphic wood features pointing to secondary woodiness. This hypothesis is further supported by the molecular phylogeny, demonstrating that these most woody species are derived from herbaceous (or less woody) species in at least five independent clades. Wood formation in H. triflora is mostly confined to the ribs of the stems and shows paedomorphic wood features as well, suggesting that the common ancestor of Balsaminaceae was probably herbaceous.

CONCLUSIONS

The terms 'herbaceousness' and 'woodiness' are notoriously difficult to use in Balsaminaceae. However, anatomical observations and molecular sequence data show that the woodier species are derived from less woody or clearly herbaceous species, demonstrating that secondary woodiness has evolved in parallel.