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

A Genomic Analysis of Factors Driving lincRNA Diversification: Lessons from Plants

Transcriptomic analyses from across eukaryotes indicate that most of the genome is transcribed at some point in the developmental trajectory of an organism. One class of these transcripts is termed long intergenic noncoding RNAs (lincRNAs). Recently, attention has focused on understanding the evolutionary dynamics of lincRNAs, particularly their conservation within genomes. Here, we take a comparative genomic and phylogenetic approach to uncover factors influencing lincRNA emergence and persistence in the plant family Brassicaceae, to which Arabidopsis thaliana belongs. We searched 10 genomes across the family for evidence of > 5000 lincRNA loci from A. thaliana. From loci conserved in the genomes of multiple species, we built alignments and inferred phylogeny. We then used gene tree/species tree reconciliation to examine the duplication history and timing of emergence of these loci. Emergence of lincRNA loci appears to be linked to local duplication events, but, surprisingly, not whole genome duplication events (WGD), or transposable elements. Interestingly, WGD events are associated with the loss of loci for species having undergone relatively recent polyploidy. Lastly, we identify 1180 loci of the 6480 previously annotated A. thaliana lincRNAs (18%) with elevated levels of conservation. These conserved lincRNAs show higher expression, and are enriched for stress-responsiveness and cis-regulatory motifs known as conserved noncoding sequences (CNSs). These data highlight potential functional pathways and suggest that CNSs may regulate neighboring genes at both the genomic and transcriptomic level. In sum, we provide insight into processes that may influence lincRNA diversification by providing an evolutionary context for previously annotated lincRNAs.

Regulation of oncogenic genes by MicroRNAs and pseudogenes in human lung cancer

Lung cancer is one of the most common mortal cancer types both for men and women. Several different biomarkers have been analyzed to reveal lung cancer prognosis pathways for developing efficient therapeutics and diagnostic agents. microRNAs (miRNAs) and pseudogenes are critical biomarkers in lung cancer and alteration of their expression levels has been identified in each step of lung cancer tumorigenesis. miRNAs and pseudogenes are crucial gene regulators in normal cells as well as in lung cancer cells, and they have both oncogenic and tumor-suppressive roles in lung cancer tumorigenesis. In this study, we have determined the relationship between lung cancer related oncogenes and miRNAs along with pseudogenes in lung cancer, and the results indicate their potential as biological markers for diagnostic and therapeutic purposes.

A microdeletion at Xq22.2 implicates a glycine receptor GLRA4 involved in intellectual disability, behavioral problems and craniofacial anomalies

Background

Among the 21 annotated genes at Xq22.2, PLP1 is the only known gene involved in Xq22.2 microdeletion and microduplication syndromes with intellectual disability. Using an atypical microdeletion, which does not encompass PLP1, we implicate a novel gene GLRA4 involved in intellectual disability, behavioral problems and craniofacial anomalies.

Case presentation

We report a female patient (DGDP084) with a de novo Xq22.2 microdeletion of at least 110 kb presenting with intellectual disability, motor delay, behavioral problems and craniofacial anomalies. While her phenotypic features such as cognitive impairment and motor delay show overlap with Pelizaeus-Merzbacher disease (PMD) caused by PLP1 mutations at Xq22.2, this gene is not included in our patient’s microdeletion and is not dysregulated by a position effect. Because the microdeletion encompasses only three genes, GLRA4, MORF4L2 and TCEAL1, we investigated their expression levels in various tissues by RT-qPCR and found that all three genes were highly expressed in whole human brain, fetal brain, cerebellum and hippocampus. When we examined the transcript levels of GLRA4, MORF4L2 as well as TCEAL1 in DGDP084′s family, however, only GLRA4 transcripts were reduced in the female patient compared to her healthy mother. This suggests that GLRA4 is the plausible candidate gene for cognitive impairment, behavioral problems and craniofacial anomalies observed in DGDP084. Importantly, glycine receptors mediate inhibitory synaptic transmission in the brain stem as well as the spinal cord, and are known to be involved in syndromic intellectual disability.

Conclusion

We hypothesize that GLRA4 is involved in intellectual disability, behavioral problems and craniofacial anomalies as the second gene identified for X-linked syndromic intellectual disability at Xq22.2. Additional point mutations or intragenic deletions of GLRA4 as well as functional studies are needed to further validate our hypothesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12883-016-0642-z) contains supplementary material, which is available to authorized users.

Expression of evolutionarily novel genes in tumors

The evolutionarily novel genes originated through different molecular mechanisms are expressed in tumors. Sometimes the expression of evolutionarily novel genes in tumors is highly specific. Moreover positive selection of many human tumor-related genes in primate lineage suggests their involvement in the origin of new functions beneficial to organisms. It is suggested to consider the expression of evolutionarily young or novel genes in tumors as a new biological phenomenon, a phenomenon of TSEEN (tumor specifically expressed, evolutionarily novel) genes.

Identification and characterization of toll-like receptors (TLRs) in the Chinese tree shrew (Tupaia belangeri chinensis)

