# Discovery of MicroRNAs in Cardamom (Elettaria cardamomum Maton) under Drought Stress

Authors: N. Anjali; F. Nadiya; Jinu Thomas; K. K. Sabu
Journal: Dataset Papers in Science
 (2017)

---

## Abstract

MicroRNAs (miRNAs) are endogenous, 19–24-nucleotide-long, noncoding RNAs found in animals, plants, and some viruses. Drought is one of the major limiting factors that negatively affect productivity of plants. To date, nothing is known about the regulatory roles of miRNAs in response to drought stress in cardamom. Cardamom collected from the natural forest area in Therakkudi in the Idamalayar forest range was selected for the study. Total RNA was isolated from leaves and stems of control and drought treated cardamom plant, which was used for library construction. Ion torrent sequencing of the two small RNA libraries prepared from plants raised under well irrigated and drought stressed treatments, respectively, created 3,938,342 and 4,083,181 primary reads. A total of 150 conserved and 20 novel microRNAs were identified from both the control and the treated libraries. Discovery of 17 differentially expressed miRNAs under drought stress suggests that these might have been involved in various biological processes to improve plant tolerance to water stress. This study is the first report of drought responsive miRNAs and their targets in cardamom. The outcome of this research could provide insights into the miRNA based regulatory mechanisms in response to drought stress in monocot plants.


---

## Body

## 1. Introduction

MicroRNAs (miRNAs), which are functionally significant, noncoding RNAs (approximately 19–24 nucleotides (nt) long), act as posttranscriptional gene expression regulators by being involved in different processes like apoptosis, stress response, and differentiation and different disease conditions [1]. Sequencing is the base for awareness about the number, diversity, expression, and probable roles of small noncoding RNAs in plants [2]. Sanger sequencing or chain termination method has been the dominant method of DNA sequencing which was developed by Frederick Sanger in 1977. Genomic research has become revolutionized with the launch of next-generation sequencing (NGS) technologies since 2005. NGS is more economical than Sanger sequencing, and hence researchers are now able to carry out many experiments which were formerly problematic. Cardamom cultivation has good production potential but the plants are vulnerable to many pests, diseases, and abiotic stress factors like droughts, floods, extreme temperatures, salinity, nutrition starvation, and oxidative and heavy metal stress. Drought is one of the major limiting factors that negatively affect the crop productivity of the plant. To resist this drought stress, plants execute several mechanisms at the molecular and physiological levels [3, 4]. Knowledge on the mechanism behind plant response to drought will be useful for better productivity [5].


## 2. Methodology

