Identification and expression profiling of microRNAs in leaf tissues of Foeniculum vulgare Mill. under salinity stress.

Foeniculum vulgare Mill. commonly known as fennel, is a globally recognized aromatic medicinal plant and culinary herb with widespread popularity due to its antimicrobial, antioxidant, carminative, and diuretic properties, among others. Although the phenotypic effects of salinity stress have been previously explored in fennel, the molecular mechanisms underlying responses to elevated salinity in this plant remain elusive. MicroRNAs (miRNAs) are tiny, endogenous, and extensively conserved non-coding RNAs (ncRNAs) typically ranging from 20 to 24 nucleotides (nt) in length that play a major role in a myriad of biological functions. In fact, a number of miRNAs have been extensively associated with responses to abiotic stress in plants. Consequently, employing computational methodologies and rigorous filtering criteria, 40 putative miRNAs belonging to 25 different families were characterized from fennel in this study. Subsequently, employing the psRNATarget tool, a total of 67 different candidate target transcripts for the characterized fennel miRNAs were predicted. Additionally, the expression patterns of six selected fennel miRNAs (i.e. fvu-miR156a, fvu-miR162a-3p, fvu-miR166a-3p, fvu-miR167a-5p, fvu-miR171a-3p, and fvu-miR408-3p) were analyzed under salinity stress conditions via qPCR. This article holds notable significance as it identifies not only 40 putative miRNAs in fennel, a non-model plant, but also pioneers the analysis of their expression under salinity stress conditions.

Principles of miRNA/miRNA* function in plant MIRNA processing.

MicroRNAs (miRNAs) are essential regulators of gene expression, defined by their unique biogenesis, which requires the precise excision of the small RNA from an imperfect fold-back precursor. Unlike their animal counterparts, plant miRNA precursors exhibit variations in sizes and shapes. Plant MIRNAs can undergo processing in a base-to-loop or loop-to-base direction, with DICER-LIKE1 (DCL1) releasing the miRNA after two cuts (two-step MIRNAs) or more (sequential MIRNAs). In this study, we demonstrate the critical role of the miRNA/miRNA* duplex region in the processing of miRNA precursors. We observed that endogenous MIRNAs frequently experience suboptimal processing in vivo due to mismatches in the miRNA/miRNA* duplex, a key region that fine-tunes miRNA levels. Enhancing the interaction energy of the miRNA/miRNA* duplex in two-step MIRNAs results in a substantial increase in miRNA levels. Conversely, sequential MIRNAs display distinct and specific requirements for the miRNA/miRNA* duplexes along their foldback structure. Our work establishes a connection between the miRNA/miRNA* structure and precursor processing mechanisms. Furthermore, we reveal a link between the biological function of miRNAs and the processing mechanism of their precursors with the evolution of plant miRNA/miRNA* duplex structures.

A novel strategy to uncover specific GO terms/phosphorylation pathways in phosphoproteomic data in Arabidopsis thaliana.

BACKGROUND: Proteins are the workforce of the cell and their phosphorylation status tailors specific responses efficiently. One of the main challenges of phosphoproteomic approaches is to deconvolute biological processes that specifically respond to an experimental query from a list of phosphoproteins. Comparison of the frequency distribution of GO (Gene Ontology) terms in a given phosphoproteome set with that observed in the genome reference set (GenRS) is the most widely used tool to infer biological significance. Yet, this comparison assumes that GO term distribution between the phosphoproteome and the genome are identical. However, this hypothesis has not been tested due to the lack of a comprehensive phosphoproteome database. RESULTS: In this study, we test this hypothesis by constructing three phosphoproteome databases in Arabidopsis thaliana: one based in experimental data (ExpRS), another based in in silico phosphorylation protein prediction (PredRS) and a third that is the union of both (UnRS). Our results show that the three phosphoproteome reference sets show default enrichment of several GO terms compared to GenRS, indicating that GO term distribution in the phosphoproteomes does not match that of the genome. Moreover, these differences overshadow the identification of GO terms that are specifically enriched in a particular condition. To overcome this limitation, we present an additional comparison of the sample of interest with UnRS to uncover GO terms specifically enriched in a particular phosphoproteome experiment. Using this strategy, we found that mRNA splicing and cytoplasmic microtubule compounds are important processes specifically enriched in the phosphoproteome of dark-grown Arabidopsis seedlings. CONCLUSIONS: This study provides a novel strategy to uncover GO specific terms in phosphoproteome data of Arabidopsis that could be applied to any other organism. We also highlight the importance of specific phosphorylation pathways that take place during dark-grown Arabidopsis development.

