Next Article in Journal
Gene Loss and Evolution of the Plastome
Previous Article in Journal
Association between the TPMT*3C (rs1142345) Polymorphism and the Risk of Death in the Treatment of Acute Lymphoblastic Leukemia in Children from the Brazilian Amazon Region
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 10
10.3390/genes11101131
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Identification and Analysis of microRNAs in Chlorella sorokiniana Using High-Throughput Sequencing

by
Siti Nor Ani Azaman
1,2,
Dilan Amila Satharasinghe
3,
Sheau Wei Tan
4,
Norio Nagao
5,
Fatimah M. Yusoff
6 and
Swee Keong Yeap
7,*
1
Centre of Foundation Studies for Agricultural Science, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
2
Aquatic Animal Health and Therapeutics Laboratory (AquaHealth), Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
3
Department of Basic Veterinary Sciences, Faculty of Veterinary Medicine and Animal Science University of Peradeniya, Peradeniya 20400, Sri Lanka
4
Laboratory of Vaccine and Biomolecules (VacBio), Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
5
102 Naname-go, Shinkamigoto-cho, Minami Matsuura-gun, Nagasaki 857-4214, Japan
6
Department of Aquaculture, Faculty of Agriculture, Universiti Putra Malaysia, Serdang, Selangor 43400, Malaysia
7
China-ASEAN College of Marine Sciences, Xiamen University Malaysia, Sepang, 43900 Selangor, Malaysia
*
Author to whom correspondence should be addressed.
Genes 2020, 11(10), 1131; https://doi.org/10.3390/genes11101131
Submission received: 11 August 2020 / Revised: 29 August 2020 / Accepted: 9 September 2020 / Published: 25 September 2020
(This article belongs to the Section Microbial Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Chlorella is a popular microalga with robust physiological and biochemical characteristics, which can be cultured under various conditions. The exploration of the small RNA content of Chlorella could improve strategies for the enhancement of metabolite production from this microalga. In this study, stress was introduced to the Chlorella sorokiniana culture to produce high-value metabolites such as carotenoids and phenolic content. The small RNA transcriptome of C. sorokiniana was sequenced, focusing on microRNA (miRNA) content. From the analysis, 98 miRNAs were identified in cultures subjected to normal and stress conditions. The functional analysis result showed that the miRNA targets found were most often involved in the biosynthesis of secondary metabolites, followed by protein metabolism, cell cycle, and porphyrin and chlorophyll metabolism. Furthermore, the biosynthesis of secondary metabolites such as carotenoids, terpenoids, and lipids was found mostly in stress conditions. These results may help to improve our understanding of regulatory mechanisms of miRNA in the biological and metabolic process of Chlorella species. It is important and timely to determine the true potential of this microalga species and to support the potential for genetic engineering of microalgae as they receive increasing focus for their development as an alternative source of biofuel, food, and health supplements.

1. Introduction

Chlorella is known as one of the best microalgae for the production of high-value metabolites such as carotenoids and lipids [1,2]. Microscopic observation of lipid-droplet and starch-granule accumulation in algal cells proved the potential of Chlorella as an alternative source for biofuel production [3]. Several researchers reported that the manipulation of different environmental factors could affect biochemical production from microalgae [4,5]. For instance, Chlorella sorokiniana showed a several fold increase in carotenoid production when cultured under stress-related conditions, regulated through a carotenoid biosynthesis pathway [6].
Although metabolic transition was identified in the accumulation of pigment, lipid, and starch in Chlorella species [7], there is still a lack of genetic information and differential expression of key metabolic pathways that hamper manipulation and improvement of these microalgae. This suggests that the factors and mechanisms responsible for the differences in metabolite deposition, which have not been elucidated completely, could aid in the development of new strategies to modulate carotenoid and lipid deposition and improve their commercial value.
Small noncoding RNA (ncRNA) genes have been found in numerous organisms, acting as key regulators for development through modulation of the processing, stability, and translation of larger RNAs. Their ability to silence specific genes affects a wide range of biological functions, ranging from gene regulation during embryological development and cell differentiation to genome rearrangement [8,9]. Of the classes of small RNAs, microRNAs (miRNAs) play a pivotal role in controlling gene expressions that are involved in developmental, physiological, or metabolic processes, or stress responses through messenger RNA (mRNA) cleavage or translational repression [10,11]. The choice of control mechanism by miRNA is dependent on the degree of sequence complementary between miRNAs and their target, in which perfect pairing results in target mRNA cleavage, whereas imperfect pairing results in repression. The other mechanism of miRNA-mediated regulation is miRNA-mediated mRNA decay [12,13].
This study is the first step toward profiling the miRNAs that are involved in regulating the production of carotenoids and lipids in Chlorella sorokiniana when cultured under moderate light intensity and glucose addition during nitrogen-limited conditions using a small RNA-sequencing technique. This in turn may elucidate the important metabolic pathways (involving carotenoid, lipid, and starch) present in C. sorokiniana, which could facilitate genetic modification to improve its productivity and reduce the reliance on other model organisms such as Chlamydomonas, since their miRNA expression profiles are different.

2. Materials and Methods

2.1. Samples

Chlorella sorokiniana (NIES-2168) used in this study was obtained from the Marine Biotechnology Lab at the Faculty of Agriculture, Universiti Putra Malaysia, which originated from the National Institute for Environmental Studies (NIES), Japan. The microalgae were cultured in Bold’s Basal Medium (BBM). At the beginning of microalga preculture preparation, the freshly cultured microalgae were first inoculated with an exponentially growing algal culture at 10% (v/v) and continuously illuminated at an intensity of approximately 10 µmol photons·m−2·s−1 with a shaking speed of 30 rpm at 27 °C. The preculture microalgae were allowed to grow until the mid-logarithmic phase, reaching approximately 2.5–3.0 × 106 cells/mL on day 15. Then, microalgal cultures were further divided into two flasks with each flask containing 100 mL of 2 × 106 cells/mL. One flask of cells was then cultured under normal conditions while the other was cultured under stress conditions.
The normal conditions for subsequent culturing were the same as the preculture conditions, whereas the stress culture conditions consisted of higher light intensity than the normal conditions (100 µmol photons·m−2·s−1) and the addition of 2% glucose [1,14]. Microalgae cultured under both conditions were allowed to grow for 7 days. The experiments were conducted in a shaker incubator, and a conical flask was used as the growth chamber; the white fluorescence light source was located above the cultures. All experiments were repeated independently in triplicate, generating six datasets for next-generation sequencing (NGS) analysis. The microalgae were harvested by separating the pellet from the medium via centrifugation at 8000× g for 10 min. The pellet was then flash-frozen using liquid nitrogen and stored at −20 °C until analysis.

2.2. Sample Preparation and Sequencing

Total RNA was extracted using Trizol® (Invitrogen, Waltham, MA, USA) followed by clean-up using the modified RNeasy®Mini Kit (Qiagen, Hilden, Germany) [2]. Briefly, 100 mg of frozen tissue was ground using a prechilled mortar and pestle into a fine powder. The ground tissue sample was mixed with 1 mL of TRIzol reagent and homogenized until no visible debris remained. The homogenized samples were then incubated at room temperature for 5 min. Chloroform (0.2 mL) was added to the homogenate, vortexed vigorously for about 15 s, and incubated at room temperature for 2 to 3 min. After centrifugation, the aqueous layer of the sample was transferred to an RNeasy column from a mini RNA isolation kit (Qiagen, Hilden, Germany) for purification. Residual DNA was eliminated by performing on-column DNase I (Qiagen, Hilden, Germany) digestion at 37 °C for 30 min. The integrity of the extracted RNA was determined by gel electrophoresis, and its concentration was measured using a biospectrometer (Eppendorf BioSpectrometer, Hamburg, Germany).
A Qubit RNA HS (High-Sensitivity) assay kit (Life Technologies, Darmstadt, Germany) was used to determine the concentration of RNA in samples. Briefly, the RNA samples were diluted to a final concentration between 250 pg/µL and 100 ng/µL using the buffer provided before loading onto the Qubit 2.0 Fluorometer (Life Technologies, Darmstadt, Germany). The extracted RNA was analyzed using the Agilent High-Sensitivity RNA assay kit on the 2100 bioanalyzer (Agilent Technologies, Waldbronn, Germany) in complement with the RNA 6000 Nano LabChip kit. The assays were conducted for all normal and stress samples with three biological replicates. Data analysis was performed in accordance with the Agilent protocol.

2.3. Library Preparation

The library for small RNA-sequencing was prepared according to the NEBNext Multiplex Small RNA Library Prep Set for Illumina manual from New England Biolabs (NEB, Herts, UK). The main steps for the preparation of small RNA libraries using this kit were as follows: (1) 3′ ligation, (2) primer hybridization, (3) 5′ ligation, (4) first-strand synthesis, (5) PCR amplification, and (6) size selection. The molarity and size of the libraries were assessed by an Agilent high-sensitivity chip on a 2100 Bioanalyzer (Agilent Technologies, Germany). An individual library was prepared for each sample. Sequencing of 51 bp single-end reads was performed on the Illumina MiSeq at the Laboratory of Immunotherapeutic and Vaccine, Universiti Putra Malaysia.

2.4. Data Analysis

Raw reads obtained from Illumina MiSeq were first examined using FastQC. Reads containing adapter contaminants were removed by the data trimming program in the CLC Genomic Workbench Version 7.5.1 [15]. Only clean reads with a length between 18 nucleotides and 27 nucleotides were selected for further analysis. The clean reads from each C. sorokiniana library were annotated according to the miRBase version 22.1 and noncoding RNA (ncRNA) of Chlamydomonas reinhardtii, Coccomyxa subellipsoidea C-169, Volvox carteri, and Ostreococcus lucimarinus.

2.5. Quantitative PCR (qPCR) Validation

miR156c, miR164a, miR396c, and miR5645d that were identified in sequencing and significantly dysregulated in the C. sorokiniana cultured in stress conditions were subjected to quantitative PCR (qPCR) validation. Primers for the miR156c, miR164a, miR396c, and miR5645d were designed using miRprimer (Table S1) [16]. The primer for α-tubulin was designed using Primer-BLAST and used as the housekeeping gene for normalization. All primers were purchased from Integrated DNA Technologies (Singapore). Three replicates of total RNA extracted from C. sorokiniana cultured under normal and stress conditions were subjected to polyA tailing and reverse transcription using the miScript II RT kit. Briefly, 10 ng of total RNA was added with 5× miScript HiSpec buffer, 10× miScript nucleics mix, miScript reverse-transcriptase mix, and RNase-free water. The final mixture was incubated for 60 min at 37 °C, followed by inactivation of the miScript reverse-transcriptase mix via incubation for 5 min at 95 °C. For housekeeping gene amplification, 1 µg of total RNA was added with Quantiscript Reverse Transcriptase, 5× Quantiscript RT buffer, and RT primer mix. The final mixture was incubated for 15 min at 42 °C, followed by inactivation of Quantiscript reverse transcriptase via incubation for 3 min at 95 °C. The complementary DNA (cDNA) products synthesized using the miScript II RT kit and Quantitect RT kit were then diluted with 90 µL of RNase-free water. Subsequently, the cDNA synthesized from miRNA and mRNA was subjected to a qPCR run using the QuantiNova SYBR Green PCR kit (Qiagen, Venlo, The Netherlands) according to the manufacturer’s protocol [16]. Briefly, cDNA converted from three replicates of the sample was added with 1× QuantiNova SYBR Green PCR master mix (Qiagen, Venlo, The Netherlands), and 0.25 µM of each forward and reverse primer. The final mixture was immediately run using a BioRad CFX96 Real-Time PCR system (BioRad, Hercules, CA, USA) using the set-up of 95 °C for 2 min followed by 40 cycles of 95 °C for 5 s and 60 °C for 10 s. The melting curve was determined using a protocol of 95 °C for 1 min, 55 °C for 1 min, and 55 °C for 10 s at an increment of 0.5 °C per cycle for 80 cycles. The efficiency of the primers for miRNAs and housekeeping genes was tested by qPCR using cDNA synthesized from the stress samples, which were diluted to fivefold serial dilutions (51, 52, 53, 54, 55) using the above qPCR set-up. A no-template control (NTC) and a no--reverse-transcription control (RTC) were included in each run of samples and standard curves. The primer efficiency was calculated using the standard curves, and the results are listed in Table S1. Differential expression (fold changes) of the miRNA between the normal and stress samples was normalized using α-tubulin as a housekeeping gene with the double delta Ct method, and the results were compared between sequencing and qPCR analyses. Statistical analysis was performed using SPSS version 17.0 (IBM Corp, Armonkk, NY, USA). The comparison between normal and stress samples was investigated using Student’s t-tests. A value of p < 0.05 was considered statistically significant.

2.6. Prediction of miRNA Targets and Functional Analysis

Predicted miRNAs were identified using psRNATarget: A Plant Small RNA Target Analysis Server [17] using the Chlamydomonas reinhardtii JGI genome project Phytozome (phytozome v10.0, internal number 281) [18]. Targets were qualified on the basis of the default scoring schema with the following criteria: (a) expectation scores between 0 and 5, (b) unpaired energy score (UPE) ranging from 0.0 to 25.0, (c) high-scoring segment pair (HSP) size of 19 nucleotides, (d) maximum number of transcripts targeted less than or equal to 200 hits, (e) seed region ranging from two to 13 nucleotides, and (f) central mismatch leading to translational inhibition centered on two nucleotides. The information pertaining to transcripts targeted by the miRNA includes the type of inhibition (cleavage or translation repression) and multiplicity (number of sites in the transcript likely to be targeted by the miRNA).
Functional analysis of miRNA targets was carried out using AlgaePath [19]. The query input obtained from psRNATarget results was submitted to the gene group analysis. Information on protein class and the pathways related to miRNA targeting genes were also generated using the same tool.

3. Results

3.1. High-Throughput Small RNA-Sequencing Analysis

Small RNA-sequencing for the normal and stress-induced C. sorokiniana libraries generated approximately 2,440,339 and 3,597,169 reads, respectively. The numbers of cleans reads, annotated samples, and unannotated samples are shown in Table 1. Approximately 96,553 and 248,139 unique reads were produced from normal and stress libraries, respectively. The length distribution of the small RNAs in the normal and stress libraries is shown in Figure 1. The size distribution indicated that the most frequent size length was 20 nt, followed by 21 nt and 19 nt, in both libraries.

3.2. Identification of miRNAs in C. sorokiniana

A Venn diagram was generated to show the number of miRNAs found in normal and stress-induced C. sorokiniana libraries (Figure 2). From 98 miRNAs annotated in the database, 20 miRNAs were identified in the normal C. sorokiniana library, while 59 miRNAs were found in the stress-induced C. sorokiniana library. A total of 19 miRNAs were found to be common in both libraries.

3.3. qPCR Validation of the Selected miRNAs between Stress-Induced and Normal C. sorokiniana

The expressions of miR156c, miR164a, miR396c, and miR5645d of C. sorokiniana cultured under normal and stress conditions were validated by qPCR. Figure 2 shows that miR156c, miR164a, and miR5645d were only detected in the C. sorokiniana cultured under stress conditions. miR396c was detected in samples cultured under both conditions and was fourfold higher in the stress sample. qPCR was able to detect all miRNAs in both normal and stress samples; the expression of these selected microRNAs was upregulated under stress conditions, indicating that these miRNAs were overexpressed under stress conditions (Figure 3).

3.4. Functional Analysis of the Putative Target Genes of miRNAs in C. sorokiniana

A total of 4813 miRNA target genes were predicted and identified for known (or conserved) miRNAs (Supplementary Data 1) using psRNAtarget tools. These target genes were then subjected to functional annotation using AlgaePath, and the results obtained are presented in Table 2. From the results obtained, 14 pathways were identified as significantly enriched (p < 0.05). Functional analysis of the miRNA transcriptome of C. sorokiniana demonstrated that the highest percentage of the miRNA target genes found with significant functional annotation were involved in the biosynthesis of secondary metabolites. Additionally, the second highest percentage of the target genes were involved in purine metabolism, followed by cell cycle, porphyrin and chlorophyll metabolism, aminoacyl-transfer RNA (tRNA) biosynthesis, RNA degradation, and the pentose phosphate pathway. The other important metabolism involved in the regulation of miRNA identified in C. sorokiniana miRNA transcriptome were ether lipid metabolism, starch and sucrose metabolism, carotenoid biosynthesis, glycerolipid metabolism, isoprenoid metabolism (ubiquitin and terpenoid biosynthesis), and others (Table S2). Pathways enriched in different cultivation conditions (normal and stress) are presented in Figure 4. From the figure, more pathways related to the production of secondary metabolites were enriched in stress-induced samples, such as carotenoid biosynthesis, terpenoid backbone biosynthesis, and ubiquinone and other terpenoid–quinone biosynthesis. Meanwhile, different pathways involved in fat and lipid syntheses were observed in C. sorokiniana cultured in normal and stress conditions. For example, ether lipid and glycerophospholipid metabolism were enriched in normal conditions, while α-linolenic acid metabolism, glycerolipid metabolism, and biosynthesis of unsaturated fatty acid pathways were enriched in the stress sample (Table S3).

4. Discussion

4.1. Characteristics of miRNA in Microalgae

Based on Qin et al. [20], most of the well-studied and functionally known miRNAs are conserved and designated with identification numbers ranging from miR156 to miR408. The identified miRNAs are commonly found in a variety of terrestrial plants and nonflowering moss. The miRNAs that are species-specific were identified as nonconserved miRNAs. Nonconserved miRNAs are usually imprecisely processed and weakly expressed, and they lack known functional targets. However, some of the nonconserved miRNAs may exhibit higher expression levels under particular conditions, such as environmental changes [20,21]. However, there is controversy regarding the homology of miRNAs in terrestrial plants and algae [22,23]. A recent review by De [24] indicated that miRNAs involved in key biotic signaling networks of streptophyte algae and terrestrial algae might have similar functions. Furthermore, there were reports on the similarity of target function in green alga Chlamydomonas and animals [16]. Previous findings demonstrated that miRNAs found in animals, plants, and algae share similar biogenesis, function, and turnover pathways in diverse lineages [16,25,26].
The average length of detected miRNA in this study (20 nt) is almost similar to the length of plant miRNA. Generally, the length of miRNA is between 16 and 24 nt, whereby plant miRNA is around 21 nt in length and animal miRNA is a little longer at 22–24 nt. The range of selected miRNA length in this study was based on the range of plant and animal miRNA [27,28]. Previous studies reported that more than 60% of all plant miRNAs are 21 nt in length [8]. After the miRNA identification step, the subsequent step was to infer the roles of the putative target genes, as well as the regulatory networks they are involved in. Basically, miRNA binds to the target mRNA and regulates gene expression by cleaving or inhibiting the translation of the target mRNA. In this study, conserved miRNAs were predicted on the basis of a comparison with the reference library. From the result obtained through identification and functional analysis of target miRNAs, it was revealed that most of the identified miRNAs in both normal and stress-induced libraries of this study were significantly involved in the regulation of secondary metabolite biosynthesis.

4.2. Functional Analysis of Identified miRNA in C. sorokiniana

Several identified miRNAs in this study might be involved in regulating flavonoid biosynthesis, terpenoid biosynthesis, and alkaloid biosynthesis [29]. For instance, miR156 in Arabidopsis thaliana was shown to regulate metabolic flux during the flavonoid biosynthetic pathway by targeting the SQUAMOSA promoter-binding protein-like (SPL) protein. According to previous findings, a higher concentration of miR156 reduced SPL activity and enhanced the anthocyanin biosynthetic genes [30]. Wei et al. [31] demonstrated that the expression of Arabidopsis miR156 in Brassica napus caused an increased level of carotenoid production and reproductive shoot branching, suggesting a link between miR156 expression and carotenoid metabolism. The ability of microalgae to produce certain carotenoids can be explored by investigating miR156 expression under stress-related conditions. In addition, miR172i was identified as targeting mRNAs coding for enzymes of the phenylpropanoid pathway, such as 4-coumarate–CoA ligase, chalcone synthase, caffeoyl-CoA O-methyl transferase, dihydroflavonol 4-reductase, and C,2-hydroxyisoflavanone dehydratase in Podophyllum hexandrum [32]. miR408 identified in the stress-induced library is involved in regulating the alkaloid biosynthetic pathway in the opium poppy, which targets the mRNA of reticuline oxidase-like protein that converts S-reticuline to (S)-scoulerine in benzylisoquinoline alkaloid biosynthesis [33].
It is worth noting that the identification of miR164 in the stress-induced library from this study might be related to the disruption of chlorophyll content. In a previous report, miR164 was found to regulate chlorophyll breakdown during leaf senescence in Arabidopsis thaliana [34,35] by targeting the no apical meristem, ATAF ½ and cup-shaped cotyledon (NAC) transcription factor involved in chlorophyll catabolism. Similarly, in this study, the chlorophyll content of C. sorokiniana grown under stress-induced conditions was reduced by 80% compared to cultures grown under normal conditions [3]. These findings demonstrate a possible correlation between the predicted function of the identified miRNA and the morphological study.
Other interesting miRNAs found in this study were predicted to have similar functions to miRNAs involved in regulating lipid biosynthesis in higher organisms, such as animals and plants. miR148a identified from the stress-induced library, which is highly expressed in hepatic tissue of mouse and human, is responsible for alteration of total cholesterol (TC), low-density-lipoprotein-cholesterol (LDL-C), and triglycerides (TAGs) in blood [36]. It was reported that miR148a is involved in regulating the expression of carnitine palmitoyltransferase 1 (CPT1A), a central regulator of mitochondrial fatty-acid β-oxidation [37]. Meanwhile, in oil-producing plants such as oil palm and rapeseed, miR156, miR159, miR166, miR167, miR172, and miR397 were identified in regulating lipid metabolism [38,39]. miR156 and miR397 in oil palm plants are known to regulate fatty-acid synthesis in plastids and the transport of fatty acyl-ACP from plastids to the endoplasmic reticulum by targeting acetyl-CoA carboxylase 1 (ACC1), long-chain acyl-CoA synthetase 9 (LACS9), LACS4, and enoyl-ACP reductase (ENR). In rapeseed plants, miR156 and miR167 are involved in the regulation of oil content during seed development and maturation by targeting SPL and the adenosine diphosphate ribosylation factor (ARF) family. miR172 from rapeseed plants targets the zinc finger protein associated with fatty-acid synthesis and metabolism. The other miRNAs, such as miR396 identified in the stress-induced library, were reported to target a conserved domain of GRF and the WRC (Trp, Arg, Cys) domain, which contains a functional nuclear localization signal and a putative zinc finger motif [40]. This miRNA plays regulatory roles in plant growth and development, as well as the responses to various environmental stresses [41,42]. This is in agreement with the current study, in which higher numbers of miR396 were identified in the stress-induced library, namely, miR396a, miR396c, miR396d, and miR396o, as compared to only miR396b identified in the normal library.
An understanding of the mechanism of reaction in microalgae to various environmental factors and of the information on RNA-mediated stress regulatory networks offers an innovative path for enhancing the production of metabolites in C. sorokiniana. Suitable stress conditions could stimulate the synthesis and accumulation of lipids and high-value products derived from this microalga. The considerable amount of available omics data from green microalgae could be channeled toward achieving efficient regulation of metabolism via the application of synthetic biology by targeting the regulatory networks rather than deletion or overexpression of enzyme-coding genes [43,44].

5. Conclusions

The present report is the first to identify and profile the expression pattern of miRNA in C. sorokiniana through high-throughput sequencing. A total of 98 miRNAs were identified in the small RNA transcriptome of C. sorokiniana grown under normal and stress conditions. Functional analysis of target miRNAs genes demonstrated that most of the identified miRNAs are involved in the regulation of secondary metabolite biosynthesis. The production of secondary metabolites, such as carotenoid and terpenoid biosynthesis, as well as lipid production, was mostly found in stress samples. These findings provide new insights into the molecular mechanisms involved in the regulation of metabolite production toward stress stimulus in microalgae.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/10/1131/s1; Supplementary Data 1: A total of 4813 miRNA target genes were predicted and identified for known (or conserved) miRNAs; Table S1: MicroRNA primer sequence for qPCR; Table S2: AlgaePath analysis; Table S3: AlgaePath analysis for Normal sample and Stress sample.

Author Contributions

S.N.A.A., S.K.Y., D.A.S., S.W.T., N.N., F.M.Y. conceptualized, designed, and performed the experiments. S.N.A.A and S.K.Y. analyzed the data, wrote the paper, prepared figures and/or tables, and reviewed the drafts of the paper. D.A.S., S.W.T., N.N., and F.M.Y. analyzed the data, contributed reagents/materials/analysis tools, and reviewed drafts of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Xiamen University Malaysia (Xiamen University Malaysia Research Fund; Grant no. XMUMRF/2018-C1/ICAM/0002) and an Institute of Bioscience, Higher Institution Center of Excellence (IBS HICoE) grant from the Ministry of Higher Education, Government of Malaysia (Vote no 6369100).

Acknowledgments

We would like to thank Boo Sok Yee, Lily Boo, Norlaily Mohd Ali, Nancy Liew Woan Charn, and Hyunh Ky for their technical support and guidance with analyzing the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, N.; Wang, Y.; Li, Y.-T.; Huang, J.-C.; Chen, F. Sugar-based growth, astaxanthin accumulation and carotenogenic transcription of heterotrophic Chlorella zofingiensis (Chlorophyta). Process Biochem. 2008, 43, 1288–1292. [Google Scholar] [CrossRef]
  2. Eldh, M.; Lötvall, J.; Malmhäll, C.; Ekström, K. Importance of RNA isolation methods for analysis of exosomal RNA: Evaluation of different methods. Mol. Immunol. 2012, 50, 278–286. [Google Scholar] [CrossRef] [PubMed]
  3. Azaman, S.N.A.; Nagao, N.; Yusoff, F.M.; Tan, S.W.; Yeap, S.K. A comparison of the morphological and biochemical characteristics of Chlorella sorokiniana and Chlorella zofingiensis cultured under photoautotrophic and mixotrophic conditions. PeerJ 2017, 5, e3473. [Google Scholar] [CrossRef]
  4. Juneja, A.; Ceballos, R.M.; Murthy, G.S. Effects of environmental factors and nutrient availability on the biochemical composition of algae for biofuels production: A review. Energies 2013, 6, 4607–4638. [Google Scholar] [CrossRef] [Green Version]
  5. Hu, Q. Environmental Effects on Cell Composition. In Handbook of Microalgal Culture: Biotechnology and Applied Phycology; Richmond, A., Ed.; Blackwell Publishing Ltd.: Oxford, UK, 2004; pp. 83–93. [Google Scholar]
  6. Cordero, B.F.; Obraztsova, I.; Couso, I.; Leon, R.; Vargas, M.A.; Rodriguez, H. Enhancement of lutein production in Chlorella sorokiniana (Chorophyta) by improvement of culture conditions and random mutagenesis. Mar. Drugs 2011, 9, 1607–1624. [Google Scholar] [CrossRef] [Green Version]
  7. Gururani, M.A.; Venkatesh, J.; Tran, L.S.P. Regulation of Photosynthesis during Abiotic Stress-Induced Photoinhibition. Mol. Plant 2015, 8, 1304–1320. [Google Scholar] [CrossRef] [Green Version]
  8. He, L.; Hannon, G.J. MicroRNAs: Small RNAs with a big role in gene regulation. Nat. Rev. Genet. 2004, 5, 522–531. [Google Scholar] [CrossRef]
  9. Cerutti, H.; Ma, X.; Msanne, J.; Repas, T. RNA-mediated silencing in algae: Biological roles and tools for analysis of gene function. Eukaryot. Cell 2011, 10, 1164–1172. [Google Scholar] [CrossRef] [Green Version]
  10. Djami-Tchatchou, A.T.; Sanan-Mishra, N.; Ntushelo, K.; Dubery, I.A. Functional Roles of microRNAs in Agronomically Important Plants—Potential as Targets for Crop Improvement and Protection. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  11. Bartel, D.P. MicroRNAs-Genomics, Biogenesis, Mechanism, and Function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  12. Wu, L.; Fan, J.; Belasco, J.G. MicroRNAs direct rapid deadenylation of mRNA. Proc. Natl. Acad. Sci. USA 2006, 103, 4034–4039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Saraiya, A.A.; Li, W.; Wang, C.C. Transition of a microRNA from Repressing to Activating Translation Depending on the Extent of Base Pairing with the Target. PLoS ONE 2013, 8, e55672. [Google Scholar] [CrossRef]
  14. Ip, P.-F.; Chen, F. Production of astaxanthin by the green microalga Chlorella zofingiensis in the dark. Process Biochem. 2005, 40, 733–738. [Google Scholar] [CrossRef]
  15. Boo, L.; Ho, W.Y.; Ali, N.M.; Yeap, S.K.; Ky, H.; Chan, K.G.; Yin, W.F.; Satharasinghe, D.A.; Liew, W.C.; Tan, S.W.; et al. MiRNA Transcriptome Profiling of Spheroid-Enriched Cells with Cancer Stem Cell Properties in Human Breast MCF-7 Cell Line. Int. J. Biol. Sci. 2016, 12, 427–445. [Google Scholar] [CrossRef]
  16. Yamasaki, T.; Voshall, A.; Kim, E.-J.; Moriyama, E.; Cerutti, H.; Ohama, T. Complementarity to an miRNA seed region is sufficient to induce moderate repression of a target transcript in the unicellular green alga Chlamydomonas reinhardtii. Plant J. 2013, 76, 1045–1056. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Dai, X.; Zhuang, Z.; Zhao, P.X. psRNATarget: A plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018, 46, W49–W54. [Google Scholar] [CrossRef] [Green Version]
  18. Krishnakumar, V.; Hanlon, M.R.; Contrino, S.; Ferlanti, E.S.; Karamycheva, S.; Kim, M.; Rosen, B.D.; Cheng, C.-Y.; Moreira, W.; Mock, S.A.; et al. Araport: The Arabidopsis Information Portal. Nucleic Acids Res. 2015, 43, D1003–D1009. [Google Scholar] [CrossRef] [Green Version]
  19. Zheng, H.-Q.; Chiang-Hsieh, Y.-F.; Chien, C.-H.; Hsu, B.-K.; Liu, T.-L.; Chen, C.-N.; Chang, W.-C. AlgaePath: Comprehensive analysis of metabolic pathways using transcript abundance data from next-generation sequencing in green algae. BMC Genom. 2014, 15, 196. [Google Scholar] [CrossRef] [Green Version]
  20. Qin, Z.; Li, C.; Mao, L.; Wu, L. Novel insights from non-conserved microRNAs in plants. Front. Plant Sci. 2014, 5. [Google Scholar] [CrossRef] [Green Version]
  21. Taylor, R.S.; Tarver, J.E.; Hiscock, S.J.; Donoghue, P.C.J. Evolutionary history of plant microRNAs. Trends Plant Sci. 2014, 19, 175–182. [Google Scholar] [CrossRef]
  22. You, C.; Cui, J.; Wang, H.; Qi, X.; Kuo, L.-Y.; Ma, H.; Gao, L.; Mo, B.; Chen, X. Conservation and divergence of small RNA pathways and microRNAs in land plants. Genome Biol. 2017, 18. [Google Scholar] [CrossRef] [PubMed]
  23. Nozawa, M.; Miura, S.; Nei, M. Origins and Evolution of MicroRNA Genes in Plant Species. Genome Biol. Evol. 2012, 4, 230–239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. de Vries, S.; de Vries, J.; Rose, L.E. The Elaboration of miRNA Regulation and Gene Regulatory Networks in Plant–Microbe Interactions. Genes 2019, 10, 310. [Google Scholar] [CrossRef] [Green Version]
  25. Cui, J.; You, C.; Chen, X. The evolution of microRNAs in plants. Curr. Opin. Plant Biol. 2017, 35, 61–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Mukherjee, K.; Campos, H.; Kolaczkowski, B. Evolution of Animal and Plant Dicers: Early Parallel Duplications and Recurrent Adaptation of Antiviral RNA Binding in Plants. Mol. Biol. Evol. 2013, 30, 627–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lopez-Gomollon, S.; Beckers, M.; Rathjen, T.; Moxon, S.; Maumus, F.; Mohorianu, I.; Moulton, V.; Dalmay, T.; Mock, T. Global discovery and characterization of small non-coding RNAs in marine microalgae. BMC Genom. 2014, 15, 697. [Google Scholar] [CrossRef] [Green Version]
  28. Lou, S.; Sun, T.; Li, H.; Hu, Z. Mechanisms of microRNA-mediated gene regulation in unicellular model alga Chlamydomonas reinhardtii. Biotechnol. Biofuels 2018, 11. [Google Scholar] [CrossRef]
  29. Gupta, O.P.; Karkute, S.G.; Banerjee, S.; Meena, N.L.; Dahuja, A. Contemporary Understanding of miRNA-Based Regulation of Secondary Metabolites Biosynthesis in Plants. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  30. Gou, J.-Y.; Felippes, F.F.; Liu, C.-J.; Weigel, D.; Wang, J.-W. Negative Regulation of Anthocyanin Biosynthesis in Arabidopsis by a miR156-Targeted SPL Transcription Factor. Plant Cell 2011, 23, 1512–1522. [Google Scholar] [CrossRef] [Green Version]
  31. Wei, S.; Gruber, M.Y.; Yu, B.; Gao, M.-J.; Khachatourians, G.G.; Hegedus, D.D.; Parkin, I.A.; Hannoufa, A. Arabidopsis mutant sk156 reveals complex regulation of SPL15 in a miR156-controlled gene network. BMC Plant Biol. 2012, 12, 169. [Google Scholar] [CrossRef] [Green Version]
  32. Biswas, S.; Hazra, S.; Chattopadhyay, S. Identification of conserved miRNAs and their putative target genes in Podophyllum hexandrum (Himalayan Mayapple). Plant Gene 2016, 6, 82–89. [Google Scholar] [CrossRef] [Green Version]
  33. Boke, H.; Ozhuner, E.; Turktas, M.; Parmaksiz, I.; Ozcan, S.; Unver, T. Regulation of the alkaloid biosynthesis by miRNA in opium poppy. Plant Biotechnol. J. 2015, 13, 409–420. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, J.H.; Woo, H.R.; Kim, J.; Lim, P.O.; Lee, I.C.; Choi, S.H.; Hwang, D.; Nam, H.G.J.S. Trifurcate feed-forward regulation of age-dependent cell death involving miR164 in Arabidopsis. Science 2009, 323, 1053–1057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Oda-Yamamizo, C.; Mitsuda, N.; Sakamoto, S.; Ogawa, D.; Ohme-Takagi, M.; Ohmiya, A. The NAC transcription factor ANAC046 is a positive regulator of chlorophyll degradation and senescence in Arabidopsis leaves. Sci. Rep. 2016, 6, 23609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Aryal, B.; Singh, A.K.; Rotllan, N.; Price, N.; Fernández-Hernando, C. MicroRNAs and lipid metabolism. Curr. Opin. Lipidol. 2017, 28, 273–280. [Google Scholar] [CrossRef] [PubMed]
  37. Wagschal, A.; Najafi-Shoushtari, S.H.; Wang, L.; Goedeke, L.; Sinha, S.; Delemos, A.S.; Black, J.C.; Ramírez, C.M.; Li, Y.; Tewhey, R.; et al. Genome-wide identification of microRNAs regulating cholesterol and triglyceride homeostasis. Nat. Med. 2015, 21, 1290–1297. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, Z.; Qiao, Y.; Zhang, J.; Shi, W.; Zhang, J. Genome wide identification of microRNAs involved in fatty acid and lipid metabolism of Brassica napus by small RNA and degradome sequencing. Gene 2017, 619, 61–70. [Google Scholar] [CrossRef] [PubMed]
  39. Zheng, Y.; Chen, C.; Liang, Y.; Sun, R.; Gao, L.; Liu, T.; Li, D. Genome-wide association analysis of the lipid and fatty acid metabolism regulatory network in the mesocarp of oil palm (Elaeis guineensis Jacq.) based on small noncoding RNA sequencing. Tree Physiol. 2019, 39, 356–371. [Google Scholar] [CrossRef]
  40. Van Der Knaap, E. A Novel Gibberellin-Induced Gene from Rice and Its Potential Regulatory Role in Stem Growth. Plant Physiol. 2000, 122, 695–704. [Google Scholar] [CrossRef] [Green Version]
  41. Liu, W.; Zhou, Y.; Li, X.; Wang, X.; Dong, Y.; Wang, N.; Liu, X.; Chen, H.; Yao, N.; Cui, X.; et al. Tissue-Specific Regulation of Gma-miR396 Family on Coordinating Development and Low Water Availability Responses. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  42. Omidbakhshfard, M.A.; Proost, S.; Fujikura, U.; Mueller-Roeber, B. Growth-Regulating Factors (GRFs): A Small Transcription Factor Family with Important Functions in Plant Biology. Mol. Plant 2015, 8, 998–1010. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Fischer, N.; Rochaix, J.D. The flanking regions of PsaD drive efficient gene expression in the nucleus of the green alga Chlamydomonas reinhardtii. Mol. Genet. Genom. 2001, 265, 888–894. [Google Scholar] [CrossRef] [PubMed]
  44. Ferrante, P.; Diener, D.R.; Rosenbaum, J.L.; Giuliano, G. Nickel and low CO2-controlled motility in Chlamydomonas through complementation of a paralyzed flagella mutant with chemically regulated promoters. BMC Plant Biol. 2011, 11, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 11 01131 g001 550
Figure 1. Length distribution of small RNAs in normal and stress-induced C. sorokiniana libraries.
Figure 1. Length distribution of small RNAs in normal and stress-induced C. sorokiniana libraries.
Genes 11 01131 g001
Genes 11 01131 g002 550
Figure 2. Common and specific miRNAs identified in normal and stress-induced C. sorokiniana libraries. Conserved miRNAs were miR156 to miR408 [20]. Others were nonconserved miRNAs. The red miRNAs were selected for quantitative PCR (qPCR) validation.
Figure 2. Common and specific miRNAs identified in normal and stress-induced C. sorokiniana libraries. Conserved miRNAs were miR156 to miR408 [20]. Others were nonconserved miRNAs. The red miRNAs were selected for quantitative PCR (qPCR) validation.
Genes 11 01131 g002
Genes 11 01131 g003 550
Figure 3. Normalized fold changes of miR156c, miR164a, miR396c and miR5645d were quantified using qPCR. * p < 0.05 when compared with the normal sample.
Figure 3. Normalized fold changes of miR156c, miR164a, miR396c and miR5645d were quantified using qPCR. * p < 0.05 when compared with the normal sample.
Genes 11 01131 g003
Genes 11 01131 g004 550
Figure 4. Predicted pathways enriched by genes targeted by the miRNA expressed in normal and stress conditions.
Figure 4. Predicted pathways enriched by genes targeted by the miRNA expressed in normal and stress conditions.
Genes 11 01131 g004
Table
Table 1. Statistic for deep-sequencing results generated from normal and stress-induced Chlorella sorokiniana libraries. miRNA: microRNA.
Table 1. Statistic for deep-sequencing results generated from normal and stress-induced Chlorella sorokiniana libraries. miRNA: microRNA.
Total ReadsPercentageUnique ReadsPercentage
NormalRaw reads2,440,339
Clean reads1,004,55110096,553100
miRNA31910.284220.44
Others2320.02590.06
Unannotated1,001,12899.796,07299.5
StressRaw reads3,597,169
Clean reads1,511,785100248,139100
miRNA66880.447880.32
Others6150.0412150.08
Unannotated1,504,48299.52247,13699.6
Table
Table 2. Functional analysis using Algaepath of known miRNA targets in C. sorokiniana. Results show only significant functional annotation with p < 0.05. tRNA: transfer RNA.
Table 2. Functional analysis using Algaepath of known miRNA targets in C. sorokiniana. Results show only significant functional annotation with p < 0.05. tRNA: transfer RNA.
Pathway IDPathway NameHit Number (Query)Percentage in Queryp-Value
map01110 Biosynthesis of secondary metabolites434.08%0.0015
map00230 Purine metabolism171.61%0.0324
map04110 Cell cycle100.95%0.0077
map00860 Porphyrin and chlorophyll metabolism90.85%0.0018
map00970 Aminoacyl-tRNA biosynthesis80.76%0.0282
map03018 RNA degradation80.76%0.0406
map00030 Pentose phosphate pathway70.66%0.0041
map03050 Proteasome70.66%0.0315
map04712 Circadian rhythm, plant50.47%0.0004
map00450 Selenocompound metabolism50.47%0.0092
map00982 Drug metabolism, cytochrome P45050.47%0.03
map00983 Drug metabolism, other enzymes50.47%0.0362
map00290 Valine, leucine, and isoleucine biosynthesis40.38%0.0155
map00770 Pantothenate and CoA biosynthesis40.38%0.0531
map00565 Ether lipid metabolism30.28%0.0056

Share and Cite

MDPI and ACS Style

Azaman, S.N.A.; Satharasinghe, D.A.; Tan, S.W.; Nagao, N.; Yusoff, F.M.; Yeap, S.K. Identification and Analysis of microRNAs in Chlorella sorokiniana Using High-Throughput Sequencing. Genes 2020, 11, 1131. https://doi.org/10.3390/genes11101131

AMA Style

Azaman SNA, Satharasinghe DA, Tan SW, Nagao N, Yusoff FM, Yeap SK. Identification and Analysis of microRNAs in Chlorella sorokiniana Using High-Throughput Sequencing. Genes. 2020; 11(10):1131. https://doi.org/10.3390/genes11101131

Chicago/Turabian Style

Azaman, Siti Nor Ani, Dilan Amila Satharasinghe, Sheau Wei Tan, Norio Nagao, Fatimah M. Yusoff, and Swee Keong Yeap. 2020. "Identification and Analysis of microRNAs in Chlorella sorokiniana Using High-Throughput Sequencing" Genes 11, no. 10: 1131. https://doi.org/10.3390/genes11101131

APA Style

Azaman, S. N. A., Satharasinghe, D. A., Tan, S. W., Nagao, N., Yusoff, F. M., & Yeap, S. K. (2020). Identification and Analysis of microRNAs in Chlorella sorokiniana Using High-Throughput Sequencing. Genes, 11(10), 1131. https://doi.org/10.3390/genes11101131

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop