The human genome has 49 cytochrome c pseudogenes, including a relic of a primordial gene that still functions in mouse

Prostate Cancer-Specific Lysine 53 Acetylation of Cytochrome c Drives Metabolic Reprogramming and Protects from Apoptosis in Intact Cells

Cytochrome c (Cytc) is important for both mitochondrial respiration and apoptosis, both of which are altered in cancer cells that switch to Warburg metabolism and manage to evade apoptosis. We earlier reported that lysine 53 (K53) of Cytc is acetylated in prostate cancer. K53 is conserved in mammals that is known to be essential for binding to cytochrome c oxidase and apoptosis protease activating factor-1 (Apaf-1). Here we report the effects of this acetylation on the main functions of cytochrome c by expressing acetylmimetic K53Q in cytochrome c double knockout cells. Other cytochrome c variants analyzed were wild-type, K53R as a control that maintains the positive charge, and K53I, which is present in some non-mammalian species. Intact cells expressing K53Q cytochrome c showed 49% decreased mitochondrial respiration and a concomitant increase in glycolytic activity (Warburg effect). Furthermore, mitochondrial membrane potential was decreased, correlating with notably reduced basal mitochondrial superoxide levels and decreased cell death upon challenge with H2O2 or staurosporine. To test for markers of cancer aggressiveness and invasiveness, cells were grown in 3D spheroid culture. K53Q cytochrome c-expressing cells showed profoundly increased protrusions compared to WT, suggesting increased invasiveness. We propose that K53 acetylation of cytochrome c is an adaptive response that mediates prostate cancer metabolic reprogramming and evasion of apoptosis, which are two hallmarks of cancer, to better promote tumor survival and metastasis.

Direct tests of cytochrome c and c1 functions in the electron transport chain of malaria parasites

The mitochondrial electron transport chain (ETC) of Plasmodium malaria parasites is a major antimalarial drug target, but critical cytochrome (cyt) functions remain unstudied and enigmatic. Parasites express two distinct cyt c homologs (c and c-2) with unusually sparse sequence identity and uncertain fitness contributions. P. falciparum cyt c-2 is the most divergent eukaryotic cyt c homolog currently known and has sequence features predicted to be incompatible with canonical ETC function. We tagged both cyt c homologs and the related cyt c1 for inducible knockdown. Translational repression of cyt c and cyt c1 was lethal to parasites, which died from ETC dysfunction and impaired ubiquinone recycling. In contrast, cyt c-2 knockdown or knockout had little impact on blood-stage growth, indicating that parasites rely fully on the more conserved cyt c for ETC function. Biochemical and structural studies revealed that both cyt c and c-2 are hemylated by holocytochrome c synthase, but UV-vis absorbance and EPR spectra strongly suggest that cyt c-2 has an unusually open active site in which heme is stably coordinated by only a single axial amino acid ligand and can bind exogenous small molecules. These studies provide a direct dissection of cytochrome functions in the ETC of malaria parasites and identify a highly divergent Plasmodium cytochrome c with molecular adaptations that defy a conserved role in eukaryotic evolution.

Direct Tests of Cytochrome Function in the Electron Transport Chain of Malaria Parasites

The mitochondrial electron transport chain (ETC) of Plasmodium malaria parasites is a major antimalarial drug target, but critical cytochrome functions remain unstudied and enigmatic. Parasites express two distinct cyt c homologs ( c and c -2) with unusually sparse sequence identity and uncertain fitness contributions. P. falciparum cyt c -2 is the most divergent eukaryotic cyt c homolog currently known and has sequence features predicted to be incompatible with canonical ETC function. We tagged both cyt c homologs and the related cyt c ₁ for inducible knockdown. Translational repression of cyt c and cyt c ₁ was lethal to parasites, which died from ETC dysfunction and impaired ubiquinone recycling. In contrast, cyt c -2 knockdown or knock-out had little impact on blood-stage growth, indicating that parasites rely fully on the more conserved cyt c for ETC function. Biochemical and structural studies revealed that both cyt c and c -2 are hemylated by holocytochrome c synthase, but UV-vis absorbance and EPR spectra strongly suggest that cyt c -2 has an unusually open active site in which heme is stably coordinated by only a single axial amino-acid ligand and can bind exogenous small molecules. These studies provide a direct dissection of cytochrome functions in the ETC of malaria parasites and identify a highly divergent Plasmodium cytochrome c with molecular adaptations that defy a conserved role in eukaryotic evolution.

