Pseudogenes as an alternative source of natural antisense transcripts

LncRNAs and Chromatin Modifications Pattern m6A Methylation at the Untranslated Regions of mRNAs

New roles for RNA in mediating gene expression are being discovered at an alarming rate. A broad array of pathways control patterning of N⁶-methyladenosine (m⁶A) methylation on RNA transcripts. This review comprehensively discusses long non-coding RNAs (lncRNAs) as an additional dynamic regulator of m⁶A methylation, with a focus on the untranslated regions (UTRs) of mRNAs. Although there is extensive literature describing m⁶A modification of lncRNA, the function of lncRNA in guiding m⁶A writers has not been thoroughly explored. The independent control of lncRNA expression, its heterogeneous roles in RNA metabolism, and its interactions with epigenetic machinery, alludes to their potential in dynamic patterning of m⁶A methylation. While epigenetic regulation by histone modification of H3K36me3 has been demonstrated to pattern RNA m⁶A methylation, these modifications were specific to the coding and 3′UTR regions. However, there are observations that 5′UTR m⁶A is distinct from that of the coding and 3′UTR regions, and substantial evidence supports the active regulation of 5′UTR m⁶A methylation. Consequently, two potential mechanisms in patterning the UTRs m⁶A methylation are discussed; (1) Anti-sense lncRNA (AS-lncRNA) can either bind directly to the UTR, or (2) act indirectly via recruitment of chromatin-modifying complexes to pattern m⁶A. Both pathways can guide the m⁶A writer complex, facilitate m⁶A methylation and modulate protein translation. Findings in the lncRNA-histone-m⁶A axis could potentially contribute to the discovery of new functions of lncRNAs and clarify lncRNA-m⁶A findings in translational medicine.

PΨFinder: a practical tool for the identification and visualization of novel pseudogenes in DNA sequencing data

BACKGROUND

Processed pseudogenes (PΨgs) are disabled gene copies that are transcribed and may affect expression of paralogous genes. Moreover, their insertion in the genome can disrupt the structure or the regulatory region of a gene, affecting its expression level. These events have been identified as occurring mutations during cancer development, thus being able to identify PΨgs and their location will improve their impact on diagnostic testing, not only in cancer but also in inherited disorders.

RESULTS

We have implemented PΨFinder (P-psy-finder), a tool that identifies PΨgs, annotates known ones and predicts their insertion site(s) in the genome. The tool screens alignment files and provides user-friendly summary reports and visualizations. To demonstrate its applicability, we scanned 218 DNA samples from patients screened for hereditary colorectal cancer. We detected 423 PΨgs distributed in 96% of the samples, comprising 7 different parent genes. Among these, we confirmed the well-known insertion site of the SMAD4-PΨg within the last intron of the SCAI gene in one sample. While for the ubiquitous CBX3-PΨg, present in 82.6% of the samples, we found it reversed inserted in the second intron of the C15ORF57 gene.

CONCLUSIONS

PΨFinder is a tool that can automatically identify novel PΨgs from DNA sequencing data and determine their location in the genome with high sensitivity (95.92%). It generates high quality figures and tables that facilitate the interpretation of the results and can guide the experimental validation. PΨFinder is a complementary analysis to any mutational screening in the identification of disease-causing mutations within cancer and other diseases.

Metastatic EMT Phenotype Is Governed by MicroRNA-200-Mediated Competing Endogenous RNA Networks

Epithelial–mesenchymal transition (EMT) is a fundamental physiologically relevant process that occurs during morphogenesis and organ development. In a pathological setting, the transition from epithelial toward mesenchymal cell phenotype is hijacked by cancer cells, allowing uncontrolled metastatic dissemination. The competing endogenous RNA (ceRNA) hypothesis proposes a competitive environment resembling a large-scale regulatory network of gene expression circuits where alterations in the expression of both protein-coding and non-coding genes can make relevant contributions to EMT progression in cancer. The complex regulatory diversity is exerted through an array of diverse epigenetic factors, reaching beyond the transcriptional control that was previously thought to single-handedly govern metastatic dissemination. The present review aims to unravel the competitive relationships between naturally occurring ceRNA transcripts for the shared pool of the miRNA-200 family, which play a pivotal role in EMT related to cancer dissemination. Upon acquiring more knowledge and clinical evidence on non-genetic factors affecting neoplasia, modulation of the expression levels of diverse ceRNAs may allow for the development of novel prognostic/diagnostic markers and reveal potential targets for the disruption of cancer-related EMT.

