Skip to main content
  • Research article
  • Open access
  • Published:

Characterization and differential expression of microRNAs elicited by sulfur deprivation in Chlamydomonas reinhardtii



microRNAs (miRNAs) have been found to play an essential role in the modulation of numerous biological processes in eukaryotes. Chlamydomonas reinhardtii is an ideal model organism for the study of many metabolic processes including responses to sulfur-deprivation. We used a deep sequencing platform to extensively profile and identify changes in the miRNAs expression that occurred under sulfur-replete and sulfur-deprived conditions. The aim of our research was to characterize the differential expression of Chlamydomonas miRNAs under sulfur-deprived conditions, and subsequently, the target genes of miRNA involved in sulfur-deprivation were further predicted and analyzed.


By using high-throughput sequencing, we characterized the microRNA transcriptomes under sulphur-replete and sulfur-deprived conditions in Chlamydomonas reinhardtii. We predicted a total of 310 miRNAs which included 85 known miRNAs and 225 novel miRNAs. 13 miRNAs were the specific to the sulfur-deprived conditions. 47 miRNAs showed significantly differential expressions responding to sulfur-deprivation, and most were up-regulated in the small RNA libraries with sulfur-deprivation. Using a web-based integrated system (Web MicroRNAs Designer 3) and combing the former information from a transcriptome of Chlamydomonas reinhardtii, 22 miRNAs and their targets involved in metabolism regulation with sulfur-deprivation were verified.


Our results indicate that sulfur-deprivation may have a significant influence on small RNA expression patterns, and the differential expressions of miRNAs and interactions between miRNA and its targets might further reveal the molecular mechanism responding to sulfur-deprivation in Chlamydomonas reinhardtii.


Sulfur is an essential trace element for all organisms, and is widely used in biochemical processes. Many enzymes and antioxidant molecules such as glutathione contains sulfur. Organically bonded sulfur is a component of all proteins, in the amino acids cysteine and methionine. Generally, sulfate is the most stable form of sulfur. The available pools of sulfate can vary significantly as environmental conditions change. Most organisms have a limited capacity to store sulfur, and thus require different strategies to optimize sulfur use for survival. The ability of microbes to acclimate to periods of nutrient insufficiency is essential to their survival in the natural environment [1]. The unicellular green alga Chlamydomonas reinhardtii is an ideal model organism for the study of many metabolic processes including response to sulfur-deprivation. Sulfur-deprived Chlamydomonas cells have been used for microarray-based RNA abundance studies [2, 3], RNA-seq analysis [4], determination of metabolite profiles [5], and sustained production of H2 [68]. Chlamydomonas reinhardtii exhibits several responses to sulfur deprivation, including changes of the photosynthetic apparatus, the synthesis of enzymes, cell wall structure, SO42- transport activity, and cell size [4].

Recently, there has been considerable interest in understanding the impacts of sulfur deprivation on miRNA. miRNAs were first found in Caenorhabditis elegans through forward genetic screens of the lin-4 and let-7 mutants [9, 10]. Since then, genetic studies on various organisms have revealed that miRNAs are universally present and are key components of various gene regulatory pathways in eukaryotes. Chlamydomonas miRNAs were discovered independently by two groups [11, 12], showing that miRNAs exist not only in multicellular systems but also in unicellular eukaryotes. Experimental approaches and bioinformatics-assisted screening have identified approximately 85 Chlamydomonas reinhardtii miRNAs, which are listed in the miRBase version15.0 To investigate the role of Chlamydomonas reinhardtii miRNAs involved in sulfur deprivation, we used a deep sequencing platform to extensively profile and identify changes in the miRNAs expression that occur during sulfur-replete and sulfur-deprived conditions.


The small RNA profile of Chlamydomonas reinhardtii in both sulfur-replete and sulfur-deprived conditions

The equal numbers of algal cells were resuspended under continuous illumination for up to 72 h in TAP (with sulfate 40.55 mg/L) or TAP-S (the sulfate concentration was less than 0.42 mg/L). Two small RNA libraries were constructed using the algal cells under sulphur-replete (+S library) and with sulfur-deprived conditions (-S library) respectively. Sequencing of the Chlamydomonas reinhardtii small RNA libraries was performed with Solexa high-throughput sequencing. We obtained 11,080,539 reads from the sulfur-deprived (-S) library and 11,284,767 reads from the sulfur-replete (+S) library. After discarding low quality and shorter than 18 (nt) sequences, 9,918,931(-S library) and 10,029,992 (+S library) clean reads ranging from 18 to 30 nt were collected (Table 1). The majority of the small RNA sequences obtained from the -S and + S libraries were 20-25 nt in size, which was the typical size range for Dicer derived products (Figure 1). Comparing the common small RNAs sequences between -S and + S libraries, although the total small RNAs sequences in both libraries reached 98.01%, the common categories of unique small RNAs were only 24.88%. It indicates that sulfur-deprivation led to a significant influence on the small RNA expression patterns in Chlamydomonas reinhardtii (Figure 2).

