Pseudogenes: pseudo-functional or key regulators in health and disease?

The persistent contributions of RNA to eukaryotic gen(om)e architecture and cellular function

Currently, the best scenario for earliest forms of life is based on RNA molecules as they have the proven ability to catalyze enzymatic reactions and harbor genetic information. Evolutionary principles valid today become apparent in such models already. Furthermore, many features of eukaryotic genome architecture might have their origins in an RNA or RNA/protein (RNP) world, including the onset of a further transition, when DNA replaced RNA as the genetic bookkeeper of the cell. Chromosome maintenance, splicing, and regulatory function via RNA may be deeply rooted in the RNA/RNP worlds. Mostly in eukaryotes, conversion from RNA to DNA is still ongoing, which greatly impacts the plasticity of extant genomes. Raw material for novel genes encoding protein or RNA, or parts of genes including regulatory elements that selection can act on, continues to enter the evolutionary lottery.

Analysis of stranded information using an automated procedure for strand specific RNA sequencing

Background

Strand specific RNA sequencing is rapidly replacing conventional cDNA sequencing as an approach for assessing information about the transcriptome. Alongside improved laboratory protocols the development of bioinformatical tools is steadily progressing. In the current procedure the Illumina TruSeq library preparation kit is used, along with additional reagents, to make stranded libraries in an automated fashion which are then sequenced on Illumina HiSeq 2000. By the use of freely available bioinformatical tools we show, through quality metrics, that the protocol is robust and reproducible. We further highlight the practicality of strand specific libraries by comparing expression of strand specific libraries to non-stranded libraries, by looking at known antisense transcription of pseudogenes and by identifying novel transcription. Furthermore, two ribosomal depletion kits, RiboMinus and RiboZero, are compared and two sequence aligners, Tophat2 and STAR, are also compared.

Results

The, non-stranded, Illumina TruSeq kit can be adapted to generate strand specific libraries and can be used to access detailed information on the transcriptome. The RiboZero kit is very effective in removing ribosomal RNA from total RNA and the STAR aligner produces high mapping yield in a short time. Strand specific data gives more detailed and correct results than does non-stranded data as we show when estimating expression values and in assembling transcripts. Even well annotated genomes need improvements and corrections which can be achieved using strand specific data.

Conclusions

Researchers in the field should strive to use strand specific data; it allows for more confidence in the data analysis and is less likely to lead to false conclusions. If faced with analysing non-stranded data, researchers should be well aware of the caveats of that approach.

Electronic supplementary material

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

Introduction to microRNAs in biological systems

MicroRNAs are 20-24-nucleotide-long noncoding RNAs that bind to the 3' UTR (untranslated region) of target mRNAs. Since their discovery, microRNAs have been gaining attention for their ability to contribute to gene expression regulation under various physiological conditions. Consequently, deregulated expression of microRNAs has been linked to different disease states. Here, a brief overview of the canonical and alternative microRNA biogenesis pathways and microRNA functions in biological systems is given based on recent developments. In addition, newly emerging regulatory mechanisms, such as alternative polyadenylation, in connection with microRNA-dependent gene expression regulation are discussed.

Relationship between gene duplicability and diversifiability in the topology of biochemical networks

Background

Selective gene duplicability, the extensive expansion of a small number of gene families, is universal. Quantitatively, the number of genes (P_(K)) with K duplicates in a genome decreases precipitously as K increases, and often follows a power law (P_(k)∝k^-α). Functional diversification, either neo- or sub-functionalization, is a major evolution route for duplicate genes.

Results

Using three lines of genomic datasets, we studied the relationship between gene duplicability and diversifiability in the topology of biochemical networks. First, we explored scenario where two pathways in the biochemical networks antagonize each other. Synthetic knockout of respective genes for the two pathways rescues the phenotypic defects of each individual knockout. We identified duplicate gene pairs with sufficient divergences that represent this antagonism relationship in the yeast S. cerevisiae. Such pairs overwhelmingly belong to large gene families, thus tend to have high duplicability. Second, we used distances between proteins of duplicate genes in the protein interaction network as a metric of their diversification. The higher a gene’s duplicate count, the further the proteins of this gene and its duplicates drift away from one another in the networks, which is especially true for genetically antagonizing duplicate genes. Third, we computed a sequence-homology-based clustering coefficient to quantify sequence diversifiability among duplicate genes – the lower the coefficient, the more the sequences have diverged. Duplicate count (K) of a gene is negatively correlated to the clustering coefficient of its duplicates, suggesting that gene duplicability is related to the extent of sequence divergence within the duplicate gene family.

Conclusion

Thus, a positive correlation exists between gene diversifiability and duplicability in the context of biochemical networks – an improvement of our understanding of gene duplicability.

The Pan-Cancer analysis of pseudogene expression reveals biologically and clinically relevant tumour subtypes

Although individual pseudogenes have been implicated in tumor biology, the biomedical significance and clinical relevance of pseudogene expression have not been assessed in a systematic way. Here we generate pseudogene expression profiles in 2,808 patient samples of seven cancer types from The Cancer Genome Atlas RNA-seq data using a newly developed computational pipeline. Supervised analysis reveals a significant number of pseudogenes differentially expressed among established tumor subtypes; and pseudogene expression alone can accurately classify the major histological subtypes of endometrial cancer. Across cancer types, the tumor subtypes revealed by pseudogene expression show extensive and strong concordance with the subtypes defined by other molecular data. Strikingly, in kidney cancer, the pseudogene-expression subtypes not only significantly correlate with patient survival, but also help stratify patients in combination with clinical variables. Our study highlights the potential of pseudogene expression analysis as a new paradigm for investigating cancer mechanisms and discovering prognostic biomarkers.

piRNAQuest: searching the piRNAome for silencers

Background

PIWI-interacting RNA (piRNA) is a novel and emerging class of small non-coding RNA (sncRNA). Ranging in length from 26-32 nucleotides, this sncRNA is a potent player in guiding the vital regulatory processes within a cellular system. Inspite of having such a wide role within cellular systems, piRNAs are not well organized and classified, so that a researcher can pool out the biologically relevant information concerning this class.

Description

