The HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENE15-HISTONE DEACETYLASE9 complex associates with HYPONASTIC LEAVES 1 to modulate microRNA expression in response to abscisic acid signaling.

The regulation of microRNA (miRNA) biogenesis is crucial for maintaining plant homeostasis under biotic and abiotic stress. The crosstalk between the RNA polymerase II (Pol-II) complex and the miRNA processing machinery has emerged as a central hub modulating transcription and cotranscriptional processing of primary miRNA transcripts (pri-miRNAs). However, it remains unclear how miRNA-specific transcriptional regulators recognize MIRNA loci. Here, we show that the Arabidopsis (Arabidopsis thaliana) HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENE15 (HOS15)-HISTONE DEACETYLASE9 (HDA9) complex is a conditional suppressor of miRNA biogenesis, particularly in response to abscisic acid (ABA). When treated with ABA, hos15/hda9 mutants show enhanced transcription of pri-miRNAs that is accompanied by increased processing, leading to overaccumulation of a set of mature miRNAs. Moreover, upon recognition of the nascent pri-miRNAs, the ABA-induced recruitment of the HOS15-HDA9 complex to MIRNA loci is guided by HYPONASTIC LEAVES 1 (HYL1). The HYL1-dependent recruitment of the HOS15-HDA9 complex to MIRNA loci suppresses expression of MIRNAs and processing of pri-miRNA. Most importantly, our findings indicate that nascent pri-miRNAs serve as scaffolds for recruiting transcriptional regulators, specifically to MIRNA loci. This indicates that RNA molecules can act as regulators of their own expression by causing a negative feedback loop that turns off their transcription, providing a self-buffering system.

Polymorphic inverted repeats near coding genes impact chromatin topology and phenotypic traits in Arabidopsis thaliana.

Transposons are mobile elements that are commonly silenced to protect eukaryotic genome integrity. In plants, transposable element (TE)-derived inverted repeats (IRs) are commonly found near genes, where they affect host gene expression. However, the molecular mechanisms of such regulation are unclear in most cases. Expression of these IRs is associated with production of 24-nt small RNAs, methylation of the IRs, and drastic changes in local 3D chromatin organization. Notably, many of these IRs differ between Arabidopsis thaliana accessions, causing variation in short-range chromatin interactions and gene expression. CRISPR-Cas9-mediated disruption of two IRs leads to a switch in genome topology and gene expression with phenotypic consequences. Our data show that insertion of an IR near a gene provides an anchor point for chromatin interactions that profoundly impact the activity of neighboring loci. This turns IRs into powerful evolutionary agents that can contribute to rapid adaptation.

exp2GO: Improving Prediction of Functions in the Gene Ontology With Expression Data.

The computational methods for the prediction of gene function annotations aim to automatically find associations between a gene and a set of Gene Ontology (GO) terms describing its functions. Since the hand-made curation process of novel annotations and the corresponding wet experiments validations are very time-consuming and costly procedures, there is a need for computational tools that can reliably predict likely annotations and boost the discovery of new gene functions. This work proposes a novel method for predicting annotations based on the inference of GO similarities from expression similarities. The novel method was benchmarked against other methods on several public biological datasets, obtaining the best comparative results. exp2GO effectively improved the prediction of GO annotations in comparison to state-of-the-art methods. Furthermore, the proposal was validated with a full genome case where it was capable of predicting relevant and accurate biological functions. The repository of this project withh full data and code is available at https://github.com/sinc-lab/exp2GO.

The underground life of homeodomain-leucine zipper transcription factors.

Roots are the anchorage organs of plants, responsible for water and nutrient uptake, exhibiting high plasticity. Root architecture is driven by the interactions of biomolecules, including transcription factors and hormones that are crucial players regulating root plasticity. Multiple transcription factor families are involved in root development; some, such as ARFs and LBDs, have been well characterized, whereas others remain less well investigated. In this review, we synthesize the current knowledge about the involvement of the large family of homeodomain-leucine zipper (HD-Zip) transcription factors in root development. This family is divided into four subfamilies (I-IV), mainly according to structural features, such as additional motifs aside from HD-Zip, as well as their size, gene structure, and expression patterns. We explored and analyzed public databases and the scientific literature regarding HD-Zip transcription factors in Arabidopsis and other species. Most members of the four HD-Zip subfamilies are expressed in specific cell types and several individuals from each group have assigned functions in root development. Notably, a high proportion of the studied proteins are part of intricate regulation pathways involved in primary and lateral root growth and development.

