piRNA Populations Research Articles

PNLDC1 catalysis and postnatal germline function are required for piRNA trimming, LINE1 silencing, and spermatogenesis in mice.

PIWI-interacting RNAs (piRNAs) play critical and conserved roles in transposon silencing and gene regulation in the animal germline. Three distinct piRNA populations are present during mouse spermatogenesis: fetal piRNAs in fetal/perinatal testes, pre-pachytene and pachytene piRNAs in postnatal testes. PNLDC1 is required for piRNA 3' end maturation in multiple species. However, whether PNLDC1 is the bona fide piRNA trimmer and the physiological role of 3' trimming of different piRNA populations in spermatogenesis in mammals remain unclear. Here, by inactivating Pnldc1 exonuclease activity in vitro and in mice, we reveal that the PNLDC1 trimmer activity is essential for spermatogenesis and male fertility. PNLDC1 catalytic activity is required for both fetal and postnatal piRNA 3' end trimming. Despite this, postnatal piRNA trimming but not fetal piRNA trimming is critical for LINE1 transposon silencing. Furthermore, conditional inactivation of Pnldc1 in postnatal germ cells causes LINE1 transposon de-repression and spermatogenic arrest in mice, indicating that germline-specific postnatal piRNA trimming is essential for transposon silencing and germ cell development. Our findings highlight the germ cell-intrinsic role of PNLDC1 and piRNA trimming in mammals to safeguard the germline genome and promote fertility.

PiRNA expression patterns in high vs. low fertility bovine sperm

PIWI-interacting RNAs (piRNAs) are 24–32 nucleotide RNA sequences primarily expressed in germ cells and developing embryos that suppress transposable element expression to protect genomic integrity during epigenetic reprogramming events. We characterized the expression of piRNA sequences and their encoding clusters in sperm samples from an idiopathic fertility model of Holstein bulls with high and low Sire Conception Rates. The piRNA populations were determined to be mostly similar between fertility conditions when investigated by principal component and differential expression analysis, suggesting that a high degree of conservation in the piRNA system is likely necessary for the production of viable sperm. Both fertility conditions demonstrated evidence of ‘ping-pong’ activity – a secondary biogenesis pathway associated with active transposable element targeting and suppression. Most sperm-borne piRNAs were between 29-30 nucleotides in length and originated from 226 clusters across the genome, with the exception of chromosome 20. Mapping analysis revealed abundant targeting of several transposable element families, suggesting a suppressive function of sperm piRNAs consistent with their established roles. Expression of genes targeted by sperm-borne piRNAs is significantly reduced throughout early embryogenesis compared to the mRNA population. Limited transposable element expression is known to be essential for spermatogenesis, thus epigenetic regulation of this pathway is likely to influence sperm quality and fertilizing capacity.

An extended Tudor domain within Vreteno interconnects Gtsf1L and Ago3 for piRNA biogenesis in Bombyx mori.

Piwi-interacting RNAs (piRNAs) direct PIWI proteins to transposons to silence them, thereby preserving genome integrity and fertility. The piRNA population can be expanded in the ping-pong amplification loop. Within this process, piRNA-associated PIWI proteins (piRISC) enter a membraneless organelle called nuage to cleave their target RNA, which is stimulated by Gtsf proteins. The resulting cleavage product gets loaded into an empty PIWI protein to form a new piRISC complex. However, for piRNA amplification to occur, the new RNA substrates, Gtsf-piRISC, and empty PIWI proteins have to be in physical proximity. In this study, we show that in silkworm cells, the Gtsf1 homolog BmGtsf1L binds to piRNA-loaded BmAgo3 and localizes to granules positive for BmAgo3 and BmVreteno. Biochemical assays further revealed that conserved residues within the unstructured tail of BmGtsf1L directly interact with BmVreteno. Using a combination of AlphaFold modeling, atomistic molecular dynamics simulations, and in vitro assays, we identified a novel binding interface on the BmVreteno-eTudor domain, which is required for BmGtsf1L binding. Our study reveals that a single eTudor domain within BmVreteno provides two binding interfaces and thereby interconnects piRNA-loaded BmAgo3 and BmGtsf1L.

