Regulatory microRNAs and phasiRNAs of paclitaxel biosynthesis in Taxus chinensis.

Paclitaxel (trade name Taxol) is a rare diterpenoid with anticancer activity isolated from Taxus. At present, paclitaxel is mainly produced by the semi-synthetic method using extract of Taxus tissues as raw materials. The studies of regulatory mechanisms in paclitaxel biosynthesis would promote the production of paclitaxel through tissue/cell culture approaches. Here, we systematically identified 990 transcription factors (TFs), 460 microRNAs (miRNAs), and 160 phased small interfering RNAs (phasiRNAs) in Taxus chinensis to explore their interactions and potential roles in regulation of paclitaxel synthesis. The expression levels of enzyme genes in cone and root were higher than those in leaf and bark. Nearly all enzyme genes in the paclitaxel synthesis pathway were significantly up-regulated after jasmonate treatment, except for GGPPS and CoA Ligase. The expression level of enzyme genes located in the latter steps of the synthesis pathway was significantly higher in female barks than in male. Regulatory TFs were inferred through co-expression network analysis, resulting in the identification of TFs from diverse families including MYB and AP2. Genes with ADP binding and copper ion binding functions were overrepresented in targets of miRNA genes. The miRNA targets were mainly enriched with genes in plant hormone signal transduction, mRNA surveillance pathway, cell cycle and DNA replication. Genes in oxidoreductase activity, protein-disulfide reductase activity were enriched in targets of phasiRNAs. Regulatory networks were further constructed including components of enzyme genes, TFs, miRNAs, and phasiRNAs. The hierarchical regulation of paclitaxel production by miRNAs and phasiRNAs indicates a robust regulation at post-transcriptional level. Our study on transcriptional and posttranscriptional regulation of paclitaxel synthesis provides clues for enhancing paclitaxel production using synthetic biology technology.

Allele-specific DNA methylation and gene expression during shoot organogenesis in tissue culture of hybrid poplar.

Plant tissue regeneration is critical for genetic transformation and genome editing techniques. During the regeneration process, changes in epigenetic modifications accompany the cell fate transition. However, how allele-specific DNA methylation in two haplotypes contributes to the transcriptional dynamics during regeneration remains elusive. Here we applied an inter-species hybrid poplar (Populus alba × P. glandulosa cv. 84 K) as a system to characterize the DNA methylation landscape during de novo shoot organogenesis at allele level. Both direct and indirect shoot organogenesis showed a reduction in genome-wide DNA methylation. At gene level, non-expressed genes were hypermethylated in comparison with expressed genes. Among the genes exhibiting significant correlations between levels of DNA methylation and gene expression, the expression patterns of 75% of genes were negatively correlated with DNA methylation in the CG context, whereas the correlation patterns in the CHH context were the reverse. The allele-biased DNA methylation was consistent during shoot organogenesis, with fewer than one-thousandth of allele-specific methylation regions shifted. Analysis of allele-specific expression revealed that there were only 1909 genes showing phase-dependent allele-biased expression in the regeneration process, among which the allele pairs with greater differences in transcription factor binding sites at promoter regions exhibited greater differences in allele expression. Our results indicated a relatively independent transcriptional regulation in two subgenomes during shoot organogenesis, which was contributed by cis-acting genomic and epigenomic variations.

Multiplex Editing of the Nucleoredoxin1 Tandem Array in Poplar: From Small Indels to Translocations and Complex Inversions.

The CRISPR-Cas9 system has been deployed for precision mutagenesis in an ever-growing number of species, including agricultural crops and forest trees. Its application to closely linked genes with extremely high sequence similarities has been less explored. In this study, we used CRISPR-Cas9 to mutagenize a tandem array of seven Nucleoredoxin1 (NRX1) genes spanning ∼100 kb in Populus tremula × Populus alba. We demonstrated efficient multiplex editing with one single guide RNA in 42 transgenic lines. The mutation profiles ranged from small insertions and deletions and local deletions in individual genes to large genomic dropouts and rearrangements spanning tandem genes. We also detected complex rearrangements including translocations and inversions resulting from multiple cleavage and repair events. Target capture sequencing was instrumental for unbiased assessments of repair outcomes to reconstruct unusual mutant alleles. The work highlights the power of CRISPR-Cas9 for multiplex editing of tandemly duplicated genes to generate diverse mutants with structural and copy number variations to aid future functional characterization.

