Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome

Understanding the dynamics of eukaryotic transcriptome is essential for studying the complexity of transcriptional regulation and its impact on phenotype. However, comprehensive studies of transcriptomes at single base resolution are rare, even for modern organisms, and lacking for rice. Here, we present the first transcriptome atlas for eight organs of cultivated rice. Using high-throughput paired-end RNA-seq, we unambiguously detected transcripts expressing at an extremely low level, as well as a substantial number of novel transcripts, exons, and untranslated regions. An analysis of alternative splicing in the rice transcriptome revealed that alternative cis-splicing occurred in similar to 33% of all rice genes. This is far more than previously reported. In addition, we also identified 234 putative chimeric transcripts that seem to be produced by trans-splicing, indicating that transcript fusion events are more common than expected. In-depth analysis revealed a multitude of fusion transcripts that might be by-products of alternative splicing. Validation and chimeric transcript structural analysis provided evidence that some of these transcripts are likely to be functional in the cell. Taken together, our data provide extensive evidence that transcriptional regulation in rice is vastly more complex than previously believed.

Statistical correlation between protein secondary structure and messenger RNA stem-loop structure

Anew integrated sequence-structure database, called IADE (Integrated ASTRAL-DSSP-EMBL), incorporating matching mRNA sequence, amino acid sequence, and protein secondary structural data, is constructed. It includes 648 protein domains. Based on the IADE database, we studied the relation between RNA stem-loop frequencies and protein secondary structure. It was found that the alpha-helices and beta-strands on proteins tend to be preferably "coded" by mRNA stem region, while the coils on proteins tend to be preferably "coded" by mRNA loop region. These tendencies are more obvious if we observe the structural words (SWs). An SW is defined by a four-amino-acid-fragment that shows the pronounced secondary structural (alpha-helix or beta-strand) propensity. It is demonstrated that the deduced correlation between protein and mRNA structure can hardly be explained as the stochastic fluctuation effect. (C) 2003 Wiley Periodicals, Inc.

Origin and evolution of the RIG-I like RNA helicase gene family

Background: The DExD/H domain containing RNA helicases such as retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5) are key cytosolic pattern recognition receptors (PRRs) for detecting nucleotide pathogen associated molecular patterns (PAMPs) of invading viruses. The RIG-I and MDA5 proteins differentially recognise conserved PAMPs in double stranded or single stranded viral RNA molecules, leading to activation of the interferon system in vertebrates. They share three core protein domains including a RNA helicase domain near the C terminus (HELICc), one or more caspase activation and recruitment domains (CARDs) and an ATP dependent DExD/H domain. The RIG-I/MDA5 directed interferon response is negatively regulated by laboratory of genetics and physiology 2 (LGP2) and is believed to be controlled by the mitochondria antiviral signalling protein (MAVS), a CARD containing protein associated with mitochondria. Results: The DExD/H containing RNA helicases including RIG-I, MDA5 and LGP2 were analysed in silico in a wide spectrum of invertebrate and vertebrate genomes. The gene synteny of MDA5 and LGP2 is well conserved among vertebrates whilst conservation of the gene synteny of RIG-I is less apparent. Invertebrate homologues had a closer phylogenetic relationship with the vertebrate RIG-Is than the MDA5/LGP2 molecules, suggesting the RIG-I homologues may have emerged earlier in evolution, possibly prior to the appearance of vertebrates. Our data suggest that the RIG-I like helicases possibly originated from three distinct genes coding for the core domains including the HELICc, CARD and ATP dependent DExD/H domains through gene fusion and gene/domain duplication. Furthermore, presence of domains similar to a prokaryotic DNA restriction enzyme III domain (Res III), and a zinc finger domain of transcription factor (TF) IIS have been detected by bioinformatic analysis. Conclusion: The RIG-I/MDA5 viral surveillance system is conserved in vertebrates. The RIG-I like helicase family appears to have evolved from a common ancestor that originated from genes encoding different core functional domains. Diversification of core functional domains might be fundamental to their functional divergence in terms of recognition of different viral PAMPs.

Isolation and characterization of Argonaute 2: A key gene of the RNA interference pathway in the rare minnow, Gobiocypris rarus

