Discovery of Potential piRNAs from Next Generation Sequences of the Sexually Mature Porcine Testes

Piwi- interacting RNAs (piRNAs), a new class of small RNAs discovered from mammalian testes, are involved in transcriptional silencing of retrotransposons and other genetic elements in germ line cells. In order to identify a full transcriptome set of piRNAs expressed in the sexually mature porcine testes, small RNA fractions were extracted and were subjected to a Solexa deep sequencing. We cloned 6,913,561 clean reads of Sus Scrofa small RNAs (18–30 nt) and performed functional characterization. Sus Scrofa small RNAs showed a bimodal length distribution with two peaks at 21 nt and 29 nt. Then from 938,328 deep-sequenced small RNAs (26–30 nt), 375,195 piRNAs were identified by a k-mer scheme and 326 piRNAs were identified by homology searches. All piRNAs predicted by the k-mer scheme were then mapped to swine genome by Short Oligonucleotide Analysis Package (SOAP), and 81.61% of all uniquely mapping piRNAs (197,673) were located to 1124 defined genomic regions (5.85 Mb). Within these regions, 536 and 501 piRNA clusters generally distributed across only minus or plus genomic strand, 48 piRNA clusters distributed on two strands but in a divergent manner, and 39 piRNA clusters distributed on two strands in an overlapping manner. Furthermore, expression pattern of 7 piRNAs identified by homology searches showed 5 piRNAs displayed a ubiquitous expression pattern, although 2 piRNAs were specifically expressed in the testes. Overall, our results provide new information of porcine piRNAs and their specific expression pattern in porcine testes suggests that piRNAs have a role in regulating spermatogenesis

Human piRNAs Are Under Selection in Africans and Repress Transposable Elements

Piwi-interacting RNAs (piRNAs) are a recently discovered class of 24- to 30-nt noncoding RNAs whose best-understood function is to repress transposable elements (TEs) in animal germ lines. In humans, TE-derived sequences comprise ∼45% of the genome and there are several active TE families, including LINE-1 and Alu elements, which are a significant source of de novo mutations and intrapopulation variability. In the “ping-pong model,” piRNAs are thought to alternatively cleave sense and antisense TE transcripts in a positive feedback loop. Because piRNAs are poorly conserved between closely related species, including human and chimpanzee, we took a population genomics approach to study piRNA function and evolution. We found strong statistical evidence that piRNA sequences are under selective constraint in African populations. We then mapped the piRNA sequences to human TE sequences and found strong correlations between the age of each LINE-1 and Alu subfamily and the number of piRNAs mapping to the subfamily. This result supports the idea that piRNAs function as repressors of TEs in humans. Finally, we observed a significant depletion of piRNA matches in the reverse transcriptase region of the consensus human LINE-1 element but not of the consensus mouse LINE-1 element. This result suggests that reverse transcriptase might have an endogenous role specific to humans. Overall, our results elucidate the function and evolution of piRNAs in humans and highlight the utility of population genomics analysis for studying this rapidly evolving genetic system

piRNN: deep learning algorithm for piRNA prediction

Piwi-interacting RNAs (piRNAs) are the largest class of small non-coding RNAs discovered in germ cells. Identifying piRNAs from small RNA data is a challenging task due to the lack of conserved sequences and structural features of piRNAs. Many programs have been developed to identify piRNA from small RNA data. However, these programs have limitations. They either rely on extracting complicated features, or only demonstrate strong performance on transposon related piRNAs. Here we proposed a new program called piRNN for piRNA identification. For our software, we applied a convolutional neural network classifier that was trained on the datasets from four different species (Caenorhabditis elegans, Drosophila melanogaster, rat and human). A matrix of k-mer frequency values was used to represent each sequence. piRNN has great usability and shows better performance in comparison with other programs. It is freely available at https://github.com/bioinfolabmu/piRNN

Characterization of the Small RNA Transcriptome of the Marine Coccolithophorid, Emiliania huxleyi