SIGNIFICANCE STATEMENT

Mitochondria are critical organelles in eukaryotic cells that drive oxidative metabolism. The mitochondrion of Plasmodium malaria parasites is a major drug target that has many differences from human cells and remains poorly studied. One key difference from humans is that malaria parasites express two cytochrome c proteins that differ significantly from each other and play untested and uncertain roles in the mitochondrial electron transport chain (ETC). Our study revealed that one cyt c is essential for ETC function and parasite viability while the second, more divergent protein has unusual structural and biochemical properties and is not required for growth of blood-stage parasites. This work elucidates key biochemical properties and evolutionary differences in the mitochondrial ETC of malaria parasites.

Tissue-specific regulation of cytochrome c by post-translational modifications: respiration, the mitochondrial membrane potential, ROS, and apoptosis

Cytochrome c (Cyt c) plays a vital role in the mitochondrial electron transport chain (ETC). In addition, it is a key regulator of apoptosis. Cyt c has multiple other functions including ROS production and scavenging, cardiolipin peroxidation, and mitochondrial protein import. Cyt c is tightly regulated by allosteric mechanisms, tissue-specific isoforms, and post-translational modifications (PTMs). Distinct residues of Cyt c are modified by PTMs, primarily phosphorylations, in a highly tissue-specific manner. These modifications downregulate mitochondrial ETC flux and adjust the mitochondrial membrane potential (ΔΨm), to minimize reactive oxygen species (ROS) production under normal conditions. In pathologic and acute stress conditions, such as ischemia-reperfusion, phosphorylations are lost, leading to maximum ETC flux, ΔΨm hyperpolarization, excessive ROS generation, and the release of Cyt c. It is also the dephosphorylated form of the protein that leads to maximum caspase activation. We discuss the complex regulation of Cyt c and propose that it is a central regulatory step of the mammalian ETC that can be rate limiting in normal conditions. This regulation is important because it maintains optimal intermediate ΔΨm, limiting ROS generation. We examine the role of Cyt c PTMs, including phosphorylation, acetylation, methylation, nitration, nitrosylation, and sulfoxidation and consider their potential biological significance by evaluating their stoichiometry.-Kalpage, H. A., Bazylianska, V., Recanati, M. A., Fite, A., Liu, J., Wan, J., Mantena, N., Malek, M. H., Podgorski, I., Heath, E. I., Vaishnav, A., Edwards, B. F., Grossman, L. I., Sanderson, T. H., Lee, I., Hüttemann, M. Tissue-specific regulation of cytochrome c by post-translational modifications: respiration, the mitochondrial membrane potential, ROS, and apoptosis.

Pseudogenes in gastric cancer pathogenesis: a review article

Cancer burden rises globally at an alarming pace. According to GLOBOCAN 2012, gastric cancer (GC) is regarded as the fifth most common malignancy in the world. Being twice as high in men as in women, GC is the third leading cause of cancer mortality in both sexes globally. Being labeled as 'junk DNA', pseudogenes were considered as nonfunctional 'trash', which contribute nothing to survival of the organism; therefore, a number of strategies have been developed to circumvent their accidental detection. Recent progresses have confirmed that pseudogenes can have broad and multifaceted spectrum of activities in human cancers in general and GC in particular. Furthermore, the mentioned functions are parental gene-dependent and/or -independent. Therefore, pseudogenes can be regarded as the emerging class of elaborate modulators of gene expression involved in pathogenesis of human cancers including gastric adenocarcinoma.

Pseudogenes and Their Genome-Wide Prediction in Plants

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

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

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

Utility of syntenic relationships of VDAC1 pseudogenes for not only an understanding of the phylogenetic divergence history of rodents, but also ascertaining possible pseudogene candidates as genuine pseudogenes