Argonaute 2 gene plays a pivotal role in RNAi in many species. Herein is the first report of the cloning and characterization of Argonaute 2 gene in fish. The full-length cDNA of Gobiocypris rarus Argonaute 2 (GrAgo2) consisted of 3073 nucleotides encoding 869 amino acid residues with a calculated molecular weight of 98.499 kDa and an estimated isoelectric point of 9.18. Analysis of the deduced amino acid sequence showed the presence of two signature domains, PAZ and Piwi. RT-PCR analysis indicated that GrAgo2 mRNA expression could be detected in widespread tissues. After infection with grass carp reovirus, GrAgo2 expression was up-regulated from 12 h post-injection (p < 0.05) and returned to control levels at 48 h post-injection (p > 0.05). These data imply that GrAgo2 is involved in antiviral defense in rare minnow. (C) 2008 Published by Elsevier Ltd.

Functional domains and the antiviral effect of the double-stranded RNA-dependent protein kinase PKR from Paralichthys olivaceus

The double-stranded RNA (dsRNA)-dependent protein kinase PKR is thought to mediate a conserved antiviral pathway by inhibiting viral protein synthesis via the phosphorylation of the alpha subunit of eukaryotic initiation factor 2 (eIF2 alpha). However, little is known about the data related to the lower vertebrates, including fish. Recently, the identification of PKR-like, or PKZ, has addressed the question of whether there is an orthologous PKR in fish. Here, we identify the first fish PKR gene from the Japanese flounder Paralichthys olivaceus (PoPKR). PoPKR encodes a protein that shows a conserved structure that is characteristic of mammalian PKRs, having both the N-terminal region for dsRNA binding and the C-terminal region for the inhibition of protein translation. The catalytic activity of PoPKR is further evidence that it is required for protein translation inhibition in vitro. PoPKR is constitutively transcribed at low levels and is highly induced after virus infection. Strikingly, PoPKR overexpression increases eIF2 alpha phosphorylation and inhibits the replication of Scophthalmus maximus rhabdovirus (SMRV) in flounder embryonic cells, whereas phosphorylation and antiviral effects are impaired in transfected cells expressing the catalytically inactive PKR-K421R variant, indicating that PoPKR inhibits virus replication by phosphorylating substrate eIF2 alpha. The interaction between PoPKR and eIF2 alpha is demonstrated by coimmunoprecipitation assays, and the transfection of PoPKR-specific short interfering RNA further reveals that the enhanced eIF2 alpha phosphorylation is catalyzed by PoPKR during SMRV infection. The current data provide significant evidence for the existence of a PKR-mediated antiviral pathway in fish and reveal considerable conservation in the functional domains and the antiviral effect of PKR proteins between fish and mammals.

The cytomegalovirus promoter-driven short hairpin RNA constructs mediate effective RNA interference in zebrafish in vivo

The ability to utilize the RNA interference (RNAi) machinery for silencing target-gene expression has created a lot of excitement in the research community. In the present study, we used a cytomegalovirus (CMV) promoter-driven DNA template approach to induce short hairpin RNA (shRNA) triggered RNAi to block exogenous Enhanced Green Fluorescent Protein (EGFP) and endogenous No Tail (NTL) gene expressions. We constructed three plasmids, pCMV-EGFP-CMV-shGFP-SV40, pCMV-EGFP-CMV-shNTL-SV40, and pCMV-EGFP-CMV-shScrambled-SV40, each containing a CMV promoter driving an EGFP reporter cDNA and DNA coding for one shRNA under the control of another CMV promoter. The three shRNA-generating plasmids and pCMV-EGFP control plasmid were introduced into zebrafish embryos by microinjection. Samples were collected at 48 h after injection. Results were evaluated by phenotype observation and real-time fluorescent quantitative reverse-transcription polymerase chain reaction (Q-PCR). The shGFP-generating plasmid significantly inhibited the EGFP expression viewed under fluorescent microscope and reduced by 70.05 +/- 1.26% of exogenous EGFP gene mRNA levels compared with controls by Q-PCR. The shRNA targeting endogenous NTL gene resulted in obvious NTL phenotype of 30 +/- 4% and decreased the level of their corresponding mRNAs up to 54.52 +/- 2.05% compared with nontargeting control shRNA. These data proved the feasibility of the CMV promoter-driven shRNA expression technique to be used to inhibit exogenous and endogenous gene expressions in zebrafish in vivo.

