Ying-Hsuan Sun

Novel and Mechanical Stress–Responsive MicroRNAs in Populus trichocarpa That Are Absent from Arabidopsis

MicroRNAs (miRNAs) are small, noncoding RNAs that can play crucial regulatory roles in eukaryotes by targeting mRNAs for silencing. To test whether miRNAs play roles in the regulation of wood development in tree species, we isolated small RNAs from the developing xylem of Populus trichocarpa stems and cloned 22 miRNAs. They are the founding members of 21 miRNA gene families for 48 miRNA sequences, represented by 98 loci in the Populus genome. A majority of these miRNAs were predicted to target developmental- and stress/defense-related genes and possible functions associated with the biosynthesis of cell wall metabolites. Of the 21 P. trichocarpa miRNA families, 11 have sequence conservation in Arabidopsis thaliana but exhibited species-specific developmental expression patterns, suggesting that even conserved miRNAs may have different regulatory roles in different species. Most unexpectedly, the remaining 10 miRNAs, for which 17 predicted targets were experimentally validated in vivo, are absent from the Arabidopsis genome, suggesting possible roles in tree-specific processes. In fact, the expression of a majority of the cloned miRNAs was upregulated or downregulated in woody stems in a manner consistent with tree-specific corrective growth against tension and compression stresses, two constant mechanical loads in trees. Our results show that plant miRNAs can be induced by mechanical stress and may function in one of the most critical defense systems for structural and mechanical fitness.

Stress-responsive microRNAs in Populus

MicroRNAs (miRNAs), a group of small non-coding RNAs, have recently become the subject of intense study. They are a class of post-transcriptional negative regulators playing vital roles in plant development and growth. However, little is known about their regulatory roles in the responses of trees to the stressful environments incurred over their long-term growth. Here, we report the cloning of small RNAs from abiotic stressed tissues of Populus trichocarpa (Ptc) and the identification of 68 putative miRNA sequences that can be classified into 27 families based on sequence homology. Among them, nine families are novel, increasing the number of the known Ptc-miRNA families from 33 to 42. A total of 346 targets was predicted for the cloned Ptc-miRNAs using penalty scores of </=2.5 for mismatched patterns in the miRNA:mRNA duplexes as the criterion. Six of the selected targets were validated experimentally. The expression of a majority of the novel miRNAs was altered in response to cold, heat, salt, dehydration, and mechanical stresses. Microarray analysis of known Ptc-miRNAs identified 19 additional cold stress-responsive Ptc-miRNAs from 14 miRNA gene families. Interestingly, we found that individual miRNAs of a family responded differentially to stress, which suggests that the members of a family may have different functions. These results reveal possible roles for miRNAs in the regulatory networks associated with the long-term growth of tree species and provide useful information for developing trees with a greater level of stress resistance.

Towards a Systems Approach for Lignin Biosynthesis in Populus trichocarpa: Transcript Abundance and Specificity of the Monolignol Biosynthetic Genes

As a step toward a comprehensive description of lignin biosynthesis in Populus trichocarpa, we identified from the genome sequence 95 phenylpropanoid gene models in 10 protein families encoding enzymes for monolignol biosynthesis. Transcript abundance was determined for all 95 genes in xylem, leaf, shoot and phloem using quantitative real-time PCR (qRT-PCR). We identified 23 genes that most probably encode monolignol biosynthesis enzymes during wood formation. Transcripts for 18 of the 23 are abundant and specific to differentiating xylem. We found evidence suggesting functional redundancy at the transcript level for phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate:CoA ligase (4CL), p-hydroxycinnamoyl-CoA:quinate shikimate p-hydroxycinnamoyltransferase (HCT), caffeoyl-CoA O-methyltransferase (CCoAOMT) and coniferyl aldehyde 5-hydroxylase (CAld5H). We carried out an enumeration-based motif identification and discriminant analysis on the promoters of all 95 genes. Five core motifs correctly discriminate the 18 xylem-specific genes from the 77 non-xylem genes. These motifs are similar to promoter elements known to regulate phenylpropanoid gene expression. This work suggests that genes in monolignol biosynthesis are regulated by multiple motifs, often related in sequence.

The Cellulose Synthase Gene Superfamily and Biochemical Functions of Xylem-Specific Cellulose Synthase-Like Genes in Populus trichocarpa