In mammals, the toll-like receptors (TLRs) play a major role in initiating innate immune responses against pathogens. Comparison of the TLRs in different mammals may help in understanding the TLR-mediated responses and developing of animal models and efficient therapeutic measures for infectious diseases. The Chinese tree shrew (Tupaia belangeri chinensis), a small mammal with a close relationship to primates, is a viable experimental animal for studying viral and bacterial infections. In this study, we characterized the TLRs genes (tTLRs) in the Chinese tree shrew and identified 13 putative TLRs, which are orthologs of mammalian TLR1-TLR9 and TLR11-TLR13, and TLR10 was a pseudogene in tree shrew. Positive selection analyses using the Maximum likelihood (ML) method showed that tTLR8 and tTLR9 were under positive selection, which might be associated with the adaptation to the pathogen challenge. The mRNA expression levels of tTLRs presented an overall low and tissue-specific pattern, and were significantly upregulated upon Hepatitis C virus (HCV) infection. tTLR4 and tTLR9 underwent alternative splicing, which leads to different transcripts. Phylogenetic analysis and TLR structure prediction indicated that tTLRs were evolutionarily conserved, which might reflect an ancient mechanism and structure in the innate immune response system. Taken together, TLRs had both conserved and unique features in the Chinese tree shrew.

Evolutionary direction of processed pseudogenes

While some pseudogenes have been reported to play important roles in gene regulation, little is known about the possible relationship between pseudogene functions and evolutionary process of pseudogenes, or about the forces responsible for the pseudogene evolution. In this study, we characterized human processed pseudogenes in terms of evolutionary dynamics. Our results show that pseudogenes tend to evolve toward: lower GC content, strong dinucleotide bias, reduced abundance of transcription factor binding motifs and short palindromes, and decreased ability to form nucleosomes. We explored possible evolutionary forces that shaped the evolution pattern of pseudogenes, and concluded that mutations in pseudogenes are likely determined, at least partially, by neighbor-dependent mutational bias and recombination-associated selection.

Differential Expression of OCT4 Pseudogenes in Pluripotent and Tumor Cell Lines

Objective

The human OCT4 gene, the most important pluripotency marker, can generate at least three different transcripts (OCT4A, OCT4B, and OCT4B1) by alternative splicing. OCT4A is the main isoform responsible for the stemness property of embryonic stem (ES) cells. There also exist eight processed OCT4 pseudogenes in the human genome with high homology to the OCT4A, some of which are transcribed in various cancers. Recent conflicting reports on OCT4 expression in tumor cells and tissues emphasize the need to discriminate the expression of OCT4A from other variants as well as OCT4 pseudogenes.

Materials and Methods

In this experimental study, DNA sequencing confirmed the authenticity of transcripts of OCT4 pseudogenes and their expression patterns were investigated in a panel of different human cell lines by reverse transcription-polymerase chain reaction (RT-PCR).

Results

Differential expression of OCT4 pseudogenes in various human cancer and pluripotent cell lines was observed. Moreover, the expression pattern of OCT4-pseudogene 3 (OCT4-pg3) followed that of OCT4A during neural differentiation of the pluripotent cell line of NTERA-2 (NT2). Although OCT4-pg3 was highly expressed in undifferentiated NT2 cells, its expression was rapidly down-regulated upon induction of neural differentiation. Analysis of protein expression of OCT4A, OCT4-pg1, OCT4-pg3, and OCT4-pg4 by Western blotting indicated that OCT4 pseudogenes cannot produce stable proteins. Consistent with a newly proposed competitive role of pseudogene microRNA docking sites, we detected miR-145 binding sites on all transcripts of OCT4 and OCT4 pseudogenes.

Conclusion

Our study suggests a potential coding-independent function for OCT4 pseudogenes during differentiation or tumorigenesis.

FTH1P3, a Novel H-Ferritin Pseudogene Transcriptionally Active, Is Ubiquitously Expressed and Regulated during Cell Differentiation

Ferritin, the major iron storage protein, performs its essential functions in the cytoplasm, nucleus and mitochondria. The variable assembly of 24 subunits of the Heavy (H) and Light (L) type composes the cytoplasmic molecule. In humans, two distinct genes code these subunits, both belonging to complex multigene families. Until now, one H gene has been identified with the coding sequence interrupted by three introns and more than 20 intronless copies widely dispersed on different chromosomes. Two of the intronless genes are actively transcribed in a tissue-specific manner. Herein, we report that FTH1P3, another intronless pseudogene, is transcribed. FTH1P3 transcript was detected in several cell lines and tissues, suggesting that its transcription is ubiquitary, as it happens for the parental ferritin H gene. Moreover, FTH1P3 expression is positively regulated during the cell differentiation process.

The Reverse Transcriptase Encoded by LINE-1 Retrotransposons in the Genesis, Progression, and Therapy of Cancer