Elucidating carbohydrate metabolism in Euglena gracilis: Reverse genetics-based evaluation of genes coding for enzymes linked to paramylon accumulation.

Euglena gracilis is a eukaryotic single-celled and photosynthetic organism grouped under the kingdom Protista. This phytoflagellate can accumulate the carbon photoassimilate as a linear ß-1,3-glucan chain called paramylon. This storage polysaccharide can undergo degradation to provide glucose units to obtain ATP and reducing power both in aerobic and anaerobic growth conditions. Our group has recently characterized an essential enzyme for accumulating the polysaccharide, the UDP-glucose pyrophosphorylase (Biochimie vol 154, 2018, 176-186), which catalyzes the synthesis of UDP-glucose (the substrate for paramylon synthase). Additionally, the identification of nucleotide sequences coding for putative UDP-sugar pyrophosphorylases suggests the occurrence of an alternative source of UDP-glucose. In this study, we demonstrate the active involvement of both pyrophosphorylases in paramylon accumulation. Using techniques of single and combined knockdown of transcripts coding for these proteins, we evidenced a substantial decrease in the polysaccharide synthesis from 39 ± 7 µg/106 cells determined in the control at day 21st of growth. Thus, the paramylon accumulation in Euglena gracilis cells decreased by 60% and 30% after a single knockdown of the expression of genes coding for UDP-glucose pyrophosphorylase and UDP-sugar pyrophosphorylase, respectively. Besides, the combined knockdown of both genes resulted in a ca. 65% reduction in the level of the storage polysaccharide. Our findings indicate the existence of a physiological dependence between paramylon accumulation and the partitioning of sugar nucleotides into other metabolic routes, including the Leloir pathway's functionality in Euglena gracilis.

Extensive Analysis of miRNA Trimming and Tailing Indicates that AGO1 Has a Complex Role in miRNA Turnover.

MicroRNAs are small regulatory RNAs involved in several processes in plants ranging from development and stress responses to defense against pathogens. In order to accomplish their molecular functions, miRNAs are methylated and loaded into one ARGONAUTE (AGO) protein, commonly known as AGO1, to stabilize and protect the molecule and to assemble a functional RNA-induced silencing complex (RISC). A specific machinery controls miRNA turnover to ensure the silencing release of targeted-genes in given circumstances. The trimming and tailing of miRNAs are fundamental modifications related to their turnover and, hence, to their action. In order to gain a better understanding of these modifications, we analyzed Arabidopsis thaliana small RNA sequencing data from a diversity of mutants, related to miRNA biogenesis, action, and turnover, and from different cellular fractions and immunoprecipitations. Besides confirming the effects of known players in these pathways, we found increased trimming and tailing in miRNA biogenesis mutants. More importantly, our analysis allowed us to reveal the importance of ARGONAUTE 1 (AGO1) loading, slicing activity, and cellular localization in trimming and tailing of miRNAs.

HASTY modulates miRNA biogenesis by linking pri-miRNA transcription and processing.

Post-transcriptional gene silencing mediated by microRNAs (miRNAs) modulates numerous developmental and stress response pathways. For the last two decades, HASTY (HST), the ortholog of human EXPORTIN 5, was considered to be a candidate protein that exports plant miRNAs from the nucleus to the cytoplasm. Here, we report that HST functions in the miRNA pathway independent of its cargo-exporting activity in Arabidopsis. We found that Arabidopsis mutants with impaired HST shuttling exhibit normal subcellular distribution of miRNAs. Interestingly, protein-protein interaction and microscopy assays showed that HST directly interacts with the microprocessor core component DCL1 through its N-terminal domain. Moreover, mass spectrometry analysis revealed that HST also interacts independently of its N-terminal domain with the mediator complex subunit MED37. Further experiments revealed that HST could act as a scaffold to facilitate the recruitment of DCL1 to genomic MIRNA loci by stabilizing the DCL1-MED37 complex, which in turn promotes the transcription and proper processing of primary miRNA transcripts (pri-miRNAs). Taken together, these results suggest that HST is likely associated with the formation of the miRNA biogenesis complex at MIRNA genes, promoting the transcription and processing of pri-miRNAs rather than the direct export of processed miRNAs from the nucleus.

