Pseudogenes: pseudo or real functional elements?

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From > 85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From >85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

Proterozoic Acquisition of Archaeal Genes for Extracellular Electron Transfer: A Metabolic Adaptation of Aerobic Ammonia-Oxidizing Bacteria to Oxygen Limitation

Many aerobic microbes can utilize alternative electron acceptors under oxygen-limited conditions. In some cases, this is mediated by extracellular electron transfer (or EET), wherein electrons are transferred to extracellular oxidants such as iron oxide and manganese oxide minerals. Here, we show that an ammonia-oxidizer previously known to be strictly aerobic, Nitrosomonas communis, may have been able to utilize a poised electrode to maintain metabolic activity in anoxic conditions. The presence and activity of multiheme cytochromes in N. communis further suggest a capacity for EET. Molecular clock analysis shows that the ancestors of β-proteobacterial ammonia oxidizers appeared after Earth's atmospheric oxygenation when the oxygen levels were >10-4pO2 (present atmospheric level [PAL]), consistent with aerobic origins. Equally important, phylogenetic reconciliations of gene and species trees show that the multiheme c-type EET proteins in Nitrosomonas and Nitrosospira lineages were likely acquired by gene transfer from γ-proteobacteria when the oxygen levels were between 0.1 and 1 pO2 (PAL). These results suggest that β-proteobacterial EET evolved during the Proterozoic when oxygen limitation was widespread, but oxidized minerals were abundant.

Cymric, a Maternal and Zygotic HTK-16-Like SHARK Family Tyrosine Kinase Gene, Is Disrupted in Molgula occulta, a Tailless Ascidian

AbstractWe describe the cloning and expression of a nonreceptor tyrosine kinase, cymric (Uro-1), a HTK-16-like (HydraTyrosineKinase-16) gene, identified in a subtractive screen for maternal ascidian cDNAs in Molgula oculata, an ascidian species with a tadpole larva. The cymric gene encodes a 4-kb mRNA expressed in gonads, eggs, and embryos in the tailed M. oculata but is not detected in eggs or embryos of the closely related tailless species Molgula occulta. There is a large insertion in cymric in the M. occulta genome, as shown by transcriptome and genome analyses, resulting in it becoming a pseudogene. The cymric amino acid sequence encodes a nonreceptor tyrosine kinase with an N-terminal region containing two SH2 domains and five ankyrin repeats, similar to the HTK-16-like gene found in other ascidians. Thus, the ascidian cymric genes are members of the SHARK (Src-homology ankyrin-repeat containing tyrosine kinase) family of nonreceptor tyrosine kinases, which are found throughout invertebrates and missing from vertebrates. We show that cymric is lacking the tyrosine kinase domain in the tailless M. occulta, although the truncated mRNA is still expressed in transcriptome data. This maternal and zygotic HTK-16-like tyrosine kinase is another described pseudogene from M. occulta and appears not to be necessary for adult development.

Anti-apoptotic genes and non-coding RNAs are potential outcome predictors for ulcerative colitis

Due to the lack of clinical, immunologic, genetic, and laboratory markers to predict remission in ulcerative colitis (UC) without relapse, there is no clear recommendation regarding withdrawal of therapy. Therefore, this study was to investigate if transcriptional analysis together with Cox survival analysis might be able to reveal molecular markers that are specific for remission duration and outcome. Mucosal biopsies from patients in remission with active treatment-naïve UC and healthy control subjects underwent whole-transcriptome RNA-seq. Principal component analysis (PCA) and Cox proportional hazards regression analysis were applied to the remission data concerning duration and status of patients. A randomly chosen remission sample set was used for validation of the applied methods and results. The analyses distinguished two different UC remission patient groups with respect to remission duration and outcome (relapse). Both groups showed that altered states of UC with quiescent microscopic disease activity are still present. The patient group with the longest remission duration and no relapse revealed specific and increased expression of antiapoptotic factors belonging to the MTRNR2-like gene family and non-coding RNAs. In summary, the expression of anti-apoptotic factors and non-coding RNAs may contribute to personalized medicine approaches in UC by improving patient stratification for different treatment regimens.

The diagnostic and prognostic value of pseudogene SIGLEC17P in lung adenocarcinoma and a preliminary functional study

Among malignant tumors, lung adenocarcinoma (LUAD) is the leading cause of death worldwide. This study explored the diagnostic, prognostic value, and preliminary functional verification of sialic acid binding Ig like lectin 17, pseudogene (SIGLEC17P) in LUAD. Prognostic lncRNAs for LUAD were identified by The Cancer Genome Atlas and quantitative real-time PCR (qRT-PCR) was used to detect the expression of SIGLEC17P in LUAD and paracarcinoma tissues. Subsequently, lentiviral vectors were used to overexpress SIGLEC17P in A549 and H1299 cells. The effects of SIGLEC17P overexpression on the proliferation, migration, and invasiveness of LUAD cells (A549 and H1299) were evaluated by Cell Counting Kit-8, wound healing, and transwell migration assays, respectively. Bioinformatics analyses were performed to reveal the potential pathways in which SIGLEC17P is involved in LUAD. qRT-PCR results revealed low SIGLEC17P expression in LUAD tissues and a significant association with the N stage, T stage, and tumor node metastasis stage. Furthermore, the receiver operating characteristic curve demonstrated a reliable diagnostic value. The proliferation, migration, and invasion of LUAD cells were inhibited by overexpression of SIGLEC17P. Bioinformatics analyses suggested that SIGLEC17P might exert antioncogenic effects in LUAD through the mir-20-3p/ADH1B or mir-4476-5p/DPYSL axis. In summary, our results revealed that SIGLEC17P acts as a prognostic biomarker, independent prognostic factor, and potential therapeutic target for patients with LUAD.

Immune system disruptions implicated in whole blood epigenome-wide association study of depression among Parkinson's disease patients