In higher eukaryotic genomes, Long Interspersed Nuclear Element 1 (LINE-1) retrotransposons represent a large family of repeated genomic elements. They transpose using a reverse transcriptase (RT), which they encode as part of the ORF2p product. RT inhibition in cancer cells, either via RNA interference-dependent silencing of active LINE-1 elements, or using RT inhibitory drugs, reduces cancer cell proliferation, promotes their differentiation and antagonizes tumor progression in animal models. Indeed, the non-nucleoside RT inhibitor efavirenz has recently been tested in a phase II clinical trial with metastatic prostate cancer patients. An in-depth analysis of ORF2p in a mouse model of breast cancer showed ORF2p to be precociously expressed in precancerous lesions and highly abundant in advanced cancer stages, while being barely detectable in normal breast tissue, providing a rationale for the finding that RT-expressing tumors are therapeutically sensitive to RT inhibitors. We summarize mechanistic and gene profiling studies indicating that abundant LINE-1-derived RT can “sequester” RNA substrates for reverse transcription in tumor cells, entailing the formation of RNA:DNA hybrid molecules and impairing the overall production of regulatory miRNAs, with a global impact on the cell transcriptome. Based on these data, LINE-1-ORF2 encoded RT has a tumor-promoting potential that is exerted at an epigenetic level. We propose a model whereby LINE1-RT drives a previously unrecognized global regulatory process, the deregulation of which drives cell transformation and tumorigenesis with possible implications for cancer cell heterogeneity.

Diverse regulation of 3' splice site usage

The regulation of splice site (SS) usage is important for alternative pre-mRNA splicing and thus proper expression of protein isoforms in cells; its disruption causes diseases. In recent years, an increasing number of novel regulatory elements have been found within or nearby the 3'SS in mammalian genes. The diverse elements recruit a repertoire of trans-acting factors or form secondary structures to regulate 3'SS usage, mostly at the early steps of spliceosome assembly. Their mechanisms of action mainly include: (1) competition between the factors for RNA elements, (2) steric hindrance between the factors, (3) direct interaction between the factors, (4) competition between two splice sites, or (5) local RNA secondary structures or longer range loops, according to the mode of protein/RNA interactions. Beyond the 3'SS, chromatin remodeling/transcription, posttranslational modifications of trans-acting factors and upstream signaling provide further layers of regulation. Evolutionarily, some of the 3'SS elements seem to have emerged in mammalian ancestors. Moreover, other possibilities of regulation such as that by non-coding RNA remain to be explored. It is thus likely that there are more diverse elements/factors and mechanisms that influence the choice of an intron end. The diverse regulation likely contributes to a more complex but refined transcriptome and proteome in mammals.

High expression of dual-specificity phosphatase 5 pseudogene 1 (DUSP5P1) is associated with poor prognosis in acute myeloid leukemia

The purpose of this study was to investigate the expression status of Dual-Specificity Phosphatase 5 Pseudogene 1 (DUSP5P1) and its clinical relevance in patients with acute myeloid leukemia (AML). Real-time quantitative PCR (RQ-PCR) was performed to detect the status of DUSP5P1 expression in 89 patients with de novo AML and 24 normal controls. The level of DUSP5P1 expression was significantly up-regulated in AML compared to controls (P=0.031). The patients with high expression of DUSP5P1 had higher percentage of blasts in bone marrow (BM) than those without high expression (P=0.027). The occurrence rate of DUSP5P1 high expression was significantly higher in M1 (2/8, 25%) and M2 subtypes (9/33, 27%) than in M3 subtype (0/17, 0%) (P=0.034). At the same time, the frequency of DUSP5P1 high expression in patients with intermediate (13/53, 24%) and poor karyotypes (5/11, 45%) was significantly higher than that in patients with favorable karyotype (0/21, 0%) (P=0.003). Meanwhile, DUSP5P1 high-expressed patients had significantly shorter overall survival (OS) than those with low expression (median 4.5 vs. 10.5 months, respectively, P=0.038). Our findings indicated that high expression of DUSP5P1 may identify high-risk AML patients and is associated with poor prognosis in AML.

Are Human Translated Pseudogenes Functional?

By definition, pseudogenes are relics of former genes that no longer possess biological functions. Operationally, they are identified based on disruptions of open reading frames (ORFs) or presumed losses of promoter activities. Intriguingly, a recent human proteomic study reported peptides encoded by 107 pseudogenes. These peptides may play currently unrecognized physiological roles. Alternatively, they may have resulted from accidental translations of pseudogene transcripts and possess no function. Comparing between human and macaque orthologs, we show that the nonsynonymous to synonymous substitution rate ratio (ω) is significantly smaller for translated pseudogenes than other pseudogenes. In particular, five of 34 translated pseudogenes amenable to evolutionary analysis have ω values significantly lower than 1, indicative of the action of purifying selection. This and other findings demonstrate that some but not all translated pseudogenes have selected functions at the protein level. Hence, neither ORF disruption nor presence of protein product disproves or proves gene functionality at the protein level.

Potential Mechanisms for Cancer Resistance in Elephants and Comparative Cellular Response to DNA Damage in Humans

IMPORTANCE

Evolutionary medicine may provide insights into human physiology and pathophysiology, including tumor biology.

OBJECTIVE

To identify mechanisms for cancer resistance in elephants and compare cellular response to DNA damage among elephants, healthy human controls, and cancer-prone patients with Li-Fraumeni syndrome (LFS).

DESIGN, SETTING, AND PARTICIPANTS

A comprehensive survey of necropsy data was performed across 36 mammalian species to validate cancer resistance in large and long-lived organisms, including elephants (n = 644). The African and Asian elephant genomes were analyzed for potential mechanisms of cancer resistance. Peripheral blood lymphocytes from elephants, healthy human controls, and patients with LFS were tested in vitro in the laboratory for DNA damage response. The study included African and Asian elephants (n = 8), patients with LFS (n = 10), and age-matched human controls (n = 11). Human samples were collected at the University of Utah between June 2014 and July 2015.

