The influences of PRG-1 on the expression of small RNAs and mRNAs

Background

In metazoans, Piwi-related Argonaute proteins play important roles in maintaining germline integrity and fertility and have been linked to a class of germline-enriched small RNAs termed piRNAs. Caenorhabditis elegans encodes two Piwi family proteins called PRG-1 and PRG-2, and PRG-1 interacts with the C. elegans piRNAs (21U-RNAs). Previous studies found that mutation of prg-1 causes a marked reduction in the expression of 21U-RNAs, temperature-sensitive defects in fertility and other phenotypic defects.

Results

In this study, we wanted to systematically demonstrate the function of PRG-1 in the regulation of small RNAs and their targets. By analyzing small RNAs and mRNAs with and without a mutation in prg-1 during C. elegans development, we demonstrated that (1) mutation of prg-1 leads to a decrease in the expression of 21U-RNAs, and causes 35 ~ 40% of miRNAs to be down-regulated; (2) in C. elegans, approximately 3% (6% in L4) of protein-coding genes are differentially expressed after mutating prg-1, and 60 ~ 70% of these substantially altered protein-coding genes are up-regulated; (3) the target genes of the down-regulated miRNAs and the candidate target genes of the down-regulated 21U-RNAs are enriched in the up-regulated protein-coding genes; and (4) PRG-1 regulates protein-coding genes by down-regulating small RNAs (miRNAs and 21U-RNAs) that target genes that participate in the development of C. elegans.

Conclusions

In prg-1-mutated C. elegans, the expression of miRNAs and 21U-RNAs was reduced, and the protein-coding targets, which were associated with the development of C. elegans, were up-regulated. This may be the mechanism underlying PRG-1 function.

Electronic supplementary material

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

An atlas and analysis of bovine skeletal muscle long noncoding RNAs

Long noncoding RNAs (lncRNAs) have various biological functions and have been extensively studied in recent years. However, the identification and characterization of bovine lncRNAs in skeletal muscle has been very limited compared with that of lncRNAs in other model organisms. In this study, 7188 bovine skeletal muscle lncRNAs were identified by RNA‐Seq and a stringent screening procedure in four different muscle tissues. These lncRNAs shared many characteristics with other mammalian lncRNAs, such as a shorter open reading frame and lower expression level than for mRNAs. Furthermore, the chromosomal locations and global expression patterns for these lncRNAs are also described in detail. More importantly, we detected the important interaction relationships of lncRNAs–miRNAs–mRNAs related to muscle development among 36 lncRNAs, 62 miRNAs and 12 mRNAs. Our results provide a global expression pattern of lncRNAs specific to bovine skeletal muscle and provide important targets for revealing the function of bovine muscle development by thoroughly studying the interaction relationships of lncRNAs–miRNAs–mRNAs.     

Identification and characterization of endogenous small interfering RNAs from rice

RNA silencing-mediated small interfering RNAs (siRNAs) and microRNAs (miRNAs) have diverse natural roles, ranging from regulation of gene expression and heterochromatin formation to genome defense against transposons and viruses. Unlike miRNAs, endogenous siRNAs are generally not conserved between species; consequently, their identification requires experimental approaches. Thus far, endogenous siRNAs have not been reported from rice, which is a model species for monocotyledonous plants. We identified a large set of putative endogenous siRNAs from root, shoot and inflorescence small RNA cDNA libraries of rice. Most of these siRNAs are from intergenic regions, although a substantial proportion (22%) originates from the introns and exons of protein-coding genes. Northern and RT–PCR analysis revealed that the expression of some of the siRNAs is tissue specific or developmental stage specific. A total of 25 transposons and 21 protein-coding genes were predicted to be cis-targets of some of the siRNAs. Based on sequence homology, we also predicted 111 putative trans-targets for 44 of the siRNAs. Interestingly, ~46% of the predicted trans-targets are transposable elements, which suggests that endogenous siRNAs may play an important role in the suppression of transposon proliferation. Using RNA ligase-mediated-5′ rapid amplification of cDNA end assays, we validated three of the predicted targets and provided evidence for both cis- and trans-silencing of target genes by siRNAs-guided mRNA cleavage.     

Identifying expression of new small RNAs by microarrays

Although a large number of small RNAs (sRNAs) have been discovered, it is very likely that the screens conducted so far have not reached saturation. Recently, many methods for predicting and identifying new sRNAs have been developed. However, it remains unclear what the total number of sRNAs within a genome is and how many types of sRNAs exist in plants and animals. In this article, combined methods of dynamic programming prediction, enrichment of sRNAs, and microarray analysis are developed to screen and evaluate a new class of sRNAs from introns of human, protein-encoding genes. The methods used by our laboratories to design capture probes and label enriched small RNAs are thoroughly described here. The microarray results show that our modified technologies are useful to enhance sensitivity and specificity of arrays, identify expression patterns within different cells, and discover differential expression of sRNAs during the differentiation process of bone marrow stem cells. Accordingly, the combination of computational prediction and microarray analysis may be a feasible and practical approach for profiling studies of both known and predicted small RNAs.     

Cloning and expression analysis of piRNA-like RNAs: adult testis-specific small RNAs in chicken