Table 1 Statistics of small RNA sequences from Chlamydomonas reinhardtii + S and -S libraries
Figure 1
figure 1

The size distribution of small RNAs in both + S (NS) and -S (WS) libraries.

Figure 2
figure 2

Summary of the common and specific tags of small RNAs in + S (sample 1) and -S (sample 3) libraries, including the summary of unique tags (a) and total tags (b). (a) sample1 specific: number of unique sample 1 (+S) specific sRNAs and the percentage; sample3 specific: number of unique sample 3 (-S) specific sRNAs and the percentage; sample1 & sample3: number of unique common sRNAs between two samples and percentage. (b) sample1 specific: number of total sample1 (+S) specific sRNAs and the percentage; sample3 specific: number of total sample3 (-S) specific sRNAs and the percentage; sample1 & sample3: number of total common sRNAs between two samples and percentage.

The 20-24 nt sequences from the + S and -S libraries were aligned to the draft Chlamydomonas reinhardtii genome using SOAP [13]. A total of 4,598,243 (+S) and 6,039,480 (-S) sequences were found to match the genome perfectly (see Additional file 1: Figure S1). These small RNAs were used for further analysis. With the alignment to Genbank, Rfam, Exon and Intron, the composition of small RNAs were annotated as siRNA, miRNA, snRNA, snoRNA etc. (Figure 3). Among all the categories of small RNAs, siRNAs expression (small interfering RNA) did not change on sulfur deprivation. However, the microRNAs expression showed more significant up-regulation in response to sulfur-deprivation (Table 2).

Figure 3
figure 3

Small RNA annotation of + S (a, b) and -S (c, d) libraries. Charts a and c show the numbers of unique sRNA tags matched to all categories of RNA; charts b and d show the number of total sRNA tags matched to all categories of RNA.

Table 2 Composition of the small RNAs in the + S and -S libraries

Identifying novel potential miRNAs in Chlamydomonas reinhardtii

To date, miRBase had a collection of 85 Chlamydomonas miRNA. The characteristic hairpin structure of miRNA precursor can be used to predict novel miRNA. We used the predictive software Mireap to predict novel miRNAs by exploring their secondary structures, and the minimum free energy of the unannotated small RNA tags which could be mapped to the genome. We predicted 225 novel miRNAs by Solexa sequencing and Mireap predictive software in Chlamydomonas reinhardtii. Some of the novel potential miRNAs (n51, n62, n84, n182, n196, Figure 4) have more than ten thousand reads but were not detected by former research [11, 12]. These results showed that high-throughput sequencing of small RNAs is also an ideal strategy to analyze small RNAs profiles and identify novel potential miRNAs in Chlamydomonas reinhardtii.

Figure 4
figure 4

The hairpin secondary structures of 5 predicted miRNA precursors (n51, n62, n84, n182, n196). The positions of mature miRNAs are highlighted in red.

Differential expression of Chlamydomonas reinhardtii miRNAs in sulfur-deprived and sulphur-replete conditions

In order to detect the effect of sulfur-deprivation on Chlamydomonas reinhardtii miRNAs, the expression abundance of 310 candidate miRNAs (minimum short RNA sequence depth above 1) in both -S and + S libraries was examined. Overall, approximately 15% (47 miRNA) of miRNAs showed significant changes in expression (fold changes > 2) upon sulfur-deprivation. Among these miRNA with altered expression, 24 known miRNAs and 23 predicted miRNAs were up-regulated with 13 miRNAs being -S library specific (Tables 3, 4).

Table 3 Known Chlamydomonas reinhardtii miRNAs that are responsive to sulfur-deprivation
Table 4 Predicted Chlamydomonas reinhardtii miRNAs that are responsive to sulfur-deprivation

