First Evidence of Thalassochory in the Ficus Genus: Seed Dispersal Using the Kuroshio Oceanic Current.

AIM: Plants distributed between southern Taiwan and the north of the Philippines are spread among numerous small islands in an area crossed by the powerful Kuroshio current. Oceanic currents can be effective seed-dispersal agents for coastal plant species. Moreover, the Luzon Strait is an area prone to tropical cyclones. The aim of this study is to look at the dispersal capability of an endangered coastal plant species, the Mearns fig (Ficus pedunculosa var. mearnsii), using both experimental and population genetics methods. LOCATION: Southern Taiwan, the Philippines, and the islands between Luzon and Taiwan Island. METHODS: This study combined two types of analysis, i.e., buoyancy experiments on syconia and double digest restriction-associated DNA sequencing (ddRAD), to analyze the population genetics of the Mearns fig. RESULTS: We first discovered that mature Mearns fig syconia could float in seawater. They have a mean float duration of 10 days to a maximum of 21 days. Germination rates varied significantly between Mearns fig seeds that had undergone different durations of flotation treatment. Population genetic analysis shows a high degree of inbreeding among various Mearns fig populations. Moreover, no isolation by distance was found between the populations and individuals. MAIN CONCLUSIONS: From our analysis of the genetic structure of the Mearns fig populations, we can clearly highlight the effect of the Kuroshio oceanic current on the seed dispersal of this fig tree. Comprehensive analysis has shown that Mearns fig seeds are still viable before the mature syconium sinks into the seawater, and so they could use the Kuroshio Current to float to the current population locations in Taiwan.

Analyzing genetic diversity in luffa and developing a Fusarium wilt-susceptible linked SNP marker through a single plant genome-wide association (sp-GWAS) study.

BACKGROUND: Luffa (Luffa spp.) is an economically important crop of the Cucurbitaceae family, commonly known as sponge gourd or vegetable gourd. It is an annual cross-pollinated crop primarily found in the subtropical and tropical regions of Asia, Australia, Africa, and the Americas. Luffa serves not only as a vegetable but also exhibits medicinal properties, including anti-inflammatory, antidiabetic, and anticancer effects. Moreover, the fiber derived from luffa finds extensive applications in various fields such as biotechnology and construction. However, luffa Fusarium wilt poses a severe threat to its production, and existing control methods have proven ineffective in terms of cost-effectiveness and environmental considerations. Therefore, there is an urgent need to develop luffa varieties resistant to Fusarium wilt. Single-plant GWAS (sp-GWAS) has been demonstrated as a promising tool for the rapid and efficient identification of quantitative trait loci (QTLs) associated with target traits, as well as closely linked molecular markers. RESULTS: In this study, a collection of 97 individuals from 73 luffa accessions including two major luffa species underwent single-plant GWAS to investigate luffa Fusarium wilt resistance. Utilizing the double digest restriction site associated DNA (ddRAD) method, a total of 8,919 high-quality single nucleotide polymorphisms (SNPs) were identified. The analysis revealed the potential for Fusarium wilt resistance in accessions from both luffa species. There are 6 QTLs identified from 3 traits, including the area under the disease progress curve (AUDPC), a putative disease-resistant QTL, was identified on the second chromosome of luffa. Within the region of linkage disequilibrium, a candidate gene homologous to LOC111009722, which encodes peroxidase 40 and is associated with disease resistance in Cucumis melo, was identified. Furthermore, to validate the applicability of the marker associated with resistance from sp-GWAS, an additional set of 21 individual luffa plants were tested, exhibiting 93.75% accuracy in detecting susceptible of luffa species L. aegyptiaca Mill. CONCLUSION: In summary, these findings give a hint of genome position that may contribute to luffa wild resistance to Fusarium and can be utilized in the future luffa wilt resistant breeding programs aimed at developing wilt-resistant varieties by using the susceptible-linked SNP marker.

Savinin Triggers Programmed Cell Death of Ray Parenchyma Cells in Heartwood Formation of Taiwania cryptomerioides Hayata.

