24-nt sRNAs Research Articles

Post-transcriptional adaptation of the aquatic plant Spirodela polyrhiza under stress and hormonal stimuli.

The Lemnaceae family comprises aquatic plants of angiosperms gaining attention due to their utility in wastewater treatment, and rapid production of biomass that can be used as feed, fuel, or food. Moreover, it can serve as a model species for neotenous growth and environmental adaptation. The latter properties are subject to post-transcriptional regulation of gene expression, meriting investigation of how miRNAs in Spirodela polyrhiza, the most basal and most thoroughly sequenced member of the family, are expressed under different growth conditions. To further scientific understanding of its capacity to adapt to environmental cues, we measured miRNA expression and processing of their target sequences under different temperatures, and in the presence of abscisic acid, copper, kinetin, nitrate, and sucrose. Using two small RNA sequencing experiments and one degradome sequencing experiment, we provide evidence for 108 miRNAs. Sequencing cleaved mRNAs validated 42 conserved miRNAs with 83 targets and 24 novel miRNAs regulating 66 targets and created a list of 575 predicted and verified targets. These analyses revealed condition-induced changes in miRNA expression and cleavage activity, and resulted in the addition of stringently reviewed miRNAs to miRBase. This combination of small RNA and degradome sequencing provided not only high confidence predictions of conserved and novel miRNAs and targets, but also a view of the post-transcriptional regulation of adaptations. A unique aspect is the role of miR156 and miR172 expression and activity in its clonal propagation and neoteny. Additionally, low levels of 24nt sRNAs were observed, despite the lack of recent retrotransposition.

Assessing the Response of Small RNA Populations to Allopolyploidy Using Resynthesized Brassica napus Allotetraploids.

Allopolyploidy, combining interspecific hybridization with whole genome duplication, has had significant impact on plant evolution. Its evolutionary success is related to the rapid and profound genome reorganizations that allow neoallopolyploids to form and adapt. Nevertheless, how neoallopolyploid genomes adapt to regulate their expression remains poorly understood. The hypothesis of a major role for small noncoding RNAs (sRNAs) in mediating the transcriptional response of neoallopolyploid genomes has progressively emerged. Generally, 21-nt sRNAs mediate posttranscriptional gene silencing by mRNA cleavage, whereas 24-nt sRNAs repress transcription (transcriptional gene silencing) through epigenetic modifications. Here, we characterize the global response of sRNAs to allopolyploidy in Brassica, using three independently resynthesized Brassica napus allotetraploids originating from crosses between diploid Brassica oleracea and Brassica rapa accessions, surveyed at two different generations in comparison with their diploid progenitors. Our results suggest an immediate but transient response of specific sRNA populations to allopolyploidy. These sRNA populations mainly target noncoding components of the genome but also target the transcriptional regulation of genes involved in response to stresses and in metabolism; this suggests a broad role in adapting to allopolyploidy. We finally identify the early accumulation of both 21- and 24-nt sRNAs involved in regulating the same targets, supporting a posttranscriptional gene silencing to transcriptional gene silencing shift at the first stages of the neoallopolyploid formation. We propose that reorganization of sRNA production is an early response to allopolyploidy in order to control the transcriptional reactivation of various noncoding elements and stress-related genes, thus ensuring genome stability during the first steps of neoallopolyploid formation.

Citrus tristeza virus: Host RNA Silencing and Virus Counteraction.