Although Parkinson's Disease (PD) is typically described in terms of motor symptoms, depression is a common feature. We explored whether depression influences blood-based genome-wide DNA methylation (DNAm) in 692 subjects from a population-based PD case-control study, using both a history of clinically diagnosed depression and current depressive symptoms measured by the geriatric depression scale (GDS). While PD patients in general had more immune activation and more accelerated epigenetic immune system aging than controls, the patients experiencing current depressive symptoms (GDS≥5) showed even higher levels of both markers than patients without current depressive symptoms (GDS<5). For PD patients with a history of clinical depression compared to those without, we found no differences in immune cell composition. However, a history of clinical depression among patients was associated with differentially methylated CpGs. Epigenome-wide association analysis (EWAS) revealed 35 CpGs associated at an FDR≤0.05 (569 CpGs at FDR≤0.10, 1718 CpGs at FDR≤0.15). Gene set enrichment analysis implicated immune system pathways, including immunoregulatory interactions between lymphoid and non-lymphoid cells (p-adj = 0.003) and cytokine-cytokine receptor interaction (p-adj = 0.004). Based on functional genomics, 25 (71%) of the FDR≤0.05 CpGs were associated with genetic variation at 45 different methylation quantitative trait loci (meQTL). Twenty-six of the meQTLs were also expression QTLs (eQTLs) associated with the abundance of 53 transcripts in blood and 22 transcripts in brain (substantia nigra, putamen basal ganglia, or frontal cortex). Notably, cg15199181 was strongly related to rs823114 (SNP-CpG p-value = 3.27E-310), a SNP identified in a PD meta-GWAS and related to differential expression of PM20D1, RAB29, SLC41A1, and NUCKS1. The entire set of genes detected through functional genomics was most strongly overrepresented for interferon-gamma-mediated signaling pathway (enrichment ratio = 18.8, FDR = 4.4e-03) and T cell receptor signaling pathway (enrichment ratio = 13.2, FDR = 4.4e-03). Overall, the current study provides evidence of immune system involvement in depression among Parkinson's patients.

•

Parkinson's disease (PD) is associated with clinical depression prior to PD onset and depressive symptoms after PD diagnosis.

•

Epigenome-wide analysis revealed CpGs related to current depressive symptoms and a history of clinical depression among PD patients.

•

Patients experiencing current depressive symptoms had the highest epigenetic-based neutrophil-to-lymphocyte ratio on average.

•

Patients with a history of clinical depression had differentially methylated CpGs in genes enriched for immune system pathways.

•

Many of the depression associated CpGs were linked to differential expression through meQTL/eQTLs, which included GWAS variants.

Comparative Cytogenetic and Sequence Analysis of U Small Nuclear RNA Genes in Three Ancistrus Species (Siluriformes: Loricariidae)

Ancistrus presents a wide karyotypic diversity, resulting from numeric and structural chromosomal rearrangements. It has been proposed that some genome-specific regions containing repetitive units could organize prone-to-break DNA sites in Loricariidae, triggering chromosomal rearrangements such as Robertsonian fusions (Rb fusions), centric fissions, translocations, and inversions. The tandemly repeats of the small nuclear RNAs (snRNAs) gene families are considered good cytogenetic markers for understanding chromosomal remodeling events among closely related species, but these snRNAs have been scarcely analyzed in Ancistrus. This study presented the nucleotide sequencing and comparative in situ location of U snRNA sequences from Ancistrus aguaboensis, Ancistrus cf. multispinis, and Ancistrus sp. (2n = 50, 52, and 50, respectively), aiming to provide information about snRNA clusters in the genome and chromosome evolution in Ancistrus. U snRNA nucleotide sequences of Ancistrus presented identity to orthologous copies and folded their secondary structures correctly. In situ localization and karyotyping of the three Ancistrus species revealed clustered copies of U2 and U5 snRNA gene families to a single chromosome site, one chromosome pair bearing U1 snRNA sequence, and one main locus of U4 snRNA sequence, besides scattered signals along the chromosomes. Previous studies related the participation of the rRNA gene families in centric fusion events, contributing to chromosome rearrangements and karyotype plasticity present in Loricariidae. In this study, homeologies in U snRNA loci chromosomal locations were detected, indicating the occurrence of conserved sites of these gene families in these three Ancistrus species with 2n = 50 or 52 chromosomes.

Major and minor U small nuclear RNAs genes characterization in a neotropical fish genome: Chromosomal remodeling and repeat units dispersion in Parodontidae

In association with many proteins, small nuclear RNAs (snRNAs) organize the spliceosomes that play a significant role in processing precursor mRNAs during gene expression. According to snRNAs genic arrangements, two kinds of spliceosomes (major and minor) can be organized into eukaryotic cells. Although in situ localization of U1 and U2 snDNAs have been performed in fish karyotypes, studies with genomic characterization and functionality of U snRNAs integrated into chromosomal changes on Teleostei are still scarce. This study aimed to achieve a genomic characterization of the U snRNAs genes in Apareiodon sp. (2n = 54, ZZ/ZW), apply these data to recognize functional/defective copies, and map chromosomal changes involving snDNAs in Parodontidae species karyotype diversification. Nine snRNA multigene families (U1, U2, U4, U5, U6, U11, U12, U4atac and U6atac) arranged in putatively functional copies in the genome were analyzed. Proximal Sequence Elements (PSE) and TATA-box promoters occurrence, besides an entire transcribed region and conserved secondary structures, qualify them for spliceosome activity. In addition, several defective copies or pseudogenes were identified for the snRNAs that make up the major spliceosome. In situ localization of snDNAs in five species of Parodontidae demonstrated that U1, U2, and U4 snDNAs were involved in chromosomal location changes or units dispersion. The U snRNAs defective/pseudogenes units dispersion could be favored by the probable occurrence of active retrotransposition enzymes in the Apareiodon genome. The U2 and U4 snDNAs sites were involved in independent events in the differentiation of sex chromosomes among Parodontidae lineages. The study characterized U snRNA genes that compose major and minor spliceosomes in the Apareiodon sp. genome and proposes that their defective copies trigger chromosome differentiation and diversification events in Parodontidae.

The World of Pseudogenes: New Diagnostic and Therapeutic Targets in Cancers or Still Mystery Molecules?