Wood from forest trees modified for more cellulose or hemicelluloses could be a major feedstock for fuel ethanol. Xylan and glucomannan are the two major hemicelluloses in wood of angiosperms. However, little is known about the genes and gene products involved in the synthesis of these wood polysaccharides. Using Populus trichocarpa as a model angiosperm tree, we report here a systematic analysis in various tissues of the absolute transcript copy numbers of cellulose synthase superfamily genes, the cellulose synthase (CesA) and the hemicellulose-related cellulose synthase-like (Csl) genes. Candidate Csl genes were characterized for biochemical functions in Drosophila Schneider 2 (S2) cells. Of the 48 identified members, 37 were found expressed in various tissues. Seven CesA genes are xylem specific, suggesting gene networks for the synthesis of wood cellulose. Four Csl genes are xylem specific, three of which belong to the CslA subfamily. The more xylem-specific CslA subfamily is represented by three types of members: PtCslA1, PtCslA3, and PtCslA5. They share high sequence homology, but their recombinant proteins produced by the S2 cells exhibited distinct substrate specificity. PtCslA5 had no catalytic activity with the substrates for xylan or glucomannan. PtCslA1 and PtCslA3 encoded mannan synthases, but PtCslA1 further encoded a glucomannan synthase for the synthesis of (1→4)-β-d-glucomannan. The expression of PtCslA1 is most highly xylem specific, suggesting a key role for it in the synthesis of wood glucomannan. The results may help guide further studies to learn about the regulation of cellulose and hemicellulose synthesis in wood.

Functional genomics and cell wall biosynthesis in loblolly pine

Loblolly pine (Pinus taeda L.) is the most widely planted tree species in the USA and an important tree in commercial forestry world-wide. The large genome size and long generation time of this species present obstacles to both breeding and molecular genetic analysis. Gene discovery by partial DNA sequence determination of cDNA clones is an effective means of building a knowledge base for molecular investigations of mechanisms governing aspects of pine growth and development, including the commercially relevant properties of secondary cell walls in wood. Microarray experiments utilizing pine cDNA clones can be used to gain additional information about the potential roles of expressed genes in wood formation. Different methods have been used to analyze data from first-generation pine microarrays, with differing degrees of success. Disparities in predictions of differential gene expression between cDNA sequencing experiments and microarray experiments arise from differences in the nature of the respective analyses, but both approaches provide lists of candidate genes which should be further investigated for potential roles in cell wall formation in differentiating pine secondary xylem. Some of these genes seem to be specific to pine, while others also occur in model plants such as Arabidopsis, where they could be more efficiently investigated.

Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa

Laccases, as early as 1959, were proposed to catalyze the oxidative polymerization of monolignols. Genetic evidence in support of this hypothesis has been elusive due to functional redundancy of laccase genes. An Arabidopsis double mutant demonstrated the involvement of laccases in lignin biosynthesis. We previously identified a subset of laccase genes to be targets of a microRNA (miRNA) ptr-miR397a in Populus trichocarpa. To elucidate the roles of ptr-miR397a and its targets, we characterized the laccase gene family and identified 49 laccase gene models, of which 29 were predicted to be targets of ptr-miR397a. We overexpressed Ptr-MIR397a in transgenic P. trichocarpa. In each of all nine transgenic lines tested, 17 PtrLACs were down-regulated as analyzed by RNA-seq. Transgenic lines with severe reduction in the expression of these laccase genes resulted in an ∼40% decrease in the total laccase activity. Overexpression of Ptr-MIR397a in these transgenic lines also reduced lignin content, whereas levels of all monolignol biosynthetic gene transcripts remained unchanged. A hierarchical genetic regulatory network (GRN) built by a bottom-up graphic Gaussian model algorithm provides additional support for a role of ptr-miR397a as a negative regulator of laccases for lignin biosynthesis. Full transcriptome–based differential gene expression in the overexpressed transgenics and protein domain analyses implicate previously unidentified transcription factors and their targets in an extended hierarchical GRN including ptr-miR397a and laccases that coregulate lignin biosynthesis in wood formation. Ptr-miR397a, laccases, and other regulatory components of this network may provide additional strategies for genetic manipulation of lignin content.

Adenylation of plant miRNAs

The modification or degradation of RNAs including miRNAs may play vital roles in regulating RNA functions. The polyadenylation- and exosome-mediated RNA decay is involved in the degradation of plant RNAs including the primary miRNA processing intermediates. However, plant miRNA levels are not affected by exosome depletion. Here, we report the cloning of a large number of 5′ and/or 3′ truncated versions of the known miRNAs from various tissues of Populus trichocarpa (black cottonwood). It suggests that plant miRNAs may be degraded through either 5′ to 3′ or 3′ to 5′ exonucleolytic digestion. We also show that a significant portion of the isolated miRNAs contains, at the 3′-end, one or a few post-transcriptionally added adenylic acid residues, which are distinct in length from the polyadenylate tail added to other plant RNAs for exosome-mediated degradation. Using an in vitro miRNA degradation system, where synthesized miRNA oligos were degraded in extracts of P. trichocarpa cells, we revealed that the adenylated miRNAs were degraded slower than others without adenylation. It indicates that addition of adenylic acid residues on the 3′-end plays a negative role in miRNA degradation. Our results provide new information for understanding the mechanism of miRNA degradation.

Complete Proteomic-Based Enzyme Reaction and Inhibition Kinetics Reveal How Monolignol Biosynthetic Enzyme Families Affect Metabolic Flux and Lignin in Populus trichocarpa

