NONCODEv4: exploring the world of long non-coding RNA genes

Exploring Consensus RNA Substructural Patterns Using Subgraph Mining

Frequently recurring RNA structural motifs play important roles in RNA folding process and interaction with other molecules. Traditional index-based and shape-based schemas are useful in modeling RNA secondary structures but ignore the structural discrepancy of individual RNA family member. Further, the in-depth analysis of underlying substructure pattern is insufficient due to varied and unnormalized substructure data. This prevents us from understanding RNAs functions and their inherent synergistic regulation networks. This article thus proposes a novel labeled graph-based algorithm RnaGraph to uncover frequently RNA substructure patterns. Attribute data and graph data are combined to characterize diverse substructures and their correlations, respectively. Further, a top-k graph pattern mining algorithm is developed to extract interesting substructure motifs by integrating frequency and similarity. The experimental results show that our methods assist in not only modelling complex RNA secondary structures but also identifying hidden but interesting RNA substructure patterns.

Proteogenomics from a bioinformatics angle: A growing field

Proteogenomics is a research area that combines areas as proteomics and genomics in a multi-omics setup using both mass spectrometry and high-throughput sequencing technologies. Currently, the main goals of the field are to aid genome annotation or to unravel the proteome complexity. Mass spectrometry based identifications of matching or homologues peptides can further refine gene models. Also, the identification of novel proteoforms is also made possible based on detection of novel translation initiation sites (cognate or near-cognate), novel transcript isoforms, sequence variation or novel (small) open reading frames in intergenic or un-translated genic regions by analyzing high-throughput sequencing data from RNAseq or ribosome profiling experiments. Other proteogenomics studies using a combination of proteomics and genomics techniques focus on antibody sequencing, the identification of immunogenic peptides or venom peptides. Over the years, a growing amount of bioinformatics tools and databases became available to help streamlining these cross-omics studies. Some of these solutions only help in specific steps of the proteogenomics studies, e.g. building custom sequence databases (based on next generation sequencing output) for mass spectrometry fragmentation spectrum matching. Over the last few years a handful integrative tools also became available that can execute complete proteogenomics analyses. Some of these are presented as stand-alone solutions, whereas others are implemented in a web-based framework such as Galaxy. In this review we aimed at sketching a comprehensive overview of all the bioinformatics solutions that are available for this growing research area. © 2015 Wiley Periodicals, Inc. Mass Spec Rev 36:584-599, 2017.

Genome-wide identification and characterization of putative lncRNAs in the diamondback moth, Plutella xylostella (L.)

Long non-coding RNAs (lncRNAs) are of particular interest because of their contributions to many biological processes. Here, we present the genome-wide identification and characterization of putative lncRNAs in a global insect pest, Plutella xylostella. A total of 8096 lncRNAs were identified and classified into three groups. The average length of exons in lncRNAs was longer than that in coding genes and the GC content was lower than that in mRNAs. Most lncRNAs were flanked by canonical splice sites, similar to mRNAs. Expression profiling identified 114 differentially expressed lncRNAs during the DBM development and found that majority were temporally specific. While the biological functions of lncRNAs remain uncharacterized, many are microRNA precursors or competing endogenous RNAs involved in micro-RNA regulatory pathways. This work provides a valuable resource for further studies on molecular bases for development of DBM and lay the foundation for discovery of lncRNA functions in P. xylostella.

The LncRNA Connectivity Map: Using LncRNA Signatures to Connect Small Molecules, LncRNAs, and Diseases

Well characterized the connections among diseases, long non-coding RNAs (lncRNAs) and drugs are important for elucidating the key roles of lncRNAs in biological mechanisms in various biological states. In this study, we constructed a database called LNCmap (LncRNA Connectivity Map), available at http://www.bio-bigdata.com/LNCmap/, to establish the correlations among diseases, physiological processes, and the action of small molecule therapeutics by attempting to describe all biological states in terms of lncRNA signatures. By reannotating the microarray data from the Connectivity Map database, the LNCmap obtained 237 lncRNA signatures of 5916 instances corresponding to 1262 small molecular drugs. We provided a user-friendly interface for the convenient browsing, retrieval and download of the database, including detailed information and the associations of drugs and corresponding affected lncRNAs. Additionally, we developed two enrichment analysis methods for users to identify candidate drugs for a particular disease by inputting the corresponding lncRNA expression profiles or an associated lncRNA list and then comparing them to the lncRNA signatures in our database. Overall, LNCmap could significantly improve our understanding of the biological roles of lncRNAs and provide a unique resource to reveal the connections among drugs, lncRNAs and diseases.

Identification of long non-coding RNAs in two anthozoan species and their possible implications for coral bleaching

Long non-coding RNAs (lncRNAs) have been shown to play regulatory roles in a diverse range of biological processes and are associated with the outcomes of various diseases. The majority of studies about lncRNAs focus on model organisms, with lessened investigation in non-model organisms to date. Herein, we have undertaken an investigation on lncRNA in two zoanthids (cnidarian): Protolpalythoa varibilis and Palythoa caribaeorum. A total of 11,206 and 13,240 lncRNAs were detected in P. variabilis and P. caribaeorum transcriptome, respectively. Comparison using NONCODE database indicated that the majority of these lncRNAs is taxonomically species-restricted with no identifiable orthologs. Even so, we found cases in which short regions of P. caribaeorum’s lncRNAs were similar to vertebrate species’ lncRNAs, and could be associated with lncRNA conserved regulatory functions. Consequently, some high-confidence lncRNA-mRNA interactions were predicted based on such conserved regions, therefore revealing possible involvement of lncRNAs in posttranscriptional processing and regulation in anthozoans. Moreover, investigation of differentially expressed lncRNAs, in healthy colonies and colonial individuals undergoing natural bleaching, indicated that some up-regulated lncRNAs in P. caribaeorum could posttranscriptionally regulate the mRNAs encoding proteins of Ras-mediated signal transduction pathway and components of innate immune-system, which could contribute to the molecular response of coral bleaching.

DES-ncRNA: A knowledgebase for exploring information about human micro and long noncoding RNAs based on literature-mining

