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 treespecific 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.

SummaryMicroRNAs (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 PtcmiRNAs 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.

Amborella trichopoda is strongly supported as the single living species of the sister lineage to all other extant flowering plants, providing a unique reference for inferring the genome content and structure of the most recent common ancestor (MRCA) of living angiosperms. Sequencing the Amborella genome, we identified an ancient genome duplication predating angiosperm diversification, without evidence of subsequent, lineage-specific genome duplications. Comparisons between Amborella and other angiosperms facilitated reconstruction of the ancestral angiosperm gene content and gene order in the MRCA of core eudicots. We identify new gene families, gene duplications, and floral protein-protein interactions that first appeared in the ancestral angiosperm. Transposable elements in Amborella are ancient and highly divergent, with no recent transposon radiations. Population genomic analysis across Amborella's native range in New Caledonia reveals a recent genetic bottleneck and geographic structure with conservation implications.

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 PtrMIR397a in transgenic P. trichocarpa. In each of all nine transgenic lines tested, 17 PtrLACs were down-regulated as analyzed by RNAseq. 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.L ignin, an abundant biological polymer affecting the ecology of the terrestrial biosphere, is vital for the integrity of plant cell walls, the strength of stems, and resistance against pests and pathogens (1). Lignin is also a major barrier in the pulping and biomass-to-ethanol processes (2-4). For extracting cellulose (pulping) or for enzymatic degradation of cellulose for bioethanol, harsh chemical or physical treatments are used to reduce interactions with lignin or other cell wall components (2-4). Reducing lignin content or altering lignin structure to reduce its recalcitrance are major goals for more efficient processing.Lignin is polymerized primarily from three monolignol precursors, p-coumaryl alcohol, coniferyl alcohol, and sinapyl alcohol (1, 5). Over five decades, efforts have been made to understand the biosynthesis of the primary monolignols and to modify the quantity or composition of lignin. The polymerization of monolignols into a lignin polymer has long been thought to occur through oxidative polymerization catalyzed by either laccases or peroxidases (6). The mechanisms and specificity of the roles of the oxidative enzymes in lignin polymerization have been controversial (7).Laccases (EC. 1.10.3.2) are multicopper oxidoreductases. Plant laccase was the first enzyme sh...

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.