Rodent and human genomes were screened to identify pseudogenes of the type 1 voltage-dependent anion channel (VDAC1) in mitochondria. In addition to the 16 pseudogenes of rat VDAC1 identified in our recent study, 15 and 13 sequences were identified as pseudogenes of VDAC1 in mouse and human genome, respectively; and 4, 2, and 1 sequences, showing lower similarities with the VDAC1 sequence, were identified as "possible pseudogene candidates" in rat, mouse, and human, respectively. No syntenic combination was observed between rodent and human pseudogenes, but 2 and 1 possible pseudogene candidates of VDAC1 of rat and mouse, respectively, were found to have syntenic counterparts in mouse and rat genome, respectively; and these syntenic counterparts were genuine VDAC1 pseudogenes. Therefore, syntenic combinations of pseudogenes of VDAC1 were useful not only for a better understanding of the phylogenetic divergence history of rodents but also for ascertaining possible pseudogene candidates as genuine pseudogenes.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

Pseudogenes and their composers: delving in the 'debris' of human genome

Pseudogenes, the nonfunctional homologs of functional genes and thus exemplified as 'genomic fossils' provide intriguing snapshots of the evolutionary history of human genome. These defunct copies generally arise by retrotransposition or duplication followed by various genetic disablements. In this study, focusing on human pseudogenes and their functional homologues we describe their characteristic features and relevance to protein sequence evolution. We recapitulate that pseudogenes harbor disease-causing degenerative sequence variations in conjunction with the immense disease gene association of their progenitors. Furthermore, we also discuss the issue of functional resurrection and the potentiality observed in some pseudogenes to regulate their functional counterparts.

Pseudogenes: newly discovered players in human cancer

Because they are generally noncoding and thus considered nonfunctional and unimportant, pseudogenes have long been neglected. Recent advances have established that the DNA of a pseudogene, the RNA transcribed from a pseudogene, or the protein translated from a pseudogene can have multiple, diverse functions and that these functions can affect not only their parental genes but also unrelated genes. Therefore, pseudogenes have emerged as a previously unappreciated class of sophisticated modulators of gene expression, with a multifaceted involvement in the pathogenesis of human cancer.

Phosphorylation of mammalian cytochrome c and cytochrome c oxidase in the regulation of cell destiny: respiration, apoptosis, and human disease

The mitochondrial oxidative phosphorylation (OxPhos) system not only generates the vast majority of cellular energy, but is also involved in the generation of reactive oxygen species (ROS), and apoptosis. Cytochrome c (Cytc) and cytochrome c oxidase (COX) represent the terminal step of the electron transport chain (ETC), the proposed rate-limiting reaction in mammals. Cytc and COX show unique regulatory features including allosteric regulation, isoform expression, and regulation through cell signaling pathways. This chapter focuses on the latter and discusses all mapped phosphorylation sites based on the crystal structures of COX and Cytc. Several signaling pathways have been identified that target COX including protein kinase A and C, receptor tyrosine kinase, and inflammatory signaling. In addition, four phosphorylation sites have been mapped on Cytc with potentially large implications due to its multiple functions including apoptosis, a pathway that is overactive in stressed cells but inactive in cancer. The role of COX and Cytc phosphorylation is reviewed in a human disease context, including cancer, inflammation, sepsis, asthma, and ischemia/reperfusion injury as seen in myocardial infarction and ischemic stroke.

Regulation of mitochondrial respiration and apoptosis through cell signaling: cytochrome c oxidase and cytochrome c in ischemia/reperfusion injury and inflammation