The purpose of this study was to investigate the relationship between lignan biosynthesis and programmed cell death (PCD) of ray parenchyma cells during the heartwood formation of Taiwania (Taiwania cryptomerioides Hayata). Since the PCD of ray parenchyma cells and the synthesis of lignans are the two main processes involved in the formation of heartwood, both of which need to be completed through gene regulation. Based on the results of genomics and bioinformatics analysis, that the PCD of tracheids are induced by genotoxic, and the PCD of ray parenchyma cells is induced by biological factors, such as fungi, bacteria, and viruses, which could induce oxidative stress. According to the results of time-of-flight secondary ion mass spectrometry (ToF-SIMS) analysis, lignans are produced in ray parenchyma cells, and the accumulation of savinin and its downstream lignans might be the cause of PCD in ray parenchyma cells. An in vitro experiment further confirmed that the accumulation of savinin could cause protoplasts of Taiwania's xylem to produce taiwanin A, which is the marker of heartwood formation in Taiwania. Resulting in an increase in reactive oxygen species (ROS) content, which could induce oxidative stress in ray parenchyma cells and potentially lead to PCD. Based on these findings, we conclude that accumulation of savinin could be induced PCD of ray parenchyma cells in heartwood formation in Taiwania.

Single-cell transcriptomics unveils xylem cell development and evolution.

BACKGROUND: Xylem, the most abundant tissue on Earth, is responsible for lateral growth in plants. Typical xylem has a radial system composed of ray parenchyma cells and an axial system of fusiform cells. In most angiosperms, fusiform cells comprise vessel elements for water transportation and libriform fibers for mechanical support, while both functions are performed by tracheids in other vascular plants such as gymnosperms. Little is known about the developmental programs and evolutionary relationships of these xylem cell types. RESULTS: Through both single-cell and laser capture microdissection transcriptomic profiling, we determine the developmental lineages of ray and fusiform cells in stem-differentiating xylem across four divergent woody angiosperms. Based on cross-species analyses of single-cell clusters and overlapping trajectories, we reveal highly conserved ray, yet variable fusiform, lineages across angiosperms. Core eudicots Populus trichocarpa and Eucalyptus grandis share nearly identical fusiform lineages, whereas the more basal angiosperm Liriodendron chinense has a fusiform lineage distinct from that in core eudicots. The tracheids in the basal eudicot Trochodendron aralioides, an evolutionarily reversed trait, exhibit strong transcriptomic similarity to vessel elements rather than libriform fibers. CONCLUSIONS: This evo-devo framework provides a comprehensive understanding of the formation of xylem cell lineages across multiple plant species spanning over a hundred million years of evolutionary history.

The complete plastid genome of Selaginella erythropus (Selaginellaceae), a species with distinctive giant chloroplasts.

The plastid genome of the deep-shade plant Selaginella erythropus, which has highly unusual chloroplasts, was characterized using Illumina pair-end sequencing. This plastome is 140,151 bp in length with a large single-copy region (LSC) of 56,133 bp, a small single-copy region (SSC) of 61,268 bp, and two direct repeats (DRs) of 11,375 bp. The overall GC content is 50.68%, while those of LSC, SSC, and DR are 48.96%, 50.3%, and 55.96%, respectively. The plastome contains 102 genes, including 76 protein-coding, 15 tRNA (12 tRNA species), and 8 rRNA genes (4 rRNA species). The phylogenetic analysis shows that S. erythropus is closely related to S. moellendorffii and S. doederleinii. This result is consistent with the previous phylogenetic relationship inferred from multiple plastid and nuclear loci. However, only S. erythropus has the two-zoned giant chloroplast, the bizonoplast. The plastome provides an excellent reference for understanding the unique chloroplast differentiation in Selaginellaceae.

Improving wood properties for wood utilization through multi-omics integration in lignin biosynthesis.