Pseudogenes were once considered as “junk DNA”, due to loss of their functions as a result of the accumulation of mutations, such as frameshift and presence of premature stop-codons and relocation of genes to inactive heterochromatin regions of the genome. Pseudogenes are divided into two large groups, processed and unprocessed, according to their primary structure and origin. Only 10% of all pseudogenes are transcribed into RNAs and participate in the regulation of parental gene expression at both transcriptional and translational levels through senseRNA (sRNA) and antisense RNA (asRNA). In this review, about 150 pseudogenes in the different types of cancers were analyzed. Part of these pseudogenes seem to be useful in molecular diagnostics and can be detected in various types of biological material including tissue as well as biological fluids (liquid biopsy) using different detection methods. The number of pseudogenes, as well as their function in the human genome, is still unknown. However, thanks to the development of various technologies and bioinformatic tools, it was revealed so far that pseudogenes are involved in the development and progression of certain diseases, especially in cancer.

Studying the Oncosuppressive Functions of PTENP1 as a ceRNA

PTENP1 is a processed pseudogene of the tumour suppressor phosphatase and tensin homolog deleted on chromosome 10 (PTEN). It functions posttranscriptionally to regulate PTEN by acting as a sponge for microRNAs that target PTEN. PTENP1 therefore functions as a competitive endogenous RNA (ceRNA), competing with PTEN for binding of microRNAs (miRNA) and thereby modulating PTEN cellular abundance. Studies of the overexpression of PTENP1 all confirm its oncosuppressive function to be mediated through the suppression of cell proliferation, induction of apoptosis, and inhibition of cell migration and invasion of cancer cells of differing types. These oncosuppressive functions are a direct consequence of miRNA binding by PTENP1 and the subsequent liberation of PTEN from miRNA induced suppression. In this chapter, we will focus initially on the description of a high efficiency transient transfection method to introduce and overexpress PTENP1 in the cell type of interest, followed by accurate methodologies to measure transfection efficiency by flow cytometry. We will then continue to describe two methods to analyze cell proliferation, namely the CCK-8 assay and Click-iT^® EdU assay. Due to commonalities in the manifestation of the oncosuppressive effects of PTENP1, mediated through its role as a ceRNA, the methods presented in this chapter will have wide applicability to a variety of different cell types.

Gauging the trends of pseudogenes in plants

Pseudogenes, the debilitated parts of ancient genes, were previously scrapped off as junk or discarded genes with no functional significance. Pseudogenes have come under scrutiny for their functionality, since recent studies have unveiled their importance in the regulation of their corresponding parent genes and various biological mechanisms. Despite the enormous occurrence of pseudogenes in plants, the lack of experimental validation has contributed toward their unresolved roles in gene regulation. Contrarily, most of the studies associated with gene regulation have been mainly reported for humans, mice, and other mammalian genomes. Consequently, in order to present a cumulative report on plant-based pseudogenes research, an attempt has been made to assemble multiple studies presenting the pseudogene classification, the prediction and the determination of comparative accuracies of various computational pipelines, and recent trends in analyzing their biological functions, and regulatory mechanisms. This review represents the classical, as well as the recent advances on pseudogene identification and their potential roles in transcriptional regulation, which could possibly invigorate the quality of genome annotation, evolutionary analysis, and complexity surrounding the regulatory pathways in plants. Thus, when the ambiguous boundary girdling the pseudogenes eventually recedes on account of their explicit orchestration role, research in flora would no longer saunter compared to that on fauna.

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.

Competing endogenous RNAs in lung cancer

Competing endogenous RNAs (ceRNAs) containing microRNA response elements can competitively interact with microRNA via miRNA response elements, which can combine non-coding RNAs with protein-coding RNAs through complex ceRNA networks. CeRNAs include non-coding RNAs (long non-coding RNAs, circular RNAs, and transcribed pseudogenes) and protein-coding RNAs (mRNAs). Molecular interactions in ceRNA networks can coordinate many biological processes; however, they may also lead to ceRNA network imbalance and thus contribute to cancer occurrence when disturbed. Recent studies indicate that many dysregulated RNAs derived from lung cancer may function as ceRNAs to regulate multitudinous biological functions for lung cancer, including tumor cell proliferation, apoptosis, growth, invasion, migration, and metastasis. This study therefore reviewed the research progress in the field of non-coding and protein-coding RNAs as ceRNAs in lung cancer, and highlighted validated ceRNAs involved in biological lung cancer functions. Furthermore, the roles of ceRNAs as novel prognostic and diagnostic biomarkers were also discussed. Interpreting the involvement of ceRNAs networks in lung cancer will provide new insight into cancer pathogenesis and treatment strategies.

lncRNA DUXAP8 Facilitates Multiple Malignant Phenotypes and Resistance to PARP Inhibitor in HCC via Upregulating FOXM1

In this study, we examined the clinical significance and molecular mechanisms of a long non-coding RNA (lncRNA), double homeobox A pseudogene 8 (DUXAP8) in hepatocellular carcinoma (HCC). DUXAP8 expression was compared using quantitative real-time PCR in HCC versus adjacent tissues and in HCC cell lines versus normal hepatic epithelial cells. The correlations between DUXAP8 level and clinicopathological features were analyzed. Assays including MTT, colony-forming analysis, Transwell assay, western blot, xenograft formation, experimental metastasis, luciferase assay, RNA pull-down, and RNA immunoprecipitation were used to examine DUXAP8-induced malignant phenotypes, its regulation on forkhead box protein M1 (FOXM1), and the importance of FOXM1 in mediating DUXAP8 phenotypes. Our results showed that DUXAP8 was significantly upregulated in HCC tissues or cell lines associated with tumors of advanced grades, tumors that were positive for lymph node metastasis, and patients with poor overall survival. DUAXP8 was essential in maintaining multiple malignant phenotypes (including resistance to olaparib) both in vitro and in vivo. Mechanistically, DUXAP8 upregulated FOXM1 expression by sponging miR-485-5p and interacting with the RNA-binding protein Fused in Sarcoma (FUS). Functionally, FOXM1 essentially mediated the oncogenic phenotypes of DUXAP8. Collectively, DUXAP8 acts through two distinct mechanisms to upregulate FOXM1 and becomes a pleotropic oncogenic lncRNA in HCC.