To dissect the host RNA silencing response incited by citrus tristeza virus (CTV, genus Closterovirus), a (+) ssRNA of ~19300nt, and the counter reaction deployed by the virus via its three RNA silencing suppressors (RSS), the small RNAs (sRNAs) of three virus-host combinations were deep sequenced. The subsequent analysis indicated that CTV sRNAs (1) constitute more than half of the total sRNAs in the susceptible Mexican lime and sweet orange, while only 3.5% in the restrictive sour orange; (2) are mostly of 21-22 nt, with those of (+) sense predominating slightly; and (3) derive from all the CTV genome, as evidenced by its entire recomposition from viral sRNA contigs but adopt an asymmetric pattern with a hotspot mapping at the 3'-terminal ~2500nt. The citrus homologues of Arabidopsis Dicer-like (DCL) 4 and 2 most likely generate the 21 and 22nt CTV sRNAs, respectively, by dicing the gRNA and the 3' co-terminal sgRNAs and, particularly, their double-stranded forms accumulating in infected cells. The plant sRNA profile, very similar and dominated by the 24 nt sRNAs in the three mock-inoculated controls, displayed a major reduction of the 24 nt sRNAs in Mexican lime and sweet orange, but not in sour orange. CTV infection also influences the levels of certain microRNAs.The high accumulation of CTV sRNAs in two of the citrus hosts examined suggests that it is not their synthesis, but their function, the target of the RSS encoded by CTV: p25 (intercellular), p23 (intracellular) and p20 (both). The two latter might block the loading of CTV sRNAs into the RNA silencing complex or interfere with it through alternative mechanisms. Of the three CTV RSS, p23 is the one that has been more thoroughly studied. It is a multifunctional RNA-binding protein with a putative Zn finger domain and basic motifs that (1) has no homologues in other closteroviruses, (2) accumulates in the nucleolus and plasmodesmata, (3) regulates the asymmetric balance of CTV (+) and (-) RNA strands, and (4) induces CTV syndromes and stimulates systemic infection in certain citrus species when expressed as a transgene ectopically or in phloem-associated cells.

Small RNA differential expression and regulation in Tuxpeño maize embryogenic callus induction and establishment

Somatic embryogenesis represents an alternative developmental process used to achieve genetic transformation and to approach key questions in maize development. It is known that embryogenic callus induction and plant regeneration are accompanied by microRNA expression changes. However, small RNA (sRNA) populations have not been explored during the proliferative callus subculture establishment and their impact on maintaining the dedifferentiated status and embryogenic potential is far from being completely understood. Here we globally tested the sRNA populations in explants (immature embryos), induced and established maize embryogenic callus from the Mexican cultivar VS-535, Tuxpeño landrace. We detected readjustments in 24 nt and 21–22 nt sRNAs during the embryogenic callus (EC) establishment and maintenance. A follow up on specific microRNAs (miRNAs) indicated that miRNAs related to stress response substantially increase upon the callus proliferation establishment, correlating with a reduction in some of their target levels. On the other hand, while 24 nt-long heterochromatic small interfering RNAs (hc-siRNAs) derived from transposable retroelements transiently decreased in abundance during the EC establishment, a population of 22 nt-hc-siRNAs increased. This was accompanied by reduction in transposon expression in the established callus subcultures. We conclude that stress- and development-related miRNAs are highly expressed upon maize EC callus induction and during maintenance of the subcultures, while miRNAs involved in hormone response only transiently increase during induction. In addition, the establishment of a proliferative status in embryogenic callus is accompanied by important readjustments in hc-siRNAs mapping to long tandem repeat (LTR) retrotransposons, and their expression regulation.

A possible role of CTV.20 gene methylation in response to Citrus tristeza virus infection

Citrus tristeza virus (CTV) is the pathogen causing tristeza diseases in several Rutaceae species and leading to significant economic damage to citrus worldwide. The Ctv locus provides broad spectrum resistance to CTV in Poncirus trifoliata L. Raf. This locus is present also in tolerant and susceptible species, so an epigenetic mechanism of Ctv expression regulation was proposed. Indeed, a difference in plant 24-nt sRNAs distribution corresponding to CTV.20 gene was previously observed in susceptible species following CTV infection. This gene, encoding for a plant virus movement-like protein, was investigated as a candidate gene for CTV susceptibility. Here, we show the presence of differences in methylation status of a specific region of CTV.20 in two susceptible species, sour and sweet orange, following CTV infection. On the contrary, no significant differences were observed in the tolerant Citrange carrizo following the infection. Moreover, a hypermethylation of the whole CTV.20 gene was observed in Citrange carrizo both healthy and infected, and in healthy sour and sweet orange. This preliminary study allows hypothesizing a possible role of methylation in regulation of CTV.20 expression involved in the CTV susceptibility.