Small RNAs (smRNAs) control a variety of cellular processes by silencing target genes at the transcriptional or post-transcription level. While extensively studied in plants, relatively little is known about smRNAs and their targets in marine phytoplankton, such as Emiliania huxleyi (E. huxleyi). Deep sequencing was performed of smRNAs extracted at different time points as E. huxleyi cells transition from logarithmic to stationary phase growth in batch culture. Computational analyses predicted 18 E. huxleyi specific miRNAs. The 18 miRNA candidates and their precursors vary in length (18-24 nt and 71-252 nt, respectively), genome copy number (3-1,459), and the number of genes targeted (2-107). Stem-loop real time reverse transcriptase (RT) PCR was used to validate miRNA expression which varied by nearly three orders of magnitude when growth slows and cells enter stationary phase. Stem-loop RT PCR was also used to examine the expression profiles of miRNA in calcifying and non-calcifying cultures, and a small subset was found to be differentially expressed when nutrients become limiting and calcification is enhanced. In addition to miRNAs, endogenous small RNAs such as ra-siRNAs, ta-siRNAs, nat-siRNAs, and piwiRNAs were predicted along with the machinery for the biogenesis and processing of si-RNAs. This study is the first genome-wide investigation smRNAs pathways in E. huxleyi. Results provide new insights into the importance of smRNAs in regulating aspects of physiological growth and adaptation in marine phytoplankton and further challenge the notion that smRNAs evolved with multicellularity, expanding our perspective of these ancient regulatory pathways

SMALL RNA EXPRESSION DURING PROGRAMMED REARRAGEMENT OF A VERTEBRATE GENOME

The sea lamprey (Petromyzon marinus) undergoes programmed genome rearrangements (PGRs) during embryogenesis that results in the deletion of ~0.5 Gb of germline DNA from the somatic lineage. The underlying mechanism of these rearrangements remains largely unknown. miRNAs (microRNAs) and piRNAs (PIWI interacting RNAs) are two classes of small noncoding RNAs that play important roles in early vertebrate development, including differentiation of cell lineages, modulation of signaling pathways, and clearing of maternal transcripts. Here, I utilized next generation sequencing to determine the temporal expression of miRNAs, piRNAs, and other small noncoding RNAs during the first five days of lamprey embryogenesis, a time series that spans the 24-32 cell stage to the formation of the neural crest. I obtained expression patterns for thousands of miRNA and piRNA species. These studies identified several thousand small RNAs that are expressed immediately before, during, and immediately after PGR. Significant sequence variation was observed at the 3’ end of miRNAs, representing template-independent covalent modifications. Patterns observed in lamprey are consistent with expectations that the addition of adenosine and uracil residues plays a role in regulation of miRNA stability during the maternal-zygotic transition. We also identified a conserved motif present in sequences without any known annotation that is expressed exclusively during PGR. This motif is similar to binding motifs of known DNA binding and nuclear export factors, and our data could represent a novel class of small noncoding RNAs operating in lamprey

Small non-coding RNAs and cancer: a bioinformatics approach

The development of Next-Generation Sequencing (NGS) technologies has added a new method to investigate the role of genes in several diseases, including breast cancer (BC). For example, RNA-Seq is mainly used to analyze the continuously changing cellular transcriptome and finding the pathways enriched/inhibited by the differently expressed genes. Small RNA-Seq is a very similar method allowing the detection of other classes of RNA such as small non-coding RNA (miRNA, piRNA, tRNA). In particular, microRNAs (miRNA) are a small non-coding RNA that regulate gene expression at the post-transcriptional level by binding to the 3'UTR (untranslated region) of target mRNAs, causing their degradation or translation inhibition. Several studies correlated an altered expression of these sncRNAs with the onset and /or progression of different diseases, including cancer. They can be used as biomarkers, playing a key role in the diagnosis, prognosis and prediction of response to specific therapies. Moreover, many studies have focused on the possibility of developing new therapeutic strategies based on microRNAs modulation and their potential use in the personalized management of cancer. Furthermore, also Piwi-interacting RNA (piRNA) is a class of small non-coding RNA molecules expressed in animal cells that are associated with both epigenetic and post-transcriptional gene silencing of retrotransposons and other genetic elements. They are distinct from miRNA in size (26–31 nt rather than 21–24 nt), lack of sequence conservation, and increased complexity. The aim of my research project was to understand the carcinogenic mechanisms and pathways modulated by these mRNAs and smallRNAs, since the clarification of the roles played by these molecules in cancer might provide new opportunities to develop novel strategies for diagnosing and treating this disease using bioinformatics tools and created it. Furthermore, a tool has been created that allows the analysis of smallRNAs by integrating various software to facilitate the use of this technology and to better explore the expected results

Identification and differential expression of piRNAs in the gonads of Amur sturgeon (Acipenser schrenckii)

