Search

BackgroundThe human genome encodes over 14,000 pseudogenes that are evolutionary relics of protein-coding genes and commonly considered as nonfunctional. Emerging evidence suggests that some pseudogenes may exert important functions. However, to what extent human pseudogenes are functionally relevant remains unclear. There has been no large-scale characterization of pseudogene function because of technical challenges, including high sequence similarity between pseudogene and parent genes, and poor annotation of transcription start sites.ResultsTo overcome these technical obstacles, we develop an integrated computational pipeline to design the first genome-wide library of CRISPR interference (CRISPRi) single-guide RNAs (sgRNAs) that target human pseudogene promoter-proximal regions. We perform the first pseudogene-focused CRISPRi screen in luminal A breast cancer cells and reveal approximately 70 pseudogenes that affect breast cancer cell fitness. Among the top hits, we identify a cancer-testis unitary pseudogene, MGAT4EP, that is predominantly localized in the nucleus and interacts with FOXA1, a key regulator in luminal A breast cancer. By enhancing the promoter binding of FOXA1, MGAT4EP upregulates the expression of oncogenic transcription factor FOXM1. Integrative analyses of multi-omic data from the Cancer Genome Atlas (TCGA) reveal many unitary pseudogenes whose expressions are significantly dysregulated and/or associated with overall/relapse-free survival of patients in diverse cancer types.ConclusionsOur study represents the first large-scale study characterizing pseudogene function. Our findings suggest the importance of nuclear function of unitary pseudogenes and underscore their underappreciated roles in human diseases. The functional genomic resources developed here will greatly facilitate the study of human pseudogene function.

CpG islands (CGIs) are aggregation of CpG dinucleotides in the promoters of mammalian genes. These CGIs are present in almost all the housekeeping genes and some tissue-specific genes in the mammalian genome. Extensive research has been done on the prevalence and role of CGIs in protein-coding genes. However, little is known about CGIs in pseudogenes. In the current research project, we focused on CGIs in three main classes of pseudogenes e.g., duplicated pseudogenes (DPGs), processed pseudogenes (PPGs), and unitary pseudogenes (UPGs). We discovered a predominant absence of CGIs in the promoters of all three pseudogenes. We also compared the CGI profile of these pseudogenes with their parent genes and found that unitary pseudogenes (UPGs) differ from the DPGs and PPGs in the sense that in the latter, lack of CGIs is a consequential event while in UPGs, this lack of CGIs in their promoters is not a result of pseudogenization process. We also discussed the implication of the results obtained from this comparison. To our knowledge, this is the first-ever study highlighting this aspect of UPGs throwing new insights into the evolution of genome in general and especially in the context of pseudogenes.

Prokaryotic genomes are generally gene dense and encode relatively few pseudogenes, or nonfunctional/inactivated remnants of genes. However, in certain contexts, such as recent ecological shifts or extreme population bottlenecks (such as those experienced by symbionts and pathogens), pseudogenes can quickly accumulate and form a substantial fraction of the genome. Identification of pseudogenes is, thus, a critical step for understanding the evolutionary forces acting upon, and the functional potential encoded within, prokaryotic genomes. Here, we present Pseudofinder, an open-source software dedicated to pseudogene identification and analysis. With Pseudofinder’s multi-pronged, reference-based approach, we demonstrate its capacity to detect a wide variety of pseudogenes, including those that are highly degraded and typically missed by gene-calling pipelines, as well newly formed pseudogenes, which can have only one or a few inactivating mutations. Additionally, Pseudofinder can detect intact genes undergoing relaxed selection, which may indicate incipient pseudogene formation. Implementation of Pseudofinder in annotation pipelines will not only clarify the functional potential of sequenced microbes, but will also generate novel insights and hypotheses regarding the evolutionary dynamics of bacterial and archaeal genomes.

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.

LINE-1-mediated retrotransposition of protein-coding mRNAs is an active process in modern humans for both germline and somatic genomes. Prior works that surveyed human data mostly relied on detecting discordant mappings of paired-end short reads, or exon junctions contained in short reads. Moreover, there have been few genome-wide comparisons between gene retrocopies in great apes and humans. In this study, we introduced a more sensitive and accurate method to identify processed pseudogenes. Our method utilizes long-read assemblies, and more importantly, is able to provide full-length retrocopy sequences as well as flanking regions which are missed by short-read based methods. From 22 human individuals, we pinpointed 40 processed pseudogenes that are not present in the human reference genome GRCh38 and identified 17 pseudogenes that are in GRCh38 but absent from some input individuals. This represents a significantly higher discovery rate than previous reports (39 pseudogenes not in the reference genome out of 939 individuals). We also provided an overview of lineage-specific retrocopies in chimpanzee, gorilla, and orangutan genomes.

The junk DNA "pseudogenes," known as genomic fossils, are characterized by their ubiquitousness and abundance within the genomic structure. These genomics sets are recognized by the potential activity of meta-regulating the parent genes; these are transcribed into interfering RNA, consequently acting on miRNA concentration, thereby shedding light on the crosstalk of the pseudogenes' miRNA, siRNA, lncRNA/tumor therapy co-relationship. Moreover, an upcoming visualization regarding pseudogenes is under investigation, which describes the potentiality of pseudogenes as a fundamental component of cancerous evolutionary processing tools. Accordingly, here is a systematic review covering pseudobirth, pseudosignatures, and functional properties of pseudogenes, concluding that these pseudogenes are hypothetically predictive tumor therapies.

As a momentous post-transcriptional regulator, microRNAs (miRNAs) are attracting more and more attention. The classical miRNAs regulated mechanism shows it binds to the targets' 3'UTR thus play the role in post-transcription. Meanwhile, single miRNA can target multiple genes, so those should compete to bind that miRNA. Vice versa, single gene can sponge mass of miRNAs as well. Thus the competitive endogenous RNAs (ceRNAs) hypothesis was put forward in 2011. The ceRNA hypothesis has made huge achievements, in particular in non-coding genes, which including long non-coding RNAs (lncRNAs), circle RNAs (circRNAs) and pseudogenes, even viral transcripts. It also contributed greatly to epigenetics development. However, an increasing number of controversies have occurred with applause. Based on this situation, this review introduces something in detail about the ceRNAs hypothesis achieved in lncRNAs, circRNAs, pseudogenes and viral transcripts, respectively. Meanwhile, it also covers controversy of the ceRNAs hypothesis.

Alternative splicing (AS) and alternative polyadenylation (APA) generate diverse transcripts in mammalian genomes during development and differentiation. Epigenetic marks such as trimethylation of histone H3 lysine 36 (H3K36me3) and DNA methylation play a role in generating transcriptome diversity. Intragenic CpG islands (iCGIs) and their corresponding host genes exhibit dynamic epigenetic and gene expression patterns during development and between different tissues. We hypothesise that iCGI-associated H3K36me3, DNA methylation and transcription can influence host gene AS and/or APA. We investigate H3K36me3 and find that this histone mark is not a major regulator of AS or APA in our model system. Genomewide, we identify over 4000 host genes that harbour an iCGI in the mammalian genome, including both previously annotated and novel iCGI/host gene pairs. The transcriptional activity of these iCGIs is tissue- and developmental stage-specific and, for the first time, we demonstrate that the premature termination of host gene transcripts upstream of iCGIs is closely correlated with the level of iCGI transcription in a DNA-methylation independent manner. These studies suggest that iCGI transcription, rather than H3K36me3 or DNA methylation, interfere with host gene transcription and pre-mRNA processing genomewide and contributes to the spatiotemporal diversification of both the transcriptome and proteome.