ADAD2 functions in spermiogenesis and piRNA biogenesis in mice.

Adenosine deaminase domain containing 2 (ADAD2) is a testis-specific protein composed of a double-stranded RNA binding domain and a non-catalytic adenosine deaminase domain. A recent study showed that ADAD2 is indispensable for the male reproduction in mice. However, the detailed functions of ADAD2 remain elusive. This study aimed to investigate the cause of male sterility in Adad2 mutant mice and to understand the molecular functions of ADAD2. Adad2 homozygous mutant mouse lines, Adad2-/- and Adad2Δ/Δ , were generated by CRISPR/Cas9. Western blotting and immunohistochemistry were used to reveal the expression and subcellular localization of ADAD2. Co-immunoprecipitation tandem mass spectrometry was employed to determine the ADAD2-interacting proteins in mouse testes. RNA-sequencing analyses were carried out to analyze the transcriptome and PIWI-interacting RNA (piRNA) populations in wildtype and Adad2 mutant testes. Adad2-/- and Adad2Δ/Δ mice exhibit male-specific sterility because of abnormal spermiogenesis. ADAD2 interacts with multiple RNA-binding proteins involved in piRNA biogenesis, including MILI, MIWI, RNF17, and YTHDC2. ADAD2 co-localizes and forms novel granules with RNF17 in spermatocytes. Ablation of ADAD2 impairs the formation of RNF17 granules, decreases the number of cluster-derived pachytene piRNAs, and increases expression of ping-pong-derived piRNAs. In collaboration with RNF17 and other RNA-binding proteins in spermatocytes, ADAD2 directly or indirectly functions in piRNA biogenesis.

The small RNA universe of Capitella teleta

RNAi is an evolutionarily fluid mechanism with dramatically different activities across animal phyla. One major group where there has been little investigation is annelid worms. Here the small RNAs of the polychaete developmental model, Capitella teleta, are profiled across development. As is seen with nearly all animals, nearly 200 hundred microRNAs were found with 58 high confidence novel species. Greater miRNA diversity was associated with later stages consistent with differentiation of tissues. Outside miRNA, a distinct composition of other small RNA pathways was found. Unlike many invertebrates, an endogenous siRNA pathway was not observed, indicating pathway loss relative to basal planarians. No processively generated siRNA‐class RNAs could be found arising from dsRNA precursors. This has significant impact on RNAi technology development for this group of animals. Unlike the apparent absence of siRNAs, a significant population of piRNAs were observed. For many piRNAs phasing and ping pong biogenesis pathways were identified. Interestingly, piRNAs were found to be highly expressed during early development, suggesting a potential role in regulation in metamorphosis. Critically, the configuration of RNAi factors in C. teletais found in other annelids and mollusks, suggesting that similar biology is likely present in the wider clade. This study is the first providing comprehensive analysis of small RNAs in annelids.

The Small RNA Universe of Capitella teleta.

RNAi is an evolutionarily fluid mechanism with dramatically different activities across animal phyla. One major group where there has been little investigation is annelid worms. Here, the small RNAs of the polychaete developmental model Capitella teleta are profiled across development. As is seen with nearly all animals, nearly 200 microRNAs were found with 58 high-confidence novel species. Greater miRNA diversity was associated with later stages consistent with differentiation of tissues. Outside miRNA, a distinct composition of other small RNA pathways was found. Unlike many invertebrates, an endogenous siRNA pathway was not observed, indicating pathway loss relative to basal planarians. No processively generated siRNA-class RNAs could be found arising from dsRNA precursors. This has a significant impact on RNAi technology development for this group of animals. Unlike the apparent absence of siRNAs, a significant population of piRNAs was observed. For many piRNAs, phasing and ping-pong biogenesis pathways were identified. Interestingly, piRNAs were found to be highly expressed during early development, suggesting a potential role in regulation in metamorphosis. Critically, the configuration of RNAi factors in C. teleta is found in other annelids and mollusks, suggesting that similar biology is likely to be present in the wider clade. This study is the first in providing comprehensive analysis of small RNAs in annelids.

