Identification and characterization of over 100 mitochondrial ribosomal protein pseudogenes in the human genome

Pseudogenes and their potential functions in hematopoiesis

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

Arginine and lysine methylation of MRPS23 promotes breast cancer metastasis through regulating OXPHOS

Mitochondrial oxidative phosphorylation (OXPHOS) is a vital regulator of tumor metastasis. However, the mechanisms governing OXPHOS to facilitate tumor metastasis remain unclear. In this study, we discovered that arginine 21(R21) and lysine 108 (K108) of mitochondrial ribosomal protein S23 (MRPS23) was methylated by the protein arginine methyltransferase 7 (PRMT7) and SET-domain-containing protein 6 (SETD6), respectively. R21 methylation accelerated the poly-ubiquitin-dependent degradation of MRPS23 to a low level. The MRPS23 degradation inhibited OXPHOS with elevated mtROS level, which consequently increased breast cancer cell invasion and metastasis. In contrast, K108 methylation increased MRPS23 stability, and K108 methylation coordinated with R21 methylation to maintain a low level of MRPS23, which was in favor of supporting breast cancer cell survival through regulating OXPHOS. Consistently, R21 and K108 methylation was correlated with malignant breast carcinoma. Significantly, our findings unveil a unique mechanism of controlling OXPHOS by arginine and lysine methylation and point to the impact of the PRMT7-SETD6-MRPS23 axis during breast cancer metastasis.

Signal Transduction in Ribosome Biogenesis: A Recipe to Avoid Disaster

Energetically speaking, ribosome biogenesis is by far the most costly process of the cell and, therefore, must be highly regulated in order to avoid unnecessary energy expenditure. Not only must ribosomal RNA (rRNA) synthesis, ribosomal protein (RP) transcription, translation, and nuclear import, as well as ribosome assembly, be tightly controlled, these events must be coordinated with other cellular events, such as cell division and differentiation. In addition, ribosome biogenesis must respond rapidly to environmental cues mediated by internal and cell surface receptors, or stress (oxidative stress, DNA damage, amino acid depletion, etc.). This review examines some of the well-studied pathways known to control ribosome biogenesis (PI3K-AKT-mTOR, RB-p53, MYC) and how they may interact with some of the less well studied pathways (eIF2α kinase and RNA editing/splicing) in higher eukaryotes to regulate ribosome biogenesis, assembly, and protein translation in a dynamic manner.

Multipronged approach to identify and validate a novel upstream regulator of Sncg in mouse retinal ganglion cells

Loss of retinal ganglion cells (RGCs) is one of the hallmarks of retinal neurodegenerative diseases, glaucoma being one of the most common. Mechanistic studies on RGCs are hindered by the lack of sufficient primary cells and consensus regarding their signature markers. Recently, γ-synuclein (SNCG) has been shown to be highly expressed in the somas and axons of RGCs. In various mouse models of glaucoma, downregulation of Sncg gene expression correlates with RGC loss. To investigate the role of Sncg in RGCs, we used a novel systems genetics approach to identify a gene that modulates Sncg expression, followed by confirmatory studies in both healthy and diseased retinae. We found that chromosome 1 harbors an expression quantitative trait locus that modulates Sncg expression in the mouse retina, and identified the prefoldin-2 (PFDN2) gene as the candidate upstream modulator of Sncg expression. Our immunohistochemical analyses revealed similar expression patterns in both mouse and human healthy retinae, with PFDN2 colocalizing with SNCG in RGCs and their axons. In contrast, in retinae from glaucoma subjects, SNCG levels were significantly reduced, although PFDN2 levels were maintained. Using a novel flow cytometry-based RGC isolation method, we obtained viable populations of murine RGCs. Knocking down Pfdn2 expression in primary murine RGCs significantly reduced Sncg expression, confirming that Pfdn2 regulates Sncg expression in murine RGCs. Gene Ontology analysis indicated shared mitochondrial function associated with Sncg and Pfdn2. These data solidify the relationship between Sncg and Pfdn2 in RGCs, and provide a novel mechanism for maintaining RGC health.

Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction?