To confirm the expression changes of Chlamydomonas miRNAs and their response to sulfur-deprived stress, we used quantitative RT-PCR analysis to validate the results of the high throughput sequencing. Fourteen miRNAs, which included 12 known miRNAs and 2 novel predicted miRNAs, were selected at random. The quantitative RT-PCR experimental results of 11 miRNAs matched these of high throughput sequencing data (Figure 5 and Tables 3, 4). However, 3 of the chosen miRNAs did not show the matched results. We deduce that this was likely due to the low quality of primers [14] or low abundance of miRNAs, and further research is needed for this problem.

Figure 5
figure 5

quantitative RT-PCR analysis of fold changes for the expression of fourteen miRNAs between + S(NS) and -S(WS).

Target predictions of Chlamydomonas reinhardtii miRNAs involved sulfur-deprivation

The target genes of Chlamydomonas. reinhardtii miRNAs were predicted by a web-base integrated system, Web MicroRNAs Designer 3 The 47 miRNA, which showed differential expressions (fold changes > 2) to sulfur-deprivation were screened for target predictions. We divided the predicted target genes into different group according to their putative functions including photosynthesis, carbon metabolism, lipid metabolism and other processes. We analyzed the differential expression correlation between the miRNA and the target genes [3, 4] in response to sulfur-deprivation; consequently, we confirmed 17 miRNAs and the target genes involved in response to sulfur-deprivation in Chlamydomonas reinhardtii (Table 5).

Table 5 General overview of Chlamydomonas reinhardtii miRNAs and their functional targets responded to sulfur-deprivation


The reliability of the method to identify Chlamydomonas miRNA

miRNAs are a group of small non-coding RNAs that play an important role in various developmental and stress response processes through negative control of gene expression [15]. miRNA identification had been reported in Chlamydomonas reinhardtii [11, 12], but only 85 miRNAs were annotated and collected to miRBase in previous studies. In this study we characterized a total of 310 Chlamydomonas reinhardtii miRNAs which included the 85 known miRNAs and 225 novel miRNAs by using Solexa sequencing technology with prediction software Mireap. We found that most known miRNAs were abundant, and that the 225 novel miRNAs were less abundant or specific to the miRNA of -S library. It further confirmed Solexa sequencing technology as a high-throughput sequencing system, which was able to produce highly accurate, reproducible and quantitative readouts of small RNAs [16, 17].

The Chlamydomonas miRNAs expression patterns for responding to sulfur-deprivation