Efficient method for isolation of high-quality RNA from Psidium guajava L. tissues.

Acquiring high-quality RNA in sufficient amounts is crucial in plant molecular biology and genetic studies. Several methods for RNA extraction from plants are available in the literature, mainly due to the great biochemical diversity present in each species and tissue, which can complicate or prevent the extraction. Psidium guajava (Myrtaceae family) is a perennial fruit tree of medicinal and economic value; nevertheless, only a few molecular studies are available for the species. One reason is the difficulty in obtaining RNA due to the content of the samples, which are rich in polyphenols, polysaccharides, and secondary metabolites. Furthermore, there are few studies available for the isolation of RNA from guava or Psidium samples, which hampers advances in the study of the genus. Here, quality and yields of RNA isolates were compared using six extraction protocols: two protocols based on the application of cetyltrimethylammonium bromide (CTAB) lysis buffer, one protocol which uses the TRIzol reagent, one which applies guanidine thiocyanate lysis buffer followed by organic phase extraction, and two commercial kits (PureLink RNA Mini Kit and RNeasy Plant Mini Kit). The CTAB-based method provided the highest RNA yields and quality for five different tissues (flower bud, immature leaf, young leaf, mature leaf, and root), genotypes, and stress conditions. For the most efficient protocol, the average yield of RNA from guava leaves was 203.06 µg/g of tissue, and the A260/A280 and A260/A230 ratios were 2.1 and 2.2, respectively. RT-qPCR analysis demonstrated that the purity of the samples was sufficient for molecular biology experiments. CTAB-based methods for RNA isolation were found to be the most efficient, providing the highest RNA yields and quality for tissues from P. guajava. Additionally, they were compatible for downstream RNA-based applications, besides being simple and cost-effective.

MicroRNA Zma-miR528 Versatile Regulation on Target mRNAs during Maize Somatic Embryogenesis.

MicroRNAs (miRNAs) are small non-coding RNAs that regulate the accumulation and translation of their target mRNAs through sequence complementarity. miRNAs have emerged as crucial regulators during maize somatic embryogenesis (SE) and plant regeneration. A monocot-specific miRNA, mainly accumulated during maize SE, is zma-miR528. While several targets have been described for this miRNA, the regulation has not been experimentally confirmed for the SE process. Here, we explored the accumulation of zma-miR528 and several predicted targets during embryogenic callus induction, proliferation, and plantlet regeneration using the maize cultivar VS-535. We confirmed the cleavage site for all tested zma-miR528 targets; however, PLC1 showed very low levels of processing. The abundance of zma-miR528 slightly decreased in one month-induced callus compared to the immature embryo (IE) explant tissue. However, it displayed a significant increase in four-month sub-cultured callus, coincident with proliferation establishment. In callus-regenerated plantlets, zma-miR528 greatly decreased to levels below those observed in the initial explant. Three of the target transcripts (MATE, bHLH, and SOD1a) showed an inverse correlation with the miRNA abundance in total RNA samples at all stages. Using polysome fractionation, zma-miR528 was detected in the polysome fraction and exhibited an inverse distribution with the PLC1 target, which was not observed at total RNA. Accordingly, we conclude that zma-miR528 regulates multiple target mRNAs during the SE process by promoting their degradation, translation inhibition or both.