Objective Sturgeons are considered living fossils, and have a very high conservation and economic value. Studies on the molecular mechanism of sturgeon gonadal development and sex differentiation would not only aid in understanding vertebrate sex determination but also benefit sturgeon aquaculture. Piwi-interacting RNAs (piRNAs) have been shown to function in germline or gonadal development. In this study, we performed small RNA deep sequencing and microarray hybridization to identify potential sturgeon piRNAs. Methods Male and female sturgeon gonads were collected and used for small RNA sequencing on an Illumina HiSeq platform with the validation of piRNA expression by microarray chip. The program Bowtie and k-mer scheme were performed to filter small RNA reads and discover potential sturgeon piRNAs. A known piRNA database, the coding sequence (CDS), 5′ and 3′ untranslated region (UTR) database of the A. Schrenckii transcriptome, Gene Ontology (GO) database and KEGG pathway database were searched subsequently to analyze the potential bio-function of sturgeon piRNAs. Results A total of 875,679 putative sturgeon piRNAs were obtained, including 93 homologous to known piRNAs and hundreds showing sex-specific and sex-biased expression. Further analysis showed that they are predominant in both the ovaries and testes and those with a sex-specific expression pattern are nearly equally distribution between sexes. This may imply a relevant role in sturgeon gonadal development. KEGG pathway and GO annotation analyses indicated that they may be related to sturgeon reproductive processes. Conclusion Our study provides the first insights into the gonadal piRNAs in a sturgeon species and should serve as a useful resource for further elucidation of the gene regulation involved in the sex differentiation of vertebrates. These results should also facilitate the technological development of early sex identification in sturgeon aquaculture

Bioinformatics and machine learning approaches to understand the regulation of mobile genetic elements

Transposable elements (TEs, or mobile genetic elements, MGEs) are ubiquitous genetic elements that make up a substantial proportion of the genome of many species. The recent growing interest in understanding the evolution and function of TEs has revealed that TEs play a dual role in genome evolution, development, disease, and drug resistance. Cells regulate TE expression against uncontrolled activity that can lead to developmental defects and disease, using multiple strategies, such as DNA chemical modification, small RNA (sRNA) silencing, chromatin modification, as well as sequence-specific repressors. Advancements in bioinformatics and machine learning approaches are increasingly contributing to the analysis of the regulation mechanisms. A plethora of tools and machine learning approaches have been developed for prediction, annotation, and expression profiling of sRNAs, for methylation analysis of TEs, as well as for genome-wide methylation analysis through bisulfite sequencing data. In this review, we provide a guided overview of the bioinformatic and machine learning state of the art of fields closely associated with TE regulation and function.peer-reviewe

Small non-coding RNA profiling and the role of piRNA pathway genes in the protection of chicken primordial germ cells

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.Background Genes, RNAs, and proteins play important roles during germline development. However, the functions of non-coding RNAs (ncRNAs) on germline development remain unclear in avian species. Recent high-throughput techniques have identified several classes of ncRNAs, including micro RNAs (miRNAs), small-interfering RNAs (siRNAs), and PIWI-interacting RNAs (piRNAs). These ncRNAs are functionally important in the genome, however, the identification and annotation of ncRNAs in a genome is challenging. The aim of this study was to identify different types of small ncRNAs particularly piRNAs, and the role of piRNA pathway genes in the protection of chicken primordial germ cells (PGCs). Results At first, we performed next-generation sequencing to identify ncRNAs in chicken PGCs, and we performed ab initio predictive analysis to identify putative piRNAs in PGCs. Then, we examined the expression of three repetitive sequence-linked piRNAs and 14 genic-transcript-linked piRNAs along with their linked genes using real-time PCR. All piRNAs and their linked genes were highly expressed in PGCs. Subsequently, we knocked down two known piRNA pathway genes of chicken, PIWI-like protein 1 (CIWI) and 2 (CILI), in PGCs using siRNAs. After knockdown of CIWI and CILI, we examined their effects on the expression of six putative piRNA-linked genes and DNA double-strand breakage in PGCs. The knockdown of CIWI and CILI upregulated chicken repetitive 1 (CR1) element and RAP2B, a member of RAS oncogene family, and increased DNA double-strand breakage in PGCs. Conclusions Our results increase the understanding of PGC-expressed piRNAs and the role of piRNA pathway genes in the protection of germ cells