The 25-26 nt Small RNAs in Phytophthora parasitica Are Associated with Efficient Silencing of Homologous Endogenous Genes.

Small RNAs (sRNAs) are important non-coding RNA regulators, playing key roles in developmental regulation, transposon suppression, environmental response, host–pathogen interaction and other diverse biological processes. However, their roles in oomycetes are poorly understood. Here, we performed sRNA sequencing and RNA sequencing of Phytophthora parasitica at stages of vegetative growth and infection of Arabidopsis roots to examine diversity and function of sRNAs in P. parasitica, a model hemibiotrophic oomycete plant pathogen. Our results indicate that there are two distinct types of sRNA-generating loci in P. parasitica genome, giving rise to clusters of 25–26 nt and 21 nt sRNAs, respectively, with no significant strand-biases. The 25–26 nt sRNA loci lie predominantly in gene-sparse and repeat-rich regions, and overlap with over 7000 endogenous gene loci. These overlapped genes are typically P. parasitica species-specific, with no homologies to the sister species P. infestans. They include approximately 40% RXLR effector genes, 50% CRN effector genes and some elicitor genes. The transcripts of most of these genes could not be detected at both the vegetative mycelium and infection stages as revealed by RNA sequencing, indicating that the 25–26 nt sRNAs are associated with efficient silencing of these genes. The 21 nt sRNA loci typically overlap with the exon regions of highly expressed genes, suggesting that the biogenesis of the 21 nt sRNAs may be dependent on the level of gene transcription and that these sRNAs do not mediate efficient silencing of homologous genes. Analyses of the published P. infestans sRNA and mRNA sequencing data consistently show that the 25–26 nt sRNAs, but not the 21 nt sRNAs, may mediate efficient gene silencing in Phytophthora.

Landscape of Fluid Sets of Hairpin-Derived 21-/24-nt-Long Small RNAs at Seed Set Uncovers Special Epigenetic Features in Picea glauca.

Conifers’ exceptionally large genome (20–30 Gb) is scattered with 60% retrotransposon (RT) components and we have little knowledge on their origin and evolutionary implications. RTs may impede the expression of flanking genes and provide sources of the formation of novel small RNA (sRNAs) populations to constrain events of transposon (TE) proliferation/transposition. Here we show a declining expression of 24-nt-long sRNAs and low expression levels of their key processing gene, pgRTL2 (RNASE THREE LIKE 2) at seed set in Picea glauca. The sRNAs in 24-nt size class are significantly less enriched in type and read number than 21-nt sRNAs and have not been documented in other species. The architecture of MIR loci generating highly expressed 24-/21-nt sRNAs is featured by long terminal repeat—retrotransposons (LTR-RTs) in families of Ty3/Gypsy and Ty1/Copia elements. This implies that the production of sRNAs may be predominantly originated from TE fragments on chromosomes. Furthermore, a large proportion of highly expressed 24-nt sRNAs does not have predictable targets against unique genes in Picea, suggestive of their potential pathway in DNA methylation modifications on, for instance, TEs. Additionally, the classification of computationally predicted sRNAs suggests that 24-nt sRNA targets may bear particular functions in metabolic processes while 21-nt sRNAs target genes involved in many different biological processes. This study, therefore, directs our attention to a possible extrapolation that lacking of 24-nt sRNAs at the late conifer seed developmental phase may result in less constraints in TE activities, thus contributing to the massive expansion of genome size.

Generation of different sizes and classes of small RNAs in barley is locus, chromosome and/or cultivar-dependent.