A multi-omics quantitative integrative analysis of lignin biosynthesis can advance the strategic engineering of wood for timber, pulp, and biofuels. Lignin is polymerized from three monomers (monolignols) produced by a grid-like pathway. The pathway in wood formation of Populus trichocarpa has at least 21 genes, encoding enzymes that mediate 37 reactions on 24 metabolites, leading to lignin and affecting wood properties. We perturb these 21 pathway genes and integrate transcriptomic, proteomic, fluxomic and phenomic data from 221 lines selected from ~2000 transgenics (6-month-old). The integrative analysis estimates how changing expression of pathway gene or gene combination affects protein abundance, metabolic-flux, metabolite concentrations, and 25 wood traits, including lignin, tree-growth, density, strength, and saccharification. The analysis then predicts improvements in any of these 25 traits individually or in combinations, through engineering expression of specific monolignol genes. The analysis may lead to greater understanding of other pathways for improved growth and adaptation.

Reciprocal cross-regulation of VND and SND multigene TF families for wood formation in Populus trichocarpa.

Secondary cell wall (SCW) biosynthesis is the biological process that generates wood, an important renewable feedstock for materials and energy. NAC domain transcription factors, particularly Vascular-Related NAC-Domain (VND) and Secondary Wall-Associated NAC Domain (SND) proteins, are known to regulate SCW differentiation. The regulation of VND and SND is important to maintain homeostasis for plants to avoid abnormal growth and development. We previously identified a splice variant, PtrSND1-A2IR , derived from PtrSND1-A2 as a dominant-negative regulator, which suppresses the transactivation of all PtrSND1 family members. PtrSND1-A2IR also suppresses the self-activation of the PtrSND1 family members except for its cognate transcription factor, PtrSND1-A2, suggesting the existence of an unknown factor needed to regulate PtrSND1-A2 Here, a splice variant, PtrVND6-C1IR , derived from PtrVND6-C1 was discovered that suppresses the protein functions of all PtrVND6 family members. PtrVND6-C1IR also suppresses the expression of all PtrSND1 members, including PtrSND1-A2, demonstrating that PtrVND6-C1IR is the previously unidentified regulator of PtrSND1-A2 We also found that PtrVND6-C1IR cannot suppress the expression of its cognate transcription factor, PtrVND6-C1PtrVND6-C1 is suppressed by PtrSND1-A2IR Both PtrVND6-C1IR and PtrSND1-A2IR cannot suppress their cognate transcription factors but can suppress all members of the other family. The results indicate that the splice variants from the PtrVND6 and PtrSND1 family may exert reciprocal cross-regulation for complete transcriptional regulation of these two families in wood formation. This reciprocal cross-regulation between families suggests a general mechanism among NAC domain proteins and likely other transcription factors, where intron-retained splice variants provide an additional level of regulation.

Accumulation of high OPDA level correlates with reduced ROS and elevated GSH benefiting white cell survival in variegated leaves.

Variegated 'Marble Queen' (Epipremnum aureum) plant has white (VMW) and green (VMG) sectors within the same leaf. The white sector cells containing undifferentiated chloroplasts are viable, but the underlying mechanism for their survival and whether these white cells would use any metabolites as signal molecules to communicate with the nucleus for maintaining their viability remain unclear. We analyzed and compared phytohormone levels with their precursors produced in chloroplasts between VMW and VMG, and further compared their transcriptomes to understand the consequences related to the observed elevated 12-oxo phytodienoic acid (OPDA), which was 9-fold higher in VMW than VMG. Transcriptomic study showed that a large group of OPDA-responsive genes (ORGs) were differentially expressed in VMW, including stress-related transcription factors and genes for reactive oxygen species (ROS) scavengers, DNA replication and repair, and protein chaperones. Induced expression of these ORGs could be verified in OPDA-treated green plants. Reduced level of ROS and higher levels of glutathione in VMW were further confirmed. Our results suggest that elevated OPDA or its related compounds are recruited by white cells as a signaling molecule(s) to up-regulate stress and scavenging activity related genes that leads to reduced ROS levels and provides survival advantages to the white cells.

Tissue and cell-type co-expression networks of transcription factors and wood component genes in Populus trichocarpa.

