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

Hypothalamic and hippocampal transcriptome changes in AppNL-G-F mice as a function of metabolic and inflammatory dysfunction

The progression of Alzheimer's disease (AD) has a silent phase that predates characteristic cognitive decline and eventually leads to active cognitive deficits. Metabolism, diet, and obesity have been correlated to the development of AD but is poorly understood. The hypothalamus is a brain region that exerts homeostatic control on food intake and metabolism and has been noted to be impacted during the active phase of Alzheimer's disease. This study, in using an amyloid overexpression AppNL-G-F mouse model under normal metabolic conditions, examines blood markers in young and old male AppNL-G-F mice (n = 5) that corresponds to the silent and active phases of AD, and bulk gene expression changes in the hypothalamus and the hippocampus. The results show a large panel of inflammatory mediators, leptin, and other proteins that may be involved in weakening the blood brain barrier, to be increased in the young AppNL-G-F mice but not in the old AppNL-G-F mice. There were also several differentially expressed genes in both the hypothalamus and the hippocampus in the young AppNL-G-F mice prior to amyloid plaque formation and cognitive decline that persisted in the old AppNL-G-F mice, including GABRa2 receptor, Wdfy1, and several pseudogenes with unknown function. These results suggests that a larger panel of inflammatory mediators may be used as blood markers to detect silent AD, and that a change in leptin and gene expression in the hypothalamus exist prior to cognitive effects, suggesting a coupling of metabolism with amyloid plaque induced cognitive decline.

Competing endogenous RNA (ceRNA) networks in Parkinson's disease: A systematic review

Parkinson's disease (PD) is a distinctive clinical syndrome with several causes and clinical manifestations. Aside from an infectious cause, PD is a rapidly developing neurological disorder with a global rise in frequency. Notably, improved knowledge of molecular pathways and the developing novel diagnostic methods may result in better therapy for PD patients. In this regard, the amount of research on ceRNA axes is rising, highlighting the importance of these axes in PD. CeRNAs are transcripts that cross-regulate one another via competition for shared microRNAs (miRNAs). These transcripts may be either coding RNAs (mRNAs) or non-coding RNAs (ncRNAs). This research used a systematic review to assess validated loops of ceRNA in PD. The Prisma guideline was used to conduct this systematic review, which entailed systematically examining the articles of seven databases. Out of 309 entries, forty articles met all criteria for inclusion and were summarized in the appropriate table. CeRNA axes have been described through one of the shared vital components of the axes, including lncRNAs such as NEAT1, SNHG family, HOTAIR, MALAT1, XIST, circRNAs, and lincRNAs. Understanding the multiple aspects of this regulatory structure may aid in elucidating the unknown causal causes of PD and providing innovative molecular therapeutic targets and medical fields.

lncRNA-mediated ceRNA network in bladder cancer

Bladder cancer is a common disease associated with high rates of morbidity and mortality. Although immunotherapy approaches such as adoptive T-cell therapy and immune checkpoint blockade have been investigated for the treatment of bladder cancer, their off-target effects and ability to affect only single targets have led to clinical outcomes that are far from satisfactory. Therefore, it is important to identify novel targets that can effectively control tumor growth and metastasis. It is well known that long noncoding RNAs (lncRNAs) are powerful regulators of gene expression. Increasing evidence has shown that dysregulated lncRNAs in bladder cancer are involved in cancer cell proliferation, migration, invasion, apoptosis, and epithelial-mesenchymal transition (EMT). In this review, we focus on the roles and underlying mechanisms of lncRNA-mediated competing endogenous RNA (ceRNA) networks in the regulation of bladder cancer progression. In addition, we discuss the potential of targeting lncRNA-mediated ceRNA networks to overcome cancer treatment resistance and its association with clinicopathological features and outcomes in bladder cancer patients. We hope this review will stimulate research to develop more effective therapeutic approaches for bladder cancer treatment.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

Pseudogene AK4P1 promotes pancreatic ductal adenocarcinoma progression through relieving miR-375-mediated YAP1 degradation

Pseudogenes have been reported to play oncogenic or tumor-suppressive roles in cancer progression. However, the molecular mechanism of most pseudogenes in pancreatic ductal adenocarcinoma (PDAC) remains unknown. Herein, we characterized a novel pseudogene-miRNA-mRNA network associated with PDAC progression using bioinformatics analysis. After screening by dreamBase and GEPIA, 12 up-regulated and 7 down-regulated differentially expressed pseudogenes (DEPs) were identified. According to survival analysis, only elevated AK4P1 indicated a poor prognosis for PDAC patients. Moreover, we found that AK4 acts as a cognate gene of AK4P1 and also predicts worse survival for PDAC patients. Furthermore, 32 miRNAs were predicted to bind to AK4P1 by starBase, among which miR-375 was identified as the most potential binding miRNA of AK4P1. A total of 477 potential target genes of miR-375 were obtained by miRNet, in which 49 hub genes with node degree ≥ 20 were identified by STRING. Subsequent analysis for hub genes demonstrated that YAP1 may be a functional downstream target of AK4P1. To confirmed the above findings, microarray, and qRT-PCR assay revealed that YAP1 was dramatically upregulated in both PDAC cells and tissues. Functional experiments showed that knockdown of YAP1 significantly suppressed PDAC cells growth, increased apoptosis, and decreased the ability of invasion. In conclusion, amplification of AK4P1 may fuel the onset and development of PDAC by targeting YAP1 through competitively binding to miR-375, and serve as a promising biomarker and therapeutic target for PDAC.

Virtual Gene Concept and a Corresponding Pragmatic Research Program in Genetical Data Science

Mendel proposed an experimentally verifiable paradigm of particle-based heredity that has been influential for over 150 years. The historical arguments have been reflected in the near past as Mendel's concept has been diversified by new types of omics data. As an effect of the accumulation of omics data, a virtual gene concept forms, giving rise to genetical data science. The concept integrates genetical, functional, and molecular features of the Mendelian paradigm. I argue that the virtual gene concept should be deployed pragmatically. Indeed, the concept has already inspired a practical research program related to systems genetics. The program includes questions about functionality of structural and categorical gene variants, about regulation of gene expression, and about roles of epigenetic modifications. The methodology of the program includes bioinformatics, machine learning, and deep learning. Education, funding, careers, standards, benchmarks, and tools to monitor research progress should be provided to support the research program.

Processed pseudogenes: A substrate for evolutionary innovation: Retrotransposition contributes to genome evolution by propagating pseudogene sequences with rich regulatory potential throughout the genome

Processed pseudogenes may serve as a genetic reservoir for evolutionary innovation. Here, we argue that through the activity of long interspersed element-1 retrotransposons, processed pseudogenes disperse coding and noncoding sequences rich with regulatory potential throughout the human genome. While these sequences may appear to be non-functional, a lack of contemporary function does not prohibit future development of biological activity. Here, we discuss the dynamic evolution of certain processed pseudogenes into coding and noncoding genes and regulatory elements, and their implication in wide-ranging biological and pathological processes. Also see the video abstract here: https://youtu.be/iUY_mteVoPI.

Genetic variations of DNA bindings of FOXA1 and co-factors in breast cancer susceptibility

Identifying transcription factors (TFs) whose DNA bindings are altered by genetic variants that regulate susceptibility genes is imperative to understand transcriptional dysregulation in disease etiology. Here, we develop a statistical framework to analyze extensive ChIP-seq and GWAS data and identify 22 breast cancer risk-associated TFs. We find that, by analyzing genetic variations of TF-DNA bindings, the interaction of FOXA1 with co-factors such as ESR1 and E2F1, and the interaction of TFs with chromatin features (i.e., enhancers) play a key role in breast cancer susceptibility. Using genetic variants occupied by the 22 TFs, transcriptome-wide association analyses identify 52 previously unreported breast cancer susceptibility genes, including seven with evidence of essentiality from functional screens in breast relevant cell lines. We show that FOXA1 and co-factors form a core TF-transcriptional network regulating the susceptibility genes. Our findings provide additional insights into genetic variations of TF-DNA bindings (particularly for FOXA1) underlying breast cancer susceptibility.

The identification of transcription factors (TFs) whose binding sites are affected by risk genetic variants remains crucial. Here, the authors develop a statistical framework to analyse ChIP-seq and GWAS data, identify 22 breast cancer risk-associated TFs and a core TF-transcriptional network for FOXA1 and co-factors.

Pseudogenes and their potential functions in hematopoiesis

Pseudogenes are DNA regions comprising defective copies of functional genes, the majority of which were generated by RNA- or DNA-level duplications. They exist across almost all forms of life and account for about one-quarter of the annotated genes in the human genome. Although these have been considered nonfunctional for decades, a growing number of pseudogenes have been found to be transcribed and to play crucial regulatory roles. Accumulating evidence indicates that they regulate gene expression through molecular interactions at the protein, RNA, and DNA levels. However, pseudogenes are often excluded in multiple genomewide analyses and functional screening, and their biological activities remain to be systematically disclosed. Here, we summarize the features of and progress of research on pseudogenes, in addition to discussing what is unknown about these genetic elements. Our previous findings, together with evidence of their poor conservation, prompted us to propose that pseudogenes may contribute to primate- or human-specific regulation, especially in hematopoiesis.

Pseudogenes: Four Decades of Discovery

A pseudogene is defined as a genomic DNA sequence that looks like a mutated or truncated version of a known functional gene. Nearly four decades since their first discovery it has been estimated that between ~12,000 and ~20,000 pseudogenes exist in the human genome. Early efforts to characterize functions for pseudogenes were unsuccessful, thus they were considered functionless relics of evolutionary selection, junk DNA or genetic fossils. Remarkably, an increasing number of pseudogenes have been reported to be expressed as RNA transcripts above and beyond levels considered accidental or spurious transcription. There is emerging evidence that some expressed pseudogene transcripts have biological functions and should be defined as a subclass of functional long noncoding RNAs (lncRNA). In this introductory chapter, I briefly summarize the history and the current knowledge of pseudogenes, and highlight the emerging functions of some pseudogenes in human biology and disease. This second iteration of Pseudogenes in Methods in Molecular Biology highlights new methodological approaches to investigate this intriguing family of lncRNAs and the extent of their biological function.

PLEKHA8P1 Promotes Tumor Progression and Indicates Poor Prognosis of Liver Cancer

Hepatocellular carcinoma (HCC) records the second-lowest 5-year survival rate despite the avalanche of research into diagnosis and therapy. One of the major obstacles in treatment is chemoresistance to drugs such as 5-fluorouracil (5-FU), making identification and elucidation of chemoresistance regulators highly valuable. As the regulatory landscape grows to encompass non-coding genes such as long non-coding RNAs (lncRNAs), a relatively new class of lncRNA has emerged in the form of pseudogene-derived lncRNAs. Through bioinformatics analyses of the TCGA LIHC dataset, we have systematically identified pseudogenes of prognostic value. Initial experimental validation of selected pseudogene-derived lncRNA (PLEKHA8P1) and its parental gene (PLEKHA8), a well-studied transport protein in Golgi complex recently implicated as an oncogene in both colorectal and liver cancer, indicates that the pseudogene/parental gene pair promotes tumor progression and that their dysregulated expression levels affect 5-FU-induced chemoresistance in human HCC cell line FT3-7. Our study has thus confirmed cancer-related functions of PLEKHA8, and laid the groundwork for identification and validation of oncogenic pseudogene-derived lncRNA that shows potential as a novel therapeutic target in circumventing chemoresistance induced by 5-FU.

Long-read cDNA sequencing identifies functional pseudogenes in the human transcriptome

Pseudogenes are gene copies presumed to mainly be functionless relics of evolution due to acquired deleterious mutations or transcriptional silencing. Using deep full-length PacBio cDNA sequencing of normal human tissues and cancer cell lines, we identify here hundreds of novel transcribed pseudogenes expressed in tissue-specific patterns. Some pseudogene transcripts have intact open reading frames and are translated in cultured cells, representing unannotated protein-coding genes. To assess the biological impact of noncoding pseudogenes, we CRISPR-Cas9 delete the nucleus-enriched pseudogene PDCL3P4 and observe hundreds of perturbed genes. This study highlights pseudogenes as a complex and dynamic component of the human transcriptional landscape.