Here we present piRNAQuest- a unified and comprehensive database of 41749 human, 890078 mouse and 66758 rat piRNAs obtained from NCBI and different small RNA sequence experiments. This database provides piRNA annotation based on their localization in gene, intron, intergenic, CDS, 5^/UTR, 3^/UTR and repetitive regions which has not been done so far. We have also annotated piRNA clusters and have elucidated characteristic motifs within them. We have looked for the presence of piRNAs and piRNA clusters in pseudogenes, which are known to regulate the expression of protein coding transcripts by generating small RNAs. All these will help researchers progress towards solving the unanswered queries on piRNA biogenesis and their mode of action. Further, expression profile for piRNA in different tissues and from different developmental stages has been provided. In addition, we have provided several tools like 'homology search’, 'dynamic cluster search’ and 'pattern search’. Overall, piRNAQuest will serve as a useful resource for exploring human, mouse and rat piRNAome. The database is freely accessible and available at http://bicresources.jcbose.ac.in/zhumur/pirnaquest/.

Conclusion

piRNAs play a remarkable role in stem cell self-renewal and various vital processes of developmental biology. Although researchers are mining different features on piRNAs, the exact regulatory mechanism is still fuzzy. Thus, understanding the true potential of these small regulatory molecules with respect to their origin, localization and mode of biogenesis is crucial. piRNAQuest will provide us with a better insight on piRNA origin and function which will help to explore the true potential of these sncRNAs.

Evaluation of coding-independent functions of the transcribed bovine aromatase pseudogene CYP19P1

Background

CYP19A1 encodes the aromatase which catalyzes the final reaction of estrogen biosynthesis. The bovine genome also contains a non-coding copy of CYP19A1, the transcribed pseudogene CYP19P1. Whereas CYP19A1 is transcribed in all estrogen-producing tissues, mainly in the placenta and gonads, the CYP19P1 transcript so far was detected in the placenta. Strikingly, one sequence segment of both transcripts exhibits an exceptional high identity of 98%, which implies selective pressure and suggests some kind of function. Only recently, indeed, coding-independent functions of several transcribed pseudogenes were reported. Therefore, we analyzed CYP19P1 and CYP19A1 transcripts with the aim to detect clues for gene–pseudogene interference.

Findings

The CYP19P1 transcript was first examined in silico for the presence of microRNA coding sequences and microRNA targets. Further, to identify tissues where CYP19P1 and CYP19A1 transcripts are co-expressed, as a pre-requisite for transcript interference, expression profiling was performed in a variety of bovine tissues. Our in silico analyses did neither reveal potential microRNA coding sequences, nor microRNA targets. Co-expression of the CYP19 loci was demonstrated in placental cotyledons and granulosa cells of dominant follicles. However, in granulosa cells of dominant follicles the concentration of CYP19P1 mRNA was very low compared to CYP19A1 mRNA.

Conclusions

CYP19P1 and CYP19A1 transcripts might interfere in placental cotyledons. However, in granulosa cells of dominant follicles relevant interference between gene and pseudogene transcripts is unlikely to occur because of the very low CYP19P1/CYP19A1 transcript ratio.

Methods to study the occurrence and the evolution of pseudogenes through a phylogenetic approach

During the last few years, the study of pseudogenes has excited enthusiasm, because it has been proven that at least some of them are involved in important biological processes. An accurate detection and analysis of pseudogenes can be achieved using comparative methods, but only the use of phylogenetic tools can provide accurate information about their birth, their evolution and their death, hence about the impact that they have on genes and genomes. Here, phylogenetic methods that allow studying pseudogene history are described.

Methods for detecting transcribed pseudogenes: PCR on regions of high sequence similarity followed by cloning and sequencing

Identifying pseudogene transcription is problematic in many cases due to the high sequence similarity between pseudogenes and their parental genes. In this chapter, we detail the procedure for the detection of pseudogene transcription using the reverse transcription polymerase chain reaction (RT-PCR) method. The protocol comprises (1) extraction of total RNA, (2) first-strand cDNA synthesis from total RNA, (3) amplification of the cDNA by PCR, and (4) cloning and sequencing of the PCR products. Technical and practical guidance is provided, and the critical points during each of the steps are discussed. In particular, the importance of designing high specific PCR primers and thoroughly eliminating genomic DNA contamination from RNA preparation is emphasized.

The rise of regulatory RNA

Discoveries over the past decade portend a paradigm shift in molecular biology. Evidence suggests that RNA is not only functional as a messenger between DNA and protein but also involved in the regulation of genome organization and gene expression, which is increasingly elaborate in complex organisms. Regulatory RNA seems to operate at many levels; in particular, it plays an important part in the epigenetic processes that control differentiation and development. These discoveries suggest a central role for RNA in human evolution and ontogeny. Here, we review the emergence of the previously unsuspected world of regulatory RNA from a historical perspective.

Pseudogenes as regulators of biological function

A pseudogene arises when a gene loses the ability to produce a protein, which can be due to mutation or inaccurate duplication. Previous dogma has dictated that because the pseudogene no longer produces a protein it becomes functionless and evolutionarily inert, being neither conserved nor removed. However, recent evidence has forced a re-evaluation of this view. Some pseudogenes, although not translated into protein, are at least transcribed into RNA. In some cases, these pseudogene transcripts are capable of influencing the activity of other genes that code for proteins, thereby altering expression and in turn affecting the phenotype of the organism. In the present chapter, we will define pseudogenes, describe the evidence that they are transcribed into non-coding RNAs and outline the mechanisms by which they are able to influence the machinery of the eukaryotic cell.

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.

The fate of tandemly duplicated genes assessed by the expression analysis of a group of Arabidopsis thaliana RING-H2 ubiquitin ligase genes of the ATL family

Gene duplication events exert key functions on gene innovations during the evolution of the eukaryotic genomes. A large portion of the total gene content in plants arose from tandem duplications events, which often result in paralog genes with high sequence identity. Ubiquitin ligases or E3 enzymes are components of the ubiquitin proteasome system that function during the transfer of the ubiquitin molecule to the substrate. In plants, several E3s have expanded in their genomes as multigene families. To gain insight into the consequences of gene duplications on the expansion and diversification of E3s, we examined the evolutionary basis of a cluster of six genes, duplC-ATLs, which arose from segmental and tandem duplication events in Brassicaceae. The assessment of the expression suggested two patterns that are supported by lineage. While retention of expression domains was observed, an apparent absence or reduction of expression was also inferred. We found that two duplC-ATL genes underwent pseudogenization and that, in one case, gene expression is probably regained. Our findings provide insights into the evolution of gene families in plants, defining key events on the expansion of the Arabidopsis Tóxicos en Levadura family of E3 ligases.

Evolution and biology of supernumerary B chromosomes