MAIN CONCLUSION: Co-expression networks based on transcriptomes of Populus trichocarpa major tissues and specific cell types suggest redundant control of cell wall component biosynthetic genes by transcription factors in wood formation. We analyzed the transcriptomes of five tissues (xylem, phloem, shoot, leaf, and root) and two wood forming cell types (fiber and vessel) of Populus trichocarpa to assemble gene co-expression subnetworks associated with wood formation. We identified 165 transcription factors (TFs) that showed xylem-, fiber-, and vessel-specific expression. Of these 165 TFs, 101 co-expressed (correlation coefficient, r > 0.7) with the 45 secondary cell wall cellulose, hemicellulose, and lignin biosynthetic genes. Each cell wall component gene co-expressed on average with 34 TFs, suggesting redundant control of the cell wall component gene expression. Co-expression analysis showed that the 101 TFs and the 45 cell wall component genes each has two distinct groups (groups 1 and 2), based on their co-expression patterns. The group 1 TFs (44 members) are predominantly xylem and fiber specific, and are all highly positively co-expressed with the group 1 cell wall component genes (30 members), suggesting their roles as major wood formation regulators. Group 1 TFs include a lateral organ boundary domain gene (LBD) that has the highest number of positively correlated cell wall component genes (36) and TFs (47). The group 2 TFs have 57 members, including 14 vessel-specific TFs, and are generally less correlated with the cell wall component genes. An exception is a vessel-specific basic helix-loop-helix (bHLH) gene that negatively correlates with 20 cell wall component genes, and may function as a key transcriptional suppressor. The co-expression networks revealed here suggest a well-structured transcriptional homeostasis for cell wall component biosynthesis during wood formation.

Identification of Ice Plant (Mesembryanthemum crystallinum L.) MicroRNAs Using RNA-Seq and Their Putative Roles in High Salinity Responses in Seedlings.

The halophyte Mesembryanthemum crystallinum (common or crystalline ice plant) is a useful model for studying molecular mechanisms of salt tolerance. The morphology, physiology, metabolism, and gene expression of ice plant have been studied and large-scale analyses of gene expression profiling have drawn an outline of salt tolerance in ice plant. A rapid root growth to a sudden increase in salinity was observed in ice plant seedlings. Using a fluorescent dye to detect Na(+), we found that ice plant roots respond to an increased flux of Na(+) by either secreting or storing Na(+) in specialized cells. High-throughput sequencing was used to identify small RNA profiles in 3-day-old seedlings treated with or without 200 mM NaCl. In total, 135 conserved miRNAs belonging to 21 families were found. The hairpin precursor of 19 conserved mcr-miRNAs and 12 novel mcr-miRNAs were identified. After 6 h of salt stress, the expression of most mcr-miRNAs showed decreased relative abundance, whereas the expression of their corresponding target genes showed increased mRNA relative abundance. The cognate target genes are involved in a broad range of biological processes: transcription factors that regulate growth and development, enzymes that catalyze miRNA biogenesis for the most conserved mcr-miRNA, and proteins that are involved in ion homeostasis and drought-stress responses for some novel mcr-miRNAs. Analyses of the functions of target genes revealed that cellular processes, including growth and development, metabolism, and ion transport activity are likely to be enhanced in roots under salt stress. The expression of eleven conserved miRNAs and two novel miRNAs were correlated reciprocally with predicted targets within hours after salt stress exposure. Several conserved miRNAs have been known to regulate root elongation, root apical meristem activity, and lateral root formation. Based upon the expression pattern of miRNA and target genes in combination with the observation of Na(+) distribution, ice plant likely responds to increased salinity by using Na(+) as an osmoticum for cell expansion and guard cell opening. Excessive Na(+) could either be secreted through the root epidermis or stored in specialized leaf epidermal cells. These responses are regulated in part at the miRNA-mediated post-transcriptional level.

Gibberellin deficiency is responsible for shy-flowering nature of Epipremnum aureum.