Review on the Computational Genome Annotation of Sequences Obtained by Next-Generation Sequencing

Next-Generation Sequencing (NGS) has made it easier to obtain genome-wide sequence data and it has shifted the research focus into genome annotation. The challenging tasks involved in annotation rely on the currently available tools and techniques to decode the information contained in nucleotide sequences. This information will improve our understanding of general aspects of life and evolution and improve our ability to diagnose genetic disorders. Here, we present a summary of both structural and functional annotations, as well as the associated comparative annotation tools and pipelines. We highlight visualization tools that immensely aid the annotation process and the contributions of the scientific community to the annotation. Further, we discuss quality-control practices and the need for re-annotation, and highlight the future of annotation.

Influence of Genetic Variance on Biomarker Levels After Occupational Exposure to 1,6-Hexamethylene Diisocyanate Monomer and 1,6-Hexamethylene Diisocyanate Isocyanurate

We evaluated the impact of genetic variance on biomarker levels in a population of workers in the automotive repair and refinishing industry who were exposed to respiratory sensitizers 1,6-hexamethylene diisocyanate (HDI) monomer and one of its trimers, HDI isocyanurate. The exposures and respective urine and plasma biomarkers 1,6-diaminohexane (HDA) and trisaminohexyl isocyanurate (TAHI) were measured in 33 workers; and genome-wide microarrays (Affymetrix 6.0) were used to genotype the workers' single-nucleotide polymorphisms (SNPs). Linear mixed model analyses have indicated that interindividual variations in both inhalation and skin exposures influenced these biomarker levels. Using exposure values as covariates and a false discovery rate < 0.10 to assess statistical significance, we observed that seven SNPs were associated with HDA in plasma, five were associated with HDA in urine, none reached significance for TAHI in plasma, and eight were associated with TAHI levels in urine. The different genotypes for the 20 significant SNPs accounted for 4- to 16-fold changes observed in biomarker levels. Associated gene functions include transcription regulation, calcium ion transport, vascular morphogenesis, and transforming growth factor beta signaling pathway, which may impact toxicokinetics indirectly by altering inflammation levels. Additionally, in an expanded analysis using a minor allele cutoff of 0.05 instead of 0.10, there were biomarker-associated SNPs within three genes that have been associated with isocyanate-induced asthma: ALK, DOCK2, and LHPP. We demonstrate that genetic variance impacts the biomarker levels in workers exposed to HDI monomer and HDI isocyanurate and that genetics can be used to refine exposure predictions in small cohorts when quantitative personal exposure and biomarker measurements are included in the models.

The genomic context of retrocopies increases their chance of functional relevancy in mammals

Described as "junk" DNA, pseudogenes are dead structures of previously active genes present in genomes. Pseudogenes are categorized into two main classes: processed pseudogenes, formed through retrotransposition, and non-processed pseudogenes, typically originated from gene decay following duplication events. The term "processed pseudogene" has changed to "retrocopy" since they are likely to evolve new functional roles and became a retrogene. Here, we surveyed 38,080 retrocopies from chimpanzee, dog, human, mouse, and rat genomes to assess their potential adaptive value. The retrocopies inserted in the same chromosome of the parental gene have higher chances of remain potentially "active" (absence of premature stop codons and frameshifts) (~26.1%), while those placed into a different chromosome have a twofold decrease chance of continuing potentially "active" (~7.52%). The genomic context of their placement seems associated with their expression. Retrocopies placed in intragenic regions and the same sense of the "host" gene have higher chances of being expressed relative to other genomic contexts. The proximity of retrocopies to their parental gene is associated with a lower decay rate, and their location likely influence their expression. Thus, despite their unclear role, retrocopies are probably involved in adaptive processes. Our results evidence natural selection acting in retrocopies.

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.

A functional role for the cancer disparity-linked genes, CRYβB2 and CRYβB2P1, in the promotion of breast cancer

Background

In the USA, the breast cancer mortality rate is 41% higher for African-American women than non-Hispanic White women. While numerous gene expression studies have classified biological features that vary by race and may contribute to poorer outcomes, few studies have experimentally tested these associations. CRYβB2 gene expression has drawn particular interest because of its association with overall survival and African-American ethnicity in multiple cancers. Several reports indicate that overexpression of the CRYβB2 pseudogene, CRYβB2P1, and not CRYβB2 is linked with race and poor outcome. It remains unclear whether either or both genes are linked to breast cancer outcomes. This study investigates CRYβB2 and CRYβB2P1 expression in human breast cancers and breast cancer cell line models, with the goal of elucidating the mechanistic contribution of CRYβB2 and CRYβB2P1 to racial disparities.

Methods

Custom scripts for CRYβB2 or CRYβB2P1 were generated and used to identify reads that uniquely aligned to either gene. Gene expression according to race and tumor subtype were assessed using all available TCGA breast cancer RNA sequencing alignment samples (n = 1221). In addition, triple-negative breast cancer models engineered to have each gene overexpressed or knocked out were developed and evaluated by in vitro, biochemical, and in vivo assays to identify biological functions.

Results

We provide evidence that CRYβB2P1 is expressed at higher levels in breast tumors compared to CRYβB2, but only CRYβB2P1 is significantly increased in African-American tumors relative to White American tumors. We show that independent of CRYβB2, CRYβB2P1 enhances tumorigenesis in vivo via promoting cell proliferation. Our data also reveal that CRYβB2P1 may function as a non-coding RNA to regulate CRYβB2 expression. A key observation is that the combined overexpression of both genes was found to suppress cell growth. CRYβB2 overexpression in triple-negative breast cancers increases invasive cellular behaviors, tumor growth, IL6 production, immune cell chemoattraction, and the expression of metastasis-associated genes. These data underscore that both CRYβB2 and CRYβB2P1 promote tumor growth, but their mechanisms for tumor promotion are likely distinct.

