Arabidopsis TETRASPANIN8 mediates exosome secretion and glycosyl inositol phosphoceramide sorting and trafficking.

Sphingolipids are components of plant membranes, and their heterogeneous distribution gives different membrane systems distinct properties. For example, glycosyl inositol phosphoceramides (GIPCs), 1 major type of sphingolipids, aggregate in the outer layer of the plasma membrane (PM), as well as in extracellular vesicles (EVs), including the small (30 to 100 nm) EVs termed exosomes. How these sphingolipids are sorted and trafficked is not clear. In this work, we report that Arabidopsis thaliana TETRASPANIN8 (TET8) acts as a sphingolipid carrier and thus regulates the export of GIPCs from the Golgi apparatus. TET8 recognized the coat protein complex I (COPI) subunit γ2-COPI and moved to its proper location in the PM; this recognition required the TET8 C-terminal tail. Deleting the C-terminal tail of TET8 largely restricted its roles in GIPC transport and endosomal trafficking. Further, we show that TET8 affects EV secretion in association with GIPCs. Thus, our findings shed light on GIPC transport and the molecular machinery involved in EV biogenesis.

Comprehensive identification and expression analysis of CRY gene family in Gossypium.

BACKGROUND: The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance. RESULTS: In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments. CONCLUSION: We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

Aromatization of natural products by a specialized detoxification enzyme.

In plants, lineage-specific metabolites can be created by activities derived from the catalytic promiscuity of ancestral proteins, although examples of recruiting detoxification systems to biosynthetic pathways are scarce. The ubiquitous glyoxalase (GLX) system scavenges the cytotoxic methylglyoxal, in which GLXI isomerizes the α-hydroxy carbonyl in the methylglyoxal-glutathione adduct for subsequent hydrolysis. We show that GLXIs across kingdoms are more promiscuous than recognized previously and can act as aromatases without cofactors. In cotton, a specialized GLXI variant, SPG, has lost its GSH-binding sites and organelle-targeting signal, and evolved to aromatize cyclic sesquiterpenes bearing α-hydroxyketones to synthesize defense compounds in the cytosol. Notably, SPG is able to transform acetylated deoxynivalenol, the prevalent mycotoxin contaminating cereals and foods. We propose that detoxification enzymes are a valuable source of new catalytic functions and SPG, a standalone enzyme catalyzing complex reactions, has potential for toxin degradation, crop engineering and design of novel aromatics.

Comparative Genomics and Functional Studies of Putative m6A Methyltransferase (METTL) Genes in Cotton.

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Comparative genomic analysis of head and body/tail of pancreatic ductal adenocarcinoma at early and late stages.

Pancreatic ductal adenocarcinoma (PDAC), one of the most lethal human cancers, can be divided into head and body/tail cancers anatomically. We previously reported a prognostic relevance of tumour location in resectable PDAC. This study aimed to further explore the mechanism underlying the molecular diversity between the head and body/tail of PDACs. We detected tumour genomes in 154 resectable (surgery) and non-resectable (biopsy) PDACs using a next-generation sequencing panel. Wilcoxon's rank test or Fisher's exact test was used for evaluating associations between clinical characteristics, mutation frequency and survival probability between the two cohorts. Compared with pancreatic head cancers, pancreatic body/tail cancers showed significantly more enriched genomic alterations in KRAS (97.1% vs 82.4%, P = 0.004) and SMAD4 (42.0% vs 21.2%, P = 0.008). At early stages (I-II), the SMAD4 mutation rate was significantly higher in pancreatic body/tail cancers than pancreatic head cancers (56.0% vs 26.5%, P = 0.021). At late stages (III-IV), pancreatic body/tail cancers presented significantly higher KRAS mutation rate (100.0% vs 75.8%, P = 0.001), higher frequency of MAPK pathway mutation (100% vs 87.8%, P = 0.040) and lower rates of druggable genomic alterations (30.8% vs 57.6%, P = 0.030) than pancreatic head cancers. Our work points out that pancreatic body/tail cancer seems to be more malignant than pancreatic head cancer at late stages.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development.

BACKGROUND: Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored. RESULTS: Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes. CONCLUSION: Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.

Mediation of JA signalling in glandular trichomes by the woolly/SlMYC1 regulatory module improves pest resistance in tomato.