B chromosomes (Bs) are dispensable components of the genome exhibiting non-Mendelian inheritance and have been widely reported on over several thousand eukaryotes, but still remain an evolutionary mystery ever since their first discovery over a century ago [1]. Recent advances in genome analysis have significantly improved our knowledge on the origin and composition of Bs in the last few years. In contrast to the prevalent view that Bs do not harbor genes, recent analysis revealed that Bs of sequenced species are rich in gene-derived sequences. We summarize the latest findings on supernumerary chromosomes with a special focus on the origin, DNA composition, and the non-Mendelian accumulation mechanism of Bs.

Transcriptional profiling reveals differential expression of a neuropeptide-like protein and pseudogenes in aryl hydrocarbon receptor-1 mutant Caenorhabditis elegans

The aryl hydrocarbon receptor (AHR) functions in higher organisms in development, metabolism and toxic responses. Its Caenorhabditis elegans (C. elegans) ortholog, AHR-1, facilitates neuronal development, growth and movement. We investigated the effect of AHR mutation on the transcriptional profile of L4 stage C. elegans using RNA-seq and quantitative real time PCR in order to understand better AHR-1 function at the genomic level. Illumina HiSeq 2000 sequencing yielded 51.1, 61.2 and 54.0 million reads from wild-type controls, ahr-1(ia03) and ahr-1(ju145) mutants, respectively, providing detection of over 18,000 transcripts in each sample. Fourteen transcripts were over-expressed and 125 under-expressed in both ahr-1 mutants when compared to wild-type. Under-expressed genes included soluble guanylate cyclase (gcy) family genes, some of which were previously demonstrated to be regulated by AHR-1. A neuropeptide-like protein gene, nlp-20, and a F-box domain protein gene fbxa-192 and its pseudogenes fbxa-191 and fbxa-193 were also under-expressed. Conserved xenobiotic response elements were identified in the 5' flanking regions of some but not all of the gcy, nlp-20, and fbxa genes. These results extend previous studies demonstrating control of gcy family gene expression by AHR-1, and furthermore suggest a role of AHR-1 in regulation of a neuropeptide gene as well as pseudogenes.

Contribution of pseudogenes to sequence diversity

Pseudogenes are very common in the genomes of a wide range of organisms and, although they were originally considered as genetic junk, now several functions have been attributed to them. One important function of pseudogenes, as discussed in this chapter, is to provide material for genetic diversity. This is most prominent in the case of immunological recognition molecules such as immunoglobulins and B- and T-cell receptors, as well as in the case of antigenic variation in intracellular pathogens. Other examples discussed are olfactory receptors, ribosomal proteins, cytochrome P450s, and pseudokinases.

Pseudogene redux with new biological significance

The study of pseudogenes, originally dismissed as genomic relics of evolutionary selection, has seen a resurgence in scientific literature, in addition to being a peculiar topic of discussion in theological debates. For a long time, pseudogenes have been touted as a beacon of natural selection and a definitive proof of evolution due to the slow mutation rate that differentiated them from their parental genes and ultimately caused their genetic demise as functional genes. It now seems that "creationists" have co-opted some recent reports identifying unheralded biological functions to pseudogens and other noncoding RNAs as evidence to undermine the existence of evolution and supporting intelligent design. This issue of Methods in Molecular Biology focused on pseudogenes will certainly not end, nor enter this debate; however, scientists who are also genomics and pseudogene enthusiasts will certainly appreciate that many scientists are thinking about these particular genetic elements in new and interesting ways. With this new interest in a biological significance and "non-junk" role for pseudogenes and other noncoding RNAs, new methods and approaches are being developed to unlock the mystery of these ancient artifacts we know as pseudogenes. In this brief introductory chapter we highlight the renewed interest in pseudogenes and review a rationale for intensification of pseudogene-related research.

Computational methods of identification of pseudogenes based on functionality: entropy and GC content

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

The OCT4 pseudogene POU5F1B is amplified and promotes an aggressive phenotype in gastric cancer

POU5F1B (POU domain class 5 transcription factor 1B), a processed pseudogene that is highly homologous to OCT4, was recently shown to be transcribed in cancer cells, but its clinical relevance and biological function have remained unclear. We now show that POU5F1B, which is located adjacent to MYC on human chromosome 8q24, is frequently amplified in gastric cancer (GC) cell lines. POU5F1B, but not OCT4, was also found to be expressed at a high level in GC cell lines and clinical specimens. In addition, the DNA copy number and mRNA abundance for POU5F1B showed a positive correlation in both cancer cell lines and GC specimens. Overexpression of POU5F1B in GC cells promoted colony formation in vitro as well as both tumorigenicity and tumor growth in vivo, and these effects were enhanced in the additional presence of MYC overexpression. Furthermore, knockdown of POU5F1B expression with a short hairpin RNA confirmed a role for the endogenous pseudogene in the promotion of cancer cell growth in vitro and tumor growth in vivo. POU5F1B overexpression induced upregulation of various growth factors in GC cells as well as exhibited mitogenic, angiogenic and antiapoptotic effects in GC xenografts. Finally, amplification of POU5F1B was detected in 17 (12%) of 145 cases of GC and was a significant predictor of poor prognosis in patients with stage IV disease. In conclusion, we found that the POU5F1B pseudogene is amplified and expressed at a high level in, as well as confers an aggressive phenotype on, GC, and that POU5F1B amplification is associated with a poor prognosis in GC patients.

Meta-analysis using a novel database, miRStress, reveals miRNAs that are frequently associated with the radiation and hypoxia stress-responses

Organisms are often exposed to environmental pressures that affect homeostasis, so it is important to understand the biological basis of stress-response. Various biological mechanisms have evolved to help cells cope with potentially cytotoxic changes in their environment. miRNAs are small non-coding RNAs which are able to regulate mRNA stability. It has been suggested that miRNAs may tip the balance between continued cytorepair and induction of apoptosis in response to stress. There is a wealth of data in the literature showing the effect of environmental stress on miRNAs, but it is scattered in a large number of disparate publications. Meta-analyses of this data would produce added insight into the molecular mechanisms of stress-response. To facilitate this we created and manually curated the miRStress database, which describes the changes in miRNA levels following an array of stress types in eukaryotic cells. Here we describe this database and validate the miRStress tool for analysing miRNAs that are regulated by stress. To validate the database we performed a cross-species analysis to identify miRNAs that respond to radiation. The analysis tool confirms miR-21 and miR-34a as frequently deregulated in response to radiation, but also identifies novel candidates as potentially important players in this stress response, including miR-15b, miR-19b, and miR-106a. Similarly, we used the miRStress tool to analyse hypoxia-responsive miRNAs. The most frequently deregulated miRNAs were miR-210 and miR-21, as expected. Several other miRNAs were also found to be associated with hypoxia, including miR-181b, miR-26a/b, miR-106a, miR-213 and miR-192. Therefore the miRStress tool has identified miRNAs with hitherto unknown or under-appreciated roles in the response to specific stress types. The miRStress tool, which can be used to uncover new insight into the biological roles of miRNAs, and also has the potential to unearth potential biomarkers for therapeutic response, is freely available at http://mudshark.brookes.ac.uk/MirStress.

