MiRNAs differentially expressed in vegetative and reproductive organs of Marchantia polymorpha - insights into their expression pattern, gene structures and function.

MicroRNAs regulate gene expression affecting a variety of plant developmental processes. The evolutionary position of Marchantia polymorpha makes it a significant model to understand miRNA-mediated gene regulatory pathways in plants. Previous studies focused on conserved miRNA-target mRNA modules showed their critical role in Marchantia development. Here, we demonstrate that the differential expression of conserved miRNAs among land plants and their targets in selected organs of Marchantia additionally underlines their role in regulating fundamental developmental processes. The main aim of this study was to characterize selected liverwort-specific miRNAs, as there is a limited knowledge on their biogenesis, accumulation, targets, and function in Marchantia. We demonstrate their differential accumulation in vegetative and generative organs. We reveal that all liverwort-specific miRNAs examined are encoded by independent transcriptional units. MpmiR11737a, MpmiR11887 and MpmiR11796, annotated as being encoded within protein-encoding genes, have their own independent transcription start sites. The analysis of selected liverwort-specific miRNAs and their pri-miRNAs often reveal correlation in their levels, suggesting transcriptional regulation. However, MpmiR11796 shows a reverse correlation to its pri-miRNA level, suggesting post-transcriptional regulation. Moreover, we identify novel targets for selected liverwort-specific miRNAs and demonstrate an inverse correlation between their expression and miRNA accumulation. In the case of one miRNA precursor, we provide evidence that it encodes two functional miRNAs with two independent targets. Overall, our research sheds light on liverwort-specific miRNA gene structure, provides new data on their biogenesis and expression regulation. Furthermore, identifying their targets, we hypothesize the potential role of these miRNAs in early land plant development and functioning.

MicroRNA172b-5p/trehalose-6-phosphate synthase module stimulates trehalose synthesis and microRNA172b-3p/AP2-like module accelerates flowering in barley upon drought stress.

MicroRNAs (miRNAs) are major regulators of gene expression during plant development under normal and stress conditions. In this study, we analyzed the expression of 150 conserved miRNAs during drought stress applied to barley ready to flower. The dynamics of miRNAs expression was also observed after rewatering. Target messenger RNA (mRNAs) were experimentally identified for all but two analyzed miRNAs, and 41 of the targets were not reported before. Drought stress applied to barley induced accelerated flowering coordinated by a pair of two differently expressed miRNAs originating from a single precursor: hvu-miR172b-3p and hvu-miR172b-5p. Increased expression of miRNA172b-3p during drought leads to the downregulation of four APETALA2(AP2)-like genes by their mRNA cleavage. In parallel, the downregulation of the miRNA172b-5p level results in an increased level of a newly identified target, trehalose-6-phosphate synthase, a key enzyme in the trehalose biosynthesis pathway. Therefore, drought-treated plants have higher trehalose content, a known osmoprotectant, whose level is rapidly dropping after watering. In addition, trehalose-6-phosphate, an intermediate of the trehalose synthesis pathway, is known to induce flowering. The hvu-miRNA172b-5p/trehalose-6-phosphate synthase and hvu-miRNA172b-3p/AP2-like create a module leading to osmoprotection and accelerated flowering induction during drought.

Transcriptome Dynamics in Triticum aestivum Genotypes Associated with Resistance against the Wheat Dwarf Virus.

