Dual RNA-seq of maize and H. seropedicae ZAE94 association, in different doses of nitrate, reveals novel insights into Plant-PGPB-environment relationship.

The interactions between plants, beneficial bacteria and their environment are profoundly shaped by various environmental factors, including light, temperature, water availability, and soil quality. Despite efforts to elucidate the molecular mechanisms involved in the association between plants and beneficial bacteria, like Plant Growth-Promoting Bacteria (PGPB), with many studies focusing on the transcriptional reprogramming in the plant, there is no report on the modulation of genetic controls from both plant and associated bacteria standpoints, in response to environment. The main goal of this study was to investigate the relationship between plant-bacteria-environment signaling, using as a model maize plants inoculated with H. seropedicae ZAE94 and cultivated with different doses of N (0.3 and 3 mM). For this purpose, we performed rRNA-depleted RNA-seq to determine the global gene expression of both maize roots and associated H. seropedicae ZAE94. Our results revealed a differential modulation of maize nitrogen metabolism, phytohormone and cell wall responses when associated with H. seropedicae ZAE94 at different N concentrations. In parallel, a modulation of the bacterial metabolism could be observed, by regulating genes involved in transport, secretion system, cell mobility, oxidoreductases, and chemotaxis, when bacteria were associated with maize roots and cultivated at different doses of N. The molecular and phenotypic data of maize plantlets suggested that different doses of N fertilization differentially regulated the beneficial effects of bacterial inoculation, as higher doses (3 mM) favored shoot elongation and lower doses (0.3 mM) favored increase in plant biomass. Our results provide a valuable integrated overview of differentially expressed genes in both maize and associated H. seropedicae ZAE94 in response to different N availability, revealing new insights into pathways involved in grass-PGPB associations.

The Impact of Non-Nodulating Diazotrophic Bacteria in Agriculture: Understanding the Molecular Mechanisms That Benefit Crops.

Agriculture is facing increasing challenges with regard to achieving sustainable growth in productivity without negatively impacting the environment. The use of bioinoculants is emerging as a sustainable solution for agriculture, especially bioinoculants based on diazotrophic bacteria. Brazil is at the forefront of studies intended to identify beneficial diazotrophic bacteria, as well as in the molecular characterization of this association on both the bacterial and plant sides. Here we highlight the main advances in molecular studies to understand the benefits brought to plants by diazotrophic bacteria. Different molecular pathways in plants are regulated both genetically and epigenetically, providing better plant performance. Among them, we discuss the involvement of genes related to nitrogen metabolism, cell wall formation, antioxidant metabolism, and regulation of phytohormones that can coordinate plant responses to environmental factors. Another important aspect in this regard is how the plant recognizes the microorganism as beneficial. A better understanding of plant-bacteria-environment interactions can assist in the future formulation of more efficient bioinoculants, which could in turn contribute to more sustainable agriculture practices.

From Trash to Luxury: The Potential Role of Plant LncRNA in DNA Methylation During Abiotic Stress.

Remarkable progress has been made in elucidating important roles of plant non-coding RNAs. Among these RNAs, long noncoding RNAs (lncRNAs) have gained widespread attention, especially their role in plant environmental stress responses. LncRNAs act at different levels of gene expression regulation, and one of these mechanisms is by recruitment of DNA methyltransferases or demethylases to regulate the target gene transcription. In this mini-review, we highlight the function of lncRNAs, including their potential role in RNA-directed DNA Methylation (RdDM) silencing pathway and their potential function under abiotic stresses conditions. Moreover, we also present and discuss studies of lncRNAs in crops. Finally, we propose a path outlook for future research that may be important for plant breeding.

A Role for Epigenetic Regulation in the Adaptation and Stress Responses of Non-model Plants.