Cellular abundance shapes function in piRNA-guided genome defense.

Defense against genome invaders universally relies on RNA-guided immunity. Prokaryotic CRISPR-Cas and eukaryotic RNA interference pathways recognize targets by complementary base-pairing, which places the sequences of their guide RNAs at the center of self/nonself discrimination. Here, we explore the sequence space of PIWI-interacting RNAs (piRNAs), the genome defense of animals, and establish functional priority among individual sequences. Our results reveal that only the topmost abundant piRNAs are commonly present in every cell, whereas rare sequences generate cell-to-cell diversity in flies and mice. We identify a skewed distribution of sequence abundance as a hallmark of piRNA populations and show that quantitative differences of more than a 1000-fold are established by conserved mechanisms of biogenesis. Finally, our genomics analyses and direct reporter assays reveal that abundance determines function in piRNA-guided genome defense. Taken together, we identify an effective sequence space and untangle two classes of piRNAs that differ in complexity and function. The first class represents the topmost abundant sequences and drives silencing of genomic parasites. The second class sparsely covers an enormous sequence space. These rare piRNAs cannot function in every cell, every individual, or every generation but create diversity with potential for adaptation in the ongoing arms race with genome invaders.

PiRNAclusterDB 2.0: update and expansion of the piRNA cluster database.