Phaseolus vulgaris MIR1511 genotypic variations differentially regulate plant tolerance to aluminum toxicity.

The common-bean (Phaseolus vulgaris), a widely consumed legume, originated in Mesoamerica and expanded to South America, resulting in the development of two geographically distinct gene pools. Poor soil condition, including metal toxicity, are often constraints to common-bean crop production. Several P. vulgaris miRNAs, including miR1511, respond to metal toxicity. The MIR1511 gene sequence from the two P. vulgaris model sequenced genotypes revealed that, as opposed to BAT93 (Mesoamerican), the G19833 (Andean) accession displays a 58-bp deletion, comprising the mature and star miR1511 sequences. Genotyping-By-Sequencing data analysis from 87 non-admixed Phaseolus genotypes, comprising different Phaseolus species and P. vulgaris populations, revealed that all the P. vulgaris Andean genotypes and part of the Mesoamerican (MW1) genotypes analyzed displayed a truncated MIR1511 gene. The geographic origin of genotypes with a complete versus truncated MIR1511 showed a distinct distribution. The P. vulgaris ALS3 (Aluminum Sensitive Protein 3) gene, known to be important for aluminum detoxification in several plants, was experimentally validated as the miR1511 target. Roots from BAT93 plants showed decreased miR1511 and increased ALS3 transcript levels at early stages under aluminum toxicity (AlT), while G19833 plants, lacking mature miR1511, showed higher and earlier ALS3 response. Root architecture analyses evidenced higher tolerance of G19833 plants to AlT. However, G19833 plants engineered for miR1511 overexpression showed lower ALS3 transcript level and increased sensitivity to AlT. Absence of miR1511 in Andean genotypes, resulting in a diminished ALS3 transcript degradation, appears to be an evolutionary advantage to high Al levels in soils with increased drought conditions.

Visualization of Transiently Expressed mRNA in Plants Using MS2.

RNA transport and localization are evolutionarily conserved processes that allow protein translation to occur at specific subcellular sites and thereby having fundamental roles in the determination of cell fates, embryonic patterning, asymmetric cell division, and cell polarity. In addition to localizing RNA molecules to specific subcellular sites, plants have the ability to exchange RNA molecules between cells through plasmodesmata (PD). Plant RNA viruses hijack the mechanisms of intracellular and intercellular RNA transport to establish localized replication centers within infected cells and then to disseminate their infectious genomes between cells and throughout the plant organism with the help of their movement proteins (MP). In this chapter, we describe the transient expression of the tobacco mosaic virus movement protein (TMV-MP) and the application of the MS2 system for the in vivo labeling of the MP-encoding mRNA. The MS2 method is based on the binding of the bacteriophage coat protein (CP) to its origin of assembly (OAS) in the phage RNA. Thus, to label a specific mRNA in vivo, a tandem repetition of a 19-nucleotide-long stem-loop (SL) sequence derived from the MS2 OAS sequence (MSL) is transcriptionally fused to the RNA under investigation. The RNA is detected by the co-expression of fluorescent protein-tagged MS2 CP (MCP), which binds to each of the MSL elements. In providing a detailed protocol for the in vivo visualization of TMV-MP mRNA tagged with the MS2 system in Nicotiana benthamiana epidermal cells, we describe (1) the specific DNA constructs, (2) Agrobacterium tumefaciens-mediated transfection for their transient expression in plants, and (3) imaging conditions required to obtain high-quality mRNA imaging data.

Identification of miRNAs linked to peanut nodule functional processes.