Almost all plants form trichomes, which protect them against insect herbivores by forming a physical barrier and releasing chemical repellents. Glandular trichomes produce a variety of specialized defensive metabolites, including volatile terpenes. Previous studies have shown that the defence hormone jasmonic acid (JA) affects trichome development and induces terpene synthases (TPSs) but the underlying molecular mechanisms remain unclear. Here, we characterized a loss-of-function allele of the HD-ZIP IV transcription factor woolly (wo) and analysed its role in mediating JA signalling in tomato. We showed that knockout of wo led to extensive trichome defects, including structural and functional changes in type VI glandular trichomes, and a dramatic reduction in terpene levels. We further found that wo directly binds to TPS gene promoters to recruit SlMYC1, a JA signalling modulator, and that together these transcription factors promote terpene biosynthesis in tomato trichomes. The wo/SlMYC1 regulatory module is inhibited by SlJAZ2 through a competitive binding mechanism, resulting in a fine-tuned JA response in tomato trichomes. Enhanced expression of SlMYC1 substantially increased terpene levels and improved tomato resistance to spider mites. Interestingly, we also found that SlMYC1 plays an additional role in glandular cell division and expansion in type VI trichomes, independent of JA. Together, our results reveal a novel, JA-mediated regulatory mechanism that promotes insect resistance in tomato.

1,10/1,11-Cyclization catalyzed by diverged plant sesquiterpene synthases is dependent on a single residue.

The exquisite chemodiversity of terpenoids is the product of the large diverse terpene synthase (TPS) superfamily. Here, by using structural and phylogenetic analyses and site-directed mutagenesis, we identified a residue (Cys440 in Nicotiana tabacum 5-epi-aristolochene synthase) proximal to an ion-binding motif common to all TPSs and named the preNSE/DTE residue, which determines the product specificity of sesquiterpene synthases from different plant species. In sesquiterpene synthases catalyzing 1,10-cyclization (1,10-cyclases) of farnesyl diphosphate, mutation of the residue in both specific and promiscuous 1,10-cyclases from different lineages leads to the accumulation of monocyclic germacrene A-11-ol, which is "short-circuited" from complex cyclization cascades, suggesting a key role of this residue in generating the first common intermediate of 1,10-cyclization. Altering this residue in a specific 1,11-cyclase results in alternative 1,10-cyclization products. Moreover, the preNSE/DTE residue can be harnessed to engineer highly specific sesquiterpene synthases for an improved proportion of high-value terpenoids, such as patchoulol, a main constituent of several traditional Chinese medicines that could treat SARS-CoV-2.

Characterization of gossypol biosynthetic pathway.

Gossypol and related sesquiterpene aldehydes in cotton function as defense compounds but are antinutritional in cottonseed products. By transcriptome comparison and coexpression analyses, we identified 146 candidates linked to gossypol biosynthesis. Analysis of metabolites accumulated in plants subjected to virus-induced gene silencing (VIGS) led to the identification of four enzymes and their supposed substrates. In vitro enzymatic assay and reconstitution in tobacco leaves elucidated a series of oxidative reactions of the gossypol biosynthesis pathway. The four functionally characterized enzymes, together with (+)-δ-cadinene synthase and the P450 involved in 7-hydroxy-(+)-δ-cadinene formation, convert farnesyl diphosphate (FPP) to hemigossypol, with two gaps left that each involves aromatization. Of six intermediates identified from the VIGS-treated leaves, 8-hydroxy-7-keto-δ-cadinene exerted a deleterious effect in dampening plant disease resistance if accumulated. Notably, CYP71BE79, the enzyme responsible for converting this phytotoxic intermediate, exhibited the highest catalytic activity among the five enzymes of the pathway assayed. In addition, despite their dispersed distribution in the cotton genome, all of the enzyme genes identified show a tight correlation of expression. Our data suggest that the enzymatic steps in the gossypol pathway are highly coordinated to ensure efficient substrate conversion.

Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium.