Genome-wide DNA methylation and RNA expression differences correlate with invasiveness in melanoma cell lines

Aims & objectives: The aim of this study was to investigate the role of DNA methylation in invasiveness in melanoma cells. Materials & methods: The authors carried out genome-wide transcriptome (RNA sequencing) and reduced representation bisulfite sequencing methylome profiling between noninvasive (n = 4) and invasive melanoma cell lines (n = 5). Results: The integration of differentially expressed genes and differentially methylated fragments (DMFs) identified 12 DMFs (two in AVPI1, one in HMG20B, two in BCL3, one in NTSR1, one in SYNJ2, one in ROBO2 and four in HORMAD2) that overlapped with either differentially expressed genes (eight DMFs and six genes) or cis-targets of lncRNAs (five DMFs associated with cis-targets and four differentially expressed lncRNAs). Conclusions: DNA methylation changes are associated with a number of transcriptional differences observed in noninvasive and invasive phenotypes in melanoma.

DiseaseLinc: Disease Enrichment Analysis of Sets of Differentially Expressed LincRNAs

Long intergenic non-coding RNAs (LincRNAs) are long RNAs that do not encode proteins. Functional evidence is lacking for most of them. Their biogenesis is not well-known, but it is thought that many lincRNAs originate from genomic duplication of coding material, resulting in pseudogenes, gene copies that lose their original function and can accumulate mutations. While most pseudogenes eventually stop producing a transcript and become erased by mutations, many of these pseudogene-based lincRNAs keep similarity to the parental gene from which they originated, possibly for functional reasons. For example, they can act as decoys for miRNAs targeting the parental gene. Enrichment analysis of function is a powerful tool to discover the functional effects of a treatment producing differential expression of transcripts. However, in the case of lincRNAs, since their function is not easy to define experimentally, such a tool is lacking. To address this problem, we have developed an enrichment analysis tool that focuses on lincRNAs exploiting their functional association, using as a proxy function that of the parental genes and has a focus on human diseases.

Pseudogenes: A Novel Source of Trans-Acting Antisense RNAs

Several recent studies support a functional role for pseudogenes, a copy of a parent gene that has lost protein-coding potential, which was for a long time thought to represent only "junk" DNA. Several hundreds of pseudogenes have now been reported as transcribed RNAs in a large variety of tissues and tumor types. Most studies have focused on pseudogenes expressed in sense direction, relative to their protein-coding gene counterpart, but some reports suggest that pseudogenes can be also transcribed as antisense RNAs (asRNAs). Key regulatory genes, such as PTEN and OCT4, have in fact been reported to be under the regulation of pseudogene-expressed asRNAs. Here, we review what is known about pseudogene-expressed asRNAs, we discuss the functional role that these transcripts may have in gene regulation and we summarize the techniques that are available to study them.

Complex Analysis of Retroposed Genes' Contribution to Human Genome, Proteome and Transcriptome

Gene duplication is a major driver of organismal evolution. One of the main mechanisms of gene duplications is retroposition, a process in which mRNA is first transcribed into DNA and then reintegrated into the genome. Most gene retrocopies are depleted of the regulatory regions. Nevertheless, examples of functional retrogenes are rapidly increasing. These functions come from the gain of new spatio-temporal expression patterns, imposed by the content of the genomic sequence surrounding inserted cDNA and/or by selectively advantageous mutations, which may lead to the switch from protein coding to regulatory RNA. As recent studies have shown, these genes may lead to new protein domain formation through fusion with other genes, new regulatory RNAs or other regulatory elements. We utilized existing data from high-throughput technologies to create a complex description of retrogenes functionality. Our analysis led to the identification of human retroposed genes that substantially contributed to transcriptome and proteome. These retrocopies demonstrated the potential to encode proteins or short peptides, act as cis- and trans- Natural Antisense Transcripts (NATs), regulate their progenitors’ expression by competing for the same microRNAs, and provide a sequence to lncRNA and novel exons to existing protein-coding genes. Our study also revealed that retrocopies, similarly to retrotransposons, may act as recombination hot spots. To our best knowledge this is the first complex analysis of these functions of retrocopies.