PIWI-interacting RNAs (piRNAs) and their partnering PIWI proteins defend the animal germline against transposable elements and play a crucial role in fertility. Numerous studies in the past have uncovered many additional functions of the piRNA pathway, including gene regulation, anti-viral defense, and somatic transposon repression. Further, comparative analyses across phylogenetic groups showed that the PIWI/piRNA system evolves rapidly and exhibits great evolutionary plasticity. However, the presence of so-called piRNA clusters as the major source of piRNAs is common to nearly all metazoan species. These genomic piRNA-producing loci are highly divergent across taxa and critically influence piRNA populations in different evolutionary lineages. We launched the initial version of the piRNA cluster database to facilitate research on regulation and evolution of piRNA-producing loci across tissues und species. In recent years the amount of small RNA sequencing data that was generated and the abundance of species that were studied has grown rapidly. To keep up with this recent progress, we have released a major update for the piRNA cluster database (https://www.smallrnagroup.uni-mainz.de/piRNAclusterDB), expanding it from 12 to a total of 51 species with hundreds of new datasets, and revised its overall structure to enable easy navigation through this large amount of data.

A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons.

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

Reprogramming of Small Noncoding RNA Populations in Peripheral Blood Reveals Host Biomarkers for Latent and Active Mycobacterium tuberculosis Infection.

In tuberculosis (TB), as in other infectious diseases, studies of small noncoding RNAs (sncRNA) in peripheral blood have focused on microRNAs (miRNAs) but have neglected the other major sncRNA classes in spite of their potential functions in host gene regulation. Using RNA sequencing of whole blood, we have therefore determined expression of miRNA, PIWI-interacting RNA (piRNA), small nucleolar RNA (snoRNA), and small nuclear RNA (snRNA) in patients with TB (n = 8), latent TB infection (LTBI; n = 21), and treated LTBI (LTBItt; n = 6) and in uninfected exposed controls (ExC; n = 14). As expected, sncRNA reprogramming was greater in TB than in LTBI, with the greatest changes seen in miRNA populations. However, substantial dynamics were also evident in piRNA and snoRNA populations. One miRNA and 2 piRNAs were identified as moderately accurate (area under the curve [AUC] = 0.70 to 0.74) biomarkers for LTBI, as were 1 miRNA, 1 piRNA, and 2 snoRNAs (AUC = 0.79 to 0.91) for accomplished LTBI treatment. Logistic regression identified the combination of 4 sncRNA (let-7a-5p, miR-589-5p, miR-196b-5p, and SNORD104) as a highly sensitive (100%) classifier to discriminate TB from all non-TB groups. Notably, it reclassified 8 presumed LTBI cases as TB cases, 5 of which turned out to have features of Mycobacterium tuberculosis infection on chest radiographs. SNORD104 expression decreased during M. tuberculosis infection of primary human peripheral blood mononuclear cells (PBMC) and M2-like (P = 0.03) but not M1-like (P = 0.31) macrophages, suggesting that its downregulation in peripheral blood in TB is biologically relevant. Taken together, the results demonstrate that snoRNA and piRNA should be considered in addition to miRNA as biomarkers and pathogenesis factors in the various stages of TB.IMPORTANCE Tuberculosis is the infectious disease with the worldwide largest disease burden and there remains a great need for better diagnostic biomarkers to detect latent and active M. tuberculosis infection. RNA molecules hold great promise in this regard, as their levels of expression may differ considerably between infected and uninfected subjects. We have measured expression changes in the four major classes of small noncoding RNAs in blood samples from patients with different stages of TB infection. We found that, in addition to miRNAs (which are known to be highly regulated in blood cells from TB patients), expression of piRNA and snoRNA is greatly altered in both latent and active TB, yielding promising biomarkers. Even though the functions of many sncRNA other than miRNA are still poorly understood, our results strongly suggest that at least piRNA and snoRNA populations may represent hitherto underappreciated players in the different stages of TB infection.

Nanopore sequencing and Hi-C scaffolding provide insight into the evolutionary dynamics of transposable elements and piRNA production in wild strains of Drosophila melanogaster.

Illumina sequencing has allowed for population-level surveys of transposable element (TE) polymorphism via split alignment approaches, which has provided important insight into the population dynamics of TEs. However, such approaches are not able to identify insertions of uncharacterized TEs, nor can they assemble the full sequence of inserted elements. Here, we use nanopore sequencing and Hi-C scaffolding to produce de novo genome assemblies for two wild strains of Drosophila melanogaster from the Drosophila Genetic Reference Panel (DGRP). Ovarian piRNA populations and Illumina split-read TE insertion profiles have been previously produced for both strains. We find that nanopore sequencing with Hi-C scaffolding produces highly contiguous, chromosome-length scaffolds, and we identify hundreds of TE insertions that were missed by Illumina-based methods, including a novel micropia-like element that has recently invaded the DGRP population. We also find hundreds of piRNA-producing loci that are specific to each strain. Some of these loci are created by strain-specific TE insertions, while others appear to be epigenetically controlled. Our results suggest that Illumina approaches reveal only a portion of the repetitive sequence landscape of eukaryotic genomes and that population-level resequencing using long reads is likely to provide novel insight into the evolutionary dynamics of repetitive elements.

Organ-specific small non-coding RNA responses in domestic (Sudani) ducks experimentally infected with highly pathogenic avian influenza virus (H5N1)

ABSTRACT The duck represents an important reservoir of influenza viruses for transmission to other avian and mammalian hosts, including humans. The increased pathogenicity of the recently emerging clades of highly pathogenic avian influenza (HPAI) viruses of the H5N1 subtype in ducks features systemic viral spread and organ-to-organ variation in viral transcription and tissue damage. We previously reported that experimental infection of Sudani ducks (Cairina moschata) with an Egyptian HPAI (H5N1) virus (clade 2.2.1.2) features high viral replication and severe tissue damage in lung, but lower viral replication and only mild histological changes in brain. Little is known about the involvement of miRNA in organ-specific responses to H5N1 viruses in ducks, and involvement of the other classes of small noncoding RNA (sncRNA) has not been investigated so far. Following RNA sequencing, we have annotated the duck sncRNome and compared global expression changes of the four major sncRNA classes (miRNAs, piRNAs, snoRNAs, snRNAs) between duck lung and brain during a 120 h time course of infection with this HPAI strain. We find major organ-specific differences in miRNA, piRNA and snoRNA populations even before infection and substantial reprogramming of all sncRNA classes throughout infection, which was less pronounced in brain. Pathway prediction analysis of miRNA targets revealed enrichment of inflammation-, infection- and apoptosis-related pathways in lung, but enrichment of metabolism-related pathways (including tryptophan metabolism) in brain. Thus, organ-specific differences in sncRNA responses may contribute to differences in viral replication and organ damage in ducks infected with isolates from this emerging HPAI clade, and likely other strains.

Rapamycin-mediated mTOR inhibition impairs silencing of sex chromosomes and the pachytene piRNA pathway in the mouse testis.

Mechanistic target of rapamycin (mTOR) controls cell growth and metabolism in response to environmental and metabolic signals. Rapamycin robustly extends the lifespan in mammals and has clinical relevance in organ transplantation and cancer therapy but side effects include male infertility. Here, we report that chronic rapamycin treatment causes spermatogenic arrest in adult male mice due to defects in sex body formation and meiotic sex chromosome inactivation (MSCI). Many sex chromosome-linked genes were up-regulated in isolated pachytene spermatocytes from rapamycin-treated mice. RNA-Seq analysis also identified mRNAs encoding the core piRNA pathway components were decreased. Furthermore, rapamycin treatment was associated with a drastic reduction in pachytene piRNA populations. The inhibitory effects of rapamycin on spermatogenesis were partially reversible, with restoration of testis mass and sperm motility within 2 months of treatment cessation. Collectively, we have defined an essential role of mTOR in MSCI and identified a novel function as a regulator of small RNA homeostasis in male germ cells.

Mitochondrial membrane-based initial separation of MIWI and MILI functions during pachytene piRNA biogenesis.

PIWI-interacting RNAs (piRNAs) engage PIWI proteins to silence transposons and promote germ cell development in animals. In diverse species, piRNA biogenesis occurs near the mitochondrial surface, and involves mitochondrial membrane-anchored factors. In mice, two cytoplasmic PIWI proteins, MIWI and MILI, receive processed pachytene piRNAs at intermitochodrial cement (IMC). However, how MIWI and MILI are initially recruited to the IMC to engage multiple steps of piRNA processing is unclear. Here, we show that mitochondria-anchored TDRKH controls multiple steps of pachytene piRNA biogenesis in mice. TDRKH specifically recruits MIWI, but not MILI, to engage the piRNA pathway. It is required for the production of the entire MIWI-bound piRNA population and enables trimming of MILI-bound piRNAs. The failure to recruit MIWI to the IMC with TDRKH deficiency results in loss of MIWI in the chromatoid body, leading to spermiogenic arrest and piRNA-independent retrotransposon LINE1 de-repression in round spermatids. Our findings identify a mitochondrial surface-based scaffolding mechanism separating the entry and actions of two critical PIWI proteins in the same piRNA pathway to drive piRNA biogenesis and germ cell development.

An integrative piRNA analysis of mouse gametes and zygotes reveals new potential origins and gene regulatory roles

Piwi-interacting RNAs (piRNAs) are a subclass of the small non-coding RNAs (sncRNAs). Their main reported function was to exert control over transposable elements (TEs) in mammalian germline. In this study undertaking a deeper bioinformatics analysis of piRNAs present in mouse oocytes, sperm cells and zygotes, we first elaborated a new piRNA database based on sequences identified as piRNAs by immunoprecipitation with PIWI proteins. Our bioinformatics analysis revealed that, at least in gametes and zygotes, piRNAs could encompass multifunctional cell-dependent regulatory molecules. Indeed, genome analysis of the piRNA mapping density (reads/kb) evidenced in all samples an enrichment of intron-derived piRNAs. Further, piRNA population was classified into sequences not associated to TEs or repeats (NRapiRNAs) and associated to repetitive genome elements (RapiRNAs). In oocytes most of the NRapiRNAs mapped to the 5′UTRs of coding mRNAs, while higher proportion of NRapiRNAs was detected in sperm cells associated to the 3′UTRs of mRNAs. This piRNA complementarity to mRNA UTRs suggests key post-transcriptional regulatory roles over mRNAs such as those encoding MHC genes. In addition, a striking association of RapiRNA with long non-coding RNAs (lncRNAs) was identified. piRNAs associated with relevant lncRNAs such as: Rab26os and GAS5 and key mRNAs, were particularly assessed.

Elucidating the Small Regulatory RNA Repertoire of the Sea Anemone Anemonia viridis Based on Whole Genome and Small RNA Sequencing.

Cnidarians harbor a variety of small regulatory RNAs that include microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs), but detailed information is limited. Here, we report the identification and expression of novel miRNAs and putative piRNAs, as well as their genomic loci, in the symbiotic sea anemone Anemonia viridis. We generated a draft assembly of the A. viridis genome with putative size of 313 Mb that appeared to be composed of about 36% repeats, including known transposable elements. We detected approximately equal fractions of DNA transposons and retrotransposons. Deep sequencing of small RNA libraries constructed from A. viridis adults sampled at a natural CO2 gradient off Vulcano Island, Italy, identified 70 distinct miRNAs. Eight were homologous to previously reported miRNAs in cnidarians, whereas 62 appeared novel. Nine miRNAs were recognized as differentially expressed along the natural seawater pH gradient. We found a highly abundant and diverse population of piRNAs, with a substantial fraction showing ping–pong signatures. We identified nearly 22% putative piRNAs potentially targeting transposable elements within the A. viridis genome. The A. viridis genome appeared similar in size to that of other hexacorals with a very high divergence of transposable elements resembling that of the sea anemone genus Exaiptasia. The genome encodes and expresses a high number of small regulatory RNAs, which include novel miRNAs and piRNAs. Differentially expressed small RNAs along the seawater pH gradient indicated regulatory gene responses to environmental stressors.