BACKGROUND: Cotton (Gossypium spp.) is the most important world-wide fiber crop but salt stress limits cotton production in coastal and other areas. Growth regulation factors (GRFs) play regulatory roles in response to salt stress, but their roles have not been studied in cotton under salt stress. RESULTS: We identified 19 GRF genes in G. raimondii, 18 in G. arboreum, 34 in G. hirsutum and 45 in G. barbadense, respectively. These GRF genes were phylogenetically analyzed leading to the recognition of seven GRF clades. GRF genes from diploid cottons (G. raimondii and G. arboreum) were largely retained in allopolyploid cotton, with subsequent gene expansion in G. barbadense relative to G. hirsutum. Most G. hirsutum GRF (GhGRF) genes are preferentially expressed in young and growing tissues. To explore their possible role in salt stress, we used qRT-PCR to study expression responses to NaCl treatment, showing that five GhGRF genes were down-regulated in leaves. RNA-seq experiments showed that seven GhGRF genes exhibited decreased expression in leaves under NaCl treatment, three of which (GhGRF3, GhGRF4, and GhGRF16) were identified by both RNA-seq and qRT-PCR. We also identified six and three GRF genes that exhibit decreased expression under salt stress in G. arboreum and G. barbadense, respectively. Consistent with its lack of leaf withering or yellowing under the salt treatment conditions, G. arboreum had better salt tolerance than G. hirsutum and G. barbadense. Our results suggest that GRF genes are involved in salt stress responses in Gossypium. CONCLUSION: In summary, we identified candidate GRF genes that were involved in salt stress responses in cotton.

Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation.

Cotton cultivars have evolved to produce extensive, long, seed-born fibers important for the textile industry, but we know little about the molecular mechanism underlying spinnable fiber formation. Here, we report how PACLOBUTRAZOL RESISTANCE 1 (PRE1) in cotton, which encodes a basic helix-loop-helix (bHLH) transcription factor, is a target gene of spinnable fiber evolution. Differential expression of homoeologous genes in polyploids is thought to be important to plant adaptation and novel phenotypes. PRE1 expression is specific to cotton fiber cells, upregulated during their rapid elongation stage and A-homoeologous biased in allotetraploid cultivars. Transgenic studies demonstrated that PRE1 is a positive regulator of fiber elongation. We determined that the natural variation of the canonical TATA-box, a regulatory element commonly found in many eukaryotic core promoters, is necessary for subgenome-biased PRE1 expression, representing a mechanism underlying the selection of homoeologous genes. Thus, variations in the promoter of the cell elongation regulator gene PRE1 have contributed to spinnable fiber formation in cotton. Overexpression of GhPRE1 in transgenic cotton yields longer fibers with improved quality parameters, indicating that this bHLH gene is useful for improving cotton fiber quality.

Successful response to the combination of immunotherapy and chemotherapy in cholangiocarcinoma with high tumour mutational burden and PD-L1 expression: a case report.

BACKGROUND: Cholangiocarcinoma, or bile duct cancer, is a gastrointestinal cancer with limited therapeutic options and a poor outcome. Studies have revealed that some major driver genes are associated with cholangiocarcinoma, but no targeted therapies have been approved. Immune checkpoint inhibitors, which are represented by inhibitors of programmed cell death 1 (PD-1)/programmed death-ligand 1 (PD-L1), have emerged as a potential therapy for multiple types of solid cancers. CASE PRESENTATION: A 53-year-old female presented with postoperative recurrence of PD-L1-positive intrahepatic cholangiocarcinoma with a high tumour mutational burden. This patient exhibited a marked response to the combination of anti-PD-1 immunotherapy and chemotherapy. CONCLUSIONS: As far as we know, this is the first case report on the success of the combination of immunotherapy and chemotherapy for advanced cholangiocarcinoma with PD-L1 positivity and a high tumour mutational burden.

Interaction between two timing microRNAs controls trichome distribution in Arabidopsis.

The miR156-targeted squamosa promoter binding protein like (SPL) transcription factors function as an endogenous age cue in regulating plant phase transition and phase-dependent morphogenesis, but the control of SPL output remains poorly understood. In Arabidopsis thaliana the spatial pattern of trichome is a hallmark of phase transition and governed by SPLs. Here, by dissecting the regulatory network controlling trichome formation on stem, we show that the miR171-targeted lost meristems 1 (LOM1), LOM2 and LOM3, encoding GRAS family members previously known to maintain meristem cell polarity, are involved in regulating the SPL activity. Reduced LOM abundance by overexpression of miR171 led to decreased trichome density on stems and floral organs, and conversely, constitutive expression of the miR171-resistant LOM (rLOM) genes promoted trichome production, indicating that LOMs enhance trichome initiation at reproductive stage. Genetic analysis demonstrated LOMs shaping trichome distribution is dependent on SPLs, which positively regulate trichome repressor genes TRICHOMELESS 1 (TCL1) and TRIPTYCHON (TRY). Physical interaction between the N-terminus of LOMs and SPLs underpins the repression of SPL activity. Importantly, other growth and developmental events, such as flowering, are also modulated by LOM-SPL interaction, indicating a broad effect of the LOM-SPL interplay. Furthermore, we provide evidence that MIR171 gene expression is regulated by its targeted LOMs, forming a homeostatic feedback loop. Our data uncover an antagonistic interplay between the two timing miRNAs in controlling plant growth, phase transition and morphogenesis through direct interaction of their targets.