Genome Sequence and Comparative Analysis of Colletotrichum gloeosporioides Isolated from Liriodendron Leaves.

Colletotrichum gloeosporioides is a hemibiotrophic pathogen causing significant losses to economically important crops and forest trees, including Liriodendron. To explore the interaction between C. gloeosporioides and Liriodendron and to identify the candidate genes determining the pathogenesis, we sequenced and assembled the whole genome of C. gloeosporioides Lc1 (CgLc1) using PacBio and Illumina next generation sequencing and performed a comparative genomic analysis between CgLc1 and Cg01, the latter being a described endophytic species of the C. gloeosporioides complex. Gene structure prediction identified 15,744 protein-coding genes and 837 noncoding RNAs. Species-specific genes were characterized using an ortholog analysis followed by a pathway enrichment analysis, which showed that genes specific to CgLc1 were enriched for the arginine biosynthetic process. Furthermore, genome synteny analysis revealed that most of the protein-coding genes fell into collinear blocks. However, two clusters of polyketide synthase genes were identified to be specific for CgLc1, suggesting that they might have an important role in virulence control. Transcriptional regulators coexpressed with polyketide synthase genes were detected through a Weighted Correlation Network Analysis. Taken together, this work provides new insight into the virulence- and pathogenesis-associated genes present in C. gloeosporioides and its possible lifestyle.

Compensatory Guaiacyl Lignin Biosynthesis at the Expense of Syringyl Lignin in 4CL1-Knockout Poplar.

The lignin biosynthetic pathway is highly conserved in angiosperms, yet pathway manipulations give rise to a variety of taxon-specific outcomes. Knockout of lignin-associated 4-coumarate:CoA ligases (4CLs) in herbaceous species mainly reduces guaiacyl (G) lignin and enhances cell wall saccharification. Here we show that CRISPR-knockout of 4CL1 in poplar (Populus tremula × alba) preferentially reduced syringyl (S) lignin, with negligible effects on biomass recalcitrance. Concordant with reduced S-lignin was downregulation of ferulate 5-hydroxylases (F5Hs). Lignification was largely sustained by 4CL5, a low-affinity paralog of 4CL1 typically with only minor xylem expression or activity. Levels of caffeate, the preferred substrate of 4CL5, increased in line with significant upregulation of caffeoyl shikimate esterase1 Upregulation of caffeoyl-CoA O-methyltransferase1 and downregulation of F5Hs are consistent with preferential funneling of 4CL5 products toward G-lignin biosynthesis at the expense of S-lignin. Thus, transcriptional and metabolic adaptations to 4CL1-knockout appear to have enabled 4CL5 catalysis at a level sufficient to sustain lignification. Finally, genes involved in sulfur assimilation, the glutathione-ascorbate cycle, and various antioxidant systems were upregulated in the mutants, suggesting cascading responses to perturbed thioesterification in lignin biosynthesis.

Tubulins, rhythms and cell walls in poplar leaves: it's all in the timing.

Plant cell walls exhibit architectural and compositional changes throughout their development and in response to external cues. While tubulins are involved in cell wall biogenesis, much remains unknown about the scope of their involvement during the orchestration of this resource-demanding process. A transgenic approach coupled with cell wall compositional analysis, RNA-seq and mining of publicly available diurnal gene expression data was used to assess the involvement of tubulins in poplar leaf cell wall biogenesis. Leaf cell walls of transgenic poplar lines with constitutive overexpression of α-tubulin (TUA) exhibited an increased abundance of homogalacturonan, along with a reduction in xylose. These changes were traced to altered expression of UDP-glucuronic acid decarboxylase (GADC) in the transgenic leaves. A model is postulated by which altered diurnal control of TUA through its constitutive overexpression led to a metabolic tradeoff affecting cellular utilization of GADC substrate UDP-glucuronic acid. While there were no effects on cellulose, hemicellulose or lignin abundance, subtle effects on hemicellulose composition and associated gene expression were noted. In addition, expression and enzymatic activity of pectin methylesterase (PME) decreased in the transgenic leaves. The change is discussed in a context of increased levels of PME substrate homogalacturonan, slow stomatal kinetics and the fate of PME product methanol. Since stomatal opening and closing depend on fundamentally contrasting microtubule dynamics, the slowing of both processes in the transgenic lines as previously reported appears to be directly related to underlying cell wall compositional changes that were caused by tubulin manipulation.