Temperature-responsive miRNAs in Drosophila orchestrate adaptation to different ambient temperatures.

The majority of Drosophila genes are expressed in a temperature-dependent manner, but the way in which small RNAs may contribute to this effect is completely unknown as we currently lack an idea of how small RNA transcriptomes change as a function of temperature. Applying high-throughput sequencing techniques complemented by quantitative real-time PCR experiments, we demonstrate that altered ambient temperature induces drastic but reversible changes in sequence composition and total abundance of both miRNA and piRNA populations. Further, mRNA sequencing reveals that the expression of miRNAs and their predicted target transcripts correlates inversely, suggesting that temperature-responsive miRNAs drive adaptation to different ambient temperatures on the transcriptome level. Finally, we demonstrate that shifts in temperature affect both primary and secondary piRNA pools, and the observed aberrations are consistent with altered expression levels of the involved Piwi-pathway factors. We further reason that enhanced ping–pong processing at 29°C is driven by dissolved RNA secondary structures at higher temperatures, uncovering target sites that are not accessible at low temperatures. Together, our results show that small RNAs are an important part of epigenetic regulatory mechanisms that ensure homeostasis and adaptation under fluctuating environmental conditions.

Transposable Element Misregulation Is Linked to the Divergence between Parental piRNA Pathways in Drosophila Hybrids.

