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

The early response of oil palm (Elaeis guineensis Jacq.) plants to water deprivation: Expression analysis of miRNAs and their putative target genes, and similarities with the response to salinity stress

Oil palm (Elaeis guineensis Jacq.) is a oilseed crop of great economic importance drastically affected by abiotic stresses. MicroRNAs (miRNAs) play crucial roles in transcription and post-transcription regulation of gene expression, being essential molecules in the response of plants to abiotic stress. To better understand the molecular mechanisms behind the response of young oil palm plants to drought stress, this study reports on the prediction and characterization of miRNAs and their putative target genes in the apical leaf of plants subjected to 14 days of water deprivation. Then, the data from this study were compared to the data from a similar study that focused on salinity stress. Both, the drought-and salt-responsive miRNAs and their putative target genes underwent correlation analysis to identify similarities and dissimilarities among them. Among the 81 identified miRNAs, 29 are specific for oil palm, including two (egu-miR28ds and egu-miR29ds) new ones - described for the first time. As for the expression profile, 62 miRNAs were significantly differentially expressed under drought stress, being five up-regulated (miR396e, miR159b, miR529b, egu-miR19sds, and egu-miR29ds) and 57 down-regulated. Transcription factors, such as MYBs, HOXs, and NF-Ys, were predicted as putative miRNA-target genes in oil palm under water deprivation; making them the most predominant group of such genes. Finally, the correlation analysis study revealed a group of putative target genes with similar behavior under salt and drought stresses. Those genes that are upregulated by these two abiotic stresses encode lncRNAs and proteins linked to stress tolerance, stress memory, modulation of ROS signaling, and defense response regulation to abiotic and biotic stresses. In summary, this study provides molecular evidence for the possible involvement of miRNAs in the drought stress response in oil palm. Besides, it shows that, at the molecular level, there are many similarities in the response of young oil palm plants to these two abiotic stresses.

Expression analysis of miRNAs and their putative target genes confirm a preponderant role of transcription factors in the early response of oil palm plants to salinity stress

BACKGROUND

Several mechanisms regulating gene expression contribute to restore and reestablish cellular homeostasis so that plants can adapt and survive in adverse situations. MicroRNAs (miRNAs) play roles important in the transcriptional and post-transcriptional regulation of gene expression, emerging as a regulatory molecule key in the responses to plant stress, such as cold, heat, drought, and salt. This work is a comprehensive and large-scale miRNA analysis performed to characterize the miRNA population present in oil palm (Elaeis guineensis Jacq.) exposed to a high level of salt stress, to identify miRNA-putative target genes in the oil palm genome, and to perform an in silico comparison of the expression profile of the miRNAs and their putative target genes.

RESULTS

A group of 79 miRNAs was found in oil palm, been 52 known miRNAs and 27 new ones. The known miRNAs found belonged to 28 families. Those miRNAs led to 229 distinct miRNA-putative target genes identified in the genome of oil palm. miRNAs and putative target genes differentially expressed under salinity stress were then selected for functional annotation analysis. The regulation of transcription, DNA-templated, and the oxidation-reduction process were the biological processes with the highest number of hits to the putative target genes, while protein binding and DNA binding were the molecular functions with the highest number of hits. Finally, the nucleus was the cellular component with the highest number of hits. The functional annotation of the putative target genes differentially expressed under salinity stress showed several ones coding for transcription factors which have already proven able to result in tolerance to salinity stress by overexpression or knockout in other plant species.

CONCLUSIONS

Our findings provide new insights into the early response of young oil palm plants to salinity stress and confirm an expected preponderant role of transcription factors - such as NF-YA3, HOX32, and GRF1 - in this response. Besides, it points out potential salt-responsive miRNAs and miRNA-putative target genes that one can utilize to develop oil palm plants tolerant to salinity stress.

Computational Identification and Comparative Analysis of Conserved miRNAs and Their Putative Target Genes in the Juglans regia and J. microcarpa Genomes