LncRNA AC093818.1 accelerates gastric cancer metastasis by epigenetically promoting PDK1 expression

Gastric cancer (GC) is a highly prevalent type of metastatic tumor. The mechanisms underlying GC metastasis are poorly understood. Some long noncoding RNAs (lncRNAs) reportedly play key roles in regulating metastasis of GC. However, the biological roles of five natural antisense lncRNAs (AC093818.1, CTD-2541M15.1, BC047644, RP11-597M12.1, and RP11-40A13.1) in GC metastasis remain unclear. In this study, the expression of these lncRNAs was measured by quantitative reverse transcription-polymerase chain reaction. Migration and invasion were evaluated by wound-healing and the Transwell assay, respectively. Stable cells were injected into the tail veins of nude mice. Sections of collected lung and liver tissues were stained using hematoxylin and eosin. Protein expression was analyzed by western blot. RNA immunoprecipitation (RIP) assay was used to verify whether the STAT3 and SP1 transcription factors bound to AC093818.1 in GC cells. Expression levels of the five lncRNAs, especially AC093818.1, were significantly upregulated in metastatic GC tissues relative to those in nonmetastatic GC tissues. AC093818.1 expression was correlated with invasion, lymphatic metastasis, distal metastasis, and tumor-node-metastasis stage. AC093818.1 expression was highly sensitive and specific in the diagnosis of metastatic or nonmetastatic GC. AC093818.1 overexpression promoted GC migration and invasion in vitro and in vivo. AC093818.1 overexpression increased PDK1, p-AKT1, and p-mTOR expression levels. AC093818.1 silencing decreased these expressions. AC093818.1 bound to transcription factors STAT3 and SP1, and SP1 or STAT3 silencing could alleviated the effect of AC093818.1 overexpression. The data demonstrate that lncRNA AC093818.1 accelerates gastric cancer metastasis by epigenetically promoting PDK1 expression. LncRNA AC093818.1 may be a potential therapeutic target for metastatic GC.

Identification of transcribed protein coding sequence remnants within lincRNAs

Long intergenic non-coding RNAs (lincRNAs) are non-coding transcripts >200 nucleotides long that do not overlap protein-coding sequences. Importantly, such elements are known to be tissue-specifically expressed and to play a widespread role in gene regulation across thousands of genomic loci. However, very little is known of the mechanisms for the evolutionary biogenesis of these RNA elements, especially given their poor conservation across species. It has been proposed that lincRNAs might arise from pseudogenes. To test this systematically, we developed a novel method that searches for remnants of protein-coding sequences within lincRNA transcripts; the hypothesis is that we can trace back their biogenesis from protein-coding genes or posterior transposon/retrotransposon insertions. Applying this method, we found 203 human lincRNA genes with regions significantly similar to protein-coding sequences. Our method provides a visualization tool to trace the evolutionary biogenesis of lincRNAs with respect to protein-coding genes by sequence divergence. Subsequently, we show the expression correlation between lincRNAs and their identified parental protein-coding genes using public RNA-seq repositories, hinting at novel gene regulatory relationships. In summary, we developed a novel computational methodology to study non-coding gene sequences, which can be applied to identify the evolutionary biogenesis and function of lincRNAs.

ASSA: Fast identification of statistically significant interactions between long RNAs

The discovery of thousands of long noncoding RNAs (lncRNAs) in mammals raises a question about their functionality. It has been shown that some of them are involved in post-transcriptional regulation of other RNAs and form inter-molecular duplexes with their targets. Sequence alignment tools have been used for transcriptome-wide prediction of RNA-RNA interactions. However, such approaches have poor prediction accuracy since they ignore RNA's secondary structure. Application of the thermodynamics-based algorithms to long transcripts is not computationally feasible on a large scale. Here, we describe a new computational pipeline ASSA that combines sequence alignment and thermodynamics-based tools for efficient prediction of RNA-RNA interactions between long transcripts. To measure the hybridization strength, the sum energy of all the putative duplexes is computed. The main novelty implemented in ASSA is the ability to quickly estimate the statistical significance of the observed interaction energies. Most of the functional hybridizations between long RNAs were classified as statistically significant. ASSA outperformed 11 other tools in terms of the Area Under the Curve on two out of four test sets. Additionally, our results emphasized a unique property of the [Formula: see text] repeats with respect to the RNA-RNA interactions in the human transcriptome. ASSA is available at https://sourceforge.net/projects/assa/.