Interspecific hybridization is a genomic stress condition that leads to the activation of transposable elements (TEs) in both animals and plants. In hybrids between Drosophila buzzatii and Drosophila koepferae, mobilization of at least 28 TEs has been described. However, the molecular mechanisms underlying this TE release remain poorly understood. To give insight on the causes of this TE activation, we performed a TE transcriptomic analysis in ovaries (notorious for playing a major role in TE silencing) of parental species and their F1 and backcrossed (BC) hybrids. We find that 15.2% and 10.6% of the expressed TEs are deregulated in F1 and BC1 ovaries, respectively, with a bias toward overexpression in both cases. Although differences between parental piRNA (Piwi-interacting RNA) populations explain only partially these results, we demonstrate that piRNA pathway proteins have divergent sequences and are differentially expressed between parental species. Thus, a functional divergence of the piRNA pathway between parental species, together with some differences between their piRNA pools, might be at the origin of hybrid instabilities and ultimately cause TE misregulation in ovaries. These analyses were complemented with the study of F1 testes, where TEs tend to be less expressed than in D. buzzatii. This can be explained by an increase in piRNA production, which probably acts as a defence mechanism against TE instability in the male germline. Hence, we describe a differential impact of interspecific hybridization in testes and ovaries, which reveals that TE expression and regulation are sex-biased.