EXPOSURES

Ionizing radiation and doxorubicin.

MAIN OUTCOMES AND MEASURES

Cancer mortality across species was calculated and compared by body size and life span. The elephant genome was investigated for alterations in cancer-related genes. DNA repair and apoptosis were compared in elephant vs human peripheral blood lymphocytes.

RESULTS

Across mammals, cancer mortality did not increase with body size and/or maximum life span (eg, for rock hyrax, 1% [95% CI, 0%-5%]; African wild dog, 8% [95% CI, 0%-16%]; lion, 2% [95% CI, 0%-7%]). Despite their large body size and long life span, elephants remain cancer resistant, with an estimated cancer mortality of 4.81% (95% CI, 3.14%-6.49%), compared with humans, who have 11% to 25% cancer mortality. While humans have 1 copy (2 alleles) of TP53, African elephants have at least 20 copies (40 alleles), including 19 retrogenes (38 alleles) with evidence of transcriptional activity measured by reverse transcription polymerase chain reaction. In response to DNA damage, elephant lymphocytes underwent p53-mediated apoptosis at higher rates than human lymphocytes proportional to TP53 status (ionizing radiation exposure: patients with LFS, 2.71% [95% CI, 1.93%-3.48%] vs human controls, 7.17% [95% CI, 5.91%-8.44%] vs elephants, 14.64% [95% CI, 10.91%-18.37%]; P < .001; doxorubicin exposure: human controls, 8.10% [95% CI, 6.55%-9.66%] vs elephants, 24.77% [95% CI, 23.0%-26.53%]; P < .001).

CONCLUSIONS AND RELEVANCE

Compared with other mammalian species, elephants appeared to have a lower-than-expected rate of cancer, potentially related to multiple copies of TP53. Compared with human cells, elephant cells demonstrated increased apoptotic response following DNA damage. These findings, if replicated, could represent an evolutionary-based approach for understanding mechanisms related to cancer suppression.

The role of mitochondrial dysfunction in age-related diseases

The aging process is accompanied by the onset of disease and a general decline in wellness. Insights into the aging process have revealed a number of cellular hallmarks of aging, among these epigenetic alterations, loss of proteostasis, mitochondrial dysfunction, cellular senescence, and stem cell exhaustion. Mitochondrial dysfunction increasingly appears to be a common factor connecting several of these hallmarks, driving the aging process and afflicting tissues throughout the body. Recent research has uncovered a much more complex involvement of mitochondria in the cell than has previously been appreciated and revealed novel ways in which mitochondrial defects feed into disease pathology. In this review we evaluate ways in which problems in mitochondria contribute to disease beyond the well-known mechanisms of oxidative stress and bioenergetic deficits, and we predict the direction that mitochondrial disease research will take in years to come.

Genome-Wide Epigenetic Regulation of Gene Transcription in Maize Seeds

Background

Epigenetic regulation is well recognized for its importance in gene expression in organisms. DNA methylation, an important epigenetic mark, has received enormous attention in recent years as it’s a key player in many biological processes. It remains unclear how DNA methylation contributes to gene transcription regulation in maize seeds. Here, we take advantage of recent technologies to examine the genome-wide association of DNA methylation with transcription of four types of DNA sequences, including protein-coding genes, pseudogenes, transposable elements, and repeats in maize embryo and endosperm, respectively.

Results

The methylation in CG, CHG and CHH contexts plays different roles in the control of gene expression. Methylation around the transcription start sites and transcription stop regions of protein-coding genes is negatively correlated, but in gene bodies positively correlated, to gene expression level. The upstream regions of protein-coding genes are enriched with 24-nt siRNAs and contain high levels of CHH methylation, which is correlated to gene expression level. The analysis of sequence content within CG, CHG, or CHH contexts reveals that only CHH methylation is affected by its local sequences, which is different from Arabidopsis.

Conclusions

In summary, we conclude that methylation-regulated transcription varies with the types of DNA sequences, sequence contexts or parts of a specific gene in maize seeds and differs from that in other plant species. Our study helps people better understand from a genome-wide viewpoint that how transcriptional expression is controlled by DNA methylation, one of the important factors influencing transcription, and how the methylation is associated with small RNAs.

The Biological Functions of Non-coding RNAs: From a Line to a Circle

Non-coding RNAs have gained increasing attention, as their physiological and pathological functions are being gradually uncovered. MicroRNAs are the most well-studied ncRNAs, which play essential roles in translational repression and mRNA degradation. In contrast, long non-coding RNAs are distinguished from other small/short non-coding RNAs by length and regulate chromatin remodeling, gene transcription and posttranscriptional modifications. Recently, circular RNAs have emerged as endogenous, abundant, conserved and stable in mammalian cells. It has been demonstrated that circular RNAs can function as miRNA sponges. Other possible biological functions of circular RNAs are still under investigation. In this review, the biogenesis and biological functions of the three major types of ncRNAs, including miRNAs, lncRNAs and circRNAs, are overviewed. In addition, the role of ncRNAs in human diseases and potential clinical applications of ncRNAs are discussed.

Differentially-Expressed Pseudogenes in HIV-1 Infection