RetrogeneDB-a database of plant and animal retrocopies

For a long time, retrocopies were considered ‘junk DNA’, but numerous studies have shown that retrocopies may gain functionality and become so-called retrogenes. Retrogenes may code fully functional proteins that coexist with parental gene products or may even replace them. Retrocopies may also function as regulatory RNAs and, for example, become a source of small interfering RNAs, act as trans natural antisense transcripts or as alternative targets for miRNAs. Numerous researchers have emphasized that retrogenes play a crucial role in various organisms’ developmental stages and diseases. Despite the ever-growing evidence of the importance of retrocopies, resources dedicated to retroposition are very limited. Here, we report an update of the RetrogeneDB, which, to the best of our knowledge, is the largest database dedicated to retrocopies. It provides annotations of 86 458 retrocopies in 62 animal and 37 plant species. The database contains information about the retrocopies’ localization, open reading frame conservation, expression, RNA Polymerase II activity and the alternative transcription start site studies. Orthologous relationships between retrogenes were also determined, which made retrocopy conservation studies much more valuable. Additionally, based on the RNA-Seq data from the Geuvadis project, the expression levels of retrocopies were estimated in a total of 50 individuals from 5 human populations. The information is now presented in a new, more user-friendly web interface, with easy access to the source data, which may be used for the downstream analysis. RetrogeneDB is freely available at http://yeti.amu.edu.pl/retrogenedb.

Database URL:http://yeti.amu.edu.pl/retrogenedb

Secondary database URL:http://rhesus.amu.edu.pl/retrogenedb

Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

Protein-Coding Genes' Retrocopies and Their Functions

Transposable elements, often considered to be not important for survival, significantly contribute to the evolution of transcriptomes, promoters, and proteomes. Reverse transcriptase, encoded by some transposable elements, can be used in trans to produce a DNA copy of any RNA molecule in the cell. The retrotransposition of protein-coding genes requires the presence of reverse transcriptase, which could be delivered by either non-long terminal repeat (non-LTR) or LTR transposons. The majority of these copies are in a state of “relaxed” selection and remain “dormant” because they are lacking regulatory regions; however, many become functional. In the course of evolution, they may undergo subfunctionalization, neofunctionalization, or replace their progenitors. Functional retrocopies (retrogenes) can encode proteins, novel or similar to those encoded by their progenitors, can be used as alternative exons or create chimeric transcripts, and can also be involved in transcriptional interference and participate in the epigenetic regulation of parental gene expression. They can also act in trans as natural antisense transcripts, microRNA (miRNA) sponges, or a source of various small RNAs. Moreover, many retrocopies of protein-coding genes are linked to human diseases, especially various types of cancer.

The distribution and evolution of Arabidopsis thaliana cis natural antisense transcripts

Background

Natural antisense transcripts (NATs) are regulatory RNAs that contain sequence complementary to other RNAs, these other RNAs usually being messenger RNAs. In eukaryotic genomes, cis-NATs overlap the gene they complement.

Results

Here, our goal is to analyze the distribution and evolutionary conservation of cis-NATs for a variety of available data sets for Arabidopsis thaliana, to gain insights into cis-NAT functional mechanisms and their significance. Cis-NATs derived from traditional sequencing are largely validated by other data sets, although different cis-NAT data sets have different prevalent cis-NAT topologies with respect to overlapping protein-coding genes. A. thaliana cis-NATs have substantial conservation (28-35% in the three substantive data sets analyzed) of expression in A. lyrata. We examined evolutionary sequence conservation at cis-NAT loci in Arabidopsis thaliana across nine sequenced Brassicaceae species (picked for optimal discernment of purifying selection), focussing on the parts of their sequences not overlapping protein-coding transcripts (dubbed ‘NOLPs’). We found significant NOLP sequence conservation for 28-34% NATs across different cis-NAT sets. This NAT NOLP sequence conservation versus A. lyrata is generally significantly correlated with conservation of expression. We discover a significant enrichment of transcription factor binding sites (as evidenced by CHIP-seq data) in NOLPs compared to randomly sampled near-gene NOLP-like DNA , that is linked to significant sequence conservation. Conversely, there is no such evidence for a general significant link between NOLPs and formation of small interfering RNAs (siRNAs), with the substantial majority of unique siRNAs arising from the overlapping portions of the cis-NATs.