Genome-wide analysis of pseudogenes reveals HBBP1's human-specific essentiality in erythropoiesis and implication in β-thalassemia

The human genome harbors 14,000 duplicated or retroposed pseudogenes. Given their functionality as regulatory RNAs and low conservation, we hypothesized that pseudogenes could shape human-specific phenotypes. To test this, we performed co-expression analyses and found that pseudogene exhibited tissue-specific expression, especially in the bone marrow. By incorporating genetic data, we identified a bone-marrow-specific duplicated pseudogene, HBBP1 (η-globin), which has been implicated in β-thalassemia. Extensive functional assays demonstrated that HBBP1 is essential for erythropoiesis by binding the RNA-binding protein (RBP), HNRNPA1, to upregulate TAL1, a key regulator of erythropoiesis. The HBBP1/TAL1 interaction contributes to a milder symptom in β-thalassemia patients. Comparative studies further indicated that the HBBP1/TAL1 interaction is human-specific. Genome-wide analyses showed that duplicated pseudogenes are often bound by RBPs and less commonly bound by microRNAs compared with retropseudogenes. Taken together, we not only demonstrate that pseudogenes can drive human evolution but also provide insights on their functional landscapes.

Functional annotation of genetic associations by transcriptome-wide association analysis provides insights into neutrophil development regulation

Genome-wide association studies (GWAS) have identified multiple genomic loci linked to blood cell traits, however understanding the biological relevance of these genetic loci has proven to be challenging. Here, we performed a transcriptome-wide association study (TWAS) integrating gene expression and splice junction usage in neutrophils (N = 196) with a neutrophil count GWAS (N = 173,480 individuals). We identified a total of 174 TWAS-significant genes enriched in target genes of master transcription factors governing neutrophil specification. Knockout of a TWAS candidate at chromosome 5q13.2, TAF9, in CD34⁺ hematopoietic and progenitor cells (HSPCs) using CRISPR/Cas9 technology showed a significant effect on neutrophil production in vitro. In addition, we identified 89 unique genes significant only for splice junction usage, thus emphasizing the importance of alternative splicing beyond gene expression underlying granulopoiesis. Our results highlight the advantages of TWAS, followed by gene editing, to determine the functions of GWAS loci implicated in hematopoiesis.

Competing Endogenous RNA Networks as Biomarkers in Neurodegenerative Diseases

Protein aggregation is classically considered the main cause of neuronal death in neurodegenerative diseases (NDDs). However, increasing evidence suggests that alteration of RNA metabolism is a key factor in the etiopathogenesis of these complex disorders. Non-coding RNAs are the major contributor to the human transcriptome and are particularly abundant in the central nervous system, where they have been proposed to be involved in the onset and development of NDDs. Interestingly, some ncRNAs (such as lncRNAs, circRNAs and pseudogenes) share a common functionality in their ability to regulate gene expression by modulating miRNAs in a phenomenon known as the competing endogenous RNA mechanism. Moreover, ncRNAs are found in body fluids where their presence and concentration could serve as potential non-invasive biomarkers of NDDs. In this review, we summarize the ceRNA networks described in Alzheimer's disease, Parkinson's disease, multiple sclerosis, amyotrophic lateral sclerosis and spinocerebellar ataxia type 7, and discuss their potential as biomarkers of these NDDs. Although numerous studies have been carried out, further research is needed to validate these complex interactions between RNAs and the alterations in RNA editing that could provide specific ceRNET profiles for neurodegenerative disorders, paving the way to a better understanding of these diseases.

FUS-dependent loading of SUV39H1 to OCT4 pseudogene-lncRNA programs a silencing complex with OCT4 promoter specificity

The resurrection of pseudogenes during evolution produced lncRNAs with new biological function. Here we show that pseudogene-evolution created an Oct4 pseudogene lncRNA that is able to direct epigenetic silencing of the parental Oct4 gene via a 2-step, lncRNA dependent mechanism. The murine Oct4 pseudogene 4 (mOct4P4) lncRNA recruits the RNA binding protein FUS to allow the binding of the SUV39H1 HMTase to a defined mOct4P4 lncRNA sequence element. The mOct4P4-FUS-SUV39H1 silencing complex holds target site specificity for the parental Oct4 promoter and interference with individual components results in loss of Oct4 silencing. SUV39H1 and FUS do not bind parental Oct4 mRNA, confirming the acquisition of a new biological function by the mOct4P4 lncRNA. Importantly, all features of mOct4P4 function are recapitulated by the human hOCT4P3 pseudogene lncRNA, indicating evolutionary conservation. Our data highlight the biological relevance of rapidly evolving lncRNAs that infiltrate into central epigenetic regulatory circuits in vertebrate cells.

Genetic and environmental determinants of variation in the plasma lipidome of older Australian twins

The critical role of blood lipids in a broad range of health and disease states is well recognised but less explored is the interplay of genetics and environment within the broader blood lipidome. We examined heritability of the plasma lipidome among healthy older-aged twins (75 monozygotic/55 dizygotic pairs) enrolled in the Older Australian Twins Study (OATS) and explored corresponding gene expression and DNA methylation associations. 27/209 lipids (13.3%) detected by liquid chromatography-coupled mass spectrometry (LC-MS) were significantly heritable under the classical ACE twin model (h2 = 0.28-0.59), which included ceramides (Cer) and triglycerides (TG). Relative to non-significantly heritable TGs, heritable TGs had a greater number of associations with gene transcripts, not directly associated with lipid metabolism, but with immune function, signalling and transcriptional regulation. Genome-wide average DNA methylation (GWAM) levels accounted for variability in some non-heritable lipids. We reveal a complex interplay of genetic and environmental influences on the ageing plasma lipidome.

Identification of Potential Diagnostic and Prognostic Pseudogenes in Hepatocellular Carcinoma Based on Pseudogene-miRNA-mRNA Competitive Network

