Discovery of potential piRNAs from next generation sequences of the sexually mature porcine testes

cDNA Cloning of Feline PIWIL1 and Evaluation of Expression in the Testis of the Domestic Cat

The PIWI clade of Argonaute proteins is essential for spermatogenesis in all species examined to date. This protein family binds specific classes of small non-coding RNAs known as PIWI-interacting RNAs (piRNAs) which together form piRNA-induced silencing complexes (piRISCs) that are recruited to specific RNA targets through sequence complementarity. These complexes facilitate gene silencing through endonuclease activity and guided recruitment of epigenetic silencing factors. PIWI proteins and piRNAs have been found to play multiple roles in the testis including the maintenance of genomic integrity through transposon silencing and facilitating the turnover of coding RNAs during spermatogenesis. In the present study, we report the first characterization of PIWIL1 in the male domestic cat, a mammalian system predicted to express four PIWI family members. Multiple transcript variants of PIWIL1 were cloned from feline testes cDNA. One isoform shows high homology to PIWIL1 from other mammals, however, the other has characteristics of a "slicer null" isoform, lacking the domain required for endonuclease activity. Expression of PIWIL1 in the male cat appears limited to the testis and correlates with sexual maturity. RNA-immunoprecipitation revealed that feline PIWIL1 binds small RNAs with an average size of 29 nt. Together, these data suggest that the domestic cat has two PIWIL1 isoforms expressed in the mature testis, at least one of which interacts with piRNAs.

The PIWI/piRNA response is relaxed in a rodent that lacks mobilizing transposable elements

Transposable elements (TEs) are genomic parasites that can propagate throughout host genomes. Mammalian genomes are typically dominated by LINE retrotransposons and their associated SINEs, and germline mobilization is a challenge to genome integrity. There are defenses against TE proliferation and the PIWI/piRNA defense is among the most well understood. However, the PIWI/piRNA system has been investigated largely in animals with actively mobilizing TEs and it is unclear how the PIWI/piRNA system functions in the absence of mobilizing TEs. The 13-lined ground squirrel provides the opportunity to examine PIWI/piRNA and TE dynamics within the context of minimal, and possibly nonexistent, TE accumulation. To do so, we compared the PIWI/piRNA dynamics in squirrels to observations from the rabbit and mouse. Despite a lack of young insertions in squirrels, TEs were still actively transcribed at higher levels compared to mouse and rabbit. All three Piwi genes were not expressed, prior to P8 in squirrel testis, and there was little TE expression change with the onset of Piwi expression. We also demonstrated there was not a major expression change in the young squirrel LINE families in the transition from juvenile to adult testis in contrast to young mouse and rabbit LINE families. These observations lead us to conclude that PIWI suppression, was weaker for squirrel LINEs and SINEs and did not strongly reduce their transcription. We speculate that, although the PIWI/piRNA system is adaptable to novel TE threats, transcripts from TEs that are no longer threatening receive less attention from PIWI proteins.

New transcriptomic tools to understand testis development and functions

The testis plays a central role in the male reproductive system - secreting several hormones including male steroids and producing male gametes. A complex and coordinated molecular program is required for the proper differentiation of testicular cell types and maintenance of their functions in adulthood. The testicular transcriptome displays the highest levels of complexity and specificity across all tissues in a wide range of species. Many studies have used high-throughput sequencing technologies to define the molecular dynamics and regulatory networks in the testis as well as to identify novel genes or gene isoforms expressed in this organ. This review intends to highlight the complementarity of these transcriptomic studies and to show how the use of different sequencing protocols contribute to improve our global understanding of testicular biology.

Identification and verification of potential piRNAs from domesticated yak testis

PIWI-interacting RNAs (piRNA) are small non-coding RNA molecules expressed in animal germ cells that interact with PIWI family proteins to form RNA–protein complexes involved in epigenetic and post-transcriptional gene silencing of retrotransposons and other genetic elements in germ line cells, including reproductive stem cell self-sustainment, differentiation, meiosis and spermatogenesis. In the present study, we performed high-throughput sequencing of piRNAs in testis samples from yaks in different stages of sexual maturity. Deep sequencing of the small RNAs (18–40 nt in length) yielded 4,900,538 unique reads from a total of 53,035,635 reads. We identified yak small RNAs (18–30 nt) and performed functional characterization. Yak small RNAs showed a bimodal length distribution, with two peaks at 22 nt and >28 nt. More than 80% of the 3,106,033 putative piRNAs were mapped to 4637 piRNA-producing genomic clusters using RPKM. 6388 candidate piRNAs were identified from clean reads and the annotations were compared with the yak reference genome repeat region. Integrated network analysis suggested that some differentially expressed genes were involved in spermatogenesis through ECM–receptor interaction and PI3K-Akt signaling pathways. Our data provide novel insights into the molecular expression and regulation similarities and diversities in spermatogenesis and testicular development in yaks at different stages of sexual maturity.