BackgroundVarious small RNA (sRNA) sizes and varieties have been identified, but their relationship as well as relationship with their origins and allocations have not been well understood or investigated.ResultsBy comparing sRNAs generated from two barley cultivars, Golden Promise (GP) and Pallas, we identified that the generation of different sizes and types of sRNAs in barley was locus-, chromosome- and/or cultivar-dependent. 20-nt sRNAs mainly comprising miRNAs and chloroplast-derived sRNAs were significantly over-expressed in Pallas vs. GP on chromosomes 3H and 6H. MiRNAs-enriched 21-nt sRNAs were significantly over-expressed in Pallas vs. GP only on chromosome 4H. On chromosome 5H this size of sRNAs was significantly under-expressed in Pallas, so were 22-nt sRNAs mainly comprising miRNAs and repeat-derived sRNAs. 24-nt sRNAs mostly derived from repeats were evenly distributed in all chromosomes and expressed similarly between GP and Pallas. Unlike other sizes of sRNAs, 24-nt sRNAs were little conserved in other plant species. Abundant sRNAs were mostly generated from 3’ terminal regions of chromosome 1H and 5’ terminal regions of chromosome 5H. Over-expressed miRNAs in GP vs. Pallas primarily function in stress responses and iron-binding.ConclusionsOur study indicates that 23−24-nt sRNAs may be linked to repressive chromatin modifications and function in genome stability while 20−21-nt sRNAs may be important for the cultivar specificity. This study provides a novel insight into the mechanism of sRNA expression and function in barley.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3023-5) contains supplementary material, which is available to authorized users.

The Tomato Spotted Wilt Virus Genome Is Processed Differentially in its Plant Host Arachis hypogaea and its Thrips Vector Frankliniella fusca

Thrips-transmitted tospoviruses are economically important viruses affecting a wide range of field and horticultural crops worldwide. Tomato spotted wilt virus (TSWV) is the type member of the Tospovirus genus with a broad host range of more than 900 plant species. Interactions between these viruses and their plant hosts and insect vectors via RNAi pathways are likely a key determinant of pathogenicity. The current investigation, for the first time, compares biogenesis of small RNAs between the plant host and insect vector in the presence or absence of TSWV. Unique viral small interfering RNA (vsiRNA) profiles are evident for Arachis hypogaea (peanut) and Frankliniella fusca (thrips vector) following infection with TSWV. Differences between vsiRNA profiles for these plant and insect species, such as the relative abundance of 21 and 22 nt vsiRNAs and locations of alignment hotspots, reflect the diverse siRNA biosynthesis pathways of their respective kingdoms. The presence of unique vsiRNAs in F. fusca samples indicates that vsiRNA generation takes place within the thrips, and not solely through uptake via feeding on vsiRNAs produced in infected A. hypogaea. The study also shows key vsiRNA profile differences for TSWV among plant families, which are evident in the case of A. hypogaea, a legume, and members of Solanaceae (S. lycopersicum and Nicotiana benthamiana). Distinctively, overall small RNA (sRNA) biogenesis in A. hypogaea is markedly affected with an absence of the 24 nt sRNAs in TSWV-infected plants, possibly leading to wide-spread molecular and phenotypic perturbations specific to this species. These findings add significant information on the host–virus–vector interaction in terms of RNAi pathways and may lead to better crop and vector specific control strategies.

Small RNA and methylation responses in susceptible and tolerant landraces of cassava infected with South African cassava mosaic virus

Endogenous small RNAs (sRNAs) associated with gene regulatory mechanisms respond to virus infection, and virus-derived small RNAs (vsRNAs) have been implicated in recovery or symptom remission in some geminivirus-host interactions. Transcriptional gene silencing (TGS) (24 nt vsRNAs) and post transcriptional gene silencing (PTGS) (21–23 nt vsRNAs) have been associated with geminivirus intergenic (IR) and coding regions, respectively. In this Illumina deep sequencing study, we compared for the first time, the small RNA response to South African cassava mosaic virus (SACMV) of cassava landrace TME3 which shows a recovery and tolerant phenotype, and T200, a highly susceptible landrace. Interestingly, different patterns in the percentage of SACMV-induced normalized total endogenous sRNA reads were observed between T200 and TME3. Notably in virus-infected T200 there was an increase in 21 nt sRNAs during the early pre-symptomatic response (12dpi) compared to mock, while in TME3, the 22 nt sRNA size class was predominant at 32days post infection with SACMV. While vsRNAs of 21–24 nt size classes mapped to the entire SACMV DNA-A and DNA-B genome components in T200 and TME3, vsRNA population counts were lower at 32 (symptomatic stage) and 67 dpi (recovery stage) in tolerant TME3 compared with T200 (non-recovery). It is suggested that the high accumulation of primary vsRNAs, which correlated with high virus titers and severe symptoms in susceptible T200, may be due to failure to target SACMV-derived mRNA. Likewise, in contrast, in TME3 low vsRNA counts may represent efficient PTGS of viral mRNA, leading to a depletion/sequestration of vsRNA populations, supporting a role for PTGS in tolerance/recovery in TME3. Notably, in TME3 at recovery (67 dpi) the percentage (expressed as a percentage of total vsRNA counts) of redundant and non-redundant (unique) 24 nt vsRNAs increased dramatically. Since methylation of the SACMV genome was not detected by bisulfite sequencing, and vsRNA counts targeting the intergenic region (where the promoters reside) were very low in both the tolerant or susceptible landraces, we could not provide conclusive evidence that 24 nt vsRNA-mediated RNA directed genome methylation plays a central role in disease phenotype in these landraces, notwithstanding recognition for a possible role in histone modification in TME3. This work represents an important step toward understanding variable roles of sRNAs in different cassava genotype-geminivirus interactions.