In recent years enormous progress has been made in understanding the role of epigenetic regulation response to environmental stimuli, especially in response to stresses. Molecular mechanisms involved in chromatin dynamics and silencing have been explained, leading to an appreciation of how new phenotypes can be generated quickly in response to environmental modifications. In some cases, it has also been shown that epigenetic modifications can be stably transmitted to the next generations. Despite this, the vast majority of studies have been carried out with model plants, particularly with Arabidopsis, and very little is known on how native plants in their natural habitat react to changes in their environment. Climate change has been affecting, sometimes drastically, the conditions of numerous ecosystems around the world, forcing populations of native species to adapt quickly. Although part of the adaptation can be explained by the preexisting genetic variation in the populations, recent studies have shown that new stable phenotypes can be generated through epigenetic modifications in few generations, contributing to the stability and survival of the plants in their natural habitat. Here, we review the recent data that suggest that epigenetic variation can help natural populations to cope to with change in their environments.

A miniature inverted-repeat transposable element, AddIn-MITE, located inside a WD40 gene is conserved in Andropogoneae grasses.

Miniature inverted-repeat transposable elements (MITEs) have been associated with genic regions in plant genomes and may play important roles in the regulation of nearby genes via recruitment of small RNAs (sRNA) to the MITEs loci. We identified eight families of MITEs in the sugarcane genome assembly with MITE-Hunter pipeline. These sequences were found to be upstream, downstream or inserted into 67 genic regions in the genome. The position of the most abundant MITE (Stowaway-like) in genic regions, which we call AddIn-MITE, was confirmed in a WD40 gene. The analysis of four monocot species showed conservation of the AddIn-MITE sequence, with a large number of copies in their genomes. We also investigated the conservation of the AddIn-MITE' position in the WD40 genes from sorghum, maize and, in sugarcane cultivars and wild Saccharum species. In all analyzed plants, AddIn-MITE has located in WD40 intronic region. Furthermore, the role of AddIn-MITE-related sRNA in WD40 genic region was investigated. We found sRNAs preferentially mapped to the AddIn-MITE than to other regions in the WD40 gene in sugarcane. In addition, the analysis of the small RNA distribution patterns in the WD40 gene and the structure of AddIn-MITE, suggests that the MITE region is a proto-miRNA locus in sugarcane. Together, these data provide insights into the AddIn-MITE role in Andropogoneae grasses.

Copaifera langsdorffii Novel Putative Long Non-Coding RNAs: Interspecies Conservation Analysis in Adaptive Response to Different Biomes.

Long non-coding RNAs (lncRNAs) are involved in multiple regulatory pathways and its versatile form of action has disclosed a new layer in gene regulation. LncRNAs have their expression levels modulated during plant development, and in response to stresses with tissue-specific functions. In this study, we analyzed lncRNA from leaf samples collected from the legume Copaifera langsdorffii Desf. (copaíba) present in two divergent ecosystems: Cerrado (CER; Ecological Station of Botanical Garden in Brasília, Brazil) and Atlantic Rain Forest (ARF; Rio de Janeiro, Brazil). We identified 8020 novel lncRNAs, and they were compared to seven Fabaceae genomes and transcriptomes, to which 1747 and 2194 copaíba lncRNAs were mapped, respectively, to at least one species. The secondary structures of the lncRNAs that were conserved and differentially expressed between the populations were predicted using in silico methods. A few selected lncRNA were confirmed by RT-qPCR in the samples from both biomes; Additionally, the analysis of the lncRNA sequences predicted that some might act as microRNA (miRNA) targets or decoys. The emerging studies involving lncRNAs function and conservation have shown their involvement in several types of biotic and abiotic stresses. Thus, the conservation of lncRNAs among Fabaceae species considering their rapid turnover, suggests they are likely to have been under functional conservation pressure. Our results indicate the potential involvement of lncRNAs in the adaptation of C. langsdorffii in two different biomes.

PlantRNA_Sniffer: A SVM-Based Workflow to Predict Long Intergenic Non-Coding RNAs in Plants.