The Intrinsically Disordered Protein CARP9 Bridges HYL1 to AGO1 in the Nucleus to Promote MicroRNA Activity.

In plants, small RNAs are loaded into ARGONAUTE (AGO) proteins to fulfill their regulatory functions. MicroRNAs (miRNAs), one of the most abundant classes of endogenous small RNAs, are preferentially loaded into AGO1. Such loading, long believed to happen exclusively in the cytoplasm, was recently proposed to also occur in the nucleus. Here, we identified CONSTITUTIVE ALTERATIONS IN THE SMALL RNAS PATHWAYS9 (CARP9), a nuclear-localized, intrinsically disordered protein, as a factor promoting miRNA activity in Arabidopsis (Arabidopsis thaliana). Mutations in the CARP9-encoding gene led to a mild reduction of miRNAs levels, impaired gene silencing, and characteristic morphological defects, including young leaf serration and altered flowering time. Intriguingly, we found that CARP9 was able to interact with HYPONASTIC LEAVES1 (HYL1), but not with other proteins of the miRNA biogenesis machinery. In the same way, CARP9 appeared to interact with mature miRNA, but not with primary miRNA, positioning it after miRNA processing in the miRNA pathway. CARP9 was also able to interact with AGO1, promoting its interaction with HYL1 to facilitate miRNA loading in AGO1. Plants deficient in CARP9 displayed reduced levels of AGO1-loaded miRNAs, partial retention of miRNA in the nucleus, and reduced levels of AGO1. Collectively, our data suggest that CARP9 might modulate HYL1-AGO1 cross talk, acting as a scaffold for the formation of a nuclear post-primary miRNA-processing complex that includes at least HYL1, AGO1, and HEAT SHOCK PROTEIN 90. In such a complex, CARP9 stabilizes AGO1 and mature miRNAs, allowing the proper loading of miRNAs in the effector complex.

Class I TCP transcription factors regulate trichome branching and cuticle development in Arabidopsis.

Trichomes and the cuticle are two specialized structures of the aerial epidermis that are important for plant organ development and interaction with the environment. In this study, we report that Arabidopsis thaliana plants affected in the function of the class I TEOSINTE BRANCHED 1, CYCLOIDEA, PCF (TCP) transcription factors TCP14 and TCP15 show overbranched trichomes in leaves and stems and increased cuticle permeability. We found that TCP15 regulates the expression of MYB106, a MIXTA-like transcription factor involved in epidermal cell and cuticle development, and overexpression of MYB106 in a tcp14 tcp15 mutant reduces trichome branch number. TCP14 and TCP15 are also required for the expression of the cuticle biosynthesis genes CYP86A4, GPAT6, and CUS2, and of SHN1 and SHN2, two AP2/EREBP transcription factors required for cutin and wax biosynthesis. SHN1 and CUS2 are also targets of TCP15, indicating that class I TCPs influence cuticle formation acting at different levels, through the regulation of MIXTA-like and SHN transcription factors and of cuticle biosynthesis genes. Our study indicates that class I TCPs are coordinators of the regulatory network involved in trichome and cuticle development.

CURLY LEAF Regulates MicroRNA Activity by Controlling ARGONAUTE 1 Degradation in Plants.

CURLY LEAF (CLF) encodes the methyltransferase subunit of the Polycomb Repressor Complex 2 (PRC2), which regulates the expression of target genes through H3K27 trimethylation. We isolated a new CLF mutant allele (clf-78) using a genetic screen designed to identify microRNA (miRNA) deficient mutants. CLF mutant plants showed impaired miRNA activity caused by increased ubiquitination and enhanced degradation of ARGONAUTE 1 (AGO1) in specific tissues. Such CLF-mediated AGO1 regulation was evident when plants were exposed to UV radiation, which caused increased susceptibility of clf mutants to some UV-induced responses. Furthermore, we showed that CLF directly regulates FBW2, which in turn triggers AGO1 degradation in the clf mutants. Interestingly, AGO1 bound to a target appeared particularly prone to degradation in the mutant plants, a process that was exacerbated when the complex bound a non-cleavable target. Thus, prolonged AGO1-target interaction seems to favor AGO1 degradation, suggesting that non-cleavable miRNA targets may overcome translation inhibition by modulating AGO1 stability in plants.