Genome-wide identification of microRNAs in pomegranate (Punica granatum L.) by high-throughput sequencing.

BackgroundMicroRNAs (miRNAs), a class of small non-coding endogenous RNAs that regulate gene expression post-transcriptionally, play multiple key roles in plant growth and development and in biotic and abiotic stress response. Knowledge and roles of miRNAs in pomegranate fruit development have not been explored.ResultsPomegranate, which accumulates a large amount of anthocyanins in skin and arils, is valuable to human health, mainly because of its antioxidant properties. In this study, we developed a small RNA library from pooled RNA samples from young seedlings to mature fruits and identified both conserved and pomegranate-specific miRNA from 29,948,480 high-quality reads. For the pool of 15- to 30-nt small RNAs, ~50 % were 24 nt. The miR157 family was the most abundant, followed by miR156, miR166, and miR168, with variants within each family. The base bias at the first position from the 5’ end had a strong preference for U for most 18- to 26-nt sRNAs but a preference for A for 18-nt sRNAs. In addition, for all 24-nt sRNAs, the nucleotide U was preferred (97 %) in the first position. Stem-loop RT-qPCR was used to validate the expression of the predominant miRNAs and novel miRNAs in leaves, male and female flowers, and multiple fruit developmental stages; miR156, miR156a, miR159a, miR159b, and miR319b were upregulated during the later stages of fruit development. Higher expression of miR156 in later fruit developmental may positively regulate anthocyanin biosynthesis by reducing SPL transcription factor. Novel miRNAs showed variation in expression among different tissues. These novel miRNAs targeted different transcription factors and hormone related regulators. Gene ontology and KEGG pathway analyses revealed predominant metabolic processes and catalytic activities, important for fruit development. In addition, KEGG pathway analyses revealed the involvement of miRNAs in ascorbate and linolenic acid, starch and sucrose metabolism; RNA transport; plant hormone signaling pathways; and circadian clock.ConclusionOur first and preliminary report of miRNAs will provide information on the synthesis of biochemical compounds of pomegranate for future research. The functions of the targets of the novel miRNAs need further investigation.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-016-0807-3) contains supplementary material, which is available to authorized users.

Identification of miRNAs and their targets through high-throughput sequencing and degradome analysis in male and female Asparagus officinalis.