Redescription of Trypanosoma siniperca Chang 1964 from freshwater fish of China based on morphological and molecular data

During the parasite fauna investigation within 2005, the freshwater fish trypanosome Trypanosoma siniperca Chang 1964 was isolated from the blood of Mandarin carp (Siniperca chuatsi) from Niushan Lake, Hubei Province, central China. Blood trypomastigotes were observed only, and the density of infection was low. Light microscopy examinations of this material made it possible to study in detail the morphology of this parasite and redescribe it according to current standards. T. siniperca is characterized also on the molecular level using the sequences of SSU rRNA gene. Phylogenetic analyses based on these sequences allowed clearer phylogenetic relationships to be established with other fish trypanosomes sequenced to date.

Efficient RNA interference in zebrafish embryos using siRNA synthesized with SP6 RNA polymerase

Double-stranded RNA (dsRNA) has been shown to be a useful tool for silencing genes in zebrafish (Danio rerio), while the blocking specificity of dsRNA is still of major concern for application. It was reported that siRNA (small interfering RNA) prepared by endoribonuclease digestion (esiRNA) could efficiently silence endogenous gene expression in mammalian embryos. To test whether esiRNA could work in zebrafish, we utilized Escherichia coli RNaseIII to digest dsRNA of zebrafish no tail (ntl), a mesoderm determinant in zebrafish and found that esi-ntl could lead to developmental defects, however, the effective dose was so close to the toxic dose that esi-ntl often led to non-specific developmental defects. Consequently, we utilized SP6 RNA polymerase to produce si-ntl, siRNA designed against ntl, by in vitro transcription. By injecting in vitro synthesized si-ntl into zebrafish zygotes, we obtained specific phenocopies of reported mutants of ntl. We achieved up to a 59%no tail phenotype when the injection concentration was as high as 4 mu g/mu L. Quantitative reverse transcription-polymerase chain reaction (RT-PCR) and whole-mount in situ hybridization analysis showed that si-ntl could largely and specifically reduce mRNA levels of the ntl gene. As a result, our data indicate that esiRNA is unable to cause specific developmental defects in zebrafish, while siRNA should be an alternative for downregulation of specific gene expression in zebrafish in cases where RNAi techniques are applied to zebrafish reverse genetics.

A key gene of the RNA interference pathway in the black tiger shrimp, Penaeus monodon: Identification and functional characterisation of Dicer-1

RNA interference (RNAi) is an evolutionarily conserved mechanism by which double-stranded RNA (dsRNA) initiates post-transcriptional silencing of homologous genes. Here we report the amplification and characterisation of a full length cDNA from black tiger shrimp (Penaeus monodon) that encodes the bidentate RNAase III Dicer, a key component of the RNAi pathway. The full length of the shrimp Dicer (Pm Dcr1) cDNA is 7629 bp in length, including a 51 untranslated region (UTR) of 130 bp, a 3' UTR of 77 bp, and an open reading frame of 7422 bp encoding a polypeptide of 2473 amino acids with an estimated molecular mass of 277.895 kDa and a predicted isoelectric point of 4.86. Analysis of the deduced amino acid sequence indicated that the mature peptide contains all the seven recognised functional domains and is most similar to the mosquito (Aedes aegypti) Dicer-1 sequence with a similarity of 34.6%. Quantitative RT-PCR analysis showed that Pm Dcr1 mRNA is most highly expressed in haemolymph and lymphoid organ tissues (P 0.05). However, there was no correlation between Pm Dcr1 mRNA levels in lymphoid organ and the viral genetic loads in shrimp naturally infected with gill-associated virus (GAV) and Mourilyan virus (P > 0.05). Treatment with synthetic dsRNA corresponding to Pm Dcr1 sequence resulted in knock-down of Pm Dcr1 mRNA expression in both uninfected shrimp and shrimp infected experimentally with GAV. Knock-down of Pm Dcr1 expression resulted in more rapid mortalities and higher viral loads. These data demonstrated that Dicer is involved in antiviral defence in shrimp. (c) 2007 Elsevier Ltd. All rights reserved.