Conclusions

Our findings provide novel data emphasizing the need to distinguish and study the biological effects of both CRYβB2 and CRYβB2P1 as both genes independently promote tumor progression. Our data demonstrate novel molecular mechanisms of two understudied, disparity-linked molecules.

Electronic supplementary material

The online version of this article (10.1186/s13058-019-1191-3) contains supplementary material, which is available to authorized users.

Hypothesis: Compulsive hoarding as an atavism

In this hypothesis paper, we propose that hoarding disorder may be an ancestral mammalian behavior, possibly hard-wired but dormant in present-day humans, that can be activated by traumatic experiences. To support this claim, we describe hoarding in ancestral mammals and provide examples of recurrence of traits, or atavisms, in humans and other mammals, highlighting its implications for neuropsychiatric diseases. We highlight the potential of using shrews as animal models for hoarding behavior and, considering the current knowledge on the molecular underpinnings of the disorder, we also point to limitations of the proposed connection.

Pseudogene RACGAP1P activates RACGAP1/Rho/ERK signalling axis as a competing endogenous RNA to promote hepatocellular carcinoma early recurrence

Accumulating evidence has indicated crucial roles for pseudogenes in human cancers. However, the roles played by pseudogenes in the pathogenesis of HCC, particularly HCC early recurrence, still incompletely elucidated. Herein, we identify a novel early recurrence related pseudogene RACGAP1P which was significantly upregulated in HCC and was associated with larger tumour size, advanced clinical stage, abnormal AFP level and shorter survival time. In vitro and in vivo experiments have shown that RACGAP1P is a prerequisite for the development of malignant characteristics of HCC cells, including cell growth and migration. Mechanistic investigations indicated that RACGAP1P elicits its oncogenic activity as a ceRNA to sequestrate miR-15-5p from its endogenous target RACGAP1, thereby leading to the upregulation of RACGAP1 and the activation of RhoA/ERK signalling. These results may provide new insights into the functional crosstalk of the pseudogene/miRNA/parent-gene genetic network during HCC early relapse and may contribute to improving the clinical intervention for this subset of HCC patients.

Long Non-Coding RNAs and the Innate Immune Response

Innate immunity provides the initial defence against infection and it is now clear that long non-coding RNAs (lncRNAs) are important regulators of this response. Following activation of the innate response, we commonly see rapid induction of these lncRNAs and this is often mediated via the pro-inflammatory transcription factor, nuclear factor-κB (NF-κB). Knockdown studies have shown that lncRNAs tend to act in trans to regulate the expression of multiple inflammatory mediators and other responses. Mechanistically, many lncRNAs have demonstrated acting through heterogeneous nuclear ribonucleoproteins, complexes that are implicated chromatin re-modelling, transcription process and translation. In addition, these lncRNAs have also been shown to interact with multiple other proteins involved in the regulation of chromatin re-modelling, as well as those proteins involved in intracellular immune signalling, which include NF-κB. In this review, we will describe the evidence that supports this emerging role of lncRNA in the innate immune response.

Pseudogenes as Functionally Significant Elements of the Genome

Pseudogene is a gene copy that has lost its original function. For a long time, pseudogenes have been considered as "junk DNA" that inevitably arises as a result of ongoing evolutionary process. However, experimental data obtained during recent years indicate this understanding of the nature of pseudogenes is not entirely correct, and many pseudogenes perform important genetic functions. In the review, we have addressed classification of pseudogenes, methods of their detection in the genome, and the problem of their evolutionary conservatism and prevalence among species belonging to different taxonomic groups in the light of modern data. The mechanisms of gene expression regulation by pseudogenes and the role of pseudogenes in pathogenesis of various human diseases are discussed.

Marfan Syndrome Variability: Investigation of the Roles of Sarcolipin and Calcium as Potential Transregulator of FBN1 Expression

Marfan syndrome (MFS) is an autosomal dominant connective tissue disorder that displays a great clinical variability. Previous work in our laboratory showed that fibrillin-1 (FBN1) messenger RNA (mRNA) expression is a surrogate endpoint for MFS severity. Therefore, an expression quantitative trait loci (eQTL) analysis was performed to identify trans-acting regulators of FBN1 expression, and a significant signal reached genome-wide significant threshold on chromosome 11. This signal delineated a region comprising one expressed gene, SLN (encoding sarcolipin), and a single pseudogene, SNX7-ps1 (CTD-2651C21.3). We first investigated the region and then looked for association between the genes in the region and FBN1 expression. For the first time, we showed that the SLN gene is weakly expressed in skin fibroblasts. There is no direct correlation between SLN and FBN1 gene expression. We showed that calcium influx modulates FBN1 gene expression. Finally, SLN gene expression is highly correlated to that of the neighboring SNX7-ps1. We were able to confirm the impact of calcium influx on FBN1 gene expression but we could not conclude regarding the role of sarcolipin and/or the eQTL locus in this regulation.

Analysis of Polymorphisms in the Mediator Complex Subunit 13-like (Med13L) Gene in the Context of Immune Function and Development of Experimental Arthritis

The Mediator complex subunit 13-like (MED13L) protein is part of the multi-protein mediator complex and plays an important role in gene transcription. Polymorphisms in the MED13L gene have been linked to congenital heart anomalies and intellectual disabilities. Despite recent evidence of indirect links of MED13L to cytokine release and inflammation, impact of genetic variations in MED13L on immune cells remains unexplored. The B10.RIII and RIIIS/J mouse strains vary in susceptibility to induced experimental autoimmune disease models. From sequencing data of the two mouse strains, we identified six polymorphisms in the coding regions of Med13L. Using congenic mice, we studied the effect of these polymorphisms on immune cell development and function along with susceptibility to collagen-induced arthritis, an animal model for rheumatoid arthritis. Combining in vivo disease data, in vitro functional data, and computational analysis of the reported non-synonymous polymorphisms, we report that genetic polymorphisms in Med13L do not affect the immune phenotype in these mice and are predicted to be non-disease associated.