Wheat dwarf virus (WDV) is one of the most important pathogens of cereal crops worldwide. To understand the molecular mechanism of resistance, here we investigated the comparative transcriptome of wheat genotypes with different levels of resistance (Svitava and Fengyou 3) and susceptibility (Akteur) to WDV. We found a significantly higher number of differentially expressed transcripts (DETs) in the susceptible genotype than in the resistant one (e.g., Svitava). The number of downregulated transcripts was also higher in the susceptible genotype than in the resistant one (Svitava) and the opposite was true for the upregulated transcripts. Further functional analysis of gene ontology (GO) enrichment identified a total of 114 GO terms for the DETs. Of these, 64 biological processes, 28 cellular components and 22 molecular function GO terms were significantly enriched. A few of these genes appear to have a specific expression pattern related to resistance or susceptibility to WDV infection. Validation of the expression pattern by RT-qPCR showed that glycosyltransferase was significantly downregulated in the susceptible genotype compared to the resistant genotypes after WDV infection, while CYCLIN-T1-3, a regulator of CDK kinases (cyclin-dependent kinase), was upregulated. On the other hand, the expression pattern of the transcription factor (TF) MYB (TraesCS4B02G174600.2; myeloblastosis domain of transcription factor) was downregulated by WDV infection in the resistant genotypes compared to the susceptible genotype, while a large number of TFs belonging to 54 TF families were differentially expressed due to WDV infection. In addition, two transcripts (TraesCS7A02G341400.1 and TraesCS3B02G239900.1) were upregulated with uncharacterised proteins involved in transport and regulation of cell growth, respectively. Altogether, our findings showed a clear gene expression profile associated with resistance or susceptibility of wheat to WDV. In future studies, we will explore the regulatory network within the same experiment context. This knowledge will broaden not only the future for the development of virus-resistant wheat genotypes but also the future of genetic improvement of cereals for resilience and WDV-resistance breeding.

Genomic, transcriptomic, and metabolomic analysis of Oldenlandia corymbosa reveals the biosynthesis and mode of action of anti-cancer metabolites.

Plants accumulate a vast array of secondary metabolites, which constitute a natural resource for pharmaceuticals. Oldenlandia corymbosa belongs to the Rubiaceae family, and has been used in traditional medicine to treat different diseases, including cancer. However, the active metabolites of the plant, their biosynthetic pathway and mode of action in cancer are unknown. To fill these gaps, we exposed this plant to eight different stress conditions and combined different omics data capturing gene expression, metabolic profiles, and anti-cancer activity. Our results show that O. corymbosa extracts are active against breast cancer cell lines and that ursolic acid is responsible for this activity. Moreover, we assembled a high-quality genome and uncovered two genes involved in the biosynthesis of ursolic acid. Finally, we also revealed that ursolic acid causes mitotic catastrophe in cancer cells and identified three high-confidence protein binding targets by Cellular Thermal Shift Assay (CETSA) and reverse docking. Altogether, these results constitute a valuable resource to further characterize the biosynthesis of active metabolites in the Oldenlandia group, while the mode of action of ursolic acid will allow us to further develop this valuable compound.

Excess nitrogen responsive HvMADS27 transcription factor controls barley root architecture by regulating abscisic acid level.

Nitrogen (N) is an important element for plant growth and development. Although several studies have examined plants' response to N deficiency, studies on plants' response to excess N, which is common in fertilizer-based agrosystems, are limited. Therefore, the aim of this study was to examine the response of barley to excess N conditions, specifically the root response. Additionally, genomic mechanism of excess N response in barley was elucidated using transcriptomic technologies. The results of the study showed that barley MADS27 transcription factor was mainly expressed in the roots and its gene contained N-responsive cis-regulatory elements in the promoter region. Additionally, there was a significant decrease in HvMADS27 expression under excess N condition; however, its expression was not significantly affected under low N condition. Phenotypic analysis of the root system of HvMADS27 knockdown and overexpressing barley plants revealed that HvMADS27 regulates barley root architecture under excess N stress. Further analysis of wild-type (WT) and transgenic barley plants (hvmads27 kd and hvmads27 c-Myc OE) revealed that HvMADS27 regulates the expression of HvBG1 ß-glucosidase, which in turn regulates abscisic acid (ABA) level in roots. Overall, the findings of this study showed that HvMADS27 expression is downregulated in barley roots under excess N stress, which induces HvBG1 expression, leading to the release of ABA from ABA-glucose conjugate, and consequent shortening of the roots.

MicroRNA biogenesis and activity in plant cell dedifferentiation stimulated by cell wall removal.

