Structural motifs in the RNA encoded by the early nodulation gene enod40 of soybean

The plant gene enod40 is highly conserved among legumes and also present in various non-legume species. It is presumed to play a central regulatory role in the Rhizobium–legume interaction, being expressed well before the initiation of cortical cell divisions resulting in nodule formation. Two small peptides encoded by enod40 mRNA as well as its secondary structure have been shown to be key elements in the signalling processes underlying nodule organogenesis. Here results concerning the secondary structure of mRNA of enod40 in soybean are presented. This study combined a theoretical approach, involving structure prediction and comparison, as well as structure probing. Our study indicates five conserved domains in enod40 mRNA among numerous leguminous species. Structure comparison suggests that some domains are also conserved in non-leguminous species and that an additional domain exists that was found only in leguminous species developing indeterminate nodules. Enzymatic and chemical probing data support the structure for three of the domains, and partially for the remaining two. The rest of the molecule appears to be less structured. Some of the domains include motifs, such as U-containing internal loops and bulges, which seem to be conserved. Therefore, they might be involved in the regulatory role of enod40 RNA.     

Translational and structural requirements of the early nodulin gene enod40, a short-open reading frame-containing RNA, for elicitation of a cell-specific growth response in the alfalfa root cortex

A diversity of mRNAs containing only short open reading frames (sORF-RNAs; encoding less than 30 amino acids) have been shown to be induced in growth and differentiation processes. The early nodulin gene enod40, coding for a 0.7-kb sORF-RNA, is expressed in the nodule primordium developing in the root cortex of leguminous plants after infection by symbiotic bacteria. Ballistic microtargeting of this gene into Medicago roots induced division of cortical cells. Translation of two sORFs (I and II, 13 and 27 amino acids, respectively) present in the conserved 5' and 3' regions of enod40 was required for this biological activity. These sORFs may be translated in roots via a reinitiation mechanism. In vitro translation products starting from the ATG of sORF I were detectable by mutating enod40 to yield peptides larger than 38 amino acids. Deletion of a Medicago truncatula enod40 region between the sORFs, spanning a predicted RNA structure, did not affect their translation but resulted in significantly decreased biological activity. Our data reveal a complex regulation of enod40 action, pointing to a role of sORF-encoded peptides and structured RNA signals in developmental processes involving sORF-RNAs.     

Reduction of cell size induced by enod40 in Arabidopsis thaliana

An extensive analysis of organ and cell size was performed in three different Arabidopsis lines transformed with the early nodulin gene enod40 under control of the CaMV35S promoter. All three transgenic lines presented a significant decrease in the mean size of both epidermal internode and leaf mesophyll cells. Flow cytometric and image analysis of enod40-transfected protoplasts prepared from wild-type Arabidopsis cell suspensions showed that transient expression of the gene resulted in reduced forward light scattering (a factor correlated with particle size) and cell size. The direct administration of ENOD40 peptide to fresh protoplasts also resulted in reduced forward scattering with respect to the control and to the administration of unrelated peptides. As far as is known this is the first report documenting a biological effect of enod40 at the cellular level in non-legume plants.     

Enod40, a short open reading frame-containing mRNA, induces cytoplasmic localization of a nuclear RNA binding protein in Medicago truncatula

In eukaryotes, diverse mRNAs containing only short open reading frames (sORF-mRNAs) are induced at specific stages of development. Their mechanisms of action may involve the RNA itself and/or sORF-encoded oligopeptides. Enod40 genes code for highly structured plant sORF-mRNAs involved in root nodule organogenesis. A novel RNA binding protein interacting with the enod40 RNA, MtRBP1 (for Medicago truncatula RNA Binding Protein 1), was identified using a yeast three-hybrid screening. Immunolocalization studies and use of a MtRBP1-DsRed2 fluorescent protein fusion showed that MtRBP1 localized to nuclear speckles in plant cells but was exported into the cytoplasm during nodule development in enod40-expressing cells. Direct involvement of the enod40 RNA in MtRBP1 relocalization into cytoplasmic granules was shown using a transient expression assay. Using a (green fluorescent protein)/MS2 bacteriophage system to tag the enod40 RNA, we detected in vivo colocalization of the enod40 RNA and MtRBP1 in these granules. This in vivo approach to monitor RNA-protein interactions allowed us to demonstrate that cytoplasmic relocalization of nuclear proteins is an RNA-mediated cellular function of a sORF-mRNA.     