Noncoding RNAs (ncRNAs), particularly microRNAs (miRNAs) and long ncRNAs (lncRNAs), are important players in diseases and emerge as novel drug targets. Thus, unraveling the relationships between ncRNAs and other biomedical entities in cells are critical for better understanding ncRNA roles that may eventually help develop their use in medicine. To support ncRNA research and facilitate retrieval of relevant information regarding miRNAs and lncRNAs from the plethora of published ncRNA-related research, we developed DES-ncRNA ( www.cbrc.kaust.edu.sa/des_ncrna ). DES-ncRNA is a knowledgebase containing text- and data-mined information from public scientific literature and other public resources. Exploration of mined information is enabled through terms and pairs of terms from 19 topic-specific dictionaries including, for example, antibiotics, toxins, drugs, enzymes, mutations, pathways, human genes and proteins, drug indications and side effects, mutations, diseases, etc. DES-ncRNA contains approximately 878,000 associations of terms from these dictionaries of which 36,222 (5,373) are with regards to miRNAs (lncRNAs). We provide several ways to explore information regarding ncRNAs to users including controlled generation of association networks as well as hypotheses generation. We show an example how DES-ncRNA can aid research on Alzheimer disease and suggest potential therapeutic role for Fasudil. DES-ncRNA is a powerful tool that can be used on its own or as a complement to the existing resources, to support research in human ncRNA. To our knowledge, this is the only knowledgebase dedicated to human miRNAs and lncRNAs derived primarily through literature-mining enabling exploration of a broad spectrum of associated biomedical entities, not paralleled by any other resource.

Distinct hepatitis B virus integration patterns in hepatocellular carcinoma and adjacent normal liver tissue

Infection by the hepatitis B virus (HBV) is one of the main etiologies of hepatocellular carcinoma (HCC). During chronic infection, HBV DNA can integrate into the human genome, and this has been postulated as a possible mechanism of HBV-induced HCC. In this study we used 2199 HBV integration sites from Dr.VIS v2.0 and mapped them to the human genome (hg19) to obtain viral integration sites (VIS) related to protein-coding and non-protein-coding genes. In total, we found 1,377 and 767 VIS within close proximity to protein coding genes and noncoding genes, respectively. Genes affected more than two times included 23.1% of protein-coding genes and 24.7% of long noncoding RNAs (lncRNA). Only 4.8% of VIS were shared between HCC and non-tumor tissues. HBV integrations were more common in chromosomes 5, 8, 10, and 19 in HCC tissue and chromosomes 1 and 2 in non-tumorous tissue. The number of integration sites on each chromosome correlated with the number of fragile sites in non-tumorous tissue but not in HCC tissue. Functional enrichment analysis of the protein-coding genes containing or in close proximity to HBV integration sites in HCC tissue showed an enrichment of cancer related gene ontology terms. Additionally, the most frequently associated lncRNA genes were related to telomere maintenance, protein modification processes, and chromosome localization. Thus, HBV may have preferred integration sites in the human genome that serve a critical role in HCC development. These results show that HCC treatment may benefit from the development of next generation anti-viral therapies.

Identification and characterization of long intergenic noncoding RNAs in bovine mammary glands

Background

Mammary glands of dairy cattle produce milk for the newborn offspring and for human consumption. Long intergenic noncoding RNAs (lincRNAs) play various functions in eukaryotic cells. However, types and roles of lincRNAs in bovine mammary glands are still poorly understood.

Results

Using computational methods, 886 unknown intergenic transcripts (UITs) were identified from five RNA-seq datasets from bovine mammary glands. Their non-coding potentials were predicted by using the combination of four software programs (CPAT, CNCI, CPC and hmmscan), with 184 lincRNAs identified. By comparison to the NONCODE2016 database and a domestic-animal long noncoding RNA database (ALDB), 112 novel lincRNAs were revealed in bovine mammary glands. Many lincRNAs were found to be located in quantitative trait loci (QTL). In particular, 36 lincRNAs were found in 172 milk related QTLs, whereas one lincRNA was within clinical mastitis QTL region. In addition, targeted genes for 10 lincRNAs with the highest fragments per kilobase of transcript per million fragments mapped (FPKM) were predicted by LncTar for forecasting potential biological functions of these lincRNAs. Further analyses indicate involvement of lincRNAs in several biological functions and different pathways.

Conclusion

Our study has provided a panoramic view of lincRNAs in bovine mammary glands and suggested their involvement in many biological functions including susceptibility to clinical mastitis as well as milk quality and production. This integrative annotation of mammary gland lincRNAs broadens and deepens our understanding of bovine mammary gland biology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3858-4) contains supplementary material, which is available to authorized users.

Identification of MyoD-Responsive Transcripts Reveals a Novel Long Non-coding RNA (lncRNA-AK143003) that Negatively Regulates Myoblast Differentiation

Myogenic differentiation factor (MyoD) is a master transcription factor in muscle development and differentiation. Although several long non-coding RNAs (lncRNAs) linked to MyoD have been found to influence muscle development, the functions of many lncRNAs have not been explored. Here we utilized lncRNA and mRNA microarray analysis to identify potential lncRNAs regulated by MyoD in muscle cells. A total of 997 differentially expressed lncRNAs (335 up-regulated and 662 down-regulated) and 1,817 differentially expressed mRNAs (148 up-regulated and 1,669 down-regulated) were identified after MyoD knockdown in C2C12 cells. Functional predictions suggested that most lncRNAs are involved in the biological pathways related to muscle differentiation and cell cycle with co-expressed genes. To gain further insight into the MyoD-mediated lncRNA expression in muscle differentiation, tissue expression profiles and MyoD overexpression were performed, and we found one of the candidate lncRNAs-AK143003 was significantly regulated by MyoD. Further analyses showed its noncoding ability and cytoplasmic localisation. Silencing of AK143003 stimulated the accumulation of myogenic marker genes, whereas AK143003 overexpression led to their decreased synthesis. This study identified a multitude of MyoD-mediated lncRNAs for further investigation and identified a novel lncRNA, lnc-AK143003, which plays a role in controlling muscle differentiation.

Unraveling long non-coding RNAs through analysis of high-throughput RNA-sequencing data

Extensive genome-wide transcriptome study mediated by high throughput sequencing technique has revolutionized the study of genetics and epigenetic at unprecedented resolution. The research has revealed that besides protein-coding RNAs, large proportions of mammalian transcriptome includes a heap of regulatory non protein-coding RNAs, the number encoded within human genome is enigmatic. Many taboos developed in the past categorized these non-coding RNAs as ''dark matter" and "junks". Breaking the myth, RNA-seq-- a recently developed experimental technique is widely being used for studying non-coding RNAs which has acquired the limelight due to their physiological and pathological significance. The longest member of the ncRNA family-- long non-coding RNAs, acts as stable and functional part of a genome, guiding towards the important clues about the varied biological events like cellular-, structural- processes governing the complexity of an organism. Here, we review the most recent and influential computational approach developed to identify and quantify the long non-coding RNAs serving as an assistant for the users to choose appropriate tools for their specific research.