Tubulin C-terminal Post-translational Modifications Do Not Occur in Wood Forming Tissue of Populus.

Cortical microtubules (MTs) are evolutionarily conserved cytoskeletal components with specialized roles in plants, including regulation of cell wall biogenesis. MT functions and dynamics are dictated by the composition of their monomeric subunits, α- (TUA) and ß-tubulins (TUB), which in animals and protists are subject to both transcriptional regulation and post-translational modifications (PTM). While spatiotemporal regulation of tubulin gene expression has been reported in plants, whether and to what extent tubulin PTMs occur in these species remain poorly understood. We chose the woody perennial Populus for investigation of tubulin PTMs in this study, with a particular focus on developing xylem where high tubulin transcript levels support MT-dependent secondary cell wall deposition. Mass spectrometry and immunodetection concurred that detyrosination, non-tyrosination and glutamylation were essentially absent in tubulins isolated from wood-forming tissues of P. deltoides and P. tremula ×alba. Label-free quantification of tubulin isotypes and RNA-Seq estimation of tubulin transcript abundance were largely consistent with transcriptional regulation. However, two TUB isotypes were detected at noticeably lower levels than expected based on RNA-Seq transcript abundance in both Populus species. These findings led us to conclude that MT composition during wood formation depends exclusively on transcriptional and, to a lesser extent, translational regulation of tubulin isotypes.

Drought response transcriptomes are altered in poplar with reduced tonoplast sucrose transporter expression.

Transgenic Populus tremula x alba (717-1B4) plants with reduced expression of a tonoplast sucrose efflux transporter, PtaSUT4, exhibit reduced shoot growth compared to wild type (WT) under sustained mild drought. The present study was undertaken to determine whether SUT4-RNAi directly or indirectly altered poplar predisposition and/or response to changes in soil water availability. While sucrose and hexose levels were constitutively elevated in shoot organs, expression responses to drought were most altered in the root tips of SUT4-RNAi plants. Prior to any drought treatment, constitutively elevated transcript levels of abscisic acid biosynthetic genes and bark/vegetative storage proteins suggested altered metabolism in root tips of RNAi plants. Stronger drought-stimulation of stress-inducible genes encoding late-embryogenesis-abundant proteins in transgenic roots was consistent with increased vulnerability to soil drying. Transcript evidence suggested an RNAi effect on intercellular water trafficking by aquaporins in stem xylem during soil drying and recovery. Co-expression network analysis predicted altered integration of abscisic acid sensing/signaling with ethylene and jasmonate sensing/signaling in RNAi compared to WT roots. The overall conclusion is that steepened shoot-root sugar gradient in RNAi plants increased sensitivity of root tips to decreasing soil water availability.

CRISPRing into the woods.

The CRISPR/Cas9 technology is a welcome breakthrough for genome editing, owing to its precision, efficiency, versatility and ease of adoption. We recently reported the first application of CRISPR/Cas9 for biallelic mutations in stably transformed Populus, extending the species range of this powerful technology to woody perennials. An underappreciated obstacle in genome editing of outcrossing species is the frequent occurrence of sequence polymorphisms that can render CRISPR/Cas9 unproductive. We discuss experimental evidence as well as genome-wide computational analysis to demonstrate the sensitivity of CRISPR/Cas9 to allelic heterozygosity, and highlight tools and strategies that can help deal with such sequence polymorphisms. With its specificity, CRISPR/Cas9 offers a less equivocal means than previous approaches for discerning functional redundancy of paralogous genes that are prevalent in plant genomes. Continuing improvements of the CRISPR/Cas9 system for multiplex genome engineering should facilitate these efforts. The paradigm shift brought about by CRISPR/Cas9 promises to accelerate not only basic research but also applied crop improvement progress.

Tubulin perturbation leads to unexpected cell wall modifications and affects stomatal behaviour in Populus.