Background

It is widely known that hepatocellular carcinoma (HCC) has high rates of morbidity and mortality. A large number of studies have indicated that pseudogenes have an important effect on the carcinogenesis of HCC. Pseudogenes can play a role through the ceRNA network. There have been numerous studies on lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA networks. However, the pseudogene-miRNA-mRNA network in HCC has rarely been researched or reported on.

Material/Methods

The Cancer Genome Atlas (TCGA) database was researched and differences between selected genes were studied. A pseudogene-miRNA-mRNA network was then constructed and clustering of pseudogenes was studied. The diagnostic value of the selected pseudogenes, their functions, and pathways were investigated using available databases to understand their possible pathogenic mechanism in HCC. The protein-protein interaction network of target genes was found and the top 10 hub genes were identified. Expression of hub genes in HCC tissues was then detected by RT-qPCR.

Results

By analyzing the gene difference and clinical data of HCC, we constructed a ceRNA network composed of 4 pseudogenes, 8 miRNAs, and 30 mRNAs. The pseudogenes AP000769.1, KRT16P1, KRT16P3, and RPLP0P2 were all correlated with the diagnosis and prognosis of HCC. Functional analyses through the Kyoto Encyclopedia of Genes and Genomes and the Gene Ontology databases indicated that pseudogenes can affect the physiological process of HCC through the p53 pathway. The top 10 hub genes identified were all highly expressed in HCC tissues and affected the patient survival rate.

Conclusions

In this study, 4 pseudogenes related to the diagnosis and prognosis of liver cancer were found through the construction of a ceRNA network. These 4 pseudogenes might constitute new therapeutic targets for liver cancer patients.

Re-recognition of pseudogenes: From molecular to clinical applications

Pseudogenes were initially regarded as "nonfunctional" genomic elements that did not have protein-coding abilities due to several endogenous inactivating mutations. Although pseudogenes are widely expressed in prokaryotes and eukaryotes, for decades, they have been largely ignored and classified as gene "junk" or "relics". With the widespread availability of high-throughput sequencing analysis, especially omics technologies, knowledge concerning pseudogenes has substantially increased. Pseudogenes are evolutionarily conserved and derive primarily from a mutation or retrotransposon, conferring the pseudogene with a "gene repository" role to store and expand genetic information. In contrast to previous notions, pseudogenes have a variety of functions at the DNA, RNA and protein levels for broadly participating in gene regulation to influence the development and progression of certain diseases, especially cancer. Indeed, some pseudogenes have been proven to encode proteins, strongly contradicting their "trash" identification, and have been confirmed to have tissue-specific and disease subtype-specific expression, indicating their own value in disease diagnosis. Moreover, pseudogenes have been correlated with the life expectancy of patients and exhibit great potential for future use in disease treatment, suggesting that they are promising biomarkers and therapeutic targets for clinical applications. In this review, we summarize the natural properties, functions, disease involvement and clinical value of pseudogenes. Although our knowledge of pseudogenes remains nascent, this field deserves more attention and deeper exploration.

The Protein-Coding Human Genome: Annotating High-Hanging Fruits

The major transcript variants of human protein-coding genes are annotated to a certain degree of accuracy combining manual curation, transcript data, and proteomics evidence. However, there is considerable disagreement on the annotation of about 2000 genes-they can be protein-coding, noncoding, or pseudogenes-and on the annotation of most of the predicted alternative transcripts. Pure transcriptome mapping approaches seem to be limited in discriminating functional expression from noise. These limitations have partially been overcome by dedicated algorithms to detect alternative spliced micro-exons and wobble splice variants. Recently, knowledge about splice mechanism and protein structure are incorporated into an algorithm to predict neighboring homologous exons, often spliced in a mutually exclusive manner. Predicted exons are evaluated by transcript data, structural compatibility, and evolutionary conservation, revealing hundreds of novel coding exons and splice mechanism re-assignments. The emerging human pan-genome is necessitating distinctive annotations incorporating differences between individuals and between populations.

The pseudogene PTENP1 regulates smooth muscle cells as a competing endogenous RNA

The long non-coding RNA (lncRNA) PTENP1 is a pseudogene of phosphatase and tensin homologue deleted on chromosome ten (PTEN), has been implicated in smooth muscle cell (SMC) proliferation and apoptosis. PTENP1 is the pseudogene of PTEN. However, it is unclear whether and how PTENP1 functions in the proliferation and apoptosis of human aortic SMCs (HASMCs). Here, we hypothesised that PTENP1 inhibits HASMC proliferation and enhances apoptosis by promoting PTEN expression. PCR analysis and Western blot assays respectively showed that both PTENP1 and PTEN were up-regulated in human aortic dissection (AD) samples. PTENP1 overexpression significantly increased the protein expression of PTEN, promoted apoptosis and inhibited the proliferation of HASMCs. PTENP1 silencing exhibited the opposite effects and mitigated H₂O₂-induced apoptosis of HASMCs. In an angiotensin II (Ang II)-induced mouse aortic aneurysm (AA) model, PTENP1 overexpression potentiated aortic SMC apoptosis, exacerbated aneurysm formation. Mechanistically, RNA pull-down assay and a series of luciferase reporter assays using miR-21 mimics or inhibitors identified PTENP1 as a molecular sponge for miR-21 to endogenously compete for the binding between miR-21 and the PTEN transcript, releasing PTEN expression. This finding was further supported by in vitro immunofluorescent evidence showing decreased cell apoptosis upon miR-21 mimic administration under baseline PTENP1 overexpression. Ex vivo rescue of PTEN significantly mitigated the SMC apoptosis induced by PTENP1 overexpression. Finally, Western blot assays showed substantially reduced Akt phosphorylation and cyclin D1 and cyclin E levels with up-regulated PTENP1 in HASMCs. Our study identified PTENP1 as a mediator of HASMC homeostasis and suggests that PTENP1 is a potential target in AD or AA intervention.

Identification of Novel Susceptibility Loci and Genes for Prostate Cancer Risk: A Transcriptome-Wide Association Study in Over 140,000 European Descendants