Not so pseudo: the evolutionary history of protein phosphatase 1 regulatory subunit 2 and related pseudogenes

Background

Pseudogenes are traditionally considered “dead” genes, therefore lacking biological functions. This view has however been challenged during the last decade. This is the case of the Protein phosphatase 1 regulatory subunit 2 (PPP1R2) or inhibitor-2 gene family, for which several incomplete copies exist scattered throughout the genome.

Results

In this study, the pseudogenization process of PPP1R2 was analyzed. Ten PPP1R2-related pseudogenes (PPP1R2P1-P10), highly similar to PPP1R2, were retrieved from the human genome assembly present in the databases. The phylogenetic analysis of mammalian PPP1R2 and related pseudogenes suggested that PPP1R2P7 and PPP1R2P9 retroposons appeared before the great mammalian radiation, while the remaining pseudogenes are primate-specific and retroposed at different times during Primate evolution. Although considered inactive, four of these pseudogenes seem to be transcribed and possibly possess biological functions. Given the role of PPP1R2 in sperm motility, the presence of these proteins was assessed in human sperm, and two PPP1R2-related proteins were detected, PPP1R2P3 and PPP1R2P9. Signatures of negative and positive selection were also detected in PPP1R2P9, further suggesting a role as a functional protein.

Conclusions

The results show that contrary to initial observations PPP1R2-related pseudogenes are not simple bystanders of the evolutionary process but may rather be at the origin of genes with novel functions.

The Janus face of the heme oxygenase/biliverdin reductase system in Alzheimer disease: it's time for reconciliation

Alzheimer disease (AD) is the most common form of dementia among the elderly and is characterized by progressive loss of memory and cognition. These clinical features are due in part to the increase of reactive oxygen and nitrogen species that mediate neurotoxic effects. The up-regulation of the heme oxygenase-1/biliverdin reductase-A (HO-1/BVR-A) system is one of the earlier events in the adaptive response to stress. HO-1/BVR-A reduces the intracellular levels of pro-oxidant heme and generates equimolar amounts of the free radical scavengers biliverdin-IX alpha (BV)/bilirubin-IX alpha (BR) as well as the pleiotropic gaseous neuromodulator carbon monoxide (CO) and ferrous iron. Two main and opposite hypotheses for a role of the HO-1/BVR-A system in AD propose that this system mediates neurotoxic and neuroprotective effects, respectively. This apparent controversy was mainly due to the fact that for over about 20years HO-1 was the only player on which all the analyses were focused, excluding the other important and essential component of the entire system, BVR. Following studies from the Butterfield laboratory that reported alterations in BVR activity along with decreased phosphorylation and increased oxidative/nitrosative post-translational modifications in the brain of subjects with AD and amnestic mild cognitive impairment (MCI) subjects, a debate was opened on the real pathophysiological and clinical significance of BVR-A. In this paper we provide a review of the main discoveries about the HO/BVR system in AD and MCI, and propose a mechanism that reconciles these two hypotheses noted above of neurotoxic and the neuroprotective aspects of this important stress responsive system.

lincRNAs: genomics, evolution, and mechanisms

Long intervening noncoding RNAs (lincRNAs) are transcribed from thousands of loci in mammalian genomes and might play widespread roles in gene regulation and other cellular processes. This Review outlines the emerging understanding of lincRNAs in vertebrate animals, with emphases on how they are being identified and current conclusions and questions regarding their genomics, evolution and mechanisms of action.

Up-regulation of SUMO1 pseudogene 3 (SUMO1P3) in gastric cancer and its clinical association

Long noncoding RNAs (lncRNAs) play crucial roles during cancer occurrence and progression. The pseudogene-expressed lncRNA is one major type of lncRNA family. However, their association with cancers is largely unknown. In this study, we focused on small ubiquitin-like modifier (SUMO) 1 pseudogene 3, SUMO1P3. Gastric cancer tissues and adjacent nontumor tissues were collected from 96 patients with gastric cancer. The SUMO1P3 levels were detected by quantitative reverse transcription-polymerase chain reaction. Then, the association between the level of SUMO1P3 in gastric cancer tissues and the clinicopathological features of patients with gastric cancer was further analyzed. A receiver operating characteristic curve was constructed for differentiating patients with gastric cancer from patients with benign gastric diseases. The results showed that SUMO1P3 was significantly up-regulated in gastric cancer tissues compared with paired-adjacent nontumorous tissues (p < 0.01). Its expression level was significantly correlated with tumor size (p = 0.003), differentiation (p = 0.002), lymphatic metastasis (p = 0.001), and invasion (p = 0.039). The area under the ROC curve of SUMO1P3 was up to 0.666. These results indicated, for the first time, that pseudogene-expressed lncRNA SUMO1P3 may be a potential biomarker in the diagnosis of gastric cancer.

Pseudogene OCT4-pg4 functions as a natural micro RNA sponge to regulate OCT4 expression by competing for miR-145 in hepatocellular carcinoma

The POU transcription factor OCT4 is a pleiotropic regulator of gene expression in embryonic stem cells. Recent studies demonstrated that OCT4 is aberrantly expressed in multiple types of human cancer; however, the underlying molecular mechanism remains largely unknown. In this study, we report that OCT4-pg4, a pseudogene of OCT4, is abnormally activated in hepatocellular carcinoma (HCC). The expression level of OCT4-pg4 is positively correlated with that of OCT4, and both gene transcripts can be directly targeted by a tumor-suppressive micro RNA miR-145. We find that the non-coding RNA OCT4-pg4 is biologically active, as it can upregulate OCT4 protein level in HCC. Mechanistic analysis revealed that OCT4-pg4 functions as a natural micro RNA sponge to protect OCT4 transcript from being inhibited by miR-145. In addition, our study also showed that OCT4-pg4 can promote growth and tumorigenicity of HCC cells, thus exerting an oncogenic role in hepatocarcinogenesis. Furthermore, survival analysis suggests that high OCT4-pg4 level is significantly correlated with poor prognosis of HCC patients. Taken together, our finding adds a new layer of post-transcriptional regulation of OCT4 and sheds new light on the treatment of human HCC.