Non-coding RNAs (ncRNAs) constitute an important set of transcripts produced in the cells of organisms. Among them, there is a large amount of a particular class of long ncRNAs that are difficult to predict, the so-called long intergenic ncRNAs (lincRNAs), which might play essential roles in gene regulation and other cellular processes. Despite the importance of these lincRNAs, there is still a lack of biological knowledge and, currently, the few computational methods considered are so specific that they cannot be successfully applied to other species different from those that they have been originally designed to. Prediction of lncRNAs have been performed with machine learning techniques. Particularly, for lincRNA prediction, supervised learning methods have been explored in recent literature. As far as we know, there are no methods nor workflows specially designed to predict lincRNAs in plants. In this context, this work proposes a workflow to predict lincRNAs on plants, considering a workflow that includes known bioinformatics tools together with machine learning techniques, here a support vector machine (SVM). We discuss two case studies that allowed to identify novel lincRNAs, in sugarcane (Saccharum spp.) and in maize (Zea mays). From the results, we also could identify differentially-expressed lincRNAs in sugarcane and maize plants submitted to pathogenic and beneficial microorganisms.

Roles of Non-Coding RNA in Sugarcane-Microbe Interaction.

Studies have highlighted the importance of non-coding RNA regulation in plant-microbe interaction. However, the roles of sugarcane microRNAs (miRNAs) in the regulation of disease responses have not been investigated. Firstly, we screened the sRNA transcriptome of sugarcane infected with Acidovorax avenae. Conserved and novel miRNAs were identified. Additionally, small interfering RNAs (siRNAs) were aligned to differentially expressed sequences from the sugarcane transcriptome. Interestingly, many siRNAs aligned to a transcript encoding a copper-transporter gene whose expression was induced in the presence of A. avenae, while the siRNAs were repressed in the presence of A. avenae. Moreover, a long intergenic non-coding RNA was identified as a potential target or decoy of miR408. To extend the bioinformatics analysis, we carried out independent inoculations and the expression patterns of six miRNAs were validated by quantitative reverse transcription-PCR (qRT-PCR). Among these miRNAs, miR408-a copper-microRNA-was downregulated. The cleavage of a putative miR408 target, a laccase, was confirmed by a modified 5'RACE (rapid amplification of cDNA ends) assay. MiR408 was also downregulated in samples infected with other pathogens, but it was upregulated in the presence of a beneficial diazotrophic bacteria. Our results suggest that regulation by miR408 is important in sugarcane sensing whether microorganisms are either pathogenic or beneficial, triggering specific miRNA-mediated regulatory mechanisms accordingly.

Sugarcane transcriptome analysis in response to infection caused by Acidovorax avenae subsp. avenae.

Sugarcane is an important tropical crop mainly cultivated to produce ethanol and sugar. Crop productivity is negatively affected by Acidovorax avenae subsp avenae (Aaa), which causes the red stripe disease. Little is known about the molecular mechanisms triggered in response to the infection. We have investigated the molecular mechanism activated in sugarcane using a RNA-seq approach. We have produced a de novo transcriptome assembly (TR7) from sugarcane RNA-seq libraries submitted to drought and infection with Aaa. Together, these libraries present 247 million of raw reads and resulted in 168,767 reference transcripts. Mapping in TR7 of reads obtained from infected libraries, revealed 798 differentially expressed transcripts, of which 723 were annotated, corresponding to 467 genes. GO and KEGG enrichment analysis showed that several metabolic pathways, such as code for proteins response to stress, metabolism of carbohydrates, processes of transcription and translation of proteins, amino acid metabolism and biosynthesis of secondary metabolites were significantly regulated in sugarcane. Differential analysis revealed that genes in the biosynthetic pathways of ET and JA PRRs, oxidative burst genes, NBS-LRR genes, cell wall fortification genes, SAR induced genes and pathogenesis-related genes (PR) were upregulated. In addition, 20 genes were validated by RT-qPCR. Together, these data contribute to a better understanding of the molecular mechanisms triggered by the Aaa in sugarcane and opens the opportunity for the development of molecular markers associated with disease tolerance in breeding programs.