PTEN/PTENP1: 'Regulating the regulator of RTK-dependent PI3K/Akt signalling', new targets for cancer therapy

Regulation of the PI-3 kinase (PI3K)/Akt signalling pathway is essential for maintaining the integrity of fundamental cellular processes, cell growth, survival, death and metabolism, and dysregulation of this pathway is implicated in the development and progression of cancers. Receptor tyrosine kinases (RTKs) are major upstream regulators of PI3K/Akt signalling. The phosphatase and tensin homologue (PTEN), a well characterised tumour suppressor, is a prime antagonist of PI3K and therefore a negative regulator of this pathway. Loss or inactivation of PTEN, which occurs in many tumour types, leads to overactivation of RTK/PI3K/Akt signalling driving tumourigenesis. Cellular PTEN levels are tightly regulated by a number of transcriptional, post-transcriptional and post-translational regulatory mechanisms. Of particular interest, transcription of the PTEN pseudogene, PTENP1, produces sense and antisense transcripts that exhibit post-transcriptional and transcriptional modulation of PTEN expression respectively. These additional levels of regulatory complexity governing PTEN expression add to the overall intricacies of the regulation of RTK/PI-3 K/Akt signalling. This review will discuss the regulation of oncogenic PI3K signalling by PTEN (the regulator) with a focus on the modulatory effects of the sense and antisense transcripts of PTENP1 on PTEN expression, and will further explore the potential for new therapeutic opportunities in cancer treatment.

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.

Deep sequencing of natural and experimental populations of Drosophila melanogaster reveals biases in the spectrum of new mutations

Mutations provide the raw material of evolution, and thus our ability to study evolution depends fundamentally on having precise measurements of mutational rates and patterns. We generate a data set for this purpose using (1) de novo mutations from mutation accumulation experiments and (2) extremely rare polymorphisms from natural populations. The first, mutation accumulation (MA) lines are the product of maintaining flies in tiny populations for many generations, therefore rendering natural selection ineffective and allowing new mutations to accrue in the genome. The second, rare genetic variation from natural populations allows the study of mutation because extremely rare polymorphisms are relatively unaffected by the filter of natural selection. We use both methods in Drosophila melanogaster, first generating our own novel data set of sequenced MA lines and performing a meta-analysis of all published MA mutations (∼2000 events) and then identifying a high quality set of ∼70,000 extremely rare (≤0.1%) polymorphisms that are fully validated with resequencing. We use these data sets to precisely measure mutational rates and patterns. Highlights of our results include: a high rate of multinucleotide mutation events at both short (∼5 bp) and long (∼1 kb) genomic distances, showing that mutation drives GC content lower in already GC-poor regions, and using our precise context-dependent mutation rates to predict long-term evolutionary patterns at synonymous sites. We also show that de novo mutations from independent MA experiments display similar patterns of single nucleotide mutation and well match the patterns of mutation found in natural populations.

rRNA Pseudogenes in Filamentous Ascomycetes as Revealed by Genome Data

The nuclear ribosomal DNA (rDNA) is considered as a paradigm of concerted evolution. Components of the rDNA tandem repeats (45S) are widely used in phylogenetic studies of different organisms and the internal transcribed spacer (ITS) region was recently selected as a fungal DNA bar code. However, rRNA pseudogenes, as one kind of escape from concerted evolution, were reported in a wide range of organisms, especially in plants and animals. Moreover, large numbers of 5S rRNA pseudogenes were identified in several filamentous ascomycetes. To study whether rDNA evolves in a strict concerted manner and test whether rRNA pseudogenes exist in more species of ascomycetes, intragenomic rDNA polymorphisms were analyzed using whole genome sequences. Divergent rDNA paralogs were found to coexist within a single genome in seven filamentous ascomycetes examined. A great number of paralogs were identified as pseudogenes according to the mutation and secondary structure analyses. Phylogenetic analyses of the three rRNA coding regions of the 45S rDNA repeats, i.e., 18S, 5.8S, and 28S, revealed an interspecies clustering pattern of those different rDNA paralogs. The identified rRNA pseudogenic sequences were validated using specific primers designed. Mutation analyses revealed that the repeat-induced point (RIP) mutation was probably responsible for the formation of those rRNA pseudogenes.

The pseudogene-derived long noncoding RNA SFTA1P is down-regulated and suppresses cell migration and invasion in lung adenocarcinoma

Pseudogenes were once considered to be genomic fossils without biological function. Interestingly, recent evidence showed that a lot of pseudogenes are transcribed in human cancers, and their alterations contribute to multiple cancer development and progression. It is apparent that many pseudogenes transcribe noncoding RNAs and contribute to the role noncoding genome plays in human cancers. On this basis, some pseudogene transcripts are currently ranked among the classes of long noncoding RNAs. In this study, we identified a new pseudogene-derived long noncoding RNA termed SFTA1P by analyzing the microarray data of non-small cell lung cancer from Gene Expression Omnibus datasets. We found that SFTA1P expression was significantly decreased in non-small cell lung cancer tissues compared with normal tissues in non-small cell lung cancer microarray data. Moreover, decreased SFTA1P expression is only correlated with lung adenocarcinoma patients' poor survival time but not with lung squamous cell carcinoma patients' survival. In addition, gain-of-function studies including growth curves, migration, invasion assays, and in vivo studies were performed to verify the tumor suppressor role of SFTA1P in non-small cell lung cancer. Finally, the potential underlying pathways involved in SFTA1P were investigated by analyzing the SFTA1P-correlated genes in The Cancer Genome Atlas lung adenocarcinoma and normal tissues RNA sequencing data. Taken together, these findings demonstrate that pseudogene-derived long noncoding RNA SFTA1P exerts the tumor suppressor functions in human lung adenocarcinoma. Our investigation reveals the novel roles of pseudogene in lung adenocarcinoma, which may serve as a new target for lung adenocarcinoma diagnosis and therapy.