BackgroundMicroRNAs (miRNAs), a class of non-coding small RNAs (sRNAs), regulate various biological processes. Although miRNAs have been identified and characterized in several plant species, miRNAs in Asparagus officinalis have not been reported. As a dioecious plant with homomorphic sex chromosomes, asparagus is regarded as an important model system for studying mechanisms of plant sex determination.ResultsTwo independent sRNA libraries from male and female asparagus plants were sequenced with Illumina sequencing, thereby generating 4.13 and 5.88 million final clean reads, respectively. Both libraries predominantly contained 24-nt sRNAs, followed by 21-nt sRNAs. Further analysis identified 154 conserved miRNAs, which belong to 26 families, and 39 novel miRNA candidates seemed to be specific to asparagus. Comparative profiling revealed that 63 miRNAs exhibited significant differential expression between male and female plants, which was confirmed by real-time quantitative PCR analysis. Among them, 37 miRNAs were significantly up-regulated in the female library, whereas the others were preferentially expressed in the male library. Furthermore, 40 target mRNAs representing 44 conserved and seven novel miRNAs were identified in asparagus through high-throughput degradome sequencing. Functional annotation showed that these target mRNAs were involved in a wide range of developmental and metabolic processes.ConclusionsWe identified a large set of conserved and specific miRNAs and compared their expression levels between male and female asparagus plants. Several asparagus miRNAs, which belong to the miR159, miR167, and miR172 families involved in reproductive organ development, were differentially expressed between male and female plants, as well as during flower development. Consistently, several predicted targets of asparagus miRNAs were associated with floral organ development. These findings suggest the potential roles of miRNAs in sex determination and reproductive developmental processes in asparagus.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-016-0770-z) contains supplementary material, which is available to authorized users.

Systemic silencing: Mobile sRNA stabilizes genomes.

Transcriptional gene silencing is a pivotal mechanism for regulating gene expression and genome stability. In Arabidopsis, combined analyses of small RNAs (sRNAs) and DNA methylation reveals that mobile 24-nt sRNAs are involved in reinforcing genome-wide silencing of transposons through DNA methylation.

Mobile small RNAs regulate genome-wide DNA methylation.

RNA silencing at the transcriptional and posttranscriptional levels regulates endogenous gene expression, controls invading transposable elements (TEs), and protects the cell against viruses. Key components of the mechanism are small RNAs (sRNAs) of 21-24 nt that guide the silencing machinery to their nucleic acid targets in a nucleotide sequence-specific manner. Transcriptional gene silencing is associated with 24-nt sRNAs and RNA-directed DNA methylation (RdDM) at cytosine residues in three DNA sequence contexts (CG, CHG, and CHH). We previously demonstrated that 24-nt sRNAs are mobile from shoot to root in Arabidopsis thaliana and confirmed that they mediate DNA methylation at three sites in recipient cells. In this study, we extend this finding by demonstrating that RdDM of thousands of loci in root tissues is dependent upon mobile sRNAs from the shoot and that mobile sRNA-dependent DNA methylation occurs predominantly in non-CG contexts. Mobile sRNA-dependent non-CG methylation is largely dependent on the DOMAINS REARRANGED METHYLTRANSFERASES 1/2 (DRM1/DRM2) RdDM pathway but is independent of the CHROMOMETHYLASE (CMT)2/3 DNA methyltransferases. Specific superfamilies of TEs, including those typically found in gene-rich euchromatic regions, lose DNA methylation in a mutant lacking 22- to 24-nt sRNAs (dicer-like 2, 3, 4 triple mutant). Transcriptome analyses identified a small number of genes whose expression in roots is associated with mobile sRNAs and connected to DNA methylation directly or indirectly. Finally, we demonstrate that sRNAs from shoots of one accession move across a graft union and target DNA methylation de novo at normally unmethylated sites in the genomes of root cells from a different accession.

Differentially expressed small RNAs in Arabidopsis galls formed by Meloidogyne javanica: a functional role for miR390 and its TAS3-derived tasiRNAs.