Cytochrome c (Cytc) and cytochrome c oxidase (COX) catalyze the terminal reaction of the mitochondrial electron transport chain (ETC), the reduction of oxygen to water. This irreversible step is highly regulated, as indicated by the presence of tissue-specific and developmentally expressed isoforms, allosteric regulation, and reversible phosphorylations, which are found in both Cytc and COX. The crucial role of the ETC in health and disease is obvious since it, together with ATP synthase, provides the vast majority of cellular energy, which drives all cellular processes. However, under conditions of stress, the ETC generates reactive oxygen species (ROS), which cause cell damage and trigger death processes. We here discuss current knowledge of the regulation of Cytc and COX with a focus on cell signaling pathways, including cAMP/protein kinase A and tyrosine kinase signaling. Based on the crystal structures we highlight all identified phosphorylation sites on Cytc and COX, and we present a new phosphorylation site, Ser126 on COX subunit II. We conclude with a model that links cell signaling with the phosphorylation state of Cytc and COX. This in turn regulates their enzymatic activities, the mitochondrial membrane potential, and the production of ATP and ROS. Our model is discussed through two distinct human pathologies, acute inflammation as seen in sepsis, where phosphorylation leads to strong COX inhibition followed by energy depletion, and ischemia/reperfusion injury, where hyperactive ETC complexes generate pathologically high mitochondrial membrane potentials, leading to excessive ROS production. Although operating at opposite poles of the ETC activity spectrum, both conditions can lead to cell death through energy deprivation or ROS-triggered apoptosis.

Silencing, positive selection and parallel evolution: busy history of primate cytochromes C

Cytochrome c (cyt c) participates in two crucial cellular processes, energy production and apoptosis, and unsurprisingly is a highly conserved protein. However, previous studies have reported for the primate lineage (i) loss of the paralogous testis isoform, (ii) an acceleration and then a deceleration of the amino acid replacement rate of the cyt c somatic isoform, and (iii) atypical biochemical behavior of human cyt c. To gain insight into the cause of these major evolutionary events, we have retraced the history of cyt c loci among primates. For testis cyt c, all primate sequences examined carry the same nonsense mutation, which suggests that silencing occurred before the primates diversified. For somatic cyt c, maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses yielded the same tree topology. The evolutionary analyses show that a fast accumulation of non-synonymous mutations (suggesting positive selection) occurred specifically on the anthropoid lineage root and then continued in parallel on the early catarrhini and platyrrhini stems. Analysis of evolutionary changes using the 3D structure suggests they are focused on the respiratory chain rather than on apoptosis or other cyt c functions. In agreement with previous biochemical studies, our results suggest that silencing of the cyt c testis isoform could be linked with the decrease of primate reproduction rate. Finally, the evolution of cyt c in the two sister anthropoid groups leads us to propose that somatic cyt c evolution may be related both to COX evolution and to the convergent brain and body mass enlargement in these two anthropoid clades.

The multiple functions of cytochrome c and their regulation in life and death decisions of the mammalian cell: From respiration to apoptosis

Cytochrome c (Cytc) is essential in mitochondrial electron transport and intrinsic type II apoptosis. Mammalian Cytc also scavenges reactive oxygen species (ROS) under healthy conditions, produces ROS with the co-factor p66(Shc), and oxidizes cardiolipin during apoptosis. The recent finding that Cytc is phosphorylated in vivo underpins a model for the pivotal role of Cytc regulation in making life and death decisions. An apoptotic sequence of events is proposed involving changes in Cytc phosphorylation, increased ROS via increased mitochondrial membrane potentials or the p66(Shc) pathway, and oxidation of cardiolipin by Cytc followed by its release from the mitochondria. Cytc regulation in respiration and cell death is discussed in a human disease context including neurodegenerative and cardiovascular diseases, cancer, and sepsis.

Identification and analysis of unitary pseudogenes: historic and contemporary gene losses in humans and other primates

Novel human pseudogenes are identified that had previous functionality and their age is estimated. The rate of loss-of-function occurred uniformly.

Background

Unitary pseudogenes are a class of unprocessed pseudogenes without functioning counterparts in the genome. They constitute only a small fraction of annotated pseudogenes in the human genome. However, as they represent distinct functional losses over time, they shed light on the unique features of humans in primate evolution.

Results