Epipremnum aureum is an extremely popular houseplant belonging to the Araceae family of angiosperms, but it does not flower either in the wild or under cultivation. We uncovered the potential causes of its shy-flowering nature by building the transcriptome using next-generation sequencing and identifying floral-related genes that are differentially expressed between vertical growth (VG, adult) and horizontal growth (HG, juvenile) plants. Transcripts of the gibberellin (GA) biosynthetic gene EaGA3ox1 and GA-responsive floral meristem identity gene EaLFY were absent in both VG and HG plants, suggesting that a deficiency of bioactive GAs may be responsible for its shy-flowering nature. This hypothesis is supported by undetectable or low levels of three bioactive GAs, and exogenous GA3 triggered flowering in both plants. Our study resolves the mystery why E. aureum fails to flower, and reveals the positive role of GAs in floral transition in perennials.

A simple improved-throughput xylem protoplast system for studying wood formation.

Isolated protoplasts serve as a transient expression system that is highly representative of stable transgenics in terms of transcriptome responses. They can also be used as a cellular system to study gene transactivation and nucleocytoplasmic protein trafficking. They are particularly useful for systems studies in which stable transgenics and mutants are unavailable. We present a protocol for the isolation and transfection of protoplasts from wood-forming tissue, the stem-differentiating xylem (SDX), in the model woody plant Populus trichocarpa. The method involves tissue preparation, digestion of SDX cell walls, protoplast isolation and DNA transfection. Our approach is markedly faster and provides better yields than previous protocols; small (milligrams)- to large (20 g)-scale SDX preparations can be achieved in ~60 s, with isolation of protoplasts and their subsequent transfection taking ~50 min. Up to ten different samples can be processed simultaneously in this time scale. Our protocol gives a high yield (~2.5 × 10(7) protoplasts per g of SDX) of protoplasts sharing 96% transcriptome identity with intact SDX.

Complete proteomic-based enzyme reaction and inhibition kinetics reveal how monolignol biosynthetic enzyme families affect metabolic flux and lignin in Populus trichocarpa.

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.

SND1 transcription factor-directed quantitative functional hierarchical genetic regulatory network in wood formation in Populus trichocarpa.

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.

Regulation of phenylalanine ammonia-lyase (PAL) gene family in wood forming tissue of Populus trichocarpa.

Phenylalanine ammonia-lyase (PAL) catalyzes the initial step of phenylpropanoid biosynthesis in plants. Five PAL genes (PtrPAL1 to 5) have been identified in Populus trichocarpa. These genes are classified into two subgroups according to their transcript sequence similarity and tissue specificity. However, the regulation of these genes and their protein functions are not well understood. In this study, enzymatic properties of each PtrPALs were characterized based on their recombinant proteins expressed in E.coli. Subcellular localizations of each PtrPALs in stem wood forming tissue were investigated and individual PtrPAL protein abundances in cytosol and membrane protein fractions were measured using protein cleavage-isotope dilution mass spectrometry (PC-IDMS). Protein/mRNA ratios of PtrPALs were further verified using RNA-Seq and gel-enhanced liquid chromatography mass spectrometry (GeLC-MS). All PtrPALs have similar catalytic properties for the deamination of L-phenylalanine, their major substrate. All PtrPALs have similar subcellular locations in stem wood forming tissue, with major amount in the cytosol (93-96 %) and less in the membrane (4-7 %). However, the protein/mRNA ratios of subgroup A (PtrPAL2, 4 and 5) are about five times that of subgroup B (PtrPAL1 and 3) in stem wood forming tissue, while all PtrPALs have similar transcript abundances. These results indicate a greater functional significance of subgroup A PtrPALs for stem wood formation, and highlight the role of gene post-transcriptional regulation.

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.

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-A2(IR), that acts as a dominant negative of SND1 transcriptional network genes in Populus trichocarpa. PtrSND1-A2(IR) 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-A2(IR) 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-A2(IR) lacks DNA binding and transactivation abilities but retains dimerization capability. PtrSND1-A2(IR) is localized exclusively in cytoplasmic foci. In the presence of any full-size PtrSND1 member, PtrSND1-A2(IR) 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-A2(IR) 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-A2(IR) may contribute to transcriptional homeostasis to avoid deleterious effects on xylogenesis and plant growth.