Long Noncoding RNA CCAT2 as a Potential Novel Biomarker to Predict the Clinical Outcome of Cancer Patients: A Meta-Analysis

Background: Colon Cancer-Associated Transcript 2 (CCAT2) has been demonstrated associated with clinical outcomes in various tumors. However, the results from each study were unfortunately insufficient and not completely consistent. Therefore, we conduct a systematic meta-analysis to evaluate the value for a feasible biomarker for metastasis and prognosis.

Methods: A meta-analysis was performed using data obtained through a systematic search of PubMed, EMBASE, Cochrane Library, China National Knowledge Infrastructure, Wanfang database and VIP database. The pooled odds ratio (OR) and hazard ratio (HR) with 95% Confidence interval (CI ) using random-effect were used to identify the relationship of CCAT2 with clinical outcome of cancer patients. Subgroup analysis and sensitivity analysis were performed.

Results: A total of 867 patients from eight studies were finally included. Patients with high CCAT2 expression underwent an increased risk of lymph node metastasis (LNM) (OR=3.09, 95% CI: 1.53-6.26) and distant metastasis (DM) (OR=7.70, 95% CI: 3.26-18.17). CCAT2 was also significantly correlated with overall survival (OS) (HR=2.19, 95%CI: 1.70-2.82) and progression-free survival (PFS) (HR=2.59, 95% CI: 1.78-3.76). Moderate heterogeneity was observed in meta-analysis for LNM. However, the results remained robust in multiple sensitivity analyses.

Conclusions: High expression of CCAT2 was linked with poor clinical outcome. CCAT2 can serve as a potential molecular marker for prognosis in different types of cancers.

LncRNA-SVUGP2 suppresses progression of hepatocellular carcinoma

Numerous studies indicate that long noncoding RNAs (lncRNAs) are dysregulated in hepatocellular carcinoma (HCC) and might serve as potential diagnostic biomarkers and therapeutic targets of HCC. Therefore, it is interesting to globally identify the lncRNAs altered in HCC. In our study, we used microarray to profile the levels of lncRNAs and mRNAs in three pairs of HCC and their adjacent noncancerous samples. We found lncRNA-SVUGP2, which is a splice variant of the UGP2 gene, was down-regulated in HCC samples and correlates with a better prognosis in patients with HCC. Overexpression of lncRNA-SVUGP2 in HepG2 and Hep3B liver cancer cells suppresses cell proliferation in vitro and tumor growth in vivo. Moreover, lncRNA-SVUGP2 suppresses the invasion ability of liver cancer cell lines and downregulates the mRNA and protein levels of MMP2 and 9. Additionally, lncRNA-SVUGP2 positively or negatively correlates with many mRNAs in liver cancer tissues, indicating it is multifunctional in regulating carcinogenesis.

The origin of Metazoa: a unicellular perspective

The first animals evolved from an unknown single-celled ancestor in the Precambrian period. Recently, the identification and characterization of the genomic and cellular traits of the protists most closely related to animals have shed light on the origin of animals. Comparisons of animals with these unicellular relatives allow us to reconstruct the first evolutionary steps towards animal multicellularity. Here, we review the results of these investigations and discuss their implications for understanding the earliest stages of animal evolution, including the origin of metazoan genes and genome function.

Significant changes in circulating microRNA by dietary supplementation of selenium and coenzyme Q10 in healthy elderly males. A subgroup analysis of a prospective randomized double-blind placebo-controlled trial among elderly Swedish citizens

Background

Selenium and coenzyme Q10 is essential for important cellular functions. A low selenium intake is reported from many European countries, and the endogenous coenzyme Q10 production is decreasing in the body with increasing age. Supplementation with selenium and coenzyme Q10 in elderly have shown reduced cardiovascular mortality and reduced levels of markers of inflammation. However, microRNA analyses could give important information on the mechanisms behind the clinical effects of supplementation.

Methods

Out of the 443 healthy elderly participants that were given supplementation with 200 μg Se/day as organic selenium yeast tablets, and 200 mg/day of coenzyme Q10 capsules, or placebo for 4 years, 25 participants from each group were randomized and evaluated regarding levels of microRNA. Isolation of RNA from plasma samples and quantitative PCR analysis were performed. Volcano- and principal component analyses (PCA)–plots were used to illustrate the differences in microRNA expression between the intervention, and the placebo groups. Serum selenium concentrations were measured before intervention.

Findings

On average 145 different microRNAs out of 172 were detected per sample. In the PCA plots two clusters could be identified indicating significant difference in microRNA expression between the two groups. The pre-treatment expression of the microRNAs did not differ between active treatment and the placebo groups. When comparing the post-treatment microRNAs in the active and the placebo groups, 70 microRNAs exhibited significant differences in expression, also after adjustment for multiple measurements. For the 20 microRNAs with the greatest difference in expression the difference was up to more than 4 fold and with a P-value that were less than 4.4e^-8.

Conclusions

Significant differences were found in expression of more than 100 different microRNAs with up to 4 fold differences as a result of the intervention of selenium and coenzyme Q10 combined. The changes in microRNA could be a part of mechanisms underlying the clinical effects earlier reported that reduced cardiovascular mortality, gave better cardiac function, and showed less signs of inflammation and oxdative stress following the intervention. However, more research is needed to understand biological mechanisms of the protective effects of selenium and Q10 supplementation.

Review: Regulation of the cancer epigenome by long non-coding RNAs

Long non-coding RNAs have emerged as highly versatile players in the regulation of gene expression in development and human disease, particularly cancer. Hundreds of lncRNAs become dysregulated across tumor types, and multiple lncRNAs have demonstrated functions as tumor-suppressors or oncogenes. Furthermore, studies have demonstrated that dysregulation of lncRNAs results in alterations of the epigenome in cancer cells, potentially providing a novel mechanism for the massive epigenomic alterations observed in many tumors. Here, we highlight and provide some illustrious examples of lncRNAs in various epigenetic regulatory processes, including coordination of chromatin dynamics, regulation of DNA methylation, modulation of other non-coding RNAs and mRNA stability, and control of epigenetic substrate availability through altered tumor metabolism. In light of all these known and emerging functions in epigenetic regulation of tumorigenesis and cancer progression, lncRNAs represent attractive targets for future therapeutic strategies in cancer.