Light Triggers the miRNA-Biogenetic Inconsistency for De-etiolated Seedling Survivability in Arabidopsis thaliana.

The shift of dark-grown seedlings into light causes enormous transcriptome changes followed by a dramatic developmental transition. Here, we show that microRNA (miRNA) biogenesis also undergoes regulatory changes during de-etiolation. Etiolated seedlings maintain low levels of primary miRNAs (pri-miRNAs) and miRNA processing core proteins, such as Dicer-like 1, SERRATE, and HYPONASTIC LEAVES 1, whereas during de-etiolation both pri-miRNAs and the processing components accumulate to high levels. However, the levels of most miRNAs do not notably increase in response to light. To reconcile this inconsistency, we demonstrated that an unknown suppressor decreases miRNA-processing activity and light-induced SMALL RNA DEGRADING NUCLEASE 1 shortens the half-life of several miRNAs in de-etiolated seedlings. Taken together, these data suggest a novel mechanism, miRNA-biogenetic inconsistency, which accounts for the intricacy of miRNA biogenesis during de-etiolation. This mechanism is essential for the survival of de-etiolated seedlings after long-term skotomorphogenesis and their optimal adaptation to ever-changing light conditions.

Dynamic regulation of chromatin topology and transcription by inverted repeat-derived small RNAs in sunflower.

Transposable elements (TEs) are extremely abundant in complex plant genomes. siRNAs of 24 nucleotides in length control transposon activity in a process that involves de novo methylation of targeted loci. Usually, these epigenetic modifications trigger nucleosome condensation and a permanent silencing of the affected loci. Here, we show that a TE-derived inverted repeat (IR) element, inserted near the sunflower HaWRKY6 locus, dynamically regulates the expression of the gene by altering chromatin topology. The transcripts of this IR element are processed into 24-nt siRNAs, triggering DNA methylation on its locus. These epigenetic marks stabilize the formation of tissue-specific loops in the chromatin. In leaves, an intragenic loop is formed, blocking HaWRKY6 transcription. While in cotyledons (Cots), formation of an alternative loop, encompassing the whole HaWRKY6 gene, enhances transcription of the gene. The formation of this loop changes the promoter directionality, reducing IR transcription, and ultimately releasing the loop. Our results provide evidence that TEs can act as active and dynamic regulatory elements within coding loci in a mechanism that combines RNA silencing, epigenetic modification, and chromatin remodeling machineries.

Alternative use of miRNA-biogenesis co-factors in plants at low temperatures.

Plants use molecular mechanisms to sense temperatures, trigger quick adaptive responses and thereby cope with environmental changes. MicroRNAs (miRNAs) are key regulators of plant development under such conditions. The catalytic action of DICER LIKE 1 (DCL1), in conjunction with HYPONASTIC LEAVES 1 (HYL1) and SERRATE (SE), produces miRNAs from double-stranded RNAs. As plants lack a stable internal temperature to which enzymatic reactions could be optimized during evolution, reactions such as miRNA processing have to be adjusted to fluctuating environmental temperatures. Here, we report that with decreasing ambient temperature, the plant miRNA biogenesis machinery becomes more robust, producing miRNAs even in the absence of the key DCL1 co-factors HYL1 and SE. This reduces the morphological and reproductive defects of se and hyl1 mutants, restoring seed production. Using small RNA-sequencing and bioinformatics analyses, we have identified specific miRNAs that become HYL1/SE independent for their production in response to temperature decrease. We found that the secondary structure of primary miRNAs is key for this temperature recovery. This finding may have evolutionary implications as a potential adaptation-driving mechanism to a changing climate.