The protein-coding genes and pseudogenes of Cuscuta australis had the diverse contribution to the formation and evolution of parasitism. The codon usage pattern analysis of these two type genes could be used to understand the gene transcription and translation. In this study, we systematically analyzed the codon usage patterns of protein-coding sequences and pseudogenes sequences in C. australis. The results showed that the high frequency codons of protein coding sequences and pseudogenes had the same A/U bias in the third position. However, these two sequences had converse bias at the third base in optimal codons: the protein coding sequences preferred G/C-ending codons while pseudogene sequences preferred A/U-ending codons. Neutrality plot and effective number of codons plot revealed that natural selection played a more important role than mutation pressure in two sequences codon usage bias. Furthermore, the gene expression level had a significant positive correlation with codon usage bias in C. australis. Highly-expressed protein coding genes exhibited a higher codon bias than lowly-expressed genes. Meanwhile, the high-expression genes tended to use G/C-ending synonymous codons. This result further verified the optimal codons usage bias and its correlation with the gene expression in C. australis.

The Ptr1 (Pseudomonas tomato race 1) locus in Solanum lycopersicoides confers resistance to strains of Pseudomonas syringae pv. tomato expressing AvrRpt2 and Ralstonia pseudosolanacearum expressing RipBN. Here we describe the identification and phylogenetic analysis of the Ptr1 gene. A single recombinant among 585 F2 plants segregating for the Ptr1 locus was discovered that narrowed the Ptr1 candidates to eight nucleotide-binding leucine-rich repeat protein (NLR)-encoding genes. From analysis of the gene models in the S. lycopersicoides genome sequence and RNA-Seq data, two of the eight genes emerged as the strongest candidates for Ptr1. One of these two candidates was found to encode Ptr1 based on its ability to mediate recognition of AvrRpt2 and RipBN when it was transiently expressed with these effectors in leaves of Nicotiana glutinosa. The ortholog of Ptr1 in tomato and in Solanum pennellii is a pseudogene. However, a functional Ptr1 ortholog exists in Nicotiana benthamiana and potato, and both mediate recognition of AvrRpt2 and RipBN. In apple and Arabidopsis, recognition of AvrRpt2 is mediated by the Mr5 and RPS2 proteins, respectively. Phylogenetic analysis places Ptr1 in a distinct clade compared with Mr5 and RPS2, and it therefore appears to have arisen by convergent evolution for recognition of AvrRpt2.

Monoclonal antibodies (mAbs) are widely utilized as therapeutic drugs for various diseases, such as cancer, autoimmune diseases, and infectious diseases. Using the avian-derived B cell line DT40, we previously developed an antibody display technology, namely, the ADLib system, which rapidly generates antigen-specific mAbs. Here, we report the development of a human version of the ADLib system and showcase the streamlined generation and optimization of functional human mAbs. Tailored libraries were first constructed by replacing endogenous immunoglobulin genes with designed human counterparts. From these libraries, clones producing full-length human IgGs against distinct antigens can be isolated, as exemplified by the selection of antagonistic mAbs. Taking advantage of avian biology, effective affinity maturation was achieved in a straightforward manner by seamless diversification of the parental clones into secondary libraries followed by single-cell sorting, quickly affording mAbs with improved affinities and functionalities. Collectively, we demonstrate that the human ADLib system could serve as an integrative platform with unique diversity for rapid de novo generation and optimization of therapeutic or diagnostic antibody leads. Furthermore, our results suggest that libraries can be constructed by introducing exogenous genes into DT40 cells, indicating that the ADLib system has the potential to be applied for the rapid and effective directed evolution and optimization of proteins in various fields beyond biomedicine.

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.

Pseudogenes have been known as non-functional molecular relics developed from inactivated genomic mutations while evolution. However, they have a various function at distinct molecular stages (DNA, RNA and protein) in different biological events including in cancer and associated with parental-gene-dependently and parental-gene independently. The interactions of pseudogenes with ancestor genes or other genes alter in their sequences and transcriptional processes. Pseudogene of RNA sequences have multiple functions in post-transcriptional activities as antisense RNAs, endogenous small-interference RNAs, and competing endogenous RNAs. Moreover, it also plays vital roles in controlling of its parent genes and other pseudogenes transcribing into RNA. Pseudogene transcripts make small interfering RNA or decline cellular miRNA level. The current review focuses on pseudogene functional signature in human genome by regulating the gene expression.

Head and neck squamous cell carcinoma (HNSCC) is the sixth malignancy, which is characterized by poor prognosis or high mortality because of the lack of predicting markers. Aberrant cancer pseudogenes have been found predictive for prognosis. We aim to identify a pseudogene-based prognosis signature for HNSCC by machine learning. RNA-seq data were downloaded from The Cancer Genome Atlas, and 700 differentially-expressed pseudogenes were identified. The survival-related pseudogenes were screened through COX-regression analysis, which includes univariate regression, least absolute shrinkage and selection operator regression, and multivariate regression, and a five-pseudogene signature was constructed. The value of prediction for the signature was validated in multiple subgroups in terms of survival. Gene set enrichment analysis (GSEA) and coexpression analysis were used to determine the underlying biological functions. Seven hundred dysregulated pseudogenes were identified, and the five-pseudogene signature can distinguish the low-risk and high-risk patients for both training and testing sets and predicted prognosis with high sensitivity and specificity. Furthermore, the signature was applicable to patients of different genders, ages, stages, and grades. Coexpression analysis revealed that the five-pseudogene is associated with immune system. GSEA showed cancer-related biological process and pathways the five-pseudogene involved in. The five-pseudogene signature is not only a novel marker for prognosis but also a promising signature for monitoring therapeutic schedule. Therefore, our findings may have potential clinical significance.

Pseudogenes are mutated copies of protein-coding genes that cannot be translated into proteins, but a small subset of pseudogenes has been detected at the protein level. Although ubiquitin pseudogenes represent one of the most abundant pseudogene families in many organisms, little is known about their expression and signaling potential. By re-analyzing public RNA-sequencing and proteomics datasets, we here provide evidence for the expression of several ubiquitin pseudogenes including UBB pseudogene 4 (UBBP4), which encodes UbKEKS (Q2K, K33E, Q49K, N60S). The functional consequences of UbKEKS conjugation appear to differ from canonical ubiquitylation. Quantitative proteomics shows that UbKEKS modifies specific proteins including lamins. Knockout of UBBP4 results in slower cell division, and accumulation of lamin A within the nucleolus. Our work suggests that a subset of proteins reported as ubiquitin targets may instead be modified by ubiquitin variants that are the products of wrongly annotated pseudogenes and induce different functional effects., Ubiquitin pseudogenes are present in many organisms but whether they encode functional proteins has remained unclear. Here, the authors show that human UBB pseudogene 4 produces ubiquitin variants with amino acid compositions and cellular functions that are distinct from canonical ubiquitin.