Not all pseudogenes are transcriptionally silent as previously thought. Pseudogene transcripts, although not translated, contribute to the non-coding RNA pool of the cell that regulates the expression of other genes. Pseudogene transcripts can also directly compete with the parent gene transcripts for mRNA stability and other cell factors, modulating their expression levels. Tissue-specific and cancer-specific differential expression of these “functional” pseudogenes has been reported. To ascertain potential pseudogene:gene interactions in HIV-1 infection, we analyzed transcriptomes from infected and uninfected T-cells and found that 21 pseudogenes are differentially expressed in HIV-1 infection. This is interesting because parent genes of one-third of these differentially-expressed pseudogenes are implicated in HIV-1 life cycle, and parent genes of half of these pseudogenes are involved in different viral infections. Our bioinformatics analysis identifies candidate pseudogene:gene interactions that may be of significance in HIV-1 infection. Experimental validation of these interactions would establish that retroviruses exploit this newly-discovered layer of host gene expression regulation for their own benefit.

Functions of noncoding sequences in mammalian genomes

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

Competing endogenous RNA: A novel posttranscriptional regulatory dimension associated with the progression of cancer

The existence of artificial sponges and antisense oligonucleotides designed to decrease the availability of microRNAs (miRNAs), a family of small non-coding RNAs that target RNA transcripts through miRNA response elements (MREs) involved in gene expression, suggests that miRNAs may also be regulated. The wide range of RNA transcripts harboring MREs, termed competing endogenous RNAs (ceRNAs), includes protein-coding messenger RNAs (mRNAs) and non-coding RNAs, for example long non-coding RNAs, pseudogenes and circular RNAs, which compete for a common pool of miRNAs as natural decoys. These ceRNAs are co-regulated and produce large, complex posttranscriptional regulatory networks, which have been implicated in numerous biological processes. The present review discusses recent discoveries that implicate natural microRNA decoys in the development of cancer.

Extensive identification and analysis of conserved small ORFs in animals

Background

There is increasing evidence that transcripts or transcript regions annotated as non-coding can harbor functional short open reading frames (sORFs). Loss-of-function experiments have identified essential developmental or physiological roles for a few of the encoded peptides (micropeptides), but genome-wide experimental or computational identification of functional sORFs remains challenging.

Results

Here, we expand our previously developed method and present results of an integrated computational pipeline for the identification of conserved sORFs in human, mouse, zebrafish, fruit fly, and the nematode C. elegans. Isolating specific conservation signatures indicative of purifying selection on amino acid (rather than nucleotide) sequence, we identify about 2,000 novel small ORFs located in the untranslated regions of canonical mRNAs or on transcripts annotated as non-coding. Predicted sORFs show stronger conservation signatures than those identified in previous studies and are sometimes conserved over large evolutionary distances. The encoded peptides have little homology to known proteins and are enriched in disordered regions and short linear interaction motifs. Published ribosome profiling data indicate translation of more than 100 novel sORFs, and mass spectrometry data provide evidence for more than 70 novel candidates.

Conclusions

Taken together, we identify hundreds of previously unknown conserved sORFs in major model organisms. Our computational analyses and integration with experimental data show that these sORFs are expressed, often translated, and sometimes widely conserved, in some cases even between vertebrates and invertebrates. We thus provide an integrated resource of putatively functional micropeptides for functional validation in vivo.

Electronic supplementary material

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

Structural Bioinformatics Inspection of neXtProt PE5 Proteins in the Human Proteome

One goal of the Human Proteome Project is to identify at least one protein product for each of the ∼20,000 human protein-coding genes. As of October 2014, however, there are 3564 genes (18%) that have no or insufficient evidence of protein existence (PE), as curated by neXtProt; these comprise 2647 PE2-4 missing proteins and 616 PE5 dubious protein entries. We conducted a systematic examination of the 616 PE5 protein entries using cutting-edge protein structure and function modeling methods. Compared to a random sample of high-confidence PE1 proteins, the putative PE5 proteins were found to be over-represented in the membrane and cell surface proteins and peptides fold families. Detailed functional analyses show that most PE5 proteins, if expressed, would belong to transporters and receptors localized in the plasma membrane compartment. The results suggest that experimental difficulty in identifying membrane-bound proteins and peptides could have precluded their detection in mass spectrometry and that special enrichment techniques with improved sensitivity for membrane proteins could be important for the characterization of the PE5 "dark matter" of the human proteome. Finally, we identify 66 high scoring PE5 protein entries and find that six of them were reported in recent mass spectrometry databases; an illustrative annotation of these six is provided. This work illustrates a new approach to examine the potential folding and function of the dubious proteins comprising PE5, which we will next apply to the far larger group of missing proteins comprising PE2-4.

Regulatory roles of LINE-1-encoded reverse transcriptase in cancer onset and progression

LINE-1 retrotransposons encode the reverse transcriptase (RT) enzyme, required for their own mobility, the expression of which is inhibited in differentiated tissues while being active in tumors. Experimental evidence indicate that the inhibition of LINE-1-derived RT restores differentiation in cancer cells, inhibits tumor progression and yields globally reprogrammed transcription profiles. Newly emerging data suggest that LINE-1-encoded RT modulates the biogenesis of miRNAs, by governing the balance between the production of regulatory double-stranded RNAs and RNA:DNA hybrid molecules, with a direct impact on global gene expression. Abnormally high RT activity unbalances the transcriptome in cancer cells, while RT inhibition restores ‘normal’ miRNA profiles and their regulatory networks. This RT-dependent mechanism can target the myriad of transcripts - both coding and non-coding, sense and antisense - in eukaryotic transcriptomes, with a profound impact on cell fates. LINE-1-encoded RT emerges therefore as a key regulator of a previously unrecognized mechanism in tumorigenesis