Mitochondrial dysfunction is central to numerous diseases of oxidative stress. Changes in mitochondrial DNA (MtDNA) content, often measured as mitochondrial genome to nuclear genome ratio (Mt/N) using real time quantitative PCR, have been reported in a broad range of human diseases, such as diabetes and its complications, obesity, cancer, HIV complications, and ageing. We propose the hypothesis that MtDNA content in body fluids and tissues could be a biomarker of mitochondrial dysfunction and review the evidence supporting this theory. Increased reactive oxygen species resulting from an external trigger such as hyperglycaemia or increased fat in conditions of oxidative stress could lead to enhanced mitochondrial biogenesis, and increased Mt/N. Altered MtDNA levels may contribute to enhanced oxidative stress and inflammation and could play a pathogenic role in mitochondrial dysfunction and disease. Changes in Mt/N are detectable in circulating cells such as peripheral blood mononuclear cells and these could be used as surrogate to predict global changes in tissues and organs. We review a large number of studies reporting changes in MtDNA levels in body fluids such as circulating blood cells, cell free serum, saliva, sperm, and cerebrospinal fluid as well as in tumour and normal tissue samples. However, the data are often conflicting as the current methodology used to measure Mt/N can give false results because of one or more of the following reasons (1) use of mitochondrial primers which co-amplify nuclear pseudogenes (2) use of nuclear genes which are variable and/or duplicated in numerous locations (3) a dilution bias caused by the differing genome sizes of the mitochondrial and nuclear genome and (4) template preparation protocols which affect the yields of nuclear and mitochondrial genomes. Development of robust and reproducible methodology is needed to test the hypothesis that MtDNA content in body fluids is biomarker of mitochondrial dysfunction.

Proteome profile of functional mitochondria from human skeletal muscle using one-dimensional gel electrophoresis and HPLC-ESI-MS/MS

Mitochondria can be isolated from skeletal muscle in a manner that preserves tightly coupled bioenergetic function in vitro. The purpose of this study was to characterize the composition of such preparations using a proteomics approach. Mitochondria isolated from human vastus lateralis biopsies were functional as evidenced by their response to carbohydrate and fat-derived fuels. Using one-dimensional gel electrophoresis and HPLC-ESI-MS/MS, 823 unique proteins were detected, and 487 of these were assigned to the mitochondrion, including the newly characterized SIRT5, MitoNEET and RDH13. Proteins detected included 9 of the 13 mitochondrial DNA-encoded proteins and 86 of 104 electron transport chain (ETC) and ETC-related proteins. In addition, 59 of 78 proteins of the 55S mitoribosome, several TIM and TOM proteins and cell death proteins were present. This study presents an efficient method for future qualitative assessments of proteins from functional isolated mitochondria from small samples of healthy and diseased skeletal muscle.

Organization and evolution of mitochondrial gene clusters in human

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average+/-SD 1.3+/-0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.

The human ABC transporter pseudogene family: Evidence for transcription and gene-pseudogene interference

Background

Pseudogenes are an integral component of the human genome. Little attention, however, has so far been paid to the phenomenon that some pseudogenes are transcriptionally active. Recently, we demonstrated that the human ortholog of the rodent testis-specific ATP-binding cassette (ABC) transporter Abca17 is a ubiquitously transcribed pseudogene (ABCA17P). The aim of the present study was to establish a complete inventory of all ABC transporter pseudogenes in the human genome and to identify transcriptionally active ABC transporter pseudogenes. Moreover, we tested the hypothesis that a regulatory interdependency exists between ABC transporter pseudogenes and their parental protein coding equivalents.

Results

Systematic bioinformatic analysis revealed the existence of 22 ABC transporter pseudogenes within the human genome. We identified two clusters on chromosomes 15 and 16, respectively, which harbor almost half of all pseudogenes (n = 10). Available information from EST and mRNA databases and RT-PCR expression profiling indicate that a large portion of the ABC transporter pseudogenes (45%, n = 10) are transcriptionally active and some of them are expressed as alternative splice variants. We demonstrate that both pseudogenes of the pseudoxanthoma elasticum gene ABCC6, ABCC6P1 and ABCC6P2, are transcribed. ABCC6P1 and ABCC6 possess near-identical promoter sequences and their tissue-specific expression profiles are strikingly similar raising the possibility that they form a gene-pseudogene dual transcription unit. Intriguingly, targeted knockdown of the transcribed pseudogene ABCC6P1 resulted in a significant reduction of ABCC6 mRNA expression levels.

Conclusion