microRNAs (miRNAs) are non-coding small RNAs that regulate gene expression at post-transcriptional level. Thousands of miRNAs have been identified in legumes, but studies about miRNAs linked to peanut nodule functionality are scarce. In this work we analyzed transcriptional changes in peanut nodules to identify miRNAs involved in functional processes of these organs. We found 32 miRNAs precursors differentially expressed in nodules compared with roots, and predicted the potential targets of their corresponding mature miRNAs. Among them, 20 belong to 14 conserved miRNAs families and 12 are Arachis hypogaea-specific miRNAs. Expression levels of 3 miRNAs (ahy-miR399, ahy-miR159 and ahy-miR3508) were confirmed experimentally by qPCR. We also demonstrated that the expression of these miRNAs was not affected by inoculation of a biocontrol bacterium or a fungal pathogen. The catalogue of differentially expressed miRNA precursors and the expression of the corresponding mature miRNA potential targets in the nodules of A. hypogaea obtained in this work is a database of strong candidates, including A. hypogaea-specific miRNAs, for the regulation of the nodule functionality. The analysis of their role in this process will certainly lead to the characterization of essential regulators in these particular aeschynomenoid nodules.

C4 and crassulacean acid metabolism within a single leaf: deciphering key components behind a rare photosynthetic adaptation.

Although biochemically related, C4 and crassulacean acid metabolism (CAM) systems are expected to be incompatible. However, Portulaca species, including P. oleracea, operate C4 and CAM within a single leaf, and the mechanisms behind this unique photosynthetic arrangement remain largely unknown. Here, we employed RNA-seq to identify candidate genes involved exclusively or shared by C4 or CAM, and provided an in-depth characterization of their transcript abundance patterns during the drought-induced photosynthetic transitions in P. oleracea. Data revealed fewer candidate CAM-specific genes than those recruited to function in C4 . The putative CAM-specific genes were predominantly involved in night-time primary carboxylation reactions and malate movement across the tonoplast. Analysis of gene transcript-abundance regulation and photosynthetic physiology indicated that C4 and CAM coexist within a single P. oleracea leaf under mild drought conditions. Developmental and environmental cues were shown to regulate CAM expression in stems, whereas the shift from C4 to C4 -CAM hybrid photosynthesis in leaves was strictly under environmental control. Moreover, efficient starch turnover was identified as part of the metabolic adjustments required for CAM operation in both organs. These findings provide insights into C4 /CAM connectivity and compatibility, contributing to a deeper understanding of alternative ways to engineer CAM into C4 crop species.

Reprogramming of Root Cells during Nitrogen-Fixing Symbiosis Involves Dynamic Polysome Association of Coding and Noncoding RNAs.

Translational control is a widespread mechanism that allows the cell to rapidly modulate gene expression in order to provide flexibility and adaptability to eukaryotic organisms. We applied translating ribosome affinity purification combined with RNA sequencing to characterize translational regulation of mRNAs at early stages of the nitrogen-fixing symbiosis established between Medicago truncatula and Sinorhizobium meliloti Our analysis revealed a poor correlation between transcriptional and translational changes and identified hundreds of regulated protein-coding and long noncoding RNAs (lncRNAs), some of which are regulated in specific cell types. We demonstrated that a short variant of the lncRNA Trans-acting small interference RNA3 (TAS3) increased its association to the translational machinery in response to rhizobia. Functional analysis revealed that this short variant of TAS3 might act as a target mimic that captures microRNA390, contributing to reduce trans acting small interference Auxin Response Factor production and modulating nodule formation and rhizobial infection. The analysis of alternative transcript variants identified a translationally upregulated mRNA encoding subunit 3 of the SUPERKILLER complex (SKI3), which participates in mRNA decay. Knockdown of SKI3 decreased nodule initiation and development, as well as the survival of bacteria within nodules. Our results highlight the importance of translational control and mRNA decay pathways for the successful establishment of the nitrogen-fixing symbiosis.

A dicistronic precursor encoding miR398 and the legume-specific miR2119 coregulates CSD1 and ADH1 mRNAs in response to water deficit.