Cortical microtubules are integral to plant morphogenesis, cell wall synthesis, and stomatal behaviour, presumably by governing cellulose microfibril orientation. Genetic manipulation of tubulins often leads to abnormal plant development, making it difficult to probe additional roles of cortical microtubules in cell wall biogenesis. Here, it is shown that expressing post-translational C-terminal modification mimics of α-tubulin altered cell wall characteristics and guard cell dynamics in transgenic Populus tremula x alba that otherwise appear normal. 35S promoter-driven transgene expression was high in leaves but unusually low in xylem, suggesting high levels of tubulin transgene expression were not tolerated in wood-forming tissues during regeneration of transformants. Cellulose, hemicellulose, and lignin contents were unaffected in transgenic wood, but expression of cell wall-modifying enzymes, and extractability of lignin-bound pectin and xylan polysaccharides were increased in developing xylem. The results suggest that pectin and xylan polysaccharides deposited early during cell wall biogenesis are more sensitive to subtle tubulin perturbation than cellulose and matrix polysaccharides deposited later. Tubulin perturbation also affected guard cell behaviour, delaying drought-induced stomatal closure as well as light-induced stomatal opening in leaves. Pectins have been shown to confer cell wall flexibility critical for reversible stomatal movement, and results presented here are consistent with microtubule involvement in this process. Taken together, the data show the value of growth-compatible tubulin perturbations for discerning microtubule functions, and add to the growing body of evidence for microtubule involvement in non-cellulosic polysaccharide assembly during cell wall biogenesis.

Bayesian phylogeny of sucrose transporters: ancient origins, differential expansion and convergent evolution in monocots and dicots.

Sucrose transporters (SUTs) are essential for the export and efficient movement of sucrose from source leaves to sink organs in plants. The angiosperm SUT family was previously classified into three or four distinct groups, Types I, II (subgroup IIB), and III, with dicot-specific Type I and monocot-specific Type IIB functioning in phloem loading. To shed light on the underlying drivers of SUT evolution, Bayesian phylogenetic inference was undertaken using 41 sequenced plant genomes, including seven basal lineages at key evolutionary junctures. Our analysis supports four phylogenetically and structurally distinct SUT subfamilies, originating from two ancient groups (AG1 and AG2) that diverged early during terrestrial colonization. In both AG1 and AG2, multiple intron acquisition events in the progenitor vascular plant established the gene structures of modern SUTs. Tonoplastic Type III and plasmalemmal Type II represent evolutionarily conserved descendants of AG1 and AG2, respectively. Type I and Type IIB were previously thought to evolve after the dicot-monocot split. We show, however, that divergence of Type I from Type III SUT predated basal angiosperms, likely associated with evolution of vascular cambium and phloem transport. Type I SUT was subsequently lost in monocots along with vascular cambium, and independent evolution of Type IIB coincided with modified monocot vasculature. Both Type I and Type IIB underwent lineage-specific expansion. In multiple unrelated taxa, the newly-derived SUTs exhibit biased expression in reproductive tissues, suggesting a functional link between phloem loading and reproductive fitness. Convergent evolution of Type I and Type IIB for SUT function in phloem loading and reproductive organs supports the idea that differential vascular development in dicots and monocots is a strong driver for SUT family evolution in angiosperms.

Condensed tannin biosynthesis and polymerization synergistically condition carbon use, defense, sink strength and growth in Populus.

The partitioning of carbon for growth, storage and constitutive chemical defenses is widely framed in terms of a hypothetical sink-source differential that varies with nutrient supply. According to this framework, phenolics accrual is passive and occurs in source leaves when normal sink growth is not sustainable due to a nutrient limitation. In assessing this framework, we present gene and metabolite evidence that condensed tannin (CT) accrual is strongest in sink leaves and sequesters carbon in a way that impinges upon foliar sink strength and upon phenolic glycoside (PG) accrual in Populus. The work was based on two Populus fremontii × angustifolia backcross lines with contrasting rates of CT accrual and growth, and equally large foliar PG reserves. However, foliar PG accrual was developmentally delayed in the high-CT, slow-growth line (SG), and nitrogen-limitation led to increased foliar PG accrual only in the low-CT, fast-growth line (FG). Metabolite profiling of developing leaves indicated comparatively carbon-limited amino acid metabolism, depletion of several Krebs cycle intermediates and reduced organ sink strength in SG. Gene profiling indicated that CT synthesis decreased as leaves expanded and PGs increased. A most striking finding was that the nitrogenous monoamine phenylethylamine accumulated only in leaves of SG plants. The potential negative impact of CT hyper-accumulation on foliar sink strength, as well as a mechanism for phenylethylamine involvement in CT polymerization in Populus are discussed. Starch accrual in source leaves and CT accrual in sink leaves of SG may both contribute to the maintenance of a slow-growth phenotype suited to survival in nutrient-poor habitats.