The human genome contains a surprisingly small number of ABC transporter pseudogenes relative to other known gene families. They are unevenly distributed across the chromosomes. Importantly, a significant portion of the ABC transporter pseudogenes is transcriptionally active. The downregulation of ABCC6 mRNA levels by targeted suppression of the expression of its pseudogene ABCC6P1 provides evidence, for the first time, for a regulatory interdependence of a transcribed pseudogene and its protein coding counterpart in the human genome.

MRPS18CP2 alleles and DEFA3 absence as putative chromosome 8p23.1 modifiers of hearing loss due to mtDNA mutation A1555G in the 12S rRNA gene

Background

Mitochondrial DNA (mtDNA) mutations account for at least 5% of cases of postlingual, nonsyndromic hearing impairment. Among them, mutation A1555G is frequently found associated with aminoglycoside-induced and/or nonsyndromic hearing loss in families presenting with extremely variable clinical phenotypes. Biochemical and genetic data have suggested that nuclear background is the main factor involved in modulating the phenotypic expression of mutation A1555G. However, although a major nuclear modifying locus was located on chromosome 8p23.1 and regardless intensive screening of the region, the gene involved has not been identified.

Methods

With the aim to gain insights into the factors that determine the phenotypic expression of A1555G mutation, we have analysed in detail different genetic and genomic elements on 8p23.1 region (DEFA3 gene absence, CLDN23 gene and MRPS18CP2 pseudogene) in a group of 213 A1555G carriers.

Results

Family based association studies identified a positive association for a polymorphism on MRPS18CP2 and an overrepresentation of DEFA3 gene absence in the deaf group of A1555G carriers.

Conclusion

Although none of the factors analysed seem to have a major contribution to the phenotype, our findings provide further evidences of the involvement of 8p23.1 region as a modifying locus for A1555G 12S rRNA gene mutation.

Pseudogene.org: a comprehensive database and comparison platform for pseudogene annotation

The Pseudogene.org knowledgebase serves as a comprehensive repository for pseudogene annotation. The definition of a pseudogene varies within the literature, resulting in significantly different approaches to the problem of identification. Consequently, it is difficult to maintain a consistent collection of pseudogenes in detail necessary for their effective use. Our database is designed to address this issue. It integrates a variety of heterogeneous resources and supports a subset structure that highlights specific groups of pseudogenes that are of interest to the research community. Tools are provided for the comparison of sets and the creation of layered set unions, enabling researchers to derive a current ‘consensus’ set of pseudogenes. Additional features include versatile search, the capacity for robust interaction with other databases, the ability to reconstruct older versions of the database (accounting for changing genome builds) and an underlying object-oriented interface designed for researchers with a minimal knowledge of programming. At the present time, the database contains more than 100 000 pseudogenes spanning 64 prokaryote and 11 eukaryote genomes, including a collection of human annotations compiled from 16 sources.

Retroposition of processed pseudogenes: the impact of RNA stability and translational control

Human processed pseudogenes are copies of cellular RNAs reverse transcribed and inserted into the nuclear genome by the enzymatic machinery of L1 (LINE1) non-LTR retrotransposons. Although it is generally accepted that germline expression is crucial for the heritable retroposition of cellular mRNAs, little is known about the influences of RNA stability, mRNA quality control and compartmentalization of translation on the retroposition of processed pseudogenes. We found that frequently retroposed human mRNAs are derived from stable transcripts with translation-competent functional reading frames that are resistant to nonsense-mediated RNA decay. They are preferentially translated on free cytoplasmic ribosomes and encode soluble proteins. Our results indicate that interactions between mRNAs and L1 proteins seem to occur at free cytoplasmic ribosomes.

Multiple independent pseudogene derivations indicate increased instability of the Mdm2 locus in Mus caroli

Under conditions of genomic stress, the Mdm locus in human and in mouse is prone to instability manifested as amplification and oncogenesis. The Mdm2 gene is a known oncogene that is amplified in approximately one-third of sarcomas and whose protein product interacts with the tumor suppressor p53. Concimitant with such gene amplification events is the activation and mobilization of endogenous retroelements, typically through the relaxation of epigenetic controlling mechanisms. Processed pseudogenes, which can be formed through endogenous LINE retroelement activity, may indicate increased genomic instability. We have isolated processed pseudogenes for Mdm2 in Mus caroli DNA, likely formed from independent events in different individuals. This is the first identification and characterization of an Mdm2 pseudogene in any organism. Multiple retrotransposition events are suggested by the variable sequence and genomic structure of the identified pseudogenes across all exons and the 3'UTR. The high degree of similarity between the gene and each pseudogene, as well as the lack of evidence for an Mdm2 pseudogene in several other species of Mus, indicate evolutionarily recent retrotransposition events leading to the formation of the Mdm2 pseudogenes in M. caroli. Previous studies on the Mdm2 locus in Mus caroli showed amplification and overexpression of this gene on double minute chromosomes in a Mus musculus x Mus caroli interspecific hybrid. The identification of an Mdm2 retropseudogene within this species further highlights the predisposition to instability for this region of the genome.