MicroRNAs (miRNAs) are important factors for the post-transcriptional regulation of protein-coding genes in plants and animals. They are discovered either by sequencing small RNAs or computationally. We employed a sequence-homology-based computational approach to identify conserved miRNAs and their target genes in Persian (English) walnut, Juglans regia, and its North American wild relative, J. microcarpa. A total of 119 miRNA precursors (pre-miRNAs) were detected in the J. regia genome and 121 in the J. microcarpa genome and miRNA target genes were predicted and their functional annotations were performed in both genomes. In the J. regia genome, 325 different genes were targets; 87.08% were regulated by transcript cleavage and 12.92% by translation repression. In the J. microcarpa genome, 316 different genes were targets; 88.92% were regulated by transcript cleavage and 11.08% were regulated by translation repression. Totals of 1.3% and 2.0% of all resistance gene analogues (RGA) and 2.7% and 2.6% of all transcription factors (TFs) were regulated by miRNAs in the J. regia and J. microcarpa genomes, respectively. Juglans genomes evolved by a whole genome duplication (WGD) and consist of eight pairs of fractionated homoeologous chromosomes. Within each pair, the chromosome that has more genes with greater average transcription also harbors more pre-miRNAs and more target genes than its homoeologue. While only minor differences were detected in pre-miRNAs between the J. regia and J. microcarpa genomes, about one-third of the pre-miRNA loci were not conserved between homoeologous chromosome within each genome. Pre-miRNA and their corresponding target genes showed a tendency to be collocated within a subgenome.

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.

Comparative Study of Withanolide Biosynthesis-Related miRNAs in Root and Leaf Tissues of Withania somnifera

Withania somnifera, popularly known as Indian ginseng, is one of the most important medicinal plants. The plant is well studied in terms of its pharmaceutical activities and genes involved in biosynthetic pathways. However, not much is known about the regulatory mechanism of genes responsible for the production of secondary metabolites. The idea was to identify miRNA transcriptome responsible for the regulation of withanolide biosynthesis, specifically of root and leaf tissues individually. The transcriptome data of in vitro culture of root and leaf tissues of the plant was considered for miRNA identification. A total of 24 and 39 miRNA families were identified in root and leaf tissues, respectively. Out of these, 15 and 27 miRNA families have shown their involvement in different biological functions in root and leaf tissues, respectively. We report here, specific miRNAs and their corresponding target genes for corresponding root and leaf tissues. The target genes have also been analyzed for their role in withanolide metabolism. Endogenous root-miR5140, root-miR159, leaf-miR477, and leaf-miR530 were reported for regulation of withanolide biosynthesis.

Genome-wide analysis of miRNAs and Tasi-RNAs in Zea mays in response to phosphate deficiency

Globally important cereal crop maize provides important nutritions and starch in dietary foods. Low phosphate (LPi) availability in the soil frequently limits the maize quality and yield across the world. Small non-coding RNAs (Snc-RNAs) play crucial roles in growth and adaptation of plants to the environment. Snc-RNAs like microRNAs (miRs) and trans-acting small interfering RNAs (Tasi-Rs) play important functions in posttranscriptional regulation of gene expression, which controls plant development, reproduction, and biotic/abiotic stress responses. In order to identify the miR and Tasi-R alterations in leaf and root of maize in response to sufficient phosphate and LPi at 3LS and 4LS, the snc-RNA population libraries for 0th, 1st, 2nd, 4th, and 8th day were constructed. These libraries were used for genome-wide alignment and RNA-fold analysis for possible prediction of potential miRs and Tasi-Rs. This study reported 174 known and conserved differentially expressed miRs of 27 miR families of maize plant. In addition, leaf and root specific potential novel miRs representing 155 new families were also discovered. Differentially expressed conserved as well as novel miR functions in root and leaf during early stage of Pi starvation were extensively discussed. Leaf and root specific miRs as well as common miRs with their target genes, participating in different biological, cellular, and metabolic processes were explored. Further, four miR390-directed Tasi-Rs which belong to TAS3 gene family along with other orthologs of Tasi-Rs were also identified. Finally, the study provides an insight into the composite regulatory mechanism of miRs in maize in response to Pi deficiency.