Constitutively elevated salicylic acid levels alter photosynthesis and oxidative state but not growth in transgenic populus.

Salicylic acid (SA) has long been implicated in plant responses to oxidative stress. SA overproduction in Arabidopsis thaliana leads to dwarfism, making in planta assessment of SA effects difficult in this model system. We report that transgenic Populus tremula × alba expressing a bacterial SA synthase hyperaccumulated SA and SA conjugates without negative growth consequences. In the absence of stress, endogenously elevated SA elicited widespread metabolic and transcriptional changes that resembled those of wild-type plants exposed to oxidative stress-promoting heat treatments. Potential signaling and oxidative stress markers azelaic and gluconic acids as well as antioxidant chlorogenic acids were strongly coregulated with SA, while soluble sugars and other phenylpropanoids were inversely correlated. Photosynthetic responses to heat were attenuated in SA-overproducing plants. Network analysis identified potential drivers of SA-mediated transcriptome rewiring, including receptor-like kinases and WRKY transcription factors. Orthologs of Arabidopsis SA signaling components NON-EXPRESSOR OF PATHOGENESIS-RELATED GENES1 and thioredoxins were not represented. However, all members of the expanded Populus nucleoredoxin-1 family exhibited increased expression and increased network connectivity in SA-overproducing Populus, suggesting a previously undescribed role in SA-mediated redox regulation. The SA response in Populus involved a reprogramming of carbon uptake and partitioning during stress that is compatible with constitutive chemical defense and sustained growth, contrasting with the SA response in Arabidopsis, which is transient and compromises growth if sustained.

Genome-wide analysis of the complex transcriptional networks of rice developing seeds.

BACKGROUND: The development of rice (Oryza sativa) seed is closely associated with assimilates storage and plant yield, and is fine controlled by complex regulatory networks. Exhaustive transcriptome analysis of developing rice embryo and endosperm will help to characterize the genes possibly involved in the regulation of seed development and provide clues of yield and quality improvement. PRINCIPAL FINDINGS: Our analysis showed that genes involved in metabolism regulation, hormone response and cellular organization processes are predominantly expressed during rice development. Interestingly, 191 transcription factor (TF)-encoding genes are predominantly expressed in seed and 59 TFs are regulated during seed development, some of which are homologs of seed-specific TFs or regulators of Arabidopsis seed development. Gene co-expression network analysis showed these TFs associated with multiple cellular and metabolism pathways, indicating a complex regulation of rice seed development. Further, by employing a cold-resistant cultivar Hanfeng (HF), genome-wide analyses of seed transcriptome at normal and low temperature reveal that rice seed is sensitive to low temperature at early stage and many genes associated with seed development are down-regulated by low temperature, indicating that the delayed development of rice seed by low temperature is mainly caused by the inhibition of the development-related genes. The transcriptional response of seed and seedling to low temperature is different, and the differential expressions of genes in signaling and metabolism pathways may contribute to the chilling tolerance of HF during seed development. CONCLUSIONS: These results provide informative clues and will significantly improve the understanding of rice seed development regulation and the mechanism of cold response in rice seed.

Rice ABI5-Like1 regulates abscisic acid and auxin responses by affecting the expression of ABRE-containing genes.