BACKGROUND: Despite the frequent use of protoplast-to-plant system in in vitro cultures of plants, the molecular mechanisms regulating the first and most limiting stages of this process, i.e., protoplast dedifferentiation and the first divisions leading to the formation of a microcallus, have not been elucidated. RESULTS: In this study, we investigated the function of miRNAs in the dedifferentiation of A. thaliana mesophyll cells in a process stimulated by the enzymatic removal of the cell wall. Leaf cells, protoplasts and CDPs (cells derived from protoplasts) cultured for 24, 72 and 120 h (first cell division). In protoplasts, a strong decrease in the amount of AGO1 in both the nucleus and the cytoplasm, as well as dicing bodies (DBs), which are considered to be sites of miRNA biogenesis, was shown. However during CDPs division, the amounts of AGO1 and DBs strongly increased. MicroRNA transcriptome studies demonstrated that lower amount of differentially expressed miRNAs are present in protoplasts than in CDPs cultured for 120 h. Then analysis of differentially expressed miRNAs, selected pri-miRNA and mRNA targets were performed. CONCLUSION: This result indicates that miRNA function is not a major regulation of gene expression in the initial but in later steps of dedifferentiation during CDPs divisions. miRNAs participate in organogenesis, oxidative stress, nutrient deficiencies and cell cycle regulation in protoplasts and CDPs. The important role played by miRNAs in the process of dedifferentiation of mesophyll cells was confirmed by the increased mortality and reduced cell division of CDPs derived from mutants with defective miRNA biogenesis and miR319b expression.

Pi-starvation induced transcriptional changes in barley revealed by a comprehensive RNA-Seq and degradome analyses.

BACKGROUND: Small RNAs (sRNAs) are 20-30 nt regulatory elements which are responsible for plant development regulation and participate in many plant stress responses. Insufficient inorganic phosphate (Pi) concentration triggers plant responses to balance the internal Pi level. RESULTS: In this study, we describe Pi-starvation-responsive small RNAs and transcriptome changes in barley (Hordeum vulgare L.) using Next-Generation Sequencing (NGS) RNA-Seq data derived from three different types of NGS libraries: (i) small RNAs, (ii) degraded RNAs, and (iii) functional mRNAs. We find that differentially and significantly expressed miRNAs (DEMs, Bonferroni adjusted p-value < 0.05) are represented by 15 molecules in shoot and 13 in root; mainly various miR399 and miR827 isomiRs. The remaining small RNAs (i.e., those without perfect match to reference sequences deposited in miRBase) are considered as differentially expressed other sRNAs (DESs, p-value Bonferroni correction < 0.05). In roots, a more abundant and diverse set of other sRNAs (DESs, 1796 unique sequences, 0.13% from the average of the unique small RNA expressed under low-Pi) contributes more to the compensation of low-Pi stress than that in shoots (DESs, 199 unique sequences, 0.01%). More than 80% of differentially expressed other sRNAs are up-regulated in both organs. Additionally, in barley shoots, up-regulation of small RNAs is accompanied by strong induction of two nucleases (S1/P1 endonuclease and 3'-5' exonuclease). This suggests that most small RNAs may be generated upon nucleolytic cleavage to increase the internal Pi pool. Transcriptomic profiling of Pi-starved barley shoots identifies 98 differentially expressed genes (DEGs). A majority of the DEGs possess characteristic Pi-responsive cis-regulatory elements (P1BS and/or PHO element), located mostly in the proximal promoter regions. GO analysis shows that the discovered DEGs primarily alter plant defense, plant stress response, nutrient mobilization, or pathways involved in the gathering and recycling of phosphorus from organic pools. CONCLUSIONS: Our results provide comprehensive data to demonstrate complex responses at the RNA level in barley to maintain Pi homeostasis and indicate that barley adapts to Pi-starvation through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.

Core spliceosomal Sm proteins as constituents of cytoplasmic mRNPs in plants.