The Pseudogene DUXAP8 Promotes Non-small-cell Lung Cancer Cell Proliferation and Invasion by Epigenetically Silencing EGR1 and RHOB

Recently, the non-protein-coding functional elements in the human genome have been identified as key regulators in postgenomic biology, and a large number of pseudogenes as well as long non-coding RNAs (lncRNAs) have been found to be transcribed in multiple human cancers. However, only a small proportion of these pseudogenes has been functionally characterized. In this study, we screened for pseudogenes associated with human non-small-cell lung cancer (NSCLC) by comparative analysis of several independent datasets from the GEO. We identified a transcribed pseudogene named DUXAP8 that is upregulated in tumor tissues. Patients with higher DUXAP8 expression exhibited shorter survival, suggesting DUXAP8 as a new candidate prognostic marker for NSCLC patients. Knockdown of DUXAP8 impairs cell growth, migration, and invasion, and induces apoptosis both in vitro and in vivo. Mechanistically, DUXAP8 represses the tumor suppressors EGR1 and RHOB by recruiting histone demethylase LSD1 and histone methyltransferase EZH2, thereby promoting cell proliferation, migration, and invasion. These findings indicate that the pseudogene DUXAP8 may act as an oncogene in NSCLC by silencing EGR1 and RHOB transcription by binding with EZH2 and LSD1, which may offer a novel therapeutic target for this disease.

The Effect of Botulinum Toxin Type A on Expression Profiling of Long Noncoding RNAs in Human Dermal Fibroblasts

Objective. This study was aimed at analyzing the expressions of long noncoding RNAs (lncRNAs) in Botulinum Toxin Type A (BoNTA) treated human dermal fibroblasts (HDFs) in vitro. Methods. We used RNA sequencing to characterize the lncRNAs and mRNAs transcriptome in the control and BoNTA treated group, in conjunction with application of GO (gene ontology) analysis and KEGG (kyoto encyclopedia of genes and genomes) analysis to delineate the alterations in gene expression. We also obtained quantitative real time polymerase chain reaction (qRT-PCR) to confirm some differentially expressed genes. Results. Numerous differentially expressed genes were observed by microarrays between the two groups. qRT-PCR confirmed the changes of six lncRNAs (RP11-517C16.2-001, FR271872, LOC283352, RP11-401E9.3, FGFR3P, and XXbac-BPG16N22.5) and nine mRNAs (NOS2, C13orf15, FOS, FCN2, SPINT1, PLAC8, BIRC5, NOS2, and COL19A1). Farther studies indicated that the downregulating effect of BoNTA on the expression of FGFR3P was time-related and the dosage of BoNTA at a range from 2.5 U/10⁶ cells to 7.5 U/10⁶ cells increased the expression of FGFR3P and COL19A1 in HDFs as well. Conclusion. The expression profiling of lncRNAs was visibly changed in BoNTA treated HDFs. Further studies should focus on several lncRNAs to investigate their functions in BoNTA treated HDFs and the underlying mechanisms.

An XIST-related small RNA regulates KRAS G-quadruplex formation beyond X-inactivation

X-inactive-specific transcript (XIST), a long non-coding RNA, is essential for the initiation of X-chromosome inactivation. However, little is known about other roles of XIST in the physiological process in eukaryotic cells. In this study, the bioinformatics approaches revealed XIST could be processed into a small non-coding RNA XPi2. The XPi2 RNA was confirmed by a northern blot assay; its expression was gender-independent, suggesting the role of XPi2 was beyond X-chromosome inactivation. The pull-down assay combined with LC-MS-MS identified two XPi2-associated proteins, nucleolin and hnRNP A1, connected to the formation of G-quadruplex. Moreover, the microarray data showed the knockdown of XPi2 down-regulated the KRAS pathway. Consistently, we tested the expression of ten genes, including KRAS, which was correlated with a G-quadruplex formation and found the knockdown of XPi2 caused a dramatic decrease in the transcription level of KRAS among the ten genes. The results of CD/NMR assay also supported the interaction of XPi2 and the polypurine-polypyrimidine element of KRAS. Accordingly, XPi2 may stimulate the KRAS expression by attenuating G-quadruplex formation. Our present work sheds light on the novel role of small RNA XPi2 in modulating the G-quadruplex formation which may play some essential roles in the KRAS- associated carcinogenesis.

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.

Expression of evolutionarily novel genes in tumors

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

Computational Identification of the Paralogs and Orthologs of Human Cytochrome P450 Superfamily and the Implication in Drug Discovery

The human cytochrome P450 (CYP) superfamily consisting of 57 functional genes is the most important group of Phase I drug metabolizing enzymes that oxidize a large number of xenobiotics and endogenous compounds, including therapeutic drugs and environmental toxicants. The CYP superfamily has been shown to expand itself through gene duplication, and some of them become pseudogenes due to gene mutations. Orthologs and paralogs are homologous genes resulting from speciation or duplication, respectively. To explore the evolutionary and functional relationships of human CYPs, we conducted this bioinformatic study to identify their corresponding paralogs, homologs, and orthologs. The functional implications and implications in drug discovery and evolutionary biology were then discussed. GeneCards and Ensembl were used to identify the paralogs of human CYPs. We have used a panel of online databases to identify the orthologs of human CYP genes: NCBI, Ensembl Compara, GeneCards, OMA ("Orthologous MAtrix") Browser, PATHER, TreeFam, EggNOG, and Roundup. The results show that each human CYP has various numbers of paralogs and orthologs using GeneCards and Ensembl. For example, the paralogs of CYP2A6 include CYP2A7, 2A13, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 2F1, 2J2, 2R1, 2S1, 2U1, and 2W1; CYP11A1 has 6 paralogs including CYP11B1, 11B2, 24A1, 27A1, 27B1, and 27C1; CYP51A1 has only three paralogs: CYP26A1, 26B1, and 26C1; while CYP20A1 has no paralog. The majority of human CYPs are well conserved from plants, amphibians, fishes, or mammals to humans due to their important functions in physiology and xenobiotic disposition. The data from different approaches are also cross-validated and validated when experimental data are available. These findings facilitate our understanding of the evolutionary relationships and functional implications of the human CYP superfamily in drug discovery.

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

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