CACX is one of the most common cancer affecting women world-wide. Here, expression microarray analysis revealed 8 over-expressed transcribed pseudogenes ( GBP1P1 , HLA-DRB6 , HLA H , SLC6A10P , NAPSB, KRT16P2 , PTTG3P and RNF126P1), down-regulated 7 lincRNAs ( H19 , MIR100HG , MEG3 , DIO3OS , HOXA11-AS , CD27-AS1 and EPB41L4A-AS ) and 6 snoRNAs ( SNORD97 , SNORD3A , SNORD3C , SNORD3D , SNORA12 and SCARNA9 ) as DEncGs (log2 fold-change ≥ ±1.0) in CACX. Consequently, down-regulation of lincRNA MEG3 and over-expression of pseudogenes, GBP1P1 and PTTG3P in the microarray analysis were found concordant with the real-time quantitative PCR results upon validation. Then, Ingenuity® Pathway analysis (IPA®) analysis with deregulated DEncGs identified functionally important gene, H19 . Further, validation ( n = 52) of expression confirmed frequent downregulation of H19 with significant association with its deletion (LOH) and promoter methylation ( n = 128) in CACX. Moreover, clinico-pathological analysis found Indian CACX patients ( n = 26) with alterations of H19 by deletion or, promoter methylation with concomitant low expression have poor prognosis.

Objective: Ovarian cancer is one of the most common and lethal cancer types in women. The molecular mechanism of ovarian cancer progression is still unclear. Results: Here, we first reported that expression levels of three genes, GJB2, S100A2 and SPOCK2, were significantly higher in advanced stage than that in early stage of ovarian cancer, and upregulation of them indicated poor prognosis of patients with ovarian cancer. Subsequently, 8, 6 and 20 miRNAs were predicted to target GJB2, S100A2 and SPOCK2, respectively. Among these miRNA-mRNA pairs, hsa-miR-363-3p-SPOCK2 axis was the most potential in suppressing progression of ovarian cancer. Mechanistically, we found that hsa-miR-363-3p-SPOCK2 axis was involved in regulation of actin cytoskeleton. Moreover, 6 pseudogenes and 8 lncRNAs were identified to potentially inhibit hsa-miR-363-3p-SPOCK2 axis in ovarian cancer. Conclusions: Collectively, we elucidate a regulatory role of pseudogene/lncRNA-hsa-miR-363-3p-SPOCK2 pathway in progression of ovarian cancer, which may provide effective therapeutic approaches and promising prognostic biomarkers for ovarian cancer. Materials and methods: Differentially expressed genes (DEGs) in ovarian cancer were first screened using {"type":"entrez-geo","attrs":{"text":"GSE12470","term_id":"12470"}}GSE12470, after which DEGs expression were validated using GEPIA. Kaplan-Meier analysis was employed to assess the prognostic values. Potential miRNAs were predicted by seven target prediction databases, and upstream lncRNAs and pseudogenes of hsa-miR-363-3p were forecasted through miRNet or starBase. UALCAN and starBase were used to obtain the co-expressed genes of SPOCK. Enrichment analysis for these co-expressed genes was performed by Enrichr.

Summary Mutations in the NADPH oxidase, which is crucial for the respiratory burst in phagocytes, result in chronic granulomatous disease (CGD). The only curative treatment option for CGD patients, who suffer from severe infections, is allogeneic bone marrow transplantation. Over 90% of patients with mutations in the p47phox subunit of the oxidase complex carry the deletion c.75_76delGT (ΔGT). This frequent mutation most likely originates via gene conversion from one of the two pseudogenes NCF1B or NCF1C, which are highly homologous to NCF1 (encodes p47phox) but carry the ΔGT mutation. We applied CRISPR/Cas9 to generate patient-like p47-ΔGT iPSCs for disease modeling. To avoid unpredictable chromosomal rearrangements by CRISPR/Cas9-mediated cleavage in the pseudogenes, we developed a gene-correction approach to specifically target NCF1 but leave the pseudogenes intact. Functional assays revealed restored NADPH oxidase activity and killing of bacteria in corrected phagocytes as well as the specificity of this approach., Graphical Abstract, Highlights • CRISPR/Cas9-mediated disease modeling of p47-CGD in iPSCs mimics pathogenesis • Targeted repair allows correction of p47-CGD without cleavage in the pseudogenes • Corrected iPSC-derived granulocytes have a restored NADPH oxidase function • Corrected iPSC-derived macrophages phagocytose and kill bacteria, Klatt, Schambach and colleagues generated patient-like p47-ΔGT iPSCs for disease modeling of chronic granulomatous disease and subsequent gene correction of NCF1 (encodes p47phox) using CRISPR/Cas9. Their strategy avoided cleavage of the highly homologous NCF1 pseudogenes, which might have a functional role. Corrected granulocytes and macrophages demonstrated restored NADPH oxidase activity, formation of neutrophil extracellular traps, and killing of phagocytosed bacteria.

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle. This article is protected by copyright. All rights reserved.

Simple Summary Competing endogenous RNAs (ceRNAs) have gained attention in cancer research owing to their involvement in microRNA-mediated gene regulation. Here, we identified a shared ceRNA network across five hormone-dependent (HD) cancers (prostate, breast, colon, rectal, and endometrial), that contain two long non-coding RNAs, nine mRNAs, and seventy-four microRNAs. Among them, two mRNAs and forty-one microRNAs were associated with at least one HD cancer survival. A similar analytical approach can be applied to identify shared ceRNAs across a group of related cancers, which will significantly contribute to understanding their shared disease biology. Abstract The discovery of microRNAs (miRNAs) has fundamentally transformed our understanding of gene regulation. The competing endogenous RNA (ceRNA) hypothesis postulates that messenger RNAs and other RNA transcripts, such as long non-coding RNAs and pseudogenes, can act as natural miRNA sponges. These RNAs influence each other’s expression levels by competing for the same pool of miRNAs through miRNA response elements on their target transcripts, thereby modulating gene expression and protein activity. In recent years, these ceRNA regulatory networks have gained considerable attention in cancer research. Several studies have identified cancer-specific ceRNA networks. Nevertheless, prior bioinformatic analyses have focused on long-non-coding RNA-associated ceRNA networks. Here, we identify an extended ceRNA network (including both long non-coding RNAs and pseudogenes) shared across a group of five hormone-dependent (HD) cancers, i.e., prostate, breast, colon, rectal, and endometrial cancers, using data from The Cancer Genome Atlas (TCGA). We performed a functional enrichment analysis for differentially expressed genes in the shared ceRNA network of HD cancers, followed by a survival analysis to determine their prognostic ability. We identified two long non-coding RNAs, nine genes, and seventy-four miRNAs in the shared ceRNA network across five HD cancers. Among them, two genes and forty-one miRNAs were associated with at least one HD cancer survival. This study is the first to investigate pseudogene-associated ceRNAs across a group of related cancers and highlights the value of this approach to understanding the shared molecular pathogenesis in a group of related diseases.

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.

