Genomics: protein fossils live on as RNA

The Elusive Trypanosoma cruzi Disperse Gene Protein Family (DGF-1)

Chagas disease, caused by Trypanosoma cruzi infections, is included in the group of neglected diseases, and efforts to develop new therapeutic or immunoprevention approaches have not been successful. After the publication of the T. cruzi genome, the number of molecular and biochemical studies on this parasite has increased considerably, many of which are focused on families of variant surface proteins, especially trans-sialidases, mucins, and mucin-associated proteins. The disperse gene protein 1 family (DGF-1) is one of the most abundant families in the T. cruzi genome; however, the large gene size, high copy numbers, and low antibody titers detected in infected humans make it an unattractive study target. However, here we argue that given the ubiquitous presence in all T. cruzi species, and physicochemical characteristics, the DGF-1 gene family may play and important role in host-parasite interactions.

Pseudogene SNRPFP1 derived long non-coding RNA facilitates hepatocellular carcinoma progress in vitro by sponging tumor-suppressive miR-126-5p

Pseudogene-derived transcripts, especially those barely transcribed in normal tissues, have been regarded as a kind of non-coding RNAs, and present potential functions in tumorigenicity and tumor development in human beings. However, their exact effects on hepatocellular carcinoma (HCC) remain largely unknown. On basis of our previous research and the constructed online database for the non-coding RNAs related to HCC, a series of pseudogene transcripts have been discovered, and SNRPFP1, the homologous pseudogene of SNRPF, was found to produce an anomalously high expression long non-coding RNA in HCC. In this study, we validated the expression of the SNRPFP1 transcript in both HCC tissues and cell lines. The adverse correlation between SNRPFP1 expression and patients' outcomes was observed. And depletion of SNRPF1 in HCC cells significantly suppressed cell proliferation and apoptosis resistance. Meanwhile, the motility of HCC cells was potently impaired. Interestingly, miR-126-5p, one of the tumor-suppressive genes commonly decreased in HCC, was found negatively expressed and correlated with SNRPF1, and a specific region of SNRPF1 transcript is directly binding to miR-126-5p in a molecular sponge way. The rescue experiment by knock-out miR-126-5p significantly reversed the cell growth suppression and a higher ratio of cell apoptosis induced by SNRPF1 depletion. Lastly, we concluded that SNRPF1 is a pseudogene active in HCC, and its abnormally over-expressed transcript is a strong promoter of HCC cell progress in vitro by sponging miR-126-5p. We believe that the findings in this study provide new strategies for HCC prevention and therapeutic treatment.

Pseudogene UBE2MP1 derived transcript enhances in vitro cell proliferation and apoptosis resistance of hepatocellular carcinoma cells through miR-145-5p/RGS3 axis

Pseudogenes are barely transcribed at normal, while the anomalous transcripts of them are mostly regarded as long non-coding RNAs (lncRNAs), which play potential functions in human tumorigenicity and development. The exact effects of pseudogene-derived transcripts on hepatocellular carcinoma (HCC) are ambiguous. According to our previous research and constructed database on the HCC-related lncRNAs, we noticed that UBE2MP1 was transcriptionally activated in HCC as a pseudogene from the ubiquitin-conjugating enzyme member UBE2M. In this study, we validated the high expression of the UBE2MP1 transcript in HCC and its adverse correlation with dismal outcomes for the patients. UBE2MP1 depletion at the transcript level significantly impaired cell proliferation and apoptosis resistance in HCC cell lines. Notably, we discovered that the UBE2MP1 transcript shared a specific sequence, binding to the miR-145-5p seed region with a typical ceRNA effect. Simultaneously, we verified an axis of miR-145-5p/RGS3 in HCC cells, which promoted cell proliferation and apoptosis resistance with significance. And modulation of UE2MP1 could remarkably affect RGS3 expression and consequentially influence HCC cell growth in vitro. And combined with the rescue experiment modulating either miR-145-5p or RGS3 furtherly indicated UBE2MP1 as an upstream regulator of the axis in promoting HCC cell growth and maintenance. Thus, our findings provide new strategies for HCC prevention and individual treatment.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

Isolation and Identification of Bioactive Compounds from Bidens spp. Using HPLC-DAD and GC-MS Analysis and Their Biological Activity as Anticancer Molecules

The genus Bidens a member of family Compositae, is widely documented as an ethno-medicinally important genus of plants. In the present study, anticancer potential of three ethno-medicinally important species i.e., B. bipinnata, B. biternata and B. pilosa were tested. For in-vitro evaluation, an MTT (Thiazolyl blue tetrazolium bromide) assay was performed against cervical cancer cells (HeLa), hepatocellular carcinoma (HepG), and adenocarcinoma human alveolar basal epithelial cells (A549). For in vivo evaluation, Artemia salina, Danio rerio, and Caenorhabditis elegans were used. Among all the tested extracts, the ethanol extract of B. biternata appeared to have highest anticancer activity, and the compounds responsible for this activity were identified to be Tris (2,4-di-tert-butylphenyl), 4-hydroxy-2,4′-dimethoxychalcone, and 2,4-di-tert-butylphenol. This is the first report of the isolation of Tris (2,4-di-tert-butylphenyl) phosphate from the genus Bidens and the first report of 4-hydroxy-2,4′-dimethoxychalcone and 2,4-di-tert-butylphenol from B. biternata. Among the isolated compounds, 4-hydroxy-2,4′-dimethoxychalcone showed the highest anticancer activity with an LD50 value of 236.7 µg/mL. Therefore, this compound carries promising potential for being established as a pharmaceutical for chemoprevention and chemotherapy.

GC content of plant genes is linked to past gene duplications