Association of jacalin-related lectins with wheat responses to stresses revealed by transcriptional profiling

Jacalin-related lectins (JRLs) are carbohydrate-binding proteins widely present in plants and have one or more jacalin domains in common. However, JRLs’ structural types and functions are still poorly understood. In the present study, a total of 67 wheat (Triticum aestivum) JRL genes were identified through an exhausted search of EST database coupling with genome walking using published 454 sequence reads of Chinese Spring. A comparison of the translated wheat JRL proteins with those from other plants showed plant JRLs generally had low sequence similarity within and between species but exhibited conserved modular domain structures. More JRL genes encoded multiple jacalin domains in Arabidopsis thaliana, whereas more genes encoded chimeric JRLs in cereal plants. Dirigent domain-containing JRL genes were Poaceae-specific and accounted for nearly half of the identified wheat JRL genes. The dirigent domains were evolutionarily significantly correlated with the covalently linked jacalin domains. A phylogenetic analysis showed JRL proteins have experienced a substantial diversification after speciation. Moreover, new structural features conserved across the taxa were identified. Digital expression analysis and RT-PCR assays showed the expression of wheat JRL genes was largely tissue specific, typically low, and mostly inducible by biotic and abiotic stresses and stress hormones. These results suggest plant JRLs are critical for plant adaptation to stressful environments.     

Diversity,expression and mRNA targeting abilities of Argonaute-targeting miRNAs among selected vascular plants

Background

Micro (mi)RNAs are important regulators of plant development. Across plant lineages, Dicer-like 1 (DCL1) proteins process long ds-like structures to produce micro (mi) RNA duplexes in a stepwise manner. These miRNAs are incorporated into Argonaute (AGO) proteins and influence expression of RNAs that have sequence complementarity with miRNAs. Expression levels of AGOs are greatly regulated by plants in order to minimize unwarranted perturbations using miRNAs to target mRNAs coding for AGOs. AGOs may also have high promoter specificity-sometimes expression of AGO can be limited to just a few cells in a plant. Viral pathogens utilize various means to counter antiviral roles of AGOs including hijacking the host encoded miRNAs to target AGOs. Two host encoded miRNAs namely miR168 and miR403 that target AGOs have been described in the model plant Arabidopsis and such a mechanism is thought to be well conserved across plants because AGO sequences are well conserved.

Results

We show that the interaction between AGO mRNAs and miRNAs is species-specific due to the diversity in sequences of two miRNAs that target AGOs, sequence diversity among corresponding target regions in AGO mRNAs and variable expression levels of these miRNAs among vascular plants. We used miRNA sequences from 68 plant species representing 31 plant families for this analysis. Sequences of miR168 and miR403 are not conserved among plant lineages, but surprisingly they differ drastically in their sequence diversity and expression levels even among closely related plants. Variation in miR168 expression among plants correlates well with secondary structures/length of loop sequences of their precursors.

Conclusions

Our data indicates a complex AGO targeting interaction among plant lineages due to miRNA sequence diversity and sequences of miRNA targeting regions among AGO mRNAs, thus leading to the assumption that the perturbations by viruses that use host miRNAs to target antiviral AGOs can only be species-specific. We also show that rapid evolution and likely loss of expression of miR168 isoforms in tobacco is related to the insertion of MITE-like transposons between miRNA and miRNA* sequences, a possible mechanism showing how miRNAs are lost in few plant lineages even though other close relatives have abundantly expressing miRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1049) contains supplementary material, which is available to authorized users.     

Discovery of Replicating Circular RNAs by RNA-Seq and Computational Algorithms

Replicating circular RNAs are independent plant pathogens known as viroids, or act to modulate the pathogenesis of plant and animal viruses as their satellite RNAs. The rate of discovery of these subviral pathogens was low over the past 40 years because the classical approaches are technical demanding and time-consuming. We previously described an approach for homology-independent discovery of replicating circular RNAs by analysing the total small RNA populations from samples of diseased tissues with a computational program known as progressive filtering of overlapping small RNAs (PFOR). However, PFOR written in PERL language is extremely slow and is unable to discover those subviral pathogens that do not trigger in vivo accumulation of extensively overlapping small RNAs. Moreover, PFOR is yet to identify a new viroid capable of initiating independent infection. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was 3 to 8 times faster than PFOR. A new computational program was further developed and incorporated into PFOR2 to allow the identification of circular RNAs by deep sequencing of long RNAs instead of small RNAs. PFOR2 analysis of the small RNA libraries from grapevine and apple plants led to the discovery of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA (AHVd-like RNA), respectively. GLVd was proposed as a new species in the genus Apscaviroid, because it contained the typical structural elements found in this group of viroids and initiated independent infection in grapevine seedlings. AHVd-like RNA encoded a biologically active hammerhead ribozyme in both polarities, and was not specifically associated with any of the viruses found in apple plants. We propose that these computational algorithms have the potential to discover novel circular RNAs in plants, invertebrates and vertebrates regardless of whether they replicate and/or induce the in vivo accumulation of small RNAs.     

Next-Generation Sequencing of a 40 Mb Linkage Interval Reveals TSPAN12 Mutations in Patients with Familial Exudative Vitreoretinopathy

Familial exudative vitreoretinopathy (FEVR) is a genetically heterogeneous retinal disorder characterized by abnormal vascularisation of the peripheral retina, often accompanied by retinal detachment. To date, mutations in three genes (FZD4, LRP5, and NDP) have been shown to be causative for FEVR. In two large Dutch pedigrees segregating autosomal-dominant FEVR, genome-wide SNP analysis identified an FEVR locus of ∼40 Mb on chromosome 7. Microsatellite marker analysis suggested similar at risk haplotypes in patients of both families. To identify the causative gene, we applied next-generation sequencing in the proband of one of the families, by analyzing all exons and intron-exon boundaries of 338 genes, in addition to microRNAs, noncoding RNAs, and other highly conserved genomic regions in the 40 Mb linkage interval. After detailed bioinformatic analysis of the sequence data, prioritization of all detected sequence variants led to three candidates to be considered as the causative genetic defect in this family. One of these variants was an alanine-to-proline substitution in the transmembrane 4 superfamily member 12 protein, encoded by TSPAN12. This protein has very recently been implicated in regulating the development of retinal vasculature, together with the proteins encoded by FZD4, LRP5, and NDP. Sequence analysis of TSPAN12 revealed two mutations segregating in five of 11 FEVR families, indicating that mutations in TSPAN12 are a relatively frequent cause of FEVR. Furthermore, we demonstrate the power of targeted next-generation sequencing technology to identify disease genes in linkage intervals.     

A Nucleolar PUF RNA-binding Protein with Specificity for a Unique RNA Sequence