The interferon-induced transmembrane (IFITM) gene family performs multiple functions in immunity, including inhibition of virus entry into cells. The IFITM repertoire varies widely between species and consists of protein-coding genes and pseudogenes. The selective forces driving pseudogenization within gene families are rarely understood. In this issue, the human pseudogene IFITM4P is characterized as a virus-induced, long noncoding RNA that contributes to restriction of influenza A virus by regulating mRNA levels of IFITM1, IFITM2, and IFITM3.

Background Pseudogenes are dysfunctional copies of protein-coding genes that showed critical regulatory roles during carcinogenesis. Plasminogen like A (PLGLA) is a transcribed unprocessed pseudogene and biasedly expressed in liver. But its function has not been studied in hepatocellular carcinoma (HCC). Aims We aimed to explore the role of PLGLA in HCC. Methods The expression of PLGLA and its association with pathological features in HCC patients was analyzed using The Cancer Genome Atlas (TCGA) datasets. Quantitative reverse transcription PCR (qRT-PCR) was used to validate PLGLA level in HCC tissue samples and cell lines. Gain-of-function experiments in vitro and in vivo were employed to assess the impact of PLGLA on HCC cell proliferation, migration and invasion. Luciferase reporter assay and RNA pull-down assay were conducted to confirm the interaction among PLGLA, miR-324–3p and GLYATL1. Results PLGLA expression was significantly downregulated in HCC tissues and cell lines. Furthermore, low PLGLA expression was positively associated with tumor progression and poor prognosis. PLGLA restoration markedly suppressed cell proliferation, migration and invasion. Mechanistically, PLGLA could competitively bind to miR-324–3p and acted as a competitive endogenous RNA (ceRNA) to enhance GLYATL1 expression. Conclusions Our results established a novel tumor suppressive role of PLGLA in HCC pathogenesis and highlighted its potential as a therapeutic target for HCC treatment.

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.

Despite more than a century of research, genetic manipulation of Treponema pallidum subsp. pallidum (T. pallidum), the causative agent of syphilis, has not been successful. The lack of genetic engineering tools has severely limited understanding of the mechanisms behind T. pallidum success as a pathogen. A recently described method for in vitro cultivation of T. pallidum, however, has made it possible to experiment with transformation and selection protocols in this pathogen. Here, we describe an approach that successfully replaced the tprA (tp0009) pseudogene in the SS14 T. pallidum strain with a kanamycin resistance (kanR) cassette. A suicide vector was constructed using the pUC57 plasmid backbone. In the vector, the kanR gene was cloned downstream of the tp0574 gene promoter. The tp0574prom-kanR cassette was then placed between two 1-kbp homology arms identical to the sequences upstream and downstream of the tprA pseudogene. To induce homologous recombination and integration of the kanR cassette into the T. pallidum chromosome, in vitro-cultured SS14 strain spirochetes were exposed to the engineered vector in a CaCl2-based transformation buffer and let recover for 24 hours before adding kanamycin-containing selective media. Integration of the kanR cassette was demonstrated by qualitative PCR, droplet digital PCR (ddPCR), and whole-genome sequencing (WGS) of transformed treponemes propagated in vitro and/or in vivo. ddPCR analysis of RNA and mass spectrometry confirmed expression of the kanR message and protein in treponemes propagated in vitro. Moreover, tprA knockout (tprAko-SS14) treponemes grew in kanamycin concentrations that were 64 times higher than the MIC for the wild-type SS14 (wt-SS14) strain and in infected rabbits treated with kanamycin. We demonstrated that genetic manipulation of T. pallidum is attainable. This discovery will allow the application of functional genetics techniques to study syphilis pathogenesis and improve syphilis vaccine development., Author summary Syphilis is still an endemic disease in many low- and middle-income countries, and it has been resurgent in high-income nations for almost two decades. In endemic areas, syphilis causes significant morbidity and mortality, particularly when its causative agent, the spirochete Treponema pallidum subsp. pallidum (T. pallidum) is transmitted to the fetus during pregnancy. A better understanding of T. pallidum biology and syphilis pathogenesis would help devise better control strategies for this infection. One of the limitations associated with working with T. pallidum was our inability to genetically alter this pathogen to evaluate the function of genes encoding virulence factors or create attenuated strains that could be informative for vaccine development when studied using the rabbit model of the disease. Here, we report a transformation protocol that allowed us to replace a specific region of the T. pallidum genome containing a pseudogene (i.e., a non-functional gene) with a stably integrated kanamycin resistance gene. To our knowledge, this is the first-ever report of a method to achieve a genetically modified T. pallidum strain.

BackgroundCanonical nonsense-mediated decay (NMD) is an important splicing-dependent process for mRNA surveillance in mammals. However, processed pseudogenes are not able to trigger NMD due to their lack of introns. It is largely unknown whether they have evolved other surveillance mechanisms.ResultsHere, we find that the RNAs of pseudogenes, especially processed pseudogenes, have dramatically higher m6A levels than their cognate protein-coding genes, associated with de novo m6A peaks and motifs in human cells. Furthermore, pseudogenes have rapidly accumulated m6A motifs during evolution. The m6A sites of pseudogenes are evolutionarily younger than neutral sites and their m6A levels are increasing, supporting the idea that m6A on the RNAs of pseudogenes is under positive selection. We then find that the m6A RNA modification of processed, rather than unprocessed, pseudogenes promotes cytosolic RNA degradation and attenuates interference with the RNAs of their cognate protein-coding genes. We experimentally validate the m6A RNA modification of two processed pseudogenes,DSTNP2andNAP1L4P1, which promotes the RNA degradation of both pseudogenes and their cognate protein-coding genesDSTNandNAP1L4. In addition, the m6A ofDSTNP2regulation of DSTN is partially dependent on the miRNA miR-362-5p.ConclusionsOur discovery reveals a novel evolutionary role of m6A RNA modification in cleaning up the unnecessary processed pseudogene transcripts to attenuate their interference with the regulatory network of protein-coding genes.

BackgroundPseudogenes are non-functional copies of protein coding genes that typically follow a different molecular evolutionary path as compared to functional genes. The inclusion of pseudogene sequences in DNA barcoding and metabarcoding analysis can lead to misleading results. None of the most widely used bioinformatic pipelines used to process marker gene (metabarcode) high throughput sequencing data specifically accounts for the presence of pseudogenes in protein-coding marker genes. The purpose of this study is to develop a method to screen for nuclear mitochondrial DNA segments (nuMTs) in large COI datasets. We do this by: (1) describing gene and nuMT characteristics from an artificial COI barcode dataset, (2) show the impact of two different pseudogene removal methods on perturbed community datasets with simulated nuMTs, and (3) incorporate a pseudogene filtering step in a bioinformatic pipeline that can be used to process Illumina paired-end COI metabarcode sequences. Open reading frame length and sequence bit scores from hidden Markov model (HMM) profile analysis were used to detect pseudogenes.ResultsOur simulations showed that it was more difficult to identify nuMTs from shorter amplicon sequences such as those typically used in metabarcoding compared with full length DNA barcodes that are used in the construction of barcode libraries. It was also more difficult to identify nuMTs in datasets where there is a high percentage of nuMTs. Existing bioinformatic pipelines used to process metabarcode sequences already remove some nuMTs, especially in the rare sequence removal step, but the addition of a pseudogene filtering step can remove up to 5% of sequences even when other filtering steps are in place.ConclusionsOpen reading frame length filtering alone or combined with hidden Markov model profile analysis can be used to effectively screen out apparent pseudogenes from large datasets. There is more to learn from COI nuMTs such as their frequency in DNA barcoding and metabarcoding studies, their taxonomic distribution, and evolution. Thus, we encourage the submission of verified COI nuMTs to public databases to facilitate future studies.