A proteomic-based predictive kinetic metabolic-flux model was developed for monolignol biosynthesis in Populus trichocarpa. Absolute quantities of all monolignol pathway proteins and 189 kinetic parameters were generated to construct the model, which was experimentally validated in transgenic P. trichocarpa and provides a comprehensive description of the monolignol biosynthetic pathway. We established a predictive kinetic metabolic-flux model for the 21 enzymes and 24 metabolites of the monolignol biosynthetic pathway using Populus trichocarpa secondary differentiating xylem. To establish this model, a comprehensive study was performed to obtain the reaction and inhibition kinetic parameters of all 21 enzymes based on functional recombinant proteins. A total of 104 Michaelis-Menten kinetic parameters and 85 inhibition kinetic parameters were derived from these enzymes. Through mass spectrometry, we obtained the absolute quantities of all 21 pathway enzymes in the secondary differentiating xylem. This extensive experimental data set, generated from a single tissue specialized in wood formation, was used to construct the predictive kinetic metabolic-flux model to provide a comprehensive mathematical description of the monolignol biosynthetic pathway. The model was validated using experimental data from transgenic P. trichocarpa plants. The model predicts how pathway enzymes affect lignin content and composition, explains a long-standing paradox regarding the regulation of monolignol subunit ratios in lignin, and reveals novel mechanisms involved in the regulation of lignin biosynthesis. This model provides an explanation of the effects of genetic and transgenic perturbations of the monolignol biosynthetic pathway in flowering plants.

Splice variant of the SND1 transcription factor is a dominant negative of SND1 members and their regulation in Populus trichocarpa

Secondary Wall-Associated NAC Domain 1s (SND1s) are transcription factors (TFs) known to activate a cascade of TF and pathway genes affecting secondary cell wall biosynthesis (xylogenesis) in Arabidopsis and poplars. Elevated SND1 transcriptional activation leads to ectopic xylogenesis and stunted growth. Nothing is known about the upstream regulators of SND1. Here we report the discovery of a stem-differentiating xylem (SDX)-specific alternative SND1 splice variant, PtrSND1-A2IR, that acts as a dominant negative of SND1 transcriptional network genes in Populus trichocarpa. PtrSND1-A2IR derives from PtrSND1-A2, one of the four fully spliced PtrSND1 gene family members (PtrSND1-A1, -A2, -B1, and -B2). Each full-size PtrSND1 activates its own gene, and all four full-size members activate a common MYB gene (PtrMYB021). PtrSND1-A2IR represses the expression of its PtrSND1 member genes and PtrMYB021. Repression of the autoregulation of a TF family by its only splice variant has not been previously reported in plants. PtrSND1-A2IR lacks DNA binding and transactivation abilities but retains dimerization capability. PtrSND1-A2IR is localized exclusively in cytoplasmic foci. In the presence of any full-size PtrSND1 member, PtrSND1-A2IR is translocated into the nucleus exclusively as a heterodimeric partner with full-size PtrSND1s. Our findings are consistent with a model in which the translocated PtrSND1-A2IR lacking DNA-binding and transactivating abilities can disrupt the function of full-size PtrSND1s, making them nonproductive through heterodimerization, and thereby modulating the SND1 transcriptional network. PtrSND1-A2IR may contribute to transcriptional homeostasis to avoid deleterious effects on xylogenesis and plant growth.

SND1 Transcription Factor–Directed Quantitative Functional Hierarchical Genetic Regulatory Network in Wood Formation in Populus trichocarpa

Novel methods were developed and demonstrated for the discovery of genetic regulatory networks in wood-forming tissues. Transfection of protoplasts from differentiating xylem with the transcription factor gene Ptr-SND1-B1 and novel computational analysis revealed a three-level hierarchical genetic regulatory network that was verified by ChIP and Ptr-SND1-B1 overexpression in transgenic plants. Wood is an essential renewable raw material for industrial products and energy. However, knowledge of the genetic regulation of wood formation is limited. We developed a genome-wide high-throughput system for the discovery and validation of specific transcription factor (TF)–directed hierarchical gene regulatory networks (hGRNs) in wood formation. This system depends on a new robust procedure for isolation and transfection of Populus trichocarpa stem differentiating xylem protoplasts. We overexpressed Secondary Wall-Associated NAC Domain 1s (Ptr-SND1-B1), a TF gene affecting wood formation, in these protoplasts and identified differentially expressed genes by RNA sequencing. Direct Ptr-SND1-B1–DNA interactions were then inferred by integration of time-course RNA sequencing data and top-down Graphical Gaussian Modeling–based algorithms. These Ptr-SND1-B1-DNA interactions were verified to function in differentiating xylem by anti-PtrSND1-B1 antibody-based chromatin immunoprecipitation (97% accuracy) and in stable transgenic P. trichocarpa (90% accuracy). In this way, we established a Ptr-SND1-B1–directed quantitative hGRN involving 76 direct targets, including eight TF and 61 enzyme-coding genes previously unidentified as targets. The network can be extended to the third layer from the second-layer TFs by computation or by overexpression of a second-layer TF to identify a new group of direct targets (third layer). This approach would allow the sequential establishment, one two-layered hGRN at a time, of all layers involved in a more comprehensive hGRN. Our approach may be particularly useful to study hGRNs in complex processes in plant species resistant to stable genetic transformation and where mutants are unavailable.