Non-coding RNAs: multi-tasking molecules in the cell

In the last years it has become increasingly clear that the mammalian transcriptome is highly complex and includes a large number of small non-coding RNAs (sncRNAs) and long noncoding RNAs (lncRNAs). Here we review the biogenesis pathways of the three classes of sncRNAs, namely short interfering RNAs (siRNAs), microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs). These ncRNAs have been extensively studied and are involved in pathways leading to specific gene silencing and the protection of genomes against virus and transposons, for example. Also, lncRNAs have emerged as pivotal molecules for the transcriptional and post-transcriptional regulation of gene expression which is supported by their tissue-specific expression patterns, subcellular distribution, and developmental regulation. Therefore, we also focus our attention on their role in differentiation and development. SncRNAs and lncRNAs play critical roles in defining DNA methylation patterns, as well as chromatin remodeling thus having a substantial effect in epigenetics. The identification of some overlaps in their biogenesis pathways and functional roles raises the hypothesis that these molecules play concerted functions in vivo, creating complex regulatory networks where cooperation with regulatory proteins is necessary. We also highlighted the implications of biogenesis and gene expression deregulation of sncRNAs and lncRNAs in human diseases like cancer.

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.

Long non-coding RNAs in haematological malignancies

Long non-coding RNAs (lncRNAs) are functional RNAs longer than 200 nucleotides in length. LncRNAs are as diverse as mRNAs and they normally share the same biosynthetic machinery based on RNA polymerase II, splicing and polyadenylation. However, lncRNAs have low coding potential. Compared to mRNAs, lncRNAs are preferentially nuclear, more tissue specific and expressed at lower levels. Most of the lncRNAs described to date modulate the expression of specific genes by guiding chromatin remodelling factors; inducing chromosomal loopings; affecting transcription, splicing, translation or mRNA stability; or serving as scaffolds for the organization of cellular structures. They can function in cis, cotranscriptionally, or in trans, acting as decoys, scaffolds or guides. These functions seem essential to allow cell differentiation and growth. In fact, many lncRNAs have been shown to exert oncogenic or tumor suppressor properties in several cancers including haematological malignancies. In this review, we summarize what is known about lncRNAs, the mechanisms for their regulation in cancer and their role in leukemogenesis, lymphomagenesis and hematopoiesis. Furthermore, we discuss the potential of lncRNAs in diagnosis, prognosis and therapy in cancer, with special attention to haematological malignancies.

Evolution of the human gastrokine locus and confounding factors regarding the pseudogenicity of GKN3

In a screen for genes expressed specifically in gastric mucous neck cells, we identified GKN3, the recently discovered third member of the gastrokine family. We present confirmatory mouse data and novel porcine data showing that mouse GKN3 expression is confined to mucous cells of the corpus neck and antrum base and is prominently expressed in metaplastic lesions. GKN3 was proposed originally to be expressed in some human populations and a pseudogene in others. To investigate that hypothesis, we studied human GKN3 evolution in the context of its paralogous genomic neighbors, GKN1 and GKN2. Haplotype analysis revealed that GKN3 mimics GKN2 in patterns of exonic SNP allocation, whereas GKN1 appeared to be more stringently selected. GKN3 showed signatures of both directional selection and population based selective sweeps in humans. One such selective sweep includes SNP rs10187256, originally identified as an ancestral tryptophan to premature STOP codon mutation. The derived (nonancestral) allele went to fixation in Asia. We show that another SNP, rs75578132, identified 5 bp downstream of rs10187256, exhibits a second selective sweep in almost all Europeans, some Latinos, and some Africans, possibly resulting from a reintroduction of European genes during African colonization. Finally, we identify a mutation that would destroy the splice donor site in the putative exon3-intron3 boundary, which occurs in all human genomes examined to date. Our results highlight a stomach-specific human genetic locus, which has undergone various selective sweeps across European, Asian, and African populations and thus reflects geographic and ethnic patterns in genome evolution.

microRNAs and ceRNAs: RNA networks in pathogenesis of cancer

microRNAs (miRNAs) are a class of endogenous, single-stranded non-coding RNAs of 20-23 nucleotides in length, functioning as negative regulators of gene expression at the post-transcriptional level. The dysregulation of miRNAs has been demonstrated to play critical roles in tumorigenesis, either through inhibiting tumor suppressor genes or activating oncogenes inappropriately. Besides their promising clinical applications in cancer diagnosis and treatment, recent studies have uncovered that miRNAs could act as a regulatory language, through which messenger RNAs, transcribed pseudogenes, and long noncoding RNAs crosstalk with each other and form a novel regulatory network. RNA transcripts involved in this network have been termed as competing endogenous RNAs (ceRNAs), since they influence each other's level by competing for the same pool of miRNAs through miRNA response elements (MREs) on their target transcripts. The discovery of miRNA-ceRNA network not only provides the possibility of an additional level of post-transcriptional regulation, but also dictates a reassessment of the existing regulatory pathways involved in cancer initiation and progression.

microRNAs and ceRNAs: RNA networks in pathogenesis of cancer

microRNAs (miRNAs) are a class of endogenous, single-stranded non-coding RNAs of 20-23 nucleotides in length, functioning as negative regulators of gene expression at the post-transcriptional level. The dysregulation of miRNAs has been demonstrated to play critical roles in tumorigenesis, either through inhibiting tumor suppressor genes or activating oncogenes inappropriately. Besides their promising clinical applications in cancer diagnosis and treatment, recent studies have uncovered that miRNAs could act as a regulatory language, through which messenger RNAs, transcribed pseudogenes, and long noncoding RNAs crosstalk with each other and form a novel regulatory network. RNA transcripts involved in this network have been termed as competing endogenous RNAs (ceRNAs), since they influence each other's level by competing for the same pool of miRNAs through miRNA response elements (MREs) on their target transcripts. The discovery of miRNA-ceRNA network not only provides the possibility of an additional level of post-transcriptional regulation, but also dictates a reassessment of the existing regulatory pathways involved in cancer initiation and progression.

Evolutionary paths of the cAMP-dependent protein kinase (PKA) catalytic subunits