The Landscape of long noncoding RNA classification

Advances in the depth and quality of transcriptome sequencing have revealed many new classes of long noncoding RNAs (lncRNAs). lncRNA classification has mushroomed to accommodate these new findings, even though the real dimensions and complexity of the noncoding transcriptome remain unknown. Although evidence of functionality of specific lncRNAs continues to accumulate, conflicting, confusing, and overlapping terminology has fostered ambiguity and lack of clarity in the field in general. The lack of fundamental conceptual unambiguous classification framework results in a number of challenges in the annotation and interpretation of noncoding transcriptome data. It also might undermine integration of the new genomic methods and datasets in an effort to unravel the function of lncRNA. Here, we review existing lncRNA classifications, nomenclature, and terminology. Then, we describe the conceptual guidelines that have emerged for their classification and functional annotation based on expanding and more comprehensive use of large systems biology-based datasets.

Aquaporin10 is a pseudogene in cattle and their relatives

Background

Although AQP10 is mainly expressed in the human GI tract, its physiological role is unclear. In fact, we previously reported that mouse AQP10 is a pseudogene. It is possible that AQP10 is also a pseudogene in other animals.

Methods

Genome databases were searched for AQP10 orthologs and the genomic DNA of each candidate pseudogene was sequenced to confirm its mutations. The expression of the AQP10 mRNA was examined by RT-PCR in the small intestine where human AQP10 is highly expressed.

Results

The genomic database of some mammals had insertions and deletions in the exons of the AQP10 gene, including cattle (Bos taurus), sheep (Ovis aries) and goats (Capra hircus). In the bovine AQP10 gene, exon 1 and 5 had deletions resulting in a frame-shift or a premature termination, respectively, which were confirmed by the direct exon sequencing of the genomic DNA. In the RT-PCR experiments, the PCR primer sets for exon 1/2 and exon 4/5 failed to detect the bands for AQP10 mRNA in the duodenum and jejunum. Similar AQP10 gene mutations were also confirmed in the genomic DNA from sheep and goats. Although these animals were derived from porcine ancestors, the exons of the swine (Sus scrofa) AQP10 gene were complete without mutations. Therefore, AQP10 gene might have turned to a pseudogene around 65 million years before when cattle evolved from porcine ancestors.

Conclusion

AQP10 of ruminantia which regurgitate and rechew their food may have lost its role possibly due to the redundant expression of other aquaglyceroporins.

A LINE-1-encoded reverse transcriptase-dependent regulatory mechanism is active in embryogenesis and tumorigenesis

LINE-1 (long interspersed nuclear elements) retrotransposons constitute a large family of retrotransposable elements, accounting for 17% of the human genome. They encode proteins required for their own mobilization, including a reverse transcriptase (RT) enzyme highly expressed in mouse embryos and mouse and human cancer cells and repressed in somatic differentiated healthy cells. We have found that reverse transcription takes place in early murine embryos, yielding an increase in LINE-1 copy number during preimplantation development, which also occurs in tumor progression. RT inhibition irreversibly arrests embryo development, reduces cancer cell proliferation, promotes differentiation, antagonizes tumor growth, and causes a global reprogramming of transcription profiles. These results strongly suggest that a previously unrecognized RT-dependent regulatory mechanism operates during preimplantation development, is repressed during differentiation to normal tissues, and, when erroneously reactivated in adult life, promotes cell transformation and cancer progression by "resurrecting" embryonic transcriptional pathways. The RT-dependent mechanism emerges as a major source of genetic and epigenetic changes with physiological, pathological, and evolutionary implications.

Pseudogene-derived endogenous siRNAs and their function

Pseudogenes were once considered genomic fossils, but recent studies indicate that they may function as gene regulators through the generation of endogenous small interfering RNAs (esiRNAs), antisense RNAs, and decoys for microRNAs. In this review, we summarize pseudogene study methods, emphasizing relevant publicly available resources, and we describe a systematic pipeline to identify pseudogene-derived esiRNAs and their targets, which can lead to a deeper understanding of pseudogene function.

Making novel proteins from pseudogenes

MOTIVATION

Recently, we made synthetic proteins from non-coding DNA of Escherichia coli. Encouraged by this, we asked: can we artificially express pseudogenes into novel and functional proteins? What kind of structures would be generated? Would these proteins be stable? How would the organism respond to the artificial reactivation of pseudogenes?

RESULTS

To answer these questions, we studied 16 full-length protein equivalents of pseudogenes. The sequence-based predictions indicated interesting molecular and cellular functional roles for pseudogene-derived proteins. Most of the proteins were predicted to be involved in the amino acid biosynthesis, energy metabolism, purines and pyrimidine biosynthesis, central intermediary metabolism, transport and binding. Interestingly, many of the pseudogene-derived proteins were predicted to be enzymes. Furthermore, proteins showed strong evidence of stable tertiary structures. The prediction scores for structure, function and stability were found to be favorable in most of the cases.