Identification of acylation products in SHAPE chemistry

SHAPE chemistry (selective 2'-hydroxyl acylation analyzed by primer extension) has been developed to specifically target flexible nucleotides (often unpaired nucleotides) independently to their purine or pyrimidine nature for RNA secondary structure determination. However, to the best of our knowledge, the structure of 2'-O-acylation products has never been confirmed by NMR or X-ray data. We have realized the acylation reactions between cNMP and NMIA under SHAPE chemistry conditions and identified the acylation products using standard NMR spectroscopy and LC-MS/MS experiments. For cAMP and cGMP, the major acylation product is the 2'-O-acylated compound (>99%). A trace amount of N-acylated cAMP has also been identified by LC-UV-MS². While for cCMP, the isolated acylation products are composed of 96% of 2'-O-acylated, 4% of N,O-diacylated, and trace amount of N-acylated compounds. In addition, the characterization of the major 2'-O-acylated compound by NMR showed slight differences in the conformation of the acylated sugar between the three cyclic nucleotides. This interesting result should be useful to explain some unexpected reactivity of the SHAPE chemistry.

Potentials of Long Noncoding RNAs (LncRNAs) in Sarcoma: From Biomarkers to Therapeutic Targets

Sarcoma includes some of the most heterogeneous tumors, which make the diagnosis, prognosis and treatment of these rare yet diverse neoplasms especially challenging. Long noncoding RNAs (lncRNAs) are important regulators of cancer initiation and progression, which implies their potential as neoteric prognostic and diagnostic markers in cancer, including sarcoma. A relationship between lncRNAs and sarcoma pathogenesis and progression is emerging. Recent studies demonstrate that lncRNAs influence sarcoma cell proliferation, metastasis, and drug resistance. Additionally, lncRNA expression profiles are predictive of sarcoma prognosis. In this review, we summarize contemporary advances in the research of lncRNA biogenesis and functions in sarcoma. We also highlight the potential for lncRNAs to become innovative diagnostic and prognostic biomarkers as well as therapeutic targets in sarcoma.

Present Scenario of Long Non-Coding RNAs in Plants

Small non-coding RNAs have been extensively studied in plants over the last decade. In contrast, genome-wide identification of plant long non-coding RNAs (lncRNAs) has recently gained momentum. LncRNAs are now being recognized as important players in gene regulation, and their potent regulatory roles are being studied comprehensively in eukaryotes. LncRNAs were first reported in humans in 1992. Since then, research in animals, particularly in humans, has rapidly progressed, and a vast amount of data has been generated, collected, and organized using computational approaches. Additionally, numerous studies have been conducted to understand the roles of these long RNA species in several diseases. However, the status of lncRNA investigation in plants lags behind that in animals (especially humans). Efforts are being made in this direction using computational tools and high-throughput sequencing technologies, such as the lncRNA microarray technique, RNA-sequencing (RNA-seq), RNA capture sequencing, (RNA CaptureSeq), etc. Given the current scenario, significant amounts of data have been produced regarding plant lncRNAs, and this amount is likely to increase in the subsequent years. In this review we have documented brief information about lncRNAs and their status of research in plants, along with the plant-specific resources/databases for information retrieval on lncRNAs.

Genome-wide screening and characterization of long non-coding RNAs involved in flowering development of trifoliate orange (Poncirus trifoliata L. Raf.)

Long non-coding RNAs (lncRNAs) have been demonstrated to play critical regulatory roles in post-transcriptional and transcriptional regulation in Arabidopsis. However, lncRNAs and their functional roles remain poorly characterized in woody plants, including citrus. To identify lncRNAs and investigate their role in citrus flowering, paired-end strand-specific RNA sequencing was performed for precocious trifoliate orange and its wild-type counterpart. A total of 6,584 potential lncRNAs were identified, 51.6% of which were from intergenic regions. Additionally, 555 lncRNAs were significantly up-regulated and 276 lncRNAs were down-regulated in precocious trifoliate orange, indicating that lncRNAs could be involved in the regulation of trifoliate orange flowering. Comparisons between lncRNAs and coding genes indicated that lncRNAs tend to have shorter transcripts and lower expression levels and that they display significant expression specificity. More importantly, 59 and 7 lncRNAs were identified as putative targets and target mimics of citrus miRNAs, respectively. In addition, the targets of Pt-miR156 and Pt-miR396 were confirmed using the regional amplification reverse-transcription polymerase chain reaction method. Furthermore, overexpression of Pt-miR156a1 and Pt-miR156a1 in Arabidopsis resulted in an extended juvenile phase, short siliques, and smaller leaves in transgenic plants compared with control plants. These findings provide important insight regarding citrus lncRNAs, thus enabling in-depth functional analyses.

Inhibition of long non-coding RNA UCA1 by CRISPR/Cas9 attenuated malignant phenotypes of bladder cancer

CRISPR/Cas9 is a novel and effective genome editing technique, but its application is not widely expanded to manipulate long non-coding RNA (lncRNA) expression. The lncRNA urothelial carcinoma-associated 1 (UCA1) is upregulated in bladder cancer and promotes the progression of bladder cancer. Here, we design gRNAs specific to UCA1 and construct CRISPR/Cas9 systems targeting UCA1. Single CRISPR/Cas9-UCA1 can effectively inhibit UCA1 expression when transfected into 5637 and T24 bladder cancer cells, while the combined transfection of the two most effective CRISPR/Cas9-UCA1s can generate more satisfied inhibitory effect. CRISPR/Cas9-UCA1s attenuate UCA1 expression via targeted genome-specific DNA cleavage, resulting in the significant inhibition of cell proliferation, migration and invasion in vitro and in vivo. The mechanisms associated with the inhibitory effect of CRISPR/Cas9-UCA1 on malignant phenotypes of bladder cancer are attributed to the induction of cell cycle arrest at G1 phase, a substantial increase of apoptosis, and an enhanced activity of MMPs. Additionally, urinary UCA1 can be used as a non-invasive diagnostic marker for bladder cancer as revealed by a meta-analysis. Collectively, our data suggest that CRISPR/Cas9 technique can be used to down-modulate lncRNA expression, and urinary UCA1 may be used as a non-invasive marker for diagnosis of bladder cancer.

Insights into the Function of Long Noncoding RNAs in Sepsis Revealed by Gene Co-Expression Network Analysis

Sepsis is a major cause of death and its incidence and mortality increase exponentially with age. Most gene expression studies in sepsis have focused in protein-coding genes and the expression patterns, and potential roles of long noncoding RNAs (lncRNAs) have not been investigated yet. In this study, we performed co-expression network analysis of protein-coding and lncRNAs measured in neutrophil granulocytes from adult and elderly septic patients, along with age-matched healthy controls. We found that the genes displaying highest network similarity are predominantly differently expressed in sepsis and are enriched in loci encoding proteins with structural or regulatory functions related to protein translation and mitochondrial energetic metabolism. A number of lncRNAs are strongly connected to genes from these pathways and may take part in regulatory loops that are perturbed in sepsis. Among those, the ribosomal pseudogenes RP11-302F12.1 and RPL13AP7 are differentially expressed and appear to have a regulatory role on protein translation in both the elderly and adults, and lncRNAs MALAT1, LINC00355, MYCNOS, and AC010970.2 display variable connection strength and inverted expression patterns between adult and elderly networks, suggesting that they are the best candidates to be further studied to understand the mechanisms by which the immune response is impaired by age. In summary, we report the expression of lncRNAs that are deregulated in patients with sepsis, including subsets that display hub properties in molecular pathways relevant to the disease pathogenesis and that may participate in gene expression regulatory circuits related to the poorer disease outcome observed in elderly subjects.

Long non-coding RNA Linc-RAM enhances myogenic differentiation by interacting with MyoD

Long non-coding RNAs (lncRNAs) are important regulators of diverse biological processes. Here we report on functional identification and characterization of a novel long intergenic non-coding RNA with MyoD-regulated and skeletal muscle-restricted expression that promotes the activation of the myogenic program, and is therefore termed Linc-RAM (Linc-RNA Activator of Myogenesis). Linc-RAM is transcribed from an intergenic region of myogenic cells and its expression is upregulated during myogenesis. Notably, in vivo functional studies show that Linc-RAM knockout mice display impaired muscle regeneration due to the differentiation defect of satellite cells. Mechanistically, Linc-RAM regulates expression of myogenic genes by directly binding MyoD, which in turn promotes the assembly of the MyoD–Baf60c–Brg1 complex on the regulatory elements of target genes. Collectively, our findings reveal the functional role and molecular mechanism of a lineage-specific Linc-RAM as a regulatory lncRNA required for tissues-specific chromatin remodelling and gene expression.

Long non-coding RNAs are important regulators of many diverse biological processes. Here the authors describe Linc-RAM, which regulates myogenesis by binding MyoD and promoting the assembly of the MyoD–Baf60c–Brg1 complex at target genes.

Distinctive Patterns of Transcription and RNA Processing for Human lincRNAs

Numerous long intervening noncoding RNAs (lincRNAs) are generated from the mammalian genome by RNA polymerase II (Pol II) transcription. Although multiple functions have been ascribed to lincRNAs, their synthesis and turnover remain poorly characterized. Here, we define systematic differences in transcription and RNA processing between protein-coding and lincRNA genes in human HeLa cells. This is based on a range of nascent transcriptomic approaches applied to different nuclear fractions, including mammalian native elongating transcript sequencing (mNET-seq). Notably, mNET-seq patterns specific for different Pol II CTD phosphorylation states reveal weak co-transcriptional splicing and poly(A) signal-independent Pol II termination of lincRNAs as compared to pre-mRNAs. In addition, lincRNAs are mostly restricted to chromatin, since they are rapidly degraded by the RNA exosome. We also show that a lincRNA-specific co-transcriptional RNA cleavage mechanism acts to induce premature termination. In effect, functional lincRNAs must escape from this targeted nuclear surveillance process.

Transcriptional profiling analysis and functional prediction of long noncoding RNAs in cancer

Long noncoding RNAs (lncRNAs), which are noncoding RNAs (ncRNAs) with length more than 200 nucleotides (nt), have been demonstrated to be involved in various types of cancer. Consequently, it has been frequently discussed that lncRNAs with aberrant expression in cancer serve as potential diagnostic biomarkers and therapeutic targets. However, one major challenge of developing cancer biomarkers is tumor heterogeneity which means that tumor cells show different cellular morphology, metastatic potential as well as gene expression. In this study, a custom designed microarray platform covering both mRNAs and lncRNAs was applied to tumor tissues of gastric, colon, liver and lung. 316 and 157 differentially expressed (DE-) protein coding genes and lncRNAs common to these four types of cancer were identified respectively. Besides, the functional roles of common DE-lncRNAs were inferred based on their expression and genomic position correlation with mRNAs. Moreover, mRNAs and lncRNAs with tissue specificity were also identified, suggesting their particular roles with regard to specific biogenesis and functions of different organs. Based on the large-scale survey of mRNAs and lncRNAs in four types of cancer, this study may offer new biomarkers common or specific for various types of cancer.

Identification and functional analysis of long intergenic noncoding RNA genes in porcine pre-implantation embryonic development

Genome-wide transcriptome studies have identified thousands of long intergenic noncoding RNAs (lincRNAs), some of which play important roles in pre-implantation embryonic development (PED). Pig is an ideal model for reproduction, however, porcine lincRNAs are still poorly characterized and it is unknown if they are associated with porcine PED. Here we reconstructed 195,531 transcripts in 122,007 loci, and identified 7,618 novel lincRNAs from 4,776 loci based on published RNA-seq data. These lincRNAs show low exon number, short length, low expression level, tissue-specific expression and cis-acting, which is consistent with previous reports in other species. By weighted co-expression network analysis, we identified 5 developmental stages specific co-expression modules. Gene ontology enrichment analysis of these specific co-expression modules suggested that many lincRNAs are associated with cell cycle regulation, transcription and metabolism to regulate the process of zygotic genome activation. Futhermore, we identified hub lincRNAs in each co-expression modules, and found two lincRNAs TCONS_00166370 and TCONS_00020255 may play a vital role in porcine PED. This study systematically analyze lincRNAs in pig and provides the first catalog of lincRNAs that might function as gene regulatory factors of porcine PED.

Comparative Transcriptomic and Epigenomic Analyses Reveal New Regulators of Murine Brown Adipogenesis

Increasing energy expenditure through brown adipocyte recruitment is a promising approach to combat obesity. We report here the comprehensive profiling of the epigenome and transcriptome throughout the lineage commitment and differentiation of C3H10T1/2 mesenchymal stem cell line into brown adipocytes. Through direct comparison to datasets from differentiating white adipocytes, we systematically identify stage- and lineage-specific coding genes, lncRNAs and microRNAs. Utilizing chromatin state maps, we also define stage- and lineage-specific enhancers, including super-enhancers, and their associated transcription factor binding motifs and genes. Through these analyses, we found that in brown adipocytes, brown lineage-specific genes are pre-marked by both H3K4me1 and H3K27me3, and the removal of H3K27me3 at the late stage is necessary but not sufficient to promote brown gene expression, while the pre-deposition of H3K4me1 plays an essential role in poising the brown genes for expression in mature brown cells. Moreover, we identify SOX13 as part of a p38 MAPK dependent transcriptional response mediating early brown cell lineage commitment. We also identify and subsequently validate PIM1, SIX1 and RREB1 as novel regulators promoting brown adipogenesis. Finally, we show that SIX1 binds to adipogenic and brown marker genes and interacts with C/EBPα, C/EBPβ and EBF2, suggesting their functional cooperation during adipogenesis.

Identification of Specific Long Non-Coding RNA Expression: Profile and Analysis of Association with Clinicopathologic Characteristics and BRAF Mutation in Papillary Thyroid Cancer

BACKGROUND

In recent years, long non-coding RNAs (lncRNAs) have been shown to play a critical regulatory role in cancer biology. However, the contribution of lncRNAs to papillary thyroid cancer (PTC) remains largely unknown.

METHODS

RNA sequencing and quantitative reverse transcription polymerase chain reaction were used to detect and verify changes to the transcriptome profile in 12 PTC tissues compared to paired normal adjacent tissues. The statistical correlation between differentially expressed lncRNAs and clinicopathologic characteristics was analyzed, and potential lncRNA functions were predicted by examining annotations for the co-expressed mRNAs. Furthermore, the specific subgroup patterns of the PTC transcriptome remodeled by BRAF mutations were also analyzed.

RESULTS

A total of 188 lncRNAs and 505 mRNAs were differentially expressed in 50% or more of the PTC tissues (fold change >2; p < 0.05) as assessed by RNA-sequencing compared with paired normal adjacent tissues. Forty-seven lncRNAs and 39 mRNAs were verified in 31 pairs of PTC specimens using quantitative reverse transcription polymerase chain reaction, and the results were consistent with the RNA sequencing data. The lncRNAs NONHSAT076747 and NONHSAT122730 were associated with lymph node metastasis, and NONHSAG051968 expression was negatively correlated with tumor size. A co-expression network between differentially expressed lncRNAs and protein-coding RNAs was constructed and analyzed, and functional analysis suggested that the differentially expressed genes mainly participate in ECM-receptor interactions and the focal adhesion pathway. Furthermore, a specific PTC transcriptome subtype pattern stratified by BRAF mutation was also uncovered. The p53 signaling pathway was the most highly enriched pathway among the BRAF mutation-related genes.

CONCLUSIONS

This study reveals specific changes to the lncRNA profile associated with PTC, and provides new insight into its pathogenesis. The PTC-associated lncRNAs NONHSAG051968, NONHSAT076747, and NONHSAT122730 might be potential diagnostic and therapeutic targets for PTC patients, and lncRNAs with subtype-specific expression stratified by BRAF mutation might be significant in individual molecular subtypes.

Long intergenic non-coding RNA 271 is predictive of a poorer prognosis of papillary thyroid cancer

A number of long non-coding RNAs (lncRNAs) have been found to play critical roles in oncogenesis and tumor progression. We aimed to investigate whether lncRNAs could act as prognostic biomarkers for papillary thyroid cancer (PTC) that may assist us in evaluating disease status and prognosis for patients. We found 220 lncRNAs with expression alteration from the annotated 2773 lncRNAs approved by the HUGO gene nomenclature committee in The Cancer Genome Atlas (TCGA) dataset, of which FAM41C, CTBP1-AS2, LINC00271, HAR1A, LINC00310 and HAS2-AS1 were associated with recurrence. After adjusting classical clinicopathogical factors and BRAF^V600E mutation, LINC00271 was found to be an independent risk factor for extrathyroidal extension, lymph node metastasis, advanced tumor stage III/IV and recurrence in multivariate analyses. Additionally, LINC00271 expression was significantly downregulated in PTCs versus adjacent normal tissues (P < 0.001). The Gene Set Enrichment Analysis (GSEA) revealed that genes associated with cell adhesion molecules, cell cycle, P53 signaling pathway and JAK/STAT signaling pathway were remarkably enriched in lower-LINC00271 versus higher-LINC00271 tumors. In conclusion, LINC00271 was identified as a possible suppressor gene in PTC in our study, and it may serve as a potential predictor of poor prognoses in PTC.

LincSNP 2.0: an updated database for linking disease-associated SNPs to human long non-coding RNAs and their TFBSs

We describe LincSNP 2.0 (http://bioinfo.hrbmu.edu.cn/LincSNP), an updated database that is used specifically to store and annotate disease-associated single nucleotide polymorphisms (SNPs) in human long non-coding RNAs (lncRNAs) and their transcription factor binding sites (TFBSs). In LincSNP 2.0, we have updated the database with more data and several new features, including (i) expanding disease-associated SNPs in human lncRNAs; (ii) identifying disease-associated SNPs in lncRNA TFBSs; (iii) updating LD-SNPs from the 1000 Genomes Project; and (iv) collecting more experimentally supported SNP-lncRNA-disease associations. Furthermore, we developed three flexible online tools to retrieve and analyze the data. Linc-Mart is a convenient way for users to customize their own data. Linc-Browse is a tool for all data visualization. Linc-Score predicts the associations between lncRNA and disease. In addition, we provided users a newly designed, user-friendly interface to search and download all the data in LincSNP 2.0 and we also provided an interface to submit novel data into the database. LincSNP 2.0 is a continually updated database and will serve as an important resource for investigating the functions and mechanisms of lncRNAs in human diseases.

Genome-wide analysis of alternative splicing during human heart development

Alternative splicing (AS) drives determinative changes during mouse heart development. Recent high-throughput technological advancements have facilitated genome-wide AS, while its analysis in human foetal heart transition to the adult stage has not been reported. Here, we present a high-resolution global analysis of AS transitions between human foetal and adult hearts. RNA-sequencing data showed extensive AS transitions occurred between human foetal and adult hearts, and AS events occurred more frequently in protein-coding genes than in long non-coding RNA (lncRNA). A significant difference of AS patterns was found between foetal and adult hearts. The predicted difference in AS events was further confirmed using quantitative reverse transcription-polymerase chain reaction analysis of human heart samples. Functional foetal-specific AS event analysis showed enrichment associated with cell proliferation-related pathways including cell cycle, whereas adult-specific AS events were associated with protein synthesis. Furthermore, 42.6% of foetal-specific AS events showed significant changes in gene expression levels between foetal and adult hearts. Genes exhibiting both foetal-specific AS and differential expression were highly enriched in cell cycle-associated functions. In conclusion, we provided a genome-wide profiling of AS transitions between foetal and adult hearts and proposed that AS transitions and deferential gene expression may play determinative roles in human heart development.

Predicting the Organelle Location of Noncoding RNAs Using Pseudo Nucleotide Compositions

Noncoding RNAs (ncRNAs) are implicated in various biological processes. Recent findings have demonstrated that the function of ncRNAs correlates with their provenance. Therefore, the recognition of ncRNAs from different organelle genomes will be helpful to understand their molecular functions. However, the weakness of experimental techniques limits the progress toward studying organellar ncRNAs and their functional relevance. As a complement of experiments, computational method provides an important choice to identify ncRNA in different organelles. Thus, a computational model was developed to identify ncRNAs from kinetoplast and mitochondrion organelle genomes. In this model, RNA sequences are encoded by "pseudo dinucleotide composition." It was observed by the jackknife test that the overall success rate achieved by the proposed model was 90.08 %. We hope that the proposed method will be helpful in predicting ncRNA organellar locations.

NONCODEv4: Annotation of Noncoding RNAs with Emphasis on Long Noncoding RNAs

The rapid development of high-throughput sequencing technologies and bioinformatics algorithms now enables detection and profiling of a large number of noncoding transcripts. Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides, are accumulating with important roles involved in biological processes and tissue physiology. In this chapter, we describe the use of NONCODEv4, a database that provide a comprehensive catalog of noncoding RNAs with particularly detailed annotations for lncRNAs. NONCODEv4 stores more than half million transcripts, of which more than 200,000 are lncRNAs. NONCODEv4 raises the concept of lncRNA genes and explores their expression and functions based on public transcriptome data. NONCODEv4 also integrated a series of online tools and have a web interface easy to use. NONCODEv4 is available at http://www.noncode.org/ http://www.bioinfo.org/ noncode.

lncScore: alignment-free identification of long noncoding RNA from assembled novel transcripts

RNA-Seq based transcriptome assembly has been widely used to identify novel lncRNAs. However, the best-performing transcript reconstruction methods merely identified 21% of full-length protein-coding transcripts from H. sapiens. Those partial-length protein-coding transcripts are more likely to be classified as lncRNAs due to their incomplete CDS, leading to higher false positive rate for lncRNA identification. Furthermore, potential sequencing or assembly error that gain or abolish stop codons also complicates ORF-based prediction of lncRNAs. Therefore, it remains a challenge to identify lncRNAs from the assembled transcripts, particularly the partial-length ones. Here, we present a novel alignment-free tool, lncScore, which uses a logistic regression model with 11 carefully selected features. Compared to other state-of-the-art alignment-free tools (e.g. CPAT, CNCI, and PLEK), lncScore outperforms them on accurately distinguishing lncRNAs from mRNAs, especially partial-length mRNAs in the human and mouse datasets. In addition, lncScore also performed well on transcripts from five other species (Zebrafish, Fly, C. elegans, Rat, and Sheep). To speed up the prediction, multithreading is implemented within lncScore, and it only took 2 minute to classify 64,756 transcripts and 54 seconds to train a new model with 21,000 transcripts with 12 threads, which is much faster than other tools. lncScore is available at https://github.com/WGLab/lncScore.

Single-cell RNA sequencing reveals dynamic changes in A-to-I RNA editome during early human embryogenesis

BACKGROUND

A-to-I RNA-editing mediated by ADAR (adenosine deaminase acting on RNA) enzymes that converts adenosine to inosine in RNA sequence can generate mutations and alter gene regulation in metazoans. Previous studies have shown that A-to-I RNA-editing plays vital roles in mouse embryogenesis. However, the RNA-editing activities in early human embryonic development have not been investigated.

RESULTS

Here, we characterized genome-wide A-to-I RNA-editing activities during human early embryogenesis by profiling 68 single cells from 29 human embryos spanning from oocyte to morula stages. We demonstrate dynamic changes in genome-wide RNA-editing during early human embryogenesis in a stage-specific fashion. In parallel with ADAR expression level changes, the genome-wide A-to-I RNA-editing levels in cells remained relatively stable until 4-cell stage, but dramatically decreased at 8-cell stage, continually decreased at morula stage. We detected 37 non-synonymously RNA-edited genes, of which 5 were frequently found in cells of multiple embryonic stages. Moreover, we found that A-to-I editings in miRNA-targeted regions of a substantial number of genes preferably occurred in one or two sequential stages.

CONCLUSIONS

Our single-cell analysis reveals dynamic changes in genome-wide RNA-editing during early human embryogenesis in a stage-specific fashion, and provides important insights into early human embryogenesis.

Regulation of PCGEM1 by p54/nrb in prostate cancer

PCGEM1 is a long non-coding RNA (lncRNA) that is often upregulated in prostate cancer. However, little is known how PCGEM1 is regulated. In the present study, we show transcriptional regulation of PCGEM1 in response to androgen deprivation by p54/nrb. While ectopic expression of p54/nrb increases, suppression of p54/nrb by RNAi or knockout (KO) reduces PCGEM1. Moreover, rescue experiments indicate that re-expression of p54/nrb in KO cells restores the ability to induce PCGEM1, leading to upregulation of the androgen receptor splice variant AR3 which has been shown to play a role in castration resistance. Finally, 3,3′-Diindolylmethane (DIM), a known chemoprevention agent, is capable of suppressing PCGEM1 expression by preventing the interaction of p54/nrb with the PCGEM1 promoter. In particular, DIM reduces tumor growth by suppression of PCGEM1 and promoting apoptosis in the castrated xenograft mouse model. Together, these results demonstrate a novel mechanism of p54/nrb-mediated expression of PCGEM1 and AR3, contributing to castration resistance in prostate cancer.

Long noncoding RNA (lincRNA), a new paradigm in gene expression control

Long intergenic non-coding RNAs (lincRNAs) are defined as RNA transcripts that are longer than 200 nucleotides. By definition, these RNAs must not have open reading frames that encode proteins. Many of these transcripts are encoded by RNA polymerase II, are spliced, and are poly-adenylated. This final fact indicates that there is a trove of information about lincRNAs in databases such as the Gene Expression Omnibus (GEO), which is a repository for RNAseq and microarray data. Recent experiments indicate that there are upwards of 15,000 lincRNAs encoded by the human genome. The term "intergenic" refers to the identification of these transcripts from regions of the genome that do not contain protein-encoding genes. These regions coincide with what was once labeled as the "junk DNA" portions of our genomes, which, upon careful examination by whole genome RNA sequencing experiments, clearly encode RNA transcripts. LincRNAs also contain promoter- or enhancer-associated RNAs that are gene proximal and can be either in the sense or antisense orientation, relative to the protein-coding gene with which they are associated. In this review, we describe the functions of lincRNAs playing roles in biological processes such as gene expression control, scaffold formation, and epigenetic control.

BmncRNAdb: a comprehensive database of non-coding RNAs in the silkworm, Bombyx mori

BACKGROUND

Long non-coding RNAs (lncRNAs) may play critical roles in a wide range of developmental processes of higher organisms. Recently, lncRNAs have been widely identified across eukaryotes and many databases of lncRNAs have been developed for human, mouse, fruit fly, etc. However, there is rare information about them in the only completely domesticated insect, silkworm (Bombyx mori).

DESCRIPTION

In this study, we systematically scanned lncRNAs using the available silkworm RNA-seq data and public unigenes. Finally, we identified and collected 6281 lncRNAs in the silkworm. Besides, we also collected 1986 microRNAs (miRNAs) from previous studies. Then, we organized them into a comprehensive and web-based database, BmncRNAdb. This database offers a user-friendly interface for data browse and online analysis as well as the three online tools for users to predict the target genes of lncRNA or miRNA.

CONCLUSIONS

We have systematically identified and collected the silkworm lncRNAs and constructed a comprehensive database of the silkworm lncRNAs and miRNAs. This work gives a glimpse into lncRNAs of the silkworm and lays foundations for the ncRNAs study of the silkworm and other insects in the future. The BmncRNAdb is freely available at http://gene.cqu.edu.cn/BmncRNAdb/index.php .

Satellite non-coding RNAs: the emerging players in cells, cellular pathways and cancer

For several decades, transcriptional inactivity was considered as one of the particular features of constitutive heterochromatin and, therefore, of its major component, satellite DNA sequences. However, more recently, succeeding evidences have demonstrated that these sequences can indeed be transcribed, yielding satellite non-coding RNAs with important roles in the organization and regulation of genomes. Since then, several studies have been conducted, trying to understand the function(s) of these sequences not only in the normal but also in cancer genomes. It is thought that the association between cancer and satncRNAs is mostly due to the influence of these transcripts in the genome instability, a hallmark of cancer. The few reports on satellite DNA transcription in cancer contexts point to its overexpression; however, this scenario may be far more complex, variable, and influenced by a number of factors and the exact role of satncRNAs in the oncogenic process remains poorly understood. The greater is the knowledge on the association of satncRNAs with cancer, the greater would be the opportunity to assist cancer treatment, either by the design of effective therapies targeting these molecules or by using them as biomarkers in cancer diagnosis, prognosis, and with predictive value.

ANGIOGENES: knowledge database for protein-coding and noncoding RNA genes in endothelial cells

Increasing evidence indicates the presence of long noncoding RNAs (lncRNAs) is specific to various cell types. Although lncRNAs are speculated to be more numerous than protein-coding genes, the annotations of lncRNAs remain primitive due to the lack of well-structured schemes for their identification and description. Here, we introduce a new knowledge database “ANGIOGENES” (http://angiogenes.uni-frankfurt.de) to allow for in silico screening of protein-coding genes and lncRNAs expressed in various types of endothelial cells, which are present in all tissues. Using the latest annotations of protein-coding genes and lncRNAs, publicly-available RNA-seq data was analyzed to identify transcripts that are expressed in endothelial cells of human, mouse and zebrafish. The analyzed data were incorporated into ANGIOGENES to provide a one-stop-shop for transcriptomics data to facilitate further biological validation. ANGIOGENES is an intuitive and easy-to-use database to allow in silico screening of expressed, enriched and/or specific endothelial transcripts under various conditions. We anticipate that ANGIOGENES serves as a starting point for functional studies to elucidate the roles of protein-coding genes and lncRNAs in angiogenesis.

A genetic algorithm-based weighted ensemble method for predicting transposon-derived piRNAs

BACKGROUND

Predicting piwi-interacting RNA (piRNA) is an important topic in the small non-coding RNAs, which provides clues for understanding the generation mechanism of gamete. To the best of our knowledge, several machine learning approaches have been proposed for the piRNA prediction, but there is still room for improvements.

RESULTS

In this paper, we develop a genetic algorithm-based weighted ensemble method for predicting transposon-derived piRNAs. We construct datasets for three species: Human, Mouse and Drosophila. For each species, we compile the balanced dataset and imbalanced dataset, and thus obtain six datasets to build and evaluate prediction models. In the computational experiments, the genetic algorithm-based weighted ensemble method achieves 10-fold cross validation AUC of 0.932, 0.937 and 0.995 on the balanced Human dataset, Mouse dataset and Drosophila dataset, respectively, and achieves AUC of 0.935, 0.939 and 0.996 on the imbalanced datasets of three species. Further, we use the prediction models trained on the Mouse dataset to identify piRNAs of other species, and the models demonstrate the good performances in the cross-species prediction.

CONCLUSIONS

Compared with other state-of-the-art methods, our method can lead to better performances. In conclusion, the proposed method is promising for the transposon-derived piRNA prediction. The source codes and datasets are available in https://github.com/zw9977129/piRNAPredictor .