The unicellular green alga Chlamydomonas reinhardtii has been used as a model organism for studying the responses of eukaryotes to sulfur-deprivation [4]. The completion of Chlamydomonas reinhardtii genome sequencing and substantially expressed sequence tag (EST) project has greatly increased its utility as a model system [18]. Although some studies have been carried out on Chlamydomonas miRNAs regulation [11, 12], this is the first report on the miRNAs which responded to sulphur-deprivation stress in Chlamydomonas reinhardtii. By comparing the expression levels of two miRNAs libraries(+S library and -S library), we found that approximately 15% of miRNAs expressions showed evident changes (fold changes > 2) with sulfur deprivation (Tables 3, 4). These results further indicate that sulfur-responsive processes are complicated and need a lot of miRNAs to participate in the regulation of gene expression at both transcriptional and post-transcriptional levels [3, 4]. Interestingly, although most miRNA which were reported by several groups [10, 11] were found in this study, we found that no Chlamydomonas miRNAs have identifiable orthologs with miR395s. miR395s play a very important roles in the sulfur-deprived response in higher plants [19, 20] where miR395 regulates sulfate distribution and metabolism in the plant cells. Further, miR395 targets a low-affinity sulfate transporter AST68 (At5g10180) and 3 enzymes in the sulfate assimilation pathway (i.e., ATP sulfurylases [(APS1: At3g22890), (APS3: At4g14680) and (APS4: At5g43780)] [19, 20]. The highly conserved miR395 and the presence of multiple copies of miR395 in diverse plant species is necessary for sulfur homeostasis. The lack of a universally conserved miR395s family in Chlamydomonas reinhardtii suggests that green algae miRNAs may have a different pathway for responding to sulfur-deprivation than higher plants.

Target genes of Chlamydomonas miRNAs involved in sulfur-deprivation

miRNAs negatively regulated their targets by cleavage-induced degradation, and the effects of miRNAs over expression were thus often reflected in decreased mRNA levels of the target gene [21]. Although we found the differential expression of multitudinous Chlamydomonas miRNA between sulfur-replete and sulfur-deprived conditions, it is difficult to confirm their target genes because of the complex mechanism of interaction between miRNAs and their target transcripts was not determined. Zhao et al (2009) selected the precursor of miRNA cre-MIR1162 as backbone to design amiRNAs silencing the MAA7 and RBCS2 genes [22]. Attila et al (2009) selected the precursor of miRNA cre-MIR 1157 as a backbone to efficiently produce amiRNAs targeting COX90, PSY and DCL1 genes [23]. James et al (2010) also used amiRNAs silencing HydA1, HydA2, and Hyd3 genes [24]. This far we have not found the specific bioinformatics tools to predict Chlamydomonas miRNAs targets, so we chose the web-base integrated system (Web MicroRNAs Designer 3) to predict Chlamydomonas miRNAs involved in sulfur-deprivation. The predicted results contained false positive rates similar to other reports on miRNA target prediction [14], so it was necessary to verify these predicted targets. However, it was more difficult to screen the miRNAs targets without the 3' UTR sequence database and abundance information of mRNA expression for response to sulfur-deprivation because most miRNA targeted the 3' UTR of mRNA and expression abundances showed the negative correlative between miRNA and their target mRNA. Fortunately, the Chlamydomonas reinhardtii trancriptome was characterized from sulfur-replete and sulfur-depleted conditions [3, 4], which provided very important information to analyze the interaction between miRNA and its target mRNA. Theoretically, with the sulfur-deprivation stress, the amount of related mRNAs expression will increased, the transcription levels of its target genes should be decreased by miRNA-mediated gene silencing. In fact, we found that expression abundances of most miRNAs showed a negative correlation with their levels of target mRNA under sulfur-deprivation among 22 miRNAs which regulated the metabolic activities (Table 5). Interestingly, the four target genes which encoded low phosphate bleaching protein (LPB), 6-phosphogluconate dehydrogenase (6-PGDH), decarboxylating amino acid transporter and serine/threonine protein kinase were up-regulated while the expressions of the their miRNA including n222, miR1158, n197 and n005 were up. LPB is important for acclimation of Chlamydomonas reinhardtii to phosphorus and sulfur deprivation [25]. 6-PGDH is the first enzyme for pentose phosphate pathway (PPP), which is common for plant responses to abiotic stresses, and serine/threonine protein kinase is required for acclimation of the alga to sulfur deprivation [26]. These four genes are important for Chlamydomonas reinhardtii to survive under sulfur-deprivation, so their miRNA should be down regulated when the algal cell is stressed by sulfur deprivation. However, our results from high-throughput sequencing demonstrate that their regulated miRNAs are increased with sulfur deprivation. These results may indicate that the expression of the 4 genes are regulated by multiple factor, the miRNA regulation may not be the major force. Further experiments are needed to verify this hypothesis.

We were particularly interested in the relative miRNA and its targets for hydrogen bioproduction. The differential expression profile of Chlamydomonas miRNAs led to massive changes in gene expression and metabolism which was closely associated with H2 photo-production. Several miRNAs targeted to genes involved photosynthesis created an anaerobic environment and induced the activity of hydrogenase. The target gene of miR1166.1 encoded the PG phosphate synthase protein. PG is synthesized from cytidyldiphosphate (CDP)-diacylglycerol and glycerol-3-phosphate by the catalytic action of PG-phosphate synthase. Being the only phospholipid in cyanobacteria, PG was required for the accumulation of chlorophyll-protein complexes in the thylakoid membrane and for the normal functioning of photosystem II (PSII) [27]. In the absence of O2, in order to generate ATP, green algae resorted to anaerobic photosynthetic metabolism by miRNA regulation and evolved H2 in the light and consumed endogenous substrates.


We have performed a deep-sequencing analysis of miRNAs in Chlamydomonas reinhardtii and provided a genome-wide, quantitative view of how sulfur-deprivation impacts the expression of small RNAs in Chlamydomonas reinhardtii. Our data confirm 24 known miRNAs and 23 predicted miRNAs with altered expression under sulfur-deprivation, most of which were up-regulated and 13 were -S library specific. Target predictions revealed that a variety of metabolic processes may be affected by changing the expression of miRNAs. Our study has delivered new insights into the role of miRNAs involved in sulfur-deprivation and provided a new approach to understand the biohydrogen production from the small RNA level in Chlamydomonas reinhardtii.


Growth of the algae

Chlamydomonas reinhardtii CC849 were obtained from Chlamydomonas Genetic Centre (c/o Dr. Elizabeth H. Harris, Department of Botany, Duke University, Durham, NC27706, USA). The algal strain was grown in a Tris-Acetate-Phosphate (TAP) medium at 25°C and under continuous cool-white fluorescent lamps(≈200 μmol photons m-2 s-1). To impose S deprivation [6, 8], the liquid cultures were grown into mid-logarithmic phase, algal cells were collected by centrifugation, were washed twice with liquid TAP medium without S (TAP-S, for 1 L of Medium: 2X Filner's Beijernicks Solution 25 ml; 1 M Potassium Phosphate 1 ml; Trace mineral solution 1 ml; Tris-Base 2.42 g; adjust pH to 7.0 by Glacial Acetic Acid. Sulfur-deprivation media (TAP-S) were prepared by replacement of the S-salts by their chloride counterparts). Equal numbers of cells were resuspended in TAP or TAP-S under continuous illumination for up to 72 h, with cell aliquots were collected for RNA isolation, and the sulfate concentration in the supernatant was determined by Dionex ICS-1100 ion chromatogram. Small RNA library construction was carried out as follows: for the + S library, RNA was isolated from the algal cells which were resuspended in sulfur-replete TAP media; for the -S library, RNA was isolated from the algal cells which were resuspended in sulfur-free media TAP-S.

Preparation of total RNA

Total RNA was extracted using Trizol reagent (Invitrogen). The 72 h cells cultured at 25°C in TAP and TAP-S were collected. Total RNA was extracted according to the manufacturer's protocol. The quality of RNA was examined by using an Agilent 2100 Bioanalyzer. The same amount of total RNA was used to construct the two libraries and the samples were prepared in a similar manner. Sequencing of the two libraries was performed on the Illumina's Solexa Sequencer and the samples were run side by side.

Small RNA library construction and high-throughput sequencing

After PAGE purification of small RNA molecules under 30 bases and ligation of a pair of Solexa adaptors to their 5'and 3'ends (Illumina, San Diego, CA. USA), the small RNA molecules were amplified using the adaptor primers for 17 cycles and the fragments around 90 bp (small RNA + adaptors) were isolated with agarose gel. The purified DNA was used directly for cluster generation and sequencing analysis using the Illumina's Solexa Sequencer according to the manufacturer's instructions. The image files generated by the sequencer were processed to produce digital-quality data. The following procedures were performed with Solexa to summarize the data: evaluation of the sequencing quality, calculation of the length distribution of small RNA reads, and filtration of the reads contaminated by rRNA, tRNA, mRNA, snRNA, and snoRNA. Finally, clean reads were compared with a miRBase database (release 15.0).

Bioinformatics analysis

Data quality and length distribution

We eliminated some contaminant reads from the fq file and to obtain the final clean reads and then summarized the length distribution of these clean reads. Normally, length of small RNA is between 18 nt and 30 nt. The length distribution analysis was helpful to see the composition of small RNA samples. For example, miRNA is normally 21 nt or 22 nt, siRNA is 24 nt, and piRNA is 30 nt. The data was processed by the following steps: 1) Elimination of low quality reads (the criteria for this was listed in the explanation of meaning of each row in the result tables); 2) Elimination of reads with 5' primer contaminants; 3) Elimination of reads without 3' primer; 4) Elimination of reads without the insert tag; 5) Elimination of reads with poly A; 6) Elimination of reads shorter than 18 nt; 7) Summarization of the length distribution of the clean reads.