Conclusions

In aggregate, our results suggest that many cis-NAT NOLPs function in the regulation of conserved promoter/regulatory elements that they ‘over-hang’.

Electronic supplementary material

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

Characterization of human pseudogene-derived non-coding RNAs for functional potential

Thousands of pseudogenes exist in the human genome and many are transcribed, but their functional potential remains elusive and understudied. To explore these issues systematically, we first developed a computational pipeline to identify transcribed pseudogenes from RNA-Seq data. Applying the pipeline to datasets from 16 distinct normal human tissues identified ∼ 3,000 pseudogenes that could produce non-coding RNAs in a manner of low abundance but high tissue specificity under normal physiological conditions. Cross-tissue comparison revealed that the transcriptional profiles of pseudogenes and their parent genes showed mostly positive correlations, suggesting that pseudogene transcription could have a positive effect on the expression of their parent genes, perhaps by functioning as competing endogenous RNAs (ceRNAs), as previously suggested and demonstrated with the PTEN pseudogene, PTENP1. Our analysis of the ENCODE project data also found many transcriptionally active pseudogenes in the GM12878 and K562 cell lines; moreover, it showed that many human pseudogenes produced small RNAs (sRNAs) and some pseudogene-derived sRNAs, especially those from antisense strands, exhibited evidence of interfering with gene expression. Further integrated analysis of transcriptomics and epigenomics data, however, demonstrated that trimethylation of histone 3 at lysine 9 (H3K9me3), a posttranslational modification typically associated with gene repression and heterochromatin, was enriched at many transcribed pseudogenes in a transcription-level dependent manner in the two cell lines. The H3K9me3 enrichment was more prominent in pseudogenes that produced sRNAs at pseudogene loci and their adjacent regions, an observation further supported by the co-enrichment of SETDB1 (a H3K9 methyltransferase), suggesting that pseudogene sRNAs may have a role in regional chromatin repression. Taken together, our comprehensive and systematic characterization of pseudogene transcription uncovers a complex picture of how pseudogene ncRNAs could influence gene and pseudogene expression, at both epigenetic and post-transcriptional levels.

Pseudogenes and their composers: delving in the 'debris' of human genome

Pseudogenes, the nonfunctional homologs of functional genes and thus exemplified as 'genomic fossils' provide intriguing snapshots of the evolutionary history of human genome. These defunct copies generally arise by retrotransposition or duplication followed by various genetic disablements. In this study, focusing on human pseudogenes and their functional homologues we describe their characteristic features and relevance to protein sequence evolution. We recapitulate that pseudogenes harbor disease-causing degenerative sequence variations in conjunction with the immense disease gene association of their progenitors. Furthermore, we also discuss the issue of functional resurrection and the potentiality observed in some pseudogenes to regulate their functional counterparts.

Natural antisense transcripts in plants: a review and identification in soybean infected with Phakopsora pachyrhizi SuperSAGE library

Natural antisense ranscripts (NAT) are RNA molecules complementary to other endogenous RNAs. They are capable of regulating the expression of target genes at different levels (transcription, mRNA stability, translation, etc.). Such a property makes them ideal for interventions in organisms' metabolism. The present study reviewed plant NAT aspects, including features, availability and genesis, conservation and distribution, coding capacity, NAT pair expression, and functions. Besides, an in silico identification of NATs pairs was presented, using deepSuperSAGE libraries of soybean infected or not with Phakopsora pachyrhizi. Results showed that around 1/3 of the 77,903 predicted trans-NATs (by PlantsNATsDB database) detected had unitags mapped in both sequences of each pair. The same 1/3 of the 436 foreseen cis-NATs showed unitags anchored in both sequences of the related pairs. For those unitags mapped in NAT pairs, a modulation expression was assigned as upregulated, downregulated, or constitutive, based on the statistical analysis (P < 0.05). As a result, the infected treatment promoted the expression of 2,313 trans-NATs pairs comprising unitags exclusively from that library (1,326 pairs had unitags only found in the mock library). To understand the regulation of these NAT pairs could be a key aspect in the ASR plant response.