Many small RNAs have been cloned from animal gonads, for example, endogenous small interfering RNAs (endo-siRNAs) were found in oocytes and piwi-interacting RNAs (piRNAs) were found in testis. Gallus gallus (chicken) is an important model organism, but few small RNAs have been identified from its gonads. In this study, we isolated and cloned 156 small RNAs from adult chicken testes. Since there is a reasonably even distribution from 22 to 33 nt, these small RNAs are slightly longer than miRNAs and endo-siRNAs. Genome mapping indicated that these small RNAs were derived from intergenic regions, exons, introns, and repetitive elements including chicken repeat 1, long terminal repeats, and simple repeats. Since they are similar with piRNAs, we named them piRNA-like RNAs (pil-RNAs). Northern blotting of 16 selected sequences showed that nine are specifically expressed in the adult testis. The vast majority of these pil-RNAs are poorly conserved between species, suggesting that they are unique to the adult chicken testis. Further analysis of the cloned pil-RNAs will improve our understanding of the function of small RNAs in animal gonad development.     

Identification of small RNAs in late developmental stage of rice anthers

Small RNAs including microRNA (miRNA) and small interfering RNA (siRNA) are known as repressors of gene expression. There are many plant proteins involved in small RNA-mediated gene silencing, such as Dicer ribonucleases and RNA-dependent RNA polymerases. However, most of these proteins have been reported to be absent in the late developmental stage of the plant male gamete, pollen. In order to clarify the existence of the small RNAs during maturation of pollen, we cloned and sequenced small RNAs from rice anthers including tricellular pollen. From fifty six candidates of small RNAs, we identified two known miRNAs (miR166 and miR167), eight potential miRNAs, and ten putative heterochromatic siRNAs (hc-siRNAs). RNA gel blot analyses clearly showed that miR166 and miR167 were accumulated in the uninuclear pollen stage of anther development and remained until the tricellular pollen stage. Our cloning and RNA gel blot analyses of small RNAs led us to propose a possible function of small RNA-mediated gene regulation for the development of male gametes in rice.     

Sex specific expression and distribution of small RNAs in papaya

Background

Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.

Results

We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.

Conclusion

By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Y^h chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.

Electronic supplementary material

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

An expression atlas of miRNAs in Arabidopsis thaliana

MicroRNAs(miRNAs) are small non-coding RNAs that regulate a variety of biological processes. miRNA expression often exhibits spatial and temporal specificity. However, genome-wide miRNA expression patterns in different organs during development of Arabidopsis thaliana have not yet been systemically investigated. In this study, we sequenced small RNA libraries generated from 27 different organ/tissue types, which cover the entire life cycle of Arabidopsis. Analysis of the sequencing data revealed that most miRNAs are ubiquitously expressed, whereas a small set of miRNAs display highly specific expression patterns. In addition, different miRNA members within the same family have distinct spatial and temporal expression patterns. Moreover, we found that some miRNAs are produced from different arms of their hairpin precursors at different developmental stages. This work provides new insights into the regulation of miRNA biogenesis and a rich resource for future investigation of miRNA functions in Arabidopsis.     

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA–target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA–target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA–target interactions.

A population of rice short interfering RNAs derived from TRANS-ACTING SIRNA3 loci direct target mRNA cleavage, revealing a broader role of TRANS-ACTING SIRNA3 genes in rice growth and stress responses.     

Characterization and expression patterns of small RNAs in synthesized Brassica hexaploids

Polyploidy has played an important role in promoting plant evolution through genomic merging and doubling. We used high-throughput sequencing to compare miRNA expression profiles between Brassica hexaploid and its parents. A total of 613, 784 and 742 known miRNAs were identified in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. We detected 618 miRNAs were differentially expressed (log₂Ratio ≥ 1, P ≤ 0.05) between Brassica hexaploid and its parents, and 425 miRNAs were non-additively expressed in Brassica hexaploid, which suggest a trend of non-additive miRNA regulation following hybridization and polyploidization. Remarkably, majority of the non-additively expressed miRNAs in the Brassica hexaploid are repressed, and there was a bias toward repression of B. rapa miRNAs, which is consistent with the progenitor-biased gene repression in the synthetic allopolyploids. In addition, we identified 653 novel mature miRNAs in Brassica hexaploid and its parents. Finally, we found that almost all the non-additive accumulation of siRNA clusters exhibited a low-parent pattern in Brassica hexaploid. Non-additive small RNA regulation is involved in a range of biological pathways, probably providing a driving force for variation and adaptation in allopolyploids.     

An expression atlas of connexin genes in the mouse

Connexin genes are involved in several human diseases such as hearing and dermatological and peripheral nerve disorders. Connexins are protein units of gap junctions and form homotypic, heterotypic, or heteromeric complexes known as connexons. Data on the expression patterns of members of this family are partial and fragmentary. We therefore cloned all the identifiable murine homologs of human CONNEXIN genes and analyzed their expression patterns in embryonic and neonatal mouse tissues. We found that connexins are preferentially expressed in tissues derived from ectoderm and/or endoderm. Our data provide a comprehensive and detailed atlas of expression of connexin genes and in some cases suggest possible interactions of proteins that are coexpressed in the same tissue. Knowledge of temporal and spatial distribution of connexins also allows the identification of candidate genes for human diseases and provides important insight into mechanisms that lead to human disorders due to mutations in CONNEXIN genes.