Whole genome association study of rheumatoid arthritis using 27 039 microsatellites

A major goal of current human genome-wide studies is to identify the genetic basis of complex disorders. However, the availability of an unbiased, reliable, cost efficient and comprehensive methodology to analyze the entire genome for complex disease association is still largely lacking or problematic. Therefore, we have developed a practical and efficient strategy for whole genome association studies of complex diseases by charting the human genome at 100 kb intervals using a collection of 27,039 microsatellites and the DNA pooling method in three successive genomic screens of independent case-control populations. The final step in our methodology consists of fine mapping of the candidate susceptible DNA regions by single nucleotide polymorphisms (SNPs) analysis. This approach was validated upon application to rheumatoid arthritis, a destructive joint disease affecting up to 1% of the population. A total of 47 candidate regions were identified. The top seven loci, withstanding the most stringent statistical tests, were dissected down to individual genes and/or SNPs on four chromosomes, including the previously known 6p21.3-encoded Major Histocompatibility Complex gene, HLA-DRB1. Hence, microsatellite-based genome-wide association analysis complemented by end stage SNP typing provides a new tool for genetic dissection of multifactorial pathologies including common diseases.

Identification of 11 pseudogenes in the DNA methyltransferase gene family in rodents and humans and implications for the functional loci

DNA (cytosine-5-)-methyltransferase genes are important for normal development in mice and humans. We describe here 11 pseudogenes spread among human, mouse, and rat belonging to this gene family, ranging from 1 pseudogene in humans to 7 in rat, all belonging to the Dnmt3 subfamily. All except 1 rat Dnmt3b pseudogene appear to be transcriptionally silent. Dnmt3a2, a transcript variant of Dnmt3a starting at an alternative promoter, had the highest number of processed pseudogenes, while none were found for the canonical Dnmt3a, suggesting the former transcript is more highly expressed in germ cells. Comparison of human, mouse, and rat Dnmt3a2 sequences also suggests that human exon 8 is a recent acquisition. Alignment of the 3'UTR of Dnmt3a2 among the functional genes and the processed pseudogenes suggested that a second polyadenylation site downstream of the RefSeq poly(A) was being used in mice, resulting in a longer 3'UTR, a finding confirmed by RT-PCR in mouse tissues. We also found conserved cytoplasmic polyadenylation elements, usually implicated in regulating translation in oocytes, in Dnmt3b and Dnmt1. Expression of DNMT3B in the mouse oocyte was confirmed by immunocytochemistry. These results clarify the structure of a number of loci in the three species examined and provide some useful insights into the structure and evolution of this gene family.

Millions of years of evolution preserved: a comprehensive catalog of the processed pseudogenes in the human genome

Processed pseudogenes were created by reverse-transcription of mRNAs; they provide snapshots of ancient genes existing millions of years ago in the genome. To find them in the present-day human, we developed a pipeline using features such as intron-absence, frame-disruption, polyadenylation, and truncation. This has enabled us to identify in recent genome drafts approximately 8000 processed pseudogenes (distributed from http://pseudogene.org). Overall, processed pseudogenes are very similar to their closest corresponding human gene, being 94% complete in coding regions, with sequence similarity of 75% for amino acids and 86% for nucleotides. Their chromosomal distribution appears random and dispersed, with the numbers on chromosomes proportional to length, suggesting sustained "bombardment" over evolution. However, it does vary with GC-content: Processed pseudogenes occur mostly in intermediate GC-content regions. This is similar to Alus but contrasts with functional genes and L1-repeats. Pseudogenes, moreover, have age profiles similar to Alus. The number of pseudogenes associated with a given gene follows a power-law relationship, with a few genes giving rise to many pseudogenes and most giving rise to few. The prevalence of processed pseudogenes agrees well with germ-line gene expression. Highly expressed ribosomal proteins account for approximately 20% of the total. Other notables include cyclophilin-A, keratin, GAPDH, and cytochrome c.