Long non-coding RNA in cancer

Long non-coding RNAs (lncRNAs) are pervasively transcribed in the genome and are emerging as new players in tumorigenesis due to their various functions in transcriptional, posttranscriptional and epigenetic mechanisms of gene regulation. LncRNAs are deregulated in a number of cancers, demonstrating both oncogenic and tumor suppressive roles, thus suggesting their aberrant expression may be a substantial contributor in cancer development. In this review, we will summarize their emerging role in human cancer and discuss their perspectives in diagnostics as potential biomarkers.

Regulatory mechanisms of long noncoding RNAs in vertebrate central nervous system development and function

Long noncoding RNAs (lncRNAs) have emerged as an important class of molecules that regulate gene expression at epigenetic, transcriptional, and post-transcriptional levels through a wide array of mechanisms. This regulation is of particular importance in the central nervous system (CNS), where precise modulation of gene expression is required for proper neuronal and glial production, connection and function. There are relatively few functional studies that characterize lncRNA mechanisms, but possible functions can often be inferred based on existing examples and the lncRNA's relative genomic position. In this review, we will discuss mechanisms of lncRNAs as predicted by genomic contexts and the possible impact on CNS development, function, and disease pathogenesis. There is no doubt that investigation of the mechanistic role of lncRNAs will open a new and exciting direction in studying CNS development and function.

The complex transcriptional landscape of the anucleate human platelet

Background

Human blood platelets are essential to maintaining normal hemostasis, and platelet dysfunction often causes bleeding or thrombosis. Estimates of genome-wide platelet RNA expression using microarrays have provided insights to the platelet transcriptome but were limited by the number of known transcripts. The goal of this effort was to deep-sequence RNA from leukocyte-depleted platelets to capture the complex profile of all expressed transcripts.

Results

From each of four healthy individuals we generated long RNA (≥40 nucleotides) profiles from total and ribosomal-RNA depleted RNA preparations, as well as short RNA (<40 nucleotides) profiles. Analysis of ~1 billion reads revealed that coding and non-coding platelet transcripts span a very wide dynamic range (≥16 PCR cycles beyond β-actin), a result we validated through qRT-PCR on many dozens of platelet messenger RNAs. Surprisingly, ribosomal-RNA depletion significantly and adversely affected estimates of the relative abundance of transcripts. Of the known protein-coding loci, ~9,500 are present in human platelets. We observed a strong correlation between mRNAs identified by RNA-seq and microarray for well-expressed mRNAs, but RNASeq identified many more transcripts of lower abundance and permitted discovery of novel transcripts.

Conclusions

Our analyses revealed diverse classes of non-coding RNAs, including: pervasive antisense transcripts to protein-coding loci; numerous, previously unreported and abundant microRNAs; retrotransposons; and thousands of novel un-annotated long and short intronic transcripts, an intriguing finding considering the anucleate nature of platelets. The data are available through a local mirror of the UCSC genome browser and can be accessed at: http://cm.jefferson.edu/platelets_2012/.

Pseudogenes: newly discovered players in human cancer

Because they are generally noncoding and thus considered nonfunctional and unimportant, pseudogenes have long been neglected. Recent advances have established that the DNA of a pseudogene, the RNA transcribed from a pseudogene, or the protein translated from a pseudogene can have multiple, diverse functions and that these functions can affect not only their parental genes but also unrelated genes. Therefore, pseudogenes have emerged as a previously unappreciated class of sophisticated modulators of gene expression, with a multifaceted involvement in the pathogenesis of human cancer.

Computational approaches to discovering noncoding RNA

New developments are being brought to the field of molecular biology with the mounting evidence that RNA transcripts not translated into protein (noncoding RNAs, ncRNAs) hold a variety of biological functions. Computational discovery of ncRNAs is one of these developments, fueled not only by the urge to characterize these sequences but also by necessity to prioritize ones with the most relevant functions for experimental verification. The heterogeneity in size and mode of activity of ncRNAs is reflected in the corresponding diversity of computational methods for their study. Sequence and structural analysis, conservation across species, and relative position to other genomic elements are being used for ncRNA detection. In addition, the recent development of techniques that allow deep sequencing of cell transcripts either globally or from isolated ncRNA-related material is leading the field toward increased use of such high-throughput data. We expect that imminent breakthroughs will include the classification of newer types of ncRNA and new insights into miRNA and piRNA biology, eventually leading toward the completion of a catalog of all human ncRNAs.