3',5'-cyclic adenosine monophosphate (cAMP) dependent protein kinase or protein kinase A (PKA) has served as a prototype for the large family of protein kinases that are crucially important for signal transduction in eukaryotic cells. The PKA catalytic subunits Cα and Cβ, encoded by the two genes PRKACA and PRKACB, respectively, are among the best understood and characterized human kinases. Here we have studied the evolution of this gene family in chordates, arthropods, mollusks and other animals employing probabilistic methods and show that Cα and Cβ arose by duplication of an ancestral PKA catalytic subunit in a common ancestor of vertebrates. The two genes have subsequently been duplicated in teleost fishes. The evolution of the PRKACG retroposon in simians was also investigated. Although the degree of sequence conservation in the PKA Cα/Cβ kinase family is exceptionally high, a small set of signature residues defining Cα and Cβ subfamilies were identified. These conserved residues might be important for functions that are unique to the Cα or Cβ clades. This study also provides a good example of a seemingly simple phylogenetic problem which, due to a very high degree of sequence conservation and corresponding weak phylogenetic signals, combined with problematic nonphylogenetic signals, is nontrivial for state-of-the-art probabilistic phylogenetic methods.

RNAseq versus genome-predicted transcriptomes: a large population of novel transcripts identified in an Illumina-454 Hydra transcriptome

Background

Evolutionary studies benefit from deep sequencing technologies that generate genomic and transcriptomic sequences from a variety of organisms. Genome sequencing and RNAseq have complementary strengths. In this study, we present the assembly of the most complete Hydra transcriptome to date along with a comparative analysis of the specific features of RNAseq and genome-predicted transcriptomes currently available in the freshwater hydrozoan Hydra vulgaris.

Results

To produce an accurate and extensive Hydra transcriptome, we combined Illumina and 454 Titanium reads, giving the primacy to Illumina over 454 reads to correct homopolymer errors. This strategy yielded an RNAseq transcriptome that contains 48’909 unique sequences including splice variants, representing approximately 24’450 distinct genes. Comparative analysis to the available genome-predicted transcriptomes identified 10’597 novel Hydra transcripts that encode 529 evolutionarily-conserved proteins. The annotation of 170 human orthologs points to critical functions in protein biosynthesis, FGF and TOR signaling, vesicle transport, immunity, cell cycle regulation, cell death, mitochondrial metabolism, transcription and chromatin regulation. However, a majority of these novel transcripts encodes short ORFs, at least 767 of them corresponding to pseudogenes. This RNAseq transcriptome also lacks 11’270 predicted transcripts that correspond either to silent genes or to genes expressed below the detection level of this study.

Conclusions

We established a simple and powerful strategy to combine Illumina and 454 reads and we produced, with genome assistance, an extensive and accurate Hydra transcriptome. The comparative analysis of the RNAseq transcriptome with genome-predicted transcriptomes lead to the identification of large populations of novel as well as missing transcripts that might reflect Hydra-specific evolutionary events.

Pseudogenes: pseudo or real functional elements?

Pseudogenes are genomic remnants of ancient protein-coding genes which have lost their coding potentials through evolution. Although broadly existed, pseudogenes used to be considered as junk or relics of genomes which have not drawn enough attentions of biologists until recent years. With the broad applications of high-throughput experimental techniques, growing lines of evidence have strongly suggested that some pseudogenes possess special functions, including regulating parental gene expression and participating in the regulation of many biological processes. In this review, we summarize some basic features of pseudogenes and their functions in regulating development and diseases. All of these observations indicate that pseudogenes are not purely dead fossils of genomes, but warrant further exploration in their distribution, expression regulation and functions. A new nomenclature is desirable for the currently called 'pseudogenes' to better describe their functions.

Long noncoding RNAs: past, present, and future

Long noncoding RNAs (lncRNAs) have gained widespread attention in recent years as a potentially new and crucial layer of biological regulation. lncRNAs of all kinds have been implicated in a range of developmental processes and diseases, but knowledge of the mechanisms by which they act is still surprisingly limited, and claims that almost the entirety of the mammalian genome is transcribed into functional noncoding transcripts remain controversial. At the same time, a small number of well-studied lncRNAs have given us important clues about the biology of these molecules, and a few key functional and mechanistic themes have begun to emerge, although the robustness of these models and classification schemes remains to be seen. Here, we review the current state of knowledge of the lncRNA field, discussing what is known about the genomic contexts, biological functions, and mechanisms of action of lncRNAs. We also reflect on how the recent interest in lncRNAs is deeply rooted in biology's longstanding concern with the evolution and function of genomes.

Carboxylesterase 1 gene duplication and mRNA expression in adipose tissue are linked to obesity and metabolic function

CONTEXT AND AIMS

Carboxylesterase 1 (CES1) appears to play an important role in the control of the metabolism of triglycerides and cholesterol in adipocytes and other cell types including hepatocytes. Therefore, it is relevant to gain insights into the genetic versus non-genetic mechanisms involved in the control of CES1 mRNA expression. Here, we investigated CES1 mRNA expression level in adipose tissue and its association with measures of adiposity and metabolic function in a population of elderly twins. Furthermore, the heritability of CES1 mRNA expression level in adipose tissue and the effect of CES1 gene duplication were assessed.

METHODOLOGY

A total of 295 monozygotic and dizygotic twin subjects (62-83 years) with (n = 48) or without (n = 247) type 2 diabetes mellitus were enrolled in the study. They were subjected to a standard oral glucose tolerance test and excision of abdominal subcutaneous fat biopsies during the fasting state. Levels of CES1 mRNA and copy number of the gene were assessed by quantitative PCR.

RESULTS

CES1 mRNA expression level in adipose tissue was positively associated with body-mass index (P<0.001), homeostasis model assessment-insulin resistance (P = 0.003) and level of fasting glucose (P = 0.002), insulin (P = 0.006), and triglycerides (P = 0.003). The heritability for the expression of CES1 mRNA in adipose tissue was high. CES1 gene duplication was positively associated with insulin sensitivity (P = 0.05) as well as glucose tolerance (P = 0.03) and negatively associated with homeostasis model assessment-insulin resistance (P = 0.02). Duplication of CES1 was not linked to mRNA level of this gene (P = 0.63).

CONCLUSION

CES1 mRNA in adipose tissue appears to be under strong genetic control and was associated with measures of metabolic function raising the possibility of a potential role of this enzyme in the development of type 2 diabetes mellitus. Further studies are needed to understand the potential effect of CES1 gene duplication on adipocyte and whole-body metabolic functions.

Regional differences in gene expression and promoter usage in aged human brains

To characterize the promoterome of caudate and putamen regions (striatum), frontal and temporal cortices, and hippocampi from aged human brains, we used high-throughput cap analysis of gene expression to profile the transcription start sites and to quantify the differences in gene expression across the 5 brain regions. We also analyzed the extent to which methylation influenced the observed expression profiles. We sequenced more than 71 million cap analysis of gene expression tags corresponding to 70,202 promoter regions and 16,888 genes. More than 7000 transcripts were differentially expressed, mainly because of differential alternative promoter usage. Unexpectedly, 7% of differentially expressed genes were neurodevelopmental transcription factors. Functional pathway analysis on the differentially expressed genes revealed an overrepresentation of several signaling pathways (e.g., fibroblast growth factor and wnt signaling) in hippocampus and striatum. We also found that although 73% of methylation signals mapped within genes, the influence of methylation on the expression profile was small. Our study underscores alternative promoter usage as an important mechanism for determining the regional differences in gene expression at old age.

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.

Evidence for conserved post-transcriptional roles of unitary pseudogenes and for frequent bifunctionality of mRNAs

Background

Recent reports have highlighted instances of mRNAs that, in addition to coding for protein, regulate the abundance of related transcripts by altering microRNA availability. These two mRNA roles - one mediated by RNA and the other by protein - are inter-dependent and hence cannot easily be separated. Whether the RNA-mediated role of transcripts is important, per se, or whether it is a relatively innocuous consequence of competition by different transcripts for microRNA binding remains unknown.

Results

Here we took advantage of 48 loci that encoded proteins in the earliest eutherian ancestor, but whose protein-coding capability has since been lost specifically during rodent evolution. Sixty-five percent of such loci, which we term 'unitary pseudogenes', have retained their expression in mouse and their transcripts exhibit conserved tissue expression profiles. The maintenance of these unitary pseudogenes' spatial expression profiles is associated with conservation of their microRNA response elements and these appear to preserve the post-transcriptional roles of their protein-coding ancestor. We used mouse Pbcas4, an exemplar of these transcribed unitary pseudogenes, to experimentally test our genome-wide predictions. We demonstrate that the role of Pbcas4 as a competitive endogenous RNA has been conserved and has outlived its ancestral gene's loss of protein-coding potential.

Conclusions

These results show that post-transcriptional regulation by bifunctional mRNAs can persist over long evolutionary time periods even after their protein coding ability has been lost.

Regulation of eukaryotic gene expression by the untranslated gene regions and other non-coding elements

There is now compelling evidence that the complexity of higher organisms correlates with the relative amount of non-coding RNA rather than the number of protein-coding genes. Previously dismissed as "junk DNA", it is the non-coding regions of the genome that are responsible for regulation, facilitating complex temporal and spatial gene expression through the combinatorial effect of numerous mechanisms and interactions working together to fine-tune gene expression. The major regions involved in regulation of a particular gene are the 5' and 3' untranslated regions and introns. In addition, pervasive transcription of complex genomes produces a variety of non-coding transcripts that interact with these regions and contribute to regulation. This review discusses recent insights into the regulatory roles of the untranslated gene regions and non-coding RNAs in the control of complex gene expression, as well as the implications of this in terms of organism complexity and evolution.

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.

Is prnt a pseudogene? Identification of ram Prt in testis and ejaculated spermatozoa

A hallmark of prion diseases or transmissible spongiform encephalopaties is the conversion of the cellular prion protein (PrP(C)), expressed by the prion gene (prnp), into an abnormally folded isoform (PrP(Sc)) with amyloid-like features that causes scrapie in sheep among other diseases. prnp together with prnd (which encodes a prion-like protein designated as Doppel), and prnt (that encodes the prion protein testis specific--Prt) with sprn (shadow of prion protein gene, that encodes Shadoo or Sho) genes, constitute the "prion gene complex". Whereas a role for prnd in the proper functioning of male reproductive system has been confirmed, the function of prnt, a recently discovered prion family gene, comprises a conundrum leading to the assumption that ruminant prnt is a pseudogene with no protein expression. The main objective of the present study was to identify Prt localization in the ram reproductive system and simultaneously to elucidate if ovine prnt gene is transcribed into protein-coding RNA. Moreover, as Prt is a prnp-related protein, the amyloid propensity was also tested for ovine and caprine Prt. Recombinant Prt was used to immunize BALB/c mice, and the anti-Prt polyclonal antibody (APPA) immune response was evaluated by ELISA and Western Blot. When tested by indirect immunofluorescence, APPA showed high avidity to the ram sperm head apical ridge subdomain, before and after induced capacitation, but did not show the same behavior against goat spermatozoa, suggesting high antibody specificity against ovine-Prt. Prt was also found in the testis when assayed by immunohistochemistry during ram spermatogenesis, where spermatogonia, spermatocytes, spermatids and spermatozoa, stained positive. These observations strongly suggest ovine prnt to be a translated protein-coding gene, pointing to a role for Prt protein in the ram reproductive physiology. Besides, caprine Prt appears to exhibit a higher amyloid propensity than ovine Prt, mostly associated with its phenylalanine residue.

Pseudogenes as weaknesses of ACTB (Actb) and GAPDH (Gapdh) used as reference genes in reverse transcription and polymerase chain reactions

The genes encoding β-actin (ACTB in human or Actb in mouse) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH in human or Gapdh in mouse) are the two most commonly used references for sample normalization in determination of the mRNA level of interested genes by reverse transcription (RT) and ensuing polymerase chain reactions (PCR). In this study, bioinformatic analyses revealed that the ACTB, Actb, GAPDH and Gapdh had 64, 69, 67 and 197 pseudogenes (PGs), respectively, in the corresponding genome. Most of these PGs are intronless and similar in size to the authentic mRNA. Alignment of several PGs of these genes with the corresponding mRNA reveals that they are highly homologous. In contrast, the hypoxanthine phosphoribosyltransferase-1 gene (HPRT1 in human or Hprt in mouse) only had 3 or 1 PG, respectively, and the mRNA has unique regions for primer design. PCR with cDNA or genomic DNA (gDNA) as templates revealed that our HPRT1, Hprt and GAPDH primers were specific, whereas our ACTB and Actb primers were not specific enough both vertically (within the cDNA) and horizontally (compared cDNA with gDNA). No primers could be designed for the Gapdh that would not mis-prime PGs. Since most of the genome is transcribed, we suggest to peers to forgo ACTB (Actb) and GAPDH (Dapdh) as references in RT-PCR and, if there is no surrogate, to use our primers with extra caution. We also propose a standard operation procedure in which design of primers for RT-PCR starts from avoiding mis-priming PGs and all primers need be tested for specificity with both cDNA and gDNA.