We have developed a pipeline to detect human unitary pseudogenes through analyzing the global inventory of orthologs between the human genome and its mammalian relatives. We focus on gene losses along the human lineage after the divergence from rodents about 75 million years ago. In total, we identify 76 unitary pseudogenes, including previously annotated ones, and many novel ones. By comparing each of these to its functioning ortholog in other mammals, we can approximately date the creation of each unitary pseudogene (that is, the gene 'death date') and show that for our group of 76, the functional genes appear to be disabled at a fairly uniform rate throughout primate evolution - not all at once, correlated, for instance, with the 'Alu burst'. Furthermore, we identify 11 unitary pseudogenes that are polymorphic - that is, they have both nonfunctional and functional alleles currently segregating in the human population. Comparing them with their orthologs in other primates, we find that two of them are in fact pseudogenes in non-human primates, suggesting that they represent cases of a gene being resurrected in the human lineage.

Conclusions

This analysis of unitary pseudogenes provides insights into the evolutionary constraints faced by different organisms and the timescales of functional gene loss in humans.

Evolutionary alkaline transition in human cytochrome c

Conformational transitions in cytochrome c (cyt c) are being realized to be responsible for its multi-functions. Among a number of conformational transitions in cyt c, the alkaline transition has attracted much attention. The cDNA of human cyt c is cloned by RT-PCR and a high-effective expression system for human cyt c has been developed in this study. The equilibrium and kinetics of the alkaline transition of human cyt c have been systematically investigated for the first time, and compared with those of yeast and horse cyt c from an evolutionary perspective. The pK(a) value for the alkaline transition of human cyt c is apparently higher than that of yeast and horse. Kinetic studies suggest that it is increasingly difficult for the alkaline transition of cyt c from yeast, horse and human. Molecular modeling of human cyt c shows that the omega loop where the lysine residue is located apparently further away from heme in human cyt c than in yeast iso-1 and horse heart cyt c. These results regarding alkaline conformational transition provide valuable information for understanding the molecular basis for the biological multi-functions of cyt c.

Identification of a new rice blast resistance gene, Pid3, by genomewide comparison of paired nucleotide-binding site--leucine-rich repeat genes and their pseudogene alleles between the two sequenced rice genomes

Rice blast, caused by Magnaporthe oryzae, is one of the most devastating diseases. The two major subspecies of Asian cultivated rice (Oryza sativa L.), indica and japonica, have shown obvious differences in rice blast resistance, but the genomic basis that underlies the difference is not clear. We performed a genomewide comparison of the major class of resistant gene family, the nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family, between 93-11 (indica) and Nipponbare (japonica) with a focus on their pseudogene members. We found great differences in either constitution or distribution of pseudogenes between the two genomes. According to this comparison, we designed the PCR-based molecular markers specific to the Nipponbare NBS-LRR pseudogene alleles and used them as cosegregation markers for blast susceptibility in a segregation population from a cross between a rice blast-resistant indica variety and a susceptible japonica variety. Through this approach, we identified a new blast resistance gene, Pid3, in the indica variety, Digu. The allelic Pid3 loci in most of the tested japonica varieties were identified as pseudogenes due to a nonsense mutation at the nucleotide position 2208 starting from the translation initiation site. However, this mutation was not found in any of the tested indica varieties, African cultivated rice varieties, or AA genome-containing wild rice species. These results suggest that the pseudogenization of Pid3 in japonica occurred after the divergence of indica and japonica.

Molecular and phylogenetic characterization of cytochromes c from Haemonchus contortus and Trichostrongylus vitrinus (Nematoda: Trichostrongylida)