PUF proteins are a conserved group of sequence specific RNA-binding proteins that bind to RNA in a modular fashion. The RNA-binding domain of PUF proteins typically consists of eight clustered Puf repeats. Plant genomes code for large families of PUF proteins that show significant variability in their predicted Puf repeat number, organization, and amino acid sequence. Here we sought to determine whether the observed variability in the RNA-binding domains of four plant PUFs results in a preference for nonclassical PUF RNA target sequences. We report the identification of a novel RNA binding sequence for a nucleolar Arabidopsis PUF protein that contains an atypical RNA-binding domain. The Arabidopsis PUM23 (APUM23) binding sequence was 10 nucleotides in length, contained a centrally located UUGA core element, and had a preferred cytosine at nucleotide position 8. These RNA sequence characteristics differ from those of other PUF proteins, because all natural PUFs studied to date bind to RNAs that contain a conserved UGU sequence at their 5′ end and lack specificity for cytosine. Gel mobility shift assays validated the identity of the APUM23 binding sequence and supported the location of 3 of the 10 predicted Puf repeats in APUM23, including the cytosine-binding repeat. The preferred 10-nucleotide sequence bound by APUM23 is present within the 18S rRNA sequence, supporting the known role of APUM23 in 18S rRNA maturation. This work also reveals that APUM23, an ortholog of yeast Nop9, could provide an advanced structural backbone for Puf repeat engineering and target-specific regulation of cellular RNAs.     

The conserved 3′ Angio‐domain defines a superfamily of short interspersed nuclear elements (SINEs) in higher plants

Repetitive sequences are ubiquitous components of eukaryotic genomes affecting genome size and evolution as well as gene regulation. Among them, short interspersed nuclear elements (SINEs) are non‐coding retrotransposons usually shorter than 1000 bp. They contain only few short conserved structural motifs, in particular an internal promoter derived from cellular RNAs and a mostly AT‐rich 3′ tail, whereas the remaining regions are highly variable. SINEs emerge and vanish during evolution, and often diversify into numerous families and subfamilies that are usually specific for only a limited number of species. In contrast, at the 3′ end of multiple plant SINEs we detected the highly conserved ‘Angio‐domain’. This 37 bp segment defines the Angio‐SINE superfamily, which encompasses 24 plant SINE families widely distributed across 13 orders within the plant kingdom. We retrieved 28 433 full‐length Angio‐SINE copies from genome assemblies of 46 plant species, frequently located in genes. Compensatory mutations in and adjacent to the Angio‐domain imply selective restraints maintaining its RNA structure. Angio‐SINE families share segmental sequence similarities, indicating a modular evolution with strong Angio‐domain preservation. We suggest that the conserved domain contributes to the evolutionary success of Angio‐SINEs through either structural interactions between SINE RNA and proteins increasing their transpositional efficiency, or by enhancing their accumulation in genes.     

Deep Sequencing and Microarray Hybridization Identify Conserved and Species-Specific MicroRNAs during Somatic Embryogenesis in Hybrid Yellow Poplar

Background

To date, several studies have indicated a major role for microRNAs (miRNAs) in regulating plant development, but miRNA-mediated regulation of the developing somatic embryo is poorly understood, especially during early stages of somatic embryogenesis in hardwood plants. In this study, Solexa sequencing and miRNA microfluidic chips were used to discover conserved and species-specific miRNAs during somatic embryogenesis of hybrid yellow poplar (Liriodendron tulipifera×L. chinense).

Methodology/Principal Findings

A total of 17,214,153 reads representing 7,421,623 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from all stages of somatic embryos. Through a combination of deep sequencing and bioinformatic analyses, we discovered 83 sequences with perfect matches to known miRNAs from 33 conserved miRNA families and 273 species-specific candidate miRNAs. MicroRNA microarray results demonstrated that many conserved and species-specific miRNAs were expressed in hybrid yellow poplar embryos. In addition, the microarray also detected another 149 potential miRNAs, belonging to 29 conserved families, which were not discovered by deep sequencing analysis. The biological processes and molecular functions of the targets of these miRNAs were predicted by carrying out BLAST search against Arabidopsis thaliana GenBank sequences and then analyzing the results with Gene Ontology.

Conclusions

Solexa sequencing and microarray hybridization were used to discover 232 candidate conserved miRNAs from 61 miRNA families and 273 candidate species-specific miRNAs in hybrid yellow poplar. In these predicted miRNAs, 64 conserved miRNAs and 177 species-specific miRNAs were detected by both sequencing and microarray hybridization. Our results suggest that miRNAs have wide-ranging characteristics and important roles during all stages of somatic embryogenesis in this economically important species.