IMPACT

To our best knowledge, this is the first such report that predicts the possibility of making functional and stable proteins from pseudogenes. In future, it would be interesting to experimentally synthesize and validate these predictions.

Noncoding RNAs in diabetes vascular complications

Diabetes mellitus is the most common metabolic disorder and is recognised as a dominant health threat of our time. Diabetes induces a widespread damage of the macro- and microvasculature in different organs and tissues and disrupts the endogenous vascular repair mechanisms, thus causing diffuse and severe complications. Moreover, diabetic patients respond poorly to surgical interventions aiming to "revascularise" (i.e., to restore blood flow supply) the ischemic myocardium or lower limbs. The molecular causes underpinning diabetes vascular complications are still underappreciated and druggable molecular targets for therapeutic interventions have not yet clearly emerged. Moreover, diabetes itself and diabetes complications are often silent killers, requiring new prognostic, diagnostic and predictive biomarkers for use in the clinical practice. Noncoding RNA (ncRNAs) are emerging as new fundamental regulators of gene expression. The small microRNAs (miRNAs, miRs) have opened the field capturing the attention of basic and clinical scientists for their potential to become new therapeutic targets and clinical biomarkers. More recently, long ncRNAs (lncRNAs) have started to be actively investigated, leading to first exciting reports, which further suggest their important and yet largely unexplored contribution to vascular physiology and disease. This review introduces the different ncRNA types and focuses at the ncRNA roles in diabetes vascular complications. Furthermore, we discuss the potential value of ncRNAs as clinical biomarkers, and we examine the possibilities for therapeutic intervention targeting ncRNs in diabetes. This article is part of a Special Issue titled: Non-coding RNAs.

Gene expression profile analysis identifies metastasis and chemoresistance-associated genes in epithelial ovarian carcinoma cells

The purpose of this study was to identify genes that associated with higher ability of metastasis and chemotherapic resistance in epithelial ovarian carcinoma (EOC) cells. An oligonucleotide microarray with probe sets complementary to 41,000⁺ unique human genes and transcripts was used to determine whether gene expression profile may differentiate three epithelial ovarian cell lines (RMG-I-C, COC1 and HO8910) from their sub-lines (RMG-I-H, COCI/DDP and HO8910/PM) with higher ability of metastasis and chemotherapic resistance. Quantitative real-time PCR and immunohistochemical staining validated the microarray results. Hierarchic cluster analysis of gene expression identified 49 genes that exhibited ≥2.0-fold change and P value ≤0.05. Highly differential expression of GCET2, NLRP4, FOXP1 and SNX29 genes was validated by quantitative PCR in all cell line samples. Finally, FOXP1 was validated at the protein level by immunohistochemistry in paraffin embedded ovarian tissues (i.e., for metastasis, 15 primary EOC and 10 omental metastasis [OM]; for chemoresistance, 13 sensitive and 13 resistant EOC). The identification of higher ability of metastasis and chemotherapic resistance-associated genes may provide a foundation for the development of new type-specific diagnostic strategies and treatment for metastasis and chemotherapic resistance in epithelial ovarian cancer.

The pluripotency factor Nanog is directly upregulated by the androgen receptor in prostate cancer cells

BACKGROUND

The Androgen Receptor (AR) is a nuclear hormone receptor that functions as a critical oncogene in all stages of prostate cancer progression, including progression to castration-resistance following androgen-deprivation therapy. Thus, identifying and targeting critical AR-regulated genes is one potential method to block castration-resistant cancer proliferation. Of particular importance are transcription factors that regulate stem cell pluripotency; many of these genes are emerging as critical oncogenes in numerous tumor cell types. Of these, Nanog has been previously shown to increase the self-renewal and stem-like properties of prostate cancer cells. Thus, we hypothesized that Nanog is a candidate AR target gene that may impart castration-resistance.

METHODS

We modulated AR signaling in LNCaP prostate cancer cells and assayed for Nanog expression. Direct AR binding to the NANOG promoter was tested using AR Chromatin Immunoprecipation (ChIP) and analyses of publically available AR ChIP-sequencing data-sets. Nanog over-expressing cells were analyzed for cell growth and cytotoxicity in response to the AR antagonist enzalutamide and the microtubule stabilizing agent docetaxel.

RESULTS

AR signaling upregulates Nanog mRNA and protein. AR binds directly to the NANOG promoter, and was not identified within 75 kb of the NANOGP8 pseudogene, suggesting the NANOG gene locus was preferentially activated. Nanog overexpression in LNCaP cells increases overall growth, but does not increase resistance to enzalutamide or docetaxel.

CONCLUSIONS

Nanog is a novel oncogenic AR target gene in prostate cancer cells, and stable expression of Nanog increases proliferation and growth of prostate cancer cells, but not resistance to enzalutamide or docetaxel.

19q13.33→qter trisomy in a girl with intellectual impairment and seizures

Rearrangements in chromosome 19 are rare. Among the 35 patients with partial 19q trisomy described, only six have a breakpoint defined by array. The 19q duplication results in a variable phenotype, including dysmorphisms, intellectual disability and seizure. In a female patient, although G-banding at 550 band-resolution was normal, multiplex ligation-dependent probe amplification (MLPA) technique and genomic array showed a 10.6 Mb terminal duplication of chromosome 19q13. Fluorescent in situ hybridization (FISH) revealed that the duplicated region was attached to the short arm of chromosome 21 and silver staining showed four small acrocentrics with nucleolar organization region (NOR) activity, suggesting that the breakpoint in chromosome 21 was at p13. This is the first de novo translocation between 19q13.33 and 21p13 described in liveborn. The chromosome 19 is known to be rich in coding and non-coding regions, and chromosomal rearrangements involving this chromosome are very harmful. Furthermore, the 19q13.33→qter region is dense in pseudogenes and microRNAs, which are potent regulators of gene expression. The trisomic level of this region may contribute to deregulation of global gene expression, and consequently, may lead to abnormal development on the carriers of these rearrangements.

•

The first patient with a de novo translocation between 19q13.33 and 21p13 reported in liveborn.

•

The patient clinical and cytogenetic analyses are reported in details.

•

Rearrangements in 19q13.33→qter region are correlated to intellectual disability and seizures.

•

Chromosomal rearrangements involving rich coding and non-coding regions appear to be very harmful.

The parasite Trichomonas vaginalis expresses thousands of pseudogenes and long non-coding RNAs independently from functional neighbouring genes

Background

The human pathogen Trichomonas vaginalis is a parabasalian flagellate that is estimated to infect 3% of the world’s population annually. With a 160 megabase genome and up to 60,000 genes residing in six chromosomes, the parasite has the largest genome among sequenced protists. Although it is thought that the genome size and unusual large coding capacity is owed to genome duplication events, the exact reason and its consequences are less well studied.

Results

Among transcriptome data we found thousands of instances, in which reads mapped onto genomic loci not annotated as genes, some reaching up to several kilobases in length. At first sight these appear to represent long non-coding RNAs (lncRNAs), however, about half of these lncRNAs have significant sequence similarities to genomic loci annotated as protein-coding genes. This provides evidence for the transcription of hundreds of pseudogenes in the parasite. Conventional lncRNAs and pseudogenes are expressed in Trichomonas through their own transcription start sites and independently from flanking genes in Trichomonas. Expression of several representative lncRNAs was verified through reverse-transcriptase PCR in different T. vaginalis strains and case studies exclude the use of alternative start codons or stop codon suppression for the genes analysed.

Conclusion

Our results demonstrate that T. vaginalis expresses thousands of intergenic loci, including numerous transcribed pseudogenes. In contrast to yeast these are expressed independently from neighbouring genes. Our results furthermore illustrate the effect genome duplication events can have on the transcriptome of a protist. The parasite’s genome is in a steady state of changing and we hypothesize that the numerous lncRNAs could offer a large pool for potential innovation from which novel proteins or regulatory RNA units could evolve.

Electronic supplementary material

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

Noncoding RNA response to xenobiotic exposure: an indicator of toxicity and carcinogenicity

INTRODUCTION

Human exposure to certain environmental and occupational chemicals is one of the major risk factors for noncommunicable diseases, including cancer. Therefore, it is desirable to take advantage of subtle exposure-related adverse cellular events for early disease detection and to identify potential dangers caused by new and currently under-evaluated drugs and chemicals. Nongenotoxic events due to carcinogen/toxicant exposure are a general hallmark of sustained cellular stress leading to tumorigenesis. These processes are globally regulated via noncoding RNAs (ncRNAs). Tumorigenesis-associated genotoxic and nongenotoxic events lead to the altered expression of ncRNAs and may provide a mechanistic link between chemical exposure and tumorigenesis. Current advances in toxicogenomics are beginning to provide valuable insight into gene-chemical interactions at the transcriptome level.

AREAS COVERED

In this review, we summarize recent information about the impact of xenobiotics on ncRNAs. Evidence highlighted in this review suggests a critical role of ncRNAs in response to carcinogen/toxicant exposure.

EXPERT OPINION

Benefits for the use of ncRNAs in carcinogenicity assessment include remarkable tissue specificity, early appearance, low baseline variability, and their presence and stability in biological fluids, which suggests that the incorporation of ncRNAs in the evaluation of cancer risk assessment may enhance substantially the efficiency of toxicity and carcinogenicity testing.

Long non-coding RNAs in the regulation of the immune response

•

Widespread changes in lncRNA expresssion are associated with the immune response.

•

lncRNAs regulate the inflammatory response following activation of innate immunity.

•

lncRNAs regulate T cell differentiation and migration.

•

The action of long non-coding RNAs is mediated via diverse mechanisms.

It is increasingly clear that long non-coding RNAs (lncRNAs) regulate a variety biological responses, and that they do so by a diverse range of mechanisms. In the field of immunology, recent publications have shown widespread changes in the expression of lncRNAs during the activation of the innate immune response and T cell development, differentiation, and activation. These lncRNAs control important aspects of immunity such as production of inflammatory mediators, differentiation, and cell migration through regulating protein–protein interactions or via their ability to basepair with RNA and DNA. We review the current understanding of the mechanism of action of these immune-related lncRNAs, discuss their impact on physiological and pathological processes, and highlight important areas of inquiry at the intersection between immunology and lncRNA biology.

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.