Targeting insect mitochondrial complex I for plant protection.

Plant engineered to express double-stranded RNA (dsRNA) can target the herbivorous insect gene for silencing. Although mounting evidence has emerged to support feasibility of this new pest control technology, field application is slow largely due to lack of potent targets. Here, we show that suppression of the gene encoding NDUFV2, a subunit of mitochondrial complex I that catalyses NADH dehydrogenation in respiratory chain, was highly lethal to insects. Feeding cotton bollworm (Helicoverpa armigera) larvae with transgenic cotton tissues expressing NDUFV2 dsRNA led to mortality up to 80% within 5 days, and almost no larvae survived after 7 days of feeding, due to the altered mitochondrial structure and activity. Transcriptome comparisons showed a drastic repression of dopa decarboxylase genes. Reciprocal assays with Asian corn borer (Ostrinia furnacalis), another lepidopteran species, revealed the sequence-specific effect of NDUFV2 suppression. Furthermore, the hemipteran lugus Apolygus lucorum was also liable to NDUFV2 repression. These data demonstrate that the mitochondrial complex I is a promising target with both sequence specificity and wide applicability for the development of new-generation insect-proof crops.

Functional characterization of a basic helix-loop-helix (bHLH) transcription factor GhDEL65 from cotton (Gossypium hirsutum).

Cotton fiber is proposed to share some similarity with the Arabidopsis thaliana leaf trichome, which is regulated by the MYB-bHLH-WD40 transcription complex. Although several MYB transcription factors and WD40 family proteins in cotton have been characterized, little is known about the role of bHLH family proteins in cotton. Here, we report that GhDEL65, a bHLH protein from cotton (Gossypium hirsutum), is a functional homologue of Arabidopsis GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) in regulating trichome development. Transcripts of GhDEL65 were detected in 0 ∼ 1 days post-anthesis (DPA) ovules and abundant in 3-DPA fibers, implying that GhDEL65 may act in early fiber development. Ectopic expression of GhDEL65 in Arabidopsis gl3 egl3 double mutant partly rescued the trichome development, and constitutive expression of GhDEL65 in wild-type plants led to increased trichome density on rosette leaves and stems, mainly by activating the transcription of two key positive regulators of trichome development, GLABRA1 (GL1) and GLABRA2 (GL2), and suppressed the expression of a R3 single-repeat MYB factor TRIPTYCHON (TRY). GhDEL65 could interact with cotton R2R3 MYB transcription factors GhMYB2 and GhMYB3, as well as the WD40 protein GhTTG3, suggesting that the MYB-bHLH-WD40 protein complex also exists in cotton fiber cell, though its function in cotton fiber development awaits further investigation.

Computational study of effect of water finger on ion transport through water-oil interface.

When an ion transports from water to oil through water-oil interface, it accompanies hydrated water molecules and transiently forms a chain of water, called "water finger." We thoroughly investigated the role of the water finger in chloride ion transport through water-dichloromethane interface by using molecular dynamics technique. We developed a proper coordinate w to describe the water finger structure and calculated the free energy landscape and the friction for the ion transport as a function of ion position z and the water finger coordinate w. It is clearly shown that the formation and break of water finger accompanies an activation barrier for the ion transport, which has been overlooked in the conventional free energy curve along the ion position z. The present analysis of the friction does not support the hypothesis of augmented local friction (reduced local diffusion coefficient) at the interface. These results mean that the experimentally observed rate constants of interfacial ion transfer are reduced from the diffusion-limited one because of the activation barrier associated to the water finger, not the anomalous local diffusion. We also found that the nascent ion just after the break of water finger has excessive hydration water than that in the oil phase.

Microscopic Barrier Mechanism of Ion Transport through Liquid-Liquid Interface.

Microscopic mechanism of ion transport through water-oil interface was investigated with molecular dynamics simulation. The formation/breaking of a water finger during the ion passage was explicitly formulated in the free energy surface. The calculated 2D free energy surface clearly revealed a hidden barrier of ion passage accompanied by the water finger. This barrier elucidates the retarded rate of interfacial ion transfer.