Mapping to genome

We mapped the small RNA tags to the genome by SOAP to analyze their expression and distribution on the genome. Program and Parameters: soap -v 0 -r 2 -s 7 -p 7 -a clean.fa -d ref_genome.fa -o match_genome.soap.

Summary of known miRNA alignment

We aligned the small RNA to the miRNA precursor of corresponding species (using mature miRNA if there was no precursor information of that species in miRBase14.0) to obtain the miRNA count as well as base bias on the first position of identified miRNAs with certain lengths and on each position of all identified miRNAs respectively.

siRNA identification

Small interfering RNA (siRNA) is a 22-24 nt long double-strand RNA, each strand of which is 2 nt longer than the other on the 3' end. According to this structural feature, we aligned tags from clean reads to each other to find sRNAs meeting this criteria. These tags might be potential siRNA candidates. Program and Parameters: Software developed by BGI-tag2siRNA.

Alignment to Genbank

We annotated the small RNA tags with rRNA, scRNA, snoRNA, snRNA and tRNA from Genbank and elimination of matched tags from unannotated tags. Program and Parameters: blastall -p blastn -FF -e 0.01

Alignment to Rfam

Annotation of the small RNA tags with sequences from Rfam and elimination of matched tags from unannotated tags. Program and Parameters: blastall -p blastn -FF -e 0.01

