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.

Identification of transcription factors that bind to the 5'-UTR of the barley PHO2 gene.

In barley and other higher plants, phosphate homeostasis is maintained by a regulatory network involving the PHO2 (PHOSPHATE2) encoding ubiquitin-conjugating (UBC) E2 enzyme, the PHR1 (PHOSPHATE STARVATION RESPONSE 1) transcription factor (TF), IPS1 (INDUCED BYPHOSPHATESTARVATION1) RNA, and miR399. During phosphate ion (Pi) deprivation, PHR1 positively regulates MIR399 expression, after transcription and processing mature miR399 guides the Ago protein to the 5'-UTR of PHO2 transcripts. Non-coding IPS1 RNA is highly expressed during Pi starvation, and the sequestration of miR399 molecules protects PHO2 mRNA from complete degradation. Here, we reveal new cis- and trans-regulatory elements that are crucial for efficient PHO2 gene expression in barley. We found that the 5'-UTR of PHO2 contains two PHR1 binding sites (P1BSs) and one Pi-responsive PHO element. Using a yeast one-hybrid (Y1H) assay, we identified two candidate proteins that might mediate this transcriptional regulation: a barley PHR1 ortholog and a TF containing an uncharacterized MYB domain. Additional results classified this new potential TF as belonging to the APL (ALTERED PHLOEM DEVELOPMENT) protein family, and we observed its nuclear localization in barley protoplasts. Pi starvation induced the accumulation of barley APL transcripts in both the shoots and roots. Interestingly, the deletion of the P1BS motif from the first intron of the barley 5'-UTR led to a significant increase in the transcription of a downstream ß-glucuronidase (GUS) reporter gene in tobacco leaves. Our work extends the current knowledge about putative cis- and trans-regulatory elements that may affect the expression of the barley PHO2 gene.

Plant PHR Transcription Factors: Put on A Map.

The phosphate starvation response (PHR) protein family exhibits the MYB and coiled-coil domains. In plants, within the either 5' untranslated regions (UTRs) or promoter regions of phosphate starvation-induced (PSI) genes are characteristic cis-regulatory elements, namely PHR1 binding sequence (P1BS). The most widely studied PHR protein family members, such as AtPHR1 in Arabidopsis thaliana (L.) and OsPHR2 in Oryza sativa (L.), may activate the gene expression of a broad range of PSI genes by binding to such elements in a phosphate (Pi) dependent manner. In Pi signaling, PHR transcription factors (TFs) can be selectively activated or deactivated by other proteins to execute the final step of signal transduction. Several new proteins have been associated with the AtPHR1/OsPHR2 signaling cascade in the last few years. While the PHR TF transcriptional role has been studied intensively, here we highlight the recent findings of upstream molecular components and other signaling pathways that may interfere with the PHR final mode of action in plants. Detailed information about transcriptional regulation of the AtPHR1 gene itself and its upstream molecular events has been reviewed.

miRNA Detection by Stem-Loop RT-qPCR in Studying microRNA Biogenesis and microRNA Responsiveness to Abiotic Stresses.

This chapter is devoted to a PCR-based method for analyzing the expression level of mature miRNAs which utilizes the TaqMan® technology. Stem-loop RT-qPCR requires preparation of separate cDNA templates for each analyzed miRNA as reverse transcription occurs in the presence of a miRNA-specific stem-loop reverse primer. In quantitative analysis, SYBR® Green is not used but the more sensitive TaqMan® probe that on 5' end contains a covalently attached fluorophore and on 3' quencher. When quencher and fluorophore are spatially separated due to nucleolytic DNA polymerase activity, the signal is released and quantified. This section provides a detailed and comprehensive protocol allowing for the successful analysis of mature miRNA levels in analyzed sample. Reverse transcription combined with classic real-time PCR as well as ddPCR™ (Droplet Digital™ PCR) will be presented.

HorTILLUS-A Rich and Renewable Source of Induced Mutations for Forward/Reverse Genetics and Pre-breeding Programs in Barley (Hordeum vulgare L.).

TILLING (Targeting Induced Local Lesions IN Genomes) is a strategy used for functional analysis of genes that combines the classical mutagenesis and a rapid, high-throughput identification of mutations within a gene of interest. TILLING has been initially developed as a discovery platform for functional genomics, but soon it has become a valuable tool in development of desired alleles for crop breeding, alternative to transgenic approach. Here we present the HorTILLUS ( Hordeum-TILLING-University of Silesia) population created for spring barley cultivar "Sebastian" after double-treatment of seeds with two chemical mutagens: sodium azide (NaN3) and N-methyl-N-nitrosourea (MNU). The population comprises more than 9,600 M2 plants from which DNA was isolated, seeds harvested, vacuum-packed, and deposited in seed bank. M3 progeny of 3,481 M2 individuals was grown in the field and phenotyped. The screening for mutations was performed for 32 genes related to different aspects of plant growth and development. For each gene fragment, 3,072-6,912 M2 plants were used for mutation identification using LI-COR sequencer. In total, 382 mutations were found in 182.2 Mb screened. The average mutation density in the HorTILLUS, estimated as 1 mutation per 477 kb, is among the highest mutation densities reported for barley. The majority of mutations were G/C to A/T transitions, however about 8% transversions were also detected. Sixty-one percent of mutations found in coding regions were missense, 37.5% silent and 1.1% nonsense. In each gene, the missense mutations with a potential effect on protein function were identified. The HorTILLUS platform is the largest of the TILLING populations reported for barley and best characterized. The population proved to be a useful tool, both in functional genomic studies and in forward selection of barley mutants with required phenotypic changes. We are constantly renewing the HorTILLUS population, which makes it a permanent source of new mutations. We offer the usage of this valuable resource to the interested barley researchers on cooperative basis.

Heat Stress Affects Pi-related Genes Expression and Inorganic Phosphate Deposition/Accumulation in Barley.

Phosphorus (P) in plants is taken from soil as an inorganic phosphate (Pi) and is one of the most important macroelements in growth and development. Plants actively react to Pi starvation by the induced expression of Pi transporters, MIR399, MIR827, and miR399 molecular sponge - IPS1 genes and by the decreased expression of the ubiquitin-conjugating enzyme E2 (PHOSPHATE2 - PHO2) and Pi sensing and transport SPX-MFS genes. The PHO2 protein is involved in the degradation of Pi transporters PHT1;1 (from soil to roots) and PHO1 (from roots to shoots). The decreased expression of PHO2 leads to Pi accumulation in shoots. In contrast, the pho1 mutant shows a decreased level of Pi concentration in shoots. Finally, Pi starvation leads to decreased Pi concentration in all plant tissues. Little is known about plant Pi homeostasis in other abiotic stress conditions. We found that, during the first hour of heat stress, Pi accumulated in barley shoots but not in the roots, and transcriptomic data analysis as well as RT-qPCR led us to propose an explanation for this phenomenon. Pi transport inhibition from soil to roots is balanced by lower Pi efflux from roots to shoots directed by the PHO1 transporter. In shoots, the PHO2 mRNA level is decreased, leading to an increased Pi level. We concluded that Pi homeostasis in barley during heat stress is maintained by dynamic changes in Pi-related genes expression.