Non-functional plastid ndh gene fragments are present in the nuclear genome of Norway spruce (Picea abies L. Karsch): insights from in silico analysis of nuclear and organellar genomes

Many genes have been lost from the prokaryote plastidial genome during the early events of endosymbiosis in eukaryotes. Some of them were definitively lost, but others were relocated and functionally integrated to the host nuclear genomes through serial events of gene transfer during plant evolution. In gymnosperms, plastid genome sequencing has revealed the loss of ndh genes from several species of Gnetales and Pinaceae, including Norway spruce (Picea abies). This study aims to trace the ndh genes in the nuclear and organellar Norway spruce genomes. The plastid genomes of higher plants contain 11 ndh genes which are homologues of mitochondrial genes encoding subunits of the proton-pumping NADH-dehydrogenase (nicotinamide adenine dinucleotide dehydrogenase) or complex I (electron transport chain). Ndh genes encode 11 NDH polypeptides forming the Ndh complex (analogous to complex I) which seems to be primarily involved in chloro-respiration processes. We considered ndh genes from the plastidial genome of four gymnosperms (Cryptomeria japonica, Cycas revoluta, Ginkgo biloba, Podocarpus totara) and a single angiosperm species (Arabidopsis thaliana) to trace putative homologs in the nuclear and organellar Norway spruce genomes using tBLASTn to assess the evolutionary fate of ndh genes in Norway spruce and to address their genomic location(s), structure, integrity and functionality. The results obtained from tBLASTn were subsequently analyzed by performing homology search for finding ndh specific conserved domains using conserved domain search. We report the presence of non-functional plastid ndh gene fragments, excepting ndhE and ndhG genes, in the nuclear genome of Norway spruce. Regulatory transcriptional elements like promoters, TATA boxes and enhancers were detected in the upstream regions of some ndh fragments. We also found transposable elements in the flanking regions of few ndh fragments suggesting nuclear rearrangements in those regions. These evidences support the hypothesis that, at least in Picea, ndh translocations from the plastid to the nuclear genome have occurred, and that there might have been a functional machinery at some time during evolution to accommodate them within a nuclear-encoded environment, or attempts to form it.

Pseudogene-expressed RNAs: a new frontier in cancers

Over the past decade, the importance of non-protein-coding functional elements in the human genome has emerged from the water and been identified as a key revelation in post-genomic biology. Since the completion of the ENCODE (Encyclopedia of DNA Elements) and FANTOM (Functional Annotation of Mammals) project, tens of thousands of pseudogenes as well as numerous long non-coding RNA (lncRNA) genes were identified. However, while pseudogenes were initially regarded as non-functional relics littering the human genome during evolution, recent studies have revealed that they play critical roles at multiple levels in diverse physiological and pathological processes, especially in cancer through parental-gene-dependent or parental-gene-independent regulation. Herein, we review the current knowledge of pseudogenes and synthesize the nascent evidence for functional properties and regulatory modalities exerted by pseudogene-transcribed RNAs in human cancers and prospect the potential as molecular signatures in cancer reclassification and tailored therapy.

Functions of noncoding sequences in mammalian genomes

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

Dissecting the genetics of the human transcriptome identifies novel trait-related trans-eQTLs and corroborates the regulatory relevance of non-protein coding loci†

Genetics of gene expression (eQTLs or expression QTLs) has proved an indispensable tool for understanding biological pathways and pathomechanisms of trait-associated SNPs. However, power of most genome-wide eQTL studies is still limited. We performed a large eQTL study in peripheral blood mononuclear cells of 2112 individuals increasing the power to detect trans-effects genome-wide. Going beyond univariate SNP-transcript associations, we analyse relations of eQTLs to biological pathways, polygenetic effects of expression regulation, trans-clusters and enrichment of co-localized functional elements. We found eQTLs for about 85% of analysed genes, and 18% of genes were trans-regulated. Local eSNPs were enriched up to a distance of 5 Mb to the transcript challenging typically implemented ranges of cis-regulations. Pathway enrichment within regulated genes of GWAS-related eSNPs supported functional relevance of identified eQTLs. We demonstrate that nearest genes of GWAS-SNPs might frequently be misleading functional candidates. We identified novel trans-clusters of potential functional relevance for GWAS-SNPs of several phenotypes including obesity-related traits, HDL-cholesterol levels and haematological phenotypes. We used chromatin immunoprecipitation data for demonstrating biological effects. Yet, we show for strongly heritable transcripts that still little trans-chromosomal heritability is explained by all identified trans-eSNPs; however, our data suggest that most cis-heritability of these transcripts seems explained. Dissection of co-localized functional elements indicated a prominent role of SNPs in loci of pseudogenes and non-coding RNAs for the regulation of coding genes. In summary, our study substantially increases the catalogue of human eQTLs and improves our understanding of the complex genetic regulation of gene expression, pathways and disease-related processes.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

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

Pseudogene-derived lncRNAs: emerging regulators of gene expression

In the more than one decade since the completion of the Human Genome Project, the prevalence of non-protein-coding functional elements in the human genome has emerged as a key revelation in post-genomic biology. Highlighted by the ENCODE (Encyclopedia of DNA Elements) and FANTOM (Functional Annotation of Mammals) consortia, these elements include tens of thousands of pseudogenes, as well as comparably numerous long non-coding RNA (lncRNA) genes. Pseudogene transcription and function remain insufficiently understood. However, the field is of great importance for human disease due to the high sequence similarity between pseudogenes and their parental protein-coding genes, which generates the potential for sequence-specific regulation. Recent case studies have established essential and coordinated roles of both pseudogenes and lncRNAs in development and disease in metazoan systems, including functional impacts of lncRNA transcription at pseudogene loci on the regulation of the pseudogenes’ parental genes. This review synthesizes the nascent evidence for regulatory modalities jointly exerted by lncRNAs and pseudogenes in human disease, and for recent evolutionary origins of these systems.

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.

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

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