In recent years, research has increasingly focused on the key role of post-transcriptional regulation of messenger ribonucleoprotein (mRNP) function and turnover. As a result of the complexity and dynamic nature of mRNPs, the full composition of a single mRNP complex remains unrevealed and mRNPs are poorly described in plants. Here we identify canonical Sm proteins as part of the cytoplasmic mRNP complex, indicating their function in the post-transcriptional regulation of gene expression in plants. Sm proteins comprise an evolutionarily ancient family of small RNA-binding proteins involved in pre-mRNA splicing. The latest research indicates that Sm could also impact on mRNA at subsequent stages of its life cycle. In this work we show that in the microsporocyte cytoplasm of Larix decidua, the European larch, Sm proteins accumulate within distinct cytoplasmic bodies, also containing polyadenylated RNA. To date, several types of cytoplasmic bodies involved in the post-transcriptional regulation of gene expression have been described, mainly in animal cells. Their role and molecular composition in plants remain less well established, however. A total of 222 mRNA transcripts have been identified as cytoplasmic partners for Sm proteins. The specific colocalization of these mRNAs with Sm proteins within cytoplasmic bodies has been confirmed via microscopic analysis. The results from this work support the hypothesis, that evolutionarily conserved Sm proteins have been adapted to perform a whole repertoire of functions related to the post-transcriptional regulation of gene expression in Eukaryota. This adaptation presumably enabled them to coordinate the interdependent processes of splicing element assembly, mRNA maturation and processing, and mRNA translation regulation, and its degradation.

A Functional Network of Novel Barley MicroRNAs and Their Targets in Response to Drought.

The regulation of mRNA (messenger RNA) levels by microRNA-mediated activity is especially important in plant responses to environmental stresses. In this work, we report six novel barley microRNAs, including two processed from the same precursor that are severely downregulated under drought conditions. For all analyzed microRNAs, we found target genes that were upregulated under drought conditions and that were known to be involved in a plethora of processes from disease resistance to chromatin-protein complex formation and the regulation of transcription in mitochondria. Targets for novel barley microRNAs were confirmed through degradome data analysis and RT-qPCR using primers flanking microRNA-recognition site. Our results show a broad transcriptional response of barley to water deficiency conditions through microRNA-mediated gene regulation and facilitate further research on drought tolerance in crops.

A stable tRNA-like molecule is generated from the long noncoding RNA GUT15 in Arabidopsis.

The Arabidopsis GUT15 RNA belongs to a class of noncoding RNAs that are expressed from the intergenic regions of protein-coding genes. We show that the RNA polymerase II transcribed GUT15 transcript serves as a precursor for two stable RNA species, a tRNA-like molecule and GUT15-tRF-F5, which are both encoded by the final intron in the GUT15 gene. The GUT15-encoded tRNA-like molecule cannot be autonomously transcribed by RNA polymerase III. However, this molecule contains a CCA motif, suggesting that it may enter the tRNA maturation pathway. The GUT15-encoded tRNA-like sequence has an inhibiting effect on the splicing of its host intron. Moreover, we demonstrate that the canonical tRNA genes nested within introns do not affect the splicing patterns of their host protein-coding transcripts.

tRex: A Web Portal for Exploration of tRNA-Derived Fragments in Arabidopsis thaliana.