Cellular and molecular responses of adult testis to changes in nutrition: novel insights from the sheep model

This review explores the cellular and molecular mechanisms that regulate spermatogenesis in the post-pubertal testis that is regressing in response to mild undernutrition, using the sexually mature male sheep as a model. Testis regression leads to reductions in daily sperm production and in the quality of ejaculated spermatozoa (poorer movement, DNA damage). There is also a reduction in spermatogenic efficiency that appears to be caused, at least partially, by increases in germ cell apoptosis. Sertoli cell number does not change with testis regression, although about 1% of Sertoli cells do appear to retain proliferative ability after puberty. On the other hand, Sertoli cell function is disrupted during testis regression, as evidenced by a disorganization of tight junctions and indications that cell differentiation and maturation are reversed. Disrupted Sertoli cell function can explain, at least partially, the increase in germ cell apoptosis and any decrease in the rate of spermatogenesis, the two major contributors to spermatogenic efficiency. These outcomes seem to be mediated by changes in two RNA-based processes: (i) the expression of small non-coding RNAs that are involved in the regulation of Sertoli cell function, spermatogenesis and germ cell apoptosis and (ii) alternative pre-mRNA splicing that affects the regulation of spermatogenesis but does not appear to affect germ cell apoptosis, at least during testis progression induced by undernutrition in the male sheep. These research outcomes can be extended to other animal models and are relevant to issues in human male fertility.

Genome-wide differential expression of genes and small RNAs in testis of two different porcine breeds and at two different ages

Some documented evidences proved small RNAs (sRNA) and targeted genes are involved in mammalian testicular development and spermatogenesis. However, the detailed molecular regulation mechanisms of them remain largely unknown so far. In this study, we obtained a total of 10,716 mRNAs, 67 miRNAs and 16,953 piRNAs which were differentially expressed between LC and LW pig breeds or between the two sexual maturity stages. Of which, we identified 16 miRNAs and 28 targeted genes possibly related to spermatogenesis; 14 miRNA and 18 targeted genes probably associated with cell adhesion related testis development. We also annotated 579 piRNAs which could potentially regulate cell death, nucleosome organization and other basic biology process, which implied that those piRNAs might be involved in sexual maturation difference. The integrated network analysis results suggested that some differentially expressed genes were involved in spermatogenesis through the ECM-receptor interaction, focal adhesion, Wnt and PI3K-Akt signaling pathways, some particular miRNAs have the negative regulation roles and some special piRNAs have the positive and negative regulation roles in testicular development. Our data provide novel insights into the molecular expression and regulation similarities and diversities of spermatogenesis and testicular development in different pig breeds at different stages of sexual maturity.

Transposable Element Targeting by piRNAs in Laurasiatherians with Distinct Transposable Element Histories

PIWI proteins and PIWI-interacting RNAs (piRNAs) are part of a cellular pathway that has evolved to protect genomes against the proliferation of transposable elements (TEs). PIWIs and piRNAs assemble into complexes that are involved in epigenetic and post-transcriptional repression of TEs. Most of our understanding of the mechanisms of piRNA-mediated TE silencing comes from fruit fly and mouse models. However, even in these well-studied animals it is unclear how piRNA responses relate to variable TE expression and whether the strength of the piRNA response affects TE content over time. Here, we assessed the evolutionary interactions between TE and piRNAs in a statistical framework using three nonmodel laurasiatherian mammals as a study system: dog, horse, and a vesper bat. These three species diverged ∼80 million years ago and have distinct genomic TE contents. By comparing species with distinct TE landscapes, we aimed to identify clear relationships among TE content, expression, and piRNAs. We found that the TE subfamilies that are the most transcribed appear to elicit the strongest “ping-pong” response. This was most evident among long interspersed elements, but the relationships between expression and ping-pong pilRNA (piRNA-like) expression were more complex among SINEs. SINE transcripts were equally abundant in the dog and horse yet new SINE insertions were relatively rare in the horse genome, where we identified a stronger piRNA response. Our analyses suggest that the piRNA response can have a strong impact on the TE composition of a genome. However, our results also suggest that the presence of a robust piRNA response is apparently not sufficient to stop TE mobilization and accumulation.