Long non-coding RNAs: versatile master regulators of gene expression and crucial players in cancer

With rapid development of sequencing technologies such as deep sequencing and whole genome high-density tiling array, we now know that most of the "junk" genomic sequences are transcribed as non-coding RNAs (ncRNAs). A large number of long ncRNA transcripts (> 200bp) have been identified, and these long ncRNAs (LncRNAs) are found to be crucial regulators for epigenetic modulation, transcription, and translation. In this review, we briefly summarize the regulatory function of LncRNAs with a particular focus on the underlying mechanisms of LncRNAs in oncogenesis, tumor metastasis and suppression.

Septin genomics: a road less travelled

The human septins are part of a gene family, that is a group of genes with similar sequences and usually but not invariably share similar functions that are descended from a common ancestor. Here we review our current knowledge of the human septin gene family and highlight areas of uncertainty. Currently 13 human septin genes are known (SEPT1 to SEPT12 and SEPT14). What was known as SEPT13 is now defined as one of many SEPT7 related pseudogenes. The family is characterized by complex genomics and extensive (but not universal) splicing, giving rise to a plethora of septin isoforms. For only a few members of the family do we have a comprehensive insight into these transcripts and isoforms. Given the formation of countless septin homotypic and heterotypic interactions our understanding of the biology and pathobiology of the septin family will require a detailed understanding of the genomics, transcriptomics and regulation of all members of this diverse and complex family.

Functional evidence of post-transcriptional regulation by pseudogenes

Pseudogenes have been mainly considered as functionless evolutionary relics since their discovery in 1977. However, multiple mechanisms of pseudogene functionality have been proposed both at the transcriptional and post-transcriptional level. This review focuses on the role of pseudogenes as post-transcriptional regulators. Two lines of research have recently presented strong evidence of their potential function as post-transcriptional regulators of the corresponding parental genes from which they originate. First, pseudogene genomic sequences can encode siRNAs. Second, pseudogene transcripts can act as indirect post-transcriptional regulators decoying ncRNA, in particular miRNAs that target the parental gene. This has been demonstrated for PTEN and KRAS, two genes involved in tumorigenesis. The role of pseudogenes in disease has not been proven and seems to be the next research landmark. In this review, we chronicle the events following the initial discovery of the 'useless' pseudogene to its breakthrough as a functional molecule with hitherto unbeknownst potential to influence human disease.

Global expression analysis revealed novel gender-specific gene expression features in the blood fluke parasite Schistosoma japonicum

Background

Schistosoma japonicum is one of the remarkable Platyhelminths that are endemic in China and Southeast Asian countries. The parasite is dioecious and can reside inside the host for many years. Rapid reproduction by producing large number of eggs and count-react host anti-parasite responses are the strategies that benefit long term survival of the parasite. Praziquantel is currently the only drug that is effective against the worms. Development of novel antiparasite reagents and immune-prevention measures rely on the deciphering of parasite biology. The decoding of the genomic sequence of the parasite has made it possible to dissect the functions of genes that govern the development of the parasite. In this study, the polyadenylated transcripts from male and female S. japonicum were isolated for deep sequencing and the sequences were systematically analysed.

Results

First, the number of genes actively expressed in the two sexes of S. japonicum was similar, but around 50% of genes were biased to either male or female in expression. Secondly, it was, at the first time, found that more than 50% of the coding region of the genome was transcribed from both strands. Among them, 65% of the genes had sense and their cognate antisense transcripts co-expressed, whereas 35% had inverse relationship between sense and antisense transcript abundance. Further, based on gene ontological analysis, more than 2,000 genes were functionally categorized and biological pathways that are differentially functional in male or female parasites were elucidated.

Conclusions

Male and female schistosomal parasites differ in gene expression patterns, many metabolic and biological pathways have been identified in this study and genes differentially expressed in gender specific manner were presented. Importantly, more than 50% of the coding regions of the S. japonicum genome transcribed from both strands, antisense RNA-mediated gene regulation might play a critical role in the parasite biology.