Functional redundancy of the two 5-hydroxylases in monolignol biosynthesis of Populus trichocarpa: LC-MS/MS based protein quantification and metabolic flux analysis.

Flowering plants have syringyl and guaiacyl subunits in lignin in contrast to the guaiacyl lignin in gymnosperms. The biosynthesis of syringyl subunits is initiated by coniferaldehyde 5-hydroxylase (CAld5H). In Populus trichocarpa there are two closely related CAld5H enzymes (PtrCAld5H1 and PtrCAld5H2) associated with lignin biosynthesis during wood formation. We used yeast recombinant PtrCAld5H1 and PtrCAld5H2 proteins to carry out Michaelis-Menten and inhibition kinetics with LC-MS/MS based absolute protein quantification. CAld5H, a monooxygenase, requires a cytochrome P450 reductase (CPR) as an electron donor. We cloned and expressed three P. trichocarpa CPRs in yeast and show that all are active with both CAld5Hs. The kinetic analysis shows both CAld5Hs have essentially the same biochemical functions. When both CAld5Hs are coexpressed in the same yeast membranes, the resulting enzyme activities are additive, suggesting functional redundancy and independence of these two enzymes. Simulated reaction flux based on Michaelis-Menten kinetics and inhibition kinetics confirmed the redundancy and independence. Subcellular localization of both CAld5Hs as sGFP fusion proteins expressed in P. trichocarpa differentiating xylem protoplasts indicate that they are endoplasmic reticulum resident proteins. These results imply that during wood formation, 5-hydroxylation in monolignol biosynthesis of P. trichocarpa requires the combined metabolic flux of these two CAld5Hs to maintain adequate biosynthesis of syringyl lignin. The combination of genetic analysis, absolute protein quantitation-based enzyme kinetics, homologous CPR specificity, SNP characterization, and ER localization provides a more rigorous basis for a comprehensive systems understanding of 5-hydroxylation in lignin biosynthesis.

Comprehensive quantification of monolignol-pathway enzymes in Populus trichocarpa by protein cleavage isotope dilution mass spectrometry.

The economic value of wood/pulp from many tree species is largely dictated by the quantity and chemical properties of lignin, which is directly related to the composition and linkages of monolignols comprising the polymer. Although much is known regarding the monolignol biosynthetic pathway, our understanding is still deficient due to the lack of quantitative information at the proteomic level. We developed an assay based on protein cleavage isotope dilution mass spectrometry (PC-IDMS) for the determination of all potential, primary enzymes involved in the biosynthesis of monolignols and the peroxidases responsible for their polymerization to form lignin in the model tree species, Populus trichocarpa. Described is the identification of quantitative surrogate peptides through shotgun analysis of native and recombinant proteins, optimization of trypsin proteolysis using fractional factorial design of experiments, and development of a liquid chromatography-selected reaction monitoring method for specific detection of all targeted peptides. Of the 25 targeted enzymes, three were undetected in the normal xylem tissues, and all but two of the detectable species showed good day-to-day precision (CV < 10%). This represents the most comprehensive assay for quantification of proteins regulating monolignol biosynthesis and will lead to a better understanding of lignin formation at a systems level.

MicroRNAs in trees.

MicroRNAs (miRNAs) are 20-24 nucleotide long molecules processed from a specific class of RNA polymerase II transcripts that mainly regulate the stability of mRNAs containing a complementary sequence by targeted degradation in plants. Many features of tree biology are regulated by miRNAs affecting development, metabolism, adaptation and evolution. MiRNAs may be modified and harnessed for controlled suppression of specific genes to learn about gene function, or for practical applications through genetic engineering. Modified (artificial) miRNAs act as dominant suppressors and are particularly useful in tree genetics because they bypass the generations of inbreeding needed for fixation of recessive mutations. The purpose of this review is to summarize the current status of information on miRNAs in trees and to guide future studies on the role of miRNAs in the biology of woody perennials and to illustrate their utility in directed genetic modification of trees.