Arabidopsis MYC2 interacts with DELLA proteins in regulating sesquiterpene synthase gene expression.

Arabidopsis thaliana flowers emit volatile terpenes, which may function in plant-insect interactions. Here, we report that Arabidopsis MYC2, a basic helix-loop-helix transcription factor, directly binds to promoters of the sesquiterpene synthase genes TPS21 and TPS11 and activates their expression. Expression of TPS21 and TPS11 can be induced by the phytohormones gibberellin (GA) and jasmonate (JA), and both inductions require MYC2. The induction of TPS21 and TPS11 results in increased emission of sesquiterpene, especially (E)-ß-caryophyllene. DELLAs, the GA signaling repressors, negatively affect sesquiterpene biosynthesis, as the sesquiterpene synthase genes were repressed in plants overaccumulating REPRESSOR OF GA1-3 (RGA), one of the Arabidopsis DELLAs, and upregulated in a penta DELLA-deficient mutant. Yeast two-hybrid and coimmunoprecipitation assays demonstrated that DELLAs, represented by RGA, directly interact with MYC2. In yeast cells, the N terminus of MYC2 was responsible for binding to RGA. MYC2 has been proposed as a major mediator of JA signaling and crosstalk with abscisic acid, ethylene, and light signaling pathways. Our results demonstrate that MYC2 is also connected to GA signaling in regulating a subset of genes. In Arabidopsis inflorescences, it integrates both GA and JA signals into transcriptional regulation of sesquiterpene synthase genes and promotes sesquiterpene production.

Rational engineering of plasticity residues of sesquiterpene synthases from Artemisia annua: product specificity and catalytic efficiency.

Most TPSs (terpene synthases) contain plasticity residues that are responsible for diversified terpene products and functional evolution, which provide a potential for improving catalytic efficiency. Artemisinin, a sesquiterpene lactone from Artemisia annua L., is widely used for malaria treatment and progress has been made in engineering the production of artemisinin or its precursors. In the present paper, we report a new sesquiterpene synthase from A. annua, AaBOS (A. annua α-bisabolol synthase), which has high sequence identity with AaADS (A. annua amorpha-4,11-diene synthase), a key enzyme in artemisinin biosynthesis. Comparative analysis of the two enzymes by domain-swapping and structure-based mutagenesis led to the identification of several plasticity residues, whose alteration changed the product profile of AaBOS to include γ-humulene as the major product. To elucidate the underlying mechanisms, we solved the crystal structures of AaBOS and a γ-humulene-producing AaBOS mutant (termed AaBOS-M2). Among the plasticity residues, position 399, located in the substrate-binding pocket, is crucial for both enzymes. In AaBOS, substitution of threonine for leucine (AaBOSL339T) is required for γ-humulene production; whereas in AaADS, replacing the threonine residue with serine (AaADST399S) resulted in a substantial increase in the activity of amorpha-4,11-diene production, probably as a result of accelerated product release. The present study demonstrates that substitution of plasticity residues has potential for improving catalytic efficiency of the enzyme.

miRNA control of vegetative phase change in trees.

After germination, plants enter juvenile vegetative phase and then transition to an adult vegetative phase before producing reproductive structures. The character and timing of the juvenile-to-adult transition vary widely between species. In annual plants, this transition occurs soon after germination and usually involves relatively minor morphological changes, whereas in trees and other perennial woody plants it occurs after months or years and can involve major changes in shoot architecture. Whether this transition is controlled by the same mechanism in annual and perennial plants is unknown. In the annual forb Arabidopsis thaliana and in maize (Zea mays), vegetative phase change is controlled by the sequential activity of microRNAs miR156 and miR172. miR156 is highly abundant in seedlings and decreases during the juvenile-to-adult transition, while miR172 has an opposite expression pattern. We observed similar changes in the expression of these genes in woody species with highly differentiated, well-characterized juvenile and adult phases (Acacia confusa, Acacia colei, Eucalyptus globulus, Hedera helix, Quercus acutissima), as well as in the tree Populus x canadensis, where vegetative phase change is marked by relatively minor changes in leaf morphology and internode length. Overexpression of miR156 in transgenic P. x canadensis reduced the expression of miR156-targeted SPL genes and miR172, and it drastically prolonged the juvenile phase. Our results indicate that miR156 is an evolutionarily conserved regulator of vegetative phase change in both annual herbaceous plants and perennial trees.