Exon and intron alignment

We aligned small RNA tags to exons and introns of mRNA to find the degraded fragments of mRNA in the small RNA tags. Program and Parameters: Software developed by BGI-overlap

Small RNA annotation

We Summarized all prior alignments and annotation before. In the previous alignment and annotation, some small RNA tags may be mapped to more than one category. To make every unique small RNA map to only one annotation, we followed the following priority rule: rRNAetc(in which Genbank > Rfam) > known miRNA > repeat > exon > intron. Program and Parameters: Software developed by BGI-tag2annotation

Known miRNA expression profile

We aligned small RNA tags to the miRNA precursor/mature miRNA of corresponding species in miRBase14.0. To show detailed information of alignment, including structure of known miRNA precursor, length and count of tags from the sample, etc, click the miRNA id in the left table to see detailed information of that miRNA. Note: Only part of the known miRNA alignment are shown in this report.

Novel miRNA prediction

The characteristic hairpin structure of miRNA precursor was used to predict novel miRNA. We used the predictive software Mireap to predict novel miRNA by exploring the secondary structure, and the minimum free energy of the unannotated small RNA tags which could be mapped to genome. Mireap can be accessed from the following link:

Program and parameters

Software developed by BGI- Mireap Minimal miRNA sequence length (18); Maximal miRNA sequence length (25); Minimal miRNA reference sequence length (20); Maximal miRNA reference sequence length (23); Maximal copy number of miRNAs on reference (20); Maximal free energy allowed for a miRNA precursor (-18 kcal/mol); Maximal space between miRNA and miRNA* (300); Minimal base pairs of miRNA and miRNA* (16); Maximal bulge of miRNA and miRNA* (4); Maximal asymmetry of miRNA/miRNA* duplex (4); Flank sequence length of miRNA precursor (20); ram and Parameters: blastall -p blastn -FF -e 0.01.

Differential expression of known miRNA

Comparison of the known miRNA expression between two samples to determine the differentially expressed miRNA. The procedures are shown as below: (1) Normalize the expression of miRNA in two samples (control and treatment) to get the expression of transcript per million(TPM). Normalization forum:Normalized expression = Actual miRNA count/Total count of clean reads*1000000; (2) Calculate fold-change and P-value from the normalized expression according the Bayesian method developed by Audic and Claverie (1997) [28]. Then generate the log2 ratio plot and scatter plot.

Fold-change forum: Fold change = log 2 (treatment/control)

P-value forum:

p ( x | y ) = N 2 N 1 ( x + y ) ! x ! y ! 1 + N 2 N 1 ( x + y + 1 ) C ( y y min | x ) = y = 0 y y min p ( y | x ) D ( y y max | x ) = y y max p ( y | x )

The N1 and x represent total counts of clean reads and normalized expression level of a given miRNA in sRNA library of + S sample, respectively. The N2 and y represent total counts of clean reads and normalized expression levels of a given miRNA in sRNA library of -S sample, respectively.

Quantitative RT-PCR validation of miRNAs