A Quick HYL1-Dependent Reactivation of MicroRNA Production Is Required for a Proper Developmental Response after Extended Periods of Light Deprivation.

Light is the most influential environmental stimulus for plant growth. In response to deficient light, plants reprogram their development to adjust their growth in search for a light source. A fine reprogramming of gene expression orchestrates this adaptive trait. Here we show that plants alter microRNA (miRNA) biogenesis in response to light transition. When plants suffer an unusual extended period of light deprivation, the miRNA biogenesis factor HYPONASTIC LEAVES 1 (HYL1) is degraded but an inactive pool of phosphorylated protein remains stable inside the nucleus. Degradation of HYL1 leads to the release of gene silencing, triggering a proper response to dark and shade. Upon light restoration, a quick dephosphorylation of HYL1 leads to the reactivation of miRNA biogenesis and a switch toward a developmental program that maximizes the light uptake. Our findings define a unique and fast regulatory mechanism controlling the plant silencing machinery during plant light response.

Robustness of Transposable Element Regulation but No Genomic Shock Observed in Interspecific Arabidopsis Hybrids.

The merging of two divergent genomes in a hybrid is believed to trigger a "genomic shock", disrupting gene regulation and transposable element (TE) silencing. Here, we tested this expectation by comparing the pattern of expression of transposable elements in their native and hybrid genomic context. For this, we sequenced the transcriptome of the Arabidopsis thaliana genotype Col-0, the A. lyrata genotype MN47 and their F1 hybrid. Contrary to expectations, we observe that the level of TE expression in the hybrid is strongly correlated to levels in the parental species. We detect that at most 1.1% of expressed transposable elements belonging to two specific subfamilies change their expression level upon hybridization. Most of these changes, however, are of small magnitude. We observe that the few hybrid-specific modifications in TE expression are more likely to occur when TE insertions are close to genes. In addition, changes in epigenetic histone marks H3K9me2 and H3K27me3 following hybridization do not coincide with TEs with changed expression. Finally, we further examined TE expression in parents and hybrids exposed to severe dehydration stress. Despite the major reorganization of gene and TE expression by stress, we observe that hybridization does not lead to increased disorganization of TE expression in the hybrid. Although our study did not examine TE transposition activity in hybrids, the examination of the transcriptome shows that TE expression is globally robust to hybridization. The term "genomic shock" is perhaps not appropriate to describe transcriptional modification in a viable hybrid merging divergent genomes.

A role for LAX2 in regulating xylem development and lateral-vein symmetry in the leaf.

Background and Aims: The symmetry of venation patterning in leaves is highly conserved within a plant species. Auxins are involved in this process and also in xylem vasculature development. Studying transgenic Arabidopsis plants ectopically expressing the sunflower transcription factor HaHB4, it was observed that there was a significant lateral-vein asymmetry in leaves and in xylem formation compared to wild type plants. To unravel the molecular mechanisms behind this phenotype, genes differentially expressed in these plants and related to auxin influx were investigated. Methods: Candidate genes responsible for the observed phenotypes were selected using a co-expression analysis. Single and multiple mutants in auxin influx carriers were characterized by morphological, physiological and molecular techniques. The analysis was further complemented by restoring the wild type (WT) phenotype by mutant complementation studies and using transgenic soybean plants ectopically expressing HaHB4 . Key Results: LAX2 , down-regulated in HaHB4 transgenic plants, was bioinformatically chosen as a candidate gene. The quadruple mutant aux1 lax1 lax2 lax3 and the single mutants, except lax1, presented an enhanced asymmetry in venation patterning. Additionally, the xylem vasculature of the lax2 mutant and the HaHB4 -expressing plants differed from the WT vasculature, including increased xylem length and number of xylem cell rows. Complementation of the lax2 mutant with the LAX2 gene restored both lateral-vein symmetry and xylem/stem area ratio in the stem, showing that auxin homeostasis is required to achieve normal vascular development. Interestingly, soybean plants ectopically expressing HaHB4 also showed an increased asymmetry in the venation patterning, accompanied by the repression of several GmLAX genes. Conclusions: Auxin influx carriers have a significant role in leaf venation pattering in leaves and, in particular, LAX2 is required for normal xylem development, probablt controlling auxin homeostasis.