Genome-wide association study-identified prostate cancer risk variants explain only a relatively small fraction of its familial relative risk, and the genes responsible for many of these identified associations remain unknown. To discover novel prostate cancer genetic loci and possible causal genes at previously identified risk loci, we performed a transcriptome-wide association study in 79,194 cases and 61,112 controls of European ancestry. Using data from the Genotype-Tissue Expression Project, we established genetic models to predict gene expression across the transcriptome for both prostate models and cross-tissue models and evaluated model performance using two independent datasets. We identified significant associations for 137 genes at P < 2.61 × 10-6, a Bonferroni-corrected threshold, including nine genes that remained significant at P < 2.61 × 10-6 after adjusting for all known prostate cancer risk variants in nearby regions. Of the 128 remaining associated genes, 94 have not yet been reported as potential target genes at known loci. We silenced 14 genes and many showed a consistent effect on viability and colony-forming efficiency in three cell lines. Our study provides substantial new information to advance our understanding of prostate cancer genetics and biology. SIGNIFICANCE: This study identifies novel prostate cancer genetic loci and possible causal genes, advancing our understanding of the molecular mechanisms that drive prostate cancer.

Rapid evolution of α-gliadin gene family revealed by analyzing Gli-2 locus regions of wild emmer wheat

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

A Transcriptome-Wide Association Study Among 97,898 Women to Identify Candidate Susceptibility Genes for Epithelial Ovarian Cancer Risk

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10-6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis.Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419-30. ©2018 AACR.

Role of Pseudogenes in Tumorigenesis

Functional genomics has provided evidence that the human genome transcribes a large number of non-coding genes in addition to protein-coding genes, including microRNAs and long non-coding RNAs (lncRNAs). Among the group of lncRNAs are pseudogenes that have not been paid attention in the past, compared to other members of lncRNAs. However, increasing evidence points the important role of pseudogenes in diverse cellular functions, and dysregulation of pseudogenes are often associated with various human diseases including cancer. Like other types of lncRNAs, pseudogenes can also function as master regulators for gene expression and thus, they can play a critical role in various aspects of tumorigenesis. In this review we discuss the latest developments in pseudogene research, focusing on how pseudogenes impact tumorigenesis through different gene regulation mechanisms. Given the high sequence homology with the corresponding parent genes, we also discuss challenges for pseudogene research.

Formation of a Family of Long Intergenic Noncoding RNA Genes with an Embedded Translocation Breakpoint Motif in Human Chromosomal Low Copy Repeats of 22q11.2-Some Surprises and Questions

A family of long intergenic noncoding RNA (lincRNA) genes, FAM230 is formed via gene sequence duplication, specifically in human chromosomal low copy repeats (LCR) or segmental duplications. This is the first group of lincRNA genes known to be formed by segmental duplications and is consistent with current views of evolution and the creation of new genes via DNA low copy repeats. It appears to be an efficient way to form multiple lincRNA genes. But as these genes are in a critical chromosomal region with respect to the incidence of abnormal translocations and resulting genetic abnormalities, the 22q11.2 region, and also carry a translocation breakpoint motif, several intriguing questions arise concerning the presence and function of the translocation breakpoint sequence in RNA genes situated in LCR22s.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Transcriptome-Wide Association Study Identifies Susceptibility Loci and Genes for Age at Natural Menopause

OBJECTIVE

To identify novel susceptibility genes for age at natural menopause (ANM).

METHODS

Using transcription data generated in tissues from normal hypothalami (n = 73) and ovaries (n = 68) and high-density genotyping data provided by the Genotype-Tissue Expression (GTEx) database, we built 16 164 genetic models to predict gene expression across the transcriptome in these tissues. We used these models and summary statistics data from genome-wide association studies (GWAS) of ANM generated in 69 360 women of European ancestry to identify genes with their predicted expression related to ANM.

RESULTS

We found the predicted expression of 34 genes to be significantly associated with ANM at a Bonferroni-corrected threshold of P < 3.09 ×10^-6. These include 4 genes located more than 1 Mb away from any previously GWAS-identified ANM-associated variants, 24 genes that reside in known GWAS-identified loci but have not been previously implicated, and 6 genes previously implicated as ANM-associated genes.

CONCLUSION

Results from this transcriptome-wide association study, which integrated Expression quantitative trait loci (eQTL) data with summary statistics of GWAS of ANM, improves our understanding of the genetics and biology of female reproductive aging.

A family of long intergenic non-coding RNA genes in human chromosomal region 22q11.2 carry a DNA translocation breakpoint/AT-rich sequence

FAM230C, a long intergenic non-coding RNA (lincRNA) gene in human chromosome 13 (chr13) is a member of lincRNA genes termed family with sequence similarity 230. An analysis using bioinformatics search tools and alignment programs was undertaken to determine properties of FAM230C and its related genes. Results reveal that the DNA translocation element, the Translocation Breakpoint Type A (TBTA) sequence, which consists of satellite DNA, Alu elements, and AT-rich sequences is embedded in the FAM230C gene. Eight lincRNA genes related to FAM230C also carry the TBTA sequences. These genes were formed from a large segment of the 3’ half of the FAM230C sequence duplicated in chr22, and are specifically in regions of low copy repeats (LCR22)s, in or close to the 22q.11.2 region. 22q11.2 is a chromosomal segment that undergoes a high rate of DNA translocation and is prone to genetic deletions. FAM230C-related genes present in other chromosomes do not carry the TBTA motif and were formed from the 5’ half region of the FAM230C sequence. These findings identify a high specificity in lincRNA gene formation by gene sequence duplication in different chromosomes.

The Unforeseen Non-Coding RNAs in Head and Neck Cancer

Previously ignored non-coding RNAs (ncRNAs) have become the subject of many studies. However, there is an imbalance in the amount of consideration that ncRNAs are receiving. Some transcripts such as microRNAs (miRNAs) or small interfering RNAs (siRNAs) have gained much attention, but it is necessary to investigate other “pieces of the RNA puzzle”. These can offer a more complete view over normal and pathological cell behavior. The other ncRNA species are less studied, either due to their recent discovery, such as stable intronic sequence RNA (sisRNA), YRNA, miRNA-offset RNAs (moRNA), telomerase RNA component (TERC), natural antisense transcript (NAT), transcribed ultraconserved regions (T-UCR), and pseudogene transcript, or because they are still largely seen as non-coding transcripts with no relevance to pathogenesis. Moreover, some are still considered housekeeping RNAs, for instance small nucleolar RNAs (snoRNAs) and TERC. Our review summarizes the biogenesis, mechanism of action and potential role of less known ncRNAs in head and neck cancer, with a particular focus on the installment and progress for this particular cancer type.

HMGA1-pseudogenes and cancer

Pseudogenes are DNA sequences with high homology to the corresponding functional gene, but, because of the accumulation of various mutations, they have lost their initial functions to code for proteins. Consequently, pseudogenes have been considered until few years ago dysfunctional relatives of the corresponding ancestral genes, and then useless in the course of genome evolution. However, several studies have recently established that pseudogenes are owners of key biological functions. Indeed, some pseudogenes control the expression of functional genes by competitively binding to the miRNAs, some of them generate small interference RNAs to negatively modulate the expression of functional genes, and some of them even encode functional mutated proteins. Here, we concentrate our attention on the pseudogenes of the HMGA1 gene, that codes for the HMGA1a and HMGA1b proteins having a critical role in development and cancer progression. In this review, we analyze the family of HMGA1 pseudogenes through three aspects: classification, characterization, and their possible function and involvement in cancer.

Ensemble Modeling Approach Targeting Heterogeneous RNA-Seq data: Application to Melanoma Pseudogenes

We studied the transcriptome landscape of skin cutaneous melanoma (SKCM) using 103 primary tumor samples from TCGA, and measured the expression levels of both protein coding genes and non-coding RNAs (ncRNAs). In particular, we emphasized pseudogenes potentially relevant to this cancer. While cataloguing the profiles based on the known biotypes, all the employed RNA-Seq methods generated just a small consensus of significant biotypes. We thus designed an approach to reconcile the profiles from all methods following a simple strategy: we selected genes that were confirmed as differentially expressed by the ensemble predictions obtained in a regression model. The main advantages of this approach are: 1) Selection of a high-confidence gene set identifying relevant pathways; 2) Use of a regression model whose covariates embed all method-driven outcomes to predict an averaged profile; 3) Method-specific assessment of prediction power and significance. Furthermore, the approach can be generalized to any biological system for which noisy RNA-Seq profiles are computed. As our analyses concerned bio-annotations of both high-quality protein coding genes and ncRNAs, we considered the associations between pseudogenes and parental genes (targets). Among the candidate targets that were validated, we identified PINK1, which is studied in patients with Parkinson and cancer (especially melanoma).

A subset of conserved mammalian long non-coding RNAs are fossils of ancestral protein-coding genes

Background

Only a small portion of human long non-coding RNAs (lncRNAs) appear to be conserved outside of mammals, but the events underlying the birth of new lncRNAs in mammals remain largely unknown. One potential source is remnants of protein-coding genes that transitioned into lncRNAs.

Results

We systematically compare lncRNA and protein-coding loci across vertebrates, and estimate that up to 5% of conserved mammalian lncRNAs are derived from lost protein-coding genes. These lncRNAs have specific characteristics, such as broader expression domains, that set them apart from other lncRNAs. Fourteen lncRNAs have sequence similarity with the loci of the contemporary homologs of the lost protein-coding genes. We propose that selection acting on enhancer sequences is mostly responsible for retention of these regions. As an example of an RNA element from a protein-coding ancestor that was retained in the lncRNA, we describe in detail a short translated ORF in the JPX lncRNA that was derived from an upstream ORF in a protein-coding gene and retains some of its functionality.

Conclusions

We estimate that ~ 55 annotated conserved human lncRNAs are derived from parts of ancestral protein-coding genes, and loss of coding potential is thus a non-negligible source of new lncRNAs. Some lncRNAs inherited regulatory elements influencing transcription and translation from their protein-coding ancestors and those elements can influence the expression breadth and functionality of these lncRNAs.

Electronic supplementary material

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

Xp11.22 deletions encompassing CENPVL1, CENPVL2, MAGED1 and GSPT2 as a cause of syndromic X-linked intellectual disability

By searching a clinical database of over 60,000 individuals referred for array-based CNV analyses and online resources, we identified four males from three families with intellectual disability, developmental delay, hypotonia, joint hypermobility and relative macrocephaly who carried small, overlapping deletions of Xp11.22. The maximum region of overlap between their deletions spanned ~430 kb and included two pseudogenes, CENPVL1 and CENPVL2, whose functions are not known, and two protein coding genes-the G1 to S phase transition 2 gene (GSPT2) and the MAGE family member D1 gene (MAGED1). Deletions of this ~430 kb region have not been previously implicated in human disease. Duplications of GSPT2 have been documented in individuals with intellectual disability, but the phenotypic consequences of a loss of GSPT2 function have not been elucidated in humans or mouse models. Changes in MAGED1 have not been associated with intellectual disability in humans, but loss of MAGED1 function is associated with neurocognitive and neurobehavioral phenotypes in mice. In all cases, the Xp11.22 deletion was inherited from an unaffected mother. Studies performed on DNA from one of these mothers did not show evidence of skewed X-inactivation. These results suggest that deletions of an ~430 kb region on chromosome Xp11.22 that encompass CENPVL1, CENPVL2, GSPT2 and MAGED1 cause a distinct X-linked syndrome characterized by intellectual disability, developmental delay, hypotonia, joint hypermobility and relative macrocephaly. Loss of GSPT2 and/or MAGED1 function may contribute to the intellectual disability and developmental delay seen in males with these deletions.

Protein-Coding Genes' Retrocopies and Their Functions

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

Last rolls of the yoyo: Assessing the human canonical protein count

In 2004, when the protein estimate from the finished human genome was only 24,000, the surprise was compounded as reviewed estimates fell to 19,000 by 2014. However, variability in the total canonical protein counts (i.e. excluding alternative splice forms) of open reading frames (ORFs) in different annotation portals persists. This work assesses these differences and possible causes. A 16-year analysis of Ensembl and UniProtKB/Swiss-Prot shows convergence to a protein number of ~20,000. The former had shown some yo-yoing, but both have now plateaued. Nine major annotation portals, reviewed at the beginning of 2017, gave a spread of counts from 21,819 down to 18,891. The 4-way cross-reference concordance (within UniProt) between Ensembl, Swiss-Prot, Entrez Gene and the Human Gene Nomenclature Committee (HGNC) drops to 18,690, indicating methodological differences in protein definitions and experimental existence support between sources. The Swiss-Prot and neXtProt evidence criteria include mass spectrometry peptide verification and also cross-references for antibody detection from the Human Protein Atlas. Notwithstanding, hundreds of Swiss-Prot entries are classified as non-coding biotypes by HGNC. The only inference that protein numbers might still rise comes from numerous reports of small ORF (smORF) discovery. However, while there have been recent cases of protein verifications from previous miss-annotation of non-coding RNA, very few have passed the Swiss-Prot curation and genome annotation thresholds. The post-genomic era has seen both advances in data generation and improvements in the human reference assembly. Notwithstanding, current numbers, while persistently discordant, show that the earlier yo-yoing has largely ceased. Given the importance to biology and biomedicine of defining the canonical human proteome, the task will need more collaborative inter-source curation combined with broader and deeper experimental confirmation in vivo and in vitro of proteins predicted in silico. The eventual closure could be well be below ~19,000.

Identification and analysis of a processed cytochrome P450 pseudogene of the disease vector Aedes aegypti

OBJECTIVE

To clone cytochrome P450 from Aedes aegypti (Ae. aegypti) and determine the characteristics using bioinformatics tools.

METHODS

Cytochrome P450 of Ae. aegypti was amplified using polymerase chain reaction, cloned and sequenced. Evolutionary relationship of the sequence was inferred and bioinformatics tools were used to predict subcellular localisation, signal peptide, transmembrane helix, phosphorylation, O-glycosylation, secondary and tertiary structures of the deduced protein.

RESULTS

Polymerase chain reaction rather amplified a cytochrome P450 pseudogene which was named CYP4H44P (GenBank accession number KF779932). The pseudogene has 1537 nucleotides and an open reading frame of 335 amino acids containing cytochrome P450 motifs except the WxxxR motif. It is highly homologous to CYP4H28 and CYP4H28v2. Phylogenetic analysis and evolutionary divergence showed strong clustering with CYP4H28 alleles and least divergence from the alleles respectively. The deduced protein was predicted to be found in the cytoplasm and likely to be phosphorylated but devoid of signal peptide, transmembrane helix and O-glycosylated sites. The secondary and tertiary structures were also generated.

CONCLUSIONS

A cytochrome P450 pseudogene, CYP4H44P was cloned from Ae. aegypti. The pseudogene is homologous with CYP4H28 alleles and seems to have recently diverged from this group. Isolating this pseudogene is an important step for evaluating its biological role in the mosquito and for the evolutionary analysis of Ae. aegypti CYPs.

Identification of Multiple Forms of RNA Transcripts Associated with Human-Specific Retrotransposed Gene Copies

The human genome contains thousands of retrocopies, mostly as processed pseudogenes, which were recently shown to be prevalently transcribed. In particular, those specifically acquired in the human lineage are able to modulate gene expression in a manner that contributed to the evolution of human-specific traits. Therefore, knowledge of the human-specific retrocopies that are transcribed or their full-length transcript structure contributes to better understand human genome evolution. In this study, we identified 16 human-specific retrocopies that harbor 5' CpG islands by in silico analysis and showed that 12 were transcribed in normal tissues and cancer cell lines with a variety of expression patterns, including cancer-specific expression. Determination of the structure of the transcripts associated with the retrocopies revealed that none were transcribed from their 5' CpG islands, but rather, from inside the 3' UTR and the nearby 5' flanking region of the retrocopies as well as the promoter of neighboring genes. The multiple forms of the transcripts, such as chimeric and individual transcripts in both the sense and antisense orientation, might have introduced novel post-transcriptional regulation into the genome during human evolution. These results shed light on the potential role of human-specific retrocopies in the evolution of gene regulation and genomic disorders.

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.

RP1-13D10.2 Is a Novel Modulator of Statin-Induced Changes in Cholesterol

BACKGROUND

Numerous genetic contributors to cardiovascular disease risk have been identified through genome-wide association studies; however, identifying the molecular mechanism underlying these associations is not straightforward. The Justification for the Use of Statins in Primary Prevention: An Intervention Trial Evaluating Rosuvastatin (JUPITER) trial of rosuvastatin users identified a sub-genome-wide association of rs6924995, a single-nucleotide polymorphism ≈10 kb downstream of myosin regulatory light chain interacting protein (MYLIP, aka IDOL and inducible degrader of low-density lipoprotein receptor [LDLR]), with LDL cholesterol statin response. Interestingly, although this signal was initially attributed to MYLIP, rs6924995 lies within RP1-13D10.2, an uncharacterized long noncoding RNA.

METHODS AND RESULTS

Using simvastatin and sham incubated lymphoblastoid cell lines from participants of the Cholesterol and Pharmacogenetics (CAP) simvastatin clinical trial, we found that statin-induced change in RP1-13D10.2 levels differed between cell lines from the tails of the white and black low-density lipoprotein cholesterol response distributions, whereas no difference in MYLIP was observed. RP1-13D10.2 overexpression in Huh7 and HepG2 increased LDLR transcript levels, increased LDL uptake, and decreased media levels of apolipoprotein B. In addition, we found a trend of slight differences in the effects of RP1-13D10.2 overexpression on LDLR transcript levels between hepatoma cells transfected with the rs6924995 A versus G allele and a suggestion of an association between rs6924995 and RP1-10D13.2 expression levels in the CAP lymphoblastoid cell lines. Finally, RP1-13D10.2 expression levels seem to be sterol regulated, consistent with its potential role as a novel lipid regulator.