tRNA-derived fragments (tRFs) constitute a new class of short regulatory RNAs that are a product of nascent or mature tRNA processing. tRF sequences have been identified in all domains of life; however, most published research pertains to human, yeast and some bacterial organisms. Despite growing interest in plant tRFs and accumulating evidence of their function in plant development and stress responses, no public, web-based repository dedicated to these molecules is currently available. Here, we introduce tRex (http://combio.pl/trex)-the first comprehensive data-driven online resource specifically dedicated to tRFs in the model plant Arabidopsis thaliana. The portal is based on verified Arabidopsis tRNA annotation and includes in-house-generated and publicly available small RNA sequencing experiments from various tissues, ecotypes, genotypes and stress conditions. The provided web-based tools are designed in a user-friendly manner and allow for seamless exploration of the data that are presented in the form of dynamic tables and cumulative coverage profiles. The tRex database is connected to external genomic and citation resources, which makes it a one-stop solution for Arabidopsis tRF-related research.

Maf1-mediated regulation of yeast RNA polymerase III is correlated with CCA addition at the 3' end of tRNA precursors.

In eukaryotic cells tRNA synthesis is negatively regulated by the protein Maf1, conserved from yeast to humans. Maf1 from yeast Saccharomyces cerevisiae mediates repression of trna transcription when cells are transferred from medium with glucose to medium with glycerol, a non-fermentable carbon source. The strain with deleted gene encoding Maf1 (maf1Δ) is viable but accumulates tRNA precursors. In this study tRNA precursors were analysed by RNA-Seq and Northern hybridization in wild type strain and maf1Δ mutant grown in glucose medium or upon shift to repressive conditions. A negative effect of maf1Δ mutant on the addition of the auxiliary CCA nucleotides to the 3' end of pre-tRNAs was observed in cells shifted to unfavourable growth conditions. This effect was reduced by overexpression of the yeast CCA1 gene encoding ATP(CTP):tRNA nucleotidyltransferase. The CCA sequence at the 3' end is important for export of tRNA precursors from the nucleus and essential for tRNA charging with amino acids. Data presented here indicate that CCA-addition to intron-containing end-processed tRNA precursors is a limiting step in tRNA maturation when there is no Maf1 mediated RNA polymerase III (Pol III) repression. The correlation between CCA synthesis and Pol III regulation by Maf1 could be important in coordination of tRNA transcription, processing and regulation of translation.

Developmental changes in barley microRNA expression profiles coupled with miRNA target analysis.

MicroRNAs are 19- to 24-nt-long single-stranded RNAs that are crucial regulators of gene expression which control plant development and response to environmental cues. We have analyzed microtranscriptomes of five barley developmental stages. Generally, during the barley development, miR168-3p and miR1432-5p levels increase while the 5'U-miR156-5p level decreases (with exception for the 2-week-old barley). We have identified two miR156-5p izomiRs (called 5'U-miR156-5p [20 nt] and 5'UU-miR156-5p [21 nt]), which were expressed differently during barley development. The 5' U-miR156-5p level decreased in 3-week-, 6-week-, and 68-day-old barley, when compared to the 1-week-old plants. Meanwhile, the 5' UU-miR156-5p level increased significantly in the 68-day-old barley plants. Moreover, only the 5' U-miR156 isomiR recognizes and guides unique transcription factor mRNAs from the Squamosa Promoter Binding Protein-Like (SPL) family. We identified many non-canonical microRNAs with changed expression levels during the barley development. Here, we present the profiles of microRNA expression characteristics for particular barley developmental stages. These analyses are accompanied by the experimental degradome analysis of miRNA targets.

Arabidopsis microRNA expression regulation in a wide range of abiotic stress responses.

Arabidopsis microRNA expression regulation was studied in a wide array of abiotic stresses such as drought, heat, salinity, copper excess/deficiency, cadmium excess, and sulfur deficiency. A home-built RT-qPCR mirEX platform for the amplification of 289 Arabidopsis microRNA transcripts was used to study their response to abiotic stresses. Small RNA sequencing, Northern hybridization, and TaqMan® microRNA assays were performed to study the abundance of mature microRNAs. A broad response on the level of primary miRNAs (pri-miRNAs) was observed. However, stress response at the level of mature microRNAs was rather confined. The data presented show that in most instances, the level of a particular mature miRNA could not be predicted based on the level of its pri-miRNA. This points to an essential role of posttranscriptional regulation of microRNA expression. New Arabidopsis microRNAs responsive to abiotic stresses were discovered. Four microRNAs: miR319a/b, miR319b.2, and miR400 have been found to be responsive to several abiotic stresses and thus can be regarded as general stress-responsive microRNA species.

The liverwort Pellia endiviifolia shares microtranscriptomic traits that are common to green algae and land plants.

Liverworts are the most basal group of extant land plants. Nonetheless, the molecular biology of liverworts is poorly understood. Gene expression has been studied in only one species, Marchantia polymorpha. In particular, no microRNA (miRNA) sequences from liverworts have been reported. Here, Illumina-based next-generation sequencing was employed to identify small RNAs, and analyze the transcriptome and the degradome of Pellia endiviifolia. Three hundred and eleven conserved miRNA plant families were identified, and 42 new liverwort-specific miRNAs were discovered. The RNA degradome analysis revealed that target mRNAs of only three miRNAs (miR160, miR166, and miR408) have been conserved between liverworts and other land plants. New targets were identified for the remaining conserved miRNAs. Moreover, the analysis of the degradome permitted the identification of targets for 13 novel liverwort-specific miRNAs. Interestingly, three of the liverwort microRNAs show high similarity to previously reported miRNAs from Chlamydomonas reinhardtii. This is the first observation of miRNAs that exist both in a representative alga and in the liverwort P. endiviifolia but are not present in land plants. The results of the analysis of the P. endivifolia microtranscriptome support the conclusions of previous studies that placed liverworts at the root of the land plant evolutionary tree of life.

Transcriptionally and post-transcriptionally regulated microRNAs in heat stress response in barley.

Heat stress is one of the major abiotic factors that can induce severe plant damage, leading to a decrease in crop plant productivity. Despite barley being a cereal of great economic importance, few data are available concerning its thermotolerance mechanisms. In this work microRNAs (miRNAs) involved in heat stress response in barley were investigated. The level of selected barley mature miRNAs was examined by hybridization. Quantitative real-time PCR (RT-qPCR) was used to monitor the changes in the expression profiles of primary miRNA (pri-miRNA) precursors, as well as novel and conserved target genes during heat stress. The miRNA-mediated cleavage sites in the target transcripts were confirmed by degradome analysis and the 5' RACE (rapid amplification of cDNA ends) approach. Four barley miRNAs (miR160a, 166a, 167h, and 5175a) were found which are heat stress up-regulated at the level of both mature miRNAs and precursor pri-miRNAs. Moreover, the splicing of introns hosting miR160a and miR5175a is also heat induced. The results demonstrate transcriptional and post-transcriptional regulation of heat-responsive miRNAs in barley. The observed induction of miRNA expression is correlated with the down-regulation of the expression level of their experimentally identified new and conservative target genes.

Identification of nutrient-responsive Arabidopsis and rapeseed microRNAs by comprehensive real-time polymerase chain reaction profiling and small RNA sequencing.

Comprehensive expression profiles of Arabidopsis (Arabidopsis thaliana) MIRNA genes and mature microRNAs (miRs) are currently not available. We established a quantitative real-time polymerase chain reaction platform that allows rapid and sensitive quantification of 177 Arabidopsis primary miR transcripts (pri-miRs). The platform was used to detect phosphorus (P) or nitrogen (N) status-responsive pri-miR species. Several pri-miR169 species as well as pri-miR398a were found to be repressed during N limitation, whereas during P limitation, pri-miR778, pri-miR827, and pri-miR399 species were induced and pri-miR398a was repressed. The corresponding responses of the biologically active, mature miRs were confirmed using specific stem-loop reverse transcription primer quantitative polymerase chain reaction assays and small RNA sequencing. Interestingly, the latter approach also revealed high abundance of some miR star strands. Bioinformatic analysis of small RNA sequences with a modified miRDeep algorithm led to the identification of the novel P limitation-induced miR2111, which is encoded by two loci in the Arabidopsis genome. Furthermore, miR2111, miR169, a miR827-like sequence, and the abundances of several miR star strands were found to be strongly dependent on P or N status in rapeseed (Brassica napus) phloem sap, flagging them as candidate systemic signals. Taken together, these results reveal the existence of complex small RNA-based regulatory networks mediating plant adaptation to mineral nutrient availability.

Yellow lupine cyclophilin interacts with nucleic acids.

To investigate properties of yellow lupine cytosolic cyclophilin, an expression vector pET15CYP was constructed. The CyP cDNA (GenBank accession no.Y16088) reveals an open reading frame of 172 amino acids with the conserved tryptophan residue at position 128 and an insertion of seven amino acids spanning positions 48-54. Yellow lupine cyclophilin, purified after expression in E. coli cells, exhibits peptidyl-prolyl cis/trans isomerase activity when assayed with a synthetic oligopeptide. We have demonstrated that the recombinant cyclophilin is able to interact with nucleic acids, both single and double stranded DNA fragments as well as RNA.

[Recombinant protein production in Escherichia coli]. / Produkcja rekombinowanych bialek w Escherichia coli.

Growing needs for efficient recombinant production pose new challenges; starting from cell growth optimization under overexpression conditions, improving vectors, gene and protein sequence to suit them to protein biosynthesis machinery of the host, through extending the knowledge of protein folding, fusion protein construction, and coexpression systems, to improvements in protein purification and renaturation technologies. Hitherto Escherichia coli is the most defined and the cheapest protein biosynthesis system. With its wealth of available mutants tested is the best suited to economically test new gene constructs and to scale up the recombinant protein production.