Significance Genome reduction is commonly observed in bacteria of several phyla engaging in obligate nutritional symbioses with insects. In Actinobacteria, however, little is known about the process of genome evolution, despite their importance as prolific producers of antibiotics and their increasingly recognized role as defensive partners of insects and other organisms. Here, we show that “Streptomyces philanthi,” a defensive symbiont of digger wasps, has a G+C-enriched genome in the early stages of erosion, with inactivating mutations in a large proportion of genes, causing dependency on its hosts for certain nutrients, which was validated in axenic symbiont cultures. Additionally, overexpressed catabolic and biosynthetic pathways of the bacteria inside the host indicate host–symbiont metabolic integration for streamlining and control of antibiotic production., Genome erosion is a frequently observed result of relaxed selection in insect nutritional symbionts, but it has rarely been studied in defensive mutualisms. Solitary beewolf wasps harbor an actinobacterial symbiont of the genus Streptomyces that provides protection to the developing offspring against pathogenic microorganisms. Here, we characterized the genomic architecture and functional gene content of this culturable symbiont using genomics, transcriptomics, and proteomics in combination with in vitro assays. Despite retaining a large linear chromosome (7.3 Mb), the wasp symbiont accumulated frameshift mutations in more than a third of its protein-coding genes, indicative of incipient genome erosion. Although many of the frameshifted genes were still expressed, the encoded proteins were not detected, indicating post-transcriptional regulation. Most pseudogenization events affected accessory genes, regulators, and transporters, but “Streptomyces philanthi” also experienced mutations in central metabolic pathways, resulting in auxotrophies for biotin, proline, and arginine that were confirmed experimentally in axenic culture. In contrast to the strong A+T bias in the genomes of most obligate symbionts, we observed a significant G+C enrichment in regions likely experiencing reduced selection. Differential expression analyses revealed that—compared to in vitro symbiont cultures—“S. philanthi” in beewolf antennae showed overexpression of genes for antibiotic biosynthesis, the uptake of host-provided nutrients and the metabolism of building blocks required for antibiotic production. Our results show unusual traits in the early stage of genome erosion in a defensive symbiont and suggest tight integration of host–symbiont metabolic pathways that effectively grants the host control over the antimicrobial activity of its bacterial partner.

Pseudogenes are gene copies presumed to mainly be functionless relics of evolution due to acquired deleterious mutations or transcriptional silencing. When transcribed, pseudogenes may encode proteins or enact RNA-intrinsic regulatory mechanisms. However, the extent, characteristics and functional relevance of the human pseudogene transcriptome are unclear. Short-read sequencing platforms have limited power to resolve and accurately quantify pseudogene transcripts owing to the high sequence similarity of pseudogenes and their parent genes. Using deep full-length PacBio cDNA sequencing of normal human tissues and cancer cell lines, we identify here hundreds of novel transcribed pseudogenes. Pseudogene transcripts are expressed in tissue-specific patterns, exhibit complex splicing patterns and contribute to the coding sequences of known genes. We survey pseudogene transcripts encoding intact open reading frames (ORFs), representing potential unannotated protein-coding genes, and demonstrate their efficient translation in cultured cells. To assess the impact of noncoding pseudogenes on the cellular transcriptome, we delete the nucleus-enriched pseudogene PDCL3P4 transcript from HAP1 cells and observe hundreds of perturbed genes. This study highlights pseudogenes as a complex and dynamic component of the transcriptional landscape underpinning human biology and disease.

Motivation Retrocopies or processed pseudogenes are gene copies resulting from mRNA retrotransposition. These gene duplicates can be fixed, somatically inserted or polymorphic in the genome. However, knowledge regarding unfixed retrocopies (retroCNVs) is still limited, and the development of computational tools for effectively identifying and genotyping them is an urgent need. Results Here, we present sideRETRO, a pipeline dedicated not only to detecting retroCNVs in whole-genome or whole-exome sequencing data but also to revealing their insertion sites, zygosity and genomic context and classifying them as somatic or polymorphic events. We show that sideRETRO can identify novel retroCNVs and genotype them, in addition to finding polymorphic retroCNVs in whole-genome and whole-exome data. Therefore, sideRETRO fills a gap in the literature and presents an efficient and straightforward algorithm to accelerate the study of bona fide retroCNVs. Availability and implementation sideRETRO is available at https://github.com/galantelab/sideRETRO Supplementary information Supplementary data are available at Bioinformatics online.

The eukaryotic ribosomal DNA (rDNA) cluster consists of multiple copies of three genes (18S, 5.8S, and 28S rDNA) and two internal transcribed spacers (ITS1 and ITS2). Typically, rDNA is thought to be evolved in a strict concerted evolution pattern, in which show little or even no variations within species, while sequence diversity between species or higher taxa. However, in the present study, we reported a high level of intra-individual polymorphism of ITS1-5.8S-ITS2 rDNA in the genome of Paraplagusia japonica (Pleuronectiformes: Cynoglossidae), suggesting it conformed a non-concerted evolutionary pattern. Three divergent type sequences (Type A, B, and C) differing in length, GC content, secondary structure stability and minimum free energy coexisted in the genome of each studied individual. According to the particular features of ribosomal pseudogenes, Type B and Type C sequences were postulated as pseudogenes. Expression analyses gave clear evidence that both the functional gene and pseudogenes were abundantly transcribed, indicating expression is not an applicable criterion for identifying ribosomal pseudogenes. These findings will enrich useful data for the studies on rDNA evolution in flatfishes, as well as provide fresh insight on the transcribed ribosomal pseudogenes and further consideration of the definition of ribosomal pseudogene.