Sulfur-deprivation procedure is the same as previously described. Both total RNA from 72 h cells cultivated in TAP and TAP-S were isolated respectively for quantitative RT-PCR using the TRIZOL reagent (Invitrogen, Life Technologies, Carlsbad, CA). Low molecular weight RNA and high molecular weight RNA were separated with 4 M LiCl. sRNAs were polyadenylated at 37°C for 60 min in a 50 μl reaction mixture with 1.5 μg of total RNA, 1 mM ATP, 2.5 mM MgCl2, and 4 U poly(A) polymerase (Takara, Japan). Poly (A)-tailed sRNA was recovered by phenol/chloroform extraction and ethanol precipitation. The sRNAs were dissolved, treated with RNase-free DnaseI (Takara, Japan) and reversely transcribed using poly (T) adapter. Real-time PCR was performed using SYBRR Green Real-time PCR Master Mix (Toyobo, Osaka, Japan) and all the primers used were as listed in Additional file 2: Table S1. For each reaction, 1 μL of diluted cDNA (equivalent to 100 pg of total RNA) was mixed with 10 μL of 2 × SYBR green reaction mix (SYBRR Green qRT-PCR Master Mix; Takara, Japan), and 5 pmol of the forward and the reverse primers were added to make a final volume of 20 μL. The conditions for the PCR amplification were as follows: polymerase activation was conducted at 95°C for 30 s; followed by 40 cycles at 95°C for 5 s, 60°C for 31 s. The specificity of the primer amplicons was tested by analysis of a melting curve. The U4 snoRNA was used as a reference gene in the real-time PCR detection of miRNAs. The data was analyzed using the 2--ΔΔCt program, all with an R2 above 0.998. For this quantitative RT-PCR analysis, 3 technical replicates and 2 biological replicates were used.


  1. Grossman A, Takahashi H: Macronutrient utilization by photosynthetic eukaryotes and the fabric of interactions. Annu Rev Plant Physiol Plant Mol Biol. 2001, 52: 163-210. 10.1146/annurev.arplant.52.1.163.

    Article  CAS  PubMed  Google Scholar 

  2. Zhang Z, Shrager J, Jain M, Chang CW, Vallon O, Grossman AR: Insights into the survival of Chlamydomonas reinhardtii during sulfur starvation based on microarray analysis of gene expression. Eukaryot Cell. 2004, 3: 1331-1348. 10.1128/EC.3.5.1331-1348.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Nguyen AV, Thomas-Hall SR, Malnoe A, Timmins M, Mussgnug JH, Rupprecht J, Kruse O, Hankamer B, Schenk PM: The transcriptome of photo-biological hydrogen production induced by sulphur deprivation in the green alga Chlamydomonas reinhardtii. Eukaryot Cell. 2008, 7: 1965-1979. 10.1128/EC.00418-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Gonzalez-Ballester D, Casero D, Cokus S, Pellegrini M, Merchant SS, Grossman AR: RNA-Seq Analysis of Sulfur-Deprived Chlamydomonas Cells Reveals Aspects of Acclimation Critical for Cell Survival. Plant Cell. 2010, 22: 2058-2084. 10.1105/tpc.109.071167.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Bolling C, Fiehn O: Metabolite profiling of Chlamydomonas reinhardtii under nutrient deprivation. Plant Physiol. 2005, 139: 1995-2005. 10.1104/pp.105.071589.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Melis A, Zhang LP, Forestier M, Ghirardi ML, Seibert M: Sustained photobiological hydrogen gas production upon reversible inactivation of oxygen evolution in the green alga Chlamydomonas reinhardtii. Plant Physiol. 2000, 122 (1): 127-135. 10.1104/pp.122.1.127.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Ghirardi ML, Zhang JP, Lee JW, Flynn T, Seibert M, Greenbaum E, Melis A: Microalgae: a green source of renewable H2. Trends Biotechnol. 2000, 18 (12): 506-511. 10.1016/S0167-7799(00)01511-0.

    Article  CAS  PubMed  Google Scholar 

  8. Ghirardi ML, Posewitz MC, Maness PC, Dubini A, Yu J, Seibert M: Hydrogenases and hydrogen photoproduction in oxygenic photosynthetic organisms. Annu Rev Plant Biol. 2007, 58: 71-91. 10.1146/annurev.arplant.58.032806.103848.

    Article  CAS  PubMed  Google Scholar 

  9. Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75: 843-854. 10.1016/0092-8674(93)90529-Y.

    Article  CAS  PubMed  Google Scholar 

  10. Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403: 901-906. 10.1038/35002607.

    Article  CAS  PubMed  Google Scholar 

  11. Zhao T, Li G, Mi S, Li S, Hannon GJ, Wang X-J, Qi Y: A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes Dev. 2007, 21: 1190-1203. 10.1101/gad.1543507.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Molnar A, Schwach F, Studholme DJ, Thuenemann EC, Baulcombe DC: miRNAs control gene expression in the single-cell alga Chlamydomonas reinhardtii. Nature. 2007, 447: 1126-1130. 10.1038/nature05903.

    Article  CAS  PubMed  Google Scholar 

  13. Li RQ, Li YR, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24: 713-714. 10.1093/bioinformatics/btn025.

    Article  CAS  PubMed  Google Scholar 

  14. Wei Y, Chen S, Yang P, Ma Z, Kang L: Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust. Genome Biol. 2009, 10: R6-10.1186/gb-2009-10-1-r6. doi:10.1186/gb-2009-10-1-r6

    Article  PubMed Central  PubMed  Google Scholar 

  15. Bonnet E, Van de Peer Y, Rouze P: The small RNA world of plants. New Phytol. 2006, 171: 451-468. 10.1111/j.1469-8137.2006.01806.x.

    Article  CAS  PubMed  Google Scholar 

  16. Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, Holoch D, Lim C, Tuschl T: Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods. 2008, 44: 3-12. 10.1016/j.ymeth.2007.09.009.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92: 255-264. 10.1016/j.ygeno.2008.07.001.

    Article  CAS  PubMed  Google Scholar 

  18. Shu L, Hu Zhangli: Small silencing RNAs in Chlamydomonas reinhardtii. Minerva Biotecnologica. 2010, 22 (1): 29-37.

    Google Scholar 

  19. Sunkar R, Chinnusamy V, Zhu J, Zhu JK: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12: 301-309. 10.1016/j.tplants.2007.05.001.

    Article  CAS  PubMed  Google Scholar 

  20. Jones-Rhoades MJ, Bartel B, Bartel DP: MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.

    Article  CAS  PubMed  Google Scholar 

  21. Jung J-H, Seo PJ, Park C-M: MicroRNA biogenesis and function in higher plants. Plant Biotechnol Rep. 2009, 3: 111-126. 10.1007/s11816-009-0085-8.

    Article  Google Scholar 

  22. Zhao T, Wang W, Bai X, Qi YJ: Gene silencing by artificial microRNAs in Chlamydomonas. Plant J. 2009, 58 (1): 157-164. 10.1111/j.1365-313X.2008.03758.x.

    Article  CAS  PubMed  Google Scholar 

  23. Molnar A, Bassett A, Thuenemann EC, Schwach F, Karkare S, Ossowski S, Weigel D, Baulcombe DC: Highly specific gene silencing by artificial microRNAs in the unicellular alga Chlamydomonas reinhardtii. Plant J. 2009, 58 (1): 165-174. 10.1111/j.1365-313X.2008.03767.x.

    Article  CAS  PubMed  Google Scholar 

  24. James EG, Attila M, David CB, Janneke B: RNA silencing of hydrogenase(-like) genes and investigation of their physiological roles in the green alga Chlamydomonas reinhardtii. Biochem J. 2010, 431 (3): 345-351.

    Article  Google Scholar 

  25. Chang CW, Moseley JL, Wykoff D, Grossman AR: The LPB1 gene is important for acclimation of Chlamydomonas reinhardtii to phosphorus and sulfur deprivation. Plant Physiol. 2005, 138: 319-329. 10.1104/pp.105.059550.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. González-Ballester D, Pollock SV, Pootakham W, Grossman AR: The central role of a SNRK2 kinase in sulfur deprivation responses. Plant Physiol. 2008, 147: 216-227. 10.1104/pp.108.116137.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Gombos Z, Várkonyi Z, Hagio M, Kovács L, Masamoto K, Itoh S, Wada H: Phosphatidylglycerol requirement for the function of electron acceptor plastoquinone QB in the photosystem II reaction center. Biochemistry. 2002, 41: 3796-3802. 10.1021/bi011884h.

    Article  CAS  PubMed  Google Scholar 

  28. Audic S, Claverie J-M: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.

    CAS  PubMed  Google Scholar 

Download references


We thank the BGI (Beijing Genomics Institute) for their Solexa sequencing. Also thank Prof. Craig Thompson and Dr. W. Lindsey White for language corrections. This work was financially supported by the National Natural Science Foundation of China (Grant No.31070323, 41176106), and Shenzhen Grant Plan for Science and Technology.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Zhangli Hu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

LS contributed to execute experiments and write the manuscript, ZH contributed to design the research plans and write the manuscript. Both authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Figure S1 Mapping of small RNAs in the + S (a) and -S (b) libraries to genome by SOAP. Y axis represents the number of small RNA tags that locate on each chromosome. The numbers of sRNAs on the sense strand of chromosome are positive (shown in blue), and those on the antisence strand of chromosome are negative (shown in red). X axis shows the chromosomes. (PPT 443 KB)

Additional file 2: Table S1 The primers used in this experiment. (DOC 36 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Shu, L., Hu, Z. Characterization and differential expression of microRNAs elicited by sulfur deprivation in Chlamydomonas reinhardtii. BMC Genomics 13, 108 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: