Construction of small RNA cDNA libraries for deep sequencing

Small RNAs (21-24 nucleotides) including microRNAs (miRNAs) and small interfering RNAs (siRNAs) are potent regulators of gene expression in both plants and animals. Several hundred genes encoding miRNAs and thousands of siRNAs have been experimentally identified by cloning approaches. New sequencing technologies facilitate the identification of these molecules and provide global quantitative expression data in a given biological sample. Here, we describe the methods used in our laboratory to construct small RNA cDNA libraries for high-throughput sequencing using technologies such as MPSS, 454 or SBS.     

Deep sequencing of virus-infected cells reveals HIV-encoded small RNAs

Small virus-derived interfering RNAs (viRNAs) play an important role in antiviral defence in plants, insects and nematodes by triggering the RNA interference (RNAi) pathway. The role of RNAi as an antiviral defence mechanism in mammalian cells has been obscure due to the lack of viRNA detection. Although viRNAs from different mammalian viruses have recently been identified, their functions and possible impact on viral replication remain unknown. To identify viRNAs derived from HIV-1, we used the extremely sensitive SOLiD(TM) 3 Plus System to analyse viRNA accumulation in HIV-1-infected T lymphocytes. We detected numerous small RNAs that correspond to the HIV-1 RNA genome. The majority of these sequences have a positive polarity (98.1%) and could be derived from miRNAs encoded by structured segments of the HIV-1 RNA genome (vmiRNAs). A small portion of the viRNAs is of negative polarity and most of them are encoded within the 3'-UTR, which may represent viral siRNAs (vsiRNAs). The identified vsiRNAs can potently repress HIV-1 production, whereas suppression of the vsiRNAs by antagomirs stimulate virus production. These results suggest that HIV-1 triggers the production of vsiRNAs and vmiRNAs to modulate cellular and/or viral gene expression.     

Diversity of endogenous small non-coding RNAs in Oryza sativa

Small non-coding RNAs play important roles in regulating cell functions by controlling mRNA turnover and translational repression in eukaryotic cells. Here we isolated 162 endogenous small RNA molecules from Oryza sativa, which ranged from 16 to 35 nt in length. Further analysis indicated that they represented a diversity of small RNA molecules, including 17 microRNAs (miRNAs), 30 tiny non-coding RNAs (tncRNAs) and 20 repeat-associated small interfering RNAs (rasiRNAs). Among 17 miRNAs, 13 were novel miRNA candidates and their potential targets were important regulatory genes in the rice genome. We also found that a cluster of small RNAs, including many rasiRNAs, matched to a nuclear DNA fragment that evolutionarily derived from chloroplast. These results demonstrate clearly the existence of distinct types of small RNAs in rice and further suggest that small RNAs may control gene regulation through diverse mechanisms.     

Fatty acid desaturases from the microalga Thalassiosira pseudonana

Analysis of a draft nuclear genome sequence of the diatom Thalassiosira pseudonana revealed the presence of 11 open reading frames showing significant similarity to functionally characterized fatty acid front-end desaturases. The corresponding genes occupy discrete chromosomal locations as determined by comparison with the recently published genome sequence. Phylogenetic analysis showed that two of the T. pseudonana desaturase (Tpdes) sequences grouped with proteobacterial desaturases that lack a fused cytochrome b5 domain. Among the nine remaining gene sequences, temporal expression analysis revealed that seven were expressed in T. pseudonana cells. One of these, TpdesN, was previously characterized as encoding a Delta11-desaturase active on palmitic acid. From the six remaining putative desaturase genes, we report here that three, TpdesI, TpdesO and TpdesK, respectively encode Delta6-, Delta5- and Delta4-desaturases involved in production of the health beneficial polyunsaturated fatty acid DHA (docosahexaenoic acid). Furthermore, we show that one of the remaining genes, TpdesB, encodes a Delta8-sphingolipid desaturase with strong preference for dihydroxylated substrates.     

Cloning and Expression Studies of Novel Small RNAs in Tetraploid Cotton