Root-knot nematodes (RKNs) induce inside the vascular cylinder the giant cells (GCs) embedded in the galls. The distinctive gene repression in early-developing GCs could be facilitated by small RNAs (sRNA) such as miRNAs, and/or epigenetic mechanisms mediated by 24nt-sRNAs, rasiRNAs and 21-22nt-sRNAs. Therefore, the sRNA-population together with the role of the miR390/TAS3/ARFs module were studied during early gall/GC formation. Three sRNA libraries from 3-d-post-inoculation (dpi) galls induced by Meloidogyne javanica in Arabidopsis and three from uninfected root segments were sequenced following Illumina-Solexa technology. pMIR390a::GUS and pTAS3::GUS lines were assayed for nematode-dependent promoter activation. A sensor line indicative of TAS3-derived tasiRNAs binding to the ARF3 sequence (pARF3:ARF3-GUS) together with a tasiRNA-resistant ARF3 line (pARF3:ARF3m-GUS) were used for functional analysis. The sRNA population showed significant differences between galls and controls, with high validation rate and correspondence with their target expression: 21-nt sRNAs corresponding mainly to miRNAs were downregulated, whilst 24-nt-sRNAs from the rasiRNA family were mostly upregulated in galls. The promoters of MIR390a and TAS3, active in galls, and the pARF3:ARF3-GUS line, indicated a role of TAS3-derived-tasiRNAs in galls. The regulatory module miR390/TAS3 is necessary for proper gall formation possibly through auxin-responsive factors, and the abundance of 24-nt sRNAs (mostly rasiRNAs) constitutes a gall hallmark.

Identification and expression profiles of sRNAs and their biogenesis and action-related genes in male and female cones of Pinus tabuliformis.

BackgroundSmall RNA (sRNA) play pivotal roles in reproductive development, and their biogenesis and action mechanisms are well characterised in angiosperm plants; however, corresponding studies in conifers are very limited. To improve our understanding of the roles of sRNA pathways in the reproductive development of conifers, the genes associated with sRNA biogenesis and action pathways were identified and analysed, and sRNA sequencing and parallel analysis of RNA ends (PARE) were performed in male and female cones of the Chinese pine (Pinus tabuliformis).ResultsBased on high-quality reference transcriptomic sequences, 21 high-confidence homologues involved in sRNA biogenesis and action in P. tabuliformis were identified, including two different DCL3 genes and one AGO4 gene. More than 75 % of genes involved in sRNA biogenesis and action have higher expression levels in female than in male cones. Twenty-six microRNA (miRNA) families and 74 targets, including 46 24-nt sRNAs with a 5’ A, which are specifically expressed in male cones or female cones and probably bind to AGO4, were identified.ConclusionsThe sRNA pathways have higher activity in female than in male cones, and the miRNA pathways are the main sRNA pathways in P. tabuliformis. The low level of 24-nt short-interfering RNAs in conifers is not caused by the absence of biogenesis-related genes or AGO-binding proteins, but most likely caused by the low accumulation of these key components. The identification of sRNAs and their targets, as well as genes associated with sRNA biogenesis and action, will provide a good starting point for investigations into the roles of sRNA pathways in cone development in conifers.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1885-6) contains supplementary material, which is available to authorized users.

Nicotiana small RNA sequences support a host genome origin of cucumber mosaic virus satellite RNA.

Satellite RNAs (satRNAs) are small noncoding subviral RNA pathogens in plants that depend on helper viruses for replication and spread. Despite many decades of research, the origin of satRNAs remains unknown. In this study we show that a β-glucuronidase (GUS) transgene fused with a Cucumber mosaic virus (CMV) Y satellite RNA (Y-Sat) sequence (35S-GUS:Sat) was transcriptionally repressed in N. tabacum in comparison to a 35S-GUS transgene that did not contain the Y-Sat sequence. This repression was not due to DNA methylation at the 35S promoter, but was associated with specific DNA methylation at the Y-Sat sequence. Both northern blot hybridization and small RNA deep sequencing detected 24-nt siRNAs in wild-type Nicotiana plants with sequence homology to Y-Sat, suggesting that the N. tabacum genome contains Y-Sat-like sequences that give rise to 24-nt sRNAs capable of guiding RNA-directed DNA methylation (RdDM) to the Y-Sat sequence in the 35S-GUS:Sat transgene. Consistent with this, Southern blot hybridization detected multiple DNA bands in Nicotiana plants that had sequence homology to Y-Sat, suggesting that Y-Sat-like sequences exist in the Nicotiana genome as repetitive DNA, a DNA feature associated with 24-nt sRNAs. Our results point to a host genome origin for CMV satRNAs, and suggest novel approach of using small RNA sequences for finding the origin of other satRNAs.

Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C