Discovery of piRNAs Pathway Associated with Early-Stage Spermatogenesis in Chicken

Piwi-interacting RNAs (piRNAs) play a key role in spermatogenesis. Here, we describe the piRNAs profiling of primordial germ cells (PGCs), spermatogonial stem cells (SSCs), and the spermatogonium (Sp) during early-stage spermatogenesis in chicken. We obtained 31,361,989 reads from PGCs, 31,757,666 reads from SSCs, and 46,448,327 reads from Sp cells. The length distribution of piRNAs in the three samples showed peaks at 33 nt. The resulting genes were subsequently annotated against the Gene Ontology (GO) database. Five genes (RPL7A, HSPA8, Pum1, CPXM2, and PRKCA) were found to be involved in cellular processes. Interactive pathway analysis (IPA) further revealed three important pathways in early-stage spermatogenesis including the FGF, Wnt, and EGF receptor signaling pathways. The gene Pum1 was found to promote germline stem cell proliferation, but it also plays a role in spermatogenesis. In conclusion, we revealed characteristics of piRNAs during early spermatogonial development in chicken and provided the basis for future research.

Abundant and Altered Expression of PIWI-Interacting RNAs during Cardiac Hypertrophy

BACKGROUND

The discovery of PIWI-interacting RNAs (piRNAs) has fundamentally changed our understanding of post transcriptional regulation of transposons and other genes. Unlike miRNA and siRNA, the piRNAs are the most abundant but least studied RNA species in mammals. Although the expression of PIWI proteins and piRNAs has long been regarded as germline specific, increasing evidences suggest the expression of piRNAs in somatic cells.

METHODS

In this study, the small RNA sequencing executed during induction of cardiac hypertrophy in both in vivo and in vitro conditions were annotated for the expression of piRNAs. The expression of piRNAs was validated by qPCR and RNA immunoprecipitation. In addition, the presence of piRNAs in circulation of myocardial infarction patients was studied by qPCR.

RESULTS

We identified an abundant and altered expression of piRNAs during cardiac hypertrophy. The differentially expressed piRNAs was validated by qPCR and RNA immunoprecipitation. The significantly and differentially expressed piRNAs were predicted to target different retrotransposons and mRNAs in the rat genome. The detection of specific piRNA in serum of myocardial infarction patients suggests the potential of piRNA for diagnosis.

CONCLUSION

Overall this study is the first to provide a whole-genome analysis of the large repertoire of piRNAs in the cardiac system and this would pave a new path to understanding the molecular aetiology of piRNA and retrotransposons in the physiology and pathology of the cardiac system.

Roles of small RNAs in the effects of nutrition on apoptosis and spermatogenesis in the adult testis

We tested whether reductions in spermatozoal quality induced by under-nutrition are associated with increased germ cell apoptosis and disrupted spermatogenesis, and whether these effects are mediated by small RNAs. Groups of 8 male sheep were fed for a 10% increase or 10% decrease in body mass over 65 days. Underfeeding increased the number of apoptotic germ cells (P < 0.05) and increased the expression of apoptosis-related genes (P < 0.05) in testicular tissue. We identified 44 miRNAs and 35 putative piRNAs that were differentially expressed in well-fed and underfed males (FDR < 0.05). Some were related to reproductive system development, apoptosis (miRNAs), and sperm production and quality (piRNAs). Novel-miR-144 (miR-98), was found to target three apoptotic genes (TP53, CASP3, FASL). The proportion of miRNAs as a total of small RNAs was greater in well-fed males than in underfed males (P < 0.05) and was correlated (r = 0.8, P < 0.05) with the proportion of piRNAs in well-fed and underfed males. In conclusion, the reductions in spermatozoal quality induced by under-nutrition are caused, at least partly, by disruptions to Sertoli cell function and increased germ cell apoptosis, mediated by changes in the expression of miRNAs and piRNAs.

piRNAs from Pig Testis Provide Evidence for a Conserved Role of the Piwi Pathway in Post-Transcriptional Gene Regulation in Mammals