Although cytochrome c genes (cyt c) and proteins (CYT C) have been relatively well studied in mammals, very little is known about them in parasitic helminths. In the present study, we investigated this group of molecules in Haemonchus contortus (barber's pole worm) and Trichostrongylus vitrinus (black scour worm), two parasitic nematodes of small ruminants. The cyt c gene (512 bp) of H. contortus had one intron and encoded a transcript of 345 nucleotides, whilst that of T. vitrinus (792 bp) had two introns and encoded a transcript of 360 nucleotides. The transcription of cyt c in T. vitrinus was substantially greater in adult males compared with females, although no such gender-enrichment was evident in adults of H. contortus. These findings were supported at the protein level by immunoblot analyses. The inferred proteins (designated Hc-CYT C and Tv-CYT C, respectively) shared nucleotide and amino acid identities of 78% and 85%, respectively. The alignment of these and other CYT C sequences from nematodes, flatworms, insects and mammals identified conserved motifs associated with CYT C oxidase- and reductase- as well as haem-binding. One residue (histidine-26) was conserved for mammals, whereas this residue was absent from all nematodes; the functional significance of this difference is not yet known. Both phylogenetic analysis and protein modelling revealed that CYT C proteins of nematodes are structurally distinct from those of mammals and other organisms, suggesting their potential as targets for parasite intervention.

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.

Metazoan OXPHOS gene families: Evolutionary forces at the level of mitochondrial and nuclear genomes

Mitochondrial and nuclear DNAs contribute to encode the whole mitochondrial protein complement. The two genomes possess highly divergent features and properties, but the forces influencing their evolution, even if different, require strong coordination. The gene content of mitochondrial genome in all Metazoa is in a frozen state with only few exceptions and thus mitochondrial genome plasticity especially concerns some molecular features, i.e. base composition, codon usage, evolutionary rates. In contrast the high plasticity of nuclear genomes is particularly evident at the macroscopic level, since its redundancy represents the main feature able to introduce genetic material for evolutionary innovations. In this context, genes involved in oxidative phosphorylation (OXPHOS) represent a classical example of the different evolutionary behaviour of mitochondrial and nuclear genomes. The simple DNA sequence of Cytochrome c oxidase I (encoded by the mitochondrial genome) seems to be able to distinguish intra- and inter-species relations between organisms (DNA Barcode). Some OXPHOS subunits (cytochrome c, subunit c of ATP synthase and MLRQ) are encoded by several nuclear duplicated genes which still represent the trace of an ancient segmental/genome duplication event at the origin of vertebrates.

Remarkably high activities of testicular cytochrome c in destroying reactive oxygen species and in triggering apoptosis

Hydrogen peroxide (H(2)O(2)) is the major reactive oxygen species (ROS) produced in sperm. High concentrations of H(2)O(2) in sperm induce nuclear DNA fragmentation and lipid peroxidation and result in cell death. The respiratory chain of the mitochondrion is one of the most productive ROS generating systems in sperm, and thus the destruction of ROS in mitochondria is critical for the cell. It was recently reported that H(2)O(2) generated by the respiratory chain of the mitochondrion can be efficiently destroyed by the cytochrome c-mediated electron-leak pathway where the electron of ferrocytochrome c migrates directly to H(2)O(2) instead of to cytochrome c oxidase. In our studies, we found that mouse testis-specific cytochrome c (T-Cc) can catalyze the reduction of H(2)O(2) three times faster than its counterpart in somatic cells (S-Cc) and that the T-Cc heme has the greater resistance to being degraded by H(2)O(2). Together, these findings strongly imply that T-Cc can protect sperm from the damages caused by H(2)O(2). Moreover, the apoptotic activity of T-Cc is three to five times greater than that of S-Cc in a well established apoptosis measurement system using Xenopus egg extract. The dramatically stronger apoptotic activity of T-Cc might be important for the suicide of male germ cells, considered a physiological mechanism that regulates the number of sperm produced and eliminates those with damaged DNA. Thus, it is very likely that T-Cc has evolved to guarantee the biological integrity of sperm produced in mammalian testis.

Evolution of ATP synthase subunit c and cytochrome c gene families in selected Metazoan classes

To investigate the integrated evolution of mitochondrial and nuclear genomes in the eukaryotic cell, we have focused our attention on OXPHOS (oxidative phosphorylation) gene families which encode proteins involved in the main mitochondrial function. The present study reports the phylogenetic analysis of two OXPHOS gene families: ATP synthase subunit c (or lipid binding protein, LBP) and Cytochrome c (Cytc). Both gene families possess a higher expansion trend than the typically low duplication rate of OXPHOS genes in Metazoa, but follow a completely different evolutionary history, especially in mammals. LBP is represented by three well conserved isoforms in all mammals (P1, P2, P3): only P3 possesses a clearly conserved isoform in all Vertebrates, P1 and P2 were already present before the bird-mammal divergence and there are preliminary evidence from the in silico analysis that P1, the most evolutionary divergent isoform, is poorly expressed and not regulated by NRF1. In contrast, Cytc family presents at least two duplicated genes in all the analysed Vertebrates, is subject to a high expansion trend, especially of processed pseudogenes in mammals, and some events of gain and loss of function can be supposed.

Evolution of the NANOG pseudogene family in the human and chimpanzee genomes

Background

The NANOG gene is expressed in mammalian embryonic stem cells where it maintains cellular pluripotency. An unusually large family of pseudogenes arose from it with one unprocessed and ten processed pseudogenes in the human genome. This article compares the NANOG gene and its pseudogenes in the human and chimpanzee genomes and derives an evolutionary history of this pseudogene family.

Results

The NANOG gene and all pseudogenes except NANOGP8 are present at their expected orthologous chromosomal positions in the chimpanzee genome when compared to the human genome, indicating that their origins predate the human-chimpanzee divergence. Analysis of flanking DNA sequences demonstrates that NANOGP8 is absent from the chimpanzee genome.

Conclusion

Based on the most parsimonious ordering of inferred source-gene mutations, the deduced evolutionary origins for the NANOG pseudogene family in the human and chimpanzee genomes, in order of most ancient to most recent, are NANOGP6, NANOGP5, NANOGP3, NANOGP10, NANOGP2, NANOGP9, NANOGP7, NANOGP1, and NANOGP4. All of these pseudogenes were fixed in the genome of the human-chimpanzee common ancestor. NANOGP8 is the most recent pseudogene and it originated exclusively in the human lineage after the human-chimpanzee divergence. NANOGP1 is apparently an unprocessed pseudogene. Comparison of its sequence to the functional NANOG gene's reading frame suggests that this apparent pseudogene remained functional after duplication and, therefore, was subject to selection-driven conservation of its reading frame, and that it may retain some functionality or that its loss of function may be evolutionarily recent.

3',5'-cyclic nucleotide phosphodiesterase 11A: localization in human tissues

3',5'-Cyclic nucleotide phosphodiesterase 11 (PDE11) is the most recently discovered family of human 3',5'-cyclic nucleotide phosphodiesterases (PDEs). This family contains one gene, PDE11A, with four splice variants (PDE11A1-PDE11A4). The physiological role of PDE11A has not been determined. Tadalafil (Cialis), a PDE5A inhibitor used for the treatment of male erectile dysfunction, has been reported to partially inhibit PDE11. It was therefore of interest to consider the pattern of expression of PDE11 in human tissues. Although four PDE11A mRNA transcripts have been reported, we detected protein corresponding to only one of them, PDE11A4, in human prostate, pituitary, heart and liver. Using immunohistochemistry, there was strong PDE11A antibody staining in the glandular epithelium of the prostate and weak staining of neuronal cells within parasympathetic ganglia in the heart. No PDE11A protein was detected in blood vessels or cardiac myocytes. None of the four potential PDE11A proteins were detected in human skeletal muscle, testis, or penis.

Large-scale analysis of pseudogenes in the human genome

Pseudogenes are considered as genomic fossils: disabled copies of functional genes that were once active in the ancient genome. Recently, whole-genome computational approaches have revealed thousands of pseudogenes in the genomes of the human and other eukaryotes. Identification of these pseudogenes can improve the accuracy of gene annotation. It also offers new insight on the evolutionary history and the stability of the genome as a whole.

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.

Comparative analysis of processed pseudogenes in the mouse and human genomes

Pseudogenes are important resources in evolutionary and comparative genomics because they provide molecular records of the ancient genes that existed in the genome millions of years ago. We have systematically identified approximately 5000 processed pseudogenes in the mouse genome, and estimated that approximately 60% are lineage specific, created after the mouse and human diverged. In both mouse and human genomes, similar types of genes give rise to many processed pseudogenes. These tend to be housekeeping genes, which are highly expressed in the germ line. Ribosomal-protein genes, in particular, form the largest sub-group. The processed pseudogenes in the mouse occur with a distinctly different chromosomal distribution than LINEs or SINEs - preferentially in GC-poor regions. Finally, the age distribution of mouse-processed pseudogenes closely resembles that of LINEs, in contrast to human, where the age distribution closely follows Alus (SINEs).

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.