Objective 21-hydroxylase deficiency (21-OHD) is the most common cause of congenital adrenal hyperplasia due to CYP21A2 gene mutation . The aim of study is to expand CYP21A2 mutational spectrum in the Chinese population and to provide novel genetic information in terms of ethnic diversity. Design and methods 95 Chinese suspected 21-OHD patients with phenotypes varying from salt-wasting (SW) to nonclassic symptoms were recruited. The clinical characteristics were retrospectively analyzed. Sanger sequencing and multiplex ligation-dependent probe amplification were used to detect point mutations and large gene deletions , respectively. Results 20 different mutant alleles were detected in 35 patients with 21-OHD. The most common variant was c.293-13A/C>G (30.0%), followed by p.I173N (20.0%), large gene conversions (14.3%), large gene deletions (11.4%), and p.R484Pfs*58 (4.3%). Remarkably, we identified a novel F450L variant, in silico predicted to be associated with the salt-wasting form. Two variants including p.R409C and p.R427H, previously considered as conserved in specific ethnicities due to a founder effect, were detected in our cohort. Further, a rare p.H63L + p.V70L variant, hitherto only observed in the Chinese population, in trans with different variants corresponding to the salt-wasting form resulted in diverse phenotypes. Conclusions One novel and four rare variants of CYP21A2 gene corresponding to severe phenotypes were identified in our cohort. Two variants including p.R409C and p.R427H have wider ethnic distributions. Therefore, the sequence of CYP21A2 gene must be analyzed carefully in case rare or novel deleterious variants exist. Our findings improve the understanding of CYP21A2 mutational spectrum in 21-OHD patients and contribute to the precise diagnosis and prenatal counseling.

An exhaustive analysis of a group of closely related parasitic plants shows a predominantly gradual reduction in plastid genome composition and provides the most reduced plastomes in the genus Cuscuta. Parasitic plants have a diminished to completely absent reliance on photosynthesis, and are characterized by sweeping morphological, physiological, and genomic changes. The plastid genome (plastome) is highly conserved in autotrophic plants but is often reduced in parasites, and provides a useful system for documenting the genomic effects of a loss of photosynthesis. Previous studies have shown a substantial degree of heterogeneity in plastome length and composition across the species of the genus Cuscuta. Specifically, species in Cuscuta sect. Ceratophorae were suspected to exhibit even more dynamic plastome evolution than the rest of the genus. This complex of eight closely related species was exhaustively sampled here, and one accession per species was sequenced via a high-throughput approach. Complete plastid genomes were assembled and annotated for each of these species and were found to be 61-87 kbp in length, representing a 45-60% reduction relative to autotrophic Convolvulaceae. The most reduced plastomes on this spectrum have lost the bulk of their photosynthetic genes and are the first fully holoparasitic plastomes described for Cuscuta. The fine-scale nature of the system introduced here allowed us to phylogenetically triangulate the locations of gene loss and pseudogenization events precisely, and to construct a step-by-step model of plastome evolution in these plants. This model reveals an intense burst of gene loss along the branch leading to the most reduced plastomes, and a few idiosyncratic changes elsewhere, allowing us to conclude that the tempo of plastid evolution in sect. Ceratophorae is a blend of gradual and punctuated mode.

Objectives In this study, we aimed to present a strategy for the detection of the RHD pseudogene (RHDψ) based on a real-time polymerase chain reaction (PCR) assay. Background The D-negative phenotype is associated with many genetic alterations. In populations with African ancestry, this phenotype commonly results from the silent variant RHDψ. The evaluation of RHDψ is essential for correct inference of the RhD phenotype in order to avoid false-positive results. Methods We utilised a new method for the simultaneous detection of RHDψ and a fragment from exon 5 of the wild-type RHD gene based on duplex real-time PCR assay. Results The PCR assay allowed specific detection of RHDψ. There was complete agreement between the results generated by the new test and the results generated by molecular analysis performed using end-point PCR methods previously described. Conclusions The assay developed is easy to execute and presents the potential for routine use at blood banks and other associated facilities where it is desired to determine the presence of RHDψ.

Pseudogenes (Ψs), nonfunctional relatives of functional genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Evolutionary analysis provides clues to Ψ origins and effects on gene regulation. However, few systematic studies of plant Ψs have been conducted, hampering comparative analyses. Here, we examined the origin, evolution, and expression patterns of Ψs and their relationships with noncoding sequences in seven angiosperm plants. We identified ∼250,000 Ψs, most of which are more lineage specific than protein-coding genes. The distribution of Ψs on the chromosome indicates that genome recombination may contribute to Ψ elimination. Most Ψs evolve rapidly in terms of sequence and expression levels, showing tissue- or stage-specific expression patterns. We found that a surprisingly large fraction of nontransposable element regulatory noncoding RNAs (microRNAs and long noncoding RNAs) originate from transcription of Ψ proximal upstream regions. We also found that transcription factor binding sites preferentially occur in putative Ψ proximal upstream regions compared with random intergenic regions, suggesting that Ψs have conditioned genome evolution by providing transcription factor binding sites that serve as promoters and enhancers. We therefore propose that rapid rewiring of Ψ transcriptional regulatory regions is a major mechanism driving the origin of novel regulatory modules.

Background The C-C motif chemokine ligand 16 (CCL16) is a potent pro-inflammatory chemokine and a chemoattractant for monocytes and lymphocytes. In normal plasma, it is present at high concentrations and elicits its effects on cells by interacting with cell surface chemokine receptors. In the European rabbit and in rodents such as mouse, rat and guinea pig, CCL16 was identified as a pseudogene, while in the thirteen-lined ground squirrel it appears to be potentially functional. To gain insight into the evolution of this gene in the superorder Glires (rodents and lagomorphs), we amplified the CCL16 gene from eleven Leporidae and seven Ochotonidae species. Results We compared our sequences with CCL16 sequences of twelve rodent species retrieved from public databases. The data show that for all leporid species studied CCL16 is a pseudogene. This is primarily due to mutations at the canonical Cys Cys motif, creating either premature stop codons, or disrupting amino acid replacements. In the Mexican cottontail, CCL16 is pseudogenized due to a frameshift deletion. Additionally, in the exon 1 (signal peptide), there are frameshift deletions present in all leporids studied. In contrast, in Ochotona species, CCL16 is potentially functional, except for an allele in Hoffmann’s pika. In rodents, CCL16 is functional in a number of species, but patterns of pseudogenization similar to those observed in lagomorphs also exist. Conclusions Our results suggest that while functional in the Glires ancestor, CCL16 underwent pseudogenization in some species. This process occurred stochastically or in specific lineages at different moments in the evolution of Glires. These observations suggest that the CCL16 had different evolutionary constrains in the Glires group that could be associated with the CCL16 biological function. Electronic supplementary material The online version of this article (10.1186/s12862-019-1390-7) contains supplementary material, which is available to authorized users.

Background The genus Anaplasma is made up of organisms characterized by small genomes that are undergoing reductive evolution. Anaplasma ovis, one of the seven recognized species in this genus, is an understudied pathogen of sheep and other ruminants. This tick-borne agent is thought to induce only mild clinical disease; however, small deficits may add to larger economic impacts due to the wide geographic distribution of this pathogen. Results In this report we present the first complete genome sequence for A. ovis and compare the genome features with other closely related species. The 1,214,674 bp A. ovis genome encodes 933 protein coding sequences, the split operon arrangement for ribosomal RNA genes, and more pseudogenes than previously recognized for other Anaplasma species. The metabolic potential is similar to other Anaplasma species. Anaplasma ovis has a small repertoire of surface proteins and transporters. Several novel genes are identified. Conclusions Analyses of these important features and significant gene families/genes with potential to be vaccine candidates are presented in a comparative context. The availability of this genome will significantly facilitate research for this pathogen.