Piwi-interacting (pi-) RNAs guide germline-expressed Piwi proteins in order to suppress the activity of transposable elements (TEs). But notably, the majority of pachytene piRNAs in mammalian testes is not related to TEs. This raises the question of whether the Piwi/piRNA pathway exerts functions beyond TE silencing. Although gene-derived piRNAs were described many times, a possible gene-regulatory function was doubted due to the absence of antisense piRNAs. Here we sequenced and analyzed piRNAs expressed in the adult testis of the pig, as this taxon possesses the full set of mammalian Piwi paralogs while their spermatozoa are marked by an extreme fitness due to selective breeding. We provide an exhaustive characterization of porcine piRNAs and genomic piRNA clusters. Moreover, we reveal that both sense and antisense piRNAs derive from protein-coding genes, while exhibiting features that clearly show that they originate from the Piwi/piRNA-mediated post-transcriptional silencing pathway, commonly referred to as ping-pong cycle. We further show that the majority of identified piRNA clusters in the porcine genome spans exonic sequences of protein-coding genes or pseudogenes, which reveals a mechanism by which primary antisense piRNAs directed against mRNA can be generated. Our data provide evidence that spliced mRNAs, derived from such loci, are not only targeted by piRNAs but are also subject to ping-pong cycle processing. Finally, we demonstrate that homologous genes are targeted and processed by piRNAs in pig, mouse and human. Altogether, this strongly suggests a conserved role for the mammalian Piwi/piRNA pathway in post-transcriptional regulation of protein-coding genes, which did not receive much attention so far.

Characterization of Sus scrofa small non-coding RNAs present in both female and male gonads

Small non-coding RNAs (sncRNAs) are indispensable for proper germ cell development, emphasizing the need for greater elucidation of the mechanisms of germline development and regulation of this process by sncRNAs. We used deep sequencing to characterize three families of small non-coding RNAs (piRNAs, miRNAs, and tRFs) present in Sus scrofa gonads and focused on the small RNA fraction present in both male and female gonads. Although similar numbers of reads were obtained from both types of gonads, the number of unique RNA sequences in the ovaries was several times lower. Of the sequences detected in the testes, 2.6% of piRNAs, 9% of miRNAs, and 10% of tRFs were also present in the ovaries. Notably, the majority of the shared piRNAs mapped to ribosomal RNAs and were derived from clustered loci. In addition, the most abundant miRNAs present in the ovaries and testes are conserved and are involved in many biological processes such as the regulation of homeobox genes, the control of cell proliferation, and carcinogenesis. Unexpectedly, we detected a novel sncRNA type, the tRFs, which are 30-36-nt RNA fragments derived from tRNA molecules, in gonads. Analysis of S. scrofa piRNAs show that testes specific piRNAs are biased for 5' uracil but both testes and ovaries specific piRNAs are not biased for adenine at the 10th nucleotide position. These observations indicate that adult porcine piRNAs are predominantly produced by a primary processing pathway or other mechanisms and secondary piRNAs generated by ping-pong mechanism are absent.

Characterization and comparative profiling of ovarian microRNAs during ovine anestrus and the breeding season

Background

Seasonal estrus is a critical limiting factor of animal fecundity, and it involves changes in both ovarian biology and hormone secretion in different seasons. Previous studies indicate that two classes of small RNAs (miRNAs and piRNAs) play important regulatory roles in ovarian biology. To understand the roles of small RNA-mediated post-transcriptional regulation in ovine seasonal estrus, the variation in expression patterns of ovarian small RNAs during anestrus and the breeding season were analyzed using Solexa sequencing technology. In addition, reproductive hormone levels were determined during ovine anestrus and the breeding season.

Results

A total of 483 miRNAs (including 97 known, 369 conserved and 17 predicated novel miRNAs), which belong to 183 different miRNA families, were identified in ovaries of Tan sheep and Small Tail Han (STH) sheep. Compared with the three stages of the breeding season, 25 shared significantly differentially expressed (including 19 up- and six down-regulated) miRNAs were identified in ovine anestrus. KEGG Pathway analysis revealed that the target genes for some of the differentially expressed miRNAs were involved in reproductive hormone related pathways (e.g. steroid biosynthesis, androgen and estrogen metabolism and GnRH signaling pathway) as well as follicular/luteal development related pathways. Moreover, the expression of the differentially expressed miRNAs and most of their target genes were negatively correlated in the above pathways. Furthermore, the levels of estrogen, progesterone and LH in ovine anestrus were significantly lower than those in the breeding season. Combining the results of pathway enrichment analysis, expression of target genes and hormone measurement, we suggest that these differentially expressed miRNAs in anestrus might participate in attenuation of ovarian activity by regulating the above pathways. Besides miRNAs, a large and unexpectedly diverse set of piRNAs were also identified.

Conclusions