Small RNAs are a group of non-coding RNAs that downregulate gene expression in a sequence-specific manner to control plant growth and development. The objective of the present study was to clone and characterize several small RNAs in cotton. To identify small RNAs that are involved in the development of cotton bolls and fibers, we generated cDNA libraries from cotton bolls at 13?days post-anthesis from two cotton cultivars, Pima Phy 76 (Gossypium bardadense) and Acala 1517?C99 (Gossypium hirsutum). Screening of these libraries identified eight small RNAs, seven of which have not been reported in other plant species and appear to be absent in the known sequences of other plant species. Their predicted target genes are known to be involved in cotton fiber development. The cloned small RNAs displayed lower and differential expression in the examined boll developmental stages using RT-PCR and quantitative RT-PCR. The genetic polymorphism of the small RNAs at the DNA level was evaluated by miRNA-amplified fragment length polymorphism (AFLP) analysis using primers designed from the small miRNA genes in combination with AFLP primers. Homologous small RNA gene sequences were further isolated using this homology-based genotyping approach, and potential hairpin structures were identified. The results represent a novel method to isolate small including miRNA genes at the RNA and DNA levels in many plant species where genome sequences are not available or expressed sequence tags are limited.     

Characterization of unique small RNA populations from rice grain

Small RNAs (approximately 20 to 24 nucleotides) function as naturally occurring molecules critical in developmental pathways in plants and animals. Here we analyze small RNA populations from mature rice grain and seedlings by pyrosequencing. Using a clustering algorithm to locate regions producing small RNAs, we classified hotspots of small RNA generation within the genome. Hotspots here are defined as 1 kb regions within which small RNAs are significantly overproduced relative to the rest of the genome. Hotspots were identified to facilitate characterization of different categories of small RNA regulatory elements. Included in the hotspots, we found known members of 23 miRNA families representing 92 genes, one trans acting siRNA (ta-siRNA) gene, novel siRNA-generating coding genes and phased siRNA generating genes. Interestingly, over 20% of the small RNA population in grain came from a single foldback structure, which generated eight phased 21-nt siRNAs. This is reminiscent of a newly arising miRNA derived from duplication of progenitor genes. Our results provide data identifying distinct populations of small RNAs, including phased small RNAs, in mature grain to facilitate characterization of small regulatory RNA expression in monocot species.     

MicroRNAs: something important between the genes

Non-coding small endogenous RNAs, of 21-24 nucleotides in length, have recently emerged as important regulators of gene expression in both plants and animals. At least three categories of small RNAs exist in plants: short interfering RNAs (siRNAs) deriving from viruses or transgenes and mediating virus resistance or transgene silencing via RNA degradation; siRNAs deriving from transposons or transgene promoters and controlling transposon and transgene silencing probably via chromatin changes; and microRNAs (miRNAs) deriving from intergenic regions of the genome and regulating the expression of endogenous genes either by mRNA cleavage or translational repression. The disruption of miRNA-mediated regulation causes developmental abnormalities in plants, demonstrating that miRNAs play an important role in the regulation of developmental decisions.     

Specific enrichment of miRNAs in Arabidopsis thaliana infected with Tobacco mosaic virus.

RNA silencing is a broadly conserved machinery and is involved in many biological events. Small RNAs are key molecules in RNA silencing pathway that guide sequence-specific gene regulations and chromatin modifications. The silencing machinery works as an anti-viral defense in virus-infected plants. It is generally accepted that virus-specific small interfering (si) RNAs bind to the viral genome and trigger its cleavage. Previously, we have cloned and obtained sequences of small RNAs from Arabidopsis thaliana infected or uninfected with crucifer Tobacco mosaic virus. MicroRNAs (miRNAs) accumulated to a higher percentage of total small RNAs in the virus-infected plants. This was partly because the viral replication protein binds to the miRNA/miRNA* duplexes. In the present study, we mapped the sequences of small RNAs other than virus-derived siRNAs to the Arabidopsis genome and assigned each small RNA. It was demonstrated that only miRNAs increased as a result of viral infection. Furthermore, some newly identified miRNAs and miRNA candidates were found from the virus-infected plants despite a limited number of examined sequences. We propose that it is advantageous to use virus-infected plants as a source for cloning and identifying new miRNAs.     