Abscisic acid (ABA) regulates plant development and is crucial for plant responses to biotic and abiotic stresses. Studies have identified the key components of ABA signaling in Arabidopsis (Arabidopsis thaliana), some of which regulate ABA responses by the transcriptional regulation of downstream genes. Here, we report the functional identification of rice (Oryza sativa) ABI5-Like1 (ABL1), which is a basic region/leucine zipper motif transcription factor. ABL1 is expressed in various tissues and is induced by the hormones ABA and indole-3-acetic acid and stress conditions including salinity, drought, and osmotic pressure. The ABL1 deficiency mutant, abl1, shows suppressed ABA responses, and ABL1 expression in the Arabidopsis abi5 mutant rescued the ABA sensitivity. The ABL1 protein is localized to the nucleus and can directly bind ABA-responsive elements (ABREs; G-box) in vitro. A gene expression analysis by DNA chip hybridization confirms that a large proportion of down-regulated genes of abl1 are involved in stress responses, consistent with the transcriptional activating effects of ABL1. Further studies indicate that ABL1 regulates the plant stress responses by regulating a series of ABRE-containing WRKY family genes. In addition, the abl1 mutant is hypersensitive to exogenous indole-3-acetic acid, and some ABRE-containing genes related to auxin metabolism or signaling are altered under ABL1 deficiency, suggesting that ABL1 modulates ABA and auxin responses by directly regulating the ABRE-containing genes.

Genome-wide analysis revealed the complex regulatory network of brassinosteroid effects in photomorphogenesis.

Light and brassinosteroids (BRs) have been proved to be crucial in regulating plant growth and development; however, the mechanism of how they synergistically function is still largely unknown. To explore the underlying mechanisms in photomorphogenesis, genome-wide analyses were carried out through examining the gene expressions of the dark-grown WT or BR biosynthesis-defective mutant det2 seedlings in the presence of light stimuli or exogenous Brassinolide (BL). Results showed that BR deficiency stimulates, while BL treatment suppresses, the expressions of light-responsive genes and photomorphogenesis, confirming the negative effects of BR in photomorphogenesis. This is consistent with the specific effects of BR on the expression of genes involved in cell wall modification, cellular metabolism and energy utilization during dark-light transition. Further analysis revealed that hormone biosynthesis and signaling-related genes, especially those of auxin, were altered under BL treatment or light stimuli, indicating that BR may modulate photomorphogenesis through synergetic regulation with other hormones. Additionally, suppressed ubiquitin-cycle pathway during light-dark transition hinted the presence of a complicated network among light, hormone, and protein degradation. The study provides the direct evidence of BR effects in photomorphogenesis and identified the genes involved in BR and light signaling pathway, which will help to elucidate the molecular mechanism of plant photomorphogenesis.

Global analysis of gene expression profiles in Brassica napus developing seeds reveals a conserved lipid metabolism regulation with Arabidopsis thaliana.

In order to study Brassica napus fatty acid (FA) metabolism and relevant regulatory networks, a systematic identification of fatty acid (FA) biosynthesis-related genes was conducted. Following gene identification, gene expression profiles during B. napus seed development and FA metabolism were performed by cDNA chip hybridization (>8000 EST clones from seed). The results showed that FA biosynthesis and regulation, and carbon flux, were conserved between B. napus and Arabidopsis. However, a more critical role of starch metabolism was detected for B. napus seed FA metabolism and storage-component accumulation when compared with Arabidopsis. In addition, a crucial stage for the transition of seed-to-sink tissue was 17-21 d after flowering (DAF), whereas FA biosynthesis-related genes were highly expressed primarily at 21 DAF. Hormone (auxin and jasmonate) signaling is found to be important for FA metabolism. This study helps to reveal the global regulatory network of FA metabolism in developing B. napus seeds.

Characterization and expression profiles of miRNAs in rice seeds.

Small RNAs (sRNAs) are common and effective modulators of gene expression in eukaryotic organisms. To characterize the sRNAs expressed during rice seed development, massively parallel signature sequencing (MPSS) was performed, resulting in the obtainment of 797,399 22-nt sequence signatures, of which 111,161 are distinct ones. Analysis on the distributions of sRNAs on chromosomes showed that most sRNAs originate from interspersed repeats that mainly consist of transposable elements, suggesting the major function of sRNAs in rice seeds is transposon silencing. Through integrative analysis, 26 novel miRNAs and 12 miRNA candidates were identified. Further analysis on the expression profiles of the known and novel miRNAs through hybridizing the generated chips revealed that most miRNAs were expressed preferentially in one or two rice tissues. Detailed comparison of the expression patterns of miRNAs and corresponding target genes revealed the negative correlation between them, while few of them are positively correlated. In addition, differential accumulations of miRNAs and corresponding miRNA*s suggest the functions of miRNA*s other than being passenger strands of mature miRNAs, and in regulating the miRNA functions.