The miRNA profiles of ovine ovaries in anestrus were presented for the first time. The identification and characterization of miRNAs that are differentially expressed between ovine anestrus and the breeding season will help understanding of the role of miRNAs in the regulation of seasonal estrus, and provides candidates for determining miRNAs which could be potentially used to regulate ovine seasonal estrus.

Electronic supplementary material

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

piClust: a density based piRNA clustering algorithm

Piwi-interacting RNAs (piRNAs) are recently discovered, endogenous small non-coding RNAs. piRNAs protect the genome from invasive transposable elements (TE) and sustain integrity of the genome in germ cell lineages. Small RNA-sequencing data can be used to detect piRNA activations in a cell under a specific condition. However, identification of cell specific piRNA activations requires sophisticated computational methods. As of now, there is only one computational method, proTRAC, to locate activated piRNAs from the sequencing data. proTRAC detects piRNA clusters based on a probabilistic analysis with assumption of a uniform distribution. Unfortunately, we were not able to locate activated piRNAs from our proprietary sequencing data in chicken germ cells using proTRAC. With a careful investigation on data sets, we found that a uniform or any statistical distribution for detecting piRNA clusters may not be assumed. Furthermore, small RNA-seq data contains many different types of RNAs which was not carefully taken into account in previous studies. To improve piRNA cluster identification, we developed piClust that uses a density based clustering approach without assumption of any parametric distribution. In previous studies, it is known that piRNAs exhibit a strong tendency of forming piRNA clusters in syntenic regions of the genome. Thus, the density based clustering approach is effective and robust to the existence of non-piRNAs or noise in the data. In experiments with piRNA data from human, mouse, rat and chicken, piClust was able to detect piRNA clusters from total small RNA-seq data from germ cell lines, while proTRAC was not successful. piClust outperformed proTRAC in terms of sensitivity and running time (up to 200 folds). piClust is currently available as a web service at http://epigenomics.snu.ac.kr/piclustweb.

Prediction of transcribed PIWI-interacting RNAs from CHO RNAseq data

Chinese hamster ovary (CHO) cells are currently the most important mammalian host for the manufacture of biopharmaceuticals. To enhance our understanding of cellular processes, pathways, and the genetic setup of CHO cell lines, we predicted PIWI interacting RNAs (piRNAs) from small RNA sequencing data. Although piRNAs are the least understood class of small non-coding RNAs that mediate RNA silencing, it is believed that they play a pivotal role in protecting genome integrity by repressing transposable elements. Since genomic integrity is the key to prolonged stability of recombinant CHO cell lines, we characterized piRNA sequences and expression in six CHO cell lines by computational analysis of an existing small RNA sequencing dataset using proTRAC and the published CHO genome as reference. Here we present the result of this analysis consisting of 25,626 piRNAs and 540 piRNA clusters. Moreover we provide first evidence for differential piRNA expression in adherent and suspension-adapted CHO-K1 and DUKXB11 host cell lines as well as their recombinant derivatives, indicating that piRNAs might be tools for cell line development and engineering.

Altered expression of porcine Piwi genes and piRNA during development

Three Sus scrofa Piwi genes (Piwil1, Piwil2 and Piwil4) encoding proteins of 861, 985 and 853 aminoacids, respectively, were cloned and sequenced. Alignment of the Piwi proteins showed the high identity between Sus scrofa and Homo sapiens. Relative transcript abundance of porcine Piwil1, Piwil2 and Piwil4 genes in testes, ovaries and oocytes derived from sexually immature and mature animals was examined using Real-Time PCR. Expression of the three Piwi mRNAs was proved to be tissue specific and restricted exclusively to the gonads. In testes of adult pigs the highest relative transcript abundance was observed for the Sus scrofa Piwil1 gene. On the other hand, in testes of neonatal pigs the Piwil1 transcript level was over 2-fold reduced while the level of Piwil2 transcript was higher. As regards the expression of the Piwil4 transcript, its level was 34-fold elevated in testes of neonatal piglet when compared to adult male. In ovaries of prepubertal and pubertal female pigs transcript abundance of the three Piwi genes was significantly reduced in comparison with testes. However, similarly to testes, in ovaries of neonatal pigs the Piwil2 gene was characterized by the highest relative transcript abundance among the three Piwi genes analysed. In prepubertal and pubertal oocytes Piwil1 transcript was the most abundant whereas the expression of Piwil4 was undetectable. We also demonstrated that expression of piRNA occurs preferentially in the gonads of adult male and female pigs. Moreover, a piRNA subset isolated from ovaries was 2-3 nucleotides longer than the piRNA from testes.