Persistent infections with hepatitis B virus (HBV) or hepatitis C virus (HCV) account for the majority of cases of hepatic cirrhosis and hepatocellular carcinoma (HCC) worldwide. Small, non-coding RNAs play important roles in virus-host interactions. We used high throughput sequencing to conduct an unbiased profiling of small (14-40 nts) RNAs in liver from Japanese subjects with advanced hepatitis B or C and hepatocellular carcinoma (HCC). Small RNAs derived from tRNAs, specifically 30–35 nucleotide-long 5′ tRNA-halves (5′ tRHs), were abundant in non-malignant liver and significantly increased in humans and chimpanzees with chronic viral hepatitis. 5′ tRH abundance exceeded microRNA abundance in most infected non-cancerous tissues. In contrast, in matched cancer tissue, 5′ tRH abundance was reduced, and relative abundance of individual 5′ tRHs was altered. In hepatitis B-associated HCC, 5′ tRH abundance correlated with expression of the tRNA-cleaving ribonuclease, angiogenin. These results demonstrate that tRHs are the most abundant small RNAs in chronically infected liver and that their abundance is altered in liver cancer.

Epigenetic transitions leading to heritable, RNA-mediated de novo silencing in Arabidopsis thaliana

In plants, RNA-directed DNA methylation (RdDM), a mechanism where epigenetic modifiers are guided to target loci by small RNAs, plays a major role in silencing of transposable elements (TEs) to maintain genome integrity. So far, two RdDM pathways have been identified: RNA Polymerase IV (PolIV)-RdDM and RNA-dependent RNA Polymerase 6 (RDR6)-RdDM. PolIV-RdDM involves a self-reinforcing feedback mechanism that maintains TE silencing, but cannot explain how epigenetic silencing is first initiated. A function of RDR6-RdDM is to reestablish epigenetic silencing of active TEs, but it is unknown if this pathway can induce DNA methylation at naïve, non-TE loci. To investigate de novo establishment of RdDM, we have used virus-induced gene silencing (VIGS) of an active flowering Wageningen epiallele. Using genetic mutants we show that unlike PolIV-RdDM, but like RDR6-RdDM, establishment of VIGS-mediated RdDM requires PolV and DRM2 but not Dicer like-3 and other PolIV pathway components. DNA methylation in VIGS is likely initiated by a process guided by virus-derived small (s) RNAs that are 21/22-nt in length and reinforced or maintained by 24-nt sRNAs. We demonstrate that VIGS-RdDM as a tool for gene silencing can be enhanced by use of mutant plants with increased production of 24-nt sRNAs to reinforce the level of RdDM.

Comparison of micro RNA (miRNA) Profiles and Some miRNA Target Gene Expression levels in Roots of Non-infected and Huanglongbing-infected Tangerine (Citrus reticulata Blanco cv.‘Sanhu’ ) Trees

Solexa sequencing was used to reveal the changes in small RNAome profile in roots of mock-inoculated (CK) and Huanglongbing bacteria-inoculated (HLB) ‘Sanhu’ tangerine (Citrus reticulata Blanco) trees. Results showed that the number of reads of both unique and total sRNAs decreased apparently in roots following infection with HLB. Distribution in length of sRNAs changed also remarkably, showing an increase in 22 nt and 21 nt small RNAs (sRNAs) and a decrease in 24 nt sRNAs in HLB-infected samples. A total of 42 known micro RNAs (miRNAs) belonging to 27 highly conserved miRNA families were identified. Comparisons showed that 33 known miRNAs exhibited a significant expression difference between CK and HLB-inoculated samples. In addition, 34 novel miRNAs, among which 24 were differentially expressed, were also identified, and their expression levels were analyzed by qRT-PCR. Three hundred and eighty five potential target genes were predicted for most of the 57 differentially expressed miRNAs. GO and KEGG annotation analysis revealed that most miRNA-target genes were those implicated in developmental process, response to stress and stimulus, transcription and protein metabolism. The characterization of the miRNAomes between the healthy and HLB infected Sanhu tangerine roots provided new insight into the involvement of miRNAs in HLB infection of citrus.