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.

Inferring the Regulatory Network of miRNAs on Terpene Trilactone Biosynthesis Affected by Environmental Conditions.

As a medicinal tree species, ginkgo (Ginkgo biloba L.) and terpene trilactones (TTLs) extracted from its leaves are the main pharmacologic activity constituents and important economic indicators of its value. The accumulation of TTLs is known to be affected by environmental stress, while the regulatory mechanism of environmental response mediated by microRNAs (miRNAs) at the post-transcriptional levels remains unclear. Here, we focused on grafted ginkgo grown in northwestern, southwestern, and eastern-central China and integrally analyzed RNA-seq and small RNA-seq high-throughput sequencing data as well as metabolomics data from leaf samples of ginkgo clones grown in natural environments. The content of bilobalide was highest among detected TTLs, and there was more than a twofold variation in the accumulation of bilobalide between growth conditions. Meanwhile, transcriptome analysis found significant differences in the expression of 19 TTL-related genes among ginkgo leaves from different environments. Small RNA sequencing and analysis showed that 62 of the 521 miRNAs identified were differentially expressed among different samples, especially the expression of miRN50, miR169h/i, and miR169e was susceptible to environmental changes. Further, we found that transcription factors (ERF, MYB, C3H, HD-ZIP, HSF, and NAC) and miRNAs (miR319e/f, miRN2, miRN54, miR157, miR185, and miRN188) could activate or inhibit the expression of TTL-related genes to participate in the regulation of terpene trilactones biosynthesis in ginkgo leaves by weighted gene co-regulatory network analysis. Our findings provide new insights into the understanding of the regulatory mechanism of TTL biosynthesis but also lay the foundation for ginkgo leaves' medicinal value improvement under global change.

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.

The Metasequoia genome and evolutionary relationships among redwoods.

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

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.

Evidences for a role of two Y-specific genes in sex determination in Populus deltoides.

Almost all plants in the genus Populus are dioecious (i.e. trees are either male or female), but it is unknown whether dioecy evolved in a common ancestor or independently in different subgenera. Here, we sequence the small peritelomeric X- and Y-linked regions of P. deltoides chromosome XIX. Two genes are present only in the Y-linked region. One is a duplication of a non-Y-linked, female-specifically expressed response regulator, which produces siRNAs that block this gene's expression, repressing femaleness. The other is an LTR/Gypsy transposable element family member, which generates long non-coding RNAs. Overexpression of this gene in A. thaliana promotes androecium development. We also find both genes in the sex-determining region of P. simonii, a different poplar subgenus, which suggests that they are both stable components of poplar sex-determining systems. By contrast, only the duplicated response regulator gene is present in the sex-linked regions of P. davidiana and P. tremula. Therefore, findings in our study suggest dioecy may have evolved independently in different poplar subgenera.

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.

Author Correction: Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation.

In the Supplementary Information file originally published with this Letter, the authors mistakenly omitted Supplementary Table 14; this has now been amended.

Salicylic acid activates poplar defense against the biotrophic rust fungus Melampsora larici-populina via increased biosynthesis of catechin and proanthocyanidins.

Poplar trees synthesize flavan-3-ols (catechin and proanthocyanidins) as a defense against foliar rust fungi, but the regulation of this defense response is poorly understood. Here, we investigated the role of hormones in regulating flavan-3-ol accumulation in poplar during rust infection. We profiled levels of defense hormones, signaling genes, and flavan-3-ol metabolites in black poplar leaves at different stages of rust infection. Hormone levels were manipulated by external sprays, genetic engineering, and drought to reveal their role in rust fungal defenses. Levels of salicylic acid (SA), jasmonic acid, and abscisic acid increased in rust-infected leaves and activated downstream signaling, with SA levels correlating closely with those of flavan-3-ols. Pretreatment with the SA analog benzothiadiazole increased flavan-3-ol accumulation by activating the MYB-bHLH-WD40 complex and reduced rust proliferation. Furthermore, transgenic poplar lines overproducing SA exhibited higher amounts of flavan-3-ols constitutively via the same transcriptional activation mechanism. These findings suggest a strong association among SA, flavan-3-ol biosynthesis, and rust resistance in poplars. Abscisic acid also promoted poplar defense against rust infection, but likely through stomatal immunity independent of flavan-3-ols. Jasmonic acid did not confer any apparent defense responses to the fungal pathogen. We conclude that SA activates flavan-3-ol biosynthesis in poplar against rust infection.

Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation.

The genus Liriodendron belongs to the family Magnoliaceae, which resides within the magnoliids, an early diverging lineage of the Mesangiospermae. However, the phylogenetic relationship of magnoliids with eudicots and monocots has not been conclusively resolved and thus remains to be determined1-6. Liriodendron is a relict lineage from the Tertiary with two distinct species-one East Asian (L. chinense (Hemsley) Sargent) and one eastern North American (L. tulipifera Linn)-identified as a vicariad species pair. However, the genetic divergence and evolutionary trajectories of these species remain to be elucidated at the whole-genome level7. Here, we report the first de novo genome assembly of a plant in the Magnoliaceae, L. chinense. Phylogenetic analyses suggest that magnoliids are sister to the clade consisting of eudicots and monocots, with rapid diversification occurring in the common ancestor of these three lineages. Analyses of population genetic structure indicate that L. chinense has diverged into two lineages-the eastern and western groups-in China. While L. tulipifera in North America is genetically positioned between the two L. chinense groups, it is closer to the eastern group. This result is consistent with phenotypic observations that suggest that the eastern and western groups of China may have diverged long ago, possibly before the intercontinental differentiation between L. chinense and L. tulipifera. Genetic diversity analyses show that L. chinense has tenfold higher genetic diversity than L. tulipifera, suggesting that the complicated regions comprising east-west-orientated mountains and the Yangtze river basin (especially near 30° N latitude) in East Asia offered more successful refugia than the south-north-orientated mountain valleys in eastern North America during the Quaternary glacial period.

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.