Background Trypanosoma conorhini and Trypanosoma rangeli, like Trypanosoma cruzi, are kinetoplastid protist parasites of mammals displaying divergent hosts, geographic ranges and lifestyles. Largely nonpathogenic T. rangeli and T. conorhini represent clades that are phylogenetically closely related to the T. cruzi and T. cruzi-like taxa and provide insights into the evolution of pathogenicity in those parasites. T. rangeli, like T. cruzi is endemic in many Latin American countries, whereas T. conorhini is tropicopolitan. T. rangeli and T. conorhini are exclusively extracellular, while T. cruzi has an intracellular stage in the mammalian host. Results Here we provide the first comprehensive sequence analysis of T. rangeli AM80 and T. conorhini 025E, and provide a comparison of their genomes to those of T. cruzi G and T. cruzi CL, respectively members of T. cruzi lineages TcI and TcVI. We report de novo assembled genome sequences of the low-virulent T. cruzi G, T. rangeli AM80, and T. conorhini 025E ranging from ~ 21–25 Mbp, with ~ 10,000 to 13,000 genes, and for the highly virulent and hybrid T. cruzi CL we present a ~ 65 Mbp in-house assembled haplotyped genome with ~ 12,500 genes per haplotype. Single copy orthologs of the two T. cruzi strains exhibited ~ 97% amino acid identity, and ~ 78% identity to proteins of T. rangeli or T. conorhini. Proteins of the latter two organisms exhibited ~ 84% identity. T. cruzi CL exhibited the highest heterozygosity. T. rangeli and T. conorhini displayed greater metabolic capabilities for utilization of complex carbohydrates, and contained fewer retrotransposons and multigene family copies, i.e. trans-sialidases, mucins, DGF-1, and MASP, compared to T. cruzi. Conclusions Our analyses of the T. rangeli and T. conorhini genomes closely reflected their phylogenetic proximity to the T. cruzi clade, and were largely consistent with their divergent life cycles. Our results provide a greater context for understanding the life cycles, host range expansion, immunity evasion, and pathogenesis of these trypanosomatids. Electronic supplementary material The online version of this article (10.1186/s12864-018-5112-0) contains supplementary material, which is available to authorized users.

Annotated sequence data are instrumental in nearly all realms of biology. However, the advent of next-generation sequencing has rapidly facilitated an imbalance between accurate sequence data and accurate annotation data. To increase the annotation accuracy of the Caulobacter vibrioides CB13b1a (CB13) genome, we compared the PGAP and RAST annotations of the CB13 genome. A total of 64 unique genes were identified in the PGAP annotation that were either completely or partially absent in the RAST annotation, and a total of 16 genes were identified in the RAST annotation that were not included in the PGAP annotation. Moreover, PGAP identified 73 frameshifted genes and 22 genes with an internal stop. In contrast, RAST annotated the larger segment of these frameshifted genes without indicating a change in reading frame may have occurred. The RAST annotation did not include any genes with internal stop codons, since it chose start codons that were after the internal stop. To confirm the discrepancies between the two annotations and verify the accuracy of the CB13 genome sequence data, we re-sequenced and re-annotated the entire genome and obtained an identical sequence, except in a small number of homopolymer regions. A genome sequence comparison between the two versions allowed us to determine the correct number of bases in each homopolymer region, which eliminated frameshifts for 31 genes annotated as frameshifted genes and removed 24 pseudogenes from the PGAP annotation. Both annotation systems correctly identified genes that were missed by the other system. In addition, PGAP identified conserved gene fragments that represented the beginning of genes, but it employed no corrective method to adjust the reading frame of frameshifted genes or the start sites of genes harboring an internal stop codon. In doing so, the PGAP annotation identified a large number of pseudogenes, which may reflect evolutionary history but likely do not produce gene products. These results demonstrate that re-sequencing and annotation comparisons can be used to increase the accuracy of genomic data and the corresponding gene annotation.

Background Pseudogenes are RNA transcripts with high homology with its parent protein-coding genes. Although pseudogenes lost the ability to produce protein, it still exert import biological function, and play important role in the pathogenesis of a wide varity of tumors; However, the role of pseudogenes in esophageal squamous cell carcinoma (ESCC) is poorly understood. Methods TUSC2P function in ESCC were explored using both in vitro and in vivo experiments cell proliferation, invasion and apoptosis assay was performed to evaluated the effect of TUSC2P on the tumor biology of ESCC. Expression of relative genes was assessed by quantitative real-time PCR (qRT-PCR) and western blotting in EC109 and TE-1 cell, as well as ESCC patients. 3’UTR luciferase assay was used to confirm the direct binding of miRNAs with TUSC2 and TUSC2P 3’UTR. Relation betweenTUSC2P, TUSC2 and ESCC prognosis was predicted by survival analysis (n = 56). Results Pseudogene TUSC2P was down regulated in ESCC tissues compared with paired normal adjacent tissues, and the expression of TUSC2P was significantly correlated with survivalof ESCC patients. Over expression of TUSC2P in EC109 and TE-1 cells resulted in altered expression of TUSC2, thus inhibited proliferation, invasion and promoted apoptosis. Dual luciferase assay demonstrated that TUSC2P 3’UTR decoyed miR-17-5p, miR-520a-3p, miR-608, miR-661 from binding to TUSC2. Conclusions TUSC2P can suppresses the tumor function of esophageal squamous cell carcinoma by regulating TUSC2 expression and may also serve as a prognostic factor for ESCC patients.

The competitive endogenous RNA (ceRNA) hypothesis is an attractively simple model to explain the biological role of many putatively functionless noncoding RNAs. Under this model, there exist transcripts in the cell whose role is to titrate out microRNAs such that the expression level of another target sequence is altered. That it is logistically possible for expression of one microRNA recognition element (MRE)-containing transcript to affect another is seen in the multiple examples of pathogenic effects of inappropriate expression of MRE-containing RNAs. However, the role, if any, of ceRNAs in normal biological processes and at physiological levels is disputed. By comparison of parent genes and pseudogenes we show, both for a specific example and genome-wide, that the pseudo-3′ untranslated regions (3′UTRs) of expressed pseudogenes are frequently retained and are under selective constraint in mammalian genomes. We found that the pseudo-3′UTR of BRAFP1, a previously described oncogenic ceRNA, has reduced substitutions relative to its pseudo-coding sequence, and we show sequence constraint on MREs shared between the parent gene, BRAF, and the pseudogene. Investigation of RNA-seq data reveals expression of BRAFP1 in normal somatic tissues in human and in other primates, consistent with biological ceRNA functionality of this pseudogene in nonpathogenic cellular contexts. Furthermore, we find that on a genome-wide scale pseudo-3′UTRs of mammalian pseudogenes (n = 1,629) are under stronger selective constraint than their pseudo-coding sequence counterparts, and are more often retained and expressed. Our results suggest that many human pseudogenes, often considered nonfunctional, may have an evolutionarily constrained role, consistent with the ceRNA hypothesis.