Copy number variations of the ATP-binding cassette transporter ABCC6 gene and its pseudogenes

BACKGROUND

The ATP-binding cassette transporter ABCC6 gene is located on chromosome 16 between its two pseudogenes (ABCC6P1 and ABCC6P2). Previously, we have shown that ABCC6P1 is transcribed and affects ABCC6 at the transcriptional level. In this study we aimed to determine copy number variations of ABCC6, ABCC6P1 and ABCC6P2 in different populations. Moreover, we sought to study the transcription pattern of ABCC6 and ABCC6 pseudogenes in 39 different human tissues.

FINDINGS

Genomic DNA from healthy individuals from five populations, Chinese (n = 24), Middle East (n = 20), Mexicans (n = 24), Caucasians (n = 50) and Africans (n = 24), were examined for copy number variations of ABCC6 and its pseudogenes by pyrosequencing and quantitative PCR. Copy number variation of ABCC6 was very rare (2/142; 1.4%). However, one or three copies of ABCC6P1 were relatively common (3% and 8%, respectively). Only one person had a single copy of ABCC6P2 while none had three copies. In Chinese, deletions or duplications of ABCC6P1 were more frequent than in any other population (9/24; 37.5%). The transcription pattern of ABCC6P2 was highly similar to ABCC6 and ABCC6P1, with highest transcription in liver and kidney. Interestingly, the total transcription level of pseudogenes, ABCC6P1 + ABCC6P2, was higher than ABCC6 in most tissues, including liver and kidney.

CONCLUSIONS

Copy number variations of the ABCC6 pseudogenes are quite common, especially in populations of Chinese ancestry. The expression pattern of ABCC6P2 in 39 human tissues was highly similar to that of ABCC6 and ABCC6P1 suggesting similar regulatory mechanisms for ABCC6 and its pseudogenes.

The expression of tubulin cofactor A (TBCA) is regulated by a noncoding antisense Tbca RNA during testis maturation

Background

Recently, long noncoding RNAs have emerged as pivotal molecules for the regulation of coding genes' expression. These molecules might result from antisense transcription of functional genes originating natural antisense transcripts (NATs) or from transcriptional active pseudogenes. TBCA interacts with β-tubulin and is involved in the folding and dimerization of new tubulin heterodimers, the building blocks of microtubules.

Methodology/Principal Findings

We found that the mouse genome contains two structurally distinct Tbca genes located in chromosomes 13 (Tbca13) and 16 (Tbca16). Interestingly, the two Tbca genes albeit ubiquitously expressed, present differential expression during mouse testis maturation. In fact, as testis maturation progresses Tbca13 mRNA levels increase progressively, while Tbca16 mRNA levels decrease. This suggests a regulatory mechanism between the two genes and prompted us to investigate the presence of the two proteins. However, using tandem mass spectrometry we were unable to identify the TBCA16 protein in testis extracts even in those corresponding to the maturation step with the highest levels of Tbca16 transcripts. These puzzling results led us to re-analyze the expression of Tbca16. We then detected that Tbca16 transcription produces sense and natural antisense transcripts. Strikingly, the specific depletion by RNAi of these transcripts leads to an increase of Tbca13 transcript levels in a mouse spermatocyte cell line.

Conclusions/Significance

Our results demonstrate that Tbca13 mRNA levels are post-transcriptionally regulated by the sense and natural antisense Tbca16 mRNA levels. We propose that this regulatory mechanism operates during spermatogenesis, a process that involves microtubule rearrangements, the assembly of specific microtubule structures and requires critical TBCA levels.

Disease-linked microRNA-21 exhibits drastically reduced mRNA binding and silencing activity in healthy mouse liver

MicroRNAs (miRNAs) bind to mRNAs and fine-tune protein output by affecting mRNA stability and/or translation. miR-21 is a ubiquitous, highly abundant, and stress-responsive miRNA linked to several diseases, including cancer, fibrosis, and inflammation. Although the RNA silencing activity of miR-21 in diseased cells has been well documented, the roles of miR-21 under healthy cellular conditions are not well understood. Here, we show that pharmacological inhibition or genetic deletion of miR-21 in healthy mouse liver has little impact on regulation of canonical seed-matched mRNAs and only a limited number of genes enriched in stress response pathways. These surprisingly weak and selective regulatory effects on known and predicted target mRNAs contrast with those of other abundant liver miRNAs such as miR-122 and let-7. Moreover, miR-21 shows greatly reduced binding to polysome-associated target mRNAs compared to miR-122 and let-7. Bioinformatic analysis suggests that reduced thermodynamic stability of seed pairing and target binding may contribute to this deficiency of miR-21. Significantly, these trends are reversed in human cervical carcinoma (HeLa) cells, where miRNAs including miR-21 show enhanced target binding within polysomes and where miR-21 triggers strong degradative activity toward target mRNAs. Taken together, our results suggest that, under normal cellular conditions in liver, miR-21 activity is maintained below a threshold required for binding and silencing most of its targets. Consequently, enhanced association with polysome-associated mRNA is likely to explain in part the gain of miR-21 function often found in diseased or stressed cells.

Type-I prenyl protease function is required in the male germline of Drosophila melanogaster

Many proteins require the addition of a hydrophobic prenyl anchor (prenylation) for proper trafficking and localization in the cell. Prenyl proteases play critical roles in modifying proteins for membrane anchorage. The type I prenyl protease has a defined function in yeast (Ste24p/Afc1p) where it modifies a mating pheromone, and in humans (Zmpste24) where it has been implicated in a disease of premature aging. Despite these apparently very different biological processes, the type I prenyl protease gene is highly conserved, encoded by a single gene in a wide range of animal and plant groups. A notable exception is Drosophila melanogaster, where the gene encoding the type I prenyl protease has undergone an unprecedented series of duplications in the genome, resulting in five distinct paralogs, three of which are organized in a tandem array, and demonstrate high conservation, particularly in the vicinity of the active site of the enzyme. We have undertaken targeted deletion to remove the three tandem paralogs from the genome. The result is a male fertility defect, manifesting late in spermatogenesis. Our results also show that the ancestral type I prenyl protease gene in Drosophila is under strong purifying selection, while the more recent replicates are evolving rapidly. Our rescue data support a role for the rapidly evolving tandem paralogs in the male germline. We propose that potential targets for the male-specific type I prenyl proteases include proteins involved in the very dramatic cytoskeletal remodeling events required for spermatid maturation.