Plant microRNAs are commonly encoded in transcripts containing a single microRNA precursor. Processing by DICER-LIKE 1 and associated factors results in the production of a small RNA, followed by its incorporation into an AGO-containing protein complex to guide silencing of an mRNA possessing a complementary target sequence. Certain microRNA loci contain more than one precursor stem-loop structure, thus encoding more than one microRNA in the same transcript. Here, we describe a unique case where the evolutionary conserved miR398a is encoded in the same transcript as the legume-specific miR2119. The dicistronic arrangement found in common bean was also observed in other legumes. In Phaseolus vulgaris, mature miR398 and miR2119 are repressed in response to water deficit, and we demonstrate that both are functional as they target the mRNAs for CSD1 and ADH1, respectively. Our results indicate that the repression of miR398 and miR2119 leads to coordinated up-regulation of CSD1 and ADH1 mRNAs in response to water deficit in common bean and possibly in other legumes. Furthermore, we show that miRNA directed CSD1 and ADH1 mRNAs up-regulation also occurs when common bean plants are exposed to flooding, suggesting that plant redox status and fermentation metabolism must be closely coordinated under different adverse conditions.

An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: insights on nitrogen starvation responses.

Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.

RNA-seq Analysis of Salt-Stressed Versus Non Salt-Stressed Transcriptomes of Chenopodium quinoa Landrace R49.

Quinoa (Chenopodium quinoa Willd.), a model halophytic crop species, was used to shed light on salt tolerance mechanisms at the transcriptomic level. An RNA-sequencing analysis of genotype R49 at an early vegetative stage was performed by Illumina paired-ends method comparing high salinity and control conditions in a time-course pot experiment. Genome-wide transcriptional salt-induced changes and expression profiling of relevant salt-responsive genes in plants treated or not with 300 mM NaCl were analyzed after 1 h and 5 days. We obtained up to 49 million pairs of short reads with an average length of 101 bp, identifying a total of 2416 differentially expressed genes (DEGs) based on the treatment and time of sampling. In salt-treated vs. control plants, the total number of up-regulated and down-regulated genes was 945 and 1471, respectively. The number of DEGs was higher at 5 days than at 1 h after salt treatment, as reflected in the number of transcription factors, which increased with time. We report a strong transcriptional reprogramming of genes involved in biological processes like oxidation-reduction, response to stress and response to abscisic acid (ABA), and cell wall organization. Transcript analyses by real-time RT- qPCR supported the RNA-seq results and shed light on the contribution of roots and shoots to the overall transcriptional response. In addition, it revealed a time-dependent response in the expression of the analyzed DEGs, including a quick (within 1 h) response for some genes, suggesting a "stress-anticipatory preparedness" in this highly salt-tolerant genotype.

Robust increase of leaf size by Arabidopsis thaliana GRF3-like transcription factors under different growth conditions.

An increase in crop yield is essential to reassure food security to meet the accelerating global demand. Several genetic modifications can increase organ size, which in turn might boost crop yield. Still, only in a few cases their performance has been evaluated under stress conditions. MicroRNA miR396 repress the expression of GROWTH-REGULATING FACTOR (GRF) genes that codes for transcription factors that promote organ growth. Here, we show that both Arabidopsis thaliana At-GRF2 and At-GRF3 genes resistant to miR396 activity (rGRF2 and rGRF3) increased organ size, but only rGRF3 can produce this effect without causing morphological defects. Furthermore, introduction of At-rGRF3 in Brassica oleracea can increase organ size, and when At-rGRF3 homologs from soybean and rice are introduced in Arabidopsis, leaf size is also increased. This suggests that regulation of GRF3 activity by miR396 is important for organ growth in a broad range of species. Plants harboring rGRF3 have larger leaves also under drought stress, a condition that stimulates miR396 accumulation. These plants also showed an increase in the resistance to virulent bacteria, suggesting that the size increment promoted by rGRF3 occurs without an obvious cost on plant defenses. Our findings indicate that rGRF3 can increase plant organ size under both normal and stress conditions and is a valuable tool for biotechnological applications.