CONCLUSIONS

RP1-13D10.2 is a long noncoding RNA that regulates LDLR and may contribute to low-density lipoprotein cholesterol response to statin treatment. These findings highlight the potential role of noncoding RNAs as determinants of interindividual variation in drug response.

Functional Impact and Evolution of a Novel Human Polymorphic Inversion That Disrupts a Gene and Creates a Fusion Transcript

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000–50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.

Since the discovery of chromosomal inversions almost 100 years ago, how they are maintained in natural populations has been a highly debated issue. One of the hypotheses is that inversion breakpoints could affect genes and modify gene expression levels, although evidence of this came only from laboratory mutants. In humans, a few inversions have been shown to associate with expression differences, but in all cases the molecular causes have remained elusive. Here, we have carried out a complete characterization of a new human polymorphic inversion and determined that it is specific to East Asian populations. In addition, we demonstrate that it disrupts the ZNF257 gene and, through the translocation of the first exon and regulatory sequences, creates a previously nonexistent fusion transcript, which together are associated to expression changes in several other genes. Finally, we investigate the potential evolutionary and phenotypic consequences of the inversion, and suggest that it is probably deleterious. This is therefore the first example of a natural polymorphic inversion that has position effects and creates a new chimeric gene, contributing to answer an old question in evolutionary biology.

piRNAs derived from ancient viral processed pseudogenes as transgenerational sequence-specific immune memory in mammals

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences within vertebrate genomes derived from reverse transcription and integration of ancient bornaviral nucleoprotein mRNA via the host retrotransposon machinery. While species with EBLNs appear relatively resistant to bornaviral disease, the nature of this association is unclear. We hypothesized that EBLNs could give rise to antiviral interfering RNA in the form of PIWI-interacting RNAs (piRNAs), a class of small RNA known to silence transposons but not exogenous viruses. We found that in both rodents and primates, which acquired their EBLNs independently some 25-40 million years ago, EBLNs are present within piRNA-generating regions of the genome far more often than expected by chance alone (ℙ = 8 × 10(-3)-6 × 10(-8)). Three of the seven human EBLNs fall within annotated piRNA clusters and two marmoset EBLNs give rise to bona fide piRNAs. In both rats and mice, at least two of the five EBLNs give rise to abundant piRNAs in the male gonad. While no EBLNs are syntenic between rodent and primate, some of the piRNA clusters containing EBLNs are; thus we deduce that EBLNs were integrated into existing piRNA clusters. All true piRNAs derived from EBLNs are antisense relative to the proposed ancient bornaviral nucleoprotein mRNA. These observations are consistent with a role for EBLN-derived piRNA-like RNAs in interfering with ancient bornaviral infection. They raise the hypothesis that retrotransposon-dependent virus-to-host gene flow could engender RNA-mediated, sequence-specific antiviral immune memory in metazoans analogous to the CRISPR/Cas system in prokaryotes.

Pseudogenes in Human Cancer

Recent advances in the analysis of RNA sequencing data have shown that pseudogenes are highly specific markers of cell identity and can be used as diagnostic and prognostic markers. Furthermore, genetically engineered mouse models have recently provided compelling support for a causal link between altered pseudogene expression and cancer. In this review, we discuss the most recent milestones reached in the pseudogene field and the use of pseudogenes as cancer classifiers.

DNA methylation of the LIN28 pseudogene family

Background

DNA methylation directs the epigenetic silencing of selected regions of DNA, including the regulation of pseudogenes, and is widespread throughout the genome. Pseudogenes are decayed copies of duplicated genes that have spread throughout the genome by transposition. Pseudogenes are transcriptionally silenced by DNA methylation, but little is known about how pseudogenes are targeted for methylation or how methylation levels are maintained in different tissues.

Results

We employed bisulfite next generation sequencing to examine the methylation status of the LIN28 gene and four processed pseudogenes derived from LIN28. The objective was to determine whether LIN28 pseudogenes maintain the same pattern of methylation as the parental gene or acquire a methylation pattern independent of the gene of origin. In this study, we determined that the methylation status of LIN28 pseudogenes does not resemble the pattern evident for the LIN28 gene, but rather these pseudogenes appear to acquire methylation patterns independent of the parental gene. Furthermore, we observed that methylation levels of the examined pseudogenes correlate to the location of insertion within the genome. LIN28 pseudogenes inserted into gene bodies were highly methylated in all tissues examined. In contrast, pseudogenes inserted into genomic regions that are not proximal to genes were differentially methylated in various tissue types.

Conclusions

Our analysis suggests that Lin28 pseudogenes do not aquire patterns of tissue-specific methylation as for the parental gene, but rather are methylated in patterns specific to the local genomic environment into which they were inserted.

Electronic supplementary material

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

Anaplastic Thyroid Carcinoma: A ceRNA Analysis Pointed to a Crosstalk between SOX2, TP53, and microRNA Biogenesis

It has been suggested that cancer stem cells (CSC) may play a central role in oncogenesis, especially in undifferentiated tumours. Anaplastic thyroid carcinoma (ATC) has characteristics suggestive of a tumour enriched in CSC. Previous studies suggested that the stem cell factor SOX2 has a preeminent hierarchical role in determining the characteristics of stem cells in SW1736 ATC cell line. In detail, silencing SOX2 in SW1736 is able to suppress the expression of the stem markers analysed, strongly sensitizing the line to treatment with chemotherapeutic agents. Therefore, in order to further investigate the role of SOX2 in ATC, a competing endogenous RNA (ceRNA) analysis was conducted in order to isolate new functional partners of SOX2. Among the interactors, of particular interest are genes involved in the biogenesis of miRNAs (DICER1, RNASEN, and EIF2C2), in the control cell cycle (TP53, CCND1), and in mitochondrial activity (COX8A). The data suggest that stemness, microRNA biogenesis and functions, p53 regulatory network, cyclin D1, and cell cycle control, together with mitochondrial activity, might be coregulated.

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.