Computational identification and comparative analysis of miRNA precursors in three palm species.

MAIN CONCLUSION: In the present study, miRNA precursors in the genomes of three palm species were identified. Analyzes of sequence conservation and biological function of their putative targets contribute to understand the roles of miRNA in palm biology. MicroRNAs are small RNAs of 20-25 nucleotides in length, with important functions in the regulation of gene expression. Recent genome sequencing of the palm species Elaeis guineensis, Elaeis oleifera and Phoenix dactylifera have enabled the discovery of miRNA genes, which can be used as biotechnological tools in palm trees breeding. The goal of this study is the identification of miRNA precursors in the genomes of these species and their possible biological roles suggested by the mature miRNA-based regulation of target genes. Mature miRNA sequences from Arabidopsis thaliana, Oryza sativa, and Zea mays available at the miRBase were used to predict microRNA precursors in the palm genomes. Three hundred and thirty-eight precursors, ranging from 76 to 220 nucleotide (nt) in size and distributed in 33 families were identified. Moreover, we also identified 266 miRNA precursors of Musa acuminata, which are phylogenetically close to palms species. To understand the biological function of palm miRNAs, 374 putative miRNA targets were identified. An enrichment analysis of target-gene function was carried out using the agriGO tool. The results showed that the targets are involved in plant developmental processes, mainly regulating root development. Our findings contribute to increase the knowledge on microRNA roles in palm biology and could help breeding programs of palm trees.

Genome-wide identification of microRNA and siRNA responsive to endophytic beneficial diazotrophic bacteria in maize.

BACKGROUND: Small RNA (sRNA) has been described as a regulator of gene expression. In order to understand the role of maize sRNA (Zea mays-hybrid UENF 506-8) during association with endophytic nitrogen-fixing bacteria, we analyzed the sRNA regulated by its association with two diazotrophic bacteria, Herbaspirillum seropedicae and Azospirillum brasilense. RESULTS: Deep sequencing analysis was done with RNA extracted from plants inoculated with H. seropedicae, allowing the identification of miRNA and siRNA. A total of 25 conserved miRNA families and 15 novel miRNAs were identified. A dynamic regulation in response to inoculation was also observed. A hypothetical model involving copper-miRNA is proposed, emphasizing the fact that the up-regulation of miR397, miR398, miR408 and miR528, which is followed by inhibition of their targets, can facilitate association with diazotrophic bacteria. Similar expression patterns were observed in samples inoculated with A. brasilense. Moreover, novel miRNA and siRNA were classified in the Transposable Elements (TE) database, and an enrichment of siRNA aligned with TE was observed in the inoculated samples. In addition, an increase in 24-nt siRNA mapping to genes was observed, which was correlated with an increase in methylation of the coding regions and a subsequent reduction in transcription. CONCLUSION: Our results show that maize has RNA-based silencing mechanisms that can trigger specific responses when plants interact with beneficial endophytic diazotrophic bacteria. Our findings suggest important roles for sRNA regulation in maize, and probably in other plants, during association with diazotrophic bacteria, emphasizing the up-regulation of Cu-miRNA.

Differential sRNA regulation in leaves and roots of sugarcane under water depletion.