QTLs and eQTLs mapping related to citrandarins' resistance to citrus gummosis disease.

BACKGROUND: Phytophthora nicotianae Breda de Haan (Phytophthora parasitica Dastur) causes severe damage to citrus crops worldwide. A population of citrandarins was created from the cross between the susceptible parent Citrus sunki Hort. Ex Tan. and the resistant parent Poncirus trifoliata (L.) Raf. cv. Rubidoux, both parents and two reference rootstocks (Rangpur lime and Swingle citrumelo) were grafted in a greenhouse on Rangpur lime. Inoculations were performed at 10 cm and 15 cm above the grafting region and the resulting lesions were evaluated by measuring the lesion length 60 days after inoculation. As control, non-inoculated plants of each genotype were used. In addition, we evaluated the expression of 19 candidate genes involved in citrus defense response 48 h after pathogen infection by quantitative Real-Time PCR (qPCR). We mapped genomic regions of Quantitative Trait Loci (QTLs) and Expression Quantitative Trait Loci (eQTLs) associated with resistance to P. parasitica in the linkage groups (LGs) of the previously constructed maps of C. sunki and P. trifoliata. RESULTS: We found disease severity differences among the generated hybrids, with lesion lengths varying from 1.15 to 11.13 mm. The heritability of the character was 65%. These results indicate that there is a great possibility of success in the selection of resistant hybrids within this experiment. The analysis of gene expression profile demonstrated a great variation of responses regarding the activation of plant defense pathways, indicating that citrandarins have several defense strategies to control oomycete infection. The information of the phenotypic and gene expression data made possible to detect genomic regions associated with resistance. Three QTLs and 84 eQTLs were detected in the linkage map of P. trifoliata, while one QTL and 110 eQTLs were detected in C. sunki. CONCLUSIONS: This is the first study to use eQTLs mapping in the Phytophthora-citrus interaction. Our results from the QTLs and eQTLs mapping allow us to conclude that the resistance of some citrandarins to the infection by P. parasitica is due to a favorable combination of QTLs and eQTLs transmitted by both parents.

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness.

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

Selection and testing of reference genes for accurate RT-qPCR in rice seedlings under iron toxicity.

Reverse Transcription quantitative PCR (RT-qPCR) is a technique for gene expression profiling with high sensibility and reproducibility. However, to obtain accurate results, it depends on data normalization by using endogenous reference genes whose expression is constitutive or invariable. Although the technique is widely used in plant stress analyzes, the stability of reference genes for iron toxicity in rice (Oryza sativa L.) has not been thoroughly investigated. Here, we tested a set of candidate reference genes for use in rice under this stressful condition. The test was performed using four distinct methods: NormFinder, BestKeeper, geNorm and the comparative ΔCt. To achieve reproducible and reliable results, Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines were followed. Valid reference genes were found for shoot (P2, OsGAPDH and OsNABP), root (OsEF-1a, P8 and OsGAPDH) and root+shoot (OsNABP, OsGAPDH and P8) enabling us to perform further reliable studies for iron toxicity in both indica and japonica subspecies. The importance of the study of other than the traditional endogenous genes for use as normalizers is also shown here.

Comment on mRNA-Sequencing Analysis Reveals Transcriptional Changes in Root of Maize Seedlings Treated with Two Increasing Concentrations of a New Biostimulant.

Overpopulation is already a reality, and the need for alternative technologies to meet a continuously increasing food demand has been much discussed around the world. In addition, soil decreasing fertility and desertification are obstacles that we will need to be overcome to increase crop productivity with a much-reduced dependence upon inorganic fertilizers. In this context, protein hydrolysates has emerged as an important strategy to reduce the use of inorganic fertilizers, whose applications as biostimulants for plant growth have shown very promising results.