The frequency of G and C nucleotides in genomes varies from species to species, and sometimes even between different genes in the same genome. The monocot grasses have a bimodal distribution of genic GC content absent in dicots. We categorized plant genes from 5 dicots and 4 monocot grasses by synteny to related species and determined that syntenic genes have significantly higher GC content than non-syntenic genes at their 5`-end in the third position within codons for all 9 species. Lower GC content is correlated with gene duplication, as lack of synteny to distantly related genomes is associated with past interspersed gene duplications. Two mutation types can account for biased GC content, mutation of methylated C to T and gene conversion from A to G. Gene conversion involves non-reciprocal exchanges between homologous alleles and is not detectable when the alleles are identical or heterozygous for presence-absence variation, both likely situations for genes duplicated to new loci. Gene duplication can cause production of siRNA which can induce targeted methylation, elevating mC→T mutations. Recently duplicated plant genes are more frequently methylated and less likely to undergo gene conversion, each of these factors synergistically creating a mutational environment favoring AT nucleotides. The syntenic genes with high GC content in the grasses compose a subset that have undergone few duplications, or for which duplicate copies were purged by selection. We propose a “biased gene duplication / biased mutation” (BDBM) model that may explain the origin and trajectory of the observed link between duplication and genic GC bias. The BDBM model is supported by empirical data based on joint analyses of 9 angiosperm species with their genes categorized by duplication status, GC content, methylation levels and functional classes.

Housekeeping Genes for Parkinson's Disease in Humans and Mice

A critical aspect of real-time PCR is the presence of housekeeping genes (HKGs) as an internal control for the normalization of expression data for genes of interest. It is necessary to select correct HKGs in the investigation of various pathologies. Thereby, we analyzed the stability of expression of the HKGs in Parkinson’s disease (PD). The work was carried out in the peripheral blood of patients with PD and in the brain tissues and peripheral blood of mice with MPTP-induced PD. As a result, Aars was the most stably expressed HKG in the mouse brain as a whole. However, different genes were more stably expressed in different parts of the brain. Polr2f was the most stably expressed in the cortex, Psmd6 was the most stably expressed in the cerebellum, and Psmd7 was the most stably expressed in the striatum and substantia nigra. HKGs were different in similar tissues of the studied organisms. Polr2f was the most stably expressed HKG in the peripheral blood of mice, whereas PSMD6 was the most stably expressed gene in humans. Thus, there is no universal HKG both for different brain tissues of one organism and for similar tissues of different organisms. Furthermore, the identified most stably expressed HKGs can be considered as such only under conditions in PD.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

An adaptable dry lab for SYBR based RT-qPCR primer design to reinforce concepts in molecular biology and nucleic acids

The real time PCR (qPCR) method provides a powerful method to assess levels of particular species of DNA. When combined with reverse transcription (RT-qPCR) it is the predominate technique to measure expression of gene transcripts. While this approach is very powerful, particular care must be taken in the design of the primers to facilitate specific and sensitive detection. Herein describes the framework for an undergraduate assignment which aims to teach primer design for SYBR based RT-qPCR. Beyond gaining direct experience with primer design, students will gain familiarity with important bioinformatic resources as well as a deeper theoretical understanding of the RT-qPCR approach and potential limitations. Moreover, as students' progress through the assignment they re-encounter many important concepts in molecular biology, gene expression, and nucleic acids, creating an opportunity for spiral learning. As this exercise only requires access to free web-based resources and does not require a laboratory it can be used in most science education settings. Despite not being a wet lab, this is a highly authentic research experience as this design process is commonplace in a molecular biology laboratory. Furthermore, the assignment is highly adaptable for different learning outcomes, time frames, and student background and ability. This article seeks to highlight connections and expanded learning outcomes for those already teaching such material, as well as a step-by-step guide for those new to teaching such content.

Brucella Genomics: Macro and Micro Evolution

Brucella organisms are responsible for one of the most widespread bacterial zoonoses, named brucellosis. The disease affects several species of animals, including humans. One of the most intriguing aspects of the brucellae is that the various species show a ~97% similarity at the genome level. Still, the distinct Brucella species display different host preferences, zoonotic risk, and virulence. After 133 years of research, there are many aspects of the Brucella biology that remain poorly understood, such as host adaptation and virulence mechanisms. A strategy to understand these characteristics focuses on the relationship between the genomic diversity and host preference of the various Brucella species. Pseudogenization, genome reduction, single nucleotide polymorphism variation, number of tandem repeats, and mobile genetic elements are unveiled markers for host adaptation and virulence. Understanding the mechanisms of genome variability in the Brucella genus is relevant to comprehend the emergence of pathogens.

Pseudogene AKR1B10P1 enhances tumorigenicity and regulates epithelial-mesenchymal transition in hepatocellular carcinoma via stabilizing SOX4

Pseudogenes exert potential functions in tumorigenicity and tumour process in human beings. In our previous research on oncogene AKR1B10 in hepatocellular carcinoma (HCC), its pseudogene, AKR1B10P1, was preliminarily noticed being anomalistic transcribed, whereas whether AKR1B10P1 plays any specific function in HCC is poorly understood. By using shRNA transfection and lentiviral infection, we regulated the expression of ARK1B10P1 transcript and the relative targets in two ways. As we discovered, pathological transcription of AKR1B10P1 in HCC cells significantly promotes cell growth and motility either in vitro or in vivo. AKR1B10P1 was correlated with relatively dismal features of HCC. The epithelial‐mesenchymal transition (EMT) was enhanced by up‐regulating AKR1B10P1. And, a potential sequence of AKR1B10P1 transcript was discovered directly interacting with miR‐138. SOX4, a pivotal promotor of EMT, was validated as the down‐streaming target of miR‐138. Mechanistically, degradation of SOX4 mRNA induced by miR‐138 was effectively abrogated by AKR1B10P1. In conclusion, pseudogene AKR1B10P1 exerts stabilizing effect on SOX4 in HCC, associated EMT process, by directly sponging miR‐138, which post‐transcriptionally modulates SOX4’s regulating gene.

Identification and characterization of long non-coding RNAs as competing endogenous RNAs in the cold stress response of Triticum aestivum

Long non-coding RNAs (lncRNAs) play important roles in plant development and stress responses. MicroRNAs (miRNAs) are involved in transcriptional and post-transcriptional gene regulation. It is not clear how lncRNA-mediated plant responses to cold stress and how lncRNAs, miRNAs and target mRNAs cooperate subject to the competing endogenous RNA (ceRNA). We interpreted the function of lncRNAs in the winter wheat cultivar Dongnongdongmai 1 (Dn1). A total of 9970 putative lncRNAs were initially identified from three Dn1 lncRNA libraries (5 °C, -10 °C and -25 °C) using high-throughput sequencing. Among the 14,626 genes detected via weighted gene co-expression network analysis, 7435 lncRNAs were co-expressed with 7191 mRNAs. We found six modules related to cold resistance in the lncRNA-mRNA weighted co-expression network, and the functions of mRNAs were similar in each module. Antioxidant systems and hormones played important roles in low-temperature responses. RNA sequencing analysis revealed that interactions between the 384 lncRNAs and 70 miRNAs were required for ceRNA activity. According to ceRNA activity, 225 lncRNAs, 60 miRNAs and 621 target mRNAs were involved in the regulatory networks of the cold stress response. Notably, a conserved region was found in the complementary regions of lncRNAs and miR164/408 but had reverse expression trends in the ceRNA network. Our results reveal possible roles of lncRNAs-mRNAs in the regulatory networks associated with tolerance to low temperature and provide useful information for more strategic use of genomic resources in wheat breeding.

Comprehensive analysis of pseudogene HSPB1P1 and its potential roles in hepatocellular carcinoma

The incidence and mortality rate of hepatocellular carcinoma (HCC) nowadays is still at high levels. The regulatory roles of pseudogene in cancers have been gradually recognized in recent years. However, comprehensive investigation of abnormally expressed pseudogene and related mechanisms in HCC remains lacking. GSE124535 dataset was used to identify differentially expressed pseudogenes in HCC tissues compared with normal tissues. Prognostic value of these differentially expressed pseudogenes was analyzed at GEPIA. StarBase used to analyze microRNAs (miRNAs) can bind with pseudogene, while the targets for these miRNAs were analyzed at miRTarBase. Protein-protein interaction (PPI) network was then established for miRNA targets, after that hub genes were selected. Expression correlation of pseudogene and hub genes was analyzed at StarBase. In total, 16 upregulated and 17 downregulated pseudogenes were identified. Pseudogene HSPB1P1 was identified abnormally expressed in 20 types of human cancers and could be used as an indicator for poorer overall survival of patients with HCC. Functional analyses showed that HSPB1P1 was strongly correlated with signaling pathways related to cancer progression. Further studied revealed that HSPB1P1 could direct regulate the EZH2 expression in HCC. In summary, our study indicated that HSPB1P1 was a predictor for poorer overall survival of patients with HCC and may be potential therapeutic target against HCC.

Overexpressed pseudogenes, DUXAP8 and DUXAP9, promote growth of renal cell carcinoma and serve as unfavorable prognostic biomarkers

BACKGROUND

Growing studies have reported that pseudogenes play key roles in multiple human cancers. However, expression and roles of pseudogenes in renal cell carcinoma remains absent.

RESULTS

31 upregulated and 16 downregulated pseudogenes were screened. Higher expression of DUXAP8 and DUXAP9 indicated poorer prognosis of kidney cancer. 33 and 5 miRNAs were predicted to potentially binding to DUXAP8 and DUXAP9, respectively. miR-29c-3p was identified as the most potential binding miRNAs of DUXAP8 and DUXAP9 based on expression, survival and correlation analyses. 254 target genes of miR-29c-3p were forecast. 47 hub genes with node degree >= 10 were identified. Subsequent analysis for the top 10 hub genes demonstrated that COL1A1 and COL1A2 may be two functional targets of DUXAP8 and DUXAP9. Expression of DUXAP8, DUXAP9, COL1A1 and COL1A2 were significantly increased in cancer samples compared to normal controls while miR-29c-3p expression was decreased. Luciferase reporter assay revealed that miR-29c-3p could directly bind to DUXAP8, DUXAP9, COL1A1 and COL1A2. Functional experiments showed that DUXAP8 and DUXAP9 enhanced but miR-29c-3p weakened growth of renal cell carcinoma.

CONCLUSIONS

In conclusion, upregulated DUXAP8 and DUXAP9 promote growth of renal cell carcinoma and serve as two promising prognostic biomarkers.

METHODS

Dysregulated pseudogenes were obtained by dreamBase and GEPIA. The binding miRNAs of pseudogene and targets of miRNA were predicted using starBase and miRNet. Kaplan-Meier plotter was utilized to perform survival analysis, and Enrichr database was introduced to conduct functional enrichment analysis. Hub genes were identified through STRING and Cytoscape. qRT-PCR, luciferase reporter assay, cell counting assay and colony formation assay were performed to validate in silico analytic results.

dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease

Although thousands of pseudogenes have been annotated in the human genome, their transcriptional regulation, expression profiles and functional mechanisms are largely unknown. In this study, we developed dreamBase (http://rna.sysu.edu.cn/dreamBase) to facilitate the investigation of DNA modification, RNA regulation and protein binding of potential expressed pseudogenes from multidimensional high-throughput sequencing data. Based on ∼5500 ChIP-seq and DNase-seq datasets, we identified genome-wide binding profiles of various transcription-associated factors around pseudogene loci. By integrating ∼18 000 RNA-seq data, we analysed the expression profiles of pseudogenes and explored their co-expression patterns with their parent genes in 32 cancers and 31 normal tissues. By combining microRNA binding sites, we demonstrated complex post-transcriptional regulation networks involving 275 microRNAs and 1201 pseudogenes. We generated ceRNA networks to illustrate the crosstalk between pseudogenes and their parent genes through competitive binding of microRNAs. In addition, we studied transcriptome-wide interactions between RNA binding proteins (RBPs) and pseudogenes based on 458 CLIP-seq datasets. In conjunction with epitranscriptome sequencing data, we also mapped 1039 RNA modification sites onto 635 pseudogenes. This database will provide insights into the transcriptional regulation, expression, functions and mechanisms of pseudogenes as well as their roles in biological processes and diseases.

LncRNA PCNAP1 modulates hepatitis B virus replication and enhances tumor growth of liver cancer

Rationale: Hepatitis B virus (HBV) is a major risk factor for liver cancer, in which HBV covalently closed circular DNA (cccDNA) plays crucial roles. However, the effect of pseudogene-derived long noncoding RNAs (lncRNAs) acting as functional regulators of their ancestral gene expression on HBV replication and hepatocellular carcinoma (HCC) remains unclear. In this study, we speculated that the pseudogene-derived lncRNA PCNAP1 and its ancestor PCNA might modulate HBV replication and promote hepatocarcinogenesis.

Methods: We investigated the roles of lncRNA PCNAP1 in contribution of HBV replication through modulating miR-154/PCNA/HBV cccDNA signaling in hepatocarcinogenesis by using CRISPR/Cas9, Southern blot analysis, confocal assays, et al. in primary human hepatocytes (PHH), HepaRG cells, HepG2-NTCP cells, hepatoma carcinoma cells, human liver-chimeric mice model, transgenetic mice model, in vitro tumorigenicity and clinical patients.

Results: Interestingly, the expression levels of PCNAP1 and PCNA were significantly elevated in the liver of HBV-infectious human liver-chimeric mice. Clinically, the mRNA levels of PCNAP1 and PCNA were increased in the liver of HBV-positive/HBV cccDNA-positive HCC patients. Mechanistically, PCNA interacted with HBV cccDNA in a HBc-dependent manner. PCNAP1 enhanced PCNA through sponging miR-154 targeting PCNA mRNA 3′UTR. Functionally, PCNAP1 or PCNA remarkably enhanced HBV replication and accelerated the growth of HCC in vitro and in vivo.

Conclusion: We conclude that lncRNA PCNAP1 enhances the HBV replication through modulating miR-154/PCNA/HBV cccDNA signaling and the PCNAP1/PCNA signaling drives the hepatocarcinogenesis. Our finding provides new insights into the mechanism by which lncRNA PCNAP1 enhances HBV replication and hepatocarcinogenesis.

High Expression of Pseudogene PTTG3P Indicates a Poor Prognosis in Human Breast Cancer

Pseudogenes play pivotal roles in tumorigenesis. Previous studies have suggested that pituitary tumor-transforming 3, pseudogene (PTTG3P), serves as an oncogene in human cancers. However, its expression pattern, biological function, and underlying mechanism in breast cancer remain unknown. In this study, we demonstrated an elevated expression of PTTG3P in breast cancer and discovered that PTTG3P expression correlated negatively with estrogen receptor (ER) and progesterone receptor (PR) status, but linked positively to basal-like status, triple-negative breast cancer status, Nottingham prognostic index (NPI), and Scarff-Bloom-Richardson grade. High expression of PTTG3P was also found to be associated with a poor prognosis of breast cancer. To explore the potential mechanisms of PTTG3P, a PTTG3P-microRNA (miRNA)-mRNA regulatory network was established. Co-expressed genes of PTTG3P were also obtained. Enrichment analysis for these co-expressed genes revealed that they were significantly enriched in mitotic nuclear division and cell cycle. Subsequent research on mechanism of PTTG3P indicated that its expression correlated positively with PTTG1 expression. However, no significant expression correlation between PTTG3P and PTTG2 was observed. Taken together, our findings suggest that increased expression of pseudogene PTTG3P may be used as a promising prognostic biomarker and novel therapeutic target for breast cancer.

Integrated analysis of pseudogene RP11-564D11.3 expression and its potential roles in hepatocellular carcinoma

AIM

We aim to identify differentially expressed pseudogenes and investigate their functional roles in carcinogenesis.

MATERIALS & METHODS

Here, we identify dysregulated pseudogenes, analyze their prognostic values and investigate their potential functions through pseudogene-miRNA-mRNA network from public -omics repositories.

RESULTS

We identified 16 frequently upregulated pseudogenes among which high expression levels of RP11-564D11.3 was significantly associated with poor overall survival in seven types of human cancers. RP11-564D11.3 was strongly correlated with pathways in cancer, PI3K-Akt signaling pathway and the neurotrophin signaling pathway. Further studies revealed that RP11-564D11.3 functions as a competitive endogenous RNA through targeting VEGFA in hepatocellular carcinoma.

CONCLUSION

Our findings suggest RP11-564D11.3 as a novel biomarker and therapeutic potential target against hepatocellular carcinoma.

Analysis of the relationship between interleukin polymorphisms within miRNA-binding regions and alcoholic liver disease

INTRODUCTION

Alcohol consumption promotes inflammation through the Toll-like receptor 4 (TLR4)/nuclear factor (NF)-?B pathway, leading to organic damage. Some micro-RNA (miRNA) molecules modulate this inflammatory response by downregulating TLR4/NF-?B pathway mediators, like interleukins (ILs). Thus, polymorphisms within IL genes located near miRNA binding sites could modify the risk of ethanol-induced damage. The present study analyzed potential relationships between alcoholism or alcoholic liver disease (ALD) and IL12B 2124 G>T (rs1368439), IL16 5000 C>T (rs1131445), IL1R1 3114 C>T (rs3917328), and NFKB1 3400 A>G (rs4648143) polymorphisms.

PATIENTS AND METHODS

The study included 301 male alcoholic patients and 156 male healthy volunteers. Polymorphisms were genotyped using TaqMan^® PCR assays for allelic discrimination. Allele and genotype frequencies were compared between groups. Logistic regression analysis was performed to analyze the inheritance model.

RESULTS

Analysis of the IL1R1 (rs3917328) polymorphism showed that the proportion of alleleT carriers (CT and TT genotypes) was higher in healthy controls (9.7%) than in alcoholic patients (6.5%; P=.042). However, multivariable logistic regression analyses did not yield a significant result. No differences between groups were found for other analyzed polymorphisms.

CONCLUSIONS

Our study describes, for the first time, the expected frequencies of certain polymorphisms within miRNA-binding sites in alcoholic patients with and without ALD. Further studies should be developed to clarify the potential relevance of these polymorphisms in alcoholism and ALD development.

Landscape and variation of novel retroduplications in 26 human populations

Retroduplications come from reverse transcription of mRNAs and their insertion back into the genome. Here, we performed comprehensive discovery and analysis of retroduplications in a large cohort of 2,535 individuals from 26 human populations, as part of 1000 Genomes Phase 3. We developed an integrated approach to discover novel retroduplications combining high-coverage exome and low-coverage whole-genome sequencing data, utilizing information from both exon-exon junctions and discordant paired-end reads. We found 503 parent genes having novel retroduplications absent from the reference genome. Based solely on retroduplication variation, we built phylogenetic trees of human populations; these represent superpopulation structure well and indicate that variable retroduplications are effective population markers. We further identified 43 retroduplication parent genes differentiating superpopulations. This group contains several interesting insertion events, including a SLMO2 retroduplication and insertion into CAV3, which has a potential disease association. We also found retroduplications to be associated with a variety of genomic features: (1) Insertion sites were correlated with regular nucleosome positioning. (2) They, predictably, tend to avoid conserved functional regions, such as exons, but, somewhat surprisingly, also avoid introns. (3) Retroduplications tend to be co-inserted with young L1 elements, indicating recent retrotranspositional activity, and (4) they have a weak tendency to originate from highly expressed parent genes. Our investigation provides insight into the functional impact and association with genomic elements of retroduplications. We anticipate our approach and analytical methodology to have application in a more clinical context, where exome sequencing data is abundant and the discovery of retroduplications can potentially improve the accuracy of SNP calling.

Pseudogenes and Their Genome-Wide Prediction in Plants

Pseudogenes are paralogs generated from ancestral functional genes (parents) during genome evolution, which contain critical defects in their sequences, such as lacking a promoter, having a premature stop codon or frameshift mutations. Generally, pseudogenes are functionless, but recent evidence demonstrates that some of them have potential roles in regulation. The majority of pseudogenes are generated from functional progenitor genes either by gene duplication (duplicated pseudogenes) or retro-transposition (processed pseudogenes). Pseudogenes are primarily identified by comparison to their parent genes. Bioinformatics tools for pseudogene prediction have been developed, among which PseudoPipe, PSF and Shiu's pipeline are publicly available. We compared these three tools using the well-annotated Arabidopsis thaliana genome and its known 924 pseudogenes as a test data set. PseudoPipe and Shiu's pipeline identified ~80% of A. thaliana pseudogenes, of which 94% were shared, while PSF failed to generate adequate results. A need for improvement of the bioinformatics tools for pseudogene prediction accuracy in plant genomes was thus identified, with the ultimate goal of improving the quality of genome annotation in plants.

Two sets of RNAi components are required for heterochromatin formation in trans triggered by truncated transgenes

Across kingdoms, RNA interference (RNAi) has been shown to control gene expression at the transcriptional- or the post-transcriptional level. Here, we describe a mechanism which involves both aspects: truncated transgenes, which fail to produce intact mRNA, induce siRNA accumulation and silencing of homologous loci in trans in the ciliate Paramecium. We show that silencing is achieved by co-transcriptional silencing, associated with repressive histone marks at the endogenous gene. This is accompanied by secondary siRNA accumulation, strictly limited to the open reading frame of the remote locus. Our data shows that in this mechanism, heterochromatic marks depend on a variety of RNAi components. These include RDR3 and PTIWI14 as well as a second set of components, which are also involved in post-transcriptional silencing: RDR2, PTIWI13, DCR1 and CID2. Our data indicates differential processing of nascent un-spliced and long, spliced transcripts thus suggesting a hitherto-unrecognized functional interaction between post-transcriptional and co-transcriptional RNAi. Both sets of RNAi components are required for efficient trans-acting RNAi at the chromatin level and our data indicates similar mechanisms contributing to genome wide regulation of gene expression by epigenetic mechanisms.

Role of non-coding sequence variants in cancer

Patients with cancer carry somatic sequence variants in their tumour in addition to the germline variants in their inherited genome. Although variants in protein-coding regions have received the most attention, numerous studies have noted the importance of non-coding variants in cancer. Moreover, the overwhelming majority of variants, both somatic and germline, occur in non-coding portions of the genome. We review the current understanding of non-coding variants in cancer, including the great diversity of the mutation types--from single nucleotide variants to large genomic rearrangements--and the wide range of mechanisms by which they affect gene expression to promote tumorigenesis, such as disrupting transcription factor-binding sites or functions of non-coding RNAs. We highlight specific case studies of somatic and germline variants, and discuss how non-coding variants can be interpreted on a large-scale through computational and experimental methods.

Tissue-specific expression of histone H3 variants diversified after species separation

Background

The selective incorporation of appropriate histone variants into chromatin is critical for the regulation of genome function. Although many histone variants have been identified, a complete list has not been compiled.

Results

We screened mouse, rat and human genomes by in silico hybridization using canonical histone sequences. In the mouse genome, we identified 14 uncharacterized H3 genes, among which 13 are similar to H3.3 and do not have human or rat counterparts, and one is similar to human testis-specific H3 variant, H3T/H3.4, and had a rat paralog. Although some of these genes were previously annotated as pseudogenes, their tissue-specific expression was confirmed by sequencing the 3′-UTR regions of the transcripts. Certain new variants were also detected at the protein level by mass spectrometry. When expressed as GFP-tagged versions in mouse C2C12 cells, some variants were stably incorporated into chromatin and the genome-wide distributions of most variants were similar to that of H3.3. Moreover, forced expression of H3 variants in chromatin resulted in alternate gene expression patterns after cell differentiation.

Conclusions

We comprehensively identified and characterized novel mouse H3 variant genes that encoded highly conserved amino acid sequences compared to known histone H3. We speculated that the diversity of H3 variants acquired after species separation played a role in regulating tissue-specific gene expression in individual species. Their biological relevance and evolutionary aspect involving pseudogene diversification will be addressed by further functional analysis.

Electronic supplementary material

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

A nuclear ribosomal DNA pseudogene in triatomines opens a new research field of fundamental and applied implications in Chagas disease

A pseudogene, designated as "ps(5.8S+ITS-2)", paralogous to the 5.8S gene and internal transcribed spacer (ITS)-2 of the nuclear ribosomal DNA (rDNA), has been recently found in many triatomine species distributed throughout North America, Central America and northern South America. Among characteristics used as criteria for pseudogene verification, secondary structures and free energy are highlighted, showing a lower fit between minimum free energy, partition function and centroid structures, although in given cases the fit only appeared to be slightly lower. The unique characteristics of "ps(5.8S+ITS-2)" as a processed or retrotransposed pseudogenic unit of the ghost type are reviewed, with emphasis on its potential functionality compared to the functionality of genes and spacers of the normal rDNA operon. Besides the technical problem of the risk for erroneous sequence results, the usefulness of "ps(5.8S+ITS-2)" for specimen classification, phylogenetic analyses and systematic/taxonomic studies should be highlighted, based on consistence and retention index values, which in pseudogenic sequence trees were higher than in functional sequence trees. Additionally, intraindividual, interpopulational and interspecific differences in pseudogene amount and the fact that it is a pseudogene in the nuclear rDNA suggests a potential relationships with fitness, behaviour and adaptability of triatomine vectors and consequently its potential utility in Chagas disease epidemiology and control.

Taking the pseudo out of pseudogenes

Pseudogenes are defined as fragments of once-functional genes that have been silenced by one or more nonsense, frameshift or missense mutations. Despite continuing increases in the speed of sequencing and annotating bacterial genomes, the identification and categorisation of pseudogenes remains problematic. Even when identified, pseudogenes are considered to be rare and tend to be ignored. On the contrary, pseudogenes are surprisingly prevalent and can persist for long evolutionary time periods, representing a record of once-functional genetic characteristics. Most importantly, pseudogenes provide an insight into prokaryotic evolutionary history as a record of phenotypic traits that have been lost. Focusing on the intracellular and symbiotic bacteria in which pseudogenes predominate, this review discusses the importance of identifying pseudogenes to fully understand the abilities of bacteria, and to understand prokaryotes within their evolutionary context.

Comparative analysis of pseudogenes across three phyla

Pseudogenes are degraded fossil copies of genes. Here, we report a comparison of pseudogenes spanning three phyla, leveraging the completed annotations of the human, worm, and fly genomes, which we make available as an online resource. We find that pseudogenes are lineage specific, much more so than protein-coding genes, reflecting the different remodeling processes marking each organism's genome evolution. The majority of human pseudogenes are processed, resulting from a retrotranspositional burst at the dawn of the primate lineage. This burst can be seen in the largely uniform distribution of pseudogenes across the genome, their preservation in areas with low recombination rates, and their preponderance in highly expressed gene families. In contrast, worm and fly pseudogenes tell a story of numerous duplication events. In worm, these duplications have been preserved through selective sweeps, so we see a large number of pseudogenes associated with highly duplicated families such as chemoreceptors. However, in fly, the large effective population size and high deletion rate resulted in a depletion of the pseudogene complement. Despite large variations between these species, we also find notable similarities. Overall, we identify a broad spectrum of biochemical activity for pseudogenes, with the majority in each organism exhibiting varying degrees of partial activity. In particular, we identify a consistent amount of transcription (∼15%) across all species, suggesting a uniform degradation process. Also, we see a uniform decay of pseudogene promoter activity relative to their coding counterparts and identify a number of pseudogenes with conserved upstream sequences and activity, hinting at potential regulatory roles.

Identification of maize long non-coding RNAs responsive to drought stress

Long non-coding RNAs (lncRNAs) represent a class of riboregulators that either directly act in long form or are processed to shorter miRNAs and siRNAs. Emerging evidence shows that lncRNAs participate in stress responsive regulation. In this study, to identify the putative maize lncRNAs responsive to drought stress, 8449 drought responsive transcripts were first uploaded to the Coding Potential Calculator website for classification as protein coding or non-coding RNAs, and 1724 RNAs were identified as potential non-coding RNAs. A Perl script was written to screen these 1724 ncRNAs and 664 transcripts were ultimately identified as drought-responsive lncRNAs. Of these 664 transcripts, 126 drought-responsive lncRNAs were highly similar to known maize lncRNAs; the remaining 538 transcripts were considered as novel lncRNAs. Among the 664 lncRNAs identified as drought responsive, 567 were upregulated and 97 were downregulated in drought-stressed leaves of maize. 8 lncRNAs were identified as miRNA precursor lncRNAs, 62 were classified as both shRNA and siRNA precursors, and 279 were classified as siRNA precursors. The remaining 315 lncRNAs were classified as other lncRNAs that are likely to function as longer molecules. Among these 315 lncRNAs, 10 are identified as antisense lncRNAs and 7 could pair with 17 CDS sequences with near-perfect matches. Finally, RT-qPCR results confirmed that all selected lncRNAs could respond to drought stress. These findings extend the current view on lncRNAs as ubiquitous regulators under stress conditions.

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.

Faster evolving Drosophila paralogs lose expression rate and ubiquity and accumulate more non-synonymous SNPs

Background

Duplicated genes can indefinately persist in genomes if either both copies retain the original function due to dosage benefit (gene conservation), or one of the copies assumes a novel function (neofunctionalization), or both copies become required to perform the function previously accomplished by a single copy (subfunctionalization), or through a combination of these mechanisms. Different models of duplication retention imply different predictions about substitution rates in the coding portion of paralogs and about asymmetry of these rates.

Results

We analyse sequence evolution asymmetry in paralogs present in 12 Drosophila genomes using the nearest non-duplicated orthologous outgroup as a reference. Those paralogs present in D. melanogaster are analysed in conjunction with the asymmetry of expression rate and ubiquity and of segregating non-synonymous polymorphisms in the same paralogs. Paralogs accumulate substitutions, on average, faster than their nearest singleton orthologs. The distribution of paralogs’ substitution rate asymmetry is overdispersed relative to that of orthologous clades, containing disproportionally more unusually symmetric and unusually asymmetric clades. We show that paralogs are more asymmetric in: a) clades orthologous to highly constrained singleton genes; b) genes with high expression level; c) genes with ubiquitous expression and d) non-tandem duplications. We further demonstrate that, in each asymmetrically evolving pair of paralogs, the faster evolving member of the pair tends to have lower average expression rate, lower expression uniformity and higher frequency of non-synonymous SNPs than its slower evolving counterpart.

Conclusions

Our findings are consistent with the hypothesis that many duplications in Drosophila are retained despite stabilising selection being more relaxed in one of the paralogs than in the other, suggesting a widespread unfinished pseudogenization. This phenomenon is likely to make detection of neo- and subfunctionalization signatures difficult, as these models of duplication retention also predict asymmetries in substitution rates and expression profiles.

Reviewers

This article has been reviewed by Dr. Jia Zeng (nominated by Dr. I. King Jordan), Dr. Fyodor Kondrashov and Dr. Yuri Wolf.

Phylogenetic analysis of the expansion of the MATH-BTB gene family in the grasses

MATH-BTB proteins are known to act as substrate-specific adaptors of cullin3 (CUL3)-based ubiquitin E3 ligases to target protein for ubiquitination. In a previous study we reported the presence of 31 MATH-BTB genes in the maize genome and determined the regulatory role of the MATH-BTB protein MAB1 during meiosis to mitosis transition. In contrast to maize, there are only 6 homologous genes in the model plant Arabidopsis, while this family has largely expanded in grasses. Here, we report a phylogenetic analysis of the MATH-BTB gene family in 9 land plant species including various mosses, eudicots, and grasses. We extend a previous classification of the plant MATH-BTB family and additionally arrange the expanded group into 5 grass-specific clades. Synteny studies indicate that expansion occurred to a large extent due to local gene duplications. Expression studies of 3 closely related MATH-BTB genes in maize (MAB1-3) indicate highly specific expression pattern. In summary, this work provides a solid base for further studies comparing genetic and functional information of the MATH-BTB family especially in the grasses.

High-throughput RNA sequencing of pseudomonas-infected Arabidopsis reveals hidden transcriptome complexity and novel splice variants

We report the results of a genome-wide analysis of transcription in Arabidopsis thaliana after treatment with Pseudomonas syringae pathovar tomato. Our time course RNA-Seq experiment uses over 500 million read pairs to provide a detailed characterization of the response to infection in both susceptible and resistant hosts. The set of observed differentially expressed genes is consistent with previous studies, confirming and extending existing findings about genes likely to play an important role in the defense response to Pseudomonas syringae. The high coverage of the Arabidopsis transcriptome resulted in the discovery of a surprisingly large number of alternative splicing (AS) events--more than 44% of multi-exon genes showed evidence for novel AS in at least one of the probed conditions. This demonstrates that the Arabidopsis transcriptome annotation is still highly incomplete, and that AS events are more abundant than expected. To further refine our predictions, we identified genes with statistically significant changes in the ratios of alternative isoforms between treatments. This set includes several genes previously known to be alternatively spliced or expressed during the defense response, and it may serve as a pool of candidate genes for regulated alternative splicing with possible biological relevance for the defense response against invasive pathogens.

Analysis of variable retroduplications in human populations suggests coupling of retrotransposition to cell division

In primates and other animals, reverse transcription of mRNA followed by genomic integration creates retroduplications. Expressed retroduplications are either “retrogenes” coding for functioning proteins, or expressed “processed pseudogenes,” which can function as noncoding RNAs. To date, little is known about the variation in retroduplications in terms of their presence or absence across individuals in the human population. We have developed new methodologies that allow us to identify “novel” retroduplications (i.e., those not present in the reference genome), to find their insertion points, and to genotype them. Using these methods, we catalogued and analyzed 174 retroduplication variants in almost one thousand humans, which were sequenced as part of Phase 1 of The 1000 Genomes Project Consortium. The accuracy of our data set was corroborated by (1) multiple lines of sequencing evidence for retroduplication (e.g., depth of coverage in exons vs. introns), (2) experimental validation, and (3) the fact that we can reconstruct a correct phylogenetic tree of human subpopulations based solely on retroduplications. We also show that parent genes of retroduplication variants tend to be expressed at the M-to-G1 transition in the cell cycle and that M-to-G1 expressed genes have more copies of fixed retroduplications than genes expressed at other times. These findings suggest that cell division is coupled to retrotransposition and, perhaps, is even a requirement for it.

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.

Human Transcriptome and Chromatin Modifications: An ENCODE Perspective

A decade-long project, led by several international research groups, called the Encyclopedia of DNA Elements (ENCODE), recently released an unprecedented amount of data. The ambitious project covers transcriptome, cistrome, epigenome, and interactome data from more than 1,600 sets of experiments in human. To make use of this valuable resource, it is important to understand the information it represents and the techniques that were used to generate these data. In this review, we introduce the data that ENCODE generated, summarize the observations from the data analysis, and revisit a computational approach that ENCODE used to predict gene expression, with a focus on the human transcriptome and its association with chromatin modifications.

Vascular calcification in diabetes: mechanisms and implications

Cardiovascular disease (CVD) remains the leading cause of death among adults with diabetes, and CVD prevention remains a major challenge. Coronary artery calcium (CAC) score measured by electron beam tomography (EBT) or multi-slice detector computed tomography correlates closely with plaque burden and coronary angiography, and predicts coronary events independently of other risk factors. Further, progression of CAC over several years has been shown to predict increased mortality. Coronary calcification is an active process strongly associated with atherosclerotic plaque evolution and is an accepted surrogate endpoint in studies of patients with diabetes older than 30. In this review, recent findings regarding the mechanisms and implications of vascular calcification in diabetes will be discussed.

Sex bias in copy number variation of olfactory receptor gene family depends on ethnicity

Gender plays a pivotal role in the human genetic identity and is also manifested in many genetic disorders particularly mental retardation. In this study its effect on copy number variation (CNV), known to cause genetic disorders was explored. As the olfactory receptor (OR) repertoire comprises the largest human gene family, it was selected for this study, which was carried out within and between three populations, derived from 150 individuals from the 1000 Genome Project. Analysis of 3872 CNVs detected among 791 OR loci, in which 307 loci showed CNV, revealed the following novel findings: Sex bias in CNV was significantly more prevalent in uncommon than common CNV variants of OR pseudogenes, in which the male genome showed more CNVs; and in one-copy number loss compared to complete deletion of OR pseudogenes; both findings implying a more recent evolutionary role for gender. Sex bias in copy number gain was also detected. Another novel finding was that the observed sex bias was largely dependent on ethnicity and was in general absent in East Asians. Using a CNV public database for sick children (International Standard Cytogenomic Array Consortium) the application of these findings for improving clinical molecular diagnostics is discussed by showing an example of sex bias in CNV among kids with autism. Additional clinical relevance is discussed, as the most polymorphic CNV-enriched OR cluster in the human genome, located on chr 15q11.2, is found near the Prader–Willi syndrome/Angelman syndrome bi-directionally imprinted region associated with two well-known mental retardation syndromes. As olfaction represents the primitive cognition in most mammals, arguably in competition with the development of a larger brain, the extensive retention of OR pseudogenes in females of this study, might point to a parent-of-origin indirect regulatory role for OR pseudogenes in the embryonic development of human brain. Thus any perturbation in the temporal regulation of olfactory system could lead to developmental delay disorders including mental retardation.

The GENCODE pseudogene resource

Background

Pseudogenes have long been considered as nonfunctional genomic sequences. However, recent evidence suggests that many of them might have some form of biological activity, and the possibility of functionality has increased interest in their accurate annotation and integration with functional genomics data.

Results

As part of the GENCODE annotation of the human genome, we present the first genome-wide pseudogene assignment for protein-coding genes, based on both large-scale manual annotation and in silico pipelines. A key aspect of this coupled approach is that it allows us to identify pseudogenes in an unbiased fashion as well as untangle complex events through manual evaluation. We integrate the pseudogene annotations with the extensive ENCODE functional genomics information. In particular, we determine the expression level, transcription-factor and RNA polymerase II binding, and chromatin marks associated with each pseudogene. Based on their distribution, we develop simple statistical models for each type of activity, which we validate with large-scale RT-PCR-Seq experiments. Finally, we compare our pseudogenes with conservation and variation data from primate alignments and the 1000 Genomes project, producing lists of pseudogenes potentially under selection.

Conclusions

At one extreme, some pseudogenes possess conventional characteristics of functionality; these may represent genes that have recently died. On the other hand, we find interesting patterns of partial activity, which may suggest that dead genes are being resurrected as functioning non-coding RNAs. The activity data of each pseudogene are stored in an associated resource, psiDR, which will be useful for the initial identification of potentially functional pseudogenes.

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.

In with the old, in with the new: the promiscuity of the duplication process engenders diverse pathways for novel gene creation

The gene duplication process has exhibited far greater promiscuity in the creation of paralogs with novel exon-intron structures than anticipated even by Ohno. In this paper I explore the history of the field, from the neo-Darwinian synthesis through Ohno's formulation of the canonical model for the evolution of gene duplicates and culminating in the present genomic era. I delineate the major tenets of Ohno's model and discuss its failure to encapsulate the full complexity of the duplication process as revealed in the era of genomics. I discuss the diverse classes of paralogs originating from both DNA- and RNA-mediated duplication events and their evolutionary potential for assuming radically altered functions, as well as the degree to which they can function unconstrained from the pressure of gene conversion. Lastly, I explore theoretical population-genetic considerations of how the effective population size (N(e)) of a species may influence the probability of emergence of genes with radically altered functions.

Expressed pseudogenes in the transcriptional landscape of human cancers

Pseudogene transcripts can provide a novel tier of gene regulation through generation of endogenous siRNAs or miRNA-binding sites. Characterization of pseudogene expression, however, has remained confined to anecdotal observations due to analytical challenges posed by the extremely close sequence similarity with their counterpart coding genes. Here, we describe a systematic analysis of pseudogene "transcription" from an RNA-Seq resource of 293 samples, representing 13 cancer and normal tissue types, and observe a surprisingly prevalent, genome-wide expression of pseudogenes that could be categorized as ubiquitously expressed or lineage and/or cancer specific. Further, we explore disease subtype specificity and functions of selected expressed pseudogenes. Taken together, we provide evidence that transcribed pseudogenes are a significant contributor to the transcriptional landscape of cells and are positioned to play significant roles in cellular differentiation and cancer progression, especially in light of the recently described ceRNA networks. Our work provides a transcriptome resource that enables high-throughput analyses of pseudogene expression.