International audience; Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using LRRPROFILER, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Evolution sometimes proceeds by loss, especially when structures and genes become dispensable after an environmental shift relaxes functional constraints. Subterranean vertebrates are outstanding models to analyze this process, and gene decay can serve as a readout. We sought to understand some general principles on the extent and tempo of the decay of genes involved in vision, circadian clock, and pigmentation in cavefishes. The analysis of the genomes of two Cuban species belonging to the genus Lucifuga provided evidence for the largest loss of eye-specific genes and nonvisual opsin genes reported so far in cavefishes. Comparisons with a recently evolved cave population of Astyanax mexicanus and three species belonging to the Chinese tetraploid genus Sinocyclocheilus revealed the combined effects of the level of eye regression, time, and genome ploidy on eye-specific gene pseudogenization. The limited extent of gene decay in all these cavefishes and the very small number of loss-of-function mutations per pseudogene suggest that their eye degeneration may not be very ancient, ranging from early to late Pleistocene. This is in sharp contrast with the identification of several vision genes carrying many loss-of-function mutations in ancient fossorial mammals, further suggesting that blind fishes cannot thrive more than a few million years in cave ecosystems.

Gene rearrangements have been found in several mitochondrial genomes of Mantodea, located in the gene blocks CR-I-Q-M-ND2, COX1-K-D-ATP8 and ND3-A-R-N-S-E-F-ND5. We have sequenced one mitogenome of Amelidae (Yersinia mexicana) and six mitogenomes of Mantidae to discuss the mitochondrial gene rearrangement and the phylogenetic relationship within Mantidae. These mitogenomes showed rearrangements of tRNA genes except for Asiadodis yunnanensis and Hierodula zhangi. These novel gene rearrangements of Mantidae were primarily concentrated in the region of CR-I-Q-M-ND2, including gene translocation, duplication and pseudogenization. For the occurrences of these rearrangements, the tandem duplication-random loss (TDRL) model and slipped-strand mispairing model were suitable to explain. Large non-coding regions (LNCRs) located in the region of CR-I-Q-M-ND2 were detected in most Mantidae species, whereas some LNCRs had high similarity to the control region (CR). Both BI and ML phylogenetic analyses supported the monophyly of Mantidae and the paraphyly of Mantinae. The phylogenetic results with the gene order and the location of NCRs acted as forceful evidence that specific gene rearrangements and special LNCRs may be synapomorphies for several groups of mantises.

International audience; The discovery that pseudogenes are involved in important biological processes has excited enthusiasm and increased the research interest on them. An accurate detection and analysis of pseudogenes can be achieved using comparative methods, but only the use of phylogenetic tools can provide accurate information about their birth, their evolution and their death, hence about the impact that they have on genes and genomes. Here, phylogenetic methods that allow for studying pseudogene history are described.

International audience; Most genomic and evolutionary comparative analyses rely on accurate multiple sequence alignments. With their underlying codon structure, protein-coding nucleotide sequences pose a specific challenge for multiple sequence alignment. Multiple Alignment of Coding Sequences (MACSE) is a multiple sequence alignment program that provided the first automatic solution for aligning protein-coding gene datasets containing both functional and nonfunctional sequences (pseudogenes). Through its unique features, reliable codon alignments can be built in the presence of frameshifts and stop codons suitable for subsequent analysis of selection based on the ratio of nonsynonymous to synonymous substitutions. Here we offer a practical overview and guidelines on the use of MACSE v2. This major update of the initial algorithm now comes with a graphical interface providing user-friendly access to different subprograms to handle multiple alignments of protein-coding sequences. We also present new pipelines based on MACSE v2 subprograms to handle large datasets and distributed as Singularity containers. MACSE and associated pipelines are available at: https://bioweb.supagro.inra.fr/macse/.

Frequent deletions impacting the FRA4F fragile site leads to the oncogenic overexpression of the corresponding CCSER1 gene, through the deletion of the TMSB4XP8 pseudogene. The ensuing CCSER1 overexpression elicits mitotic instability, targetable with aurora kinase inhibitors., The oncogenic role of common fragile sites (CFS), focal and pervasive gaps in the cancer genome arising from replicative stress, remains controversial. Exploiting the TCGA dataset, we found that in most CFS the genes residing within the associated focal deletions are down-regulated, including proteins involved in tumour immune recognition. In a subset of CFS, however, the residing genes are surprisingly overexpressed. Within the most frequent CFS in this group, FRA4F, which is deleted in up to 18% of cancer cases and harbours the CCSER1 gene, we identified a region which includes an intronic, antisense pseudogene, TMSB4XP8. TMSB4XP8 focal ablation or transcriptional silencing elicits the overexpression of CCSER1, through a cis-acting mechanism. CCSER1 overexpression increases proliferation and triggers centrosome amplifications, multinuclearity, and aberrant mitoses. Accordingly, FRA4F is associated in patient samples to mitotic genes deregulation and genomic instability. As a result, cells overexpressing CCSER1 become sensitive to the treatment with aurora kinase inhibitors. Our findings point to a novel tumourigenic mechanism where focal deletions increase the expression of a new class of “dormant” oncogenes.

Pseudogenes are frequently encountered noncoding sequences with a high sequence similarity to their protein-coding paralogue. For this reason, their presence is often considered troublesome in molecular diagnostics. In pseudoxanthoma elasticum (PXE), a disease predominantly caused by mutations in ATP-binding cassette family C member 6 (ABCC6), the presence of two pseudogenes complicates the analysis of sequence data. With whole-exome sequencing (WES) becoming the standard of care in molecular diagnostics, we wanted to evaluate whether this technique is as reliable as gene-specific targeted enrichment analysis for the analysis of ABCC6. We established a PCR-based targeted enrichment and next-generation sequencing testing approach and demonstrated that the ABCC6-specific enrichment combined with the applied mapping algorithm overcomes the complication of ABCC6 pseudogene aspecificities, contrary to WES. We propose a time- and cost-efficient diagnostic strategy for comprehensive and accurate molecular genetic testing of PXE, which is highly automatable.

This is the first study in which the Daphnia magna (D. magna) nuclear genome (nDNA) obtained from the GenBank database was analyzed for pseudogene sequences of mitochondrial origin. To date, there is no information about pseudogenes localized in D. magna genome. This study aimed to identify NUMTs, their length, homology, and location for potential use in evolutionary studies and to check whether their occurrence causes co-amplification during mitochondrial genome (mtDNA) analyses. Bioinformatic analysis showed 1909 fragments of the mtDNA of D. magna, of which 1630 were located in ten linkage groups (LG) of the nDNA. The best-matched NUMTs covering &gt, 90% of the gene sequence have been identified for two mt-tRNA genes, and they may be functional nuclear RNA molecules. Isolating the total DNA in mtDNA studies, co-amplification of nDNA fragments is unlikely in the case of amplification of the whole tRNA genes as well as fragments of other genes. It was observed that TRNA-MET fragments had the highest level of sequence homology, thus they could be evolutionarily the youngest. The lowest homology was found in the D-loop-derived pseudogene. It may probably be the oldest NUMT incorporated into the nDNA, however, further analysis is necessary.