Chloroplast genomes of the diatoms <Emphasis Type="Italic">Phaeodactylum tricornutum</Emphasis> and <Emphasis Type="Italic">Thalassiosira pseudonana</Emphasis>: comparison with other plastid genomes of the red lineage

The chloroplast genomes of the pennate diatom Phaeodactylum tricornutum and the centric diatom Thalassiosira pseudonana have been completely sequenced and are compared with those of other secondary plastids of the red lineage: the centric diatom Odontella sinensis, the haptophyte Emiliania huxleyi, and the cryptophyte Guillardia theta. All five chromist genomes are compact, with small intergenic regions and no introns. The three diatom genomes are similar in gene content with 127-130 protein-coding genes, and genes for 27 tRNAs, three ribosomal RNAs and two small RNAs (tmRNA and signal recognition particle RNA). All three genomes have open-reading frames corresponding to ORFs148, 355 and 380 of O. sinensis, which have been assigned the names ycf88, ycf89 and ycf90. Gene order is not strictly conserved, but there are a number of conserved gene clusters showing remnants of red algal origin. The acpP, tsf and psb28 genes appear to be on the way from the plastid to the host nucleus, indicating that endosymbiotic gene transfer is a continuing process.     

Draft Sequencing and Analysis of the Genome of Pufferfish Takifugu flavidus

The pufferfish Takifugu flavidus is an important economic species due to its outstanding flavour and high market value. It has been regarded as an excellent model of genetic study for decades as well. In the present study, three mate-pair libraries of T. flavidus genome were sequenced by the SOLiD 4 next-generation sequencing platform, and the draft genome was constructed with the short reads using an assisted assembly strategy. The draft consists of 50,947 scaffolds with an N50 value of 305.7 kb, and the average GC content was 45.2%. The combined length of repetitive sequences was 26.5 Mb, which accounted for 6.87% of the genome, indicating that the compactness of T. flavidus genome was approximative with that of T. rubripes genome. A total of 1,253 non-coding RNA genes and 30,285 protein-encoding genes were assigned to the genome. There were 132,775 and 394 presumptive genes playing roles in the colour pattern variation, the relatively slow growth and the lipid metabolism, respectively. Among them, genes involved in the microtubule-dependent transport system, angiogenesis, decapentaplegic pathway and lipid mobilization were significantly expanded in the T. flavidus genome. This draft genome provides a valuable resource for understanding and improving both fundamental and applied research with pufferfish in the future.     

Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs

Deep sequencing technologies such as Illumina, SOLiD, and 454 platforms have become very powerful tools in discovering and quantifying small RNAs in diverse organisms. Sequencing small RNA fractions always identifies RNAs derived from abundant RNA species such as rRNAs, tRNAs, snRNA, and snoRNA, and they are widely considered to be random degradation products. We carried out bioinformatic analysis of deep sequenced HeLa RNA and after quality filtering, identified highly abundant small RNA fragments, derived from mature tRNAs that are likely produced by specific processing rather than from random degradation. Moreover, we showed that the processing of small RNAs derived from tRNA^Gln is dependent on Dicer in vivo and that Dicer cleaves the tRNA in vitro.     

Diversity of the cadmium-containing carbonic anhydrase in marine diatoms and natural waters

A recent report of a novel carbonic anhydrase (CDCA1) with Cd as its metal centre in the coastal diatom Thalassiosira weissflogii has led us to search for the occurrence of this Cd enzyme (CDCA) in other marine phytoplankton and in the environment. Using degenerate primers designed from the published sequences from T. weissflogii and a putative sequence in the genome of Thalassiosira pseudonana, we show that CDCA is widespread in diatom species and ubiquitous in the environment. All detected genes share more than 64% amino acid identity with the CDCA of T. pseudonana. Analysis of the amino acid sequence of CDCA shows that the putative Cd binding site resembles that of beta-class carbonic anhydrases (CAs). The prevalence of CAs in diatoms that presumably contain Cd at their active site probably reflects the very low concentration of Zn in the marine environment and the difficulty in acquiring inorganic carbon for photosynthesis. The cdca primers developed in this study should be useful for detecting cdca genes in the field, and studying the conditions under which they are expressed.