Plants have developed multiple regulatory mechanisms to respond and adapt to stress. Drought stress is one of the major constraints to agricultural productivity worldwide and recent reports have highlighted the importance of plant sRNA in the response and adaptation to water availability. In order to increase our understanding of the roles of sRNA in response to water depletion, cultivars of sugarcane were submitted to treatment of ceasing drip irrigation for 24 hours. Deep sequencing analysis was carried out to identify the sRNA regulated in leaves and roots of sugarcane cultivars with different drought sensitivities. The pool of sRNA selected allowed the analysis of different sRNA classes (miRNA and siRNA). Twenty-eight and 36 families of conserved miRNA were identified in leaf and root libraries, respectively. Dynamic regulation of miRNA was observed and the expression profiles of eight miRNA were verified in leaf samples from three biological replicates by stem-loop qRT-PCR assay using the cultivars: SP90-1638--sensitive cultivar--and SP83-2847 and SP83-5073--tolerant cultivars. Altered miRNA regulation was correlated with changes in mRNA levels of specific targets. Two leaf libraries from individual sugarcane cultivars with contrasting drought-tolerance properties were also analyzed. An enrichment of 22-nt sRNA species was observed in leaf libraries. 22-nt miRNA triggered siRNA production by cleavage of their targets in response to water depletion. A number of genes of the sRNA biogenesis pathway were down-regulated in tolerant genotypes and up-regulated in sensitive in response to water depletion treatment. Our analysis contributes to increase the knowledge on the roles of sRNA in sugarcane submitted to water depletion.

High-throughput sequencing of small RNA transcriptome reveals salt stress regulated microRNAs in sugarcane.

Salt stress is a primary cause of crop losses worldwide, and it has been the subject of intense investigation to unravel the complex mechanisms responsible for salinity tolerance. MicroRNA is implicated in many developmental processes and in responses to various abiotic stresses, playing pivotal roles in plant adaptation. Deep sequencing technology was chosen to determine the small RNA transcriptome of Saccharum sp cultivars grown on saline conditions. We constructed four small RNAs libraries prepared from plants grown on hydroponic culture submitted to 170 mM NaCl and harvested after 1 h, 6 hs and 24 hs. Each library was sequenced individually and together generated more than 50 million short reads. Ninety-eight conserved miRNAs and 33 miRNAs* were identified by bioinformatics. Several of the microRNA showed considerable differences of expression in the four libraries. To confirm the results of the bioinformatics-based analysis, we studied the expression of the 10 most abundant miRNAs and 1 miRNA* in plants treated with 170 mM NaCl and in plants with a severe treatment of 340 mM NaCl. The results showed that 11 selected miRNAs had higher expression in samples treated with severe salt treatment compared to the mild one. We also investigated the regulation of the same miRNAs in shoots of four cultivars grown on soil treated with 170 mM NaCl. Cultivars could be grouped according to miRNAs expression in response to salt stress. Furthermore, the majority of the predicted target genes had an inverse regulation with their correspondent microRNAs. The targets encode a wide range of proteins, including transcription factors, metabolic enzymes and genes involved in hormone signaling, probably assisting the plants to develop tolerance to salinity. Our work provides insights into the regulatory functions of miRNAs, thereby expanding our knowledge on potential salt-stressed regulated genes.

Computational identification and analysis of novel sugarcane microRNAs.

BACKGROUND: MicroRNA-regulation of gene expression plays a key role in the development and response to biotic and abiotic stresses. Deep sequencing analyses accelerate the process of small RNA discovery in many plants and expand our understanding of miRNA-regulated processes. We therefore undertook small RNA sequencing of sugarcane miRNAs in order to understand their complexity and to explore their role in sugarcane biology. RESULTS: A bioinformatics search was carried out to discover novel miRNAs that can be regulated in sugarcane plants submitted to drought and salt stresses, and under pathogen infection. By means of the presence of miRNA precursors in the related sorghum genome, we identified 623 candidates of new mature miRNAs in sugarcane. Of these, 44 were classified as high confidence miRNAs. The biological function of the new miRNAs candidates was assessed by analyzing their putative targets. The set of bona fide sugarcane miRNA includes those likely targeting serine/threonine kinases, Myb and zinc finger proteins. Additionally, a MADS-box transcription factor and an RPP2B protein, which act in development and disease resistant processes, could be regulated by cleavage (21-nt-species) and DNA methylation (24-nt-species), respectively. CONCLUSIONS: A large scale investigation of sRNA in sugarcane using a computational approach has identified a substantial number of new miRNAs and provides detailed genotype-tissue-culture miRNA expression profiles. Comparative analysis between monocots was valuable to clarify aspects about conservation of miRNA and their targets in a plant whose genome has not yet been sequenced. Our findings contribute to knowledge of miRNA roles in regulatory pathways in the complex, polyploidy sugarcane genome.

Regulation of miR319 during cold stress in sugarcane.

MicroRNAs (miRNAs) are part of a novel mechanism of gene regulation that is active in plants under abiotic stress conditions. In the present study, 12 miRNAs were analysed to identify miRNAs differentially expressed in sugarcane subjected to cold stress (4 °C). The expression of miRNAs assayed by stem-loop RT-PCR showed that miR319 is up-regulated in sugarcane plantlets exposed to 4 °C for 24 h. The induction of miR319 expression during cold stress was observed in both roots and shoots. Sugarcane miR319 was also regulated by treatment with abscisic acid. Putative targets of this miRNA were identified and their expression levels were decreased in sugarcane plantlets exposed to cold. The cleavage sites of two targets were mapped using a 5' RACE PCR assay confirming the regulation of these genes by miR319. When sugarcane cultivars contrasting in cold tolerance were subjected to 4 °C, we observed up-regulation of miR319 and down-regulation of the targets in both varieties; however, the changes in expression were delayed in the cold-tolerant cultivar. These results suggest that differences in timing and levels of the expression of miR319 and its targets could be tested as markers for selection of cold-tolerant sugarcane cultivars.

In silico identification of coffee genome expressed sequences potentially associated with resistance to diseases.

Sequences potentially associated with coffee resistance to diseases were identified by in silico analyses using the database of the Brazilian Coffee Genome Project (BCGP). Keywords corresponding to plant resistance mechanisms to pathogens identified in the literature were used as baits for data mining. Expressed sequence tags (ESTs) related to each of these keywords were identified with tools available in the BCGP bioinformatics platform. A total of 11,300 ESTs were mined. These ESTs were clustered and formed 979 EST-contigs with similarities to chitinases, kinases, cytochrome P450 and nucleotide binding site-leucine rich repeat (NBS-LRR) proteins, as well as with proteins related to disease resistance, pathogenesis, hypersensitivity response (HR) and plant defense responses to diseases. The 140 EST-contigs identified through the keyword NBS-LRR were classified according to function. This classification allowed association of the predicted products of EST-contigs with biological processes, including host defense and apoptosis, and with molecular functions such as nucleotide binding and signal transducer activity. Fisher's exact test was used to examine the significance of differences in contig expression between libraries representing the responses to biotic stress challenges and other libraries from the BCGP. This analysis revealed seven contigs highly similar to catalase, chitinase, protein with a BURP domain and unknown proteins. The involvement of these coffee proteins in plant responses to disease is discussed.

In silico identification of coffee genome expressed sequences potentially associated with resistance to diseases

Sequences potentially associated with coffee resistance to diseases were identified by in silico analyses using the database of the Brazilian Coffee Genome Project (BCGP). Keywords corresponding to plant resistance mechanisms to pathogens identified in the literature were used as baits for data mining. Expressed sequence tags (ESTs) related to each of these keywords were identified with tools available in the BCGP bioinformatics platform. A total of 11,300 ESTs were mined. These ESTs were clustered and formed 979 EST-contigs with similarities to chitinases, kinases, cytochrome P450 and nucleotide binding site-leucine rich repeat (NBS-LRR) proteins, as well as with proteins related to disease resistance, pathogenesis, hypersensitivity response (HR) and plant defense responses to diseases. The 140 EST-contigs identified through the keyword NBS-LRR were classified according to function. This classification allowed association of the predicted products of EST-contigs with biological processes, including host defense and apoptosis, and with molecular functions such as nucleotide binding and signal transducer activity. Fisher's exact test was used to examine the significance of differences in contig expression between libraries representing the responses to biotic stress challenges and other libraries from the BCGP. This analysis revealed seven contigs highly similar to catalase, chitinase, protein with a BURP domain and unknown proteins. The involvement of these coffee proteins in plant responses to disease is discussed.