Wild cardamom (accession number: TBG-C75) collected from the natural forest area in Therakkudi in the Idamalayar forest range (N 10°13′13.20′′  and E 76°47′03.7′′) in Kerala state of India was selected for the study. The collected plants are maintained in a greenhouse with daily watering (voucher specimen deposited in JNTBGRI Herbarium as TBGT86201). Out of those plants, one group was labeled as “control” and another as “treated.” Water was withheld to the plants labeled as “treated” in order to begin the experiment. Stress indications like leaf rolling appeared three weeks after water withholding and at this drought condition the absolute moisture content of the soil was calculated to be <4.5%. Normal watering was given only to the control plant.Small RNA library preparation and next-generation sequencing were as follows. Total RNA was isolated from leaves and stems of control and drought treated cardamom using the combined miRNeasy Mini Kit and CTAB method [6]. Ion total RNA Seq. Kit V2 was used for small RNA sequence library construction following the manufacturer’s instructions after pooling the total RNA from leaf and stem tissues in both control and treated samples. The purified libraries were used for sequencing analysis with the ion torrent sequencer by the Centre for Cellular and Molecular Platforms (C-CAMP, Bangalore).miRNA identification and target prediction were as follows. miRNAs were predicted from the raw small RNA reads using various bioinformatics tools and the pipeline used is depicted in Figure1. Adaptors were removed from the raw sequencing reads using Cutadapt with error rate (-e) set to 0.1 [7]. The remaining reads were checked against snRNAs, snoRNAs, rRNAs, and tRNAs from NCBI database and the perfect matches were eliminated using Bowtie alignment tool [8]. Reads with 17–27 nt length were kept for further analysis. The identical sequences were collapsed into a single sequence by FASTQ/A Collapser tool available in the FASTX-Toolkit. A Blastn search was performed against plant mature miRNAs from miRBase database to identify conserved miRNAs in cardamom small RNA libraries. The remaining reads were mapped onto the transcriptome sequence ofCurcuma longa, the closest relative of cardamom, for the identification of novel miRNAs, as the whole genome or transcriptome sequence of cardamom is not available. The aligned reads were used as input to predict novel miRNAs with the software miRDeep-P (Dataset Item 1 (Table)) [9]. Targets of all the cardamom miRNAs determined in this study were predicted using the psRNATarget (http://bioinfo.noble.org/psRNATarget/) software with default parameters (Dataset Item 2 (Table)) [10]. The model plantArabidopsis thaliana was selected from the input list of cDNA library in the psRNATarget tool as the reference genome. Singular Enrichment Analysis (SEA) tool from agriGO toolkit (http://bioinfo.cau.edu.cn/agriGO/) was used for gene ontology enrichment. Target genes and functional annotation of differentially expressed genes are given in Dataset Item 3 (Table).Figure 1
The pipeline used for predicting miRNAs from raw small RNA reads of cardamom using various bioinformatics tools.Differential expression analysis was performed as follows. The frequency of miRNAs in control and treatment libraries was normalized to transcripts per million (TPM) by the following formula: normalized expression = (actual miRNA count/total count of clean reads) × 1,000,000 (Dataset Item 4 (Table)). Fold change of miRNA expression between treatment and control libraries was determined and miRNAs with higher than 1.5-fold difference in expression levels were considered to be differentially expressed under drought stress. If the normalized expression of miRNA in both the control and the treatment libraries showed a value of less than one, it is removed due to the very low level of expression. Fold change with positive values indicates upregulation and negative values represent downregulation of miRNAs in the drought treated library.


## 3. Dataset Description

The dataset associated with this Dataset Paper consists of 4 items which are described as follows.1
Details of novel miRNAs identified from control and treated small RNA libraries using miRDeep-P tool.
25
36
10
Star Score
Minimum Folding Energy Score
Frequency Score
Score
Flank First End
Flank First Sequence
Flank First Structure
Flank Second Beginning
Flank Second Sequence
Flank Second Structure
Frequency
Loop Beginning
Loop End
Loop Sequence
Loop Structure
Mature Arm
Mature Beginning
Mature End
Mature Query
Mature Sequence
Mature Strand
Mature Structure
Precursor Sequence
Precursor Structure
Primary Beginning
Primary End
Primary Sequence ID
Primary Minimum Folding Energy
Primary Sequence
Primary Structure
Star Arm
Star Beginning
Star End
Star Read
Star Sequence
Star Structure2
Targets of all the cardamom miRNAs determined in this study predicted using the psRNATarget software with default parameters. The column Expectation presents the expectation to score the complementarity between miRNA and their target transcript; the column UPE, the maximum energy to unpair the target site; the column miRNA Start, the sequence number of miRNA from which complementarity between miRNA and target transcript starts; the column miRNA End, the sequence number of miRNA in which complementarity between miRNA and target transcript ends; the column Target Start, the sequence number of target transcript from which complementarity between miRNA and target transcript starts; the column Target End, the sequence number of target transcript in which complementarity between miRNA and target transcript ends; the column Inhibition, the mechanism of gene silencing; and the column Multiplicity, the multiplicity of target site.
2442
13
0
miRNA Sequence
Target mRNA ID
Expectation
UPE
miRNA Start
miRNA End
Target Start
Target End
miRNA Aligned Fragment
Target Aligned Fragment
Inhibition
Target Description
Multiplicity3
Functional annotation of target genes of differentially expressed miRNA families under drought stress. The column miRBase ID presents the miRBase ID of miRNA family and the column Targets presents the target mRNAs and their functional annotation.
13
2
0
miRBase ID
Targets4
The frequency of miRNAs in control and treatment libraries was normalized to transcripts per million (TPM) and fold change of miRNA expression between treatment and control library was determined. The column Fold Change presents the fold change of miRNA expression between treated library and control library.
124
6
0
Name of miRNA
Count of miRNA in Control Library
Count of miRNA in Treated Library
Normalized Reads in Control Library
Normalized Reads in Treated Library
Fold Change


## 4. Concluding Remarks

Small cardamom is affected by many pests, diseases, and abiotic stress factors that have been identified as major production constraints. If these negative factors can be cut down, cardamom production can still go high. This dataset helps in the identification and characterization of drought responsive miRNAs in cardamom, which would eventually be useful in the development of varieties with a lower requirement of water. This dataset allows researchers to develop microsatellite markers from noncoding regions of cardamom genome.



---

*Source: 9507485-2017-07-31.xml*