A uORF Represses the Transcription Factor AtHB1 in Aerial Tissues to Avoid a Deleterious Phenotype.

AtHB1 is an Arabidopsis (Arabidopsis thaliana) homeodomain-leucine zipper transcription factor that participates in hypocotyl elongation under short-day conditions. Here, we show that its expression is posttranscriptionally regulated by an upstream open reading frame (uORF) located in its 5' untranslated region. This uORF encodes a highly conserved peptide (CPuORF) that is present in varied monocot and dicot species. The Arabidopsis uORF and its maize (Zea mays) homolog repressed the translation of the main open reading frame in cis, independent of the sequence of the latter. Published ribosome footprinting results and the analysis of a frame-shifted uORF, in which the repression capability was lost, indicated that the uORF causes ribosome stalling. The regulation exerted by the CPuORF was tissue specific and did not act in the absence of light. Moreover, a photosynthetic signal is needed for the CPuORF action, since plants with uncoupled chloroplasts did not show uORF-dependent repression. Plants transformed with the native AtHB1 promoter driving AtHB1 expression did not show differential phenotypes, whereas those transformed with a construct in which the uORF was mutated exhibited serrated leaves, compact rosettes, and, most significantly, short nondehiscent anthers and siliques containing fewer or no seeds. Thus, we propose that the uncontrolled expression of AtHB1 is deleterious for the plant and, hence, finely repressed by a translational mechanism.

The Footprint of Polygenic Adaptation on Stress-Responsive Cis-Regulatory Divergence in the Arabidopsis Genus.

Adaptation of a complex trait often requires the accumulation of many modifications to finely tune its underpinning molecular components to novel environmental requirements. The investigation of cis-acting regulatory modifications can be used to pinpoint molecular systems partaking in such complex adaptations. Here, we identify cis-acting modifications with the help of an interspecific crossing scheme designed to distinguish modifications derived in each of the two sister species, Arabidopsis halleri and A. lyrata Allele-specific expression levels were assessed in three environmental conditions chosen to reflect interspecific ecological differences: cold exposure, dehydration, and standard conditions. The functions described by Gene Ontology categories enriched in cis-acting mutations are markedly different in A. halleri and A. lyrata, suggesting that polygenic adaptation reshaped distinct polygenic molecular functions in the two species. In the A. halleri lineage, an excess of cis-acting changes affecting metal transport and homeostasis was observed, confirming that the well-known heavy metal tolerance of this species is the result of polygenic selection. In A. lyrata, we find a marked excess of cis-acting changes among genes showing a transcriptional response to cold stress in the outgroup species A. thaliana The adaptive relevance of these changes will have to be validated. We finally observed that polygenic molecular functions enriched in derived cis-acting changes are more constrained at the amino acid level. Using the distribution of cis-acting variation to tackle the polygenic basis of adaptation thus reveals the contribution of mutations of small effect to Darwinian adaptation.

The cytochrome c oxidase biogenesis factor AtCOX17 modulates stress responses in Arabidopsis.

COX17 is a soluble protein from the mitochondrial intermembrane space that participates in the transfer of copper for cytochrome c oxidase (COX) assembly in eukaryotic organisms. In this work, we studied the function of both Arabidopsis thaliana AtCOX17 genes using plants with altered expression levels of these genes. Silencing of AtCOX17-1 in a cox17-2 knockout background generates plants with smaller rosettes and decreased expression of genes involved in the response of plants to different stress conditions, including several genes that are induced by mitochondrial dysfunctions. Silencing of either of the AtCOX17 genes does not affect plant development or COX activity but causes a decrease in the response of genes to salt stress. In addition, these plants contain higher reactive oxygen and lipid peroxidation levels after irrigation with high NaCl concentrations and are less sensitive to abscisic acid. In agreement with a role of AtCOX17 in stress and abscisic acid responses, both AtCOX17 genes are induced by several stress conditions, abscisic acid and mutation of the transcription factor ABI4. The results indicate that AtCOX17 is required for optimal expression of a group of stress-responsive genes, probably as a component of signalling pathways that link stress conditions to gene expression responses.