Bovine piRNA-like RNAs are associated with both transposable elements and mRNAs.

PIWI proteins and their associated piRNAs have been the focus of intensive research in the past decade; therefore, their participation in the maintenance of genomic integrity during spermatogenesis has been well established. Recent studies have suggested important roles for the PIWI/piRNA system outside of gametogenesis, based on the presence of piRNAs and PIWI proteins in several somatic tissues, cancers, and the early embryo. Here, we investigated the small RNA complement present in bovine gonads, gametes, and embryos through next-generation sequencing. A distinct piRNA population was present in the testis as expected. However, we also found a large population of slightly shorter, 24-27 nt piRNA-like RNA (pilRNAs) in pools of oocytes and zygotes. These oocyte and embryo pilRNAs exhibited many of the canonical characteristics of piRNAs including a 1U bias, the presence of a 'ping-pong' signature, genomic clustering, and transposable element targeting. Some of the major transposons targeted by oocyte and zygote pilRNA were from the LINE RTE and ERV1 classes. We also identified pools of pilRNA potentially derived from, or targeted at, specific mRNA sequences. We compared the frequency of these gene-associated pilRNAs to the fold change in the expression of respective mRNAs from two previously reported transcriptome datasets. We observed significant negative correlations between the number of pilRNAs targeting mRNAs, and their fold change in expression between the 4-8 cell and 8-16 cell stages. Together, these results represent one of the first characterizations of the PIWI/piRNA pathway in the translational bovine model, and in the novel context of embryogenesis.

Germ cell and tumor associated piRNAs in the medaka and Xiphophorus melanoma models.

BackgroundA growing number of studies report an abnormal expression of Piwi-interacting RNAs (piRNAs) and the piRNA processing enzyme Piwi in many cancers. Whether this finding is an epiphenomenon of the chaotic molecular biology of the fast dividing, neoplastically transformed cells or is functionally relevant to tumorigenesisis is difficult to discern at present. To better understand the role of piRNAs in cancer development small laboratory fish models can make a valuable contribution. However, little is known about piRNAs in somatic and neoplastic tissues of fish.ResultsTo identify piRNA clusters that might be involved in melanoma pathogenesis, we use several transgenic lines of medaka, and platyfish/swordtail hybrids, which develop various types of melanoma. In these tumors Piwi, is expressed at different levels, depending on tumor type. To quantify piRNA levels, whole piRNA populations of testes and melanomas of different histotypes were sequenced. Because no reference piRNA cluster set for medaka or Xiphophorus was yet available we developed a software pipeline to detect piRNA clusters in our samples and clusters were selected that were enriched in one or more samples. We found several loci to be overexpressed or down-regulated in different melanoma subtypes as compared to hyperpigmented skin. Furthermore, cluster analysis revealed a clear distinction between testes, low-grade and high-grade malignant melanoma in medaka.ConclusionsOur data imply that dysregulation of piRNA expression may be associated with development of melanoma. Our results also reinforce the importance of fish as a suitable model system to study the role of piRNAs in tumorigenesis.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2697-z) contains supplementary material, which is available to authorized users.